We employ Monte Carlo simulations to study the nonequilibrium relaxation of driven Ising lattice gases in two dimensions. Whereas the temporal scaling of the density autocorrelation function in the nonequilibrium steady state does not allow a precise measurement of the critical exponents, these can be accurately determined from the aging scaling of the two-time autocorrelations and the order parameter evolution following a quench to the critical point. We obtain excellent agreement with renormalization group predictions based on the standard Langevin representation of driven Ising lattice gases.
We employ Monte Carlo simulations to study the nonequilibrium relaxation of driven Ising lattice gases in two dimensions. Whereas the temporal scaling of the density autocorrelation function in the nonequilibrium steady state does not allow a precise measurement of the critical exponents, these can be accurately determined from the aging scaling of the two-time autocorrelations and the order parameter evolution following a quench to the critical point. We obtain excellent agreement with renormalization group predictions based on the standard Langevin representation of driven Ising lattice gases.
We numerically investigate the long-time behavior of the density-density auto-correlation function in driven lattice gases, with particle exclusion and periodic boundary conditions in one, two, and three dimensions using precise Monte Carlo simulations of larger system sizes than previous studies. 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, correlations are weak in three dimensions and consequently the mean-field description is adequate. We also investigate the relaxation towards the non-equilibrium steady state in the two-time density-density auto-correlations, starting from strongly correlated initial conditions. We obtain simple aging scaling behavior in one, two, and three dimensions, with the expected power laws. We numerically investigate the behavior of driven lattice gases with nearest neighbor interactions at half-filling with periodic boundary conditions below and at the critical temperature using Monte Carlo simulations of very large lattices in two dimensions. This work is one of few that explores the relaxation to a non-equilibrium steady state. We obtain data collapse for the finite-size scaling form of density-density auto-correlation function at the critical point. We achieve data collapse using finite-size scaling of the time-dependent order parameter during the transient regime starting from strongly correlated initial conditions. We present simple aging scaling of the density-density auto-correlation function at the critical point starting from strongly correlated initial conditions using Monte Carlo simulations of two different lattice anisotropies. We thus unambiguously confirm the critical exponents determined by renormalization group methods using measurement of dynamic quantities in the transient regime. Measuring these dynamic quantities in the transient regime provides more conclusive measurements of the critical exponents than previous studies measuring static quantities in the stationary state. We provide qualitative arguments that the lattice anisotropy determines the steady-state for sub-critical quenches.
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.
With dynamic Monte Carlo simulations,aging phenomena of two-dimensional Ising model and dynamic XY model were studied numerically.Systems were quenched to their critical point from a disordered initial state.The autocorrelation function was calculated by heat-bath algorithm.Critical exponents were also calculated.It was confirmed that for the 2D XY model,there is a logarithmic correction to scaling.The universality proposal and their scale relation in aging phenomena were predicted.The scale relation agreed with renormalization group calculations.
With a large-scale Monte Carlo simulation, the aging of the two-dimensional dynamic Ising and XY models is numerically studied. Starting from a disordered state, the systems are quenched to critical temperature Tc and then evolve according to dynamics of model A. By calculating the autocorrelation function A(t,t′) we observe the aging phenomenon for both spin systems. Then, relatively accurate values of the dynamic critical exponent λc and z are estimated. Especially, we confirm that for the two-dimensional XY model, there is a strong logarithmic correction to scaling in the dynamic relaxation starting from a completely disordered.
The crossover between quasiequilibrium and nonequilibrium spin-glass dynamics has been studied in Monte Carlo simulations of two- and three-dimensional short-range Ising spin-glass systems. The spin system was quenched in zero field to a low temperature T. After equilibrating for a time tw, a weak magnetic field was applied and the time dependence of the magnetization M(t) and the spin autocorrelation function q(t) was studied. Our results show that the relaxation rate S(t)=\ensuremath{\partial}M(t)/\ensuremath{\partial} log10(t) exhibits a maximum at t\ensuremath{\approxeq}tw. If the temperature is lowered or raised, immediately prior to the field application, the apparent age (ta) of the system, as defined by the maximum in relaxation rate, is shifted to longer or shorter time scales, respectively. No significant differences between the aging behavior in two and three dimensions are found. The relation \ensuremath{\chi}(t)=[1-q(t)]/T holds only for short time scales (t\ensuremath{\ll}tw), clear deviations are found at longer time scales (t>tw). Since this relation should according to the fluctuation-dissipation theorem hold for a system in equilibrium, this result is interpreted as an indication of a crossover between equilibrium and nonequilibrium dynamics. Our results are qualitatively in agreement with experiments as well as with theoretical models for the nonequilibrium spin-glass dynamics.
We employ Monte Carlo simulations to investigate the nonequilibrium relaxation properties of the two- and three-dimensional Coulomb glass with different long-range repulsive interactions. Specifically, we explore the aging scaling laws in the two-time density autocorrelation function. We find that, in the time window and parameter range accessible to us, the scaling exponents are not universal, depending on the filling fraction and temperature: As either the temperature decreases or the filling fraction deviates more from half filling, the exponents reflect markedly slower relaxation kinetics. In comparison with a repulsive Coulomb potential, appropriate for impurity states in strongly disordered semiconductors, we observe that, for logarithmic interactions, the soft pseudogap in the density of states is considerably broader, and the dependence of the scaling exponents on external parameters is much weaker. The latter situation is relevant for flux creep in the disorder-dominated Bose glass phase of type-II superconductors subject to columnar pinning centers.