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:2305.14121 · Physical review. E

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

Sora Shiratani, Synge Todo

Published Dec 3, 2024• 2 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreve.110.064106

Related papers

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

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 +2 more

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.

arXiv:2403.05223•SciPost Physics•Aug 23, 2024•5 citations

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

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

Finite-temperature critical behavior of long-range quantum Ising models

Eduardo González Lazo, Markus Heyl, Marcello Dalmonte +1 more

We study the phase diagram and critical properties of quantum Ising chains with long-range ferromagnetic interactions decaying in a power-law fashion with exponent α \\alpha , in regimes of direct interest for current trapped ion experiments. Using large-scale path integral Monte Carlo simulations, we investigate both the ground-state and the nonzero-temperature regimes. We identify the phase boundary of the ferromagnetic phase and obtain accurate estimates for the ferromagnetic-paramagnetic transition temperatures. We further determine the critical exponents of the respective transitions. Our results are in agreement with existing predictions for interaction exponents α > 1 \\alpha > 1 up to small deviations in some critical exponents. We also address the elusive regime α < 1 \\alpha < 1 , where we find that the universality class of both the ground-state and nonzero-temperature transition is consistent with the mean-field limit at α = 0 \\alpha = 0 . Our work not only contributes to the understanding of the equilibrium properties of long-range interacting quantum Ising models, but can also be important for addressing fundamental dynamical aspects, such as issues concerning the open question of thermalization in such models.

arXiv:2104.15070•SciPost Physics•Oct 13, 2021•34 citations

Numerical Analysis of Quantum Phase Transitions with Dynamic Control of Anisotropy

Shinya Yasuda, Synge Todo

We present a method that optimizes the aspect ratio of a lattice during the Monte Carlo simulations and realizes the virtually isotropic lattice automatically. Based on the correlation lengths, the system linear lengths, including the imaginary-time direction for the quantum models, are tuned by the Robbins-Monro algorithm. The method enables us to compare directly the value of critical amplitude among different anisotropic models, and helps with identifying the universality class. We apply the present method to anisotropic and transverse field Ising models and two kinds of dimerized Heisenberg antiferromagnets, the staggered and columnar dimer models, and discuss the universality and the corrections to scaling.

Asia Pacific Physics Conference (APPC)•Mar 19, 2014•0 citations
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.
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.
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​.
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.
We study the phase diagram and critical properties of quantum Ising chains with long-range ferromagnetic interactions decaying in a power-law fashion with exponent α \\alpha , in regimes of direct interest for current trapped ion experiments. Using large-scale path integral Monte Carlo simulations, we investigate both the ground-state and the nonzero-temperature regimes. We identify the phase boundary of the ferromagnetic phase and obtain accurate estimates for the ferromagnetic-paramagnetic transition temperatures. We further determine the critical exponents of the respective transitions. Our results are in agreement with existing predictions for interaction exponents α > 1 \\alpha > 1 up to small deviations in some critical exponents. We also address the elusive regime α < 1 \\alpha < 1 , where we find that the universality class of both the ground-state and nonzero-temperature transition is consistent with the mean-field limit at α = 0 \\alpha = 0 . Our work not only contributes to the understanding of the equilibrium properties of long-range interacting quantum Ising models, but can also be important for addressing fundamental dynamical aspects, such as issues concerning the open question of thermalization in such models.
We present a method that optimizes the aspect ratio of a lattice during the Monte Carlo simulations and realizes the virtually isotropic lattice automatically. Based on the correlation lengths, the system linear lengths, including the imaginary-time direction for the quantum models, are tuned by the Robbins-Monro algorithm. The method enables us to compare directly the value of critical amplitude among different anisotropic models, and helps with identifying the universality class. We apply the present method to anisotropic and transverse field Ising models and two kinds of dimerized Heisenberg antiferromagnets, the staggered and columnar dimer models, and discuss the universality and the corrections to scaling.