Ferroelectric quantum critical point in superconducting hydrides: The case of H$_3$S

By employing path integral molecular dynamics with a machine-learned potential, this study reveals that the high-temperature superconductivity in H$_3$S occurs in a paraelectric region dominated by large nuclear quantum fluctuations above a ferroelectric quantum critical point at approximately 134 GPa, which belongs to the 4D Ising universality class.

The Gist

Imagine you have a tiny, super-dense ball of sulfur and hydrogen atoms. Under extreme pressure, this ball becomes a superconductor—a material that conducts electricity with zero resistance. For a long time, scientists were puzzled by why this happens in a specific material called H3S. They knew it worked best at a certain pressure (about 155 GPa), but the map of how the atoms behave was missing.

This paper is like drawing that missing map. The researchers used a powerful computer simulation to track how the atoms dance, not just as solid balls, but as "fuzzy clouds" of probability (a quantum effect). Here is what they found, explained simply:

1. The "Fuzzy" Atoms and the Magic Pressure

In the world of tiny atoms, things aren't solid; they wiggle and jitter. The researchers found that at a specific pressure (around 134 GPa), the hydrogen atoms in H3S hit a "tipping point."

• The Analogy: Imagine a ball sitting in a bowl. If the bowl is deep, the ball stays in the middle. If you shake the bowl (heat) or squeeze it (pressure), the ball might start rolling around.
• The Discovery: At this tipping point, called a Quantum Critical Point (QCP), the atoms are in a state of maximum confusion. They aren't settled in one spot, but they aren't totally random either. They are "fluctuating" wildly, like a crowd of people trying to decide which way to turn.

2. The Phase Change: From "Symmetrical" to "Lopsided"

The material can exist in two main shapes (phases):

• The "Perfectly Balanced" Phase (Paraelectric): The hydrogen atoms sit right in the middle between sulfur atoms. It's symmetrical, like a perfectly balanced seesaw.
• The "Lopsided" Phase (Ferroelectric): The hydrogen atoms get pushed to one side. The seesaw tips over.

The paper shows that the transition from "balanced" to "lopsided" doesn't happen exactly where the superconductivity is strongest. Instead, the superconductivity peak happens in the "balanced" zone, but right next to the tipping point where the atoms are wiggling the most.

3. The "Superconducting Sweet Spot"

Here is the big surprise:

• Old Idea: Scientists thought the superconductivity peak happened because the material was switching from balanced to lopsided.
• New Finding: The paper shows the peak actually happens in the balanced zone, but right next to the chaos.
• The Analogy: Think of a surfer. The best waves aren't the calm, flat water, nor are they the chaotic, crashing surf. The best waves are right where the ocean is starting to get rough. The "roughness" (quantum fluctuations) helps the electrons pair up and flow without resistance. The paper suggests that the wild wiggling of the hydrogen atoms near the tipping point acts like a boost for the superconductivity.

4. The "4D Ising" Rulebook

The researchers analyzed the math behind this tipping point and found it follows a very specific rulebook known as the 4D Ising universality class.

• The Analogy: Imagine different games (like chess, checkers, or Go). Even though they look different, they might all follow the same underlying logic of how pieces move. The researchers found that the way these atoms behave follows the same "logic" as a specific, complex mathematical model used to describe how things change state in four dimensions. This confirms that their discovery is a fundamental law of physics, not just a fluke.

5. Why This Matters for the "Map"

Before this study, the map of H3S was blurry. The researchers used a new type of "AI brain" (a machine learning potential) to run simulations that were too expensive to do with old methods.

• They found that if you ignore the quantum "fuzziness" of the atoms (treating them like solid billiard balls), you get the wrong map. The quantum wiggles shift the transition pressure by a huge amount (about 50 GPa).
• By including these wiggles, they finally located the "tipping point" (QCP) and showed that the superconducting peak sits in a region of strong quantum fluctuations, just above the tipping point.

Summary

The paper reveals that the superconducting magic in H3S isn't caused by the material simply changing shape. Instead, it happens because the material is hovering right next to a "quantum tipping point" where the atoms are vibrating wildly. These wild vibrations act like a catalyst, helping the electricity flow perfectly. The researchers have now mapped exactly where this happens and proved it follows a specific, universal mathematical rule.

Technical Summary

Problem Statement
Sulfur hydride (H$_3$S) is a prominent high-temperature superconductor with a critical temperature ($T_c$) of 203 K at approximately 155 GPa. While the superconducting dome exhibits a peak similar to cuprates and doped SrTiO$_3$, the precise nature of the phase diagram, particularly the relationship between the superconducting maximum and the structural phase transition, remains unresolved. Previous theoretical studies suggested a ferroelectric phase transition from the body-centered cubic $Im\bar{3}m$ phase to the trigonal $R3m$ phase, driven by hydrogen displacement. However, these theories placed the transition pressure significantly lower (by over 40 GPa) than the experimental $T_c$ peak. Furthermore, a comprehensive understanding of the H$_3$S phase diagram across an extended temperature and pressure range, including the interplay of thermal and nuclear quantum effects (NQE), has been lacking.

Methodology
To address these gaps, the authors employed Path Integral Molecular Dynamics (PIMD) to treat nuclear quantum effects exactly. Given the prohibitive computational cost of combining PIMD with ab initio electronic structure methods (such as Density Functional Theory, DFT) for the required system sizes and temperatures, the authors utilized a Machine Learning Interatomic Potential (MLIP). Specifically, they trained a MACE (Higher Order Equivariant Message Passing Neural Network) potential on a dataset of DFT-BLYP configurations spanning pressures from 40 to 220 GPa at 100 K and 200 K.

The simulations were conducted using the i-PI package with NVT ensembles for system sizes of $L \times L \times L$ primitive cells ($L \in \{2, 3, 4\}$) to enable finite-size scaling analysis. The study covered temperatures from 50 K to 300 K. Key observables included:

• Order Parameters: Global proton displacement ($\Delta$) and its absolute value ($\Delta_{abs}$) to distinguish between paraelectric and ferroelectric phases.
• Local Observables: Distribution of local proton displacements and local dipole moments.
• Dynamical Properties: Imaginary-time phonon Green's functions ($G_{ij}(\tau)$) to extract anharmonic phonon frequencies and study retardation effects.
• Scaling Analysis: Finite-size scaling of the order parameter near the critical point to determine the universality class.

Key Results

1. Phase Diagram and Quantum Critical Point (QCP): The study reveals a ferroelectric transition line that shifts to lower pressures as temperature decreases. By extrapolating to the thermodynamic limit and $T=0$ K, the authors identify a ferroelectric Quantum Critical Point (QCP) at $p_{QCP} \approx 134 \pm 2$ GPa.
2. Universality Class: Finite-size scaling analysis of the order parameter near the QCP yields a critical exponent $\nu \approx 0.483 \pm 0.014$. This value is consistent with the 4D Ising universality class (where $\nu = 1/2$), assuming Lorentz invariance ($z=1$).
3. Discrepancy with Experimental $T_c$: The calculated ferroelectric transition line (e.g., at 124 GPa for 200 K) is significantly lower than the experimental superconducting peak at 155 GPa. Consequently, the experimental $T_c$ maximum falls within the paraelectric $Im\bar{3}m$ phase, not the ferroelectric phase.
4. Role of Quantum Fluctuations: The region surrounding the QCP and extending into the paraelectric phase is characterized by large nuclear quantum fluctuations. The authors quantify this using the ratio of the variance of the local Green's function at imaginary time $\beta/2$ to that at $\tau=0$ ($\sigma^2_g(\beta/2)/\sigma^2_g(0)$). This ratio peaks in the paraelectric phase near the QCP, indicating strong retardation effects and fluctuating dipole moments, which are suppressed in the ferroelectric phase where protons freeze.
5. Comparison with Classical Dynamics: Classical Molecular Dynamics (MD) simulations show that the ferroelectric transition occurs at much higher pressures (shifted by $\approx 50$ GPa at 200 K) compared to the quantum case, highlighting the critical role of NQE in stabilizing the paraelectric phase and lowering the transition pressure.
6. Phonon Softening: Optical phonon modes driving the transition exhibit softening as the system approaches the transition from the paraelectric side, followed by a sharp frequency jump upon entering the ferroelectric phase, consistent with a second-order phase transition.

Significance and Claims
The paper establishes that the high-$T_c$ superconductivity in H$_3$S does not occur at the ferroelectric phase transition itself, but rather in the paraelectric region dominated by strong quantum fluctuations and retardation effects near the QCP. The authors suggest that these fluctuations, potentially mediated by the proximity to the ferroelectric QCP, may enhance superconductivity. This scenario implies that standard Migdal-Eliashberg theory, which often neglects full frequency dependence and vertex corrections, may be insufficient for describing H$_3$S. The presence of a ferroelectric QCP and strong optical phonon softening supports the possibility of electron pairing mechanisms beyond linear order, similar to those proposed for dilute metals like SrTiO$_3$. The work provides a detailed, quantum-mechanically accurate phase diagram that reconciles the location of the superconducting dome with the underlying structural instabilities, emphasizing the necessity of including nuclear quantum effects and non-perturbative electron-phonon coupling evaluations in future theoretical treatments of superconducting hydrides.