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.