Non-equilibrium dynamics of the Ising model for T less-than/equal-toTc

Aging following a zero-temperature quench in the d=3 Ising model

Aging in phase-ordering kinetics of the d=3 Ising model following a quench from infinite to zero temperature is studied by means of Monte Carlo simulations. In this model the two-time spin-spin autocorrelator C_{ag} is expected to obey dynamical scaling and to follow asymptotically a power-law decay with the autocorrelation exponent λ. Previous work indicated that the lower Fisher-Huse bound of λ≥d/2=1.5 is violated in this model. Using much larger systems than previously studied, the instantaneous exponent for λ we obtain at late times does not disagree with this bound. By conducting systematic fits to the data of C_{ag} using different Ansätze for the leading correction term, we find λ=1.58(14), with most of the error attributed to the systematic uncertainty regarding the Ansätze. This result is in contrast to the recent report that below the roughening transition universality might be violated.

Automated characterization of spatial and dynamical heterogeneity in supercooled liquids via implementation of machine learning

A computational approach by an implementation of the principle component analysis (PCA) withK-means and Gaussian mixture (GM) clustering methods from machine learning algorithms to identify structural and dynamical heterogeneities of supercooled liquids is developed. In this method, a collection of the average weighted coordination numbers (WCNs‾) of particles calculated from particles' positions are used as an order parameter to build a low-dimensional representation of feature (structural) space forK-means clustering to sort the particles in the system into few meso-states using PCA. Nano-domains or aggregated clusters are also formed in configurational (real) space from a direct mapping using associated meso-states' particle identities with some misclassified interfacial particles. These classification uncertainties can be improved by a co-learning strategy which utilizes the probabilistic GM clustering and the information transfer between the structural space and configurational space iteratively until convergence. A final classification of meso-states in structural space and domains in configurational space are stable over long times and measured to have dynamical heterogeneities. Armed with such a classification protocol, various studies over the thermodynamic and dynamical properties of these domains indicate that the observed heterogeneity is the result of liquid-liquid phase separation after quenching to a supercooled state.

Phase-ordering kinetics in the Allen-Cahn (Model A) class: Universal aspects elucidated by electrically induced transition in liquid crystals

The two-dimensional (2D) Ising model is the statistical physics textbook example for phase transitions and their kinetics. Quenched through the Curie point with Glauber rates, the late-time description of the ferromagnetic domain coarsening finds its place at the scalar sector of the Allen-Cahn (or Model A) class, which encompasses phase-ordering kinetics endowed with a nonconserved order parameter. Resisting exact results sought for theoreticians since Lifshitz's first account in 1962, the central quantities of 2D Model A-most scaling exponents and correlation functions-remain known up to approximate theories whose disparate outcomes urge experimental assessment. Here we perform such assessment based on a comprehensive study of the coarsening of 2D twisted nematic liquid crystals whose kinetics is induced by a superfast electrical switching from a spatiotemporally chaotic (disordered) state to a two-phase concurrent, equilibrium one. Tracking the dynamics via optical microscopy, we first show the sharp evidence of well-established Model A aspects, such as the dynamic exponent z=2 and the dynamic scaling hypothesis, to then move forward. We confirm the Bray-Humayun theory for Porod's regime describing intradomain length scales of the two-point spatial correlators and show that their nontrivial decay beyond the Porod's scale can be captured in a free-from-parameter fashion by Gaussian theories, namely the Ohta-Jasnow-Kawasaki (OJK) and Mazenko theories. Regarding time-related statistics, we corroborate the aging hypothesis in Model A systems, which includes the collapse of two-time correlators into a master curve whose format is, actually, best accounted for by a solution of the local scaling invariance theory: the same solution that fits the 2D nonconserved Ising model correlator along with the Fisher-Huse conjecture. We also suggest the true value for the local persistence exponent in Model A class, in disfavor of the exact outcome for the diffusion and OJK equations. Finally, we observe a fractal morphology for persistence clusters and extract their universal dimension. Given its accuracy and possibilities, this experimental setup may work as a prototype to address further universality issues in the realm of nonequilibrium systems.

High-degeneracy Potts coarsening

I examine the fate of a kinetic Potts ferromagnet with a high ground-state degeneracy that undergoes a deep quench to zero temperature. I consider single spin-flip dynamics on triangular lattices of linear dimension 8≤L≤128 and set the number of spin states q equal to the number of lattice sites L×L. The ground state is the most abundant final state, and is reached with probability ≈0.71. Three-hexagon states occur with probability ≈0.26, and hexagonal tessellations with more than three clusters form with probabilities of O(10^{-3}) or less. Spanning stripe states-where the domain walls run along one of the three lattice directions-appear with probability ≈0.03. "Blinker" configurations, which contain perpetually flippable spins, also emerge, but with a probability that is vanishingly small with the system size.

Aging exponents for nonequilibrium dynamics following quenches from critical points

Via Monte Carlo simulations we study nonequilibrium dynamics in the nearest-neighbor Ising model, following quenches to points inside the ordered region of the phase diagram. With the broad objective of quantifying the nonequilibrium universality classes corresponding to spatially correlated and uncorrelated initial configurations, in this paper we present results for the decay of the order-parameter autocorrelation function for quenches from the critical point. This autocorrelation is an important probe for the aging dynamics in far-from-equilibrium systems and typically exhibits power-law scaling. From the state-of-the-art analysis of the simulation results, we quantify the corresponding exponents (λ) for both conserved and nonconserved (order-parameter) dynamics of the model in space dimension d=3. Via structural analysis we demonstrate that the exponents satisfy a bound. We also revisit the d=2 case to obtain more accurate results. It appears that irrespective of the dimension, λ is approximately the same for both conserved and nonconserved dynamics.

Initial correlation dependence of aging in phase separating solid binary mixtures and ordering ferromagnets

Following quenches of initial configurations having long range spatial correlations, prepared at the demixing critical point, to points inside the miscibility gap, we study aging phenomena in solid binary mixtures. Results on the decay of the two-time order-parameter autocorrelation functions, obtained from Monte Carlo simulations of the two-dimensional Ising model, with Kawasaki exchange kinetics, are analyzed via state-of-the art methods. The outcome is compared with that obtained for the ordering in uniaxial ferromagnets. For the latter, we have performed Monte Carlo simulations of the same model using the Glauber mechanism. For both types of systems we provide comparative discussion of our results with reference to those concerning quenches with configurations having no spatial correlation. We also discuss the role of structure on the decay of these correlations.

Growth kinetics of the two-dimensional Ising model with finite cooling rates

We investigate, via spin-flip kinetic Monte Carlo simulations, the entire growth kinetic process of the two-dimensional kinetic Ising model cooled to zero temperature with finite cooling rates. A new slow dynamics, suppressed, and unseen in the instantaneous quench emerges and eventually crosses over to the asymptotic standard coarsening behavior. We present a numerical observation of the excess defect number density that provides direct information on the generated defects due to the critical slowing down. We find that the excess defect density reveals the dynamics of defect generation, which is shown to be in good agreement with the scaling theory pioneered by Kibble and Zurek (KZ). The dynamic spin correlation function reveals that the impulse regime alluded by KZ is characterized by a unique critical coarsening process with the growth law dictated by the Kibble-Zurek mechanism. We determine a new dynamic exponent governing the growth kinetics at the onset times of the zero temperature. The proposed scaling scheme for the excess defect density leads to an analytic expression for this dynamic exponent involving the KZ exponents, indicating the extended influence of the KZ mechanism even down to the kinetics at the onset times of zero temperature. We also perform dynamic simulations of critical heating with finite rates from zero temperature where the power law relaxation (with respect to the inverse cooling rates) of the magnetization can be explained clearly by the KZ exponent.

Direct test of the Gaussian auxiliary field ansatz in nonconserved order parameter phase ordering dynamics

The assumption that the local order parameter is related to an underlying spatially smooth auxiliary field, u(r[over ⃗],t), is a common feature in theoretical approaches to non-conserved order parameter phase separation dynamics. In particular, the ansatz that u(r[over ⃗],t) is a Gaussian random field leads to predictions for the decay of the autocorrelation function which are consistent with observations, but distinct from predictions using alternative theoretical approaches. In this paper, the auxiliary field is obtained directly from simulations of the time-dependent Ginzburg-Landau equation in two and three dimensions. The results show that u(r[over ⃗],t) is equivalent to the distance to the nearest interface. In two dimensions, the probability distribution, P(u), is well approximated as Gaussian except for small values of u/L(t), where L(t) is the characteristic length-scale of the patterns. The behavior of P(u) in three dimensions is more complicated; the non-Gaussian region for small u/L(t) is much larger than that in two dimensions but the tails of P(u) begin to approach a Gaussian form at intermediate times. However, at later times, the tails of the probability distribution appear to decay faster than a Gaussian distribution.

Role of initial state and final quench temperature on aging properties in phase-ordering kinetics

We study numerically the two-dimensional Ising model with nonconserved dynamics quenched from an initial equilibrium state at the temperature T_{i}≥T_{c} to a final temperature T_{f} below the critical one. By considering processes initiating both from a disordered state at infinite temperature T_{i}=∞ and from the critical configurations at T_{i}=T_{c} and spanning the range of final temperatures T_{f}∈[0,T_{c}[ we elucidate the role played by T_{i} and T_{f} on the aging properties and, in particular, on the behavior of the autocorrelation C and of the integrated response function χ. Our results show that for any choice of T_{f}, while the autocorrelation function exponent λ_{C} takes a markedly different value for T_{i}=∞ [λ_{C}(T_{i}=∞)≃5/4] or T_{i}=T_{c} [λ_{C}(T_{i}=T_{c})≃1/8] the response function exponents are unchanged. Supported by the outcome of the analytical solution of the solvable spherical model we interpret this fact as due to the different contributions provided to autocorrelation and response by the large-scale properties of the system. As changing T_{f} is considered, although this is expected to play no role in the large-scale and long-time properties of the system, we show important effects on the quantitative behavior of χ. In particular, data for quenches to T_{f}=0 are consistent with a value of the response function exponent λ_{χ}=1/2λ_{C}(T_{i}=∞)=5/8 different from the one [λ_{χ}∈(0.5-0.56)] found in a wealth of previous numerical determinations in quenches to finite final temperatures. This is interpreted as due to important preasymptotic corrections associated to T_{f}>0.

Fractality in persistence decay and domain growth during ferromagnetic ordering: Dependence upon initial correlation

The dynamics of ordering in the Ising model, following quench to zero temperature, has been studied via Glauber spin-flip Monte Carlo simulations in space dimensions d=2 and 3. One of the primary objectives has been to understand phenomena associated with the persistent spins, viz., time decay in the number of unaffected spins, growth of the corresponding pattern, and its fractal dimensionality for varying correlation length in the initial configurations, prepared at different temperatures, at and above the critical value. It is observed that the fractal dimensionality and the exponent describing the power-law decay of persistence probability are strongly dependent upon the relative values of nonequilibrium domain size and the initial equilibrium correlation length. Via appropriate scaling analyses, these quantities have been estimated for quenches from infinite and critical temperatures. The above-mentioned dependence is observed to be less pronounced in the higher dimension. In addition to these findings for the local persistence, we present results for the global persistence as well. Furthermore, important observations on the standard domain growth problem are reported. For the latter, a controversy in d=3, related to the value of the exponent for the power-law growth of the average domain size with time, has been resolved.

Short-time critical dynamics of damage spreading in the two-dimensional Ising model

The short-time critical dynamics of propagation of damage in the Ising ferromagnet in two dimensions is studied by means of Monte Carlo simulations. Starting with equilibrium configurations at T=∞ and magnetization M=0 , an initial damage is created by flipping a small amount of spins in one of the two replicas studied. In this way, the initial damage is proportional to the initial magnetization M0 in one of the configurations upon quenching the system at T C, the Onsager critical temperature of the ferromagnetic-paramagnetic transition. It is found that, at short times, the damage increases with an exponent θ D=1.915(3) , which is much larger than the exponent θ=0.197 characteristic of the initial increase of the magnetization M(t). Also, an epidemic study was performed. It is found that the average distance from the origin of the epidemic (R2(t)) grows with an exponent z∗ ≈ η ≈ 1.9, which is the same, within error bars, as the exponent θ D. However, the survival probability of the epidemics reaches a plateau so that δ=0. On the other hand, by quenching the system to lower temperatures one observes the critical spreading of the damage at T D ≃ 0.51TC, where all the measured observables exhibit power laws with exponents θ D=1.026(3), δ=0.133(1), and z∗=1.74(3).

Curvature-driven coarsening in the two-dimensional Potts model

We study the geometric properties of polymixtures after a sudden quench in temperature. We mimic these systems with the q -states Potts model on a square lattice with and without weak quenched disorder, and their evolution with Monte Carlo simulations with nonconserved order parameter. We analyze the distribution of hull-enclosed areas for different initial conditions and compare our results with recent exact and numerical findings for q=2 (Ising) case. Our results demonstrate the memory of the presence or absence of long-range correlations in the initial state during the coarsening regime and exhibit superuniversality properties.

Logarithmic correction to scaling in domain-wall dynamics at Kosterlitz-Thouless phase transitions

With Monte Carlo simulations, we investigate the relaxation dynamics of domain walls at the Kosterlitz-Thouless phase transition, taking the two-dimensional XY model as an example. The dynamic scaling behavior is carefully analyzed, and a domain-wall roughening process is observed. Two-time correlation functions are calculated and aging phenomena are investigated. Inside the domain interface, a strong logarithmic correction to scaling is detected.

Nonequilibrium critical relaxation of the order parameter and energy in the two-dimensional ferromagnetic Potts model

The static and dynamic critical properties of the ferromagnetic q -state Potts models on a square lattice with q=2 and 3 are numerically studied via the nonequilibrium relaxation method. The relaxation behavior of both the order parameter and energy as well as that of the second moments are investigated, from which static and dynamic critical exponents can be obtained. We find that the static exponents thus obtained from the relaxation of the order parameter and energy together with the second moments of the order parameter exhibit a close agreement with the exact exponents, especially for the case of the q=2 (Ising) model, when care is taken in the choice of the initial states for the relaxation of the second moments. As for the case of q=3 , the estimates for the static exponents become less accurate, but still exhibit reasonable agreement with the exactly known static exponents. The dynamic critical exponent for the q=2 (Ising) model is estimated from the relaxation of the second moments of the order parameter with mixed initial conditions to give z(q=2) approximately 2.1668(19) .

Nonequilibrium critical dynamics with domain wall and surface

With Monte Carlo simulations, we investigate the relaxation dynamics with a domain wall for magnetic systems at the critical temperature. The dynamic scaling behavior is carefully analyzed, and a dynamic roughening process is observed. For comparison, similar analysis is applied to the relaxation dynamics with a free or disordered surface.

Domain growth morphology in curvature-driven two-dimensional coarsening

We study the distribution of domain areas, areas enclosed by domain boundaries ("hulls"), and perimeters for curvature-driven two-dimensional coarsening, employing a combination of exact analysis and numerical studies, for various initial conditions. We show that the number of hulls per unit area, n_{h}(A,t)dA , with enclosed area in the interval (A,A+dA) , is described, for a disordered initial condition, by the scaling function n_{h}(A,t)=2c_{h}(A+lambda_{h}t);{2} , where c_{h}=18pi sqrt[3] approximately 0.023 is a universal constant and lambda_{h} is a material parameter. For a critical initial condition, the same form is obtained, with the same lambda_{h} but with c_{h} replaced by c_{h}2 . For the distribution of domain areas, we argue that the corresponding scaling function has, for random initial conditions, the form n_{d}(A,t)=2c_{d}(lambda_{d}t);{tau'-2}(A+lambda_{d}t);{tau'} , where c_{d} and lambda_{d} are numerically very close to c_{h} and lambda_{h} , respectively, and tau'=18791 approximately 2.055 . For critical initial conditions, one replaces c_{d} by c_{d}2 and the exponent is tau=379187 approximately 2.027 . These results are extended to describe the number density of the length of hulls and domain walls surrounding connected clusters of aligned spins. These predictions are supported by extensive numerical simulations. We also study numerically the geometric properties of the boundaries and areas.

Short-time critical dynamics and aging phenomena in the two-dimensional XY model

With Monte Carlo methods, we simulate dynamic relaxation processes of the two-dimensional XY model at the Berezinskii-Kosterlitz-Thouless phase transition temperature and below, starting from both ordered and disordered initial states. The two-time correlation function A(t',t) is measured, and aging phenomena are investigated. The power-law correction in the spatial correlation length xi(t) for relaxation with an ordered initial state and the logarithmic correction for relaxation with a disordered initial state are carefully analyzed. The scaling functions of A(t',t) are then extracted.

Dynamic critical approach to self-organized criticality

A dynamic scaling ansatz for the approach to the self-organized critical (SOC) regime is proposed and tested by means of extensive simulations applied to the Bak-Sneppen model (BS), which exhibits robust SOC behavior. Considering the short-time scaling behavior of the density of sites [rho(t)] below the critical value, it is shown that (i) starting the dynamics with configurations such that rho(t=0)-->0 one observes an initial increase of the density with exponent theta=0.12(2); (ii) using initial configurations with rho(t=0)-->1, the density decays with exponent delta=0.47(2). It is also shown that the temporal autocorrelation decays with exponent Ca=0.35(2). Using these dynamically determined critical exponents and suitable scaling relationships, all known exponents of the BS model can be obtained, e.g., the dynamical exponent z=2.10(5), the mass dimension exponent D=2.42(5), and the exponent of all returns of the activity tauALL=0.39(2), in excellent agreement with values already accepted and obtained within the SOC regime.

Dynamic Monte Carlo simulations of the three-dimensional random-bond Potts model

The effect of random bonds on the phase transitions of the three-dimensional three-state Potts model is investigated with extensive dynamic Monte Carlo simulations. In the weakly disordered regime, the phase diagram is obtained with a recently suggested nonequilibrium reweighting method. The tricritical point separating the first- and second-order transitions is determined, and the critical exponents of the continuous phase transition induced by quenched randomness are estimated.

Critical behavior of an Ising system on the Sierpinski carpet: a short-time dynamics study

The short-time dynamic evolution of an Ising model embedded in an infinitely ramified fractal structure with noninteger Hausdorff dimension was studied using Monte Carlo simulations. Completely ordered and disordered spin configurations were used as initial states for the dynamic simulations. In both cases, the evolution of the physical observables follows a power-law behavior. Based on this fact, the complete set of critical exponents characteristic of a second-order phase transition was evaluated. Also, the dynamic exponent theta of the critical initial increase in magnetization, as well as the critical temperature, were computed. The exponent theta exhibits a weak dependence on the initial (small) magnetization. On the other hand, the dynamic exponent z shows a systematic decrease when the segmentation step is increased, i.e., when the system size becomes larger. Our results suggest that the effective noninteger dimension for the second-order phase transition is noticeably smaller than the Hausdorff dimension. Even when the behavior of the magnetization (in the case of the ordered initial state) and the autocorrelation (in the case of the disordered initial state) with time are very well fitted by power laws, the precision of our simulations allows us to detect the presence of a soft oscillation of the same type in both magnitudes that we attribute to the topological details of the generating cell at any scale.

Autocorrelation functions for phase separation in ternary mixtures

We present numerical and analytical results for the autocorrelation functions which characterize domain growth in ternary mixtures. The numerical results are obtained from Monte Carlo simulations of the spin-1 Blume-Emery-Griffiths model with spin-exchange kinetics. Further, we model the autocorrelation functions using an approach based on the continuous-time random walk formalism. The aging property of these functions is related to the time dependence of the domain-size distribution. Our analytical results are found to be in good agreement with the numerical data.

Critical behavior of the two-dimensional random-bond Potts model: a short-time dynamic approach

The short-time critical dynamics of the two-dimensional eight-state random-bond Potts model is investigated with large-scale Monte Carlo simulations. Dynamic relaxation starting from a disordered and an ordered state is carefully analyzed. The continuous phase transition induced by disorder is studied, and both the dynamic and static critical exponents are estimated. The static exponent beta/nu shows little dependence on the disorder amplitude r, while the dynamic exponent z and static exponent 1/nu vary with the strength of disorder.

Scaling of the linear response in simple aging systems without disorder

The time-dependent scaling of the thermoremanent and zero-field-cooled susceptibilities in ferromagnetic spin systems undergoing aging after a quench to a temperature at or below criticality is studied. A recent debate on their interpretation is resolved by showing that for systems with a short-ranged equilibrium spin-spin correlator and above their roughening temperature, the field-cooled susceptibility chi(FC)(t)-chi(0) approximately t(-A), where chi(0) is related to the equilibrium magnetization and the exponent A is related to the time-dependent scaling of the interface width between ordered domains. The same effect also dominates the scaling of the zero-field-cooled susceptibility chi(ZFC)(t,s), but does not enter into the thermoremanent susceptibility rho(TRM)(t,s). However, there may be large finite-time corrections to the scaling of rho(TRM)(t,s) which are explicitly derived and may be needed in order to extract reliable aging exponents. Consistency with the predictions of local scale invariance is confirmed in the Glauber-Ising and spherical models.

Experimental study of spinodal decomposition in a 1D conserved order parameter system

We study the zigzag instability coarsening of splay-bend walls formed in a nematic liquid crystal under external fields. The vertexes of zigzag can be considered as kinks in a one-dimensional order parameter system and the geometrical constraints associated with the necessary equal length sum of zig and zag segments impose a conserved quantity in this Cahn-Hilliard-type problem. In the late stage of coarsening, the characteristic length of the system L(t) shows a logarithmic increase in time and the dynamical scaling law holds. We then try to extract the nontrivial asymptotic scaling exponent lambda of the two-time correlation function, defined by lim(<phi(0,t)phi(0,t('))> approximately [L(t)/L(t('))](-lambda). The scaling exponents with respective time references, t(')=32 and 64 s, after quench are found to be lambda approximately 2 which is larger than the value with respective time reference t(')=0, predicted by numerical simulation.

Scaling anomalies in the coarsening dynamics of fractal viscous fingering patterns

We analyze a recent experiment of Sharon et al. (2003) on the coarsening, due to surface tension, of fractal viscous fingering patterns (FVFPs) grown in a radial Hele-Shaw cell. We argue that an unforced Hele-Shaw model, a natural model for that experiment, belongs to the same universality class as model B of phase ordering. Two series of numerical simulations with model B are performed, with the FVFPs grown in the experiment and with diffusion limited aggregates as the initial conditions. We observed Lifshitz-Slyozov scaling t(1/3) at intermediate distances and slow convergence to this scaling at small distances. Dynamic scale invariance breaks down at large distances.

Local scale invariance as dynamical space-time symmetry in phase-ordering kinetics

The scaling of the spatiotemporal response of coarsening systems is studied through simulations of the two-dimensional (2D) and 3D Ising model with Glauber dynamics. The scaling functions are consistent with the prediction of local scale invariance, thereby suggesting the extension of dynamical scaling to a space-time dynamical symmetry.

Corrections to scaling in two-dimensional dynamic XY and fully frustrated XY models

With large-scale Monte Carlo simulations, we investigate the two-dimensional dynamic XY and fully frustrated XY models. Dynamic relaxation starting from a disordered or an ordered state is carefully analyzed. It is confirmed that there is a logarithmic correction to scaling for a disordered start, but a power-law correction for an ordered start. Rather accurate values of the static exponent eta and the dynamic exponent z are estimated.

Short-time dynamics of a random Ising model with long-range interaction

Short-time critical dynamics of a random Ising model (model A) with long-range interaction decaying as r(-(d+sigma)) (where sigma is the parameter controlling the range of the interaction), is studied by the theoretic renormalization-group approach. In dimensions d<2sigma, the initial slip exponents theta(') describing the initial increase of the order parameter, and theta for the growth of the response function, which govern the short-time scaling behaviors, are calculated to the second order in sqrt[epsilon] with epsilon=2sigma-d. The crossover between the long-range interaction and the short-range interaction, which occurs at some sigma>2, is also discussed.

Phase ordering with a global conservation law: Ostwald ripening and coalescence

Globally conserved phase ordering dynamics is investigated in systems with short range correlations at t=0. A Ginzburg-Landau equation with a global conservation law is employed as the phase field model. The conditions are found under which the sharp-interface limit of this equation is reducible to the area-preserving motion by curvature. Numerical simulations show that, for both critical and off-critical quench, the equal-time pair correlation function exhibits dynamic scaling, and the characteristic coarsening length obeys l(t) approximately t(1/2). For the critical quench, our results are in excellent agreement with earlier results. For off-critical quench (Ostwald ripening) we investigate the dynamics of the size distribution function of the minority phase domains. The simulations show that, at large times, this distribution function has a self-similar form with growth exponent 1/2. The scaled distribution, however, strongly differs from the classical Wagner distribution. We attribute this difference to coalescence of domains. A theory of Ostwald ripening is developed that takes into account binary coalescence events. The theoretical scaled distribution function agrees well with that obtained in the simulations.

Normal scaling in globally conserved interface-controlled coarsening of fractal clusters

We find that globally conserved interface-controlled coarsening of diffusion-limited aggregates exhibits dynamic scale invariance (DSI) and normal scaling. This is demonstrated by a numerical solution of the Ginzburg-Landau equation with a global conservation law. The general sharp-interface limit of this equation is introduced and reduced to volume preserving motion by mean curvature. A simple example of globally conserved interface-controlled coarsening system: the sublimation/deposition dynamics of a solid and its vapor in a small closed vessel, is presented in detail. The results of the numerical simulations show that the scaled form of the correlation function has a power-law tail accommodating the fractal initial condition. The coarsening length exhibits normal dynamic scaling. A decrease of the cluster radius with time, predicted by DSI, is observed. The difference between global and local conservation is discussed.

Short-time critical dynamics of the two-dimensional random-bond Ising model

With Monte Carlo simulations we investigate the nonequilibrium critical dynamic behavior of the two-dimensional random-bond Ising model. Based on the short-time dynamic scaling form, we estimate all the static and dynamic exponents from dynamic processes starting with both disordered and ordered states. Corrections to scaling are carefully considered.

Monte Carlo simulations of short-time critical dynamics with a conserved quantity

With Monte Carlo simulations, we investigate short-time critical dynamics of the three-dimensional antiferromagnetic Ising model with a globally conserved magnetization m(s) (not the order parameter). From the power law behavior of the staggered magnetization (the order parameter), its second moment and the autocorrelation, we determine all static and dynamic critical exponents as well as the critical temperature. The universality class of m(s)=0 is the same as that without a conserved quantity, but the universality class of nonzero m(s) is different.

Corrections to scaling for the two-dimensional dynamic XY model

With large-scale Monte Carlo simulations, we confirm that for the two-dimensional XY model, there is a logarithmic correction to scaling in the dynamic relaxation starting from a completely disordered state, while only an inverse power law correction in the case of starting from an ordered state. The dynamic exponent z is z=2.04(1).

Dynamic approach to weak first-order phase transitions

A short-time dynamic approach to weak first-order phase transitions is proposed. Taking the two-dimensional Potts models as examples, from short-time behavior of nonequilibrium relaxational processes starting from high temperature and zero temperature states, pseudo-critical-points K* and K** are determined. A clear difference of the values for K* and K** distinguishes a weak first-order transition from a second-order one.

Deterministic equations of motion and phase ordering dynamics

We numerically solve microscopic deterministic equations of motion for the two-dimensional straight phi(4) theory with random initial states. Phase ordering dynamics is investigated. Dynamic scaling is found and it is dominated by a fixed point corresponding to the minimum energy of random initial states.