Field-selective criticality in 2D melting revealed by multi-field Lee-Yang zeros

By applying Lee-Yang zeros to study bilayer water under nanoconfinement, this paper reveals that 2D melting exhibits field-selective criticality where the solid-hexatic and hexatic-liquid transitions respond differently to temperature and lateral pressure, thereby resolving long-standing contradictions between simulations, models, and experiments by identifying which thermodynamic channel each probe observes.

The Gist

Imagine you are trying to figure out exactly when a block of ice turns into water. In the everyday world, this seems simple: heat it up, and it melts. But in the microscopic world of two-dimensional materials (like a single sheet of water molecules trapped between two walls), scientists have been arguing for 60 years about how this happens.

Some say it melts smoothly, like butter softening. Others say it snaps suddenly, like a glass shattering. Some even say it happens in two distinct steps. The problem is that different scientists have been looking at the ice through different "lenses," and each lens showed a slightly different picture.

This paper, by researchers at Peking University, acts like a master key that unlocks the confusion. They didn't just look at the ice; they looked at how the ice reacts to two different forces at the same time: heat (temperature) and squeeze (pressure).

Here is the story of what they found, explained simply:

1. The "Two-Lens" Problem

Imagine you are watching a magician pull a rabbit out of a hat.

• Lens A (The Heat Lens): You watch the rabbit's temperature.
• Lens B (The Squeeze Lens): You watch how much space the rabbit takes up.

In most situations, if the rabbit changes, both lenses see the change at the exact same moment. But the researchers discovered that for this trapped water, the lenses can disagree. Sometimes, the rabbit looks like it's changing size before it changes temperature, or vice versa.

In the past, scientists only looked through one lens. If they watched the heat, they saw a smooth transition. If they watched the squeeze, they saw a sudden jump. This led to the argument: "Is it smooth or sudden?" The answer, this paper says, is: "It depends on which lens you are using."

2. The "Field-Selective" Melting

The team used a sophisticated mathematical tool called Lee-Yang zeros. Think of this as a super-sensitive radar that can detect the exact moment a phase change happens, even if it's blurry.

They found two types of melting behavior in their trapped water:

• The "Split" Melting (Field-Selective):
  Imagine a crowd of people (water molecules) trying to leave a room.

  • When you look at how much space they take up (density), they seem to leave the room gradually, one by one, like a slow stream.
  • But when you look at how much energy they have (enthalpy), they all rush out at once, like a stampede.
  • The Discovery: The researchers found that for certain types of ice, the "space" lens sees a smooth transition, while the "energy" lens sees a sudden jump. This is called field-selective criticality. It means the transition is "sudden" for one observer and "smooth" for another.

• The "Two-Step" Melting:
  For other conditions, the ice doesn't melt all at once. It goes through a weird middle stage called a "hexatic" phase.

  • Think of this like a dance floor. First, the dancers (molecules) are frozen in a rigid grid (Solid).
  • Then, they break the grid but still hold hands in a circle, moving loosely (Hexatic).
  • Finally, they let go completely and run around wildly (Liquid).
  • Previous studies argued about whether the step from "Grid" to "Circle" was smooth or sudden. The researchers found that if you use a small "camera" (a small simulation), the jump looks blurry and smooth. But if you use a huge camera (a much larger simulation with 1,000+ molecules), the jump becomes crystal clear. It turns out the first step is actually a sudden jump, just a very subtle one that gets hidden in small experiments.

3. Why This Matters

The paper solves a decades-old mystery by showing that there is no single "truth" about how 2D ice melts unless you specify how you are measuring it.

• The Confusion: Previous experiments and computer simulations seemed to contradict each other. One said "smooth," another said "sudden."
• The Resolution: They were all right, but they were looking at different things. The "smooth" observers were looking at density (space), and the "sudden" observers were looking at energy.
• The New Picture: The researchers mapped out a new "weather map" for this water. They showed exactly where the "split" happens (where heat and pressure disagree) and where the "two-step" dance occurs.

The Takeaway

This paper is like realizing that a chameleon isn't just "green" or "brown." It changes color depending on the background. Similarly, 2D ice doesn't have a single way of melting. It has a dual personality: it can melt smoothly if you watch its size, but suddenly if you watch its energy.

By using advanced math to look at both "lenses" simultaneously, the authors finally organized the conflicting stories into one clear, unified picture. They didn't just find where the ice melts; they explained why everyone saw it differently.

Technical Summary

Technical Summary: Field-Selective Criticality in 2D Melting

Problem Statement
The nature of two-dimensional (2D) melting has remained a subject of debate for six decades, with conflicting scenarios proposed by theory, model systems, simulations, and atomic-resolution experiments. A central ambiguity arises because the same phase transition can appear continuous or abrupt depending on the thermodynamic observable used to diagnose it. This is particularly acute in confined systems, such as bilayer water, where different studies probe different thermodynamic responses: some rely on isotherms or density/volume changes (probing lateral pressure, $p_L$), while others track potential energy evolution with temperature (probing the thermal axis). The paper argues that if responses to temperature ($T$) and lateral pressure ($p_L$) are not locked, a single observable may yield an incomplete or misleading classification of the transition order. Furthermore, finite-size effects in simulations can round weak first-order transitions, making them appear continuous.

Methodology
To resolve these ambiguities, the authors employ a field-resolved approach using Lee-Yang zeros (LYZs) combined with enhanced sampling.

• System: The study focuses on bilayer water confined between two hydrophobic walls (separation 0.8 nm), modeled using the TIP4P/2005 rigid potential. Simulations were conducted in the isothermal-isolateral pressure ($Np_LT$) ensemble.
• Sampling: The authors utilized the OPES (On-the-fly Probability Enhanced Sampling) method to obtain well-converged densities of states (DOS) for enthalpy ($H$) and volume ($V$). Both standalone OPES-explore and combined multithermal–multibaric (MTMB) simulations were used to map the phase diagram.
• Lee-Yang Analysis: Partition functions were constructed as functions of complex temperature (at fixed real $p_L$) and complex lateral pressure (at fixed real $T$). The zeros of these partition functions (LYZs) were calculated. The real parts of the leading zeros (Lee-Yang edges) locate the phase boundaries, while their imaginary parts characterize the approach to singularity in finite systems.
• Finite-Size Scaling: To address rounding effects, the study compared a 256-molecule system with a 1024-molecule supercell. The extensivity of entropy was used to scale the DOS for larger systems ($m > 1$) to examine convergence toward the thermodynamic limit.

Key Results
The study constructs a $T$-$p_L$ phase diagram for bilayer nanoconfined water containing liquid, hexatic, square ice, and hexagonal ice phases. The analysis reveals that phase boundaries can be field-selective, meaning the thermodynamic character of a transition depends on whether it is probed via the thermal ($T$) or mechanical ($p_L$) channel.

1. One-Step Melting (Routes I & II):

   • Route I (Hexagonal Ice–Liquid): Near the maximum of the coexistence line (around $p_L \approx 70$ MPa), the authors observe a bifurcation between the $T$-edges and $p_L$-edges. While the enthalpy distribution ($H$) remains bimodal (indicating a discontinuous thermal response), the volume distribution ($V$) becomes unimodal due to the overlap of peaks ($\Delta V \to 0$). Consequently, the density response appears continuous, while the enthalpy response remains first-order. This explains why different probes yield conflicting classifications.
   • Route II (Square Ice–Liquid): This transition behaves conventionally, with both $H$ and $V$ distributions remaining bimodal and the $T$ and $p_L$ edges coinciding.

2. Two-Step Melting (Route III):

   • Solid–Hexatic Transition: In smaller systems ($N=256$), the density response appears continuous due to strong finite-size rounding of the volume distribution. However, the LYZ analysis reveals a bifurcation between $T$ and $p_L$ edges, indicating a field-selective first-order transition. The enthalpy jump is distinct even at $N=256$, while the density jump only becomes resolvable in the $N=1024$ supercell.
   • Hexatic–Liquid Transition: This transition is not field-selective; both thermal and mechanical responses are rounded similarly. In the $N=256$ system, both appear continuous, but the $N=1024$ system reveals a clear bimodal distribution in both $H$ and $V$, establishing this as a conventional first-order transition.

Significance and Claims
The paper claims to provide a unified framework that organizes apparent disagreements among previous confined-water simulations, hard-disk models, and AgI experiments.

• Reconciling Conflicts: The study posits that previous discrepancies arose because different probes accessed different thermodynamic channels (e.g., density vs. energy) and because finite-size effects obscured weak discontinuities.
• Field-Selective Criticality: The primary contribution is the identification of "field-selective criticality," where a phase boundary does not possess a single thermodynamic character. Near specific points (like the maximum of a coexistence line or in the solid-hexatic transition), the response to temperature and lateral pressure decouples.
• Diagnostic Power: The authors argue that Lee-Yang edge bifurcation serves as a robust criterion to diagnose when thermal and mechanical responses are decoupled or when finite-size rounding is unequal across channels. This allows for a more accurate classification of transition orders than relying on ensemble averages of single observables.
• Deviation from KTHNY: The results indicate that confined water departs from the original KTHNY scenario of two continuous transitions, instead exhibiting a sequence of field-selective first-order and conventional first-order transitions.

The paper concludes that distinguishing between thermal and mechanical responses is essential in low-dimensional and confined systems, as experiments and simulations often access only a subset of observables, potentially leading to misclassification of phase transition orders.