Pressure-Temperature Phase Diagram and $\lambda$-Transition in Liquid Sulfur

Using machine-learned molecular dynamics simulations, this study provides a microscopic explanation of sulfur's $\lambda$-transition by demonstrating how temperature-induced formation of non-S$_8$ rings triggers polymerization, ultimately mapping a pressure-temperature phase diagram that reveals a critical point where polymerization merges with the melting line.

The Gist

The Shape-Shifting Mystery of Sulfur: A Story of Rings and Chains

Imagine you are looking at a massive ballroom filled with dancers. Most of them are dancing in perfect, tiny circles of eight people, holding hands tightly. This is how sulfur usually behaves: it likes to form little "rings" (specifically $S_8$ rings) that float around like tiny, stable hula hoops.

But sulfur is a bit of a rebel. If you turn up the heat or squeeze it hard enough, the "dance" changes completely. This paper uses super-advanced computer simulations to act like a high-speed camera, capturing exactly how those tiny hula hoops break apart and turn into long, tangled ropes.

Here is the breakdown of what the scientists discovered:

1. The "Lambda Transition": The Great Breakup

At normal pressure, if you heat sulfur up, it first melts into a liquid of those little $S_8$ rings. But then, you hit a specific temperature (called the $\lambda$-transition), and suddenly, everything changes.

The Analogy: Imagine a room full of people spinning in small, independent circles. Suddenly, someone turns up the music (the heat), and instead of staying in their little circles, people start grabbing hands and forming massive, winding "conga lines" that stretch across the entire room. This is polymerization. The liquid goes from being "thin and watery" (like a collection of rings) to "thick and syrupy" (like a mess of long chains).

2. The "Troublemakers" (The Secret to the Chaos)

For a long time, scientists weren't sure how the rings started breaking. Did they just snap?

The researchers found that the process isn't random. Before the big "conga line" starts, some of the $S_8$ rings start turning into "imperfect" rings—maybe circles of 7 or 9 instead of 8.

The Analogy: Think of these imperfect rings as "unstable dancers." Because they aren't perfectly balanced, they are much more likely to trip or let go of their hands. These "tripping" rings act as the spark that triggers the rest of the room to break their circles and start forming the long chains.

3. The High-Pressure Twist: Melting vs. Growing

The most exciting part of the paper is what happens when you squeeze the sulfur (increasing the pressure).

At low pressure, the sulfur melts first, and then it turns into chains. But as you squeeze it harder, the melting and the chain-forming happen almost at the same time. Even weirder, the researchers found that at very high pressures, the sulfur starts forming those long chains while it is still a solid crystal!

The Analogy: Imagine a crowd of people standing in a very strict, organized grid (the crystal). Usually, they only start dancing once they break formation and start moving (melting). But if you squeeze the room tightly, the people start grabbing hands and forming lines while they are still standing in their neat rows. The "conga line" begins to grow inside the grid before the grid even breaks apart.

4. How did they do it? (The Digital Time Machine)

Studying this in a real lab is incredibly hard because these changes happen so fast and are very sensitive to how quickly you heat things up.

To solve this, the scientists used Machine Learning. They taught a computer program the "rules" of how sulfur atoms interact using complex math (called Machine-Learned Interatomic Potentials). Once the computer "learned" the personality of sulfur, they ran a digital simulation that acted like a super-powered microscope, allowing them to watch every single atom move, break, and reconnect in real-time.

Summary

In short, this paper provides a "microscopic movie" of sulfur's transformation. It shows that sulfur isn't just a simple element; it's a master of disguise that uses "imperfect" molecules to trigger a massive structural makeover, shifting from tiny rings to giant chains, especially when under pressure.

Technical Summary

Technical Summary: Pressure–Temperature Phase Diagram and $\lambda$-Transition in Liquid Sulfur

1. Problem Statement

Elemental sulfur possesses an exceptionally complex phase diagram characterized by a competition between cyclic molecular rings (e.g., $S_8$) and linear polymeric chains. A central phenomenon is the $\lambda$-transition—a temperature-induced polymerization in the liquid phase. While experimental evidence (DSC, Raman, ESR) has long established the existence of this transition, several fundamental questions remain unanswered at the atomistic level:

• Microscopic Mechanism: What specific molecular species trigger the onset of polymerization?
• Phase Sequence: How does the system transition from crystalline $\alpha$-$S_8$ through melting to the $\lambda$-transition?
• High-Pressure Behavior: How does the polymerization line interact with the melting line, and what is the mechanism when they coincide?
• Modeling Challenges: Traditional Density Functional Theory (DFT) is too computationally expensive for the large system sizes and long timescales required to simulate polymeric liquids, while empirical potentials often fail to accurately capture bond-breaking/forming events.

2. Methodology

The authors employ large-scale Molecular Dynamics (MD) simulations driven by a high-fidelity Machine-Learned Interatomic Potential (MLIP).

• MLIP Development: They trained an Allegro model (a local equivariant message-passing architecture) using a dataset generated via spin-polarized DFT (GGA-PBE).
• Active Learning: To ensure the potential could handle diverse configurations, they used a "Query by Committee" (QBC) active learning strategy, selecting structures based on the disagreement between an ensemble of neural network potentials.
• Simulation Protocol: They performed temperature-ramp (T-ramp) simulations at various pressures (ambient to 1.0 GPa) and heating rates (5 K/ns to 500 K/ns). This allowed them to observe spontaneous melting and polymerization within a reasonable computational window.
• Structural Analysis: They utilized topological analysis (atomic coordination), SOAP (Smooth Overlap of Atomic Position) descriptors to measure crystalline similarity, and lineage tracking to follow the "fate" of specific atoms.

3. Key Contributions and Results

• Identification of Polymerization "Seeds": The study provides direct evidence that non-$S_8$ rings (specifically $S_7$ and other small cycles $S_\pi$) act as reactive centers. These species have lower dissociation energies than $S_8$ and serve as the precursors that trigger the opening of chains and subsequent catenation. This challenges the traditional Tobolsky–Eisenberg model.
• Reconstruction of the Phase Diagram: The authors successfully reconstructed the $P$-$T$ phase diagram. They observed that the polymerization temperature ($T_\lambda$) decreases moderately with pressure until it merges with the melting line at a critical point.
• Mechanism of High-Pressure Transition: At pressures $\geq 0.5$ GPa, where melting and polymerization are coupled, the study reveals that polymerization actually initiates within the crystalline lattice before melting occurs. They observed non-$S_8$ rings and polymeric chains "conforming" to the lattice sites (a process akin to topochemical polymerization).
• Validation of Experimental Signatures:
  • Heat Capacity: The simulations reproduced the characteristic $\lambda$-shaped peak in heat capacity ($C_s$), aligning more closely with non-mean-field theoretical predictions.
  • Heating Rate Dependence: They captured the strong dependence of $T_\lambda$ on the heating rate, allowing them to extrapolate a $T_\lambda$ value ($\approx 500$ K) close to experimental observations.
  • Physical Properties: The simulations accurately reflected changes in potential energy, volume, and diffusivity associated with the transition.

4. Significance

This work represents a significant advancement in computational materials science by providing the first complete atomistic picture of the sulfur phase diagram from the solid state through the $\lambda$-transition. By bridging the gap between quantum mechanical accuracy (DFT) and large-scale classical MD, the authors have:

1. Resolved long-standing debates regarding the microscopic triggers of sulfur polymerization.
2. Demonstrated the utility of equivariant MLIPs in simulating complex, reactive, and phase-changing systems.
3. Provided a theoretical framework for understanding how structural transformations can be "pre-initiated" within a crystal lattice under pressure, offering insights applicable to other complex elemental systems.