Tracking metastable phases by complex Lee-Yang zeros

This paper demonstrates that metastable phases, typically suppressed in equilibrium, can be tracked and characterized as regions in the complex plane of thermal fields delineated by Lee-Yang zeros, offering a new framework for understanding and engineering non-equilibrium collective states in periodically driven systems.

The Gist

Imagine you are looking at a map of a country. Usually, this map only shows you the major cities where people live and thrive—these are the "stable phases" of matter, like ice, water, or steam. But hidden beneath the surface, in the deep valleys and foggy mountains, there are other places where people could live, but they are usually too unstable to stay there for long. These are "metastable phases." In the old way of doing science, these hidden places were like an "iceberg" below the waterline: we knew they might be there, but our standard maps couldn't show them until they suddenly popped up on the surface.

This paper introduces a new kind of "super-map" that can see the hidden parts of the iceberg before they reach the surface.

The Problem: The "Hidden" States

Think of matter like a ball rolling down a hill. It naturally settles in the deepest valley (the stable state). Sometimes, the ball gets stuck in a shallow dip halfway down the hill. It's not at the bottom, but it's not rolling away either. This is a metastable phase.

• The Old View: Standard maps (equilibrium phase diagrams) only show the deepest valleys. If the ball is stuck in the shallow dip, the map says, "Nothing here, just a slope." The shallow dip is invisible until the ball finally rolls into it and becomes a permanent city.
• The Challenge: Scientists want to find and control these shallow dips because they often have special, exotic properties that the deep valleys don't have. But finding them is like trying to find a ghost; they are hard to spot and hard to keep.

The Solution: The "Ghost Map" (Lee-Yang Zeros)

The authors propose using a mathematical tool called Lee-Yang zeros.

• The Analogy: Imagine the standard map is a 2D drawing on a flat piece of paper. The Lee-Yang method adds a third dimension: a "depth" axis.
• In this new 3D space, the "ghosts" (metastable phases) aren't invisible. They appear as specific patterns or "fences" in the complex, deeper part of the map.
• Even when the shallow dip is too unstable to exist on the flat paper (the real world), the "fences" in the 3D depth are already there, outlining exactly where that hidden state lives.

How They Proved It: The Three-Hill Model

To test this, the scientists built a simple computer model (a "toy model") with three hills:

1. Hill A and Hill C are the big, stable cities.
2. Hill B is the small, shaky hill in the middle (the metastable phase).

What happened in the simulation:

• Step 1: They started with Hill B being very weak. On the flat map, you only saw a transition from A to C. Hill B was invisible.
• Step 2: They slowly made Hill B stronger (by adjusting a knob).
• The Magic: While Hill B was still too weak to be seen on the flat map, the "Ghost Map" (the complex plane) showed a new fence appearing deep in the 3D space. As they turned the knob, this fence moved closer to the surface.
• The Result: The moment Hill B became strong enough to be a real city on the flat map, the fence from the 3D depth finally touched the surface and split, creating a clear boundary for the new city.

The Takeaway: The "Ghost Map" didn't just show the city after it appeared; it tracked the entire journey of the city from being a hidden ghost to becoming a real place.

The Real-World Test: Shaking with Light

The scientists then tried this on a more realistic system using terahertz light (a type of high-frequency vibration).

• Imagine shaking a box of marbles. If you shake it just right, you can make the marbles settle into a pattern they wouldn't normally choose.
• They used light to "shake" the material, effectively changing the landscape of the hills.
• They found that the strength of the shaking (the drive) was directly linked to the position of the "fences" in their Ghost Map.
• The Connection: By treating the shaking light as a "complex temperature," they could predict exactly how to make the hidden, metastable state appear and stay stable.

Why This Matters

This paper claims that we don't have to wait for a material to accidentally become stable to study it.

• The New Perspective: Stable phases are actually the "accidents" that happen to land on the flat surface. The "real" world of matter is the complex 3D space where all these states exist.
• The Benefit: By looking at the "Ghost Map" (the complex Lee-Yang zeros), scientists can proactively design materials. They can see the hidden states, understand how to stabilize them, and engineer new materials with special properties before they even exist in the real world.

In short, the paper says: Stop looking only at the surface. If you look deep enough into the mathematical "fog," you can see the hidden cities of matter long before they arrive.

Technical Summary

Technical Summary: Tracking Metastable Phases by Complex Lee-Yang Zeros

Problem Statement
Metastable phases (MPs) represent energetically unfavorable states that are typically suppressed in equilibrium phase diagrams (EPDs). While stable phases are reliably mapped by EPDs, the exploration of MPs has historically relied on chemical intuition or computationally exhaustive structure predictions. A fundamental challenge arises because MPs inherently correlate with dissipation, fluctuating orders, and finite lifetimes, rendering standard equilibrium thermodynamic hypotheses often incongruous or misleading when applied to them. Specifically, the thermodynamic limit tends to exaggerate free energy differences, effectively overlooking MPs entirely. Furthermore, existing theoretical frameworks (e.g., nucleation theory) focus primarily on the kinetics of decay rather than the intrinsic existence of MPs as reproducible states. A critical dilemma exists in phase diagram analysis: a single point in an EPD may correspond to multiple distinct states, necessitating additional degrees of freedom (DOFs) to distinguish them. Current simulation methods (e.g., kinetic Monte Carlo, quenched molecular dynamics) often depend heavily on prior knowledge of pathways or initial conditions, limiting proactive and systematic probing.

Methodology
The authors propose a framework to uncover "hidden" metastable phases by analyzing Lee-Yang zeros (LYZs) within a complex phase diagram. This approach utilizes the partition function, which retains full state statistics including MP signatures, and converts this information into a readable phase diagram structure via LYZs. The study employs two primary models:

1. Three-Gaussian-Peak Toy Model: A minimal model where the reference density of states (DOS), $P(E, \beta_0)$, consists of three Gaussian peaks representing three potential phases (I, II, and III). An artificial parameter, $A_2$, tunes the strength of the intermediate peak (Phase II). The authors analytically extend the inverse temperature $\beta$ to a complex variable $\tilde{\beta} = \text{Re}(\tilde{\beta}) + i\text{Im}(\tilde{\beta})$ and solve for the roots of the partition function $Z(\tilde{\beta}) = 0$.
2. Periodically Driven System: A more realistic model involving a three-well potential subjected to an oscillatory field (terahertz-like manipulation). The system is simulated using Langevin molecular dynamics (MD). The drive amplitude ($M_1$) is tuned to alter the relative heights of the DOS peaks, replacing the artificial parameter $A_2$ with an experimentally controllable field. The LYZs are computed by discretizing the sampled DOS and decomposing the partition function into a polynomial.

Key Contributions and Results
The paper demonstrates that LYZ structures in the complex plane can faithfully track MPs, delineating metastable regions even before they stabilize on the real axis.

• Complex Phase Diagram Structure: In the toy model, as the parameter $A_2$ increases, the LYZs bounding the metastable Phase II approach the real axis. Initially, for low $A_2$, the LYZs form a single boundary dividing the complex plane. As $A_2$ increases, this boundary branches, creating a new domain at large imaginary $\tilde{\beta}$ values. Eventually, the branching point approaches the real axis, forming a triple point, after which the boundary splits into two separated branches. This splitting signals the emergence and stabilization of Phase II between Phases I and III.
• Tracking Metastability: The study quantifies MP stability using $d_{BP}$, the distance from the real axis to the nearest complex branching point. A finite $d_{BP}$ indicates a metastable regime, which decreases monotonically to zero as the phase stabilizes. The complex phase diagram captures the entire stabilization process, revealing Phase II as a distinct region in the complex plane long before it appears in the real-axis EPD.
• Driven System Correlation: In the periodically driven system, increasing the drive strength ($M_1$) reproduces the stabilization of the hidden Phase II. The MD simulations show that as $M_1$ increases, the central peak of the DOS (Phase II) grows and eventually prevails. The LYZ analysis confirms this: the region associated with Phase II in the complex plane approaches and touches the real axis, while the dominant $\beta$-ranges of Phases I and III narrow.
• Quantitative Link: The study establishes a quantitative correlation between the drive strength and the LYZ structure. The imaginary part of the LYZs correlates with the drive strength, and the gap between branches ($\Delta\beta$) tracks the stability window of the MP. This links Lee-Yang theory directly to terahertz matter manipulation.

Significance and Claims
The authors claim that this work provides a scheme to describe MPs in phase diagram analysis by viewing periodic drives as complex thermal fields. Key claims include:

• Unified Framework: The complex phase diagram places stable phases and MPs on the same map, treating them as connected domains delineated by LYZs. This suggests that stable phases may be the exception rather than the rule, as only a special subset of LYZ-delineated regions "coincidentally" claims the real axis.
• Proactive Engineering: By viewing the complex plane as a landscape where MPs occupy well-defined locations, researchers can proactively engineer MPs. This involves predicting and scanning MP candidates via DOS simulations and LYZ analysis before stabilization occurs.
• Theory-Experiment Connection: The correspondence between the drive strength and the complex thermal field offers a statistical interpretation of oscillatory driving. This establishes a concrete connection between Lee-Yang theory and experimental approaches, such as terahertz manipulation, for accessing and stabilizing non-equilibrium collective states.
• Beyond Equilibrium: The findings extend the utility of LYZ analysis beyond traditional equilibrium thermodynamics (where zeros pinch the real axis) to dynamical phase transitions and non-equilibrium phenomena, utilizing imaginary thermal fields as the necessary extra degrees of freedom.

The paper concludes that this framework offers a practical route to understand and detect driven metastable phases, potentially enabling the design of novel metastable materials for advanced electronics and quantum technologies.