Abstract To understand the distribution of the Yang–Lee zeros in quantum integrable field theories we analyse the simplest of these systems given by the 2D Yang–Lee model. The grand-canonical partition function of this quantum field theory, as a function of the fugacity z and the inverse temperature β , can be computed in terms of the thermodynamics Bethe Ansatz based on its exact S -matrix. We extract the Yang–Lee zeros in the complex plane by using a sequence of polynomials of increasing order N in z which converges to the grand-canonical partition function. We show that these zeros are distributed along curves which are approximate circles as it is also the case of the zeros for purely free theories. There is though an important difference between the interactive theory and the free theories, for the radius of the zeros in the interactive theory goes continuously to zero in the high-temperature limit β → 0 while in the free theories it remains close to 1 even for small values of β , jumping to 0 only at β = 0 .
Abstract To understand the distribution of the Yang–Lee zeros in quantum integrable field theories we analyse the simplest of these systems given by the 2D Yang–Lee model. The grand-canonical partition function of this quantum field theory, as a function of the fugacity z and the inverse temperature β , can be computed in terms of the thermodynamics Bethe Ansatz based on its exact S -matrix. We extract the Yang–Lee zeros in the complex plane by using a sequence of polynomials of increasing order N in z which converges to the grand-canonical partition function. We show that these zeros are distributed along curves which are approximate circles as it is also the case of the zeros for purely free theories. There is though an important difference between the interactive theory and the free theories, for the radius of the zeros in the interactive theory goes continuously to zero in the high-temperature limit β → 0 while in the free theories it remains close to 1 even for small values of β , jumping to 0 only at β = 0 .
The distribution of Yang-Lee zeros in the ferromagnetic Ising model in both two and three dimensions is studied on the complex field plane directly in the thermodynamic limit via the tensor network methods. The partition function is represented as a contraction of a tensor network and is efficiently evaluated with an iterative tensor renormalization scheme. The free-energy density and the magnetization are computed on the complex field plane. Via the discontinuity of the magnetization, the density of the Yang-Lee zeros is obtained to lie on the unit circle, consistent with the Lee-Yang circle theorem. Distinct features are observed at different temperatures---below, above, and at the critical temperature. Application of the tensor network approach is also made to the q-state Potts models in both two and three dimensions and a previous debate on whether, in the thermodynamic limit, the Yang-Lee zeros lie on a unit circle for q>2 is resolved: they clearly do not lie on a unit circle except at the zero temperature. For the Potts models (q=3,4,5,6) investigated in two dimensions, as the temperature is lowered the radius of the zeros at a fixed angle from the real axis shrinks exponentially towards unity with the inverse temperature.
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.
We revisit the somewhat less studied problem of Yang-Lee zeros of the Ising antiferromagnet. For this purpose, we study two models, the nearest-neighbor model on a square lattice and the more tractable mean-field model corresponding to infinite-ranged coupling between all sites. In the high-temperature limit, we show that the logarithm of the Yang-Lee zeros can be written as a series in half odd integer powers of the inverse temperature, k, with the leading term ∼k^{1/2}. This result is true in any dimension and for arbitrary lattices. We also show that the coefficients of the expansion satisfy simple identities (akin to sum rules) for the nearest-neighbor case. These identities are verified numerically by computing the exact partition function for a two-dimensional square lattice of size 16×16. For the mean-field model, we write down the partition function (termed the mean-field polynomials) for the ferromagnetic (FM) and antiferromagnetic (AFM) cases and derive from them the mean-field equations. We analytically show that at high temperatures the zeros of the AFM mean-field polynomial scale as ∼k^{1/2} as well. Using a simple numerical method, we find the roots lie on certain curves (the root curves), in the thermodynamic limit for the mean-field polynomials for the AFM case as well as for the FM one. Our results show a new root curve that was not found earlier. Our results also clearly illustrate the phase transition expected for the FM and AFM cases, in the language of Yang-Lee zeros. Moreover, for the AFM case, we observe that the root curves separate two distinct phases of zero and nonzero complex staggered magnetization, and thus depict a complex phase boundary.
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
Lee-Yang zeros at O(3) and deconfined quantum critical points
Lee-Yang theory, based on the study of zeros of the partition function, is widely regarded as a powerful and complimentary approach to the study of critical phenomena and forms a foundational part of the theory of phase transitions. Its widespread use, however, is complicated by the fact that it requires introducing complex-valued fields that create an obstacle for many numerical methods, especially in the quantum case where very limited studies exist beyond one dimension. Here we present a simple and statistically exact method to compute partition function zeros with general complex-valued external fields in the context of large-scale quantum Monte Carlo simulations. We demonstrate the power of this approach by extracting critical exponents from the leading Lee-Yang zeros of 2D quantum antiferromagnets with a complex staggered field, focusing on the Heisenberg bilayer and square-lattice J-Q models. The method also allows us to introduce a complex field that couples to valence bond solid order, where we observe extended rings of zeros in the J-Q model with purely imaginary staggered and valence bond solid fields.
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.