Abstract:We use Monte Carlo simulations to demonstrate generic scaling aspects of classical phase transitions approached through a quench (or annealing) protocol where the temperature changes as a function of time with velocity $v$. Using a generalized Kibble-Zurek ansatz, we demonstrate dynamic scaling for different types of stochastic dynamics (Metropolis, Swendsen-Wang, and Wolff) on Ising models in two and higher dimensions. We show that there are dual scaling functions governing the dynamic scaling, which together describe the scaling behavior in the entire velocity range $v \in [0, \infty)$. These functions have asymptotics corresponding to the adiabatic and diabatic limit, and close to these limits they are perturbative in $v$ and $1/v$, respectively. Away from their perturbative domains, both functions cross over into the same universal power-law scaling form governed by the static and dynamic critical exponents (as well as an exponent characterizing the quench protocol). As a by-product of the scaling studies, we obtain high-precision estimates of the dynamic exponent $z$ for the two-dimensional Ising model subject to the three variants of Monte Carlo dynamics; for single-spin Metropolis updates $z_M=2.1767(5)$, for Swendsen-Wang multi-cluster updates $z_{\rm SW}=0.297(3)$, and for Wolff single-cluster updates $z_{W}=0.30(2)$. For Wolff dynamics, we find an interesting behavior with a non-analytic break-down of the quasi-adiabatic and diabatic scaling, instead of the generic smooth cross-over described by a power law. We interpret this disconnect between the two scaling regimes as a dynamic phase transition of the Wolff algorithm, caused by an effective sudden loss of ergodicity at high velocity.

Inchworm Monte Carlo for exact non-adiabatic dynamics I. Theory and algorithms

Worm-type Monte Carlo Simulation of the Ashkin-Teller Model on the Triangular Lattice

Inchworm quasi Monte Carlo for quantum impurities

Solving Caldeira-Leggett Model by Inchworm Method with Frozen Gaussian Approximation

Real-Time Simulation of Open Quantum Spin Chains with Inchworm Method

Simulation of Spin Chains with off-diagonal Coupling Using Inchworm Method

Simulation of Spin Chains with Off-Diagonal Coupling Using the Inchworm Method

Numerical and Monte Carlo Bethe Ansatz Method: 1d Heisenberg Model

Steady-state properties of multi-orbital systems using quantum Monte Carlo

Worm algorithms for classical statistical models

Nonadiabatic Dynamics Via the Symmetrical Quasi-Classical Method in the Presence of Anharmonicity.

A path integral methodology for obtaining thermodynamic properties of nonadiabatic systems using Gaussian mixture distributions

Dynamic scaling at classical phase transitions approached through non-equilibrium quenching

Monte Carlo simulations for Ising spins with spin greater than 1/2 applied to the square and triangular lattices with antiferromagnetic interactions and comparing results using Kawasaki and Glauber dynamics

Dynamical Mean Field Theory of the Bilayer Hubbard Model with Inchworm Monte Carlo

A new algorithm for Monte Carlo simulation of Ising spin systems

An Exact Monte Carlo Method for Continuum Fermion Systems

Semiclassical Dynamics in Wigner Phase Space I : Adiabatic Hybrid Wigner Dynamics

Quench dynamics in higher-dimensional Holstein models: Insights from Truncated Wigner Approaches

Diagrammatic Monte Carlo for Dissipative Quantum Impurity Models

Nonequilibrium candidate Monte Carlo: A new tool for efficient equilibrium simulation