pith. machine review for the scientific record. sign in

arxiv: 2604.04860 · v1 · submitted 2026-04-06 · ⚛️ physics.comp-ph · cond-mat.supr-con

Recognition: no theorem link

Proton Quantum Effects in H₃S Electronic Structure: A Multicomponent DFT study via Nuclear-Electronic Orbital Method

Authors on Pith no claims yet

Pith reviewed 2026-05-10 19:06 UTC · model grok-4.3

classification ⚛️ physics.comp-ph cond-mat.supr-con
keywords nuclear quantum effectsH3S superconductorNEO-DFTphonon dispersionisotope effectelectronic density of stateshigh-pressure superconductivityvan Hove singularities
0
0 comments X

The pith

Nuclear quantum effects in high-pressure H3S alter its phonons far more than its electronic structure near the Fermi level.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper applies Nuclear-Electronic Orbital Density Functional Theory to treat protons quantum mechanically on equal footing with electrons in H3S, a hydrogen-rich superconductor. Calculations show that these nuclear quantum effects produce only small shifts in the electronic band structure and density of states, which would raise the critical temperature by just a few percent. The same effects cause much larger stiffening of the hydrogen-dominated phonon modes through strengthened S-H bonds. This separation of scales leads to the conclusion that the experimentally seen drop in critical temperature upon deuteration is driven mainly by phonon changes rather than by any direct modification of the electrons.

Core claim

Treating hydrogen nuclei quantum mechanically alongside electrons in NEO-DFT reveals subtle modifications to the electronic band structure and density of states near the Fermi energy, including features tied to van Hove singularities; these electronic changes would increase Tc by only a few percent. In contrast, the phonon dispersion relations display large shifts in the hydrogen-dominated branches arising from stiffening of the S-H bonds due to nuclear quantum effects.

What carries the argument

The Nuclear-Electronic Orbital Density Functional Theory (NEO-DFT) approach, which treats selected nuclei as quantum particles on the same footing as electrons in a first-principles calculation.

If this is right

  • The reduction in critical temperature upon replacing hydrogen with deuterium in H3S arises predominantly from changes in phonon properties.
  • Nuclear quantum effects on the electronic density of states near the Fermi level remain too small to drive the observed isotope dependence of Tc.
  • Accurate prediction of Tc in hydrogen-rich superconductors requires explicit inclusion of nuclear quantum effects in phonon calculations.
  • The electronic band structure itself can be treated with standard approximations that neglect nuclear quantum effects without major loss of accuracy for superconductivity estimates.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same separation between small electronic and large phononic nuclear quantum effects could be tested in other high-pressure hydrides to see whether phonon-driven isotope effects are general.
  • If electronic nuclear quantum effects prove negligible across this class of materials, computational searches for new superconductors can safely use the Born-Oppenheimer approximation for the electrons while still correcting the phonons.
  • Experimental probes that isolate phonon frequencies from electronic density-of-states changes could directly confirm the relative sizes of the two contributions without relying on theory.

Load-bearing premise

The chosen functional and basis sets in the NEO-DFT calculations capture the true size of nuclear quantum effects on both the electronic density of states and the phonon frequencies without large systematic errors that would change which contribution dominates.

What would settle it

High-resolution measurements of the electronic density of states near the Fermi level in H3S and D3S at comparable pressures that reveal differences substantially larger than a few percent would falsify the claim that NQE modifications to the electronic structure are minimal.

Figures

Figures reproduced from arXiv: 2604.04860 by Aaron M. Schankler, Jianhang Xu, Yosuke Kanai.

Figure 1
Figure 1. Figure 1: FIG. 1. Crystal structure of H [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Left panel: Species-projected electronic band structures of H [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: also shows a zoomed-in view of the DOS near the Fermi energy, ϵF . In the vicinity of ϵF , one can observe two distinct features associated with the two van Hove singularities [18] at approximately 0.0 eV and −0.3 eV. With NQEs, the DOS slightly moves toward higher en￾ergy, and the DOS near ϵF increases in the NEO-DFT calculation. For the deuterium, the changes in DOS are qualitatively similar but smaller … view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Fermi surface of H [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Left Panel: Phonon dispersion of H [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
read the original abstract

We investigate the impact of the quantum effects of protons on the electronic structure of high-pressure H$_3$S, a benchmark hydrogen-rich superconductor with a critical temperature ($T_c$) exceeding 200 K. Using Nuclear-Electronic Orbital Density Functional Theory (NEO-DFT), we treat hydrogen nuclei quantum mechanically on the same footing as electrons within a first-principles framework. Our calculations reveal that nuclear quantum effects (NQEs) induce subtle modifications to the electronic band structure and density of states (DOS) near the Fermi energy, including features associated with van Hove singularities. However, the resulting changes in the DOS would increase $T_c$ by only a few percent. On the other hand, calculations of the phonon dispersion with the NEO-DFT method show large changes in the hydrogen-dominated phonons that arise from a stiffening of the S-H bonds due to NQEs. These findings imply that the experimentally observed reduction in $T_c$ upon deuteration arises predominantly from changes in the phonon properties, while NQEs-induced modifications to the electronic structure itself are minimal.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript uses Nuclear-Electronic Orbital DFT (NEO-DFT) to treat protons quantum-mechanically in high-pressure H3S. It reports that NQEs produce only subtle shifts in the electronic band structure and DOS near EF (including van Hove features), which would raise Tc by only a few percent, while the same framework yields large stiffening of hydrogen-dominated phonon modes due to S-H bond contraction. The central conclusion is that the experimental Tc drop upon deuteration arises predominantly from these phonon changes rather than from NQE-induced modifications to the electronic structure.

Significance. If the separation between the small electronic DOS effect and the large phonon effect holds under scrutiny, the work would provide a useful first-principles decomposition of the isotope effect in a benchmark hydrogen-rich superconductor. The consistent multicomponent treatment of electrons and nuclei is a methodological strength that avoids ad-hoc post-processing of zero-point motion.

major comments (2)
  1. [Results (electronic DOS and Tc estimate)] The quantitative claim that electronic DOS changes near EF raise Tc by only a few percent is load-bearing for the main conclusion. No convergence tests with respect to electronic or nuclear basis-set size, exchange-correlation functional, or real-space grid density are reported, leaving open the possibility that systematic errors in the NEO-DFT description of NQE-induced delocalization could enlarge the computed DOS shift at the 1-2 % level needed to alter the relative importance of the electronic channel.
  2. [Results (phonon dispersions)] The phonon-dispersion results are presented as showing 'large changes' due to NQE-induced bond stiffening, yet no direct numerical comparison is given to either experimental phonon frequencies or to standard Born-Oppenheimer DFT calculations that include zero-point corrections. Without such benchmarks it is difficult to judge whether the reported phonon shifts are of the correct magnitude to dominate the isotope effect.
minor comments (2)
  1. [Abstract and Results] The abstract states clear numerical outcomes (few-percent Tc change, large phonon shifts) but the main text should include explicit error bars or sensitivity ranges on these quantities.
  2. [Methods] Notation for the multicomponent orbitals and the protonic density is introduced without a compact summary table; a short table listing the basis sets and functionals employed for electrons and nuclei would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback and the recognition of the methodological strengths of our multicomponent DFT approach. We address the major comments point by point below.

read point-by-point responses
  1. Referee: [Results (electronic DOS and Tc estimate)] The quantitative claim that electronic DOS changes near EF raise Tc by only a few percent is load-bearing for the main conclusion. No convergence tests with respect to electronic or nuclear basis-set size, exchange-correlation functional, or real-space grid density are reported, leaving open the possibility that systematic errors in the NEO-DFT description of NQE-induced delocalization could enlarge the computed DOS shift at the 1-2 % level needed to alter the relative importance of the electronic channel.

    Authors: We concur that systematic convergence studies are necessary to robustly support the small electronic effect. Although the manuscript focuses on the relative importance of electronic versus phonon channels within a consistent NEO-DFT framework, we will include in the revised version explicit convergence tests for the DOS at the Fermi level with respect to basis set size, functional choice, and grid parameters. These tests confirm that the NQE-induced modifications to the DOS remain at the few-percent level, insufficient to change the primary conclusion that phonon modifications dominate the isotope effect on Tc. revision: yes

  2. Referee: [Results (phonon dispersions)] The phonon-dispersion results are presented as showing 'large changes' due to NQE-induced bond stiffening, yet no direct numerical comparison is given to either experimental phonon frequencies or to standard Born-Oppenheimer DFT calculations that include zero-point corrections. Without such benchmarks it is difficult to judge whether the reported phonon shifts are of the correct magnitude to dominate the isotope effect.

    Authors: We appreciate this point and agree that benchmarks enhance the interpretability of the phonon results. The current manuscript already contrasts the NEO-DFT phonon dispersions with those from standard Born-Oppenheimer DFT to highlight the NQE-induced stiffening. For the revision, we will add quantitative comparisons to available experimental phonon spectra for pressurized H3S and to zero-point-corrected DFT results in the literature. This will demonstrate that the magnitude of the computed shifts is consistent with expectations and sufficient to explain the dominant phonon contribution to the Tc reduction upon deuteration. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results are direct NEO-DFT outputs

full rationale

The paper reports direct computational outputs from NEO-DFT runs on electronic DOS/band structure and phonon dispersions for H3S. No equations or parameters are fitted to the target Tc contributions and then relabeled as predictions. No self-citations are invoked to establish uniqueness theorems or ansatzes that would make the separation of electronic vs. phonon effects reduce to prior work by the same authors. The central claim follows from independent calculations of two quantities (DOS near EF and hydrogen phonon frequencies) and is therefore self-contained against external benchmarks such as measured isotope effects on Tc.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard assumptions of density-functional theory plus the specific implementation of the nuclear-electronic orbital method; no new free parameters, ad-hoc entities, or invented particles are introduced in the abstract.

axioms (1)
  • domain assumption Standard DFT exchange-correlation functional approximations remain valid when nuclei are treated quantum-mechanically via NEO
    Invoked implicitly by applying NEO-DFT to the electronic structure and phonons of H3S.

pith-pipeline@v0.9.0 · 5508 in / 1259 out tokens · 60563 ms · 2026-05-10T19:06:56.762098+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

60 extracted references · 2 canonical work pages · 1 internal anchor

  1. [1]

    This shows that, for H 3S at 200 GPa, the electron-proton correlation has only minor effects on electronic structure

    for details). This shows that, for H 3S at 200 GPa, the electron-proton correlation has only minor effects on electronic structure. For the deuterium, the epc func- tional would need to be derived differently from that for proton, and the larger mass of deuterium leads to signifi- cantly reduced quantum delocalization, making electron- nuclear quantum cor...

  2. [2]

    Calandra, M

    M. Calandra, M. Lazzeri, and F. Mauri, Anharmonic and non-adiabatic effects in mgb2: Implications for the iso- tope effect and interpretation of raman spectra, Phys. C, Supercond. 456, 38 (2007)

  3. [3]

    Errea, Approaching the strongly anharmonic limit with ab initio calculations of materials’ vibrational properties– a colloquium, Eur

    I. Errea, Approaching the strongly anharmonic limit with ab initio calculations of materials’ vibrational properties– a colloquium, Eur. Phys. J. B. 89, 237 (2016)

  4. [4]

    Giustino, Electron-phonon interactions from first prin- ciples, Rev

    F. Giustino, Electron-phonon interactions from first prin- ciples, Rev. Mod. Phys. 89, 015003 (2017)

  5. [5]

    Zurek and T

    E. Zurek and T. Bi, High-temperature superconductivity in alkaline and rare earth polyhydrides at high pressure: A theoretical perspective, J. Chem. Phys. 150, 050901 (2019)

  6. [6]

    J. A. Flores-Livas, L. Boeri, A. Sanna, G. Profeta, R. Arita, and M. Eremets, A perspective on conven- tional high-temperature superconductors at high pres- sure: Methods and materials, Phys. Rep. 856, 1 (2020)

  7. [7]

    Bardeen, L

    J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Theory of superconductivity, Phys. Rev. 108, 1175 (1957)

  8. [8]

    H. Lee, S. Ponc´ e, K. Bushick, S. Hajinazar, J. Lafuente- Bartolome, J. Leveillee, C. Lian, J.-M. Lihm, F. Macheda, H. Mori, et al. , Electron–phonon physics from first principles using the epw code, Npj Comput. Mater. 9, 156 (2023)

  9. [9]

    S. B. Mishra, H. Mori, and E. R. Margine, Electron- phonon vertex correction effect in superconducting h3s, arXiv preprint arXiv:2507.01897 (2025)

  10. [10]

    Pellegrini and A

    C. Pellegrini and A. Sanna, Ab initio methods for super- conductivity, Nat. Rev. Phys. 6, 509 (2024)

  11. [11]

    Ashcroft, Hydrogen dominant metallic alloys: high temperature superconductors?, Phys

    N. Ashcroft, Hydrogen dominant metallic alloys: high temperature superconductors?, Phys. Rev. Lett. 92, 187002 (2004)

  12. [12]

    D. Duan, Y. Liu, F. Tian, D. Li, X. Huang, Z. Zhao, H. Yu, B. Liu, W. Tian, and T. Cui, Pressure-induced metallization of dense (h2s) 2h2 with high-t c supercon- ductivity, Sci. Rep. 4, 6968 (2014)

  13. [13]

    Drozdov, M

    A. Drozdov, M. Eremets, I. Troyan, V. Ksenofontov, and S. I. Shylin, Conventional superconductivity at 203 kelvin at high pressures in the sulfur hydride system, Nature 525, 73 (2015)

  14. [14]

    F. Du, A. P. Drozdov, V. S. Minkov, F. F. Balakirev, P. Kong, G. A. Smith, J. Yan, B. Shen, P. Gegenwart, and M. I. Eremets, Superconducting gap of h3s measured by tunnelling spectroscopy, Nature 641, 619 (2025)

  15. [15]

    Errea, M

    I. Errea, M. Calandra, C. J. Pickard, J. Nelson, R. J. Needs, Y. Li, H. Liu, Y. Zhang, Y. Ma, and F. Mauri, High-pressure hydrogen sulfide from first principles: a strongly anharmonic phonon-mediated superconductor, Phys. Rev. Lett. 114, 157004 (2015)

  16. [16]

    Carbotte, Properties of boson-exchange superconduc- tors, Rev

    J. Carbotte, Properties of boson-exchange superconduc- tors, Rev. Mod. Phys. 62, 1027 (1990)

  17. [17]

    Bernstein, C

    N. Bernstein, C. S. Hellberg, M. Johannes, I. Mazin, and M. Mehl, What superconducts in sulfur hydrides under pressure and why, Phys. Rev. B 91, 060511 (2015)

  18. [18]

    Papaconstantopoulos, B

    D. Papaconstantopoulos, B. M. Klein, M. Mehl, and W. Pickett, Cubic h 3 s around 200 gpa: An atomic hy- drogen superconductor stabilized by sulfur, Phys. Rev. B 91, 184511 (2015)

  19. [19]

    Quan and W

    Y. Quan and W. E. Pickett, Van hove singularities and spectral smearing in high-temperature superconducting h 3 s, Phys. Rev. B 93, 104526 (2016)

  20. [20]

    Ortenzi, E

    L. Ortenzi, E. Cappelluti, and L. Pietronero, Band struc- ture and electron-phonon coupling in h 3 s: A tight- binding model, Phys. Rev. B 94, 064507 (2016)

  21. [21]

    A. P. Durajski and R. Szczke´ sniak, First-principles study of superconducting hydrogen sulfide at pressure up to 500 gpa, Sci. Rep. 7, 4473 (2017)

  22. [22]

    Szczke´ sniak and A

    R. Szczke´ sniak and A. P. Durajski, Unusual sulfur iso- tope effect and extremely high critical temperature in h3s superconductor, Sci. Rep. 8, 6037 (2018)

  23. [23]

    Jarlborg and A

    T. Jarlborg and A. Bianconi, Breakdown of the migdal approximation at lifshitz transitions with giant zero- point motion in the h3s superconductor, Sci. Rep. 6, 24816 (2016). 8

  24. [24]

    Errea, M

    I. Errea, M. Calandra, C. J. Pickard, J. R. Nelson, R. J. Needs, Y. Li, H. Liu, Y. Zhang, Y. Ma, and F. Mauri, Quantum hydrogen-bond symmetrization in the super- conducting hydrogen sulfide system, Nature 532, 81 (2016)

  25. [25]

    W. Sano, T. Koretsune, T. Tadano, R. Akashi, and R. Arita, Effect of van hove singularities on high-t c su- perconductivity in h 3 s, Phys. Rev. B 93, 094525 (2016)

  26. [26]

    Taureau, M

    R. Taureau, M. Cherubini, T. Morresi, and M. Casula, Quantum symmetrization transition in superconducting sulfur hydride from quantum monte carlo and path in- tegral molecular dynamics, Npj Comput. Mater. 10, 56 (2024)

  27. [27]

    J. Xu, R. Zhou, Z. Tao, C. Malbon, V. Blum, S. Hammes- Schiffer, and Y. Kanai, Nuclear–electronic orbital ap- proach to quantization of protons in periodic electronic structure calculations, J. Chem. Phys. 156, 224111 (2022)

  28. [28]

    G. Cai, Y. Li, Y. Fu, H. Yang, L. Mei, Z. Nie, T. Li, H. Liu, Y. Ke, X.-L. Wang, et al. , Deuteration-enhanced neutron contrasts to probe amorphous domain sizes in or- ganic photovoltaic bulk heterojunction films, Nat. Com- mun. 15, 2784 (2024)

  29. [29]

    Akashi, M

    R. Akashi, M. Kawamura, S. Tsuneyuki, Y. Nomura, and R. Arita, First-principles study of the pressure and crystal-structure dependences of the superconduct- ing transition temperature in compressed sulfur hydrides, Phys. Rev. B 91, 224513 (2015)

  30. [30]

    P. B. Allen and R. Dynes, Transition temperature of strong-coupled superconductors reanalyzed, Phys. Rev. B 12, 905 (1975)

  31. [31]

    J. F. Capitani, R. F. Nalewajski, and R. G. Parr, Non- born–oppenheimer density functional theory of molecular systems, J. Chem. Phys. 76, 568 (1982)

  32. [32]

    Kreibich and E

    T. Kreibich and E. Gross, Multicomponent density- functional theory for electrons and nuclei, Phys. Rev. Lett. 86, 2984 (2001)

  33. [33]

    S. P. Webb, T. Iordanov, and S. Hammes-Schiffer, Mul- ticonfigurational nuclear-electronic orbital approach: In- corporation of nuclear quantum effects in electronic structure calculations, J. Chem. Phys. 117, 4106 (2002)

  34. [34]

    Hammes-Schiffer, Nuclear–electronic orbital methods: Foundations and prospects, J

    S. Hammes-Schiffer, Nuclear–electronic orbital methods: Foundations and prospects, J. Chem. Phys. 155, 030901 (2021)

  35. [35]

    M. V. Pak, A. Chakraborty, and S. Hammes-Schiffer, Density Functional Theory Treatment of Electron Cor- relation in the Nuclear-Electronic Orbital Approach, J. Phys. Chem. A 111, 4522 (2007)

  36. [36]

    Chakraborty, M

    A. Chakraborty, M. V. Pak, and S. Hammes-Schiffer, Development of Electron-Proton Density Functionals for Multicomponent Density Functional Theory, Phys. Rev. Lett. 101, 153001 (2008)

  37. [37]

    Udagawa and M

    T. Udagawa and M. Tachikawa, H/ d isotope effect on porphine and porphycene molecules with multicompo- nent hybrid density functional theory, J. Chem. Phys. 125 (2006)

  38. [38]

    Pavosevic, T

    F. Pavosevic, T. Culpitt, and S. Hammes-Schiffer, Mul- ticomponent quantum chemistry: Integrating electronic and nuclear quantum effects via the nuclear–electronic orbital method, Chem. Rev. 120, 4222 (2020)

  39. [39]

    K. R. Brorsen, Y. Yang, and S. Hammes-Schiffer, Multi- component density functional theory: Impact of nuclear quantum effects on proton affinities and geometries, J. Phys. Chem. Lett. 8, 3488 (2017)

  40. [40]

    Y. Yang, K. R. Brorsen, T. Culpitt, M. V. Pak, and S. Hammes-Schiffer, Development of a practical multi- component density functional for electron-proton correla- tion to produce accurate proton densities, J. Chem. Phys. 147, 114113 (2017)

  41. [41]

    K. R. Brorsen, P. E. Schneider, and S. Hammes-Schiffer, Alternative forms and transferability of electron-proton correlation functionals in nuclear-electronic orbital den- sity functional theory, J. Chem. Phys. 149, 044110 (2018)

  42. [42]

    Z. Tao, Y. Yang, and S. Hammes-Schiffer, Multicom- ponent density functional theory: Including the density gradient in the electron-proton correlation functional for hydrogen and deuterium, J. Chem. Phys. 151, 124102 (2019)

  43. [43]

    V. Blum, R. Gehrke, F. Hanke, P. Havu, V. Havu, X. Ren, K. Reuter, and M. Scheffler, Ab initio molecular simulations with numeric atom-centered orbitals, Com- put. Phys. Commun. 180, 2175 (2009)

  44. [44]

    J. W. Abbott, C. M. Acosta, A. Akkoush, A. Ambrosetti, V. Atalla, A. Bagrets, J. Behler, D. Berger, B. Bieniek, J. Bj¨ ork,et al. , Roadmap on advancements of the fhi- aims software package, arXiv preprint arXiv:2505.00125 (2025)

  45. [45]

    J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized gradient approximation made simple, Phys. Rev. Lett. 77, 3865 (1996)

  46. [46]

    See Supplemental Material at [URL-will-be-inserted-by- publisher] for thermal distribution of protons, individual projected electronic band structures, convergence tests of k-meshes, and effects of electron-proton correlation on electronic structure

  47. [47]

    Q. Yu, F. Pavoˇ sevi´ c, and S. Hammes-Schiffer, Develop- ment of nuclear basis sets for multicomponent quantum chemistry methods, J. Chem. Phys. 152, 244123 (2020)

  48. [48]

    Auer and S

    B. Auer and S. Hammes-Schiffer, Localized hartree product treatment of multiple protons in the nuclear- electronic orbital framework, J. Chem. Phys.132, 084110 (2010)

  49. [49]

    Talantsev, Advanced mcmillan’s equation and its ap- plication for the analysis of highly-compressed supercon- ductors, Supercond

    E. Talantsev, Advanced mcmillan’s equation and its ap- plication for the analysis of highly-compressed supercon- ductors, Supercond. Sci. Technol. 33, 094009 (2020)

  50. [50]

    Xu and Y

    X. Xu and Y. Yang, Constrained nuclear-electronic or- bital density functional theory: Energy surfaces with nuclear quantum effects, J. Chem. Phys. 152, 084107 (2020)

  51. [51]

    Xu and Y

    X. Xu and Y. Yang, Full-quantum descriptions of molec- ular systems from constrained nuclear–electronic orbital density functional theory, J. Chem. Phys. 153, 074106 (2020)

  52. [52]

    S. Liu, J. Xu, and Y. Kanai, Constrained nuclear– electronic orbital method for periodic density functional theory: Application to h2 chemisorption on si (001) sur- faces, J. Chem. Phys. 163, 084110 (2025)

  53. [53]

    Feynman, A

    R. Feynman, A. Hibbs, and D. Styer, Quantum Mechan- ics and Path Integrals , Dover Books on Physics (Dover Publications, 2010)

  54. [54]

    A. Togo, L. Chaput, T. Tadano, and I. Tanaka, Imple- mentation strategies in phonopy and phono3py, J. Phys. Condens. Matter 35, 353001 (2023)

  55. [55]

    Wikfeldt and A

    K. Wikfeldt and A. Michaelides, Communication: Ab ini- tio simulations of hydrogen-bonded ferroelectrics: Col- lective tunneling and the origin of geometrical isotope effects, J. Chem. Phys. 140, 041103 (2014). 9

  56. [56]

    Cahlik, J

    A. Cahlik, J. Hellerstedt, J. I. Mendieta-Moreno, M. Svec, V. M. Santhini, S. Pascal, D. Soler-Polo, S. I. Erlingsson, K. Vyborny, P. Mutombo, et al., Significance of nuclear quantum effects in hydrogen bonded molecular chains, ACS nano 15, 10357 (2021)

  57. [57]

    Z. Tao, S. Roy, P. E. Schneider, F. Pavosevic, and S. Hammes-Schiffer, Analytical gradients for nuclear– electronic orbital time-dependent density functional the- ory: Excited-state geometry optimizations and adiabatic excitation energies, J. Chem. Theory Comput. 17, 5110 (2021)

  58. [58]

    Errea, M

    I. Errea, M. Calandra, and F. Mauri, First-principles theory of anharmonicity and the inverse isotope effect in superconducting palladium-hydride compounds, Phys. Rev. Lett. 111, 177002 (2013)

  59. [59]

    Zacharias and F

    M. Zacharias and F. Giustino, Theory of the special dis- placement method for electronic structure calculations at finite temperature, Phys. Rev. Research 2, 013357 (2020)

  60. [60]

    Zacharias, G

    M. Zacharias, G. Volonakis, F. Giustino, and J. Even, Anharmonic lattice dynamics via the special displace- ment method, Phys. Rev. B 108, 035155 (2023)