Symmetry-based determination of space-time functions in nonequilibrium growth processes

Aging in the Long-Range Ising Model

The current understanding of aging phenomena is mainly confined to the study of systems with short-ranged interactions. Little is known about the aging of long-ranged systems. Here, the aging in the phase-ordering kinetics of the two-dimensional Ising model with power-law long-range interactions is studied via Monte Carlo simulations. The dynamical scaling of the two-time spin-spin autocorrelator is well described by simple aging for all interaction ranges studied. The autocorrelation exponents are consistent with λ=1.25 in the effectively short-range regime, while for stronger long-range interactions the data are consistent with λ=d/2=1. For very long-ranged interactions, strong finite-size effects are observed. We discuss whether such finite-size effects could be misinterpreted phenomenologically as subaging.

Dynamics of the critical Casimir force for a conserved order parameter after a critical quench

Fluctuation-induced forces occur generically when long-range correlations (e.g., in fluids) are confined by external bodies. In classical systems, such correlations require specific conditions, e.g., a medium close to a critical point. On the other hand, long-range correlations appear more commonly in certain nonequilibrium systems with conservation laws. Consequently, a variety of nonequilibrium fluctuation phenomena, including fluctuation-induced forces, have been discovered and explored recently. Here we address a long-standing problem of nonequilibrium critical Casimir forces emerging after a quench to the critical point in a confined fluid with order-parameter-conserving dynamics and non-symmetry-breaking boundary conditions. The interplay of inherent (critical) fluctuations and dynamical nonlocal effects (due to density conservation) gives rise to striking features, including correlation functions and forces exhibiting oscillatory time dependences. Complex transient regimes arise, depending on initial conditions and the geometry of the confinement. Our findings pave the way for exploring a wealth of nonequilibrium processes in critical fluids (e.g., fluctuation-mediated self-assembly or aggregation). In certain regimes, our results are applicable to active matter.

Pathways to equilibrium orientation fluctuations in finite stripe-forming systems

Small-angle orientation fluctuations in ordered stripe-forming systems free of topological defects can exhibit aging and anisotropic growth of two length scales. In infinitely extended systems, the stripe orientation field develops a dominant modulation length λ_{∥}^{*}(t) in the direction parallel to the stripes, which increases with time t as λ_{∥}^{*}(t)∼t^{1/4}. Simultaneously, the orientation correlation length ξ_{⊥}(t) in the direction perpendicular to the stripes increases as ξ_{⊥}(t)∼t^{1/2} [Riesch et al., Interface Focus 7, 20160146 (2017)2042-889810.1098/rsfs.2016.0146]. Here we show that finite systems of size L_{⊥}×L_{∥} with periodic boundary conditions reach equilibrium when the dominant modulation length λ_{∥}^{*}(t) reaches the system size L_{∥} in the stripe direction. The equilibration time τ_{eq}^{∥} is solely determined by L_{∥}, with τ_{eq}^{∥}∼L_{∥}^{4}. In systems with L_{⊥}<L_{∥}^{2}/2πλ_{p}, where λ_{p} is the undulation penetration length, the initial aging and coarsening dynamics changes at the crossover time τ_{C}^{⊥}∼L_{⊥}^{2} to an aging and coarsening dynamics described by the one-dimensional Mullins-Herring equation, before reaching equilibrium at τ_{∥}^{eq}. Our work reveals the two pathways to equilibrium in stripe phases with periodic boundary conditions, the finite-size scaling behavior of equilibrium orientation fluctuations, and the characteristic exponents associated with the influence of a finite system size.

Scaling properties of ageing orientation fluctuations in stripe phases

We investigate the non-equilibrium dynamics of an ordered stripe-forming system free of topological defects. In particular, we study the ageing and the coarsening of orientation fluctuations parallel and perpendicular to the stripes via computer simulations based on a minimal phase-field model (model B with Coulomb interactions). Under the influence of noise, the stripe orientation field develops fluctuations parallel to the stripes, with the dominant modulation length λ*∥ increasing with time t as λ*∥ ∼ t1/4 and the correlation length perpendicular to the stripes ξ⊥θ increasing as ξ⊥θ ∼ t1/2. We explain these anisotropic coarsening dynamics with an analytic theory based on the linear elastic model for stripe displacements first introduced by Landau and Peierls. We thus obtain the scaling forms and the scaling exponents characterizing the correlation functions and the structure factor of the stripe orientation field. Our results reveal how the coarsening of orientation fluctuations prevents a periodically modulated phase free of topological defects from reaching equilibrium.

Relaxation dynamics of vortex lines in disordered type-II superconductors following magnetic field and temperature quenches

We study the effects of rapid temperature and magnetic field changes on the nonequilibrium relaxation dynamics of magnetic vortex lines in disordered type-II superconductors by employing an elastic line model and performing Langevin molecular dynamics simulations. In a previously equilibrated system, either the temperature is suddenly changed or the magnetic field is instantaneously altered which is reflected in adding or removing flux lines to or from the system. The subsequent aging properties are investigated in samples with either randomly distributed pointlike or extended columnar defects, which allows us to distinguish the complex relaxation features that result from either type of pinning centers. One-time observables such as the radius of gyration and the fraction of pinned line elements are employed to characterize steady-state properties, and two-time correlation functions such as the vortex line height autocorrelations and their mean-square displacement are analyzed to study the nonlinear stochastic relaxation dynamics in the aging regime.

Phenomenology of aging in the Kardar-Parisi-Zhang equation

We study aging during surface growth processes described by the one-dimensional Kardar-Parisi-Zhang equation. Starting from a flat initial state, the systems undergo simple aging in both correlators and linear responses, and its dynamical scaling is characterized by the aging exponents a=-1/3, b=-2/3, λ(C)=λ(R)=1, and z=3/2. The form of the autoresponse scaling function is well described by the recently constructed logarithmic extension of local scale invariance.

Statistics of interfacial fluctuations of radially growing clusters

The dynamics of fluctuating radially growing interfaces is approached using the formalism of stochastic growth equations on growing domains. This framework reveals a number of dynamic features arising during surface growth. For fast growth, dilution, which spatially reorders the incoming matter, is responsible for the transmission of correlations. Its effects include the erasing of memory with respect to the initial condition, a partial attenuation of geometrically originated instabilities, and the restoration of universality in some special cases in which the critical exponents depend on the parameters of the equation of motion. In this sense, dilution rends the dynamics more similar to the usual one of planar systems. This fast growth regime is also characterized by the spatial decorrelation of the interface, which, in the case of radially growing interfaces, naturally originates rapid roughening and scale-dependent fractality, and suggests the advent of a self-similar fractal dimension. The center-of-mass fluctuations of growing clusters are also studied, and our analysis suggests the possible nonapplicability of usual scalings to the long-range surface fluctuations of the radial Eden model. In fact, our study points to the fact that this model belongs to a dilution-free universality class.

Relaxation of surface steps after thermal quenches: a numerical study within the terrace-step-kink model

We study the out-of-equilibrium relaxation of surface steps after thermal quenches using numerical simulations of the terrace-step-kink model for a vicinal surface. We analyze both single and interacting steps in a situation where the temperature is suddenly changed at a given quench time. We focus on a physically relevant range of temperatures and show that the relaxation of the roughness is compatible with a power-law behavior with an effective relaxation exponent close to γ = 1/2 in all cases. This value is consistent with a one-dimensional Edwards-Wilkinson equation. In particular, this means that, although the case of interacting steps is effectively a two-dimensional system, its relaxation is dominated by short length-scale fluctuations, where steps are not interacting.

Slow relaxation and aging kinetics for the driven lattice gas

We numerically investigate the long-time behavior of the density-density autocorrelation function in driven lattice gases with particle exclusion and periodic boundary conditions in one, two, and three dimensions using precise Monte Carlo simulations. In the one-dimensional asymmetric exclusion process on a ring with half the lattice sites occupied, we find that correlations induce extremely slow relaxation to the asymptotic power law decay. We compare the crossover functions obtained from our simulations with various analytic results in the literature and analyze the characteristic oscillations that occur in finite systems away from half-filling. As expected, in three dimensions correlations are weak and consequently the mean-field description is adequate. We also investigate the relaxation toward the nonequilibrium steady state in the two-time density-density autocorrelations, starting from strongly correlated initial conditions. We obtain simple aging scaling behavior in one, two, and three dimensions, with the expected power laws.

Dynamic scaling in thin-film growth with irreversible step-edge attachment

We study dynamic scaling in a model with collective diffusion (CD) of isolated atoms in terraces and irreversible aggregation at step edges. Simulations are performed in two-dimensional substrates with several diffusion to deposition ratios R identical with D/F. Data collapse of scaled roughness distributions confirms that this model is in the class of the fourth-order nonlinear growth equation by Villain, Lai, and Das Sarma (VLDS) with negligible finite-size effects, while estimates of scaling exponents show some discrepancies. This result is consistent with the prediction of a recent renormalization group approach and improves previous numerical works on related models. The roughness follows dynamic scaling as W=Lalpha/R1/2f(xi/L), with correlation length xi=(Rt)1/z, where z is the dynamic exponent. We also propose a limited mobility (LM) model where the incident atom executes up to G steps before a new atom is adsorbed, and irreversibly aggregates at step edges. This model is also shown to belong to the VLDS class. The size of the plateaus in the film surface increases as G1/2 and the lateral correlation scales as G1/2t1/z. The time evolution of the roughness reproduces that of the CD model if an equivalent parameter G approximately R2/z is chosen. This suggests the possibility of using LM models with tunable diffusion length to simulate processes with simultaneous diffusion of many atoms. A scaling approach is used to justify exponent values and dynamic relations for both models, including the significant decrease of surface roughness as R or G increases.

Changing growth conditions during surface growth

Motivated by a series of experiments that revealed a temperature dependence of the dynamic scaling regime of growing surfaces, we investigate theoretically how a nonequilibrium growth process reacts to a sudden change of system parameters. We discuss quenches between correlated regimes through exact expressions derived from the stochastic Edwards-Wilkinson equation with a variable diffusion constant. Our study reveals that a sudden change of the diffusion constant leads to remarkable changes in the surface roughness. Different dynamic regimes, characterized by a power-law or by an exponential relaxation, are identified, and a dynamic phase diagram is constructed. We conclude that growth processes provide one of the rare instances where quenches between correlated regimes yield a power-law relaxation.

Renormalization of stochastic lattice models: epitaxial surfaces

We present the application of a method [C. A. Haselwandter and D. D. Vvedensky, Phys. Rev. E 76, 041115 (2007)] for deriving stochastic partial differential equations from atomistic processes to the morphological evolution of epitaxial surfaces driven by the deposition of new material. Although formally identical to the one-dimensional (1D) systems considered previously, our methodology presents substantial additional technical issues when applied to two-dimensional (2D) surfaces. Once these are addressed, subsequent coarse-graining is accomplished as before by calculating renormalization-group (RG) trajectories from initial conditions determined by the regularized atomistic models. Our applications are to the Edwards-Wilkinson (EW) model [S. F. Edwards and D. R. Wilkinson, Proc. R. Soc. London, Ser. A 381, 17 (1982)], the Wolf-Villain (WV) model [D. E. Wolf and J. Villain, Europhys. Lett. 13, 389 (1990)], and a model with concurrent random deposition and surface diffusion. With our rules for the EW model no appreciable crossover is obtained for either 1D or 2D substrates. For the 1D WV model, discussed previously, our analysis reproduces the crossover sequence known from kinetic Monte Carlo (KMC) simulations, but for the 2D WV model, we find a transition from smooth to unstable growth under repeated coarse-graining. Concurrent surface diffusion does not change this behavior, but can lead to extended transient regimes with kinetic roughening. This provides an explanation of recent experiments on Ge(001) with the intriguing conclusion that the same relaxation mechanism responsible for ordered structures during the early stages of growth also produces an instability at longer times that leads to epitaxial breakdown. The RG trajectories calculated for concurrent random deposition and surface diffusion reproduce the crossover sequences observed with KMC simulations for all values of the model parameters, and asymptotically always approach the fixed point corresponding to the equation proposed by Villain [J. Phys. I 1, 19 (1991)] and by Lai and Das Sarma [Phys. Rev. Lett. 66, 2899 (1991)]. We conclude with a discussion of the application of our methodology to other growth settings.