Universal dynamic scaling in three-dimensional Ising spin glasses

We use a nonequilibrium Monte Carlo simulation method and dynamical scaling to study the phase transition in three-dimensional Ising spin glasses. The transition point is repeatedly approached at finite velocity v (temperature change versus time) in Monte Carlo simulations starting at a high temperature. This approach has the advantage that the equilibrium limit does not have to be strictly reached for a scaling analysis to yield critical exponents. For the dynamic exponent we obtain z = 5.85(9) for bimodal couplings distribution and z = 6.00(10) for the Gaussian case. Assuming universal dynamic scaling, we combine the two results and obtain z = 5.93 +/- 0.07 for generic 3D Ising spin glasses.