Pith. sign in

REVIEW 3 major objections 4 minor 74 references

This paper locates a displacive quantum critical point in superconducting H3S at about 134 GPa and argues that the superconducting dome peaks in the paraelectric region of strong nuclear quantum fluctuations, rather than at the structural p

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-03 05:52 UTC pith:SUUQK6MX

load-bearing objection Solid PIMD phase diagram for H3S with a plausible QCP; the 4D Ising claim is an assumed consistency check, not an independent determination, but the paper deserves refereeing. the 3 major comments →

arxiv 2602.00833 v2 pith:SUUQK6MX submitted 2026-01-31 cond-mat.supr-con

Displacive quantum critical point in superconducting hydrides: The case of H₃S

Marco Cherubini , Abhishek Raghav , Michele Casula This is my paper
classification cond-mat.supr-con
keywords H3Ssulfur hydridequantum critical pointferroelectric transitionpath integral molecular dynamicsnuclear quantum effects4D Ising universality classsuperconductivity dome
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper tries to establish that the structural phase diagram of H3S is governed by a ferroelectric transition whose phase boundary ends at a quantum critical point near 134 GPa, and that the measured superconducting critical temperature maximum around 155 GPa falls on the paraelectric side of this QCP, where hydrogen nuclei fluctuate strongly. If correct, this reframes the explanation for high-temperature superconductivity in H3S: the superconducting dome is shaped by proximity to a quantum critical point and by retarded, beyond-Migdal phonon dynamics, not by the structural transition itself. The argument uses path integral molecular dynamics with a machine-learned interatomic potential trained on density functional theory, together with finite-size scaling that assigns the transition to the 4D Ising universality class. A sympathetic reader would care because it turns a celebrated superconductor into a quantum-critical system, with concrete implications for how its superconductivity should be calculated and for other hydrogen-rich superconductors.

Core claim

The authors claim that in H3S, the protons shuttling between flanking sulfur atoms drive a displacive ferroelectric transition from a centrosymmetric cubic phase (Im-3m) to a polar trigonal phase (R3m), and that this transition line terminates at a quantum critical point at p_QCP ≈ 134 ± 2 GPa, where the phase boundary reaches zero temperature and the transition is driven purely by quantum fluctuations. A finite-size scaling analysis of the proton-displacement order parameter places the transition in the 4D Ising universality class, with correlation-length exponent ν ≈ 0.483 ± 0.014. The central physical claim is that the experimentally observed Tc peak near 155 GPa lies not at the structura

What carries the argument

The central object is the local proton displacement along the S–S direction, which serves as the order parameter for the ferroelectric transition; its average, variance, and distribution distinguish the paraelectric from the polar phase. The paper also uses the imaginary-time-resolved phonon Green function g(τ), specifically the ratio σ²_g(β/2)/σ²_g(0), to quantify retarded quantum fluctuations of local dipole moments, and a finite-size scaling collapse of the order parameter (with the 4D Ising form, z = 1, x_Δ = 1) to extract the critical exponent ν and extrapolate the QCP pressure using logarithmic corrections at the upper critical dimension.

Load-bearing premise

The load-bearing premise is that the transition belongs to the 4D Ising universality class with a fixed time–space scaling (z = 1) and a unit scaling dimension of the order parameter; the paper assumes this because its simulation grid cannot determine z directly, and if the true critical behavior deviates, the extrapolated quantum critical point at 134 GPa and the claim that the superconducting peak sits in a quantum-fluctuation region do not follow.

What would settle it

Run path integral molecular dynamics at lower temperature, e.g. 25 K, with supercells large enough to measure the imaginary-time decay of the order-parameter correlations, and extract the dynamical exponent z without assuming z = 1; if the data require z ≠ 1, or if the order-parameter collapse fails for ν = 1/2, the claimed universality class and the extrapolated QCP pressure are not correct. Experimentally, detecting a sharp structural transition at pressures above about 155 GPa in H3S would contradict the claim that the superconducting peak sits in the paraelectric phase.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • The maximum of the superconducting Tc dome falls on the paraelectric side of the quantum critical point, so superconductivity is not caused by the ferroelectric ordering itself but by the quantum fluctuations surrounding the QCP.
  • Nuclear quantum effects shift the structural phase boundary down by roughly 50 GPa relative to classical nuclei at 200 K, making explicit quantum treatment of the protons essential to place the phase diagram next to the experimentally observed dome.
  • Because the transition belongs to the 4D Ising universality class, its critical fluctuations are only marginally irrelevant and extend over a broad pressure–temperature region, consistent with a wide region of enhanced retardation effects.
  • The softening of the optical hydrogen shuttling modes and the peak in the two-body phonon Green function near the transition indicate that two-phonon and other beyond-Migdal pairing channels are amplified close to the QCP.
  • An accurate calculation of Tc in H3S must go beyond the standard Migdal–Eliashberg framework and include the full frequency dependence of the phonon Green function and vertex corrections.

Where Pith is reading between the lines

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

  • If the QCP picture is correct, isotope substitution H → D should shift the superconducting dome and the quantum-fluctuation ridge together in a way that differs from a simple phonon-frequency renormalization, giving a testable signature.
  • The paper constrains the simulation cell to cubic symmetry, thereby excluding coupling of the proton order parameter to trigonal lattice strain; if that coupling matters, the transition could become weakly first order and the 4D Ising classification would need to be revisited.
  • Showing that two-phonon correlations peak near the transition is not the same as showing that they raise Tc; an actual non-perturbative computation of the superconducting critical temperature in this regime is the natural next step to confirm the proposed mechanism.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The manuscript uses path-integral molecular dynamics (PIMD) with a MACE machine-learned potential trained on DFT-BLYP data to map the temperature-pressure phase diagram of H3S. It reports a quantum displacive/ferroelectric transition between the Im-3m paraelectric and R3m ferroelectric phases, with a quantum critical point at p_QCP ≈ 134 ± 2 GPa, and identifies the transition as belonging to the 4D Ising universality class. The authors argue that the experimental superconducting Tc dome, peaked near 155 GPa, lies in the paraelectric region where nuclear quantum fluctuations and retarded phonon correlations are strongest, suggesting a possible connection between the QCP and superconductivity. The paper also provides anharmonic phonon dispersions, local proton displacement distributions, and imaginary-time Green-function analyses to support the phase diagram and fluctuation picture.

Significance. If correct, the claimed QCP at ~134 GPa and its 4D Ising character would substantially reframe the superconducting mechanism in H3S, moving the focus from the static structural boundary to quantum-critical fluctuations and non-Migdal electron-phonon effects. The methodological package is strong: the MLIP validation is reproducible and accurate (MAE 0.363 meV/atom and 5.25 meV/Å), the finite-temperature phase boundaries are obtained by direct PIMD simulation rather than by reconstruction from harmonic or perturbative approaches, and the open data policy is commendable. The paper also offers a concrete, falsifiable prediction (p_QCP ≈ 134 GPa) and a specific macroscopic picture (Tc peak in the paraelectric regime) that can be tested by future experiment and by independent calculations.

major comments (3)
  1. [Eq. (4) and SM Eq. (1)] The finite-size scaling analysis assumes the 4D Ising universality class from the outset: z=1, x_Δ=1, and the logarithmic L-dependence of Eq. (1) in the SM. The extracted ν ≈ 0.483 ± 0.014 is therefore a consistency check, not an independent determination of the universality class. Moreover, the data collapse in Fig. 4 uses exactly three points (L=4 at T=50 K, L=3 at T=67 K, L=2 at T=100 K), all on the z=1 ray. Such a collapse cannot constrain z or the anomalous dimension, and a different z or nonzero η would shift the critical volume, and hence p_QCP. This is a load-bearing issue because both the 4D Ising claim and the location of the QCP rest on this scaling form. I recommend either softening the universality-class claim to 'consistent with 4D Ising' or adding simulations at additional sizes and temperatures, including points off the z=1 ray.
  2. [SM Methods and Sec. II] The cell is kept fixed to the cubic Im-3m shape in all PIMD simulations. The order parameter is a polar proton-displacement mode, which in a real crystal couples to acoustic strain. Such strain coupling is known to alter the universality class of a ferroelectric quantum phase transition (e.g., driving it first order or changing the exponents). The manuscript dismisses the trigonal distortion as statically small (≤0.05°), but this does not address dynamic strain fluctuations at the QCP. Since the universality class assignment is a headline result, the neglect of strain fluctuations is a major limitation that needs to be confronted, either by relaxing the cell shape in at least a subset of simulations or by a quantitative estimate of the strain-order-parameter coupling.
  3. [QCP extrapolation (Fig. 1 and SM Fig. 3)] The QCP is obtained by extrapolating the finite-size-corrected transition line down to T=0. The lowest temperature simulated is 50 K, and the functional form used for the T→0 extrapolation is not specified in the main text. The reported error bar of ±2 GPa appears to be a fit error, but it does not include the systematic uncertainty in the assumed extrapolation function. Given that the 50 K transition already lies within a few GPa of the claimed QCP, the T→0 extrapolation is a nontrivial step that should be documented explicitly (e.g., the fitting form, the included temperatures, and a sensitivity check). Without this, the quoted p_QCP value is not fully supported.
minor comments (4)
  1. [Throughout] There are several typos in the Supplemental Material: 'denisty' should be 'density' (SM Methods), 'anhamonic' → 'anharmonic' (SM PIMD phonons), 'ferromagnetic phase transition' → 'ferroelectric phase transition' (the sentence before Fig. 10), and 'distrbution' → 'distribution' (SM Fig. 6 caption).
  2. [Fig. 1 and Fig. 2(d)] The heatmap of σ²_g(β/2)/σ²_g(0) is visually useful, but the pressure-temperature interpolation is done with radially symmetric basis functions without a stated smoothing parameter. The authors may want to state that the interpolated boundaries are guides, not measurements, especially in the gray region noted in the SM.
  3. [Main text around Eq. (2)] The definition of the local displacement Δ_i^(j) is given for centroids, while the order parameter Δ in Eq. (1) is defined for centroid positions as well. It is useful that the SM shows bead-resolved results are consistent, but the main text would benefit from a sentence explaining that bead resolution does not change the phase-boundary estimate.
  4. [References] The paper cites the relevant literature on hydrides, quantum criticality, and non-Migdal superconductivity. A direct reference to the finite-size scaling cost function and to the logarithmic corrections in 4D Ising models would be helpful; the current citations are adequate but the reader may struggle to find the precise forms.

Circularity Check

2 steps flagged

The 4D-Ising classification and the extrapolated QCP pressure rest on an explicitly assumed scaling form; the fit is then cited as evidence for the same assumed class.

specific steps
  1. self definitional [Main text, Eq. (4) and the 'Finite-size scaling of the ferroelectric transition line' paragraph]
    "we assume Lorentz invariance (z=1) and vanishing anomalous dimension (xΔ=1), fulfilled by the Ising universality class... Therefore, we assume the validity of the Ising universality class and, in order to check whether our data are compatible with this scenario, we focus on the critical exponent ν. Given z=1, we fix the Lτ/L ratio to be a constant, such that the scaling function in Eq. 4 depends only on one variable."

    The claimed 4D-Ising classification is an input, not an output: z=1 and xΔ=1 are fixed to the Ising values and only ν is fitted. The three collapse points are deliberately placed on the z=1 ray (Lτ/L constant), so the data cannot independently determine z or the anomalous dimension. The later statement that the transition 'belongs to the 4D Ising universality class' is therefore a consistency check with the assumed class, not a derivation of it.

  2. fitted input called prediction [Supplemental Material, Eq. (1) and the 'Finite-size scaling of the ferroelectric transition line' section; main text QCP extrapolation paragraph]
    "According to the 4D Ising universality class, the finite-size scaling of the critical coupling is given by: pQCP(L,T)=pQCP(∞,T)+α(T)L−2 ln(L)−1/6... The fits based on the functional form in Eq. 1 are also shown... Thus, Eq. 1 captures very well the spatial dependence of the critical couplings. It is another evidence of the validity of the 4D Ising universality class."

    The QCP pressure is obtained by fitting the 4D-Ising logarithmic finite-size form, so the extrapolated p_QCP≈134±2 GPa is conditional on the assumed class. The same successful fit is then cited as 'another evidence' for the 4D-Ising class. The fit cannot independently validate the form that was used to perform the extrapolation; this makes the central numerical QCP claim partly forced by the input universality hypothesis.

full rationale

The simulation pipeline itself (PIMD with a MACE potential trained on DFT-BLYP data, validated against an independent test set) is self-contained and provides direct data for the order parameter, phonon softening, and finite-size transition pressures; those raw results are not circular. The MLIP training set reuse of the authors' earlier trajectories is not load-bearing because the model is checked against held-out DFT configurations. However, the two headline claims—the 4D-Ising classification and the extrapolated QCP location—are established only through an explicitly assumed scaling form. The paper states that it assumes z=1 and xΔ=1, checks only ν, and then uses the 4D-Ising logarithmic correction to extrapolate p_QCP to the thermodynamic limit; the quality of that same fit is presented as evidence for the 4D-Ising class. This is a legitimate consistency check but not an independent determination, so the universality/QCP claim is partially circular. The experimental-Tc comparison and the discussion of quantum fluctuations are separate and not circular. Score 5 reflects this partial reduction of the central claim to its input assumption, while acknowledging the paper's honest disclosure and the independent raw simulation data.

Axiom & Free-Parameter Ledger

2 free parameters · 5 axioms · 0 invented entities

The central claims rest on the fidelity of the BLYP-MLIP potential and on the assumed 4D Ising finite-size scaling. No new physical entities are introduced, but the QCP location and universality class are derived under these structural assumptions.

free parameters (2)
  • α0 and α1 in α(T)=α0+α1T for the finite-size scaling fit = not reported numerically
    These coefficients are fitted to the size-dependent transition pressures in SM Eq. (1) and are used to extrapolate p_QCP(L=∞, T). The QCP position depends on this fit.
  • Critical exponent ν = 0.483±0.014
    ν is obtained by minimizing the data-collapse cost function C(ν); this fitted value is then used to argue for the 4D Ising universality class.
axioms (5)
  • domain assumption The BLYP exchange-correlation functional accurately describes the potential energy surface of H3S over the studied P-T range.
    The MLIP is trained on BLYP DFT data; all conclusions inherit the accuracy of BLYP for this system. No cross-check with other functionals is provided.
  • domain assumption The MACE MLIP trained on 14640 configurations faithfully reproduces DFT energies and forces in all simulated conditions, including the 50 K and 300 K extrapolations.
    The MLIP is validated on a test set, but the training data come from T=100 and 200 K from a previous work; the behavior at the extrapolated low- and high-temperature points is assumed to remain accurate.
  • ad hoc to paper The transition belongs to the 4D Ising universality class with z=1 and vanishing anomalous dimension.
    The finite-size scaling analysis assumes this class to extract ν and to extrapolate p_QCP; it is a hypothesis tested for compatibility, not independently determined.
  • domain assumption The trigonal lattice distortion accompanying the Im-3m to R3m transition is negligible and can be safely constrained to zero.
    The PIMD simulations constrain the cell vectors to the cubic shape, based on DFT estimates of ≤0.05° distortion. Strain fluctuations, which can change the universality class, are thereby suppressed.
  • domain assumption The order parameter defined by projecting proton positions onto the eight degenerate molecular orientations correctly describes the symmetry breaking.
    The order parameter is constructed by maximizing the global displacement over eight orientations; this procedure is standard but assumes the proton displacements are the primary order parameter.

pith-pipeline@v1.3.0-alltime-deepseek · 18093 in / 14336 out tokens · 156019 ms · 2026-08-03T05:52:34.935869+00:00 · methodology

0 comments
read the original abstract

H$_3$S sulfur hydride has been widely investigated for its high superconducting critical temperature $T_c$ of 203 K at about $p_c = 155$ GPa. Despite being the precursor of superconducting hydrides, a detailed picture of its structural phase diagram in an extended temperature and pressure range is still missing. To determine it with inclusion of both thermal and quantum effects, we carry out path integral molecular dynamics combined to a MACE neural network potential trained on BLYP density functional theory configurations. The resulting H$_3$S phase diagram is characterized by the displacive transition between the centrosymmetric Im$\bar{3}$m and polar R3m phases, which originates from a quantum critical point (QCP) located at $p_\mathrm{QCP} \approx 134$ GPa. We show that the experimental $T_c$ peak falls into a centrosymmetric region of large nuclear quantum fluctuations above the displacive QCP, as measured by local phonon Green's functions resolved in imaginary time, where fluctuating moments are at play. We study the critical behavior of the system in the proximity of the QCP by a finite-size scaling analysis, showing that it belongs to the 4D Ising universality class. We finally discuss its implications for the superconducting state.

Figures

Figures reproduced from arXiv: 2602.00833 by Abhishek Raghav, Marco Cherubini, Michele Casula.

Figure 1
Figure 1. Figure 1: H [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Order parameter and local observables for [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Soft optical mode frequencies as a function of pres [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Finite-size scaling of ∆ at T = 50 K (black trian￾gles), 67 K (blue circles) and 100 K (red stars), and supercell sizes chosen to keep Lτ /L constant. We use g[uL1/ν] as a functional form, where u = V −Vc Vc , with Vc the critical vol￾ume. Inset: the critical exponent ν is obtained by minimizing the cost function C(ν) defined in the text. corrections appearing at the upper critical dimension[37– 39], such … view at source ↗
Figure 1
Figure 1. Figure 1: Evaluation of the MLIP on a 3660 configuration test set. Panel (a): Comparison between the energies predicted [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Ferroelectric critical pressures as a function of 1 [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Ferroelectric transition lines computed at different sizes, as reported in the legend. The extrapolated values ( [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Panel (a): Order parameter for the classical ferroelectric transition. Full and empty symbols indicate the absolute [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Proton centroid PIMD distributions along the S-S direction for several pressures at [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Classical proton MD distributions along the S-S direction for several pressures at [PITH_FULL_IMAGE:figures/full_fig_p012_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Quantum and classical phase diagram of sulfur hydride for [PITH_FULL_IMAGE:figures/full_fig_p013_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Order parameter at 50 K as a function of pressure. Full and empty symbols correspond to the real values of the [PITH_FULL_IMAGE:figures/full_fig_p014_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Proton bead distribution along the S-S direction for several pressures. PIMD simulations at 50 K and the 200 K are [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Anharmonic phonon dispersions computed by PIMD simulations at [PITH_FULL_IMAGE:figures/full_fig_p016_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Static limit of the local 2-body phonon Green function [PITH_FULL_IMAGE:figures/full_fig_p017_11.png] view at source ↗

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

74 extracted references · 4 canonical work pages

  1. [1]

    Onnes, The Superconductivity of Mercury, Comm

    H. Onnes, The Superconductivity of Mercury, Comm. Phys. Lab. Univ., Leiden , 122 (1911)

    1911

  2. [2]

    N. W. Ashcroft, Metallic Hydrogen: A High- Temperature Superconductor?, Physical Review Letters 21, 1748–1749 (1968)

    1968

  3. [3]

    L. P. Gor’kov and V. Z. Kresin, Pressure and high-Tc superconductivity in sulfur hydrides, Scientific Reports 6, 10.1038/srep25608 (2016)

    work page doi:10.1038/srep25608 2016

  4. [4]

    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)

    2015

  5. [6]

    Mozaffari, D

    S. Mozaffari, D. Sun, V. S. Minkov, A. P. Drozdov, D. Knyazev, J. B. Betts, M. Einaga, K. Shimizu, M. I. Eremets, L. Balicas, and F. F. Balakirev, Superconduct- ing phase diagram of H3S under high magnetic fields, Nature Communications10, 2522 (2019)

    2019

  6. [7]

    V. S. Minkov, V. B. Prakapenka, E. Greenberg, and M. I. Eremets, A Boosted Critical Temperature of 166 K in Su- perconducting D3S Synthesized from Elemental Sulfur and Hydrogen, Angewandte Chemie International Edi- tion59, 18970 (2020)

    2020

  7. [8]

    Osmond, O

    I. Osmond, O. Moulding, S. Cross, T. Muramatsu, A. Brooks, O. Lord, T. Fedotenko, J. Buhot, and S. Friedemann, Clean-limit superconductivity inIm 3m H3Ssynthesized from sulfur and hydrogen donor ammo- nia borane, Phys. Rev. B105, L220502 (2022)

    2022

  8. [9]

    J. G. Bednorz and K. A. Muller, Possible hightc super- conductivity in the Ba-La-Cu-O system, Zeitschrift fur Physik B Condensed Matter64, 189–193 (1986)

    1986

  9. [10]

    L. Gao, Y. Y. Xue, F. Chen, Q. Xiong, R. L. Meng, D. Ramirez, C. W. Chu, J. H. Eggert, and H. K. Mao, Superconductivity up to 164 K inHgBa 2Ca m−1Cu mO2m+2+δ (m=1, 2, and 3) un- 6 der quasihydrostatic pressures, Physical Review B50, 4260–4263 (1994)

    1994

  10. [11]

    Bianconi and T

    A. Bianconi and T. Jarlborg, Superconductivity above the lowest earth temperature in pressurized sulfur hy- dride, EPL (Europhysics Letters)112, 37001 (2015)

    2015

  11. [12]

    Keimer, S

    B. Keimer, S. A. Kivelson, M. R. Norman, S. Uchida, and J. Zaanen, From quantum matter to high-temperature superconductivity in copper oxides, Nature518, 179 (2015)

    2015

  12. [13]

    Michon, C

    B. Michon, C. Girod, S. Badoux, J. Kačmarčík, Q. Ma, M. Dragomir, H. Dabkowska, B. Gaulin, J.-S. Zhou, S. Pyon,et al., Thermodynamic signatures of quantum criticality in cuprate superconductors, Nature567, 218 (2019)

    2019

  13. [14]

    K. L. Ngai, Two-Phonon Deformation Potential and Su- perconductivity in Degenerate Semiconductors, Physical Review Letters32, 215–218 (1974)

    1974

  14. [16]

    van der Marel, F

    D. van der Marel, F. Barantani, and C. W. Rischau, Possible mechanism for superconductivity in doped SrTiO3, Physical Review Research1, 10.1103/physrevre- search.1.013003 (2019)

    doi:10.1103/physrevre- 2019

  15. [17]

    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-Tc supercon- ductivity, Scientific reports4, 6968 (2014)

    2014

  16. [18]

    J. A. Flores-Livas, A. Sanna, and E. Gross, High temper- ature superconductivity in sulfur and selenium hydrides at high pressure, The European Physical Journal B89, 63 (2016)

    2016

  17. [19]

    A. F. Goncharov, S. S. Lobanov, V. B. Prakapenka, and E. Greenberg, Stable high-pressure phases in the H-S sys- tem determined by chemically reacting hydrogen and sul- fur, Phys. Rev. B95, 140101 (2017)

    2017

  18. [20]

    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, Nature532, 81 (2016)

    2016

  19. [21]

    Bianco, I

    R. Bianco, I. Errea, M. Calandra, and F. Mauri, High- pressure phase diagram of hydrogen and deuterium sul- fides from first principles: Structural and vibrational properties including quantum and anharmonic effects, Phys. Rev. B97, 214101 (2018)

    2018

  20. [22]

    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 Computational Materials 10, 10.1038/s41524-024-01239-0 (2024)

    work page doi:10.1038/s41524-024-01239-0 2024

  21. [23]

    A. D. Becke, Density-functional exchange-energy approx- imation with correct asymptotic behavior, Phys. Rev. A 38, 3098 (1988)

    1988

  22. [24]

    C. Lee, W. Yang, and R. G. Parr, Development of the Colle-Salvetti correlation-energy formula into a func- tional of the electron density, Physical review B37, 785 (1988)

    1988

  23. [25]

    Batatia, D

    I. Batatia, D. P. Kovacs, G. N. C. Simm, C. Ortner, and G. Csanyi, MACE: Higher Order Equivariant Message Passing Neural Networks for Fast and Accurate Force Fields, inAdvances in Neural Information Processing Systems, edited by A. H. Oh, A. Agarwal, D. Belgrave, and K. Cho (2022)

    2022

  24. [27]

    Ceriotti, J

    M. Ceriotti, J. More, and D. E. Manolopoulos, i-PI: A Python interface for ab initio path integral molecular dy- namics simulations, Computer Physics Communications 185, 1019–1026 (2014)

    2014

  25. [28]

    See Supplemental Material at [URL will be inserted by publisher] for additional information about the compu- tational details of PIMD and classical MD simulations, the MLIP training and validation, the bead-resolved and imaginary-time resolved estimators, and the 2-body Green function calculations (see references [4, 5, 15, 21– 27, 29, 30, 37–39, 55–59] therein)

  26. [29]

    Morresi, L

    T. Morresi, L. Paulatto, R. Vuilleumier, and M. Casula, Probing anharmonic phonons by quantum correlators: A path integral approach, The Journal of Chemical Physics 154, 224108 (2021)

    2021

  27. [30]

    Morresi, R

    T. Morresi, R. Vuilleumier, and M. Casula, Hydrogen phase-IV characterization by full account of quantum an- harmonicity, Physical Review B106, 054109 (2022)

    2022

  28. [31]

    Sachdev,Quantum Phase Transitions(Cambridge University Press, 2011)

    S. Sachdev,Quantum Phase Transitions(Cambridge University Press, 2011)

    2011

  29. [32]

    Kim, Y.-C

    D.-H. Kim, Y.-C. Lin, and H. Rieger, Path inte- gral Monte Carlo study of the interacting quantum double-well model: Quantum phase transition and phase diagram, Physical Review E75, 10.1103/phys- reve.75.016702 (2007)

    doi:10.1103/phys- 2007

  30. [33]

    Mondaini, S

    R. Mondaini, S. Tarat, and R. T. Scalettar, Universality and critical exponents of the fermion sign problem, Phys- ical Review B107, 10.1103/physrevb.107.245144 (2023)

    work page doi:10.1103/physrevb.107.245144 2023

  31. [34]

    F. F. Assaad and I. F. Herbut, Pinning the Order: The Nature of Quantum Criticality in the Hubbard Model on Honeycomb Lattice, Physical Review X3, 031010 (2013)

    2013

  32. [35]

    Šuntajs, J

    J. Šuntajs, J. Bonča, T. Prosen, and L. Vidmar, Ergod- icity breaking transition in finite disordered spin chains, Physical Review B102, 064207 (2020)

    2020

  33. [36]

    A. S. Aramthottil, T. Chanda, P. Sierant, and J. Za- krzewski, Finite-size scaling analysis of the many-body localization transition in quasiperiodic spin chains, Phys. Rev. B104, 214201 (2021)

    2021

  34. [37]

    Brézin, An investigation of finite size scaling, Journal de Physique43, 15 (1982)

    E. Brézin, An investigation of finite size scaling, Journal de Physique43, 15 (1982)

    1982

  35. [38]

    Aktekin, The finite-size scaling functions of the four- dimensional Ising model, Journal of Statistical Physics 104, 1397 (2001)

    N. Aktekin, The finite-size scaling functions of the four- dimensional Ising model, Journal of Statistical Physics 104, 1397 (2001)

    2001

  36. [40]

    Dangić, L

    Ð. Dangić, L. Monacelli, R. Bianco, F. Mauri, and I. Er- rea, Large impact of phonon lineshapes on the supercon- ductivity of solid hydrogen, Communications Physics7, 150 (2024)

    2024

  37. [41]

    Marsiglio, Eliashberg theory: A short review, Annals of Physics417, 168102 (2020)

    F. Marsiglio, Eliashberg theory: A short review, Annals of Physics417, 168102 (2020)

    2020

  38. [42]

    Marsiglio, Phonon Self-Energy Effects in Migdal- Eliashberg Theory, inElectron–Phonon Interaction in Oxide Superconductors, edited by R

    F. Marsiglio, Phonon Self-Energy Effects in Migdal- Eliashberg Theory, inElectron–Phonon Interaction in Oxide Superconductors, edited by R. Baquero (World Sci- entific, Singapore, 1991) p. 167

    1991

  39. [43]

    Esterlis, B

    I. Esterlis, B. Nosarzewski, E. W. Huang, B. Moritz, T. P. Devereaux, D. J. Scalapino, and S. A. Kivelson, Breakdown of the Migdal-Eliashberg theory: A determi- nant quantum Monte Carlo study, Physical Review B97, 140501 (2018)

    2018

  40. [44]

    W. Sano, T. Koretsune, T. Tadano, R. Akashi, and R. Arita, Effect of Van Hove singularities on high-Tc su- perconductivity inH 3S, Physical Review B93, 094525 (2016)

    2016

  41. [45]

    Lucrezi, P

    R. Lucrezi, P. P. Ferreira, S. Hajinazar, H. Mori, H. Paudyal, E. R. Margine, and C. Heil, Full-bandwidth anisotropic Migdal-Eliashberg theory and its application to superhydrides, Communications Physics7, 33 (2024)

    2024

  42. [46]

    Pietronero, S

    L. Pietronero, S. Strässler, and C. Grimaldi, Nonadi- abatic superconductivity. I. Vertex corrections for the electron-phonon interactions, Physical Review B52, 10516 (1995)

    1995

  43. [47]

    Grimaldi, L

    C. Grimaldi, L. Pietronero, and S. Strässler, Nonadia- batic superconductivity. II. Generalized Eliashberg equa- tions beyond Migdal’s theorem, Physical Review B52, 10530 (1995)

    1995

  44. [48]

    S. B. Mishra, H. Mori, and E. R. Margine, Electron– phonon vertex correction effect in superconducting h3s, npj Computational Materials11, 342 (2025)

    2025

  45. [49]

    Bianco and I

    R. Bianco and I. Errea, Non-perturbative theory of the electron-phonon coupling and its first-principles imple- mentation, arXiv preprint arXiv:2303.02621 (2023)

    arXiv 2023

  46. [50]

    Collignon, X

    C. Collignon, X. Lin, C. W. Rischau, B. Fauqué, and K. Behnia, Metallicity and superconductivity in doped strontium titanate, Annual Review of Condensed Matter Physics10, 25 (2019)

    2019

  47. [51]

    M. N. Gastiasoro, J. Ruhman, and R. M. Fernandes, Su- perconductivity in dilute SrTiO3: A review, Annals of Physics417, 168107 (2020)

    2020

  48. [52]

    P. A. Volkov, P. Chandra, and P. Coleman, Superconduc- tivity from energy fluctuations in dilute quantum critical polar metals, Nature communications13, 4599 (2022)

    2022

  49. [53]

    Capitani, B

    F. Capitani, B. Langerome, J.-B. Brubach, P. Roy, A. Drozdov, M. Eremets, E. Nicol, J. Carbotte, and T. Timusk, Spectroscopic evidence of a new energy scale for superconductivity inH 3S, Nature physics13, 859 (2017)

    2017

  50. [54]

    Ferroelectric quantum critical point in super- conducting hydrides: The case ofH3S

    M. Cherubini, A. Raghav, and M. Casula, Additional data for "Ferroelectric quantum critical point in super- conducting hydrides: The case ofH3S" (2026)

    2026

  51. [56]

    Giannozzi, S

    P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- cioni, I. Dabo, A. Dal Corso, S. de Gironcoli, S. Fabris, G. Fratesi, R. Gebauer, U. Gerstmann, C. Gougoussis, A. Kokalj, M. Lazzeri, L. Martin-Samos, N. Marzari, F. Mauri, R. Mazzarello, S. Paolini, A. Pasquarello, L. Paulatto, C. Sbraccia, S. S...

    2009

  52. [57]

    Giannozzi, O

    P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. Buongiorno Nardelli, M. Calandra, R. Car, C. Cavaz- zoni, D. Ceresoli, M. Cococcioni, N. Colonna, I. Carn- imeo, A. Dal Corso, S. de Gironcoli, P. Delugas, R. A. DiStasio, A. Ferretti, A. Floris, G. Fratesi, G. Fugallo, R. Gebauer, U. Gerstmann, F. Giustino, T. Gorni, J. Jia, M. Kawamura, H.-Y. Ko, A. Ko...

    2017

  53. [58]

    Gilmer, S

    J. Gilmer, S. S. Schoenholz, P. F. Riley, O. Vinyals, and G. E. Dahl, Neural Message Passing for Quantum Chem- istry, inProceedings of the 34th International Conference on Machine Learning, Proceedings of Machine Learning Research, Vol. 70, edited by D. Precup and Y. W. Teh (PMLR, 2017) pp. 1263–1272

    2017

  54. [59]

    beyond-Migdal

    M. M. Bronstein, J. Bruna, T. Cohen, and P. Veličković, Geometric Deep Learning: Grids, Groups, Graphs, Geodesics, and Gauges (2021). Supplemental Material: Ferroelectric quantum critical point in superconducting hydrides: The case ofH 3S Marco Cherubini, 1,∗ Abhishek Raghav, 1, 2 and Michele Casula 1,† 1Institut de Min´ eralogie, de Physique des Mat´ eri...

    2021

  55. [60]

    Taureau, M

    R. Taureau, M. Cherubini, T. Morresi, and M. Casula, Quantum symmetrization transition in superconducting sulfur hy- dride from quantum Monte Carlo and path integral molecular dynamics, npj Computational Materials10, 10.1038/s41524- 024-01239-0 (2024)

    doi:10.1038/s41524- 2024

  56. [61]

    Ceriotti, M

    M. Ceriotti, M. Parrinello, T. E. Markland, and D. E. Manolopoulos, Efficient stochastic thermostatting of path integral molecular dynamics, The Journal of Chemical Physics133, 10.1063/1.3489925 (2010)

    doi:10.1063/1.3489925 2010

  57. [62]

    Ceriotti, J

    M. Ceriotti, J. More, and D. E. Manolopoulos, i-PI: A Python interface for ab initio path integral molecular dynamics simulations, Computer Physics Communications185, 1019–1026 (2014)

    2014

  58. [63]

    Bianco, I

    R. Bianco, I. Errea, M. Calandra, and F. Mauri, High-pressure phase diagram of hydrogen and deuterium sulfides from first principles: Structural and vibrational properties including quantum and anharmonic effects, Phys. Rev. B97, 214101 (2018)

    2018

  59. [64]

    A. D. Becke, Density-functional exchange-energy approximation with correct asymptotic behavior, Phys. Rev. A38, 3098 (1988)

    1988

  60. [65]

    C. Lee, W. Yang, and R. G. Parr, Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density, Physical review B37, 785 (1988)

    1988

  61. [66]

    Giannozzi, S

    P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococcioni, I. Dabo, A. Dal Corso, S. de Gironcoli, S. Fabris, G. Fratesi, R. Gebauer, U. Gerstmann, C. Gougoussis, A. Kokalj, M. Lazzeri, L. Martin-Samos, N. Marzari, F. Mauri, R. Mazzarello, S. Paolini, A. Pasquarello, L. Paulatto, C. Sbraccia, S. Sca...

    2009

  62. [67]

    Giannozzi, O

    P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. Buongiorno Nardelli, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, M. Cococcioni, N. Colonna, I. Carnimeo, A. Dal Corso, S. de Gironcoli, P. Delugas, R. A. DiStasio, A. Fer- retti, A. Floris, G. Fratesi, G. Fugallo, R. Gebauer, U. Gerstmann, F. Giustino, T. Gorni, J. Jia, M. Kawamura, H.-Y. Ko, A. Koka...

    2017

  63. [68]

    Batatia, D

    I. Batatia, D. P. Kovacs, G. N. C. Simm, C. Ortner, and G. Csanyi, MACE: Higher Order Equivariant Message Passing Neural Networks for Fast and Accurate Force Fields, inAdvances in Neural Information Processing Systems, edited by 11 A. H. Oh, A. Agarwal, D. Belgrave, and K. Cho (2022)

    2022

  64. [69]

    Batatia, S

    I. Batatia, S. Batzner, D. P. Kov´ acs, A. Musaelian, G. N. C. Simm, R. Drautz, C. Ortner, B. Kozinsky, and G. Cs´ anyi, The Design Space of E(3)-Equivariant Atom-Centered Interatomic Potentials (2022), arXiv:2205.06643

    Pith/arXiv arXiv 2022

  65. [70]

    Gilmer, S

    J. Gilmer, S. S. Schoenholz, P. F. Riley, O. Vinyals, and G. E. Dahl, Neural Message Passing for Quantum Chemistry, inProceedings of the 34th International Conference on Machine Learning, Proceedings of Machine Learning Research, Vol. 70, edited by D. Precup and Y. W. Teh (PMLR, 2017) pp. 1263–1272

    2017

  66. [71]

    M. M. Bronstein, J. Bruna, T. Cohen, and P. Veliˇ ckovi´ c, Geometric Deep Learning: Grids, Groups, Graphs, Geodesics, and Gauges (2021)

    2021

  67. [72]

    Br´ ezin, An investigation of finite size scaling, Journal de Physique43, 15 (1982)

    E. Br´ ezin, An investigation of finite size scaling, Journal de Physique43, 15 (1982)

    1982

  68. [73]

    Aktekin, The finite-size scaling functions of the four-dimensional Ising model, Journal of Statistical Physics104, 1397 (2001)

    N. Aktekin, The finite-size scaling functions of the four-dimensional Ising model, Journal of Statistical Physics104, 1397 (2001)

    2001

  69. [74]

    Kenna, Finite size scaling forO(N)φ4-theory at the upper critical dimension, Nuclear Physics B691, 292 (2004)

    R. Kenna, Finite size scaling forO(N)φ4-theory at the upper critical dimension, Nuclear Physics B691, 292 (2004)

    2004

  70. [75]

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

    2015

  71. [76]

    Einaga, M

    M. Einaga, M. Sakata, T. Ishikawa, K. Shimizu, M. I. Eremets, A. P. Drozdov, I. A. Troyan, N. Hirao, and Y. Ohishi, Crystal structure of the superconducting phase of sulfur hydride, Nature Physics12, 835 (2016)

    2016

  72. [77]

    Morresi, L

    T. Morresi, L. Paulatto, R. Vuilleumier, and M. Casula, Probing anharmonic phonons by quantum correlators: A path integral approach, The Journal of Chemical Physics154, 224108 (2021)

    2021

  73. [78]

    Morresi, R

    T. Morresi, R. Vuilleumier, and M. Casula, Hydrogen phase-IV characterization by full account of quantum anharmonicity, Physical Review B106, 054109 (2022)

    2022

  74. [79]

    D. E. Kiselov and M. V. Feigel’man, Theory of superconductivity due to Ngai’s mechanism in lightly doped SrTiO 3, Physical Review B104, 10.1103/physrevb.104.l220506 (2021)

    work page doi:10.1103/physrevb.104.l220506 2021