Strong-coupling anisotropic superconductivity in hexagonal HfRuAs from anisotropic Migdal-Eliashberg theory

This study employs anisotropic Migdal-Eliashberg theory combined with *ab initio* calculations to demonstrate that hexagonal HfRuAs is a strong-coupling, phonon-mediated superconductor characterized by a single anisotropic $s$-wave gap and a critical temperature consistent with experimental observations.

The Gist

Imagine a world where electricity flows without any resistance at all. This is the magic of superconductivity. For a long time, scientists have been trying to understand exactly how certain materials pull off this trick.

This paper is a deep dive into a specific material called hexagonal HfRuAs (a crystal made of Hafnium, Ruthenium, and Arsenic). The researchers used powerful computer simulations to figure out why this material becomes a superconductor and how it behaves.

Here is the story of their findings, broken down into simple concepts:

1. The "Dance Floor" and the "Music"

In this material, the electrons are like dancers on a crowded floor. Normally, they bump into each other and lose energy (resistance). But when the material gets cold enough, they start pairing up to dance perfectly in sync.

• The Music (Phonons): The paper explains that the "music" that gets these electrons dancing is actually the vibration of the atoms themselves. Think of the atoms as people jumping on a trampoline. When they jump, they create waves.
• The Strong Connection: The researchers found that the connection between the dancing electrons and the jumping atoms is incredibly strong. It's not a gentle tap; it's a firm handshake. In scientific terms, they call this "strong-coupling." The strength of this connection is measured at about 1.56, which is much higher than what you see in standard superconductors.

2. The "Heavy" and "Light" Dancers

The material has different "sheets" or layers of electrons (called Fermi surfaces). The paper discovered that the music isn't played equally everywhere:

• The Low Notes: The most important vibrations are the slow, low-frequency ones. These are mostly caused by the heavy Hafnium and Ruthenium atoms shaking.
• The Anisotropy (The Lopsided Dance): The dance isn't the same in every direction. On some parts of the electron "floor," the connection to the music is very strong, while on others, it's weaker. It's like a dance floor where the music is loud and clear in the center but gets muffled at the edges. This unevenness is called anisotropy.

3. The "Gap" in the Energy

To become a superconductor, the electrons need to open a "gap" in their energy levels—a protective barrier that keeps them from getting disrupted.

• A Single, Wobbly Shield: The paper found that this material has one main shield (a single gap), not multiple different ones. However, because of the "lopsided" dance mentioned earlier, this shield isn't a perfect, uniform circle. It's more like a slightly squashed or wobbly circle.
• No Holes: Crucially, the shield is fully closed. There are no holes or gaps in the shield itself. This means the superconductivity is very stable and follows a classic "s-wave" pattern (a standard, safe type of superconductivity).

4. The Temperature Puzzle

The researchers calculated that this material should become a superconductor at a temperature around 16 Kelvin (very cold, but not that cold).

• The Discrepancy: Real-world experiments have shown this material becoming superconductive at lower temperatures (between 4 K and 7 K).
• Why the difference? The paper suggests that the computer model represents a "perfect" crystal with no flaws. Real-world samples might have tiny impurities, defects, or mixed phases that act like "speed bumps," slowing down the superconductivity and lowering the temperature at which it happens.

5. The Big Conclusion

The main takeaway is that hexagonal HfRuAs is a "strong-coupling" superconductor.

• Analogy: If a weak-coupling superconductor is like two people holding hands lightly while walking, a strong-coupling superconductor is like two people locked in a tight embrace, moving as one unit.
• The Evidence: The ratio of the energy gap to the temperature is much higher than the standard limit for weak superconductors, proving that the "embrace" between the electrons and the vibrating atoms is very tight.

In summary: The paper uses advanced math to show that HfRuAs is a robust superconductor driven by strong vibrations of its own atoms. While the real-world samples aren't quite as perfect as the computer model predicts, the fundamental physics reveals a material where electrons and atoms dance together with surprising intensity.

Technical Summary

Technical Summary: Strong-coupling anisotropic superconductivity in hexagonal HfRuAs from anisotropic Migdal-Eliashberg theory

Problem and Motivation
Ternary equiatomic intermetallic compounds of the form $TT'X$ (where $T, T'$ are transition metals and $X$ is a pnictogen or silicon) often exhibit superconductivity, particularly in the hexagonal $ZrNiAl$-type ($h$-phase) structure. While members of this family like $HfRuP$ and $ZrRu(As, Si, P)$ display transition temperatures ($T_c$) exceeding 10 K, recent experimental reports on hexagonal $HfRuAs$($h$-$HfRuAs$) show a wide variation in $T_c$ (ranging from 4.3 K to 7.9 K depending on synthesis conditions). Previous theoretical studies on related compounds have similarly yielded disparate $T_c$ predictions. Furthermore, while experimental specific heat data suggests a weak-coupling, fully gapped $s$-wave state with a reduced specific heat jump ($\Delta C/\gamma T_c \approx 1.03$), a comprehensive theoretical understanding of the momentum-dependent electron-phonon coupling (EPC) and the anisotropic superconducting (SC) gap structure in $h$-$HfRuAs$ remains lacking. Specifically, fully anisotropic Migdal–Eliashberg (ME) calculations are required to resolve the multiband SC properties and Fermi surface (FS) dependent pairing interactions that standard isotropic approximations may miss.

Methodology
The authors performed a comprehensive first-principles investigation of $h$-$HfRuAs$ using the following workflow:

1. Electronic and Lattice Properties: Electronic structure, phonon dispersion, and electron-phonon coupling matrix elements were calculated using Density Functional Theory (DFT) and Density Functional Perturbation Theory (DFPT) within the Quantum ESPRESSO package. The PBEsol exchange-correlation functional and optimized norm-conserving Vanderbilt pseudopotentials were employed.
2. Wannier Interpolation: To achieve the dense sampling required for accurate EPC calculations, electron-phonon matrix elements computed on coarse Brillouin zone meshes were interpolated onto fine $k$- and $q$-point grids (32 $\times$ 32 $\times$ 48) using maximally localized Wannier functions via the EPW code.
3. Anisotropic Migdal-Eliashberg Formalism: The SC properties were investigated by solving the fully anisotropic ME equations on the imaginary-frequency axis. The calculations utilized the anisotropic EPC kernel and a screened Coulomb pseudopotential ($\mu^*_c = 0.15$), chosen based on the dominance of $Hf$-$d$ and $Ru$-$d$ orbitals at the Fermi level. The SC gap function on the real-frequency axis was obtained via analytic continuation to calculate the quasiparticle density of states (DOS).

Key Results

• Strong Electron-Phonon Coupling: The calculated Eliashberg spectral function ($\alpha^2F(\omega)$) reveals a total EPC constant $\lambda \approx 1.56$. This value places $h$-$HfRuAs$ in the strong-coupling regime, significantly higher than related compounds like $h$-$ZrRuAs$($\lambda \approx 0.79\text{--}1.32$) and $h$-$HfRuP$($\lambda \approx 0.87$). The coupling is dominated by low-frequency phonon modes ($< 10$ meV) associated primarily with $Hf$ and $Ru$ vibrations, which contribute approximately 82% of the total EPC.
• Anisotropic Multiband Superconductivity: The Fermi surface consists of three sheets: two hole-like bands ($H1, H2$) and one electron-like band ($E1$). The momentum-resolved EPC ($\lambda_{nk}$) and the SC gap ($\Delta_{nk}$) exhibit pronounced anisotropy across these sheets, with the largest variations occurring on the hole-like $H2$ band.
• Gap Structure and Symmetry: Despite the strong anisotropy, the SC state is characterized by a single anisotropic gap with overall $s$-wave symmetry. The gap is fully open (no nodes) across all Fermi surface sheets, evidenced by a U-shaped low-energy quasiparticle DOS. The gap magnitude centers around $\Delta \approx 2.9$ meV with a spread of $\sim 0.8$ meV.
• Transition Temperature and Gap Ratio: The calculated $T_c$ is approximately 16.06 K (using the fitted BCS-type temperature dependence of the averaged gap), which is higher than most experimental reports but consistent in order of magnitude with the upper bounds of previous theoretical predictions. The gap ratio $2\Delta(0)/k_B T_c \approx 4.25$ significantly exceeds the weak-coupling BCS limit of 3.53, further confirming the strong-coupling nature of the superconductivity.
• Coulomb Pseudopotential Sensitivity: The study systematically varied $\mu^*_c$ from 0.11 to 0.30. While increasing $\mu^*_c$ suppresses $T_c$ (from 13.9 K to 6.5 K) and reduces the gap magnitude, the momentum-space distribution and the anisotropic character of the gap remain robust, indicating that the qualitative conclusions regarding strong coupling and anisotropy are not dependent on the specific choice of $\mu^*_c$.

Significance and Claims
The paper establishes $h$-$HfRuAs$ as a phonon-mediated, strongly coupled anisotropic superconductor. The primary contribution is the demonstration that the superconductivity in this material is driven by strong coupling between electrons and low-frequency lattice vibrations of $Hf$ and $Ru$, leading to a large EPC constant ($\lambda \approx 1.56$). The authors clarify that while the system exhibits significant multiband anisotropy in both the EPC and the SC gap, it retains a fully gapped, single-gap $s$-wave symmetry. The discrepancy between the calculated $T_c$ (16 K) and experimental values (4.3–7.25 K) is attributed to extrinsic factors in experimental samples (disorder, phase coexistence, stoichiometry deviations) rather than a failure of the intrinsic theoretical model. The work provides a detailed microscopic description of the role of momentum-dependent electron-phonon interactions in determining the SC properties of this family of compounds, offering a consistent theoretical baseline for future investigations on high-quality single crystals.