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:2508.20422 · Journal of the Physical Society of Japan

Lee–Yang-zero Ratio Method in Three-dimensional Ising Model

Tatsuya Wada, Masakiyo Kitazawa, K. Kanaya

Published Jan 20, 2026• 3 citations
arXivHTMLPDF

Abstract

By performing Monte Carlo simulations of the three-dimensional Ising model, we apply the recently proposed Lee-Yang-zero ratio (LYZR) method to determine the location of the critical point in this model. We demonstrate that the LYZR method is as powerful as the conventional Binder-cumulant method in studying the critical point, while the LYZR method has the advantage of suppressing the violation of the finite-size scaling and non-linearity near the critical point. We also achieve a precise determination of the values of the LYZRs at the critical point, which are universal numbers. In addition, we propose an alternative method that uses only a single Lee-Yang zero and show that it is also useful for the search for the critical point.

DOI
10.7566/jpsj.95.024002

Related papers

Locating Critical Points Using Ratios of Lee-Yang Zeros

Tatsuya Wada, Masakiyo Kitazawa, K. Kanaya

We propose a method to numerically determine the location of a critical point in general systems using the finite-size scaling of Lee-Yang zeros. This method makes use of the fact that the ratios of Lee-Yang zeros on various spatial volumes intersect at the critical point. While the method is similar to the Binder-cumulant analysis, it is advantageous in suppressing the finite-volume effects arising from the mixing of variables in general systems. We show that the method works successfully for numerically locating the critical point in the three-dimensional three-state Potts model with a nonzero external field.

arXiv:2410.19345•Physical Review Letters•Apr 24, 2025•9 citations

Critical and tricritical behavior of the d = 3 Blume-Capel model: Results from small-scale Monte Carlo simulations

Leïla Moueddene, Nikolaos G. Fytas, Bertrand Berche

We investigate the location of the critical and tricritical points of the three-dimensional Blume-Capel model by analyzing the behavior of the first Lee-Yang zero, the density of partition function zeros, and higher-order cumulants of the magnetization. Our analysis is conducted through Monte-Carlo simulations, intentionally using only small system sizes. We demonstrate that this approach yields excellent results for studying the critical behavior of the model. Our findings indicate that, at the tricritical point, where logarithmic corrections are anticipated, the numerical results align closely with the theoretical exponents describing these corrections. These expected values are then employed to accurately determine the coordinates of the tricritical point. At the model's critical point, the corrections correspond to those of the three-dimensional Ising model criticality, which we also use to precisely ascertain the critical temperature at zero crystal field. Additionally, we utilize more traditional thermodynamic quantities to validate the self-consistency of our analysis.

arXiv:2410.01710•Physical review. E•Dec 23, 2024•12 citations

Lee-Yang theory, high cumulants, and large-deviation statistics of the magnetization in the Ising model

Aydin Deger, Fredrik Brange, Christian Flindt

The authors determine the Lee-Yang zeros of an Ising spin lattice in one, two, and three dimensions from the thermodynamic fluctuations that are inherent to critical regions. They can thereby predict the critical behavior of the lattice from remarkably small system sizes and, furthermore, express the large-deviation statistics of the magnetization in terms of the Lee-Yang zeros. The method used by the authors can be applied to a wide range of critical systems both in theory and experiment.

arXiv:2006.15125•Physical review. B./Physical review. B•Nov 12, 2020•34 citations

Scaling behavior of the Heisenberg model in three dimensions

A. Gordillo-Guerrero, Ralph Kenna B, Juan Jesus Ruiz-Lorenzo

We report on extensive numerical simulations of the three-dimensional Heisenberg model and its analysis through finite-size scaling of Lee-Yang zeros. Besides the critical regime, we also investigate scaling in the ferromagnetic phase. We show that, in this case of broken symmetry, the corrections to scaling contain information on the Goldstone modes. We present a comprehensive Lee-Yang analysis, including the density of zeros, and confirm recent numerical estimates for critical exponents.

arXiv:1307.5226•Physical Review E•Dec 10, 2013•5 citations

Scaling and density of Lee-Yang zeros in the four-dimensional Ising model

Ralph Kenna B, Christian B. Lang

All of the information on the behavior of the four-dimensional Ising model is contained in the distribution and density of its partition function zeros. This model is believed to belong to the same universality class as the ${\mathrm{\ensuremath{\varphi}}}_{4}^{4}$ model which plays a central role in relativistic quantum field theory. Here the scaling behavior of the edge of the distribution of zeros and the asymptotic form for the density of zeros are determined. The finite-size dependency of the density of zeros---or the distance between zeros---at the infinite volume critical point is found using both analytic and numerical approaches. As with a previous analysis of the lowest lying zero, emphasis is laid on the multiplicative logarithmic corrections to mean field scaling behavior which are related to the triviality of the Ising and ${\mathrm{\ensuremath{\varphi}}}^{4}$ models in four dimensions.

arXiv:hep-lat/9311029•Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics•Jun 1, 1994•41 citations

Verification of Ising phase transitions in the three-dimensional Ashkin-Teller model using Monte Carlo simulations

G. Szukowski, G. Kamieniarz, Grzegorz Musiał

Monte Carlo simulations in the three-dimensional (3D) Ashkin-Teller model on a cubic lattice are performed in the regions of the two-parameter space diagram where Ising-type phase transitions are expected. The scaling behavior of the Binder cumulant Q and the magnetic susceptibility in the critical region are exploited. In simulations the periodic boundary conditions and the Metropolis algorithm are used. Starting from Ising critical exponents and applying the finite-size-scaling analysis of the cumulant Q with nonlinear corrections, the accurate positions of the critical couplings on the continuous phase transition lines are calculated. For these couplings the critical exponent y(h) is calculated analyzing the magnetic susceptibility in the framework of the finite-size scaling with corrections. The values of y(h) agree with the 3D Ising model value along two lines determined by the order parameter [ssigma] .

Physical Review E•Mar 21, 2008•12 citations
By performing Monte Carlo simulations of the three-dimensional Ising model, we apply the recently proposed Lee-Yang-zero ratio (LYZR) method to determine the location of the critical point in this model. We demonstrate that the LYZR method is as powerful as the conventional Binder-cumulant method in studying the critical point, while the LYZR method has the advantage of suppressing the violation of the finite-size scaling and non-linearity near the critical point. We also achieve a precise determination of the values of the LYZRs at the critical point, which are universal numbers. In addition, we propose an alternative method that uses only a single Lee-Yang zero and show that it is also useful for the search for the critical point.
We propose a method to numerically determine the location of a critical point in general systems using the finite-size scaling of Lee-Yang zeros. This method makes use of the fact that the ratios of Lee-Yang zeros on various spatial volumes intersect at the critical point. While the method is similar to the Binder-cumulant analysis, it is advantageous in suppressing the finite-volume effects arising from the mixing of variables in general systems. We show that the method works successfully for numerically locating the critical point in the three-dimensional three-state Potts model with a nonzero external field.
We investigate the location of the critical and tricritical points of the three-dimensional Blume-Capel model by analyzing the behavior of the first Lee-Yang zero, the density of partition function zeros, and higher-order cumulants of the magnetization. Our analysis is conducted through Monte-Carlo simulations, intentionally using only small system sizes. We demonstrate that this approach yields excellent results for studying the critical behavior of the model. Our findings indicate that, at the tricritical point, where logarithmic corrections are anticipated, the numerical results align closely with the theoretical exponents describing these corrections. These expected values are then employed to accurately determine the coordinates of the tricritical point. At the model's critical point, the corrections correspond to those of the three-dimensional Ising model criticality, which we also use to precisely ascertain the critical temperature at zero crystal field. Additionally, we utilize more traditional thermodynamic quantities to validate the self-consistency of our analysis.
All of the information on the behavior of the four-dimensional Ising model is contained in the distribution and density of its partition function zeros. This model is believed to belong to the same universality class as the ${\mathrm{\ensuremath{\varphi}}}_{4}^{4}$ model which plays a central role in relativistic quantum field theory. Here the scaling behavior of the edge of the distribution of zeros and the asymptotic form for the density of zeros are determined. The finite-size dependency of the density of zeros---or the distance between zeros---at the infinite volume critical point is found using both analytic and numerical approaches. As with a previous analysis of the lowest lying zero, emphasis is laid on the multiplicative logarithmic corrections to mean field scaling behavior which are related to the triviality of the Ising and ${\mathrm{\ensuremath{\varphi}}}^{4}$ models in four dimensions.
Monte Carlo simulations in the three-dimensional (3D) Ashkin-Teller model on a cubic lattice are performed in the regions of the two-parameter space diagram where Ising-type phase transitions are expected. The scaling behavior of the Binder cumulant Q and the magnetic susceptibility in the critical region are exploited. In simulations the periodic boundary conditions and the Metropolis algorithm are used. Starting from Ising critical exponents and applying the finite-size-scaling analysis of the cumulant Q with nonlinear corrections, the accurate positions of the critical couplings on the continuous phase transition lines are calculated. For these couplings the critical exponent y(h) is calculated analyzing the magnetic susceptibility in the framework of the finite-size scaling with corrections. The values of y(h) agree with the 3D Ising model value along two lines determined by the order parameter [ssigma] .