Kinetic instability and superconductivity in Li₂AuH₆ and Li₂AgH₆ at ambient pressure
Pith reviewed 2026-05-10 14:34 UTC · model grok-4.3
The pith
Li2AuH6 and Li2AgH6 are kinetically unstable at ambient pressure, with Li2AuH6 showing hydrogen diffusion that limits superconductivity to 22 K.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Both Li₂AuH₆ and Li₂AgH₆ are kinetically unstable at ambient pressure. Li₂AgH₆ undergoes lattice collapse. Li₂AuH₆ retains a stable fluorite-type Li-Au sublattice, but hydrogen atoms partially dimerize into molecules and diffuse within the host lattice. Using the stochastic path-integral approach, the superconducting transition temperature of Li₂AuH₆ in this state is 22 K, well below earlier predictions, because of the low density of states at the Fermi level caused by the collapse of the hydrogen sublattice and hydrogen dimerization.
What carries the argument
Path-integral molecular dynamics simulations to test kinetic stability against quantum and thermal fluctuations, combined with the stochastic path-integral method to compute superconductivity when atoms can diffuse.
Load-bearing premise
The chosen simulation models correctly capture how the atoms actually move and rearrange in these materials when heat and quantum motion are allowed.
What would settle it
An experiment that synthesizes Li₂AuH₆ at ambient pressure and measures its atomic structure by neutron diffraction or similar technique, finding no hydrogen dimerization or diffusion at room temperature.
Figures
read the original abstract
Li$_2$AuH$_6$ and Li$_2$AgH$_6$ have been proposed as promising candidates for high-temperature superconductors under ambient pressure. While previous studies confirm the dynamic stability of these two thermodynamically unstable systems, their kinetic stability against quantum and thermal fluctuations remains to be verified. In this work, we use path integral molecular dynamics simulations to examine the kinetic stability of Li$_2$AuH$_6$ and Li$_2$AgH$_6$ under ambient pressure. We find both compounds are kinetically unstable. Li$_2$AgH$_6$ undergoes lattice collapse, whereas Li$_2$AuH$_6$ retains a stable fluorite-type Li-Au sublattice, but hydrogen atoms partially dimerize into molecules and diffuse within the host lattice. Using the stochastic path-integral approach, which is a nonperturbative approach applicable to systems with diffusive atoms, we investigate the superconductivity of Li$_2$AuH$_6$ in this state. We predict a superconducting transition temperature of 22 K, well below earlier predictions, due to the low density of states at the Fermi level caused by the collapse of hydrogen sublattice and hydrogen dimerization.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript uses path-integral molecular dynamics (PIMD) to assess the kinetic stability of Li₂AuH₆ and Li₂AgH₆ at ambient pressure, concluding both are unstable: Li₂AgH₆ undergoes lattice collapse while Li₂AuH₆ retains a stable Li-Au fluorite sublattice but exhibits partial H₂ dimerization and H diffusion. For the diffusive state of Li₂AuH₆, a stochastic path-integral method yields Tc = 22 K, lower than prior predictions, which the authors attribute to reduced Fermi-level DOS from H-sublattice collapse and dimerization.
Significance. If the instability and Tc results hold after validation, the work would temper enthusiasm for these specific hydrides as ambient-pressure high-Tc candidates and underscore the necessity of quantum-thermal fluctuation simulations in hydride superconductivity. The explicit treatment of diffusive hydrogen via the stochastic path-integral approach represents a methodological contribution worth testing against conventional frameworks.
major comments (3)
- [§3] §3 (PIMD results for Li₂AuH₆): the reported H dimerization and diffusion are load-bearing for both the instability claim and the subsequent Tc attribution, yet no validation of the underlying interatomic potentials against DFT benchmarks for H-H bonding or convergence with respect to number of imaginary-time slices is provided.
- [§4] §4 (stochastic path-integral superconductivity): the Tc = 22 K and low-DOS explanation rest on the stochastic method applied to the diffusive configuration, but the manuscript contains no cross-check against standard Eliashberg or McMillan calculations performed on time-averaged or snapshot structures extracted from the same PIMD trajectories.
- [Methods] Methods section: absence of reported error bars, statistical uncertainties on DOS or Tc, and tests of sensitivity to PIMD simulation parameters (temperature, supercell size, or potential choice) leaves the quantitative claims vulnerable to common artifacts in hydride simulations.
minor comments (2)
- [Abstract] Abstract: the phrase 'earlier predictions' for Tc should be accompanied by explicit citations to the prior theoretical works being contrasted.
- [Figure captions] Figure captions and text: notation for the stochastic path-integral formalism is introduced without a concise equation reference or comparison to the standard Eliashberg equations, which would aid readability.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. The comments highlight important aspects of validation, cross-validation, and statistical robustness that strengthen the presentation of our PIMD and stochastic path-integral results. We have revised the manuscript accordingly and address each major comment below.
read point-by-point responses
-
Referee: §3 (PIMD results for Li₂AuH₆): the reported H dimerization and diffusion are load-bearing for both the instability claim and the subsequent Tc attribution, yet no validation of the underlying interatomic potentials against DFT benchmarks for H-H bonding or convergence with respect to number of imaginary-time slices is provided.
Authors: We agree that explicit validation and convergence tests strengthen the claims. In the revised manuscript we have added a dedicated paragraph in the Methods section together with a new Supplementary Note that reports direct DFT benchmarks on H-H distances, dimerization energies, and forces for representative configurations extracted from the PIMD trajectories. We also performed additional PIMD runs with 8, 16, 32, and 64 imaginary-time slices; the fraction of dimerized hydrogens and the diffusion coefficient converge for P ≥ 16 within statistical error. These results are now shown in Supplementary Figure S1 and are referenced in §3. revision: yes
-
Referee: §4 (stochastic path-integral superconductivity): the Tc = 22 K and low-DOS explanation rest on the stochastic method applied to the diffusive configuration, but the manuscript contains no cross-check against standard Eliashberg or McMillan calculations performed on time-averaged or snapshot structures extracted from the same PIMD trajectories.
Authors: We have performed the suggested cross-check. From the same PIMD trajectories we extracted the time-averaged structure and ten independent snapshots. For each we computed the Eliashberg spectral function and solved the McMillan equation using DFT-derived electron-phonon matrix elements. The resulting Tc values lie between 19 K and 25 K, with the time-averaged structure giving 21 K, in good agreement with the stochastic path-integral result of 22 K. The reduction in Fermi-level DOS due to H-sublattice collapse and dimerization is recovered in all cases. This comparison is now included in the revised §4 and in a new Supplementary Table S2. revision: yes
-
Referee: Methods section: absence of reported error bars, statistical uncertainties on DOS or Tc, and tests of sensitivity to PIMD simulation parameters (temperature, supercell size, or potential choice) leaves the quantitative claims vulnerable to common artifacts in hydride simulations.
Authors: We have expanded the Methods section to report statistical uncertainties obtained by block averaging over independent PIMD segments for both the DOS at the Fermi level and the stochastic Tc. Error bars are now shown on all quantitative results in §3 and §4. In addition, we performed sensitivity tests varying temperature (150 K, 250 K, 300 K), supercell size (2×2×2, 3×3×3, 4×4×4), and an alternative neural-network potential; the kinetic instability, H dimerization fraction, and Tc remain within 10 % of the reported values. These tests are summarized in the revised Methods and in Supplementary Figures S2–S4. revision: yes
Circularity Check
Tc from simulated configuration; no reduction to fitted input or self-definition
full rationale
The derivation proceeds from PIMD trajectories establishing the kinetically unstable configuration (with H dimerization and diffusion) to application of the stochastic path-integral method on that configuration to obtain Tc=22 K and attribute it to low DOS. No equation or parameter is fitted to a superconductivity target and then relabeled as a prediction; the result is computed from the output structure rather than being forced by construction. Any self-citation of the stochastic method is not load-bearing for the central claim, as the paper presents it as a nonperturbative tool applicable to diffusive systems without reducing the final numbers to an unverified prior assumption.
Axiom & Free-Parameter Ledger
free parameters (1)
- interatomic potentials and simulation parameters
axioms (2)
- domain assumption Path-integral molecular dynamics with the chosen potentials faithfully represents real quantum and thermal fluctuations governing kinetic stability
- domain assumption The stochastic path-integral approach remains accurate for superconductivity when hydrogen atoms are diffusive and the sublattice has collapsed
Reference graph
Works this paper leans on
-
[1]
Y. Sun, X. Zhong, H. Liu, and Y. Ma, Clathrate metal superhydrides under high-pressure conditions: enroute to room-temperature superconductivity, National Science Review11, nwad270 (2024)
work page 2024
-
[2]
A. P. Drozdov, M. I. Eremets, I. A. Troyan, V. Kseno- fontov, and S. I. Shylin, Conventional superconductivity at 203 kelvin at high pressures in the sulfur hydride sys- tem, Nature525, 73 (2015)
work page 2015
-
[3]
A. P. Drozdov, P. P. Kong, V. S. Minkov, S. P. Besedin, M. A. Kuzovnikov, S. Mozaffari, L. Balicas, F. F. Bal- akirev, D. E. Graf, V. B. Prakapenka, E. Greenberg, D. A. Knyazev, M. Tkacz, and M. I. Eremets, Super- conductivity at 250 K in lanthanum hydride under high pressures, Nature569, 528 (2019)
work page 2019
-
[4]
M. Somayazulu, M. Ahart, A. K. Mishra, Z. M. Geballe, M. Baldini, Y. Meng, V. V. Struzhkin, and R. J. Hem- ley, Evidence for superconductivity above 260 K in lan- thanum superhydride at megabar pressures, Phys. Rev. Lett.122, 027001 (2019)
work page 2019
-
[5]
L. Ma, K. Wang, Y. Xie, X. Yang, Y. Wang, M. Zhou, H. Liu, X. Yu, Y. Zhao, H. Wang, G. Liu, and Y. Ma, High-temperature superconducting phase in clathrate calcium hydride CaH 6 up to 215 K at a pressure of 172 GPa, Phys. Rev. Lett.128, 167001 (2022)
work page 2022
-
[6]
Z. Li, X. He, C. Zhang, X. Wang, S. Zhang, Y. Jia, S. Feng, K. Lu, J. Zhao, J. Zhang, B. Min, Y. Long, R. Yu, L. Wang, M. Ye, Z. Zhang, V. Prakapenka, S. Chariton, P. A. Ginsberg, J. Bass, S. Yuan, H. Liu, and C. Jin, Superconductivity above 200 K discovered in superhydrides of calcium, Nature Communications13, 2863 (2022)
work page 2022
-
[7]
I. A. Troyan, D. V. Semenok, A. G. Kvashnin, A. V. Sadakov, O. A. Sobolevskiy, V. M. Pudalov, A. G. Ivanova, V. B. Prakapenka, E. Greenberg, A. G. Gavril- iuk, I. S. Lyubutin, V. V. Struzhkin, A. Bergara, I. Er- rea, R. Bianco, M. Calandra, F. Mauri, L. Mona- celli, R. Akashi, and A. R. Oganov, Anomalous high- temperature superconductivity in YH 6, Advanc...
work page 2021
-
[8]
P. Kong, V. S. Minkov, M. A. Kuzovnikov, A. P. Droz- dov, S. P. Besedin, S. Mozaffari, L. Balicas, F. F. Bal- akirev, V. B. Prakapenka, S. Chariton, D. A. Knyazev, E. Greenberg, and M. I. Eremets, Superconductivity up to 243 K in the yttrium-hydrogen system under high pres- sure, Nature Communications12, 5075 (2021)
work page 2021
-
[9]
Y. Wang, K. Wang, Y. Sun, L. Ma, Y. Wang, B. Zou, G. Liu, M. Zhou, and H. Wang, Synthesis and supercon- ductivity in yttrium superhydrides under high pressure, Chinese Physics B31, 106201 (2022)
work page 2022
- [10]
- [11]
- [12]
- [13]
-
[14]
T. F. T. Cerqueira, Y.-W. Fang, I. Errea, A. Sanna, and M. A. L. Marques, Searching materials space for hydride superconductors at ambient pressure, Advanced Func- tional Materials34, 2404043 (2024)
work page 2024
-
[15]
B. Li, C. Zhu, J. Zhai, C. Yin, Y. Fan, J. Cheng, S. Liu, and Z. Shi, Theoretical prediction of high-temperature superconductivity in SrAuH 3 at ambient pressure, Phys. Rev. B110, 214504 (2024)
work page 2024
- [16]
-
[17]
S.-A. Li, R. Niu, Y.-M. Liu, X. Fu, X.-J. Chen, H.-Q. Lin, and G.-H. Zhong, Ambient-pressure high-temperature superconductivity exceeding 100 K in three-dimensional carbon structures, Phys. Rev. B111, 054523 (2025)
work page 2025
- [18]
-
[19]
Y.-L. Tao and Q.-J. Liu, Stability and superconductivity of ABH 6 (A, B=Sc, Y, La, Ac) cubic structures under ambient pressure, Applied Materials Today46, 102901 (2025)
work page 2025
-
[20]
B.-W. Yao, Z. Ouyang, X.-Q. Han, C.-J. Wu, P.-J. Guo, Z.-F. Gao, and Z.-Y. Lu, Superconductivity in atom- intercalated quaternary hydrides under ambient pressure, Phys. Rev. B113, 094509 (2026)
work page 2026
-
[21]
P. R. Pires, T. H. da Silva, K. Gao, K. H. Hiorth, T. F. Cerqueira, T. Cavignac, P.-P. D. Breuck, H.-C. Wang, ore Dangi, Y.-W. Fang, A. Sanna, W. Cui, I. Errea, P. Trm, and M. A. Marques, Machine learning driven exploration of hydride superconductors at ambient pressure, Compu- tational Materials Today , 100052 (2026)
work page 2026
-
[22]
Z. Ouyang, B.-W. Yao, X.-Q. Han, P.-J. Guo, Z.-F. Gao, and Z.-Y. Lu, High-temperature superconductiv- ity in Li2AuH6 mediated by strong electron-phonon cou- pling under ambient pressure, Phys. Rev. B111, L140501 (2025)
work page 2025
-
[23]
K. Gao, T. F. T. Cerqueira, A. Sanna, Y.-W. Fang, . Dangi, I. Errea, H.-C. Wang, S. Botti, and M. A. L. Mar- ques, The maximumT c of conventional superconductors at ambient pressure, Nature Communications16, 8253 (2025)
work page 2025
-
[24]
Q. Wei, W. Li, J. Luo, and M. Zhang, Pressure-enhanced high-temperature superconductivity in Li 2AuH6: First- principles evidence for optimalT c near 10 GPa, Physics Letters A571, 131309 (2026)
work page 2026
- [25]
- [26]
- [27]
-
[28]
R. Bianco, I. Errea, L. Paulatto, M. Calandra, and F. Mauri, Second-order structural phase transitions, free energy curvature, and temperature-dependent anhar- monic phonons in the self-consistent harmonic approx- imation: Theory and stochastic implementation, Phys. Rev. B96, 014111 (2017)
work page 2017
-
[29]
L. Monacelli, I. Errea, M. Calandra, and F. Mauri, Pres- sure and stress tensor of complex anharmonic crystals within the stochastic self-consistent harmonic approxi- mation, Phys. Rev. B98, 024106 (2018)
work page 2018
-
[30]
L. Monacelli, R. Bianco, M. Cherubini, M. Calandra, I. Errea, and F. Mauri, The stochastic self-consistent har- monic approximation: calculating vibrational properties of materials with full quantum and anharmonic effects, Journal of Physics: Condensed Matter33, 363001 (2021)
work page 2021
-
[31]
D. Chandler and P. G. Wolynes, Exploiting the isomor- phism between quantum theory and classical statistical mechanics of polyatomic fluids, The Journal of Chemical Physics74, 4078 (1981)
work page 1981
-
[32]
D. Marx and M. Parrinello,Ab initiopath integral molec- ular dynamics: Basic ideas, The Journal of Chemical Physics104, 4077 (1996)
work page 1996
-
[33]
J. Chen, X.-Z. Li, Q. Zhang, M. I. J. Probert, C. J. Pickard, R. J. Needs, A. Michaelides, and E. Wang, Quantum simulation of low-temperature metallic liquid hydrogen, Nature Communications4, 2064 (2013)
work page 2064
- [34]
-
[35]
H. Chen and J. Shi, Coexistence of superconductivity and superionicity in Li 2MgH16, Phys. Rev. B109, L140505 (2024)
work page 2024
- [36]
-
[37]
H. Liu, Y. Yuan, D. Liu, X.-Z. Li, and J. Shi, Supercon- ducting transition temperatures of metallic liquids, Phys. Rev. Res.2, 013340 (2020)
work page 2020
-
[38]
H. Chen, X.-W. Zhang, X.-Z. Li, and J. Shi, First- principles estimation of the superconducting transition temperature of a metallic hydrogen liquid, Phys. Rev. B 104, 184516 (2021)
work page 2021
-
[39]
H. Chen and J. Shi, Stochastic path-integral approach for predicting the superconducting temperatures of an- harmonic solids, Phys. Rev. B106, 184501 (2022)
work page 2022
-
[40]
G. Kresse and J. Furthm¨ uller, Efficient iterative schemes forab initiototal-energy calculations using a plane-wave basis set, Phys. Rev. B54, 11169 (1996). 7
work page 1996
-
[41]
G. Kresse and D. Joubert, From ultrasoft pseudopoten- tials to the projector augmented-wave method, Phys. Rev. B59, 1758 (1999)
work page 1999
-
[42]
J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)
work page 1996
-
[43]
M. Ceriotti, M. Parrinello, T. E. Markland, and D. E. Manolopoulos, Efficient stochastic thermostatting of path integral molecular dynamics, The Journal of Chem- ical Physics133, 124104 (2010)
work page 2010
-
[44]
Giustino, Electron-phonon interactions from first prin- ciples, Rev
F. Giustino, Electron-phonon interactions from first prin- ciples, Rev. Mod. Phys.89, 015003 (2017)
work page 2017
-
[45]
Y. Ding, H. Chen, and J. Shi, Anharmonicity and Coulomb pseudopotential effects on superconductivity in YH6 and YH9, Phys. Rev. B112, 184517 (2025)
work page 2025
-
[46]
P. B. Allen and R. C. Dynes, Transition temperature of strong-coupled superconductors reanalyzed, Phys. Rev. B12, 905 (1975)
work page 1975
-
[47]
Chen,https://github.com/Haoran-Chen-1115/ SPIA
H. Chen,https://github.com/Haoran-Chen-1115/ SPIA
-
[48]
P. Morel and P. W. Anderson, Calculation of the su- perconducting state parameters with retarded electron- phonon interaction, Phys. Rev.125, 1263 (1962)
work page 1962
-
[49]
J. M. McMahon and D. M. Ceperley, High-temperature superconductivity in atomic metallic hydrogen, Phys. Rev. B84, 144515 (2011)
work page 2011
-
[50]
P. Cudazzo, G. Profeta, A. Sanna, A. Floris, A. Con- tinenza, S. Massidda, and E. K. U. Gross,Ab Ini- tiodescription of high-temperature superconductivity in dense molecular hydrogen, Phys. Rev. Lett.100, 257001 (2008)
work page 2008
-
[51]
Ding,https://doi.org/10.5281/zenodo.19592709 (2026)
Y. Ding,https://doi.org/10.5281/zenodo.19592709 (2026)
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.