We study how the finite-sized n-component model A with periodic boundary conditions relaxes near its bulk critical point from an initial nonequilibrium state with short-range correlations. Particular attention is paid to the universal long-time traces that the initial condition leaves. An approach based on renormalization-group improved perturbation theory in 4−ǫ space dimensions and a nonperturbative treatment of the q = 0 mode of the fluctuating orderparameter field is developed. This leads to a renormalized effective stochastic equation for this mode in the background of the other, q = 0 modes; we explicitly derive it to one-loop order, show that it takes the expected finitesize scaling form at the fixed point, and solve it numerically. Our results confirm for general n that the amplitude of the magnetization density m(t) in the linear relaxation-time regime depends on the initial magnetization in the universal fashion originally found in our large-n analysis [J. Stat. Phys. 73 (1993) 1]. The anomalous short-time power-law increase of m(t) also is recovered. For n = 1, our results are in fair agreement with recent Monte Carlo simulations by Li, Ritschel, and Zheng [J. Phys. A 27 (1994) L837] for the three-dimensional Ising model.
We study how the finite-sized n-component model A with periodic boundary conditions relaxes near its bulk critical point from an initial nonequilibrium state with short-range correlations. Particular attention is paid to the universal long-time traces that the initial condition leaves. An approach based on renormalization-group improved perturbation theory in 4−ǫ space dimensions and a nonperturbative treatment of the q = 0 mode of the fluctuating orderparameter field is developed. This leads to a renormalized effective stochastic equation for this mode in the background of the other, q = 0 modes; we explicitly derive it to one-loop order, show that it takes the expected finitesize scaling form at the fixed point, and solve it numerically. Our results confirm for general n that the amplitude of the magnetization density m(t) in the linear relaxation-time regime depends on the initial magnetization in the universal fashion originally found in our large-n analysis [J. Stat. Phys. 73 (1993) 1]. The anomalous short-time power-law increase of m(t) also is recovered. For n = 1, our results are in fair agreement with recent Monte Carlo simulations by Li, Ritschel, and Zheng [J. Phys. A 27 (1994) L837] for the three-dimensional Ising model.
We investigate, analytically near the dimension d(uc) =4 and numerically in d=3 , the nonequilibrium relaxational dynamics of the randomly diluted Ising model at criticality. Using the exact renormalization-group method to one loop, we compute the two times t, t(w) correlation function and fluctuation dissipation ratio (FDR) for any Fourier mode of the order parameter, of finite wave vector q . In the large time separation limit, the FDR is found to reach a nontrivial value X(infinity) independently of (small) q and coincide with the FDR associated to the total magnetization obtained previously. Explicit calculations in real space show that the FDR associated to the local magnetization converges, in the asymptotic limit, to this same value X(infinity). Through a Monte Carlo simulation, we compute the autocorrelation function in three dimensions, for different values of the dilution fraction p at T(c) (p) . Taking properly into account the corrections to scaling, we find, according to the renormalization-group predictions, that the autocorrelation exponent lambda(c) is independent of p . The analysis is complemented by a study of the nonequilibrium critical dynamics following a quench from a completely ordered state.
We present a method to calculate short-time nonequilibrium universal exponents within the functional-renormalization-group scheme. As an example, we consider the classical critical dynamics of the relaxational model A after a quench of the temperature of the system and calculate the initial-slip exponent which characterizes the nonequilibrium universal short-time behavior of both the order parameter and correlation functions. The value of this exponent is found to be consistent with the result of a perturbative dimensional expansion and of Monte Carlo simulations in three spatial dimensions.
We investigate the critical dynamics of the n-component relaxational models C and D, which incorporate the coupling of a nonconserved and conserved order parameter S, respectively, to the conserved energy density rho, under nonequilibrium conditions by means of the dynamical renormalization group. Detailed balance violations can be implemented isotropically by allowing for different effective temperatures for the heat baths coupling to the slow modes. In the case of model D with conserved order parameter, the energy density fluctuations can be integrated out, leaving no trace of the nonequilibrium perturbations in the asymptotic regime. For model C with scalar order parameter, in equilibrium governed by strong dynamic scaling (z(S)=z(rho)), we find no genuine nonequilibrium fixed point either. The nonequilibrium critical dynamics of model C with n=1 thus follows the behavior of other systems with nonconserved order parameter wherein detailed balance becomes effectively restored at the phase transition. For n> or =4, the energy density generally decouples from the order parameter. However, for n=2 and n=3, in the weak dynamic scaling regime (z(S)< or =z(rho)) entire lines of genuine nonequilibrium model C fixed points emerge to one-loop order, which are characterized by continuously varying static and dynamic critical exponents. Similarly, the nonequilibrium model C with spatially anisotropic noise and n<4 allows for continuously varying exponents, yet with strong dynamic scaling. Subjecting model D to anisotropic nonequilibrium perturbations leads to genuinely different critical behavior with softening only in subsectors of momentum space and correspondingly anisotropic scaling exponents. Similar to the two-temperature model B (randomly driven diffusive systems) the effective theory at criticality can be cast into an equilibrium model D dynamics, albeit incorporating long-range interactions of the uniaxial dipolar or ferroelastic type.
We investigate the short-time quantum critical dynamics in the imaginary-time relaxation processes of finite-size systems. Universal scaling behaviors exist in the imaginary-time evolution. In particular, the system undergoes a critical initial slip stage characterized by an exponent $\ensuremath{\theta}$, in which an initial power-law increase emerges in the imaginary-time correlation function when the initial state has a zero order parameter and a vanishing correlation length. Under different initial conditions, the quantum critical point and critical exponents can be determined from the universal scaling behaviors. We apply the method to the one- and two-dimensional transverse field Ising models using quantum Monte Carlo (QMC) simulations. In the one-dimensional case, we locate the quantum critical point at (h/J)c=1.00003(8) in the thermodynamic limit, and we estimate the critical initial slip exponent $\ensuremath{\theta}=0.3734(2)$ and the static exponent $\ensuremath{\beta}/\ensuremath{\nu}=0.1251(2)$ by analyzing data on chains of length L=32−−256 and 48--256, respectively. For the two-dimensional square-lattice system, the critical coupling ratio is given by 3.04451(7) in the thermodynamic limit, while the critical exponents are $\ensuremath{\theta}=0.209(4)$ and $\ensuremath{\beta}/\ensuremath{\nu}=0.518(1)$ estimated by data on systems of size L=24--64 and 32--64, respectively. Remarkably, the critical initial slip exponents obtained in both models are notably distinct from their classical counterparts due to the essential differences between classical and quantum dynamics. The short-time critical dynamics and the imaginary-time relaxation QMC approach can be readily adapted to various models.
Extensive Monte Carlo simulations are employed in order to study the dynamic critical behaviour of the one-dimensional Ising magnet, with algebraically decaying long-range interactions of the form 1/rd + σ, with σ = 0.75. The critical temperature, as well as the critical exponents, are evaluated from the power-law behaviour of suitable physical observables when the system is quenched from uncorrelated states, corresponding to infinite temperature, to the critical point. These results are compared with those obtained from the dynamic evolution of the system when it is annealed at the critical point from the ordered state. Also, the critical temperature in the infinite interaction limit is obtained by means of a finite-range scaling analysis of data measured with different truncated interaction ranges. All the estimated static critical exponents (γ/ν, β/ν, and 1/ν) are in good agreement with renormalization group (RG) results and previously reported numerical data obtained under equilibrium conditions. On the other hand, the dynamic exponent of the initial increase of the magnetization (θ) was close to RG predictions. However, the dynamic exponent z of the time correlation length is slightly different to the RG results probably due to the fact that it may depend on the specific dynamics used or because the two-loop expansion used in the RG analysis may be insufficient.
We study the short-time dynamics of systems that develop quasi-long-range order after a quench to the Kosterlitz-Thouless phase. With the working hypothesis that the universal short-time behavior, previously found in Ising-like systems, also occurs in the Kosterlitz-Thouless phase, we explore the scaling behavior of thermodynamic variables during the relaxational process following the quench. As a concrete example, we investigate the two-dimensional six-state clock model by Monte Carlo simulation. The exponents governing the magnetization, the second moment, and the autocorrelation function are calculated. From them, by means of scaling relations, estimates for the equilibrium exponents z and $\ensuremath{\eta}$ are derived. In particular, our estimates for the temperature-dependent anomalous dimension $\ensuremath{\eta}$ that governs the static correlation function are consistent with existing analytical and numerical results, and, thus, confirm our working hypothesis.