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:2403.05223 · SciPost Physics

Quantum-critical properties of the one- and two-dimensional random transverse-field Ising model from large-scale quantum Monte Carlo simulations

Calvin Krämer, Jan Alexander Koziol, Anja Langheld, Max Hörmann, Kai Phillip Schmidt

Published Aug 23, 2024• 5 citations
arXivHTMLPDF

Abstract

We study the ferromagnetic transverse-field Ising model with quenched disorder at T = 0 T = 0 in one and two dimensions by means of stochastic series expansion quantum Monte Carlo simulations using a rigorous zero-temperature scheme. Using a sample-replication method and averaged Binder ratios, we determine the critical shift and width exponents \nu_\mathrm{s} ν s and \nu_\mathrm{w} ν w as well as unbiased critical points by finite-size scaling. Further, scaling of the disorder-averaged magnetisation at the critical point is used to determine the order-parameter critical exponent \beta β and the critical exponent \nu_{\mathrm{av}} ν a v of the average correlation length. The dynamic scaling in the Griffiths phase is investigated by measuring the local susceptibility in the disordered phase and the dynamic exponent z' z ′ is extracted. By applying various finite-size scaling protocols, we provide an extensive and comprehensive comparison between the different approaches on equal footing. The emphasis on effective zero-temperature simulations resolves several inconsistencies in existing literature.

DOI
10.21468/scipostphys.17.2.061

Related papers

Quantum Monte Carlo study of the bond- and site-diluted transverse-field Ising model

Calvin Krämer, Max Hörmann, Kai Phillip Schmidt

We study the transverse-field Ising model on a square lattice with bond- and site-dilution at zero temperature by stochastic series expansion quantum Monte Carlo simulations. Tuning the transverse field $h$ and the dilution $p$, the quantum phase diagram of both models is explored. Both quantum phase diagrams show long-range order for small $h$ and small $p$. The ordered phase of each is separated from the disordered (quantum) Griffiths phase by second-order phase transitions on two critical lines touching at a multi-critical point. Using Binder ratios we locate quantum critical points with high accuracy. The order-parameter critical exponent $β$ and the average correlation-length exponent $ν_{\mathrm{av}}$ are determined along the critical lines and at the multi-critical points via finite-size scaling. We find three internally consistent sets of critical exponents and compare them with potentially connected universality classes. The quantum Griffiths phase in the vicinity of the phase transition lines is analyzed through the local susceptibility. Our results indicate that activated scaling occurs not only at the percolation transition, but also at the phase transition line for $p$ smaller than the percolation threshold $p_c$.

arXiv:2505.07627•Physical review. B./Physical review. B•Jul 8, 2025•1 citation

Renormalization group study of the two-dimensional random transverse-field Ising model

István A. Kovács, Ferenc Iglói

The infinite-disorder fixed point of the random transverse-field Ising model is expected to control the critical behavior of a large class of random quantum and stochastic systems having an order parameter with discrete symmetry. Here we study the model on the square lattice with a very efficient numerical implementation of the strong disorder renormalization group method, which makes us possible to treat finite samples of linear size up to $L=2048$. We have calculated sample dependent pseudocritical points and studied their distribution, which is found to be characterized by the same shift and width exponent: $\ensuremath{\nu}=1.24(2)$. For different types of disorder the infinite-disorder fixed point is shown to be characterized by the same set of critical exponents, for which we have obtained improved estimates: $x=0.982(15)$ and $\ensuremath{\psi}=0.48(2)$. We have also studied the scaling behavior of the magnetization in the vicinity of the critical point as well as dynamical scaling in the ordered and disordered Griffiths phases.

arXiv:1005.4740•Physical Review B•Aug 31, 2010•105 citations

Random transverse field Ising model in dimensiond> 1: scaling analysis in the disordered phase from the directed polymer model

Cécile Monthus, Thomas Garel

For the quantum Ising model with ferromagnetic random couplings J i , j > 0 and random transverse fields h i > 0 at zero temperature in finite dimensions d > 1, we consider the lowest order contributions in perturbation theory in ( J i , j / h i ) to obtain some information on the statistics of various observables in the disordered phase. We find that the two-point correlation scales as , where ξ typ is the typical correlation length, u is a random variable and ω coincides with the droplet exponent ω DP ( D = d − 1) of the directed polymer with D = ( d − 1) transverse directions. Our main conclusions are as follows. (i) Whenever ω > 0, the quantum model is governed by an infinite-disorder fixed point: there are two distinct correlation length exponents related by ν typ = (1 − ω)ν av ; the distribution of the local susceptibility χ loc presents the power-law tail P (χ loc ) ∼ 1/χ 1 + μ loc , where μ vanishes as ξ −ω av so that the averaged local susceptibility diverges in a finite neighborhood 0 < μ < 1 before criticality (Griffiths phase); the dynamical exponent z diverges near criticality as z = d /μ ∼ ξ ω av . (ii) In dimensions d ⩽ 3, any infinitesimal disorder flows toward this infinite-disorder fixed point with ω( d ) > 0 (for instance ω( d = 2) = 1/3 and ω( d = 3) ∼ 0.24). (iii) In finite dimensions d > 3, a finite disorder strength is necessary to flow toward the infinite-disorder fixed point with ω( d ) > 0 (for instance ω( d = 4) ≃ 0.19), whereas a finite-disorder fixed point remains possible for a small enough disorder strength. For the Cayley tree of effective dimension d = ∞, where ω = 0, we discuss the similarities and differences with the case of finite dimensions.

arXiv:1110.3145•Journal of Physics A Mathematical and Theoretical•Feb 16, 2012•26 citations

Finite-size scaling analysis of the two-dimensional random transverse-field Ising ferromagnet

Jiwon Choi, Seung Ki Baek

The random transverse-field Ising ferromagnet (RTFIF) is a highly disordered quantum system which contains randomness in the coupling strengths as well as in the transverse-field strengths. In one dimension, the critical properties are governed by an infinite-randomness fixed point (IRFP), and renormalization-group studies argue that the two-dimensional (2D) model is also governed by an IRFP. However, even the location of the critical point remains unsettled among quantum Monte Carlo (QMC) studies. In this work, we perform extensive QMC simulations to locate the quantum critical point and attempt a finite-size scaling analysis to observe the critical behavior. We estimate the critical field strength of the 2D RTFIF as ${\mathrm{\ensuremath{\Gamma}}}_{c}=7.52(2)$, together with critical exponents such as $\ensuremath{\beta}=1.5(3), \ensuremath{\nu}=1.6(3)$, and $z=3.3(3)$ or $\ensuremath{\psi}=0.50(3)$. We have also considered the McCoy-Wu model, which has randomness in the ferromagnetic coupling strengths but not in the transverse-field strength. Our QMC calculation shows that the critical behavior of the 2D McCoy-Wu model is closer to that of the 2D transverse-field Ising spin glass than to that of the 2D RTFIF. These numerical findings enhance our understanding of disordered 2D quantum systems.

arXiv:2310.14756•Physical review. B./Physical review. B•Oct 20, 2023•3 citations

Quantum-critical properties of the long-range transverse-field Ising model from quantum Monte Carlo simulations

Jan Alexander Koziol, Anja Langheld, Sebastian C. Kapfer +1 more

The quantum-critical properties of the transverse-field Ising model with algebraically decaying interactions are investigated by means of stochastic series expansion quantum Monte Carlo, on both the one-dimensional linear chain and the two-dimensional square lattice. We extract the critical exponents $\ensuremath{\nu}$ and $\ensuremath{\beta}$ as a function of the decay exponent of the long-range interactions. For ferromagnetic Ising interactions, we resolve the limiting regimes known from field theory, ranging from the nearest-neighbor Ising to the long-range Gaussian universality classes, as well as the intermediate regime with continuously varying critical exponents. In the long-range Gaussian regime, we treat the effect of dangerous irrelevant variables on finite-size scaling forms. For antiferromagnetic and therefore competing Ising interactions, the stochastic series expansion algorithm displays growing autocorrelation times leading to a reduced performance. Nevertheless, our results are consistent with the nearest-neighbor Ising universality for all investigated interaction ranges both on the linear chain and the square lattice.

arXiv:2103.09469•Physical review. B./Physical review. B•Jun 23, 2021•47 citations

Stochastic parameter optimization analysis of dynamical quantum critical phenomena in the long-range transverse-field Ising chain

Sora Shiratani, Synge Todo

The quantum phase transition of the one-dimensional long-range transverse-field Ising model is explored by combining the quantum Monte Carlo method and stochastic parameter optimization, specifically achieved by tuning correlation ratios so that space and imaginary time are isotropic. In our simulations, the simulator automatically determines the parameters to sample from, even without prior knowledge of the critical point and universality class. The leading-order finite-size corrections are eliminated by comparing two systems with different sizes; this procedure is also performed automatically. Varying the decay exponent of the long-range interaction, σ, we investigate σ dependence of the dynamical exponent and the other critical exponents precisely in the mean-field, nonuniversal, and two-dimensional classical Ising universality regimes. We successfully obtained numerical evidence supporting σ=7/4 as the universality boundary between the latter two.

arXiv:2305.14121•Physical review. E•Dec 3, 2024•2 citations
We study the ferromagnetic transverse-field Ising model with quenched disorder at T = 0 T = 0 in one and two dimensions by means of stochastic series expansion quantum Monte Carlo simulations using a rigorous zero-temperature scheme. Using a sample-replication method and averaged Binder ratios, we determine the critical shift and width exponents \nu_\mathrm{s} ν s and \nu_\mathrm{w} ν w as well as unbiased critical points by finite-size scaling. Further, scaling of the disorder-averaged magnetisation at the critical point is used to determine the order-parameter critical exponent \beta β and the critical exponent \nu_{\mathrm{av}} ν a v of the average correlation length. The dynamic scaling in the Griffiths phase is investigated by measuring the local susceptibility in the disordered phase and the dynamic exponent z' z ′ is extracted. By applying various finite-size scaling protocols, we provide an extensive and comprehensive comparison between the different approaches on equal footing. The emphasis on effective zero-temperature simulations resolves several inconsistencies in existing literature.
We study the transverse-field Ising model on a square lattice with bond- and site-dilution at zero temperature by stochastic series expansion quantum Monte Carlo simulations. Tuning the transverse field hhh and the dilution ppp, the quantum phase diagram of both models is explored. Both quantum phase diagrams show long-range order for small hhh and small ppp. The ordered phase of each is separated from the disordered (quantum) Griffiths phase by second-order phase transitions on two critical lines touching at a multi-critical point. Using Binder ratios we locate quantum critical points with high accuracy. The order-parameter critical exponent βββ and the average correlation-length exponent νavν_{\mathrm{av}}νav​ are determined along the critical lines and at the multi-critical points via finite-size scaling. We find three internally consistent sets of critical exponents and compare them with potentially connected universality classes. The quantum Griffiths phase in the vicinity of the phase transition lines is analyzed through the local susceptibility. Our results indicate that activated scaling occurs not only at the percolation transition, but also at the phase transition line for ppp smaller than the percolation threshold pcp_cpc​.
The infinite-disorder fixed point of the random transverse-field Ising model is expected to control the critical behavior of a large class of random quantum and stochastic systems having an order parameter with discrete symmetry. Here we study the model on the square lattice with a very efficient numerical implementation of the strong disorder renormalization group method, which makes us possible to treat finite samples of linear size up to L=2048L=2048L=2048. We have calculated sample dependent pseudocritical points and studied their distribution, which is found to be characterized by the same shift and width exponent: $\ensuremath{\nu}=1.24(2)$. For different types of disorder the infinite-disorder fixed point is shown to be characterized by the same set of critical exponents, for which we have obtained improved estimates: x=0.982(15)x=0.982(15)x=0.982(15) and $\ensuremath{\psi}=0.48(2)$. We have also studied the scaling behavior of the magnetization in the vicinity of the critical point as well as dynamical scaling in the ordered and disordered Griffiths phases.
For the quantum Ising model with ferromagnetic random couplings J i , j > 0 and random transverse fields h i > 0 at zero temperature in finite dimensions d > 1, we consider the lowest order contributions in perturbation theory in ( J i , j / h i ) to obtain some information on the statistics of various observables in the disordered phase. We find that the two-point correlation scales as , where ξ typ is the typical correlation length, u is a random variable and ω coincides with the droplet exponent ω DP ( D = d − 1) of the directed polymer with D = ( d − 1) transverse directions. Our main conclusions are as follows. (i) Whenever ω > 0, the quantum model is governed by an infinite-disorder fixed point: there are two distinct correlation length exponents related by ν typ = (1 − ω)ν av ; the distribution of the local susceptibility χ loc presents the power-law tail P (χ loc ) ∼ 1/χ 1 + μ loc , where μ vanishes as ξ −ω av so that the averaged local susceptibility diverges in a finite neighborhood 0 < μ < 1 before criticality (Griffiths phase); the dynamical exponent z diverges near criticality as z = d /μ ∼ ξ ω av . (ii) In dimensions d ⩽ 3, any infinitesimal disorder flows toward this infinite-disorder fixed point with ω( d ) > 0 (for instance ω( d = 2) = 1/3 and ω( d = 3) ∼ 0.24). (iii) In finite dimensions d > 3, a finite disorder strength is necessary to flow toward the infinite-disorder fixed point with ω( d ) > 0 (for instance ω( d = 4) ≃ 0.19), whereas a finite-disorder fixed point remains possible for a small enough disorder strength. For the Cayley tree of effective dimension d = ∞, where ω = 0, we discuss the similarities and differences with the case of finite dimensions.
The random transverse-field Ising ferromagnet (RTFIF) is a highly disordered quantum system which contains randomness in the coupling strengths as well as in the transverse-field strengths. In one dimension, the critical properties are governed by an infinite-randomness fixed point (IRFP), and renormalization-group studies argue that the two-dimensional (2D) model is also governed by an IRFP. However, even the location of the critical point remains unsettled among quantum Monte Carlo (QMC) studies. In this work, we perform extensive QMC simulations to locate the quantum critical point and attempt a finite-size scaling analysis to observe the critical behavior. We estimate the critical field strength of the 2D RTFIF as ${\mathrm{\ensuremath{\Gamma}}}_{c}=7.52(2)$, together with critical exponents such as $\ensuremath{\beta}=1.5(3), \ensuremath{\nu}=1.6(3)$, and z=3.3(3)z=3.3(3)z=3.3(3) or $\ensuremath{\psi}=0.50(3)$. We have also considered the McCoy-Wu model, which has randomness in the ferromagnetic coupling strengths but not in the transverse-field strength. Our QMC calculation shows that the critical behavior of the 2D McCoy-Wu model is closer to that of the 2D transverse-field Ising spin glass than to that of the 2D RTFIF. These numerical findings enhance our understanding of disordered 2D quantum systems.
The quantum-critical properties of the transverse-field Ising model with algebraically decaying interactions are investigated by means of stochastic series expansion quantum Monte Carlo, on both the one-dimensional linear chain and the two-dimensional square lattice. We extract the critical exponents $\ensuremath{\nu}$ and $\ensuremath{\beta}$ as a function of the decay exponent of the long-range interactions. For ferromagnetic Ising interactions, we resolve the limiting regimes known from field theory, ranging from the nearest-neighbor Ising to the long-range Gaussian universality classes, as well as the intermediate regime with continuously varying critical exponents. In the long-range Gaussian regime, we treat the effect of dangerous irrelevant variables on finite-size scaling forms. For antiferromagnetic and therefore competing Ising interactions, the stochastic series expansion algorithm displays growing autocorrelation times leading to a reduced performance. Nevertheless, our results are consistent with the nearest-neighbor Ising universality for all investigated interaction ranges both on the linear chain and the square lattice.
The quantum phase transition of the one-dimensional long-range transverse-field Ising model is explored by combining the quantum Monte Carlo method and stochastic parameter optimization, specifically achieved by tuning correlation ratios so that space and imaginary time are isotropic. In our simulations, the simulator automatically determines the parameters to sample from, even without prior knowledge of the critical point and universality class. The leading-order finite-size corrections are eliminated by comparing two systems with different sizes; this procedure is also performed automatically. Varying the decay exponent of the long-range interaction, σ, we investigate σ dependence of the dynamical exponent and the other critical exponents precisely in the mean-field, nonuniversal, and two-dimensional classical Ising universality regimes. We successfully obtained numerical evidence supporting σ=7/4 as the universality boundary between the latter two.