Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0110294 · Physical review. B, Condensed matter

Self-adapting method for the localization of quantum critical points using quantum Monte Carlo techniques

Fabien Alet, Erik Schwartz Sørensen

Published Feb 11, 2002• 3 citations
arXivHTMLPDF

Abstract

A generalization to the quantum case of a recently introduced algorithm [Tomita and Okabe, Phys. Rev. Lett. 86, 572 (2001)] for the determination of the critical temperature of classical spin models is proposed. We describe a simple method to automatically locate critical points in (Quantum) Monte Carlo simulations. The algorithm assumes the existence of a finite correlation length in at least one of the two phases surrounding the quantum critical point. We illustrate these ideas on the example of the critical interchain coupling for which coupled antiferromagnetic $S=1$ spin chains order at $T=0.$ Finite-size scaling relations are used to determine the exponents, $\ensuremath{\nu}=0.72(2)$ and $\ensuremath{\eta}=0.038(3)$ in agreement with previous estimates.

DOI
10.1103/physrevb.65.092408

Related papers

New quantum Monte Carlo study of quantum critical phenomena with Trotter-number-dependent finite-size scaling and non-equilibrium relaxation

Yoshihiko Nonomura

We propose a new efficient scheme for the quantum Monte Carlo study of quantum critical phenomena in quantum spin systems. Rieger and Young's Trotter-number-dependent finite-size scaling in quantum spin systems and Ito et al 's evaluation of the transition point with the non-equilibrium relaxation in classical spin systems are combined and generalized. That is, only one Trotter number and one inverse temperature proportional to system sizes are taken for each system size, and the transition point of the transformed classical spin model is estimated as the point at which the order parameter shows power-law decay. The present scheme is confirmed by the determination of the critical phenomenon of the one-dimensional asymmetric XY model in the ground state. The estimates of the transition point and the critical exponents , and are in good agreement with the exact solutions. The dynamical critical exponent is also estimated as , which is consistent with that of the two-dimensional Ising model.

arXiv:cond-mat/9808337•Journal of Physics A Mathematical and General•Oct 2, 1998•9 citations

Universal short-time quantum critical dynamics of finite-size systems

Yu-Rong Shu, Shuai Yin, Dao‐Xin Yao

We investigate the short-time quantum critical dynamics in the imaginary-time relaxation processes of finite-size systems. Universal scaling behaviors exist in the imaginary-time evolution. In particular, the system undergoes a critical initial slip stage characterized by an exponent $\ensuremath{\theta}$, in which an initial power-law increase emerges in the imaginary-time correlation function when the initial state has a zero order parameter and a vanishing correlation length. Under different initial conditions, the quantum critical point and critical exponents can be determined from the universal scaling behaviors. We apply the method to the one- and two-dimensional transverse field Ising models using quantum Monte Carlo (QMC) simulations. In the one-dimensional case, we locate the quantum critical point at ${(h/J)}_{c}=1.000\phantom{\rule{0.16em}{0ex}}03(8)$ in the thermodynamic limit, and we estimate the critical initial slip exponent $\ensuremath{\theta}=0.3734(2)$ and the static exponent $\ensuremath{\beta}/\ensuremath{\nu}=0.1251(2)$ by analyzing data on chains of length $L=32--256$ and 48--256, respectively. For the two-dimensional square-lattice system, the critical coupling ratio is given by $3.044\phantom{\rule{0.16em}{0ex}}51(7)$ in the thermodynamic limit, while the critical exponents are $\ensuremath{\theta}=0.209(4)$ and $\ensuremath{\beta}/\ensuremath{\nu}=0.518(1)$ estimated by data on systems of size $L=24$--64 and 32--64, respectively. Remarkably, the critical initial slip exponents obtained in both models are notably distinct from their classical counterparts due to the essential differences between classical and quantum dynamics. The short-time critical dynamics and the imaginary-time relaxation QMC approach can be readily adapted to various models.

arXiv:1705.05931•Physical review. B./Physical review. B•Sep 12, 2017•17 citations

Quantum Monte Carlo simulation method for spin systems

Anders W. Sandvik, J. Kurkijärvi

A quantum Monte Carlo simulation scheme for spin systems is presented. The method is a generalization of Handscomb's method but applicable to any length of the spin, i.e., when the spin traces cannot be evaluated analytically. The Monte Carlo sampling is extended to the space of spin vectors in addition to the usual operator-index sequences. An important technical point is that the index sequences are augmented with the aid of unit operators to a constant, self-consistently determined length. The scheme is applied to the one-dimensional antiferromagnetic spin-S Heisenberg model. Results at low temperatures are reported for S=1 and S=3/2 and system sizes up to N=64. The computed magnetic structure factor in the S=1 chain is in agreement with earlier ground-state calculations. For S=3/2 we find the exponent \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=0.49\ifmmode\pm\else\textpm\fi{}0.04 for the divergence of the antiferromagnetic structure factor. Further, the susceptibility as a function of the wave number is computed. For S=1 the staggered susceptibility \ensuremath{\chi}(\ensuremath{\pi}) at T=0 is found to take the value 20.0\ifmmode\pm\else\textpm\fi{}1.5 in units such that \ensuremath{\chi}(q)\ensuremath{\rightarrow}${\mathit{T}}^{\mathrm{\ensuremath{-}}1}$ at high temperatures (with the temperature scale defined by ${\mathit{k}}_{\mathit{B}}$=1). For S=3/2 we obtain the exponent \ensuremath{\gamma}=1.45\ifmmode\pm\else\textpm\fi{}0.05 for the divergence of the staggered susceptibility.

Physical review. B, Condensed matter•Mar 1, 1991•528 citations

New Quantum Monte Carlo Approach to Ground-State Phase Transitions in Quantum Spin Systems

Yoshihiko Nonomura

A new efficient scheme for the finite-size-scaling study of ground-state phase transitions in quantum spin systems is proposed on the basis of the quantum Monte Carlo method. That is, Trotter numbers and simulated temperatures are scaled in accordance with system sizes, and only one simulation is required for each system size. The transition point is estimated independently of the critical exponents using the non-equilibrium relaxation of the order parameter in a very large system. This scheme is confirmed by the estimation of the critical phenomenon of the one-dimensional S =1/2 asymmetric XY model.

Journal of the Physical Society of Japan•Jan 15, 1998•20 citations

Critical exponents of the classical three-dimensional Heisenberg model: A single-cluster Monte Carlo study

Christian Holm, Wolfhard Janke E

We have simulated the three-dimensional Heisenberg model on simple cubic lattices, using the single-cluster Monte Carlo update algorithm. The expected pronounced reduction of critical slowing down at the phase transition is verified. This allows simulations on significantly larger lattices than in previous studies and consequently a better control over systematic errors. In one set of simulations we employ the usual finite-size scaling methods to compute the critical exponents \ensuremath{\nu},\ensuremath{\alpha},\ensuremath{\beta},\ensuremath{\gamma},\ensuremath{\eta} from a few measurements in the vicinity of the critical point, making extensive use of histogram reweighting and optimization techniques. In another set of simulations we report measurements of improved estimators for the spatial correlation length and the susceptibility in the high-temperature phase, obtained on lattices with up to ${100}^{3}$ spins. This enables us to compute independent estimates of \ensuremath{\nu} and \ensuremath{\gamma} from power-law fits of their critical divergencies.

arXiv:hep-lat/9301002•Physical review. B, Condensed matter•Jul 1, 1993•232 citations

Stochastic approximation of dynamical exponent at quantum critical point

Shinya Yasuda, Hidemaro Suwa, Synge Todo

We have developed a unified finite-size scaling method for quantum phase transitions that requires no prior knowledge of the dynamical exponent $z$. During a quantum Monte Carlo simulation, the temperature is automatically tuned by the Robbins-Monro stochastic approximation method, being proportional to the lowest gap of the finite-size system. The dynamical exponent is estimated in a straightforward way from the system-size dependence of the temperature. As a demonstration of our novel method, the two-dimensional $S=1/2$ quantum $XY$ model in uniform and staggered magnetic fields is investigated in the combination of the world-line quantum Monte Carlo worm algorithm. In the absence of a uniform magnetic field, we obtain the fully consistent result with the Lorentz invariance at the quantum critical point, $z=1$, i.e., the three-dimensional classical $XY$ universality class. Under a finite uniform magnetic field, on the other hand, the dynamical exponent becomes two, and the mean-field universality with effective dimension $(2+2)$ governs the quantum phase transition.

arXiv:1506.04837•Physical Review B•Sep 14, 2015•10 citations
A generalization to the quantum case of a recently introduced algorithm [Tomita and Okabe, Phys. Rev. Lett. 86, 572 (2001)] for the determination of the critical temperature of classical spin models is proposed. We describe a simple method to automatically locate critical points in (Quantum) Monte Carlo simulations. The algorithm assumes the existence of a finite correlation length in at least one of the two phases surrounding the quantum critical point. We illustrate these ideas on the example of the critical interchain coupling for which coupled antiferromagnetic S=1S=1S=1 spin chains order at T=0.T=0.T=0. Finite-size scaling relations are used to determine the exponents, \ensuremathν=0.72(2)\ensuremath{\nu}=0.72(2)\ensuremathν=0.72(2) and \ensuremathη=0.038(3)\ensuremath{\eta}=0.038(3)\ensuremathη=0.038(3) in agreement with previous estimates.
We propose a new efficient scheme for the quantum Monte Carlo study of quantum critical phenomena in quantum spin systems. Rieger and Young's Trotter-number-dependent finite-size scaling in quantum spin systems and Ito et al 's evaluation of the transition point with the non-equilibrium relaxation in classical spin systems are combined and generalized. That is, only one Trotter number and one inverse temperature proportional to system sizes are taken for each system size, and the transition point of the transformed classical spin model is estimated as the point at which the order parameter shows power-law decay. The present scheme is confirmed by the determination of the critical phenomenon of the one-dimensional asymmetric XY model in the ground state. The estimates of the transition point and the critical exponents , and are in good agreement with the exact solutions. The dynamical critical exponent is also estimated as , which is consistent with that of the two-dimensional Ising model.
We investigate the short-time quantum critical dynamics in the imaginary-time relaxation processes of finite-size systems. Universal scaling behaviors exist in the imaginary-time evolution. In particular, the system undergoes a critical initial slip stage characterized by an exponent $\ensuremath{\theta}$, in which an initial power-law increase emerges in the imaginary-time correlation function when the initial state has a zero order parameter and a vanishing correlation length. Under different initial conditions, the quantum critical point and critical exponents can be determined from the universal scaling behaviors. We apply the method to the one- and two-dimensional transverse field Ising models using quantum Monte Carlo (QMC) simulations. In the one-dimensional case, we locate the quantum critical point at (h/J)c=1.00003(8){(h/J)}_{c}=1.000\phantom{\rule{0.16em}{0ex}}03(8)(h/J)c​=1.00003(8) in the thermodynamic limit, and we estimate the critical initial slip exponent $\ensuremath{\theta}=0.3734(2)$ and the static exponent $\ensuremath{\beta}/\ensuremath{\nu}=0.1251(2)$ by analyzing data on chains of length L=32−−256L=32--256L=32−−256 and 48--256, respectively. For the two-dimensional square-lattice system, the critical coupling ratio is given by 3.04451(7)3.044\phantom{\rule{0.16em}{0ex}}51(7)3.04451(7) in the thermodynamic limit, while the critical exponents are $\ensuremath{\theta}=0.209(4)$ and $\ensuremath{\beta}/\ensuremath{\nu}=0.518(1)$ estimated by data on systems of size L=24L=24L=24--64 and 32--64, respectively. Remarkably, the critical initial slip exponents obtained in both models are notably distinct from their classical counterparts due to the essential differences between classical and quantum dynamics. The short-time critical dynamics and the imaginary-time relaxation QMC approach can be readily adapted to various models.
A quantum Monte Carlo simulation scheme for spin systems is presented. The method is a generalization of Handscomb's method but applicable to any length of the spin, i.e., when the spin traces cannot be evaluated analytically. The Monte Carlo sampling is extended to the space of spin vectors in addition to the usual operator-index sequences. An important technical point is that the index sequences are augmented with the aid of unit operators to a constant, self-consistently determined length. The scheme is applied to the one-dimensional antiferromagnetic spin-S Heisenberg model. Results at low temperatures are reported for S=1 and S=3/2 and system sizes up to N=64. The computed magnetic structure factor in the S=1 chain is in agreement with earlier ground-state calculations. For S=3/2 we find the exponent \ensuremath{\gamma}\ifmmode\bar\else\textasciimacron\fi{}=0.49\ifmmode\pm\else\textpm\fi{}0.04 for the divergence of the antiferromagnetic structure factor. Further, the susceptibility as a function of the wave number is computed. For S=1 the staggered susceptibility \ensuremath{\chi}(\ensuremath{\pi}) at T=0 is found to take the value 20.0\ifmmode\pm\else\textpm\fi{}1.5 in units such that \ensuremath{\chi}(q)\ensuremath{\rightarrow}${\mathit{T}}^{\mathrm{\ensuremath{-}}1}$ at high temperatures (with the temperature scale defined by kB{\mathit{k}}_{\mathit{B}}kB​=1). For S=3/2 we obtain the exponent \ensuremath{\gamma}=1.45\ifmmode\pm\else\textpm\fi{}0.05 for the divergence of the staggered susceptibility.
A new efficient scheme for the finite-size-scaling study of ground-state phase transitions in quantum spin systems is proposed on the basis of the quantum Monte Carlo method. That is, Trotter numbers and simulated temperatures are scaled in accordance with system sizes, and only one simulation is required for each system size. The transition point is estimated independently of the critical exponents using the non-equilibrium relaxation of the order parameter in a very large system. This scheme is confirmed by the estimation of the critical phenomenon of the one-dimensional S =1/2 asymmetric XY model.
We have simulated the three-dimensional Heisenberg model on simple cubic lattices, using the single-cluster Monte Carlo update algorithm. The expected pronounced reduction of critical slowing down at the phase transition is verified. This allows simulations on significantly larger lattices than in previous studies and consequently a better control over systematic errors. In one set of simulations we employ the usual finite-size scaling methods to compute the critical exponents \ensuremath{\nu},\ensuremath{\alpha},\ensuremath{\beta},\ensuremath{\gamma},\ensuremath{\eta} from a few measurements in the vicinity of the critical point, making extensive use of histogram reweighting and optimization techniques. In another set of simulations we report measurements of improved estimators for the spatial correlation length and the susceptibility in the high-temperature phase, obtained on lattices with up to 1003{100}^{3}1003 spins. This enables us to compute independent estimates of \ensuremath{\nu} and \ensuremath{\gamma} from power-law fits of their critical divergencies.
We have developed a unified finite-size scaling method for quantum phase transitions that requires no prior knowledge of the dynamical exponent zzz. During a quantum Monte Carlo simulation, the temperature is automatically tuned by the Robbins-Monro stochastic approximation method, being proportional to the lowest gap of the finite-size system. The dynamical exponent is estimated in a straightforward way from the system-size dependence of the temperature. As a demonstration of our novel method, the two-dimensional S=1/2S=1/2S=1/2 quantum XYXYXY model in uniform and staggered magnetic fields is investigated in the combination of the world-line quantum Monte Carlo worm algorithm. In the absence of a uniform magnetic field, we obtain the fully consistent result with the Lorentz invariance at the quantum critical point, z=1z=1z=1, i.e., the three-dimensional classical XYXYXY universality class. Under a finite uniform magnetic field, on the other hand, the dynamical exponent becomes two, and the mean-field universality with effective dimension (2+2)(2+2)(2+2) governs the quantum phase transition.