Lee–Yang theory. In statistical mechanics, Lee–Yang theory, sometimes also known as Yang–Lee theory, is a scientific theory which seeks to describe phase transitions in large physical systems in the thermodynamic limit based on the properties of small, finite-size systems. The theory revolves around the complex zeros of partition functions of finite-size systems and how these may reveal the existence of phase transitions in the thermodynamic limit.

Lee–Yang theory constitutes an indispensable part of the theories of phase transitions. Originally developed for the Ising model, the theory has been extended and applied to a wide range of models and phenomena, including protein folding, percolation, complex networks, and molecular zippers.

The theory is named after the Nobel laureates Tsung-Dao Lee and Yang Chen-Ning, who were awarded the 1957 Nobel Prize in Physics for their unrelated work on parity non-conservation in weak interaction.

Introduction 
For an equilibrium system in the canonical ensemble, all statistical information about the system is encoded in the partition function,

where the sum runs over all possible microstates, and  is the inverse temperature,  is the Boltzmann constant and  is the energy of a microstate. The moments  of the energy statistics are obtained by differentiating the partition function with respect to the inverse temperature multiple times,

From the partition function, we may also obtain the free energy

Analogously to how the partition function generates the moments, the free energy generates the cumulants of the energy statistics

More generally, if the microstate energies  depend on a control parameter  and a fluctuating conjugate variable  (whose value may depend on the microstate), the moments of  may be obtained as

and the cumulants as

For instance, for a spin system, the control parameter may be an external magnetic field, , and the conjugate variable may be the total magnetization, .

Phase transitions and Lee–Yang theory 

The partition function and the free energy are intimately linked to phase transitions, for which there is a sudden change in the properties of a physical system. Mathematically, a phase transition occurs when the partition function vanishes and the free energy is singular (non-analytic). For instance, if the first derivative of the free energy with respect to the control parameter is non-continuous, a jump may occur in the average value of the fluctuating conjugate variable, such as the magnetization, corresponding to a first-order phase transition.

Importantly, for a finite-size system,  is a finite sum of exponential functions and is thus always positive for real values of . Consequently,  is always well-behaved and analytic for finite system sizes. By contrast, in the thermodynamic limit,  may exhibit a non-analytic behavior.

Using that  is an entire function for finite system sizes, Lee–Yang theory takes advantage of the fact that the partition function can be fully characterized by its zeros in the complex plane of . These zeros are often known as Lee–Yang zeros or, in the case of inverse temperature as control parameter, Fisher zeros. The main idea of Lee–Yang theory is to mathematically study how the positions and the behavior of the zeros change as the system size grows. If the zeros move onto the real axis of the control parameter in the thermodynamic limit, it signals the presence of a phase transition at the corresponding real value of .

In this way, Lee–Yang theory establishes a connection between the properties (the zeros) of a partition function for a finite size system and phase transitions that may occur in the thermodynamic limit (where the system size goes to infinity).

Examples

Molecular zipper 
The molecular zipper is a toy model which may be used to illustrate the Lee–Yang theory. It has the advantage that all quantities, including the zeros, can be computed analytically. The model is based on a double-stranded macromolecule with  links that can be either open or closed. For a fully closed zipper, the energy is zero, while for each open link the energy is increased by an amount . A link can only be open if the preceding one is also open.

For a number  of different ways that a link can be open, the partition function of a zipper with  links reads
.
This partition function has the complex zeros

where we have introduced the critical inverse temperature , with . We see that in the limit , the zeros closest to the real axis approach the critical value . For , the critical temperature is infinite and no phase transition takes place for finite temperature. By contrast, for , a phase transition takes place at the finite temperature .

To confirm that the system displays a non-analytic behavior in the thermodynamic limit, we consider the free energy

or, equivalently, the dimensionless free energy per link 
In the thermodynamic limit, one obtains
.
Indeed, a cusp develops at  in the thermodynamic limit. In this case, the first derivative of the free energy is discontinuous, corresponding to a first-order phase transition.

Ising model 

The Ising model is the original model that Lee and Yang studied when they developed their theory on partition function zeros. The Ising model consists of spin lattice with  spins , each pointing either up, , or down, . Each spin may also interact with its closest spin neighbors with a strength . In addition, an external magnetic field  may be applied (here we assume that it is uniform and thus independent of the spin indices). The Hamiltonian of the system for a certain spin configuration  then reads
 
In this case, the partition function reads

The zeros of this partition function cannot be determined analytically, thus requiring numerical approaches.

Lee–Yang theorem 

For the ferromagnetic Ising model, for which  for all , Lee and Yang showed that all zeros of  lie on the unit circle in the complex plane of the parameter . This statement is known as the Lee–Yang theorem, and has later been generalized to other models, such as the Heisenberg model.

Dynamical phase transitions 

A similar approach can be used to study dynamical phase transitions. These transitions are characterized by the Loschmidt amplitude, which plays the analogue role of a partition function.

Connections to fluctuations 
The Lee–Yang zeros may be connected to the cumulants of the conjugate variable  of the control variable . For brevity, we set  in the following. Using that the partition function is an entire function for a finite-size system, one may expand it in terms of its zeros as

where  and  are constants, and  is the :th zero in the complex plane of . The corresponding free energy then reads

Differentiating this expression  times with respect to , yields the :th order cumulant

Furthermore, using that the partition function is a real function, the Lee–Yang zeros have to come in complex conjugate pairs, allowing us to express the cumulants as

where the sum now runs only over each pair of zeros. This establishes a direct connection between cumulants and Lee–Yang zeros.

Moreover, if  is large, the contribution from zeros lying far away from  is strongly suppressed, and only the closest pair  of zeros plays an important role. One may then write

This equation may be solved as a linear system of equations, allowing for the Lee–Yang zeros to be determined directly from higher-order cumulants of the conjugate variable:

Experiments 
Being complex numbers of a physical variable, Lee–Yang zeros have traditionally been seen as a purely theoretical tool to describe phase transitions, with little or none connection to experiments. However, in a series of experiments in the 2010s, various kinds of Lee–Yang zeros have been determined from real measurements. In one experiment in 2015, the Lee–Yang zeros were extracted experimentally by measuring the quantum coherence of a spin coupled to an Ising-type spin bath. In another experiment in 2017, dynamical Lee–Yang zeros were extracted from Andreev tunneling processes between a normal-state island and two superconducting leads. Furthermore, in 2018, there was an experiment determining the dynamical Fisher zeros of the Loschmidt amplitude, which may be used to identify dynamical phase transitions.

See also 
 Lee–Yang theorem

References 

Phase transitions