Experimental Determination of Dynamical Lee-Yang Zeros

Yang-Lee zeros of certain antiferromagnetic models

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.

Experimental Observation of the Yang-Lee Quantum Criticality in Open Quantum Systems

The Yang-Lee edge singularity was originally studied from the standpoint of mathematical foundations of phase transitions. However, direct observation of anomalous scaling with the negative scaling dimension has remained elusive due to an imaginary magnetic field required for the nonunitary criticality. We experimentally implement an imaginary magnetic field with an open quantum system of heralded single photons, directly measure the partition function, and demonstrate the Yang-Lee edge singularity via the quantum-classical correspondence. We also demonstrate unconventional scaling laws for finite-temperature quantum dynamics.

Complex phase diagram and supercritical matter

The supercritical region is often described as uniform with no definite transitions. The distinct behaviors of the matter therein, e.g., as liquidlike and gaslike, however, suggest "supercritical boundaries." Here we provide a mathematical description of these phenomena by revisiting the Yang-Lee theory and introducing a complex phase diagram, specifically a four-dimensional (4D) one with complex T and p. While the traditional 2D phase diagram with real temperature T and pressure p values (the physical plane) lacks Lee-Yang (LY) zeros beyond the critical point, preventing the occurrence of criticality, the off-plane zeros in this 4D scenario still induce critical anomalies in various physical properties. This relationship is evidenced by the correlation between the Widom line and LY edges in van der Waals, 2D Ising model, and water. The diverged supercritical boundaries manifest the high-dimensional feature of the phase diagram: e.g., when LY zeros of complex T or p are projected onto the physical plane, boundaries defined by isobaric heat capacity C_{p} or isothermal compression coefficient K_{T} emanates. These results demonstrate the incipient phase transition nature of the supercritical matter.

Proposal for Observing Yang-Lee Criticality in Rydberg Atomic Arrays

Yang-Lee edge singularities (YLES) are the edges of the partition function zeros of an interacting spin model in the space of complex control parameters. They play an important role in understanding non-Hermitian phase transitions in many-body physics, as well as characterizing the corresponding nonunitary criticality. Even though such partition function zeroes have been measured in dynamical experiments where time acts as the imaginary control field, experimentally demonstrating such YLES criticality with a physical imaginary field has remained elusive due to the difficulty of physically realizing non-Hermitian many-body models. We provide a protocol for observing the YLES by detecting kinked dynamical magnetization responses due to broken PT symmetry, thus enabling the physical probing of nonunitary phase transitions in nonequilibrium settings. In particular, scaling analyses based on our nonunitary time evolution circuit with matrix product states accurately recover the exponents uniquely associated with the corresponding nonunitary CFT. We provide an explicit proposal for observing YLES criticality in Floquet quenched Rydberg atomic arrays with laser-induced loss, which paves the way towards a universal platform for simulating non-Hermitian many-body dynamical phenomena.

Lee-Yang zeros and quantum Fisher information matrix in a nonlinear system

The distribution of Lee-Yang zeros not only matters in thermodynamics and quantum mechanics, but also in mathematics. Hereby we propose a nonlinear quantum toy model and discuss the distribution of corresponding Lee-Yang zeros. Utilizing the coupling between a probe qubit and the nonlinear system, all Lee-Yang zeros can be detected in the dynamics of the probe qubit by tuning the coupling strength and linear coefficient of the nonlinear system. Moreover, the analytical expression of the quantum Fisher information matrix at the Lee-Yang zeros is provided and an interesting phenomenon is discovered. Both the coupling strength and temperature can simultaneously attain their precision limits at the Lee-Yang zeros. However, the probe qubit cannot work as a thermometer at a Lee-Yang zero if it sits on the unit circle.

Complex analysis of divergent perturbation theory at finite temperature

We investigate the convergence properties of finite-temperature perturbation theory by considering the mathematical structure of thermodynamic potentials using complex analysis. We discover that zeros of the partition function lead to poles in the internal energy and logarithmic singularities in the Helmholtz free energy that create divergent expansions in the canonical ensemble. Analyzing these zeros reveals that the radius of convergence increases at higher temperatures. In contrast, when the reference state is degenerate, these poles in the internal energy create a zero radius of convergence in the zero-temperature limit. Finally, by showing that the poles in the internal energy reduce to exceptional points in the zero-temperature limit, we unify the two main mathematical representations of quantum phase transitions.

Exact Anomalous Current Fluctuations in a Deterministic Interacting Model

We analytically compute the full counting statistics of charge transfer in a classical automaton of interacting charged particles. Deriving a closed-form expression for the moment generating function with respect to a stationary equilibrium state, we employ asymptotic analysis to infer the structure of charge current fluctuations for a continuous range of timescales. The solution exhibits several unorthodox features. Most prominently, on the timescale of typical fluctuations the probability distribution of the integrated charge current in a stationary ensemble without bias is distinctly non-Gaussian despite diffusive behavior of dynamical charge susceptibility. While inducing a charge imbalance is enough to recover Gaussian fluctuations, we find that higher cumulants grow indefinitely in time with different exponents, implying singular scaled cumulants. We associate this phenomenon with the lack of a regularity condition on moment generating functions and the onset of a dynamical critical point. In effect, the scaled cumulant generating function does not, irrespectively of charge bias, represent a faithful generating function of the scaled cumulants, yet the associated Legendre dual yields the correct large-deviation rate function. Our findings hint at novel types of dynamical universality classes in deterministic many-body systems.

Dynamical Lee-Yang zeros for continuous-time and discrete-time stochastic processes

We describe classical stochastic processes by using dynamical Lee-Yang zeros. The system is in contact with external leads and the time evolution is described by the two-state classical master equation. The cumulant generating function is written in a factorized form and the current distribution is characterized by the dynamical Lee-Yang zeros. We show that a continuous distribution of zeros is obtained by discretizing the time variable. When the transition probability is a periodically oscillating function of time, the distribution of zeros splits into many parts. We study the geometric property of the current by comparing the result with that of the adiabatic approximation. We also use the Floquet-Magnus expansion in the continuous-time case to study dynamical effects on the current at the fast-driving regime.

Moment-generating function zeros in the study of phase transitions

Partition function zeros play a central role in the study of phase transitions. Recently, energy probability distribution (EPD) zeros were proposed as an alternative approach that solves some of the implementation issues present in the Fisher zeros method by allowing drastic reduction of the polynomial. Here, a formulation based on the EPD zeros that can reduce even more the polynomial degree while maintaining its accuracy is presented. This method has shown to be computationally cheaper than the EPD zeros, allowing the study of systems by using partition function zeros that would be unfeasible otherwise. In addition, the method can be easily extended to study phase transitions in external fields while maintaining all of its improvements.

Many-body thermodynamics on quantum computers via partition function zeros

Partition functions are ubiquitous in physics: They are important in determining the thermodynamic properties of many-body systems and in understanding their phase transitions. As shown by Lee and Yang, analytically continuing the partition function to the complex plane allows us to obtain its zeros and thus the entire function. Moreover, the scaling and nature of these zeros can elucidate phase transitions. Here, we show how to find partition function zeros on noisy intermediate-scale trapped-ion quantum computers in a scalable manner, using the XXZ spin chain model as a prototype, and observe their transition from XY-like behavior to Ising-like behavior as a function of the anisotropy. While quantum computers cannot yet scale to the thermodynamic limit, our work provides a pathway to do so as hardware improves, allowing the future calculation of critical phenomena for systems beyond classical computing limits.

Classification of dark solitons via topological vector potentials

Dark solitons are some of the most interesting nonlinear excitations and are considered to be the one-dimensional topological analogs of vortices. However, in contrast to their two-dimensional vortex counterparts, the topological characteristics of a dark soliton are far from fully understood because the topological charge defined according to the phase jump cannot reflect its essential property. Here, similar to the complex extension used in the exploration of the partition-function zeros to depict thermodynamic states, we extend the complex coordinate space to explore the density zeros of dark solitons. Surprisingly we find that these zeros constitute some pointlike magnetic fields, each of which has a quantized magnetic flux of elementary π. The corresponding vector potential fields demonstrate the topology of the Wess-Zumino term and can depict the essential characteristics of dark solitons. Then we classify the dark solitons according to the Euler characteristic of the topological manifold of the vector potential fields. Our study not only reveals the topological features of dark solitons but can also be applied to explore and identify new dark solitons with high topological complexity.

Building Continuous Time Crystals from Rare Events

Symmetry-breaking dynamical phase transitions (DPTs) abound in the fluctuations of nonequilibrium systems. Here, we show that the spectral features of a particular class of DPTs exhibit the fingerprints of the recently discovered time-crystal phase of matter. Using Doob's transform as a tool, we provide a mechanism to build classical time-crystal generators from the rare event statistics of some driven diffusive systems. An analysis of the Doob's smart field in terms of the order parameter of the transition then leads to the time-crystal lattice gas (TCLG), a model of driven fluid subject to an external packing field, which presents a clear-cut steady-state phase transition to a time-crystalline phase characterized by a matter density wave, which breaks continuous time-translation symmetry and displays rigidity and long-range spatiotemporal order, as required for a time crystal. A hydrodynamic analysis of the TCLG transition uncovers striking similarities, but also key differences, with the Kuramoto synchronization transition. Possible experimental realizations of the TCLG in colloidal fluids are also discussed.

Pseudo-Yang-Lee Edge Singularity Critical Behavior in a Non-Hermitian Ising Model

The quantum phase transition of a one-dimensional transverse field Ising model in an imaginary longitudinal field is studied. A new order parameter M is introduced to describe the critical behaviors in the Yang-Lee edge singularity (YLES). The M does not diverge at the YLES point, a behavior different from other usual parameters. We term this unusual critical behavior around YLES as the pseudo-YLES. To investigate the static and driven dynamics of M, the (1+1) dimensional ferromagnetic-paramagnetic phase transition ((1+1) D FPPT) critical region, (0+1) D YLES critical region and the (1+1) D YLES critical region of the model are selected. Our numerical study shows that the (1+1) D FPPT scaling theory, the (0+1) D YLES scaling theory and (1+1) D YLES scaling theory are applicable to describe the critical behaviors of M, demonstrating that M could be a good indicator to detect the phase transition around YLES. Since M has finite value around YLES, it is expected that M could be quantitatively measured in experiments.

Transverse-size critical exponent of directed percolation from Yang-Lee zeros of survival probability

By using transfer-matrix method we compute survival probabilities for the directed percolation problem on strips of a square lattice, and get very precise estimates of their Yang-Lee zeros lying closest to the real axis in the complex plane of occupation probability. This allows us to get accurate values for transverse-size critical exponent and percolation threshold.

Zeros of partition functions in the NPT ensemble

Lee-Yang and Fisher zeros are crucial for the study of phase transitions in the grand canonical and the canonical ensembles, respectively. However, these powerful methods do not cover the isothermal-isobaric ensemble (NPT ensemble), which reflects the conditions of many experiments. In this work we present a theory of the phase transitions in terms of the zeros of the NPT-ensemble partition functions in the complex plane. The proposed theory provides an approach to calculate all the partition function zeros in the NPT ensemble, which form certain curves in the thermodynamic limit. To verify the theory we consider Tonks gas and van der Waals fluid in the NPT ensemble. In the case of Tonks gas, similarly to the Lee-Yang circle theorem, we obtain an exact equation for the zero limit curve. We also derive an approximated limit curve equation for van der Waals fluid in terms of the Szegö curve. This curve fits numerically calculated zeros and correctly describes how the phenomenon of phase transition depends on the temperature.

Exact large deviation statistics and trajectory phase transition of a deterministic boundary driven cellular automaton

We study the statistical properties of the long-time dynamics of the rule 54 reversible cellular automaton (CA), driven stochastically at its boundaries. This CA can be considered as a discrete-time and deterministic version of the Fredrickson-Andersen kinetically constrained model (KCM). By means of a matrix product ansatz, we compute the exact large deviation cumulant generating functions for a wide range of time-extensive observables of the dynamics, together with their associated rate functions and conditioned long-time distributions over configurations. We show that for all instances of boundary driving the CA dynamics occurs at the point of phase coexistence between competing active and inactive dynamical phases, similar to what happens in more standard KCMs. We also find the exact finite size scaling behavior of these trajectory transitions, and provide the explicit "Doob-transformed" dynamics that optimally realizes rare dynamical events.

Sampling rare events across dynamical phase transitions

Interacting particle systems with many degrees of freedom may undergo phase transitions to sustain atypical fluctuations of dynamical observables such as the current or the activity. In some cases, this leads to symmetry-broken space-time trajectories which enhance the probability of such events due to the emergence of ordered structures. Despite their conceptual and practical importance, these dynamical phase transitions (DPTs) at the trajectory level are difficult to characterize due to the low probability of their occurrence. However, during the last decade, advanced computational techniques have been developed to measure rare events in simulations of many-particle systems that allow the direct observation and characterization of these DPTs. Here we review the application of a particular rare-event simulation technique, based on cloning Monte Carlo methods, to characterize DPTs in paradigmatic stochastic lattice gases. In particular, we describe in detail some tricks and tips of the trade, paying special attention to the measurement of order parameters capturing the physics of the different DPTs, as well as to the finite-size effects (both in the system size and in the number of clones) that affect the measurements. Overall, we provide a consistent picture of the phenomenology associated with DPTs and their measurement.

Simulating Dynamic Quantum Phase Transitions in Photonic Quantum Walks

Signaled by nonanalyticities in the time evolution of physical observables, dynamic quantum phase transitions (DQPTs) emerge in quench dynamics of topological systems and possess an interesting geometric origin captured by dynamic topological order parameters (DTOPs). In this Letter, we report the experimental study of DQPTs using discrete-time quantum walks of single photons. We simulate quench dynamics between distinct Floquet topological phases using quantum-walk dynamics and experimentally characterize DQPTs and the underlying DTOPs through interference-based measurements. The versatile photonic quantum-walk platform further allows us to experimentally investigate DQPTs for mixed states and in parity-time-symmetric nonunitary dynamics for the first time. Our experiment directly confirms the relation between DQPTs and DTOPs in quench dynamics of topological systems and opens up the avenue of simulating emergent topological phenomena using discrete-time quantum-walk dynamics.

Dynamical quantum phase transitions: a review

Quantum theory provides an extensive framework for the description of the equilibrium properties of quantum matter. Yet experiments in quantum simulators have now opened up a route towards the generation of quantum states beyond this equilibrium paradigm. While these states promise to show properties not constrained by equilibrium principles, such as the equal a priori probability of the microcanonical ensemble, identifying the general properties of nonequilibrium quantum dynamics remains a major challenge, especially in view of the lack of conventional concepts such as free energies. The theory of dynamical quantum phase transitions attempts to identify such general principles by lifting the concept of phase transitions to coherent quantum real-time evolution. This review provides a pedagogical introduction to this field. Starting from the general setting of nonequilibrium dynamics in closed quantum many-body systems, we give the definition of dynamical quantum phase transitions as phase transitions in time with physical quantities becoming nonanalytic at critical times. We summarize the achieved theoretical advances as well as the first experimental observations, and furthermore provide an outlook to major open questions as well as future directions of research.

Lee-Yang zeros and large-deviation statistics of a molecular zipper

The complex zeros of partition functions were originally investigated by Lee and Yang to explain the behavior of condensing gases. Since then, Lee-Yang zeros have become a powerful tool to describe phase transitions in interacting systems. Today, Lee-Yang zeros are no longer just a theoretical concept; they have been determined in recent experiments. In one approach, the Lee-Yang zeros are extracted from the high cumulants of thermodynamic observables at finite size. Here we employ this method to investigate a phase transition in a molecular zipper. From the energy fluctuations in small zippers, we can predict the temperature at which a phase transition occurs in the thermodynamic limit. Even when the system does not undergo a sharp transition, the Lee-Yang zeros carry important information about the large-deviation statistics and its symmetry properties. Our work suggests an interesting duality between fluctuations in small systems and their phase behavior in the thermodynamic limit. These predictions may be tested in future experiments.

Quantum work relations and response theory in parity-time-symmetric quantum systems

In this work, we show that a universal quantum work relation for a quantum system driven arbitrarily far from equilibrium extends to a parity-time- (PT-) symmetric quantum system with unbroken PT symmetry, which is a consequence of microscopic reversibility. The quantum Jarzynski equality, linear response theory, and Onsager reciprocal relations for the PT-symmetric quantum system are recovered as special cases of the universal quantum work relation in a PT-symmetric quantum system. In the regime of broken PT symmetry, the universal quantum work relation does not hold because the norm is not preserved during the dynamics.

Partition function zeros of the p-state clock model in the complex temperature plane

We investigate the partition function zeros of the two-dimensional p-state clock model in the complex temperature plane by using the Wang-Landau method. For p=5, 6, 8, and 10, we propose a modified energy representation to enumerate exact irregular energy levels for the density of states without any binning artifacts. Comparing the leading zeros between different p's, we provide strong evidence that the upper transition at p=6 is indeed of the Berezinskii-Kosterlitz-Thouless (BKT) type in contrast to the claim of the previous Fisher zero study [Phys. Rev. E 80, 042103 (2009)10.1103/PhysRevE.80.042103]. We find that the leading zeros of p=6 at the upper transition collapse onto the zero trajectories of the larger p's including the XY limit while the finite-size behavior of p=5 differs from the converged behavior of p≥6 within the system sizes examined. In addition, we argue that the nondivergent specific heat in the BKT transition is responsible for the small partition function magnitude that decreases exponentially with increasing system size near the leading zero, fundamentally limiting access to large systems in search for zeros with an estimator under finite statistical fluctuations.

Order and Symmetry Breaking in the Fluctuations of Driven Systems

Dynamical phase transitions (DPTs) in the space of trajectories are one of the most intriguing phenomena of nonequilibrium physics, but their nature in realistic high-dimensional systems remains puzzling. Here we observe for the first time a DPT in the current vector statistics of an archetypal two-dimensional (2D) driven diffusive system and characterize its properties using the macroscopic fluctuation theory. The complex interplay among the external field, anisotropy, and vector currents in 2D leads to a rich phase diagram, with different symmetry-broken fluctuation phases separated by lines of first- and second-order DPTs. Remarkably, different types of 1D order in the form of jammed density waves emerge to hinder transport for low-current fluctuations, revealing a connection between rare events and self-organized structures which enhance their probability.