Recently we showed that the critical nonequilibrium relaxation in the Swendsen-Wang algorithm is widely described by the stretched-exponential relaxation of physical quantities in the Ising or Heisenberg models. Here we make a similar analysis in the Berezinsky-Kosterlitz-Thouless phase transition in the two-dimensional (2D) XY model and in the first-order phase transition in the 2D q=5 Potts model and find that these phase transitions are described by the simple exponential relaxation and power-law relaxation of physical quantities, respectively. We compare the relaxation behaviors of these phase transitions with those of the second-order phase transition in the three- and four-dimensional XY models and in the 2D q-state Potts models for 2≤q≤4 and show that the species of phase transitions can be clearly characterized by the present analysis. We also compare the size dependence of relaxation behaviors of the first-order phase transition in the 2D q=5 and 6 Potts models and propose a quantitative criterion on "weakness" of the first-order phase transition.
Recently we showed that the critical nonequilibrium relaxation in the Swendsen-Wang algorithm is widely described by the stretched-exponential relaxation of physical quantities in the Ising or Heisenberg models. Here we make a similar analysis in the Berezinsky-Kosterlitz-Thouless phase transition in the two-dimensional (2D) XY model and in the first-order phase transition in the 2D q=5 Potts model and find that these phase transitions are described by the simple exponential relaxation and power-law relaxation of physical quantities, respectively. We compare the relaxation behaviors of these phase transitions with those of the second-order phase transition in the three- and four-dimensional XY models and in the 2D q-state Potts models for 2≤q≤4 and show that the species of phase transitions can be clearly characterized by the present analysis. We also compare the size dependence of relaxation behaviors of the first-order phase transition in the 2D q=5 and 6 Potts models and propose a quantitative criterion on "weakness" of the first-order phase transition.
Recently, it was shown [Y. Nonomura, J. Phys. Soc. Jpn. 83, 113001 (2014)JUPSAU0031-901510.7566/JPSJ.83.113001] that the nonequilibrium critical relaxation of the two-dimensional (2D) Ising model from a perfectly ordered state in the Wolff algorithm is described by stretched-exponential decay, and a universal scaling scheme was found to connect nonequilibrium and equilibrium behaviors. In the present study we extend these findings to vector spin models, and the 3D Heisenberg model could be a typical example. To evaluate the critical temperature and critical exponents precisely using the above scaling scheme, we calculate nonequilibrium ordering from the perfectly disordered state in the Swendsen-Wang algorithm, and we find that the critical ordering process is described by stretched-exponential growth with a comparable exponent to that of the 3D XY model. The critical exponents evaluated in the present study are consistent with those in previous studies.
Although the nonequilibrium-relaxation (NER) method has been widely used in Monte Carlo studies on phase transitions in classical spin systems, such studies have been quite limited in quantum phase transitions. The reason is that the relaxation process based on cluster-update quantum Monte Carlo (QMC) algorithms, which are now standards in Monte Carlo studies on quantum systems, has been considered "too fast" for such analyses. Recently, the present authors revealed that the NER process in classical spin systems based on cluster-update algorithms is characterized by stretched-exponential critical relaxation, rather than conventional power-law relaxation in local-update algorithms. In the present article, we show that this is also the case in quantum phase transitions analyzed with the cluster-update QMC. As the simplest example of isotropic quantum spin models that exhibit quantum phase transitions, we investigate the Néel-dimer quantum phase transition in the two-dimensional S=1/2 columnar-dimerized antiferromagnetic Heisenberg model with the continuous-time loop algorithm, and we confirm stretched-exponential critical relaxation consistent with the three-dimensional classical Heisenberg model in the Swendsen-Wang algorithm.
Recently the present authors proposed the nonequilibrium-to-equilibrium scaling (NE-ES) scheme for the critical Monte Carlo relaxation process [Nonomura, J. Phys. Soc. Jpn. 83, 113001 (2014)10.7566/JPSJ.83.113001], which scales relaxation data in the whole simulation-time regions regardless of functional forms, namely, both for the stretched-exponential critical relaxation in cluster algorithms and for the power-law critical relaxation in local-update algorithms. In the present study, we generalize this scheme to off-critical relaxation process and scale relaxation data for various temperatures in the whole simulation-time regions. This proposal of the off-critical scaling in cluster algorithms cannot be described by the dynamical finite-size scaling theory based on the power-law critical relaxation. As an example, we investigate the three-dimensional classical Heisenberg model previously analyzed with the NE-ES [Nonomura and Tomita, Phys. Rev. E 93, 012101 (2016)10.1103/PhysRevE.93.012101] in the Swendsen-Wang and Metropolis algorithms.
A simple and efficient numerical analysis is proposed for the Kosterlitz-Thouless (KT) phase transition. The nonequilibrium relaxation method is applied to it. The two-dimensional ferromagnetic XY models are investigated to show the efficiency. At the KT transition point as well as inside the KT phase, the nonequilibrium relaxation of magnetization from the all-aligned state shows an asymptotic power-law decay, m(t) approximately t(-lambda(T)). Only outside the KT phase, an asymptotic single exponential decay is observed. Using a standard scaling form m(t)=tau(-lambda)(-)m(t/tau) in this regime, where tau is the relaxation time at each temperature, we find a simple and efficient numerical estimation of the KT transition point and dynamical exponent. This method can be applied to various kinds of models which show the KT-like behavior.
Recently the present authors showed that the stretched-exponential critical relaxation is widely observed in classical spin systems in the nonequilibrium relaxation process of cluster algorithms. In the present article we find that this behavior can be generalized to quantum phase transitions analyzed with the continuous-time loop algorithm, where the cluster update is introduced aprior. As an example, we consider the Neel-dimer quantum phase transition in the two-dimensional S=1/2 columnar-dimerized antiferromagnetic Heisenberg model on a square lattice, and confirm the stretched-exponential critical relaxation consistent with the three-dimensional classical Heisenberg model in the Swendsen-Wang algorithm.
Using a nonequilibrium relaxation method, we calculate the dynamic critical exponent z of the two-dimensional Ising model for the Swendsen–Wang and Wolff algorithms. We examine dynamic relaxation processes following a quench from a disordered or an ordered initial state to the critical temperature T c , and measure the exponential relaxation time of the system energy. For the Swendsen–Wang algorithm with an ordered or a disordered initial state, and for the Wolff algorithm with an ordered initial state, the exponential relaxation time fits well to a logarithmic size dependence up to a lattice size L = 8192. For the Wolff algorithm with a disordered initial state, we obtain an effective dynamic exponent z exp = 1.19(2) up to L = 2048. For comparison, we also compute the effective dynamic exponents through the integrated correlation times. In addition, an exact result of the Swendsen–Wang dynamic spectrum of a one-dimensional Ising chain is derived.