We present a numerical calculation of the Lee-Yang and Fisher zeros of the 2D Ising model using multipoint Padé approximants. We perform simulations for the 2D Ising model with ferromagnetic couplings both in the absence and in the presence of a magnetic field using a cluster spin-flip algorithm. We show that it is possible to extract genuine signature of Lee-Yang and Fisher zeros of the theory through the poles of magnetization and specific heat, using the multipoint Padé method. We extract the poles of magnetization using Padé approximants and compare their scaling with known results. We verify the circle theorem associated to the well known behavior of Lee-Yang zeros. We present our finite volume scaling analysis of the zeros done at T=Tc for a few lattice sizes, extracting to a good precision the (combination of) critical exponents βδ . The computation at the critical temperature is performed after the latter has been determined via the study of Fisher zeros, thus extracting both βc and the critical exponent ν . Results already exist for extracting the critical exponents for the Ising model in two and three dimensions making use of Fisher and Lee-Yang zeros. In this work, multipoint Padé is shown to be competitive with this respect and thus a powerful tool to study phase transitions. Published by the American Physical Society 2024
We present a numerical calculation of the Lee-Yang and Fisher zeros of the 2D Ising model using multipoint Padé approximants. We perform simulations for the 2D Ising model with ferromagnetic couplings both in the absence and in the presence of a magnetic field using a cluster spin-flip algorithm. We show that it is possible to extract genuine signature of Lee-Yang and Fisher zeros of the theory through the poles of magnetization and specific heat, using the multipoint Padé method. We extract the poles of magnetization using Padé approximants and compare their scaling with known results. We verify the circle theorem associated to the well known behavior of Lee-Yang zeros. We present our finite volume scaling analysis of the zeros done at T=Tc for a few lattice sizes, extracting to a good precision the (combination of) critical exponents βδ . The computation at the critical temperature is performed after the latter has been determined via the study of Fisher zeros, thus extracting both βc and the critical exponent ν . Results already exist for extracting the critical exponents for the Ising model in two and three dimensions making use of Fisher and Lee-Yang zeros. In this work, multipoint Padé is shown to be competitive with this respect and thus a powerful tool to study phase transitions. Published by the American Physical Society 2024
Abstract Since the landmark work of Lee and Yang, locating the zeros of the partition function in the complex magnetic-field plane has become a powerful method for studying phase transitions. Fisher later extended this approach to complex temperatures, and subsequent generalizations introduced other control parameters, such as the crystal field. In previous works (Moueddene et al 2024 J. Stat. Mech. 023206; 2024 Phys. Rev. E 110 064144) we applied this framework to the two- and three-dimensional Blume–Capel model–a system with a rich phase structure where a second-order critical line meets a first-order line at a tricritical point. We showed that the scaling of Lee-Yang, Fisher, and crystal-field zeros yields accurate critical exponents even for modest lattice sizes. In the present study, we extend this analysis and demonstrate that simulations need not be performed exactly at the nominal transition point to obtain reliable exponent estimates. Strikingly, small system sizes are sufficient, which not only improves methodological efficiency but also advances the broader goal of reducing the carbon footprint of large-scale computational studies.
There exists the famous circle theorem on the Yang-Lee zeros of the ferromagnetic Ising model. However, the Yang-Lee zeros of the antiferromagnetic Ising model are much less well understood than those of the ferromagnetic model. The precise distribution of the Yang-Lee zeros of the antiferromagnetic Ising model only with nearest-neighbor interaction J on LxL square lattices is determined as a function of temperature a=e(2betaJ) (J<0), and its relation to the phase transitions is investigated. In the thermodynamic limit (L-->infinity), the distribution of the Yang-Lee zeros of the antiferromagnetic Ising model cuts the positive real axis in the complex x=e(-2betaH) plane, resulting in the critical magnetic field +/-H(c)(a), where H(c)>0 below the critical temperature a(c)=square root of 2-1. The results suggest that the value of the scaling exponent y(h) is 1 along the critical line for a<a(c).
Lee-Yang theory is central to the analysis of thermal phase transitions. However, the underlying mechanism of the theory and the nature of Lee-Yang zeros in quantum many-body systems remains elusive. Here, we develop a unified framework for understanding quantum phase transitions from fidelity zeros induced by symmetry breaking. These zeros, arising from transitions between symmetry sectors, obey the Lee-Yang theorem and give rise to fidelity edges near critical points. Quantum criticality is further characterized through the finite-size scaling of fidelity zeros. As concrete examples, we analytically and numerically investigate fidelity zeros in one- and two-dimensional ferromagnetic Ising models under a complex magnetic field. Our results provide new insights into the mechanism of Lee-Yang theory and open avenues for exploring unexplored landscapes of phase transitions in quantum many-body systems.
The densities of Yang-Lee zeros for the Ising ferromagnet on the L x L square lattice are evaluated from the exact grand partition functions (L=3 approximately 16). The properties of the density of Yang-Lee zeros are discussed as a function of temperature T and system size L. The three different classes of phase transitions for the Ising ferromagnet--first-order phase transition, second-order phase transition, and Yang-Lee edge singularity--are clearly distinguished by estimating the magnetic scaling exponent yh from the densities of zeros for finite-size systems. The divergence of the density of zeros at Yang-Lee edge in high temperatures (Yang-Lee edge singularity), which has been detected only by the series expansion until now for the square-lattice Ising ferromagnet, is obtained from the finite-size data. The identification of the orders of phase transitions in small systems is also discussed using the density of Yang-Lee zeros.
We study Yang-Lee zeros as they approach the thermodynamic limit of the two-dimensional (2D) nearest-neighbor antiferromagnetic Ising model on square and triangular lattices. We employ the numerical linked cluster expansion (NLCE) equipped with exact enumeration (EE) of the partition function to estimate the Laplacian of the free energy, which is proportional to the zeros density. Using a modified NLCE, where the expansion can be carried directly on the Yang-Lee zeros of the involved clusters, we estimate the density of Yang-Lee zeros in the thermodynamic limit. NLCE gives significantly more zeros than EE in the complex field plane providing more insights on how the root curves look in the thermodynamic limit. For the square lattice at $T\ensuremath{\ll}{T}_{c}$, the results suggest that two vertical lines at $\ifmmode\pm\else\textpm\fi{}{h}_{c}(T)$ in the complex field plane (i.e., two concentric circles in the complex fugacity plane) are the thermodynamic root curves. A similar picture is found for the triangular lattice for phase transitions at large values of magnetic field while further study is needed for phase transitions at smaller values of magnetic field. The convergence of the NLCE and (EE) calculations of the partition function to the thermodynamic limit is studied in both lattices and the temperature-field phase diagram is obtained from Yang-Lee zeros using both methods. This NLCE-based approach will facilitate the study of different types of phase transitions using Yang-Lee zeros in future research.
Yang-Lee edge singularities are the branch point of free energy on the complex plane of physical parameters and have been shown to be the simplest universality class of phase transitions.However, Yang-Lee edge singularities have not been regarded as experimentally observable since they occur at complex physical parameters which are unphysical.A recent discovery regarding the relation between partition functions and probe spin coherence makes it experimentally feasible to access the complex plane of physical parameters.However, how to extract the critical point and the critical exponent of Yang-Lee edge singularities in many-body systems, which occurs only at the thermodynamic limit, has still been elusive.Here we show that the quantum coherence of a probe spin coupled to finite-size Ising-type spin systems presents universal scaling behavior near the Yang-Lee edge singularity.The finite-size scaling behavior of the quantum coherence of the probe spin predicts that one can extract the critical point and the critical exponent of the Yang-Lee edge singularity of an Ising-type spin system in the thermodynamic limit from the spin coherence measurement of the probe spin coupled to finite-size Ising-type spin systems.This finding provides a practical approach to studying the nature of Yang-Lee edge singularities of many-body systems.In 1952, Yang and Lee [1, 2] initiated the examination of phase transitions by studying the partition function zeros, termed Lee-Yang zeros, on the complex plane of a magnetic field h (or fugacity z for fluids).They found that for a system above its critical temperature T c , the partition function must be nonzero throughout some neighborhood of the real axis of the magnetic field within a gap determined by two edges at ( ) h T YL .Therefore the free energy is an analytic function of the real magnetic field.On the other hand, below the critical temperature T c , the Lee-Yang zeros will come arbitrarily to the real axis as the thermodynamic limit is taken, destroying the analyticity of the free energy in h for those real fields at which the Lee-Yang zeros accumulate.Kortman and Griffiths [3] pointed out that the edges of Lee-Yang zeros of a many-body system in the thermodynamic limit are singularity points, termed Yang-Lee edge singularities [4].Fisher [4] proposed that the Yang-Lee edge singularity could be considered as a new second order phase transition point with associated critical exponents.Renormalization group analyses [4] show that the Yang-Lee edge singularity is associated with an j i 3 theory and the crossover dimensionality of the Yang-Lee edge singularity is d c = 6.Later, the investigation of the Yang-Lee edge singularity has been extended to many other models, such as the classical nvector model [5], the quantum Heisenberg model [5], the spherical model [6], the quantum one-dimensional (1D) transverse Ising model [7], the hierarchical model [8], branched polymers [9], the directed-site animalsenumeration problem [10], Ising models on fractal lattices [11], Ising systems with correlated disorder [12], fluid models with repulsive-core interactions [13], the antiferromagnetic Ising model [14] etc.A comprehensive review of the study of Yang-Lee edge singularities can be found in [15, 16].Recently, Kibble-Zurek scaling was shown to appear at the Yang-Lee edge singularity [17], which supports that the Yang-Lee edge singularity is similar to conventional second order phase transitions.Yang-Lee edge singularity occurs in the ferromagnetic Ising models above its critical temperatures in a purely imaginary magnetic field h i (figure 1).For > ( ) h h T YL , the partition function acquires zeros, which becomes dense on the line with > ( ) Ih h T YL .The density of Lee-Yang zeros ( ) g T h , in the thermodynamic limit behaves as ~ -s ( ) | ( )| g T h h h T , YL when ( ) h h T YL from above, with σ being the universal critical