🔬 condensed matter

Anharmonicity and Nonadiabaticity in Hydride Superconductors

This paper investigates superconductivity in various hydrides using advanced theoretical frameworks that incorporate anharmonicity and nonadiabaticity, introducing specific diagnostic metrics to determine when corrections beyond standard Migdal-Eliashberg theory are necessary to accurately predict critical temperatures.

Original authors: Shashi B. Mishra, Francesco Belli, Eva Zurek, Elena R. Margine

Published 2026-08-18

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Shashi B. Mishra, Francesco Belli, Eva Zurek, Elena R. Margine

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

In the quest to understand how electricity can flow without resistance, scientists have long looked to materials called superconductors. These are special substances that, when cooled to very low temperatures, allow electric current to move through them without losing any energy as heat. For decades, the challenge has been finding materials that do this at temperatures we can easily reach, rather than requiring extreme cold. A major breakthrough in recent years came with the discovery that certain hydrogen-rich compounds, when squeezed under immense pressure, become superconductors at temperatures far warmer than previously thought possible, some even approaching the warmth of a summer day. To predict which materials will work and why, researchers rely on a set of mathematical tools that describe how electrons interact with the vibrating atoms of a crystal. These tools assume that the atoms vibrate in a simple, predictable way and that the electrons move much faster than the atoms can react. However, in materials packed with light hydrogen atoms, these assumptions can break down. The atoms may vibrate wildly and unpredictably, and the electrons may not be able to keep up, forcing scientists to look for more complex ways to describe the physics.

A team of researchers set out to test these complex effects in a variety of hydrogen-based superconductors to see where the standard tools succeed and where they fail. They focused on several specific materials, including hydrogen sulfide under high pressure, yttrium hydrides, a compound containing lanthanum, beryllium, and hydrogen, and a series of palladium hydrides at normal pressure. Using powerful computer simulations, they calculated how the atoms in these materials actually move, accounting for the fact that light atoms like hydrogen do not just vibrate in a simple back-and-forth motion but can explore a wider, more chaotic range of positions. They also checked whether the electrons and atoms were interacting in a way that required a more advanced description than the standard models allowed. By comparing their detailed calculations with real-world experimental data, the researchers aimed to build a reliable guide for predicting the temperature at which these materials become superconductors.

The study revealed that the behavior of these materials falls into three distinct categories, depending on how the atoms move and how they interact with electrons. For the high-pressure materials like hydrogen sulfide and the yttrium hydrides, the researchers found that both the wild vibrations of the atoms and the complex timing of the electron-atom interactions were crucial. When they ignored these effects, their predictions for the superconducting temperature were too high. However, once they included the corrections for the atoms' chaotic motion and the delayed reaction of the electrons, their calculations matched the experimental measurements almost perfectly. For instance, in hydrogen sulfide at 200 gigapascals of pressure, the standard model predicted a superconducting temperature of 244 kelvin, but the refined calculation brought this down to 184 kelvin, which is exactly what experiments had observed. Similarly, for the yttrium hydrides, the refined approach brought the predicted temperatures into close alignment with what was measured in the lab.

In contrast, the palladium hydrides, which exist at normal atmospheric pressure, told a different story. In these materials, the atoms still vibrate in a complex, anharmonic way that significantly changes the superconducting temperature, but the electrons move fast enough that the standard timing assumptions still hold true. Here, the researchers found that accounting for the atoms' wild vibrations was enough to explain the experimental results, including a curious phenomenon where heavier isotopes of hydrogen actually led to higher superconducting temperatures. The complex electron-atom timing corrections were so small in these materials that they could be safely ignored. This distinction helped the team create a clear map of which materials need which level of theoretical detail. They introduced two simple measures to classify any hydride: one that tracks how much the atoms' wild motion changes the material's properties, and another that measures how much the electron-atom timing matters.

One material, however, remained an outlier. The compound containing lanthanum, beryllium, and hydrogen showed very little of the wild atomic motion seen in the other materials, and the timing corrections were also small. Even after applying all the advanced corrections, the researchers' calculations still predicted a superconducting temperature of about 159 kelvin, which is significantly higher than the 110 kelvin observed in experiments. This suggests that while the team's methods work well for most of the materials they studied, there is still something missing in the description of this specific compound. It is possible that even more subtle effects, not yet included in their models, are at play.

The work provides a practical framework for scientists to decide how much detail is needed when studying a new superconducting material. By using the two measures they developed, researchers can quickly determine if a material requires the full, complex treatment or if a simpler approach will suffice. This clarity is essential as the field moves toward designing new materials that could one day operate at room temperature without the need for extreme pressure. The study confirms that for the most promising high-pressure hydrides, the path to accurate prediction requires acknowledging both the chaotic dance of the atoms and the subtle delays in how electrons respond to them. While the mystery of the lanthanum-beryllium-hydrogen compound remains, the new tools offer a reliable way forward for understanding the vast family of hydrogen-rich superconductors.

Technical Summary: Anharmonicity and Nonadiabaticity in Hydride Superconductors

Problem Statement
Hydrogen-rich compounds under high pressure represent the most prominent class of high-temperature superconductors, with critical temperatures ( $T_c$ ) exceeding 200 K in materials like H $_3$ S and YH $_6$ , and near-room-temperature values in LaH $_{10}$ . While the discovery and rationalization of these materials have been guided by Migdal-Eliashberg (ME) theory, quantitative agreement between calculated and experimental $T_c$ values is not uniform. In several benchmarks, state-of-the-art harmonic and adiabatic ME calculations overestimate $T_c$ .

Two primary limitations of conventional ME theory are identified as potential sources of these discrepancies:

Harmonic Approximation: Conventional theory neglects quantum anharmonic lattice dynamics, which are essential when light hydrogen atoms undergo large zero-point motion. This motion strongly renormalizes phonon spectra and electron-phonon (e-ph) coupling.
Adiabatic Approximation (Migdal's Theorem): Conventional theory neglects e-ph vertex corrections, justified only when phonon energies are small compared to the electronic energy scale. In hydrides, hydrogen-derived phonons can reach 100–250 meV, becoming comparable to the energy scale of electronic structure variations near the Fermi level ( $\varepsilon_F$ ), particularly in the presence of van Hove singularities.

The central question addressed is determining which approximation (harmonic vs. anharmonic, adiabatic vs. nonadiabatic) limits the accuracy of ME predictions for specific hydrides and whether both corrections are simultaneously required.

Methodology
The authors employ a first-principles framework combining several advanced computational techniques:

Anharmonic Phonons: Calculated using the Stochastic Self-Consistent Harmonic Approximation (SSCHA) accelerated by bespoke machine learning interatomic potentials (ML/SSCHA). This approach accounts for quantum nuclear effects and renormalizes phonon frequencies and eigenvectors.
Electron-Phonon Coupling: Matrix elements are interpolated in the Wannier representation using the EPW code.
Superconducting Properties: The isotropic Eliashberg equations are solved using Full-Bandwidth (FBW) theory, which retains the energy dependence of the electronic density of states (DOS) rather than approximating it at the Fermi level.
Vertex Corrections: First-order e-ph vertex corrections are included to account for nonadiabatic effects, modifying the pairing kernel.
Materials Studied: The framework is applied to a representative set of hydrides spanning different bonding environments: H $_3$ S (covalent S–H network), YH $_6$ and YH $_9$ (clathrate-like H networks), LaBeH $_8$ (ternary hydride with BeH $_8$ units), and the ambient-pressure PdH/PdD/PdT series.

Key Contributions and Descriptors
The paper introduces two material-specific, dimensionless descriptors to classify hydrides based on the necessity of corrections beyond standard ME theory:

Anharmonic Renormalization ( $A_\lambda$ ): Defined as $|\lambda_{anh} - \lambda_{har}| / \lambda_{har}$ , this quantifies the relative change in e-ph coupling induced by anharmonic lattice dynamics.
Vertex Ratio ( $R_V$ ): Defined as $|\lambda_V^{anh}| / \lambda_{anh}$ , this quantifies the magnitude of the first-order vertex contribution relative to the anharmonic coupling.

These descriptors allow for the classification of materials into distinct regimes regarding the breakdown of the harmonic and adiabatic approximations.

Results
The study reveals three distinct regimes among the hydrides investigated:

Strong Anharmonicity and Nonadiabaticity (H $_3$ S, YH $_6$ , YH $_9$ ):
- These high-pressure hydrides require both anharmonic and nonadiabatic corrections.
- Anharmonicity hardens hydrogen-derived optical branches, suppressing the e-ph coupling ( $\lambda$ ) by 25–40% and lowering $T_c$ .
- Vertex corrections further reduce $T_c$ . For example, in H $_3$ S at 200 GPa, the harmonic $T_c$ (244 K) is reduced to 196 K by anharmonicity and further to 184 K by vertex corrections, matching the experimental value of 184 K.
- In YH $_6$ and YH $_9$ , the combined corrections bring calculated $T_c$ values (226 K and 231 K, respectively) into close agreement with experimental measurements (224 K and 239 K).
Strong Anharmonicity but Adiabatic (PdH, PdD, PdT):
- The ambient-pressure Pd hydrides exhibit strong anharmonicity, which removes harmonic instabilities and hardens optical rattling modes.
- However, they remain adiabatic due to low phonon energies and a smooth DOS near $\varepsilon_F$ , resulting in negligible vertex corrections ( $R_V \approx 0.08$ ).
- The anharmonic treatment alone reproduces the inverse isotope effect (where $T_c$ increases with isotope mass: PdH < PdD < PdT), yielding $T_c$ values of 8.0, 10.5, and 11 K, respectively, in excellent agreement with experiment.
Weak Anharmonicity with Finite Vertex Correction (LaBeH $_8$ ):
- LaBeH $_8$ exhibits weak anharmonicity ( $A_\lambda \approx 0.02$ ), with only a slight reduction in $\lambda$ .
- Vertex corrections are finite ( $R_V \approx 0.11$ ) and reduce $T_c$ by approximately 4 K (from 163 K to 159 K).
- Despite these corrections, the predicted $T_c$ remains substantially higher than the experimental value of $\approx$ 110 K at 80 GPa. The authors note that nonlinear e-ph corrections or anisotropic treatments are not obvious remedies as they might further increase $T_c$ .

Significance
The paper establishes a first-principles framework for identifying the hierarchy of corrections required for quantitative predictions of $T_c$ in hydride superconductors. By introducing the descriptors $A_\lambda$ and $R_V$ , the authors provide a practical diagnostic tool to determine when the harmonic and adiabatic approximations break down. The framework successfully reproduces experimental trends for H $_3$ S, YH $_6$ , YH $_9$ , and the Pd hydrides, demonstrating that the simultaneous inclusion of anharmonicity and nonadiabatic vertex corrections is crucial for high-pressure hydrides with sharp electronic features. The persistent discrepancy in LaBeH $_8$ suggests that for certain materials, even these advanced corrections may be insufficient, pointing to the need for further investigation into higher-order effects or alternative mechanisms.

Technical Summary: Anharmonicity and Nonadiabaticity in Hydride Superconductors

More like this