REVIEW 3 major objections 5 minor 43 references
This paper introduces ah-SCDFT, a self-consistent scheme that adds anharmonic lattice corrections to superconducting density functional theory, and shows that for MgB2 it raises the predicted critical temperature from 36 K to 39 K, matching
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 →
A self-consistent loop between SCDFT and quasi-harmonically renormalized structures raises the predicted Tc of MgB2 from 36 K to 39 K and shifts the pressure turnaround to ~200 GPa.
T0 review reviewed 2026-08-01 challenge →
load-bearing objection A useful fixed-point extension of SCDFT with QHA-derived thermal geometries; the MgB2 numbers are plausible, but the 'anharmonic' label oversells what is a quasi-harmonic correction. the 3 major comments →
Ah-SCDFT:A general approach for superconductivity with an-harmonic corrections
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
Core claim
The central claim is that the temperature-dependent equilibrium geometry, generated by minimizing the Helmholtz free energy (electronic total energy plus harmonic phonon free energy) over strained lattices, is sufficient to capture the anharmonic lattice renormalization that matters for superconductivity in MgB2. Iterating between superconducting density functional theory and this renormalized geometry converges at Tc = 39 K, versus 36 K for the 0 K harmonic geometry. The mechanism is identified: the renormalized lattice raises the electronic density of states at the Fermi level and softens the strongly coupled in-plane boron phonon modes near the A point, increasing the electron-phonon coup
What carries the argument
The load-bearing object is a self-consistent iteration loop: (1) compute Tc with standard superconducting density functional theory at the 0 K harmonic geometry; (2) at that temperature, minimize the quasi-harmonic Helmholtz free energy over symmetric strains to obtain a thermally expanded, phonon-renormalized geometry; (3) recompute Tc on this new geometry; repeat until consecutive Tc values differ by less than 0.01 K. The quasi-harmonic approximation used here includes thermal expansion and volume-dependent harmonic phonon frequencies, but not explicit phonon-phonon interactions.
Load-bearing premise
Everything rests on the assumption that quasi-harmonic free-energy minimization (thermal expansion plus harmonic phonon frequencies at each volume) captures the anharmonic effects that matter for superconductivity, while leaving out explicit phonon-phonon interactions and zero-point anharmonicity.
What would settle it
Measure or fully anharmonically calculate the temperature dependence of the strongly coupled in-plane boron phonons in MgB2 near the A point between 0 K and 39 K; if the softening and the accompanying density-of-states increase are absent or much smaller than the quasi-harmonic prediction, the Tc enhancement to 39 K is not caused by the mechanism claimed. Equivalently, a high-pressure experiment mapping Tc between 100 and 250 GPa that finds no turning point near 200 GPa would falsify the pressure claim.
If this is right
- If the renormalization loop is correct, harmonic-only SCDFT systematically underestimates Tc in materials with soft phonons or sizable thermal expansion, and the error is fixable without new physics.
- The method transfers directly to other phonon-mediated superconductors, including hydrides under pressure where anharmonic effects are known to be large, at modest extra cost.
- The predicted pressure dependence for MgB2 implies that the nonmonotonic Tc(p) is real but occurs at pressures beyond the previously claimed 100 GPa turning point, offering a target for high-pressure experiments.
- The mechanism (DOS increase plus phonon softening) makes concrete predictions for temperature-dependent phonon spectra and electronic structure that can be tested by neutron scattering or angle-resolved photoemission.
- Since the loop converges in a few iterations, the anharmonic correction can be added to any existing SCDFT workflow without algorithmic changes.
Where Pith is reading between the lines
- The term 'anharmonic' in this scheme really means quasi-harmonic thermal renormalization; explicit anharmonic phonon-phonon interactions and zero-point anharmonicity are left out, so if those effects are substantial the agreement at 39 K could be fortuitous rather than systematic.
- The same iteration could be applied to materials that are dynamically unstable in the harmonic approximation, where the 0 K starting geometry is ill-defined; the method as stated may need a self-consistent phonon starting point to generalize.
- A sharper test would be to run the loop on a superconductor like H3S, where anharmonicity reportedly changes Tc substantially, and compare the quasi-harmonic-only result to fully anharmonic calculations.
- If the thermal renormalization explanation holds, temperature-dependent lattice constant measurements on MgB2 should show the same expansion and phonon softening used in the calculation, a straightforward experimental check.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper proposes a computational workflow, 'ah-SCDFT', that iterates between standard SCDFT superconductivity calculations on the 0 K harmonic geometry and a temperature-dependent equilibrium geometry obtained from the quasi-harmonic approximation (QHA) at the currently predicted critical temperature. The authors apply this to MgB2: the harmonic SCDFT gives Tc = 36 K, while using the QHA-renormalized geometry at 39 K raises Tc to 39 K, matching the experimental value. Under pressure, the same procedure yields a nonmonotonic Tc(P) with a turning point shifted from ~100 GPa (previous harmonic calculations) to ~200 GPa, and gives Tc = 20 K at 25 GPa in agreement with experiment. The paper claims this establishes a general, computationally inexpensive way to include anharmonic effects in first-principles superconductivity predictions.
Significance. If the claims hold, the proposed workflow is practically useful: it requires no empirical μ*, uses widely available QHA codes and SCDFT implementations, and the fixed-point iteration between Tc and lattice geometry is a clear and easily implemented idea. The reproducible pipeline (QUANTUM ESPRESSO, thermo_pw, superconducting toolkit) and the explicit convergence criterion are strengths. However, the significance is currently limited by the mismatch between the title/abstract claim of 'anharmonic corrections' and the actual content of Eqs. (1)-(2), which is the quasi-harmonic approximation. The numerical improvement over harmonic SCDFT is modest (3 K) and is presented without error bars or sensitivity analysis. The pressure prediction at ~200 GPa is not benchmarked against experiment or explicit anharmonic calculations. The central idea is defensible, but the paper overstates its novelty and the evidence for it.
major comments (3)
- The anharmonic content of the proposed method is the quasi-harmonic approximation (QHA): Eq. (2) is the harmonic phonon free energy summed over phonon wavevectors and branches, with no cubic/quartic anharmonic terms, no phonon self-energy, and no explicit phonon-phonon interaction. The geometry at T_cn is obtained by minimizing this harmonic free energy over symmetric strains, which is thermal expansion within the QHA, not a calculation of anharmonic corrections. The abstract and Introduction claim 'systematically incorporates anharmonic corrections' and 'explicitly integrates anharmonic corrections into SCDFT'; these statements are not supported by Eqs. (1)-(2). The method should be described as SCDFT with a QHA-renormalized lattice, or the authors must benchmark the QHA against explicit anharmonic calculations (e.g., SSCHA) and show that the missing phonon-phonon terms are small for Mg
- The central numerical demonstration is a 3 K increase of Tc from 36 to 39 K caused by the QHA-renormalized geometry. No numerical uncertainty or sensitivity analysis is reported. SCDFT with the Sanna2020 kernel is already accurate to a few K for many materials, and the difference between 36 and 39 K is comparable to typical methodological error. The attribution 'in excellent agreement with experiment' therefore requires error bars or a sensitivity study (pseudopotentials, k/q grids, exchange-correlation functional, QHA strain sampling). Without this, the 39 K result may be coincidental. At minimum, the paper should report the convergence of Tc with respect to numerical settings and the effect of using different pseudopotentials, since the Methods section itself uses SG15 for superconductivity and PseudoDojo for geometry.
- The pressure claim—the turning point shifts from ~100 GPa (Singh) to ~200 GPa—is made without any experimental data above 40 GPa and without validating the QHA equation of state against measured lattice constants under pressure. The text also states MgB2 becomes structurally unstable above 340 GPa from imaginary phonons, but this is not connected to any experimental structural transition. Given that the method is intended to resolve discrepancies among prior reports, the authors should compare their computed P-V/V0 curve to experiment and assess QHA validity under strong compression (e.g., anharmonic frequency shifts versus volume, or explicit anharmonic calculations at representative pressures). Otherwise, the 'turning point' is a prediction that is currently untestable from the data shown.
minor comments (5)
- Equation (2) is garbled in the manuscript: the equals sign is missing and the sum over modes and the plus sign separating zero-point and thermal parts are unclear. Please rewrite it cleanly.
- Reference [41] is cited as 'Ma et al.' but the first author is Y. Wang. Please correct the citation to match the reference list.
- The x-axis label is missing and the definition of T_cn (n=0,1,2,...) is not stated. Please add an axis label and define the iteration index.
- The data availability statement says the data are not publicly available because it is 'not technically feasible'—this is unusual for a computational paper. Please deposit the relevant input/output files or provide a more specific justification.
- The title and abstract use 'an-harmonic' with a hyphen inconsistently; the standard term 'anharmonic' should be used throughout. Also, phrases such as '12^3 Monkhorst-Pack k-point mesh' should be written as 12×12×12 for clarity.
Circularity Check
No circularity: the geometry-Tc fixed-point loop is a self-consistency condition, not a fitted or definitional reduction.
full rationale
The paper's central loop—compute SCDFT Tc at the harmonic 0 K geometry, re-evaluate the QHA geometry at that Tc, recompute Tc, and iterate—is a legitimate fixed-point equation T = F(T), not a circular derivation. The geometry entering F is obtained independently by minimizing the harmonic phonon free energy (Eqs. 1–2); it contains no superconducting order parameter and no information about the target Tc beyond the temperature argument in the Bose occupation factors. The SCDFT gap equation is an independent functional of that geometry, using the Sanna2020 kernel and external codes (QUANTUM ESPRESSO, thermo_pw, SG15/PseudoDojo pseudopotentials). No parameter such as μ* is fitted to reproduce the experimental 39 K; experiment is used only as a posteriori comparison. The pressure-dependent result follows from the same parameter-free loop. The only substantive concern is that QHA is not full anharmonicity (no explicit phonon–phonon self-energy), but that is a physical-completeness and validation issue, not circularity: the derivation does not reduce to its input by construction, and the paper does not invoke any self-citation as load-bearing evidence. Hence no circular step is identified.
Axiom & Free-Parameter Ledger
axioms (4)
- domain assumption Quasi-harmonic free energy (Eqs. 1-2), evaluated with harmonic phonons under symmetric strain, approximates the full anharmonic free energy relevant to superconductivity.
- domain assumption Electronic contribution to the Helmholtz free energy is negligible for the geometry optimization.
- domain assumption SCDFT with the Sanna2020 exchange-correlation kernel reliably predicts Tc when given the correct geometry.
- domain assumption The fixed-point iteration between Tc and the QHA geometry at T=Tc converges to the physical critical temperature.
Cite this review
Pith. "Pith review of Ah-SCDFT:A general approach for superconductivity with an-harmonic corrections." pith.science (2026). https://pith.science/paper/BED666GB
@misc{pith2026260717759,
author = {Pith},
title = {Pith review of: Ah-SCDFT:A general approach for superconductivity with an-harmonic corrections},
year = {2026},
howpublished = {\url{https://pith.science/paper/BED666GB}},
note = {Machine review of arXiv:2607.17759}
}
read the original abstract
First-principles studies of superconductivity often neglect anharmonic effects (AHE), despite their crucial role in achieving quantitative accuracy in many materials. To bridge this gap, we introduce a general computational approach, termed anharmonic superconducting density functional theory (ah-SCDFT) which systematically incorporates anharmonic corrections into standard SCDFT. This approach allows for high-fidelity predictions of superconducting properties with only a modest increase in computational cost for a limited number of superconducting calculation convergence steps. We demonstrate the effectiveness and reliability of ah-SCDFT by applying it to the prototypical superconductor MgB2, accurately reproducing its superconducting behavior under both ambient conditions and applied pressure in excellent agreement with experiment. Our results establish ah-SCDFT as a powerful, efficient, and broadly applicable approach for quantitatively reliable studies of superconductivity and a promising tool for the prediction of new superconducting materials.
Figures
Reference graph
Works this paper leans on
-
[1]
Bardeen, L
J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Theory of superconductivity, Phys. Rev. 108, 1175 (1957)
1957
-
[2]
J. P. Carbotte, Properties of boson- exchange superconductors, Rev. Mod. Phys. 62, 1027 (1990)
1990
-
[3]
E. R. Margine and F. Giustino, Anisotropic migdal - eliashberg theory using wannier functions, Phys. Rev. B 87, 024505 (2013)
2013
-
[4]
W. L. McMillan, Transition temperature of strong - coupled superconductors, Phys. Rev. 167, 331 (1968)
1968
-
[5]
R. C. Dynes, McMillan's equation and the T c of superconductors, Solid State Commun. 10, 615 (1972)
1972
-
[6]
P. B. Allen and R. C. Dynes, Transition temperature of strong-coupled superconductors reanalyzed, Phys. Rev. B 12, 905 (1975)
1975
-
[7]
L. N. Oliveira, E. K. U. Gross, and W. Kohn, Density - functional theory for superconductors, Phys. Rev. Lett. 60, 2430 (1988)
1988
-
[8]
Lüders, M
M. Lüders, M. A. L. Marques, N. N. Lathiotakis, A. Floris, G. Profeta, L. Fast, A. Continenza, S. Massidda, and E. K. U. Gross, Ab initio theory of superconductivity. I. Density functional formalism and approximate functionals, Phys. Rev. B 72, 024545 (2005)
2005
-
[9]
M. A. L. Marques, M. Lüders, N. N. Lathiotakis, G. Profeta, A. Floris, L. Fast, A. Continenza, E. K. U. Gross, and S. Massidda, Ab initio theory of superconductivity. II. Application to elemental metals, Phys. Rev. B 72, 024546 (2005)
2005
-
[10]
Sanna, C
A. Sanna, C. Pellegrini, and E. K. U. Gross, Combining eliashberg theory with density functional theory for the accurate prediction of superconducting transition temperatures and gap functions, Phys. Rev. Lett. 125, 057001 (2020)
2020
-
[11]
Kawamura, Y
M. Kawamura, Y . Hizume, and T. Ozaki, Benchmark of density functional theory for superconductors in elemental materials, Phys. Rev. B 101, 134511 (2020)
2020
-
[12]
Leibfried and W
G. Leibfried and W. Ludwig, Theory of anharmonic effects in crystals, in Solid State Physics (Academic Press, New York, 1961), V ol. 12, p. 275
1961
-
[13]
Errea, M
I. Errea, M. Calandra, and F. Mauri, Anharmonic free energies and phonon dispersions from the stochastic self-consistent harmonic approximation: Application to platinum and palladium hydrides, Phys. Rev. B 89, 064302 (2014)
2014
-
[14]
Bianco, I
R. Bianco, I. Errea, L. Paulatto, M. Calandra, and F. Mauri, Second -order structural phase transitions, free energy curvature, and temperature- dependent anharmonic phonons in the self -consistent harmonic approximation: Theory and stochastic implementation, Phys. Rev. B 96, 014111 (2017)
2017
-
[15]
Monacelli, I
L. Monacelli, I. Errea, M. Calandra, and F. Mauri, Pressure and stress tensor of complex anharmonic crystals within the stochastic self -consistent harmonic approximation, Phys. Rev. B 98, 024106 (2018)
2018
-
[16]
Monacelli, R
L. Monacelli, R. Bianco, M. Cherubini, M. Calandra, I. Errea, and F. Mauri, The stochastic self -consistent harmonic approximation: calculating vibrational properties of materials with full quantum and anharmonic effects, J. Phys.: Condens. Matter 33, 363001 (2021)
2021
-
[17]
S. Qiao, Y . Han, K. Jiang, N. Jiao, C. Si, H. Lu, and P . Zhang, Diverse and tunable charge density waves, superconductivity, and band topology in two- dimensional 1T-M 2S (M = Zr, Hf), Adv. Funct. Mater. 35, 2502257 (2025)
2025
-
[18]
X. Wang, Z. Tang, J. Li, C. He, M. Chen, C. Tang, and T. Ouyang, Phonon thermal transport in MTeX 4 (M = Zr, Hf; X = S, Se) monolayers: Role of antibonding states and higher -order phonon anharmonicity, Phys. Rev. Appl. 24, 044023 (2025)
2025
-
[19]
Jiang, Y
X. Jiang, Y . Li, W. Guo, Q. Ren, G. Tang, Z.-Z. Luo, and J. Hong, Low and anisotropic thermal conductivity in mixed-valent Sn 2S3, Phys. Rev. B 112, 144304 (2025)
2025
-
[20]
Z h o u , C
W . Z h o u , C . K e , a n d S . L i u , T w o-dimensional ferroelectric crystal with temperature-invariant ultralow thermal conductivity, Phys. Rev. B 112, 134307 (2025)
2025
-
[21]
G. He, L. Peis, E. F. Cuddy, Z. Zhao, D. Li, Y . Zhang, R. Stumberger, B. Moritz, H. Yang, H. Gao, et al. , Anharmonic strong-coupling effects at the origin of the charge density wave in CsV 3Sb5, Nat. Commun. 15, 1895 (2024)
2024
-
[22]
S. B. Mishra, H. Mori, and E. R. Margine, Electron – phonon vertex correction effect in superconducting H3S, npj Comput. Mater. 11, 342 (2025)
2025
-
[23]
S. Lu, R. Akashi, M. Kawamura, and S. Tsuneyuki, Assessing the possible superconductivity in the doped perovskite hydride KMgH 3: Effects of lattice anharmonicity and spin fluctuations, Phys. Rev. B 111 , 134516 (2025)
2025
-
[24]
Floris, G
A. Floris, G. Profeta, N. N. Lathiotakis, M. Lüders, M. A. L. Marques, C. Franchini, E. K. U. Gross, A. Continenza, and S. Massidda, Superconducting properties of MgB2 from first principles, Phys. Rev. Lett. 94, 037004 (2005)
2005
-
[25]
Floris, A
A. Floris, A. Sanna, M. Lüders, G. Profeta, N. N. Lathiotakis, M. A. L. Marques, C. Franchini, E. K. U. Gross, A. Continenza, and S. Massidda, Superconducting properties of MgB 2 from first principles, Physica C 456, 45 (2007)
2007
-
[26]
J. A. Flores -Livas and A. Sanna, Superconductivity in intercalated group-IV honeycomb structures, Phys. Rev. B 91, 054508 (2015)
2015
-
[27]
Giannozzi, S
P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococcioni, I. Dabo, et al., QUANTUM ESPRESSO: a modular and open -source software project for quantum simulations of materials, J. Phys.: Condens. Matter 21, 395502 (2009)
2009
-
[28]
Giannozzi, O
P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. 7 Buongiorno Nardelli, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, M. Cococcioni, et al. , Advanced capabilities for materials modelling with Quantum ESPRESSO, J. Phys.: Condens. Matter 29, 465901 (2017)
2017
-
[29]
Giannozzi, O
P. Giannozzi, O. Baseggio, P. Bonfà, D. Brunato, R. Car, I. Carnimeo, C. Cavazzoni, S. d. Gironcoli, P. Delugas, F. F. Ruffino, et al., Quantum ESPRESSO toward the exascale, J. Chem. Phys. 152, 154105 (2020)
2020
-
[30]
Kawamura, R
M. Kawamura, R. Akashi, and S. Tsuneyuki, Anisotropic superconducting gaps in YNi 2B2C: A first- principles investigation, Phys. Rev. B 95, 054506 (2017)
2017
-
[31]
Thermo_pw: thermoelastic properties calculator, https://dalcorso.github.io/thermo_pw/
-
[32]
D. R. Hamann, Optimized norm -conserving Vanderbilt pseudopotentials, Phy. Rev. B 88, 085117 (2013)
2013
-
[33]
Schlipf and F
M. Schlipf and F. Gygi, Optimization algorithm for the generation of ONCV pseudopotentials, Comput. Phys. Commun. 196, 36 (2015)
2015
-
[34]
Prandini, A
G. Prandini, A. Marrazzo, I. E. Castelli, N. Mounet, and N. Marzari, Precision and efficiency in solid -state pseudopotential calculations, npj Comput. Mater. 4 , 72 (2018)
2018
-
[35]
Kawamura, Y
M. Kawamura, Y . Gohda, and S. Tsuneyuki, Improved tetrahedron method for the Brillouin- zone integration applicable to response functions, Phys. Rev. B 89, 094515 (2014)
2014
-
[36]
PseudoDojo: pseudopotential database, https://www.pseudo-dojo.org/
-
[37]
Nagamatsu, N
J. Nagamatsu, N. Nakagawa, T. Muranaka, Y . Zenitani, and J. Akimitsu, Superconductivity at 39 K in magnesium diboride, Nature (London) 410, 63 (2001)
2001
-
[38]
L. Liu, F. Peng, P. Song, X. Liu, L. Zhang, X. Huang, C. Niu, C. Liu, W. Zhang, Y . Jia, et al., Generic rules for achieving room -temperature superconductivity in ternary hydrides with clathrate structures, P hys. Rev. B 107, L020504 (2023)
2023
-
[39]
Zhang, H
P. Zhang, H. Guan, H. Li, X. Zhong, R. J. Hemley, and H. Liu, High -temperature superconductivity in quinary clathrate hydrides under pressure, Phys. Rev. B 112 , 054516 (2025)
2025
-
[40]
P. P. Singh, From E 2g to other modes: Effects of pressure on electron-phonon interaction in MgB2, Phys. Rev. Lett. 97, 247002 (2006)
2006
-
[41]
Y . Wang, J. Lv, Y . Ma, T. Cui, and G. Zou, Superconductivity of MgB2 under ultrahigh pressure: A first-principles study, Phys. Rev. B 80, 092505 (2009)
2009
-
[42]
Monteverde, M
M. Monteverde, M. Núñez-Regueiro, N. Rogado, K. A. Regan, M. A. Hayward, T. He, S. M. Loureiro, and R. J. Cava, Pressure dependence of the superconducting transition temperature of magnesium diboride, Science 292, 75 (2001)
2001
-
[43]
Z. D. Gerbi and P. Singh, Effect of hydrostatic pressure on superconductivity of MgB 2, AIP Adv. 9, 095038 (2019)
2019
This paper was first reviewed by deepseek-v4-flash on August 1, 2026.
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.