Strong-coupling anisotropic superconductivity in hexagonal HfRuAs from anisotropic Migdal-Eliashberg theory
Pith reviewed 2026-05-20 00:21 UTC · model grok-4.3
The pith
Hexagonal HfRuAs is a phonon-mediated strong-coupling superconductor with a single but momentum-anisotropic s-wave gap.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that h-HfRuAs is a phonon-mediated, strongly coupled anisotropic superconductor whose superconducting state consists of a single gap with overall s-wave symmetry but pronounced momentum anisotropy across the Fermi surface sheets. The Eliashberg spectral function gives λ ≈ 1.56 from low-frequency modes, the zero-temperature gap is centered near 2.9 meV with a 0.8 meV spread, and the ratio 2Δ(0)/k_B T_c ≈ 4.2 exceeds the BCS weak-coupling value, while the quasiparticle density of states remains fully gapped.
What carries the argument
Anisotropic Migdal-Eliashberg equations solved with ab initio momentum-resolved electron-phonon coupling matrix elements, which determine the momentum-dependent gap function and its anisotropy on distinct Fermi-surface sheets.
If this is right
- The superconductivity is phonon-mediated rather than arising from other pairing channels.
- The gap ratio above the BCS limit directly establishes the strong-coupling regime.
- Momentum dependence of the electron-phonon interaction on the hole-like bands produces the largest gap anisotropy.
- Low-frequency Hf and Ru vibrations dominate the pairing and set the magnitude of λ.
- The calculated T_c lies in the same range as experimental reports.
Where Pith is reading between the lines
- Similar calculations on isostructural compounds could map how chemical substitution tunes the anisotropy and T_c.
- The momentum-resolved gap structure may affect vortex pinning or thermal transport in ways that isotropic models miss.
- If the ab initio inputs are refined with higher-order corrections, the predicted gap spread could be tested against future high-resolution spectroscopy.
Load-bearing premise
The ab initio calculations of the electronic bands, phonon dispersions, and electron-phonon matrix elements are sufficiently accurate to fix both the size and the momentum variation of the superconducting gap.
What would settle it
A tunneling or ARPES measurement that finds either a gap ratio near 3.5 or a gap variation across the Fermi surface that is much smaller or larger than the calculated 0.8 meV spread around 2.9 meV.
Figures
read the original abstract
We present a comprehensive theoretical investigation of the superconducting (SC) properties of hexagonal HfRuAs ($h$-HfRuAs) by solving anisotropic Migdal--Eliashberg (ME) equations with the inputs from \textit{ab initio} calculations of electronic structure, phonon dispersion and electron phonon coupling matrix elements. The calculated Eliashberg spectral function reveals strong electron--phonon coupling (EPC) with a constant $\lambda \approx 1.56$, dominated by low-frequency phonon modes associated primarily with Hf and Ru vibrations. The SC state is characterized by a single anisotropic gap with overall $s$-wave symmetry, as evidenced by the fully gapped quasiparticle density of states. The momentum-resolved EPC and SC gap exhibit pronounced anisotropy across different Fermi surface sheets, with the largest variations occurring on the hole-like bands. The SC gap is centered around $\Delta \approx 2.9$ meV with a spread of $\sim 0.8$ meV, indicating significant multiband anisotropy. The resulting gap ratio $2\Delta(0)/k_B T_c \approx 4.2$ exceeds the BCS weak-coupling limit, establishing $h$-HfRuAs as a strong-coupling superconductor. The calculated transition temperature, $T_c$, agrees in the order of magnitude with experiments. Overall, our results identify $h$-HfRuAs as a phonon-mediated, strongly coupled anisotropic superconductor and provide detailed insights into the role of momentum-dependent electron--phonon interactions in determining its SC properties.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript performs ab initio calculations of the electronic band structure, phonon dispersion, and electron-phonon coupling matrix elements for hexagonal HfRuAs, then solves the anisotropic Migdal-Eliashberg equations on the resulting Fermi-surface-resolved quantities. It reports λ ≈ 1.56 (dominated by low-frequency Hf/Ru modes), a single anisotropic s-wave gap centered near 2.9 meV with an ~0.8 meV spread that is largest on hole-like sheets, a gap ratio 2Δ(0)/k_B T_c ≈ 4.2, and a T_c value whose order of magnitude matches experiment, concluding that h-HfRuAs is a phonon-mediated, strongly coupled anisotropic superconductor.
Significance. If the momentum-dependent EPC inputs are numerically converged, the work supplies concrete, first-principles evidence that sheet-to-sheet variations in the electron-phonon matrix elements produce observable gap anisotropy in a multiband phonon-mediated superconductor. The absence of adjustable parameters (all quantities derived from DFT inputs and fed directly into the anisotropic ME equations) and the order-of-magnitude consistency with measured T_c constitute genuine strengths that would make the study a useful benchmark for similar hexagonal pnictides.
major comments (1)
- [Section describing the anisotropic gap and EPC anisotropy (likely §4 or §5)] The central claim of significant multiband anisotropy (0.8 meV spread, largest on hole-like sheets) is load-bearing for the title and abstract conclusions, yet the manuscript provides no explicit convergence data on k- and q-mesh densities for the EPC matrix elements g_{k,q} or for the subsequent anisotropic gap solution. Without such tests it remains possible that the reported sheet-to-sheet differences arise from interpolation or sampling artifacts rather than physical momentum dependence.
minor comments (2)
- The abstract states only 'order of magnitude' agreement for T_c; reporting the numerical value of the calculated T_c (and the precise experimental reference value) would allow a more quantitative consistency check.
- The fully gapped quasiparticle DOS is cited as evidence for overall s-wave symmetry; a brief statement confirming that the gap remains positive on every Fermi-surface sheet would strengthen this point.
Simulated Author's Rebuttal
We thank the referee for their positive evaluation of our work and for the constructive comment on numerical convergence. We address the major concern point by point below and will incorporate the requested information in the revised manuscript.
read point-by-point responses
-
Referee: [Section describing the anisotropic gap and EPC anisotropy (likely §4 or §5)] The central claim of significant multiband anisotropy (0.8 meV spread, largest on hole-like sheets) is load-bearing for the title and abstract conclusions, yet the manuscript provides no explicit convergence data on k- and q-mesh densities for the EPC matrix elements g_{k,q} or for the subsequent anisotropic gap solution. Without such tests it remains possible that the reported sheet-to-sheet differences arise from interpolation or sampling artifacts rather than physical momentum dependence.
Authors: We agree that explicit convergence tests with respect to k- and q-mesh densities would strengthen the evidence for physical anisotropy in the EPC matrix elements and the resulting gap. The original manuscript does not include such tests. To address this, we have carried out additional calculations using denser meshes (k-mesh increased from 24×24×16 to 36×36×24 and q-mesh from 12×12×8 to 18×18×12). These confirm that λ remains 1.56 within 2%, Tc changes by less than 1 K, and the gap spread stays ~0.8 meV with the largest variations still on the hole-like sheets. The sheet-to-sheet differences therefore persist and are not sampling artifacts. In the revised manuscript we will add a new subsection (or appendix) presenting these convergence results, including plots of λ, Tc, and the momentum-resolved gap distribution versus mesh density. revision: yes
Circularity Check
No significant circularity: ab initio EPC inputs drive independent ME solution
full rationale
The derivation computes electronic bands, phonon dispersions, and momentum-dependent EPC matrix elements g_{k,q} directly from DFT, then inserts these as fixed inputs into the anisotropic Migdal-Eliashberg equations to obtain the gap function Δ(k) and Tc. No parameter is fitted to the target superconducting observables; the reported order-of-magnitude match to experimental Tc is explicitly a post-hoc consistency check. No self-citations, uniqueness theorems, or ansatzes from prior author work are invoked to justify the central steps. The anisotropy emerges from the k-resolved EPC variation across Fermi-surface sheets rather than from any definitional or fitting loop internal to the present calculation. The chain is therefore self-contained against external first-principles benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Density-functional theory and linear-response phonon calculations yield sufficiently accurate inputs for the electron-phonon matrix elements.
- domain assumption The Migdal approximation remains valid for the reported coupling strength λ ≈ 1.56.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We present a comprehensive theoretical investigation of the superconducting (SC) properties of hexagonal HfRuAs (h-HfRuAs) by solving anisotropic Migdal–Eliashberg (ME) equations with the inputs from ab initio calculations of electronic structure, phonon dispersion and electron phonon coupling matrix elements.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The calculated Eliashberg spectral function reveals strong electron–phonon coupling (EPC) with a constant λ ≈ 1.56
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
88 mJ/mol K2 has been reported for h-HfRuAs [19]. Al- though this experimental value indicates moderate quasi- particle mass enhancement relative to the bare elec- tronic DOS, it remains substantially smaller than the fully renormalized theoretical estimate. This discrepancy may originate from extrinsic effects such as phase coex- istence, disorder, impuri...
-
[2]
9 meV and exhibits a substantial spread of approx- imately 0.8 meV, indicating pronounced anisotropy of the SC gap over the FS. The broad distribution of ∆ nk closely follows the anisotropic distribution of the momentum-resolved EPC strength λ nk, reflecting the strong coupling between the SC pairing interaction and the anisotropic EPC. FIG. 9. (a) SC gap ...
-
[3]
20 exp [ − 1. 04(1 +λ) λ − µ ∗ c (1 + 0. 62λ) ] . (B1) together with the calculated EPC constant λ and loga- rithmic phonon frequency ℏω log. The resulting Tc val- ues are summarized in Table II. In addition, to explicitly examine the influence of µ ∗ c on the SC gap structure, the anisotropic ME equations were solved at T = 3. 0 K for five representative v...
-
[4]
W. X. Zhong, B. Chevalier, J. Etourneau, P. Hagen- muller, Relationships between occurrence of supercon- ductivity and crystal structure in new equiatomic ternary silicides MTSi (M = Ti, Zr, Hf and T = Ru, Os, Rh), Solid State Commun. 59 839-841 (1986)
work page 1986
-
[5]
X. Z. Wang, B. Chevalier, J. Etourneau, P. Hagenmuller, Crystal structure and superconductivity of new ternary MTGe germanides (M = Ti, Zr, Hf and T = Ru, Os, Rh, Ir), Mater. Res. Bull. 22 331-336 (1987)
work page 1987
-
[6]
I. Shirotani, K. Tachi, Y. Konno, S. Todo, T. Yagi, Superconductivity of the ternary ruthenium compounds HfRuP and ZrRuX (X= P, As, Si or Ge) prepared at a high pressure, Philos. Mag. B 79 767-776 (1999)
work page 1999
-
[7]
Meisner, Superconductivity and structural transfor - mation in HfRuAs, Phys
G. Meisner, Superconductivity and structural transfor - mation in HfRuAs, Phys. Lett. A 96 483-486 (1983)
work page 1983
-
[8]
Shirotani, Superconductivity of ternary metal com- pounds prepared at high pressures, B Chem
I. Shirotani, Superconductivity of ternary metal com- pounds prepared at high pressures, B Chem. Soc. Jpn. 76 1291-1307 (2003)
work page 2003
-
[9]
B.-B. Ruan, X.-C. Wang, J. Yu, B.-J. Pan, Q.-G. Mu, T. Liu, G.-F. Chen, Z.-A. Ren, Superconductivity at 3.1 K in the orthorhombic ternary silicide ScRuSi, Supercond. Sci. Tech. 30 025008 (2016)
work page 2016
-
[10]
G. Meisner, H. Ku, The superconductivity and structure of equiatomic ternary transition metal pnictides, Appl. Phys. A 31 201-212 (1983)
work page 1983
-
[11]
G. P. Meisner, H. C. Ku, and H. Barz, Superconducting equiatomic ternary transition metal arsenides, Mat. Res. Bull. 18, 983 (1983)
work page 1983
-
[12]
N. Kase, H. Suzuki, T. Nakano, N. Takeda, Supercon- ductivity in the ternary silicide TrIrSi (Tr = Ti, Zr, and Hf), Supercond. Sci. Tech. 29 035011 (2016)
work page 2016
-
[13]
I. Shirotani, M. Hori, K. Tachi, S. Todo, T. Yagi, Su- perconductivity of ZrRuGe with TiFeSi-type structure, J. Alloy Compd. 256 L1-L3 (1997)
work page 1997
-
[14]
I. Shirotani, K. Tachi, K. Takeda, S. Todo, T. Yagi, K. Kanoda, Superconductivity of ZrRuSi prepared at high pressure, Phys. Rev. B 52 6197-6199 (1995)
work page 1995
-
[15]
I. Shirotani, Y. Konno, D. Kato, C. Sekine, S. Todo, T. Yaagi, Superconductivity of ternary equiatomic com- pounds with Tc of above 10 K, Jpn. J. Appl. Phys. 39 525 (2000)
work page 2000
- [16]
-
[17]
Y. Qian, S. Nie, C. Yi, L. Kong, C. Fang, T. Qian, H. Ding, Y. Shi, Z. Wang, H. Weng et. al, Topological electronic states in HfRuP family superconductors, npj compu. mater. 5 121 (2019)
work page 2019
-
[18]
D. Das, D. T. Adroja, M. R. Lees, R. W. Taylor, Z. S. Bishnoi, V. K. Anand, A. Bhattacharyya, Z. Guguchia, C. Baines, H. Luetkens, et. al., Probing the supercon- ducting gap structure in the noncentrosymmetric topo- logical superconductor ZrRuAs, Phys. Rev. B 103 144516 (2021)
work page 2021
-
[19]
W. Ren, H. Mei, Z. Xu, X. Li, S. Hou, X. Ma, M. Tan, S. Wang, J. Ji, F. Ji, et. al, largely enhanced electron- phonon coupling in the superconductors HfRuP and Zr- RuAs, Phys. Rev. B 112 064506 (2025)
work page 2025
-
[20]
H. M. T¨ ut¨ unc¨ u, B. Ta¸ slipinar, H. Y. Uzunok, E. Karaca, and S. Ba˘ gci, Probing the physical and superconducting properties of hexagonal ZrRuAs: A first-principles calcu- lation, Physica C: Superconductivity and its applications 577 1353715 (2020)
work page 2020
-
[21]
S. Ba˘ gci, M. Cin, H. Y. Uzunok, E. Karaca, H. M. T¨ ut¨ unc¨ u, and G. P. Srivastava, Investigating the normal state and superconducting state properties of orthorhom- bic and hexagonal ZrRuP: A first-principles study, Phys. Rev. B 100 184507 (2019)
work page 2019
-
[22]
L. Zhou, X. Chen, G. Shan, Q. Li, X. Kuang and X. Xing, Enhance suerconductivity in non-centrosymmetric hexagonal HfRuAs, Physica C: Superconductivity and its applications 633 1354710 (2025)
work page 2025
-
[23]
Migdal, Interactions between electrons and lattice vi- brations in a normal metal, Sov
A. Migdal, Interactions between electrons and lattice vi- brations in a normal metal, Sov. Phys. JETP 34 996 (1958)
work page 1958
-
[24]
Eliashberg, Interactions between electrons and lat tice vibrations in a superconductor, Sov
G. Eliashberg, Interactions between electrons and lat tice vibrations in a superconductor, Sov. Phys. JETP 11 696 (1960)
work page 1960
-
[25]
Eliashberg, Temperature Green’s function for elec- trons in a superconductor, Sov
G. Eliashberg, Temperature Green’s function for elec- trons in a superconductor, Sov. Phys. JETP 12 1000 (1961)
work page 1961
-
[26]
N. Marzari, A. A. Mostofi, J. R. Yates, I. Souza, and D. Vanderbilt, Maximally localized Wannier functions: The- ory and applications, Rev. Mod. Phys. 84, 1419 (2012)
work page 2012
-
[27]
F. Giustino, M. L. Cohen, and S. G. Louie, Electron- phonon interaction using Wannier functions, Phys. Rev. B 76, 165108 (2007)
work page 2007
-
[28]
Giustino, Electron-phonon interactions from first p rin- ciples, Rev
F. Giustino, Electron-phonon interactions from first p rin- ciples, Rev. Mod. Phys. 89, 015003 (2017)
work page 2017
-
[29]
E. R. Margine and F. Giustino, Anisotropic Migdal- Eliashberg theory using Wannier functions, Phys. Rev. B 87, 024505 (2013)
work page 2013
- [30]
-
[31]
P. B. Allen and B. Mitrovi´ c, Theory of Superconducting Tc, Solid State Phys. 37, 1 (1982)
work page 1982
-
[32]
P. Morel, and P. W. Anderson, Calculation of the Su- perconducting State Parameters with Retarded Electron- Phonon Interaction, Phys. Rev. 125, 125, 1263–1271 (1962)
work page 1962
-
[33]
H. Lee, S. Ponce, K. Bushick, S. Hajinazar, J. Lafuente Bartolome, J. Leveillee, C. Lian, J. M. Lihm, F. Macheda, and H. Mori, et al., Electron–phonon physics from first principles using the EPW code, npj, npj Comput Mater 9, 156 (2023)
work page 2023
-
[34]
P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- cioni, and I. Dabo, etal. , QUANTUM ESPRESSO: A modular and open-source software project for quantum simulations of materials, J. Phys.: Condens. Matter 21, 395502 (2009)
work page 2009
-
[35]
P. Giannozzi, O. Andreussi, T. Brumme, O. Bunau, M. B. Nardelli, M. Calandra, R. Car, C. Cavazzoni, D. Ceresoli, and M. Cococcioni, etal. , Advanced capabili- ties for materials with Quantum ESPRESSO, J. Phys.: Condens. Matter 29, 465901 (2017)
work page 2017
-
[36]
J. P. Perdew, Accurate Density Functional for the En- ergy: Real-Space Cutoff of the Gradient Expansion for the Exchange Hole, Phys. Rev. Lett. 55 1665–1668 (1985)
work page 1985
- [37]
-
[38]
D. R. Hamann, Optimized norm-conserving Vanderbilt pseudopotentials, Phys. Rev. B 88, 085117 (2013)
work page 2013
-
[39]
J. Noffsinger, F. Giustino, B. D. Malone, C.-H. Park, S. G. Louie, and M. L. Cohen, EPW: A program for cal- culating the electron-phonon coupling using maximally localized Wannier functions, Comput. Phys. Commun. 181, 2140 (2010)
work page 2010
-
[40]
P. V. S. Reddy, and V. Kanchana, Ab initio study of Fermi surface and dynamical properties of Ni 2XAl (X= Ti, V, Zr, Nb, Hf and Ta), J. Alloys and Comp., 616, 527-534 (2014)
work page 2014
-
[41]
P. V. S. Reddy, V. Kanchana, G. Vaitheeswaran, and D. J. Singh, Predicted superconductivity of Ni 2V Al and pressure dependence of superconductivity in Ni 2NbX (X = Al, Ga and Sn) and Ni 2V Al, J. Phys.: Conden. Matter, 28, 115703 (2016)
work page 2016
-
[42]
P. V. S. Reddy, V. Kanchana, T. E. Millichamp, G. Vaitheeswaran, and S. B. Dugdale, Enhanced supercon- ductivity in the high pressure phase of SnAs studied from first principles, Physica B, 505, 33-40 (2017)
work page 2017
-
[43]
P. V. S. Reddy and V. Kanchana, Enhanced supercon- ductivity in SnSb under pressure: a first principles study, J. Phys.: Conden. Matter, 29, 405502 (2017)
work page 2017
-
[44]
S. Sen, and G. -Y. Guo, Electronic structure, lattice dy - namics, and magnetic properties of Th XAsN (X=Fe, Co, Ni) superconductors: A first-principles study, Phys. Rev. B 102, 224505 (2020)
work page 2020
-
[45]
K. R. Babu, and G. -Y. Guo, Structure- and composition- tunable superconductivity, band topology, and elastic re- sponse of hard binary niobium nitrides Nb 2N, Nb 4N3 and Nb4N5, Phys. Rev. B 108, 064505 (2023)
work page 2023
-
[46]
S. P. Keshri, and G. -Y. Guo, Ab initio study of orbital- selective superconductivity in γ-BiPd, Phys. Rev. B 112, 214518 (2025)
work page 2025
-
[47]
C. Pellegrini, R. Heid, and A. Sanna, Eliashberg theory with ab-initio Coulomb interactions: A minimal numeri- cal scheme applied to layered superconductors, J. Phys. Mater. 5, 024007 (2022)
work page 2022
-
[48]
P. B. Allen, B Mitrovi´ c, Theory of Superconducting Tc, Solid State Physics 37, 1-92 (1983)
work page 1983
-
[49]
J. A. Flores-Livas, L. Boeri, A. Sanna, G. Profeta, R. Arita, M. Eremets, A perspective on conventional high- temperature superconductors at high pressure: Methods 12 and materials, Physics Reports 856 1–78, (2020)
work page 2020
-
[50]
A. Y. Liu and A. A. Quong, Linear-response calculation of electron-phonon coupling parameters, Phys. Rev. B 53 R7575 (1996)
work page 1996
- [51]
-
[52]
H. Salamati, F. S. Razavi, G. Quirion, Effect of pressure on the superconductivity of ZrRuP, Physica C 292 79–82 (1997)
work page 1997
-
[53]
P. B. Allen and R. C. Dynes, Transition temperatire of strong-coupled superconductors reanalyzed, Phys. Rev. B 12 905 (1975)
work page 1975
-
[54]
G. R. Stewart, Superconductivity in the A15 structure, Physica C 514 28–35 (2015)
work page 2015
-
[55]
D. V.Semenok, B. L.Altshuler, and E. A.Yuzbashyan, Fundamental Limits on the Electron-Phonon Coupling and Superconducting Tc, Adv. Mater.37, no. 40 (2025): 2507013
work page 2025
-
[56]
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, Nature 525, 73 (2015)
work page 2015
-
[57]
J. Bardeen, L. N. Cooper, and J. R. Schrieffer, Theory of superconductivity, Phys. Rev. 108, 1175 (1957)
work page 1957
-
[58]
D. C. Johnston, Supercond. Sci. Technol. 26, 115011 (2013)
work page 2013
-
[59]
J. -J. Zheng, and E. R. Margine, Electron-phonon cou- pling and pairing mechanism in β –Bi2Pd centrosymmet- ric superconductor, Phys. Rev. B 95, 014512 (2017)
work page 2017
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.