Showing 1–50 of 65 results for author: Cuerno, R

Kardar-Parisi-Zhang universality class in the synchronization of oscillator lattices with time-dependent noise

Authors: Ricardo Gutierrez, Rodolfo Cuerno

Abstract: Systems of oscillators subject to time-dependent noise typically achieve synchronization for long times when their mutual coupling is sufficiently strong. The dynamical process whereby synchronization is reached can be thought of as a growth process in which an interface formed by the local phase field gradually roughens and eventually saturates. Such a process is here shown to display the generic… ▽ More Systems of oscillators subject to time-dependent noise typically achieve synchronization for long times when their mutual coupling is sufficiently strong. The dynamical process whereby synchronization is reached can be thought of as a growth process in which an interface formed by the local phase field gradually roughens and eventually saturates. Such a process is here shown to display the generic scale invariance of the one-dimensional Kardar-Parisi-Zhang universality class, including a Tracy-Widom probability distribution for phase fluctuations around their mean. This is revealed by numerical explorations of a variety of oscillator systems: rings of generic phase oscillators and rings of paradigmatic limit-cycle oscillators, like Stuart-Landau and van der Pol. It also agrees with analytical expectations derived under conditions of strong mutual coupling. The nonequilibrium critical behavior that we find is robust and transcends the details of the oscillators considered. Hence, it may well be accessible to experimental ensembles of oscillators in the presence of e.g. thermal noise. △ Less

Submitted 22 July, 2024; originally announced July 2024.

Comments: 9 pages, 5 figures

Universal fluctuations of global geometrical measurements in planar clusters

Authors: Silvia N. Santalla, Iván Álvarez Domenech, Daniel Villarrubia, Rodolfo Cuerno, Javier Rodríguez-Laguna

Abstract: We characterize universal features of the sample-to-sample fluctuations of global geometrical observables, such as the area, width, length, and center-of-mass position, in random growing planar clusters. Our examples are taken from simulations of both continuous and discrete models of kinetically rough interfaces, including several universality classes, such as Kardar-Parisi-Zhang. We mostly focus… ▽ More We characterize universal features of the sample-to-sample fluctuations of global geometrical observables, such as the area, width, length, and center-of-mass position, in random growing planar clusters. Our examples are taken from simulations of both continuous and discrete models of kinetically rough interfaces, including several universality classes, such as Kardar-Parisi-Zhang. We mostly focus on the scaling behavior with time of the sample-to-sample deviation for those global magnitudes, but we have also characterized their histograms and correlations. △ Less

Submitted 26 November, 2023; originally announced November 2023.

Kardar-Parisi-Zhang fluctuations in the synchronization dynamics of limit-cycle oscillators

Authors: Ricardo Gutiérrez, Rodolfo Cuerno

Abstract: The space-time process whereby one-dimensional systems of self-sustained oscillators synchronize is shown to display generic scale invariance, with scaling properties characteristic of the Kardar-Parisi-Zhang equation with columnar noise, and phase fluctuations that follow a Tracy-Widom probability distribution. This is revealed by a numerical exploration of rings of Stuart-Landau oscillators (the… ▽ More The space-time process whereby one-dimensional systems of self-sustained oscillators synchronize is shown to display generic scale invariance, with scaling properties characteristic of the Kardar-Parisi-Zhang equation with columnar noise, and phase fluctuations that follow a Tracy-Widom probability distribution. This is revealed by a numerical exploration of rings of Stuart-Landau oscillators (the universal representation of an oscillating system close to a Hopf bifurcation) and rings of van der Pol oscillators, both of which are paradigms of self-sustained oscillators. The critical behavior is very well-defined for limit-cycle oscillations near the bifurcation point, and still dominates the behavior comparatively far from the bifurcation. In particular, the Tracy-Widom fluctuation distribution seems to be an extremely robust feature of the synchronization process. The nonequilibrium criticality here described appears to transcend the details of the coupled dynamical systems that synchronize, making plausible its experimental observation. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 11 pages, 5 figures

Shape effects in the fluctuations of random isochrones on a square lattice

Authors: Iván Álvarez Domenech, Javier Rodríguez-Laguna, Rodolfo Cuerno, Pedro Córdoba-Torres, Silvia N. Santalla

Abstract: We consider the isochrone curves in first-passage percolation on a 2D square lattice, i.e. the boundary of the set of points which can be reached in less than a given time from a certain origin. The occurrence of an instantaneous average shape is described in terms of its Fourier components, highlighting a crossover between a diamond and a circular geometry as the noise level is increased. General… ▽ More We consider the isochrone curves in first-passage percolation on a 2D square lattice, i.e. the boundary of the set of points which can be reached in less than a given time from a certain origin. The occurrence of an instantaneous average shape is described in terms of its Fourier components, highlighting a crossover between a diamond and a circular geometry as the noise level is increased. Generally, these isochrones can be understood as fluctuating interfaces with an inhomogeneous local width which reveals the underlying lattice structure. We show that once these inhomogeneities have been taken into account, the fluctuations fall into the Kardar-Parisi-Zhang (KPZ) universality class with very good accuracy, where they reproduce the Family-Vicsek Ansatz with the expected exponents and the Tracy-Widom histogram for the local radial fluctuations. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Universal interface fluctuations in the contact process

Authors: B. G. Barreales, J. J. Meléndez, R. Cuerno, J. J. Ruiz-Lorenzo

Abstract: We study the interface representation of the contact process (CP) at its directed-percolation critical point, where the scaling properties of the interface can be related to those of the original particle model. Interestingly, such a behavior happens to be intrinsically anomalous and more complex than that described by the standard Family-Vicsek dynamic scaling Ansatz of surface kinetic roughening… ▽ More We study the interface representation of the contact process (CP) at its directed-percolation critical point, where the scaling properties of the interface can be related to those of the original particle model. Interestingly, such a behavior happens to be intrinsically anomalous and more complex than that described by the standard Family-Vicsek dynamic scaling Ansatz of surface kinetic roughening. We expand on a previous numerical study by Dickman and Muñoz [Phys. Rev. E 62, 7632 (2000)] to fully characterize the kinetic roughening universality class for interface dimensions $d=1, 2$, and 3. Beyond obtaining scaling exponent values, we characterize the interface fluctuations via their probability density function (PDF) and covariance, seen to display universal properties which are qualitatively similar to those recently assessed for the Kardar-Parisi-Zhang (KPZ) and other important universality classes of kinetic roughening. Quantitatively, while for $d=1$ the interface covariance seems to be well described by the KPZ, Airy$_1$ covariance, no such agreement occurs in terms of the fluctuation PDF nor the scaling exponents. △ Less

Submitted 21 April, 2023; originally announced April 2023.

Comments: 12 pages, 9 figures

Journal ref: Physical Review E 108, 044801 (2023)

Nonequilibrium criticality driven by Kardar-Parisi-Zhang fluctuations in the synchronization of oscillator lattices

Authors: Ricardo Gutiérrez, Rodolfo Cuerno

Abstract: The synchronization of oscillator ensembles is pervasive throughout nonlinear science, from classical or quantum mechanics to biology, to human assemblies. Traditionally, the main focus has been the identification of threshold parameter values for the transition to synchronization as well as the nature of such transition. Here, we show that considering an oscillator lattice as a discrete growing i… ▽ More The synchronization of oscillator ensembles is pervasive throughout nonlinear science, from classical or quantum mechanics to biology, to human assemblies. Traditionally, the main focus has been the identification of threshold parameter values for the transition to synchronization as well as the nature of such transition. Here, we show that considering an oscillator lattice as a discrete growing interface provides unique insights into the dynamical process whereby the lattice reaches synchronization for long times. Working on a generalization of the Kuramoto model that allows for odd or non-odd couplings, we elucidate synchronization of oscillator lattices as an instance of generic scale invariance, whereby the system displays space-time criticality, largely irrespective of parameter values. The critical properties of the system (like scaling exponent values and the dynamic scaling Ansatz which is satisfied) fall into universality classes of kinetically rough interfaces with columnar disorder, namely, those of the Edwards-Wilkinson (equivalently, the Larkin model of an elastic interface in a random medium) or the Kardar-Parisi-Zhang (KPZ) equations, for Kuramoto (odd) coupling and generic (non-odd) couplings, respectively. From the point of view of kinetic roughening, the critical properties we find turn out to be quite innovative, especially concerning the statistics of the fluctuations as characterized by their probability distribution function (PDF) and covariance. While the latter happens to be that of the Larkin model irrespective of the symmetry of the coupling, in the generic non-odd coupling case the PDF turns out to be the Tracy-Widom distribution associated with the KPZ nonlinearity. This brings the synchronization of oscillator lattices into a remarkably large class of strongly-correlated, low-dimensional (classical and quantum) systems with strong universal fluctuations. △ Less

Submitted 21 April, 2023; v1 submitted 7 October, 2022; originally announced October 2022.

Comments: 26 pages, 18 figures

Journal ref: Phys. Rev. Res. 5, 023047 (2023)

Anomalous dynamic scaling of Ising interfaces

Authors: E. Rodriguez-Fernandez, S. N. Santalla, M. Castro, R. Cuerno

Abstract: Until very recently, the asymptotic occurrence of intrinsic anomalous scaling has been expected to require concomitant effects for kinetically rough interfaces, like quenched disorder or morphological instabilities. However, counterexamples have been recently reported for simpler situations dominated by time-dependent noise, as in the discrete growth system associated with an Ising model proposed… ▽ More Until very recently, the asymptotic occurrence of intrinsic anomalous scaling has been expected to require concomitant effects for kinetically rough interfaces, like quenched disorder or morphological instabilities. However, counterexamples have been recently reported for simpler situations dominated by time-dependent noise, as in the discrete growth system associated with an Ising model proposed by H.\ Dashti-Naserabadi {\em et al.}\ [Phys.\ Rev.\ E {\bf 100}, 060101(R) (2019)], who assessed the equilibrium behavior of the model. Here we revisit this system to characterize its time-dependent behavior in two and three dimensions (one-and two-dimensional interfaces, respectively). While the 3D case seems dominated by a fast evolution beyond critical dynamics, in the 2D case numerical simulations of an associated time-dependent Ginzburg-Landau equation retrieve the same static (roughness) exponents and the same intrinsic anomalous scaling Ansatz as in the equilibrium case, throughout the full time evolution. However, the dynamic exponent is seen to cross over between two different values, none of which enables identification with previously known universality classes of kinetic roughening. Simulations for larger system sizes moreover suggest breakdown of scaling behavior at the largest scales, suggesting that the previously reported scaling behavior may be effective and restricted to relatively small systems. △ Less

Submitted 14 May, 2024; v1 submitted 1 August, 2022; originally announced August 2022.

Comments: 13 pages, 12 figures

Large-scale kinetic roughening behavior of coffee-ring fronts

Authors: B. G. Barreales, J. J. Melendez, R. Cuerno, J. J. Ruiz-Lorenzo

Abstract: We have studied the kinetic roughening behavior of the fronts of coffee-ring aggregates via extensive numerical simulations of the off-lattice model considered for this context in [C.\ S.\ Dias {\it et al.}, Soft Matter {\bf 14}, 1903 (2018)]. This model describes ballistic aggregation of patchy colloids and depends on a parameter $r_\mathrm{AB}$ which controls the affinity of the two patches, A a… ▽ More We have studied the kinetic roughening behavior of the fronts of coffee-ring aggregates via extensive numerical simulations of the off-lattice model considered for this context in [C.\ S.\ Dias {\it et al.}, Soft Matter {\bf 14}, 1903 (2018)]. This model describes ballistic aggregation of patchy colloids and depends on a parameter $r_\mathrm{AB}$ which controls the affinity of the two patches, A and B. Suitable boundary conditions allow us to elucidate a discontinuous pinning-depinning transition at $r_\mathrm{AB}=0$, with the front displaying intrinsic anomalous scaling, but with unusual exponent values $α\simeq 1.2$, $α_{\rm loc} \simeq 0.5$, $β\simeq 1$, and $z\simeq 1.2$. For $0<r_\mathrm{AB}\le 1$, comparison with simulations of standard off-lattice ballistic deposition indicates the occurrence of a morphological instability induced by the patch structure. As a result, we find that the asymptotic morphological behavior is dominated by macroscopic shapes. The intermediate time regime exhibits one-dimensional KPZ exponents for $r_\mathrm{AB}> 0.01$ and the system suffers a strong crossover dominated by the $r_\mathrm{AB}=0$ behavior for $r_\mathrm{AB}\le 0.01$. A detailed analysis of correlation functions shows that the aggregate fronts are always in the moving phase for $0<r_\mathrm{AB}\le 1$ and that their kinetic roughening behavior is intrinsically anomalous for $r_\mathrm{AB}\le 0.01$. △ Less

Submitted 11 October, 2022; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: 16 pages and 20 figures. Final version published in PRE

Journal ref: Phys. Rev. E 106, 044801 (2022)

Anomalous ballistic scaling in the tensionless or inviscid Kardar-Parisi-Zhang equation

Authors: Enrique Rodriguez-Fernandez, Silvia N. Santalla, Mario Castro, Rodolfo Cuerno

Abstract: The one-dimensional Kardar-Parisi-Zhang (KPZ) equation is becoming an overarching paradigm for the scaling of nonequilibrium, spatially extended, classical and quantum systems with strong correlations. Recent analytical solutions have uncovered a rich structure regarding its scaling exponents and fluctuation statistics. However, the zero surface tension or zero viscosity case eludes such analytica… ▽ More The one-dimensional Kardar-Parisi-Zhang (KPZ) equation is becoming an overarching paradigm for the scaling of nonequilibrium, spatially extended, classical and quantum systems with strong correlations. Recent analytical solutions have uncovered a rich structure regarding its scaling exponents and fluctuation statistics. However, the zero surface tension or zero viscosity case eludes such analytical solutions and has remained ill-understood. Using numerical simulations, we elucidate a well-defined universality class for this case that differs from that of the viscous case, featuring intrinsically anomalous kinetic roughening, despite previous expectations for systems with local interactions and time-dependent noise and ballistic dynamics. The latter may be relevant to recent quantum spin chain experiments which measure KPZ and ballistic relaxation under different conditions. We identify the ensuing set of scaling exponents in previous discrete interface growth models related with isotropic percolation, and show it to describe the fluctuations of additional continuum systems related with the noisy Korteweg-de Vries equation. Along this process, we additionally elucidate the universality class of the related inviscid stochastic Burgers equation. △ Less

Submitted 25 July, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 9 pages, 8 figures, Physical Review E (in press, 2022)

Spreading fronts of wetting liquid droplets: microscopic simulations and universal fluctuations

Authors: J. M. Marcos, P. Rodríguez-López, J. J. Melendez, R. Cuerno, J. J. Ruiz-Lorenzo

Abstract: We have used kinetic Monte Carlo (kMC) simulations of a lattice gas to study front fluctuations in the spreading of a non-volatile liquid droplet onto a solid substrate. Our results are consistent with a diffusive growth law for the radius of the precursor layer, $R \sim t^δ$, with $δ\approx 1/2$ in all the conditions considered for temperature and substrate wettability, in good agreement with pre… ▽ More We have used kinetic Monte Carlo (kMC) simulations of a lattice gas to study front fluctuations in the spreading of a non-volatile liquid droplet onto a solid substrate. Our results are consistent with a diffusive growth law for the radius of the precursor layer, $R \sim t^δ$, with $δ\approx 1/2$ in all the conditions considered for temperature and substrate wettability, in good agreement with previous studies. The fluctuations of the front exhibit kinetic roughening properties with exponent values which depend on temperature $T$, but become $T$-independent for sufficiently high $T$. Moreover, strong evidences of intrinsic anomalous scaling have been found, characterized by different values of the roughness exponent at short and large length scales. Although such a behavior differs from the scaling properties of the one-dimensional Kardar-Parisi-Zhang (KPZ) universality class, the front covariance and the probability distribution function of front fluctuations found in our kMC simulations do display KPZ behavior, agreeing with simulations of a continuum height equation proposed in this context. However, this equation does not feature intrinsic anomalous scaling, at variance with the discrete model. △ Less

Submitted 22 April, 2022; originally announced April 2022.

Comments: 15 pages and 14 figures. To be published in PRE

Journal ref: Physical Review E 105, 054801 (2022)

Transition between chaotic and stochastic universality classes of kinetic roughening

Authors: E. Rodriguez-Fernandez, R Cuerno

Abstract: The dynamics of non-equilibrium spatially extended systems are often dominated by fluctuations, due to e.g.\ deterministic chaos or to intrinsic stochasticity. This reflects into generic scale invariant or kinetic roughening behavior that can be classified into universality classes defined by critical exponent values and by the probability distribution function (PDF) of field fluctuations. Suitabl… ▽ More The dynamics of non-equilibrium spatially extended systems are often dominated by fluctuations, due to e.g.\ deterministic chaos or to intrinsic stochasticity. This reflects into generic scale invariant or kinetic roughening behavior that can be classified into universality classes defined by critical exponent values and by the probability distribution function (PDF) of field fluctuations. Suitable geometrical constraints are known to change secondary features of the PDF while keeping the values of the exponents unchanged, inducing universality subclasses. Working on the Kuramoto-Sivashinsky equation as a paradigm of spatiotemporal chaos, we show that the physical nature of the prevailing fluctuations (chaotic or stochastic) can also change the universality class while respecting the exponent values, as the PDF is substantially altered. This transition takes place at a non-zero value of the stochastic noise amplitude and may be suitable for experimental verification. △ Less

Submitted 24 September, 2020; originally announced September 2020.

Comments: 12 pages, 9 figures, submitted

Journal ref: Phys. Rev. Research 3, 012020 (2021)

Kardar-Parisi-Zhang universality class for the critical dynamics of reaction-diffusion fronts

Authors: B. G. Barreales, J. J. Melendez, R. Cuerno, J. J. Ruiz-Lorenzo

Abstract: We have studied front dynamics for the discrete $A+A \leftrightarrow A$ reaction-diffusion system, which in the continuum is described by the (stochastic) Fisher-Kolmogorov-Petrovsky-Piscunov equation. We have revisited this discrete model in two space dimensions by means of extensive numerical simulations and an improved analysis of the time evolution of the interface separating the stable and un… ▽ More We have studied front dynamics for the discrete $A+A \leftrightarrow A$ reaction-diffusion system, which in the continuum is described by the (stochastic) Fisher-Kolmogorov-Petrovsky-Piscunov equation. We have revisited this discrete model in two space dimensions by means of extensive numerical simulations and an improved analysis of the time evolution of the interface separating the stable and unstable phases. In particular, we have measured the full set of critical exponents which characterize the spatio-temporal fluctuations of such front for different lattice sizes, focusing mainly in the front width and correlation length. These exponents are in very good agreement with those computed in [E. Moro, Phys. Rev. Lett. 87, 238303 (2001)] and correspond to those of the Kardar-Parisi-Zhang (KPZ) universality class for one-dimensional interfaces. Furthermore, we have studied the one-point statistics and the covariance of rescaled front fluctuations, which had remained thus far unexplored in the literature and allows for a further stringent test of KPZ universality. △ Less

Submitted 14 October, 2019; originally announced October 2019.

Comments: 19 pages and 9 figures

Journal ref: J. Stat. Mech. (2020) 023203

The derivative of the Kardar-Parisi-Zhang equation is not in the KPZ universality class

Authors: E. Rodriguez-Fernandez, R. Cuerno

Abstract: The Kardar-Parisi-Zhang (KPZ) equation is a paradigmatic model of nonequilibrium low-dimensional systems with spatiotemporal scale invariance, recently highlighting universal behavior in fluctuation statistics. Its space derivative, namely the noisy Burgers equation, has played a very important role in its study, predating the formulation of the KPZ equation proper, and being frequently held as an… ▽ More The Kardar-Parisi-Zhang (KPZ) equation is a paradigmatic model of nonequilibrium low-dimensional systems with spatiotemporal scale invariance, recently highlighting universal behavior in fluctuation statistics. Its space derivative, namely the noisy Burgers equation, has played a very important role in its study, predating the formulation of the KPZ equation proper, and being frequently held as an equivalent system. We show that, while differences in the scaling exponents for the two equations are indeed due to a mere space derivative, the fluctuations behave in a remarkably different way: while KPZ displays Tracy-Widom statistics, its derivative displays Gaussian behavior, hence being in a different universality class. We reach this conclusion via direct numerical simulations of the equations, supported by a dynamic renormalization group analysis of field statistics. △ Less

Submitted 2 March, 2020; v1 submitted 30 April, 2019; originally announced May 2019.

Comments: 13 pages, 11 figures

Journal ref: Phys. Rev. E 101, 052126 (2020)

Gaussian statistics as an emergent symmetry of the stochastic Burgers equation

Authors: Enrique Rodriguez-Fernandez, Rodolfo Cuerno

Abstract: Symmetries play a conspicuous role in the large-scale behavior of critical systems. While in equilibrium they allow to classify asymptotics into different universality classes, out of equilibrium they can emerge, some times unexpectedly, as collective properties which are not explicit in the 'bare' interactions. Here we elucidate the emergence of an up-down symmetry in the asymptotic behavior of t… ▽ More Symmetries play a conspicuous role in the large-scale behavior of critical systems. While in equilibrium they allow to classify asymptotics into different universality classes, out of equilibrium they can emerge, some times unexpectedly, as collective properties which are not explicit in the 'bare' interactions. Here we elucidate the emergence of an up-down symmetry in the asymptotic behavior of the stochastic scalar Burgers equation in one and two dimensions, manifested by the occurrence of Gaussian fluctuations for the physical field. This robustness of Gaussian behavior contradicts naive expectations, due to the detailed relation ---including the same set of symmetries--- between Burgers equation and the Kardar-Parisi-Zhang equation, which paradigmatically displays non-Gaussian fluctuations described by Tracy-Widom distributions. We reach our conclusions via a dynamic renormalization group study of the field statistics, confirmed by direct evaluation of the field probability distribution function from numerical simulations of the dynamical equation. All odd-order cumulants cancel exactly, while the excess kurtosis vanishes for large systems, indeed consistent with Gaussian behavior. △ Less

Submitted 25 April, 2019; v1 submitted 6 September, 2018; originally announced September 2018.

Comments: 6 pages, 3 figures

Journal ref: Phys. Rev. E 99, 042108 (2019)

Kardar-Parisi-Zhang Universality in First-Passage Percolation: the Role of Geodesic Degeneracy

Authors: Pedro Córdoba-Torres, Silvia N. Santalla, Rodolfo Cuerno, Javier Rodríguez-Laguna

Abstract: We have characterized the scaling behavior of the first-passage percolation (FPP) model on two types of discrete networks, the regular square lattice and the disordered Delaunay lattice, thereby addressing the effect of the underlying topology. Several distribution functions for the link-times were considered. The asymptotic behavior of the fluctuations for both the minimal arrival time and the la… ▽ More We have characterized the scaling behavior of the first-passage percolation (FPP) model on two types of discrete networks, the regular square lattice and the disordered Delaunay lattice, thereby addressing the effect of the underlying topology. Several distribution functions for the link-times were considered. The asymptotic behavior of the fluctuations for both the minimal arrival time and the lateral deviation of the geodesic path are in perfect agreement with the Kardar-Parisi-Zhang (KPZ) universality class regardless of the type of the link-time distribution and of the lattice topology. Pre-asymptotic behavior, on the other hand, is found to depend on the uniqueness of geodesics in absence of disorder in the local crossing times, a topological property of lattice directions that we term geodesic degeneracy. This property has important consequences on the model, as for example the well-known anisotropic growth in regular lattices. In this work we provide a framework to understand its effect as well as to characterize its extent. △ Less

Submitted 14 February, 2018; originally announced February 2018.

Journal ref: J. Stat. Mech. (2018) 063212

Non-universality of front fluctuations for compact colonies of non-motile bacteria

Authors: Silvia N. Santalla, Javier Rodríguez-Laguna, José P. Abad, Irma Marín, María del Mar Espinosa, Javier Muñoz-García, Luis Vázquez, Rodolfo Cuerno

Abstract: The front of a compact bacterial colony growing on a Petri dish is a paradigmatic instance of non-equilibrium fluctuations in the celebrated Eden, or Kardar-Parisi-Zhang (KPZ), universality class. While in many experiments the scaling exponents crucially differ from the expected KPZ values, the source of this disagreement has remained poorly understood. We have performed growth experiments with B.… ▽ More The front of a compact bacterial colony growing on a Petri dish is a paradigmatic instance of non-equilibrium fluctuations in the celebrated Eden, or Kardar-Parisi-Zhang (KPZ), universality class. While in many experiments the scaling exponents crucially differ from the expected KPZ values, the source of this disagreement has remained poorly understood. We have performed growth experiments with B. subtilis 168 and E. coli ATCC 25922 under conditions leading to compact colonies in the classically-alleged Eden regime, where individual motility is suppressed. Non-KPZ scaling is indeed observed for all accessible times, KPZ asymptotics being ruled out for our experiments due to the monotonic increase of front branching with time. Simulations of an effective model suggest the occurrence of transient non-universal scaling due to diffusive morphological instabilities, agreeing with expectations from detailed models of the relevant biological reaction-diffusion processes. △ Less

Submitted 5 July, 2018; v1 submitted 13 February, 2017; originally announced February 2017.

Journal ref: Phys. Rev. E 98, 012407 (2018)

Topology and the Kardar-Parisi-Zhang universality class

Authors: Silvia N. Santalla, Javier Rodriguez-Laguna, Alessio Celi, Rodolfo Cuerno

Abstract: We study the role of the topology of the background space on the one-dimensional Kardar-Parisi-Zhang (KPZ) universality class. To do so, we study the growth of balls on disordered 2D manifolds with random Riemannian metrics, generated by introducing random perturbations to a base manifold. As base manifolds we consider cones of different aperture angles $θ$, including the limiting cases of a cylin… ▽ More We study the role of the topology of the background space on the one-dimensional Kardar-Parisi-Zhang (KPZ) universality class. To do so, we study the growth of balls on disordered 2D manifolds with random Riemannian metrics, generated by introducing random perturbations to a base manifold. As base manifolds we consider cones of different aperture angles $θ$, including the limiting cases of a cylinder ($θ=0$, which corresponds to an interface with periodic boundary conditions) and a plane ($θ=π/2$, which corresponds to an interface with circular geometry). We obtain that in the former case the radial fluctuations of the ball boundaries follow the Tracy-Widom (TW) distribution of the largest eigenvalue of random matrices in the Gaussian orthogonal ensemble (TW-GOE), while on cones with any aperture angle $θ\neq 0$ fluctuations correspond to the TW-GUE distribution related with the Gaussian unitary ensemble. We provide a topological argument to justify the relevance of TW-GUE statistics for cones, and state a conjecture which relates the KPZ universality subclass with the background topology. △ Less

Submitted 26 April, 2016; v1 submitted 16 April, 2016; originally announced April 2016.

Ion-beam nanopatterning of silicon surfaces under co-deposition of non-silicide-forming impurities

Authors: B. Moon, S. Yoo, J. -S. Kim, S. J. Kang, J. Muñoz-García, R. Cuerno

Abstract: We report experiments on surface nanopatterning of Si targets which are irradiated with 2 keV Ar + ions impinging at near-glancing incidence, under concurrent co-deposition of Au impurities simultaneously extracted from a gold target by the same ion beam. Previous recent experiments by a number of groups suggest that silicide formation is a prerequisite for pattern formation in the presence of met… ▽ More We report experiments on surface nanopatterning of Si targets which are irradiated with 2 keV Ar + ions impinging at near-glancing incidence, under concurrent co-deposition of Au impurities simultaneously extracted from a gold target by the same ion beam. Previous recent experiments by a number of groups suggest that silicide formation is a prerequisite for pattern formation in the presence of metallic impurities. In spite of the fact that Au is known not to form stable compounds with the Si atoms, ripples nonetheless emerge in our experiments with nanometric wavelengths and small amplitudes, and with an orientation that changes with distance to the Au source. We provide results of sample analysis through Auger electron and energy-dispersive X-ray spectroscopies for their space-resolved chemical composition, and through atomic force, scanning transmission electron, and high-resolution transmission microscopies for their morphological properties. We discuss these findings in the light of current continuum models for this class of systems. The composition of and the dynamics within the near-surface amorphized layer that ensues is expected to play a relevant role to account for the unexpected formation of these surface structures. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Comments: 11 pages, 9 figures; submitted

Journal ref: Phys. Rev. B 93, 115430 (2016)

Symmetry of surface nanopatterns induced by ion-beam sputtering: the role of anisotropic surface diffusion

Authors: Javier Renedo, Javier Muñoz-García, Mario Castro, Rodolfo Cuerno

Abstract: Ion Beam Sputtering (IBS) is a cost-effective technique able to produce ordered nanopatterns on the surfaces of different materials. To date, most theoretical studies of this process have focused on systems which become amorphous under irradiation, e.g. semiconductors at room temperature. Thus, in spite of the large amount of experimental work on metals, or more recently on semiconductors at high… ▽ More Ion Beam Sputtering (IBS) is a cost-effective technique able to produce ordered nanopatterns on the surfaces of different materials. To date, most theoretical studies of this process have focused on systems which become amorphous under irradiation, e.g. semiconductors at room temperature. Thus, in spite of the large amount of experimental work on metals, or more recently on semiconductors at high temperatures, such experimental contexts have received relatively little theoretical attention. These systems are characterized by transport mechanisms, e.g. surface diffusion, which are anisotropic as a reflection of the crystalline structure not being overruled by the irradiation. Here, we generalize a previous continuum theory of IBS at normal incidence, in order to account for anisotropic surface diffusion. We explore systematically our generalized model in order to understand the role of anisotropy in the space ordering properties of the resulting patterns. In particular, we derive a height equation which predicts morphological transitions among hexagonal and rectangular patterns as a function of system parameters and employ an angular correlation function to assess these pattern symmetries. By suitably choosing experimental conditions, it is found that one might be able to experimentally control the type of order displayed by the patterns produced. △ Less

Submitted 11 January, 2016; originally announced January 2016.

Comments: 13 pages, 6 figures; submitted

Universal behavior of crystalline membranes: crumpling transition and Poisson ratio of the flat phase

Authors: R. Cuerno, R Gallardo Caballero, A. Gordillo-Guerrero, P. Monroy, J. J. Ruiz-Lorenzo

Abstract: We revisit the universal behavior of crystalline membranes at and below the crumpling transition, which pertains to the mechanical properties of important soft and hard matter materials, such as the cytoskeleton of red blood cells or graphene. Specifically, we perform large-scale Monte Carlo simulations of a triangulated two-dimensional phantom network which is freely fluctuating in three-dimensio… ▽ More We revisit the universal behavior of crystalline membranes at and below the crumpling transition, which pertains to the mechanical properties of important soft and hard matter materials, such as the cytoskeleton of red blood cells or graphene. Specifically, we perform large-scale Monte Carlo simulations of a triangulated two-dimensional phantom network which is freely fluctuating in three-dimensional space. We obtain a continuous crumpling transition characterized by critical exponents which we estimate accurately through the use of finite-size techniques. By controlling the scaling corrections, we additionally compute with high accuracy the asymptotic value of the Poisson ratio in the flat phase, thus characterizing the auxetic properties of this class of systems. We obtain agreement with the value which is universally expected for polymerized membranes with a fixed connectivity. △ Less

Submitted 5 February, 2016; v1 submitted 27 November, 2015; originally announced November 2015.

Comments: 10 pages, 10 figures. Final version to be published in PRE

Journal ref: Phys. Rev. E 93, 022111 (2016)

Role of Thermal Fluctuations in Nonlinear Thin Film Dewetting

Authors: S. Nesic, R. Cuerno, E. Moro, L. Kondic

Abstract: The spontaneous formation of droplets via dewetting of a thin fluid film from a solid substrate allows for materials nanostructuring, under appropriate experimental control. While thermal fluctuations are expected to play a role in this process, their relevance has remained poorly understood, particularly during the nonlinear stages of evolution. Within a stochastic lubrication framework, we show… ▽ More The spontaneous formation of droplets via dewetting of a thin fluid film from a solid substrate allows for materials nanostructuring, under appropriate experimental control. While thermal fluctuations are expected to play a role in this process, their relevance has remained poorly understood, particularly during the nonlinear stages of evolution. Within a stochastic lubrication framework, we show that thermal noise speeds up and substantially influences the formation and evolution of the droplet arrangement. As compared with their deterministic counterparts, for a fixed spatial domain, stochastic systems feature a smaller number of droplets, with a larger variability in sizes and space distribution. Finally, we discuss the influence of stochasticity on droplet coarsening for very long times. △ Less

Submitted 20 May, 2015; originally announced May 2015.

Comments: 6 pages, 5 figures, 2 movies available upon request

Random geometry and the Kardar-Parisi-Zhang universality class

Authors: Silvia N. Santalla, Javier Rodriguez-Laguna, Tom LaGatta, Rodolfo Cuerno

Abstract: We consider a model of a quenched disordered geometry in which a random metric is defined on ${\mathbb R}^2$, which is flat on average and presents short-range correlations. We focus on the statistical properties of balls and geodesics, i.e., circles and straight lines. We show numerically that the roughness of a ball of radius $R$ scales as $R^χ$, with a fluctuation exponent $χ\simeq 1/3$, while… ▽ More We consider a model of a quenched disordered geometry in which a random metric is defined on ${\mathbb R}^2$, which is flat on average and presents short-range correlations. We focus on the statistical properties of balls and geodesics, i.e., circles and straight lines. We show numerically that the roughness of a ball of radius $R$ scales as $R^χ$, with a fluctuation exponent $χ\simeq 1/3$, while the lateral spread of the minimizing geodesic between two points at a distance $L$ grows as $L^ξ$, with wandering exponent value $ξ\simeq 2/3$. Results on related first-passage percolation (FPP) problems lead us to postulate that the statistics of balls in these random metrics belong to the Kardar-Parisi-Zhang (KPZ) universality class of surface kinetic roughening, with $ξ$ and $χ$ relating to critical exponents characterizing a corresponding interface growth process. Moreover, we check that the one-point and two-point correlators converge to the behavior expected for the Airy-2 process characterized by the Tracy-Widom (TW) probability distribution function of the largest eigenvalue of large random matrices in the Gaussian unitary ensemble (GUE). Nevertheless extreme-value statistics of ball coordinates are given by the TW distribution associated with random matrices in the Gaussian orthogonal ensemble. Furthermore, we also find TW-GUE statistics with good accuracy in arrival times. △ Less

Submitted 11 February, 2015; v1 submitted 1 July, 2014; originally announced July 2014.

Journal ref: New J. Phys. 17 033018 (2015)

The circular Kardar-Parisi-Zhang equation as an inflating, self-avoiding ring polymer

Authors: Silvia N. Santalla, Javier Rodriguez-Laguna, Rodolfo Cuerno

Abstract: We consider the Kardar-Parisi-Zhang (KPZ) equation for a circular interface in two dimensions, unconstrained by the standard small-slopes and no-overhang approximations. Numerical simulations using an adaptive scheme allow us to elucidate the complete time evolution as a crossover between a short-time regime with the interface fluctuations of a self-avoiding ring or 2D vesicle, and a long-time reg… ▽ More We consider the Kardar-Parisi-Zhang (KPZ) equation for a circular interface in two dimensions, unconstrained by the standard small-slopes and no-overhang approximations. Numerical simulations using an adaptive scheme allow us to elucidate the complete time evolution as a crossover between a short-time regime with the interface fluctuations of a self-avoiding ring or 2D vesicle, and a long-time regime governed by the Tracy-Widom distribution expected for this geometry. For small noise amplitudes, scaling behavior is only of the latter type. Large noise is also seen to renormalize the bare physical parameters of the ring, akin to analogous parameter renormalization for equilibrium 3D membranes. Our results bear particular importance on the relation between relevant universality classes of scale-invariant systems in two dimensions. △ Less

Submitted 30 December, 2013; originally announced December 2013.

Comments: Accepted at Phys. Rev. E (rapid communications)

Journal ref: Phys. Rev. E 89, 010401(R) (2014)

Strong anisotropy in two-dimensional surfaces with generic scale invariance: Non-linear effects

Authors: Edoardo Vivo, Matteo Nicoli, Rodolfo Cuerno

Abstract: We expand a previous study [Phys. Rev. E 86, 051611 (2012)] on the conditions for occurrence of strong anisotropy (SA) in the scaling properties of two-dimensional surfaces displaying generic scale invariance. There, a natural Ansatz was proposed for SA, which arises naturally when analyzing data from e.g. thin-film production experiments. The Ansatz was tested in Gaussian (linear) models of surfa… ▽ More We expand a previous study [Phys. Rev. E 86, 051611 (2012)] on the conditions for occurrence of strong anisotropy (SA) in the scaling properties of two-dimensional surfaces displaying generic scale invariance. There, a natural Ansatz was proposed for SA, which arises naturally when analyzing data from e.g. thin-film production experiments. The Ansatz was tested in Gaussian (linear) models of surface dynamics and in non-linear models, like the Hwa-Kardar (HK) equation [Phys. Rev. Lett. 62, 1813 (1989)], which are susceptible of accurate approximations through the former. In contrast, here we analyze non-linear equations for which such type of approximations fail. Working within generically-scale-invariant situations, and as representative case studies, we formulate and study a generalization of the HK equation for conserved dynamics, and reconsider well-known systems, such as the conserved and the non-conserved anisotropic Kardar-Parisi-Zhang equations. Through the combined use of Dynamic Renormalization Group analysis and direct numerical simulations, we conclude that the occurrence of SA in two-dimensional surfaces requires dynamics to be conserved. We find that, moreover, SA is not generic in parameter space but requires, rather, specific shapes from the terms appearing in the equation of motion, whose justification needs detailed information on the dynamical process that is being modeled in each particular case. △ Less

Submitted 29 November, 2013; originally announced November 2013.

Dimensional fragility of the Kardar-Parisi-Zhang universality class

Authors: Matteo Nicoli, Rodolfo Cuerno, Mario Castro

Abstract: We assess the dependence on substrate dimensionality of the asymptotic scaling behavior of a whole family of equations that feature the basic symmetries of the Kardar-Parisi-Zhang (KPZ) equation. Even for cases in which, as expected from universality arguments, these models display KPZ values for the critical exponents and limit distributions, their behavior deviates from KPZ scaling for increasin… ▽ More We assess the dependence on substrate dimensionality of the asymptotic scaling behavior of a whole family of equations that feature the basic symmetries of the Kardar-Parisi-Zhang (KPZ) equation. Even for cases in which, as expected from universality arguments, these models display KPZ values for the critical exponents and limit distributions, their behavior deviates from KPZ scaling for increasing system dimensions. Such a fragility of KPZ universality contradicts naive expectations, and questions straightforward application of universality principles for the continuum description of experimental systems. △ Less

Submitted 3 August, 2013; v1 submitted 22 April, 2013; originally announced April 2013.

Journal ref: J. Stat. Mech. Theor. Exp., P11001 (2013)

Strong anisotropy in surface kinetic roughening: analysis and experiments

Authors: Edoardo Vivo, Matteo Nicoli, Martin Engler, Thomas Michely, Luis Vázquez, Rodolfo Cuerno

Abstract: We report an experimental assessment of surface kinetic roughening properties that are anisotropic in space. Working for two specific instances of silicon surfaces irradiated by ion-beam sputtering under diverse conditions (with and without concurrent metallic impurity codeposition), we verify the predictions and consistency of a recently proposed scaling Ansatz for surface observables like the tw… ▽ More We report an experimental assessment of surface kinetic roughening properties that are anisotropic in space. Working for two specific instances of silicon surfaces irradiated by ion-beam sputtering under diverse conditions (with and without concurrent metallic impurity codeposition), we verify the predictions and consistency of a recently proposed scaling Ansatz for surface observables like the two-dimensional (2D) height Power Spectral Density (PSD). In contrast with other formulations, this Ansatz is naturally tailored to the study of two-dimensional surfaces, and allows to readily explore the implications of anisotropic scaling for other observables, such as real-space correlation functions and PSD functions for 1D profiles of the surface. Our results confirm that there are indeed actual experimental systems whose kinetic roughening is strongly anisotropic, as consistently described by this scaling analysis. In the light of our work, some types of experimental measurements are seen to be more affected by issues like finite space resolution effects, etc. that may hinder a clear-cut assessment of strongly anisotropic scaling in the present and other practical contexts. △ Less

Submitted 20 December, 2012; v1 submitted 10 September, 2012; originally announced September 2012.

Journal ref: Phys. Rev. B 86, 245427 (2012)

Strong anisotropy in two-dimensional surfaces with generic scale invariance: Gaussian and related models

Authors: Edoardo Vivo, Matteo Nicoli, Rodolfo Cuerno

Abstract: Among systems that display generic scale invariance, those whose asymptotic properties are anisotropic in space (strong anisotropy, SA) have received a relatively smaller attention, specially in the context of kinetic roughening for two-dimensional surfaces. This is in contrast with their experimental ubiquity, e.g. in the context of thin film production by diverse techniques. Based on exact resul… ▽ More Among systems that display generic scale invariance, those whose asymptotic properties are anisotropic in space (strong anisotropy, SA) have received a relatively smaller attention, specially in the context of kinetic roughening for two-dimensional surfaces. This is in contrast with their experimental ubiquity, e.g. in the context of thin film production by diverse techniques. Based on exact results for integrable (linear) cases, here we formulate a SA Ansatz that, albeit equivalent to existing ones borrowed from equilibrium critical phenomena, is more naturally adapted to the type of observables that are measured in experiments on the dynamics of thin films, such as one and two-dimensional height structure factors. We test our Ansatz on a paradigmatic nonlinear stochastic equation displaying strong anisotropy like the Hwa-Kardar equation [Phys. Rev. Lett. 62, 1813 (1989)], that was initially proposed to describe the interface dynamics of running sand piles. A very important role to elucidate its SA properties is played by an accurate (Gaussian) approximation through a non-local linear equation that shares the same asymptotic properties. △ Less

Submitted 27 November, 2012; v1 submitted 10 September, 2012; originally announced September 2012.

Journal ref: Phys. Rev. E 86, 051611 (2012)

Universality of cauliflower-like fronts: from nanoscale thin films to macroscopic plants

Authors: Mario Castro, Rodolfo Cuerno, Matteo Nicoli, Luis Vázquez, Josephus G. Buijnsters

Abstract: Chemical vapor deposition (CVD) is a widely used technique to grow solid materials with accurate control of layer thickness and composition. Under mass-transport-limited conditions, the surface of thin films thus produced grows in an unstable fashion, developing a typical motif that resembles the familiar surface of a cauliflower plant. Through experiments on CVD production of amorphous hydrogenat… ▽ More Chemical vapor deposition (CVD) is a widely used technique to grow solid materials with accurate control of layer thickness and composition. Under mass-transport-limited conditions, the surface of thin films thus produced grows in an unstable fashion, developing a typical motif that resembles the familiar surface of a cauliflower plant. Through experiments on CVD production of amorphous hydrogenated carbon films leading to cauliflower-like fronts, we provide a quantitative assessment of a continuum description of CVD interface growth. As a result, we identify non-locality, non-conservation, and randomness as the main general mechanisms controlling the formation of these ubiquitous shapes. We also show that the surfaces of actual cauliflower plants and combustion fronts obey the same scaling laws, proving the validity of the theory over seven orders of magnitude in length scales. Thus, a theoretical justification is provided, that had remained elusive thus far, for the remarkable similarity between the textures of surfaces found for systems that differ widely in physical nature and typical scales. △ Less

Submitted 10 September, 2012; originally announced September 2012.

Journal ref: New. J. Phys. 14, 103039 (2012)

Kardar-Parisi-Zhang asymptotics for the two-dimensional noisy Kuramoto-Sivashinsky equation

Authors: Matteo Nicoli, Edoardo Vivo, Rodolfo Cuerno

Abstract: We study numerically the Kuramoto-Sivashinsky (KS) equation forced by external white noise in two space dimensions, that is a generic model for e.g. surface kinetic roughening in the presence of morphological instabilities. Large scale simulations using a pseudospectral numerical scheme allow us to retrieve Kardar-Parisi-Zhang (KPZ) scaling as the asymptotic state of the system, as in the 1D case.… ▽ More We study numerically the Kuramoto-Sivashinsky (KS) equation forced by external white noise in two space dimensions, that is a generic model for e.g. surface kinetic roughening in the presence of morphological instabilities. Large scale simulations using a pseudospectral numerical scheme allow us to retrieve Kardar-Parisi-Zhang (KPZ) scaling as the asymptotic state of the system, as in the 1D case. However, this is only the case for sufficiently large values of the coupling and/or system size, so that previous conclusions on non-KPZ asymptotics are demonstrated as finite size effects. Crossover effects are comparatively stronger for the 2D case than for the 1D system. △ Less

Submitted 12 March, 2012; originally announced March 2012.

Comments: 5 pages, 3 figures; supplemental material available at journal web page and/or on request

Journal ref: Phys. Rev. E 82, 045202(R) (2010)

Dynamic effects induced by renormalization in anisotropic pattern forming systems

Authors: Adrian Keller, Matteo Nicoli, Stefan Facsko, Rodolfo Cuerno

Abstract: The dynamics of patterns in large two-dimensional domains remains a challenge in non-equilibrium phenomena. Often it is addressed through mild extensions of one-dimensional equations. We show that full 2D generalizations of the latter can lead to unexpected dynamical behavior. As an example we consider the anisotropic Kuramoto-Sivashinsky equation, that is a generic model of anisotropic pattern fo… ▽ More The dynamics of patterns in large two-dimensional domains remains a challenge in non-equilibrium phenomena. Often it is addressed through mild extensions of one-dimensional equations. We show that full 2D generalizations of the latter can lead to unexpected dynamical behavior. As an example we consider the anisotropic Kuramoto-Sivashinsky equation, that is a generic model of anisotropic pattern forming systems and has been derived in different instances of thin film dynamics. A rotation of a ripple pattern by $90^{\circ}$ occurs in the system evolution when nonlinearities are strongly suppressed along one direction. This effect originates in non-linear parameter renormalization at different rates in the two system dimensions, showing a dynamical interplay between scale invariance and wavelength selection. Potential experimental realizations of this phenomenon are identified. △ Less

Submitted 12 March, 2012; originally announced March 2012.

Comments: 5 pages, 3 figures; supplemental material available at journal web page and/or on request

Journal ref: Phys. Rev. E 84, 015202(R) (2011); Phys. Rev. E 85, 029905(E) (2012)

Solid flow drives surface nanopatterning by ion-beam irradiation

Authors: M. Castro, R. Gago, L. Vázquez, J. Muñoz-García, R. Cuerno

Abstract: Ion Beam Sputtering (IBS) is known to produce surface nanopatterns over macroscopic areas on a wide range of materials. However, in spite of the technological potential of this route to nanostructuring, the physical process by which these surfaces self-organize remains poorly under- stood. We have performed detailed experiments of IBS on Si substrates that validate dynamical and morphological pred… ▽ More Ion Beam Sputtering (IBS) is known to produce surface nanopatterns over macroscopic areas on a wide range of materials. However, in spite of the technological potential of this route to nanostructuring, the physical process by which these surfaces self-organize remains poorly under- stood. We have performed detailed experiments of IBS on Si substrates that validate dynamical and morphological predictions from a hydrodynamic description of the phenomenon. Our results elucidate flow of a nanoscopically thin and highly viscous surface layer, driven by the stress created by the ion-beam, as a description of the system. This type of slow relaxation is akin to flow of macroscopic solids like glaciers or lead pipes, that is driven by defect dynamics. △ Less

Submitted 6 March, 2012; originally announced March 2012.

Comments: 12 pages, 4 figures

Hydrodynamic approach to surface pattern formation by ion beams

Authors: Mario Castro, Rodolfo Cuerno

Abstract: On the proper timescale, amorphous solids can flow. Solid flow can be observed macroscopically in glaciers or lead pipes, but it can also be artificially enhanced by creating defects. Ion Beam Sputtering (IBS) is a technique in which ions with energies in the 0.1-10 keV range impact against a solid target inducing defect creation and dynamics, and eroding its surface leading to formation of ordere… ▽ More On the proper timescale, amorphous solids can flow. Solid flow can be observed macroscopically in glaciers or lead pipes, but it can also be artificially enhanced by creating defects. Ion Beam Sputtering (IBS) is a technique in which ions with energies in the 0.1-10 keV range impact against a solid target inducing defect creation and dynamics, and eroding its surface leading to formation of ordered nanostructures. Despite its technological interest, a basic understanding of nanopattern formation processes occurring under IBS of amorphizable targets has not been clearly established, recent experiments on Si having largely questioned knowledge accumulated during the last two decades. A number of interfacial equations have been proposed in the past to describe these phenomena, typically by adding together different contributions coming from surface diffusion, ion sputtering or mass redistribution, etc. in a non-systematic way. Here, we exploit the general idea of solids flowing due to ion impacts in order to establish a general framework into which different mechanisms (such as viscous flow, stress, diffusion, or sputtering) can be incorporated, under generic physical conservation laws. As opposed to formulating phenomenological interfacial equations, this approach allows to assess systematically the relevance and interplay of different physical mechanisms influencing surface pattern formation by IBS. △ Less

Submitted 18 October, 2011; originally announced October 2011.

Comments: Applied Surface Science (in press). Some errata have been corrected with respect to the published version

Journal ref: Applied Surface Science Volume 258, Issue 9, 15 February 2012, Pages 4171--4178

Dynamical Renormalization Group Study for a Class of Non-local Interface Equations

Authors: Matteo Nicoli, Rodolfo Cuerno, Mario Castro

Abstract: We provide a detailed Dynamic Renormalization Group study for a class of stochastic equations that describe non-conserved interface growth mediated by non-local interactions. We consider explicitly both the morphologically stable case, and the less studied case in which pattern formation occurs, for which flat surfaces are linearly unstable to periodic perturbations. We show that the latter leads… ▽ More We provide a detailed Dynamic Renormalization Group study for a class of stochastic equations that describe non-conserved interface growth mediated by non-local interactions. We consider explicitly both the morphologically stable case, and the less studied case in which pattern formation occurs, for which flat surfaces are linearly unstable to periodic perturbations. We show that the latter leads to non-trivial scaling behavior in an appropriate parameter range when combined with the Kardar-Parisi-Zhang (KPZ) non-linearity, that nevertheless does not correspond to the KPZ universality class. This novel asymptotic behavior is characterized by two scaling laws that fix the critical exponents to dimension-independent values, that agree with previous reports from numerical simulations and experimental systems. We show that the precise form of the linear stabilizing terms does not modify the hydrodynamic behavior of these equations. One of the scaling laws, usually associated with Galilean invariance, is shown to derive from a vertex cancellation that occurs (at least to one loop order) for any choice of linear terms in the equation of motion and is independent on the morphological stability of the surface, hence generalizing this well-known property of the KPZ equation. Moreover, the argument carries over to other systems like the Lai-Das Sarma-Villain equation, in which vertex cancellation is known {\em not to} imply an associated symmetry of the equation. △ Less

Submitted 3 October, 2011; originally announced October 2011.

Comments: 34 pages, 9 figures. Journal of Statistical Mechanics: Theory and Experiments (in press)

Journal ref: J. Stat. Mech. Theor. Exp. P10030 (2011)

Intrinsic geometry approach to surface kinetic roughening

Authors: Javier Rodriguez-Laguna, Silvia N. Santalla, Rodolfo Cuerno

Abstract: A model for kinetic roughening of one-dimensional interfaces is presented within an intrinsic geometry framework that is free from the standard small-slope and no-overhang approximations. The model is meant to probe the consequences of the latter on the Kardar-Parisi-Zhang (KPZ) description of non-conserved, irreversible growth. Thus, growth always occurs along the local normal direction to the in… ▽ More A model for kinetic roughening of one-dimensional interfaces is presented within an intrinsic geometry framework that is free from the standard small-slope and no-overhang approximations. The model is meant to probe the consequences of the latter on the Kardar-Parisi-Zhang (KPZ) description of non-conserved, irreversible growth. Thus, growth always occurs along the local normal direction to the interface, with a rate that is subject to fluctuations and depends on the local curvature. Adaptive numerical techniques have been designed that are specially suited to the study of fractal morphologies and can support interfaces with large slopes and overhangs. Interface self-intersections are detected, and the ensuing cavities removed. After appropriate generalization of observables such as the global and local surface roughness functions, the interface scaling is seen in our simulations to be of the Family-Vicsek type for arbitrary curvature dependence of the growth rate, KPZ scaling appearing for large sytems sizes and sufficiently large noise amplitudes. △ Less

Submitted 1 June, 2011; v1 submitted 9 May, 2011; originally announced May 2011.

Journal ref: J. Stat. Mech. (2011) P05032

One-dimensional pattern of Au nanodots by ion-beam-sputtering: formation and mechanism

Authors: J-H Kim, N-B Ha, J-S Kim, M Joe, K-R Lee, R. Cuerno

Abstract: Highly ordered one-dimensional arrays of nanodots, or nanobeads, are fabricated by forming nanoripples and nanodots in sequence, entirely by ion-beamsputtering (IBS) of Au(001). This demonstrates the capability of IBS for the fabrication of sophisticated nanostructures via hierarchical self-assembly. The intricate nanobead pattern ideally serves to identify the governing mechanisms for the pattern… ▽ More Highly ordered one-dimensional arrays of nanodots, or nanobeads, are fabricated by forming nanoripples and nanodots in sequence, entirely by ion-beamsputtering (IBS) of Au(001). This demonstrates the capability of IBS for the fabrication of sophisticated nanostructures via hierarchical self-assembly. The intricate nanobead pattern ideally serves to identify the governing mechanisms for the pattern formation: Non-linear effects, especially local redeposition and surface-confined transport, are essential both for the formation and the preservation of the pattern order. △ Less

Submitted 10 February, 2011; originally announced February 2011.

Ion induced solid flow

Authors: Mario Castro, Rodolfo Cuerno

Abstract: Amorphous solids can flow over very long periods of time. Solid flow can also be artificially enhanced by creating defects, as by Ion Beam Sputtering (IBS) in which collimated ions with energies in the 0.1 to 10 keV range impact a solid target, eroding its surface and inducing formation of nanometric structures. Recent experiments have challenged knowledge accumulated during the last two decades s… ▽ More Amorphous solids can flow over very long periods of time. Solid flow can also be artificially enhanced by creating defects, as by Ion Beam Sputtering (IBS) in which collimated ions with energies in the 0.1 to 10 keV range impact a solid target, eroding its surface and inducing formation of nanometric structures. Recent experiments have challenged knowledge accumulated during the last two decades so that a basic understanding of self-organized nano-pattern formation under IBS is still lacking. We show that considering the irradiated solid to flow like a highly viscous liquids can account for the complex IBS morphological phase diagram, relegating erosion to a subsidiary role and demonstrating a controllable instance of solid flow at the nanoscale. This new perspective can allow for a full harnessing of this bottom-up route to nanostructuring. △ Less

Submitted 11 July, 2010; originally announced July 2010.

Comments: 17 pages, 5 figures

Kinetic roughening in a realistic model of non-conserved interface growth

Authors: Matteo Nicoli, Mario Castro, Rodolfo Cuerno

Abstract: We provide a quantitative picture of non-conserved interface growth from a diffusive field making special emphasis on two main issues, the range of validity of the effective small-slopes (interfacial) theories and the interplay between the emergence of morphologically instabilities in the aggregate dynamics, and its kinetic roughening properties. Taking for definiteness electrochemical depositio… ▽ More We provide a quantitative picture of non-conserved interface growth from a diffusive field making special emphasis on two main issues, the range of validity of the effective small-slopes (interfacial) theories and the interplay between the emergence of morphologically instabilities in the aggregate dynamics, and its kinetic roughening properties. Taking for definiteness electrochemical deposition as our experimental field of reference, our theoretical approach makes use of two complementary approaches: interfacial effective equations and a phase-field formulation of the electrodeposition process. Both descriptions allow us to establish a close quantitative connection between theory and experiments. Moreover, we are able to correlate the anomalous scaling properties seen in some experiments with the failure of the small slope approximation, and to assess the effective re-emergence of standard kinetic roughening properties at very long times under appropriate experimental conditions. △ Less

Submitted 22 December, 2008; originally announced December 2008.

Comments: Journal of Statistical Mechanics: Theory & Experiment, in press

Unified moving boundary model with fluctuations for unstable diffusive growth

Authors: Matteo Nicoli, Mario Castro, Rodolfo Cuerno

Abstract: We study a moving boundary model of non-conserved interface growth that implements the interplay between diffusive matter transport and aggregation kinetics at the interface. Conspicuous examples are found in thin film production by chemical vapor deposition and electrochemical deposition. The model also incorporates noise terms that account for fluctuations in the diffusive and in the attachmen… ▽ More We study a moving boundary model of non-conserved interface growth that implements the interplay between diffusive matter transport and aggregation kinetics at the interface. Conspicuous examples are found in thin film production by chemical vapor deposition and electrochemical deposition. The model also incorporates noise terms that account for fluctuations in the diffusive and in the attachment processes. A small slope approximation allows us to derive effective interface evolution equations (IEE) in which parameters are related to those of the full moving boundary problem. In particular, the form of the linear dispersion relation of the IEE changes drastically for slow or for instantaneous attachment kinetics. In the former case the IEE takes the form of the well-known (noisy) Kuramoto-Sivashinsky equation, showing a morphological instability at short times that evolves into kinetic roughening of the Kardar-Parisi-Zhang class. In the instantaneous kinetics limit, the IEE combines Mullins-Sekerka linear dispersion relation with a KPZ nonlinearity, and we provide a numerical study of the ensuing dynamics. In all cases, the long preasymptotic transients can account for the experimental difficulties to observe KPZ scaling. We also compare our results with relevant data from experiments and discrete models. △ Less

Submitted 22 December, 2008; originally announced December 2008.

Journal ref: Phys. Rev. E 78, 021601 (2008)

Coupling of morphology to surface transport in ion-beam irradiated surfaces. I. Oblique incidence

Authors: Javier Muñoz-García, Rodolfo Cuerno, Mario Castro

Abstract: We propose and study a continuum model for the dynamics of amorphizable surfaces undergoing ion-beam sputtering (IBS) at intermediate energies and oblique incidence. After considering the current limitations of more standard descriptions in which a single evolution equation is posed for the surface height, we overcome (some of) them by explicitly formulating the dynamics of the species that tran… ▽ More We propose and study a continuum model for the dynamics of amorphizable surfaces undergoing ion-beam sputtering (IBS) at intermediate energies and oblique incidence. After considering the current limitations of more standard descriptions in which a single evolution equation is posed for the surface height, we overcome (some of) them by explicitly formulating the dynamics of the species that transport along the surface, and by coupling it to that of the surface height proper. In this we follow recent proposals inspired by ``hydrodynamic'' descriptions of pattern formation in aeolian sand dunes and ion-sputtered systems. From this enlarged model, and by exploiting the time-scale separation among various dynamical processes in the system, we derive a single height equation in which coefficients can be related to experimental parameters. This equation generalizes those obtained by previous continuum models and is able to account for many experimental features of pattern formation by IBS at oblique incidence, such as the evolution of the irradiation-induced amorphous layer, transverse ripple motion with non-uniform velocity, ripple coarsening, onset of kinetic roughening and other. Additionally, the dynamics of the full two-field model is compared with that of the effective interface equation. △ Less

Submitted 26 May, 2008; originally announced May 2008.

Comments: 23 pages, 14 figures. Movies of figures 6, 7, and 8 available at http://gisc.uc3m.es/~javier/Movies/

Self-organized surface nanopatterning by ion beam sputtering

Authors: Javier Muñoz-García, Luis Vázquez, Rodolfo Cuerno, José A. Sánchez-García, Mario Castro, Raúl Gago

Abstract: The production of nanopatterns on the surfaces of targets irradiated by ion beams at low and intermediate energies has developed during the present decade to a salient degree of control over the main pattern features. However, there is still a wide experimental scatter in the type and relevance of various dynamic interfacial properties induced by this technique as a function of system type and p… ▽ More The production of nanopatterns on the surfaces of targets irradiated by ion beams at low and intermediate energies has developed during the present decade to a salient degree of control over the main pattern features. However, there is still a wide experimental scatter in the type and relevance of various dynamic interfacial properties induced by this technique as a function of system type and parameters. In parallel, diverse theoretical models exist that differ in their capabilities to reproduce such a wide range of experimental features. We provide an overview of the most recent studies of nanoripple and dot production by ion-beam sputtering, with special attention to the comparison between experiments and (continuum) models, and with a focus on those issues that remain open or, at least, ambiguous. These are perhaps more evident for the case of nanodot patterns on amorphizable targets than for the case of nanoripples due, possibly, to the relative novelty of the former. The pattern properties to be considered are those of potential increased technological importance, such as the variation of nanodot size with parameters like ion energy, target temperature and sputtering time (i.e., fluence), as well as the conditions under which the quality of the pattern order is enhanced, issues such as wavelength coarsening and order enhancement becoming relevant. △ Less

Submitted 18 June, 2007; originally announced June 2007.

Comments: To appear in Lecture Notes on Nanoscale Science and Technology, edited by Z. Wang (Springer, Heidelberg)

Generic equations for pattern formation in evolving interfaces

Authors: M. Castro, J. Muñoz-García, R. Cuerno, M. García Hernández, L. Vázquez

Abstract: We present a general formalism which allows us to derive the evolution equations describing one-dimensional (1D) and isotropic 2D interfacelike systems, that is based on symmetries, conservation laws, multiple scale arguments, and exploits the relevance of coarsening dynamics. Our approach becomes especially significant in the presence of surface morphological instabilities and allows us to clas… ▽ More We present a general formalism which allows us to derive the evolution equations describing one-dimensional (1D) and isotropic 2D interfacelike systems, that is based on symmetries, conservation laws, multiple scale arguments, and exploits the relevance of coarsening dynamics. Our approach becomes especially significant in the presence of surface morphological instabilities and allows us to classify the most relevant nonlinear terms in the continuum description of these systems. The formalism applies to systems ranging from eroded nanostructures to macroscopic pattern formation. In particular, we show the validity of the theory for novel experiments on ion plasma erosion. △ Less

Submitted 10 May, 2007; originally announced May 2007.

Journal ref: New Journal of Physics. Vol 9. 102. (2007)

Short-range stationary patterns and long-range disorder in an evolution equation for one-dimensional interfaces

Authors: Javier Muñoz-García, Rodolfo Cuerno, Mario Castro

Abstract: A novel local evolution equation for one-dimensional interfaces is derived in the context of erosion by ion beam sputtering. We present numerical simulations of this equation which show interrupted coarsening in which an ordered cell pattern develops with constant wavelength and amplitude at intermediate distances, while the profile is disordered and rough at larger distances. Moreover, for a wi… ▽ More A novel local evolution equation for one-dimensional interfaces is derived in the context of erosion by ion beam sputtering. We present numerical simulations of this equation which show interrupted coarsening in which an ordered cell pattern develops with constant wavelength and amplitude at intermediate distances, while the profile is disordered and rough at larger distances. Moreover, for a wide range of parameters the lateral extent of ordered domains ranges up to tens of cells. This behavior is new in the context of dynamics of surfaces or interfaces with morphological instabilities. We also provide analytical estimates for the stationary pattern wavelength and mean growth velocity. △ Less

Submitted 23 January, 2007; originally announced January 2007.

Journal ref: Physical Review E, Vol. 74. 050103(R) (2006)

Single tensionless transition in the Laplacian roughening model

Authors: Juan Jesus Ruiz-Lorenzo, Esteban Moro, Rodolfo Cuerno

Abstract: We report large scale Monte Carlo simulations of the equilibrium discrete Laplacian roughening (dLr) model, originally introduced as the simplest one accommodating the hexatic phase in two-dimensional melting. The dLr model is also relevant to surface roughening in molecular beam epitaxy (MBE). Our data suggest a single phase transition, possibly of the Kosterlitz-Thouless type, between a flat l… ▽ More We report large scale Monte Carlo simulations of the equilibrium discrete Laplacian roughening (dLr) model, originally introduced as the simplest one accommodating the hexatic phase in two-dimensional melting. The dLr model is also relevant to surface roughening in molecular beam epitaxy (MBE). Our data suggest a single phase transition, possibly of the Kosterlitz-Thouless type, between a flat low-temperature phase and a rough, tensionless, high-temperature phase. Thus, earlier conclusions on the order of the phase transition and on the existence of a hexatic phase are seen as due to finite size effects, the phase diagram of the dLr model being similar to that of a continuum analog previously formulated in the context of surface growth by MBE. △ Less

Submitted 23 January, 2006; originally announced January 2006.

Comments: 4 pages, 4 figures. Accepted as a Rapid Communication in Phys. Rev. E

Journal ref: Phys. Rev. E 73, 015103 (2006)

Non-linear ripple dynamics on amorphous surfaces patterned by ion-beam sputtering

Authors: Javier Muñoz-García, Mario Castro, Rodolfo Cuerno

Abstract: Erosion by ion-beam sputtering (IBS) of amorphous targets at off-normal incidence frequently produces a (nanometric) rippled surface pattern, strongly resembling macroscopic ripples on aeolian sand dunes. Suitable generalization of continuum descriptions of the latter allows us to describe theoretically for the first time the main nonlinear features of ripple dynamics by IBS, namely, wavelength… ▽ More Erosion by ion-beam sputtering (IBS) of amorphous targets at off-normal incidence frequently produces a (nanometric) rippled surface pattern, strongly resembling macroscopic ripples on aeolian sand dunes. Suitable generalization of continuum descriptions of the latter allows us to describe theoretically for the first time the main nonlinear features of ripple dynamics by IBS, namely, wavelength coarsening and non-uniform propagation velocity, that agree with similar results in experiments and discrete models. These properties are seen to be the anisotropic counterparts of in-plane ordering and (interrupted) pattern coarsening in IBS experiments on rotating substrates and at normal incidence. △ Less

Submitted 17 May, 2006; v1 submitted 18 June, 2005; originally announced June 2005.

Comments: 5 pages, 3 figures

Journal ref: Phys. Rev. Lett. 96, 086101 (2006)

Self-Organized Ordering of Nanostructures Produced by Ion-Beam Sputtering

Authors: Mario Castro, Rodolfo Cuerno, Luiz Vazquez, Raul Gago

Abstract: We study the self-organized ordering of nanostructures produced by ion-beam sputtering (IBS) of targets amorphizing under irradiation. By introducing a model akin to models of pattern formation in aeolian sand dunes, we extend consistently the current continuum theory of erosion by IBS. We obtain new non-linear effects responsible for the in-plane ordering of the structures, whose strength corre… ▽ More We study the self-organized ordering of nanostructures produced by ion-beam sputtering (IBS) of targets amorphizing under irradiation. By introducing a model akin to models of pattern formation in aeolian sand dunes, we extend consistently the current continuum theory of erosion by IBS. We obtain new non-linear effects responsible for the in-plane ordering of the structures, whose strength correlates with the degree of ordering found in experiments. Our results highlight the importance of redeposition and surface viscous flow to this nanopattern formation process. △ Less

Submitted 17 June, 2005; originally announced June 2005.

Comments: 4 pages, 2 figures

Journal ref: Phys. Rev. Lett. 94, 016102 (2005)

Influence of Collision Cascade Statistics on Pattern Formation of Ion-Sputtered Surfaces

Authors: M. Feix, A. K. Hartmann, R. Kree, J. Munoz-Garcia, R. Cuerno

Abstract: Theoretical continuum models that describe the formation of patterns on surfaces of targets undergoing ion-beam sputtering, are based on Sigmund's formula, which describes the spatial distribution of the energy deposited by the ion. For small angles of incidence and amorphous or polycrystalline materials, this description seems to be suitable, and leads to the classic BH morphological theory [R.… ▽ More Theoretical continuum models that describe the formation of patterns on surfaces of targets undergoing ion-beam sputtering, are based on Sigmund's formula, which describes the spatial distribution of the energy deposited by the ion. For small angles of incidence and amorphous or polycrystalline materials, this description seems to be suitable, and leads to the classic BH morphological theory [R.M. Bradley and J.M.E. Harper, J. Vac. Sci. Technol. A 6, 2390 (1988)]. Here we study the sputtering of Cu crystals by means of numerical simulations under the binary-collision approximation. We observe significant deviations from Sigmund's energy distribution. In particular, the distribution that best fits our simulations has a minimum near the position where the ion penetrates the surface, and the decay of energy deposition with distance to ion trajectory is exponential rather than Gaussian. We provide a modified continuum theory which takes these effects into account and explores the implications of the modified energy distribution for the surface morphology. In marked contrast with BH's theory, the dependence of the sputtering yield with the angle of incidence is non-monotonous, with a maximum for non-grazing incidence angles. △ Less

Submitted 9 July, 2004; originally announced July 2004.

Comments: 12 pages, 13 figures, RevTex

Journal ref: Phys. Rev. B 71, 125407 (2005)

Universality issues in surface kinetic roughening of thin solid films

Authors: Rodolfo Cuerno, Luis Vazquez

Abstract: Since publication of the main contributions on the theory of kinetic roughening more than fifteen years ago, many works have been reported on surface growth or erosion that employ the framework of dynamic scaling. This interest was mainly due to the predicted existence of just a few universality classes to describe the statistical properties of the morphology of growing surfaces and interfaces t… ▽ More Since publication of the main contributions on the theory of kinetic roughening more than fifteen years ago, many works have been reported on surface growth or erosion that employ the framework of dynamic scaling. This interest was mainly due to the predicted existence of just a few universality classes to describe the statistical properties of the morphology of growing surfaces and interfaces that appear in a wide range of physical systems. Nowadays, this prediction seems to be inaccurate. This situation has caused a clear detriment of these studies in spite of the undeniable existence of kinetic roughening in many different real systems, and without a clear understanding of the reasons behind the mismatch between theoretical expectations and experimental observations. In this chapter we aim to explore existing problems and shortcomings of both the theoretical and experimental approaches, focusing mainly on growth of thin solid films. Our analysis suggests that the theoretical framework as yet is not complete, while more systematic and consistent experiments need to be performed. Once these issues are taken into account, a more consistent and useful theory of kinetic roughening might develop. △ Less

Submitted 25 February, 2004; originally announced February 2004.

Comments: Review article to appear in ``Advances in Condensed Matter and Statistical Mechanics", ed. E. Korutcheva and R. Cuerno. To be published by Nova Science Publishers. 22 pages. 4 eps figures

A microscopic model for thin film spreading

Authors: Douglas B. Abraham, Rodolfo Cuerno, Esteban Moro

Abstract: A microscopic, driven lattice gas model is proposed for the dynamics and spatio-temporal fluctuations of the precursor film observed in spreading experiments. Matter is transported both by holes and particles, and the distribution of each can be described by driven diffusion with a moving boundary. This picture leads to a stochastic partial differential equation for the shape of the boundary, wh… ▽ More A microscopic, driven lattice gas model is proposed for the dynamics and spatio-temporal fluctuations of the precursor film observed in spreading experiments. Matter is transported both by holes and particles, and the distribution of each can be described by driven diffusion with a moving boundary. This picture leads to a stochastic partial differential equation for the shape of the boundary, which agrees with the simulations of the lattice gas. Preliminary results for flow in a thermal gradient are discussed. △ Less

Submitted 14 January, 2002; originally announced January 2002.

Comments: 4 pages, 3 figures. Submitted

Journal ref: Physical Review Letters 88, 206101 (2002)

Transients due to instabilities hinder Kardar-Parisi-Zhang scaling: a unified derivation for surface growth by electrochemical and chemical vapor deposition

Authors: Rodolfo Cuerno, Mario Castro

Abstract: We propose a unified moving boundary problem for surface growth by electrochemical and chemical vapor deposition, which is derived from constitutive equations into which stochastic forces are incorporated. We compute the coefficients in the interface equation of motion as functions of phenomenological parameters. The equation features the Kardar-Parisi-Zhang (KPZ) non-linearity and instabilities… ▽ More We propose a unified moving boundary problem for surface growth by electrochemical and chemical vapor deposition, which is derived from constitutive equations into which stochastic forces are incorporated. We compute the coefficients in the interface equation of motion as functions of phenomenological parameters. The equation features the Kardar-Parisi-Zhang (KPZ) non-linearity and instabilities which, depending on surface kinetics, can hinder the asymptotic KPZ scaling. Our results account for the universality and the experimental scarcity of KPZ scaling in the growth processes considered. △ Less

Submitted 9 November, 2001; v1 submitted 29 October, 2001; originally announced October 2001.

Comments: 4 double-column pages, Revtex4, no figures; to appear in Physical Review Letters

Journal ref: Physical Review Letters 87, 236103 (2001)

Dynamic renormalization group study of a generalized continuum model of crystalline surfaces

Authors: Rodolfo Cuerno, Esteban Moro

Abstract: We apply the Nozieres-Gallet dynamic renormalization group (RG) scheme to a continuum equilibrium model of a d-dimensional surface relaxing by linear surface tension and linear surface diffusion, and which is subject to a lattice potential favoring discrete values of the height variable. The model thus interpolates between the overdamped sine-Gordon model and a related continuum model of crystal… ▽ More We apply the Nozieres-Gallet dynamic renormalization group (RG) scheme to a continuum equilibrium model of a d-dimensional surface relaxing by linear surface tension and linear surface diffusion, and which is subject to a lattice potential favoring discrete values of the height variable. The model thus interpolates between the overdamped sine-Gordon model and a related continuum model of crystalline tensionless surfaces. The RG flow predicts the existence of an equilibrium roughening transition only for d = 2 dimensional surfaces, between a flat low-temperature phase and a rough high-temperature phase in the Edwards-Wilkinson (EW) universality class. The surface is always in the flat phase for any other substrate dimensions d > 2. For any value of d, the linear surface diffusion mechanism is an irrelevant perturbation of the linear surface tension mechanism, but may induce long crossovers within which the scaling properties of the linear molecular-beam epitaxy equation are observed, thus increasing the value of the sine-Gordon roughening temperature. This phenomenon originates in the non-linear lattice potential, and is seen to occur even in the absence of a bare surface tension term. An important consequence of this is that a crystalline tensionless surface is asymptotically described at high temperatures by the EW universality class. △ Less

Submitted 25 October, 2001; originally announced October 2001.

Comments: 22 pages, 5 figures. Accepted for publication in Physical Review E

Journal ref: Physical Review E 65, 016110 (2002)