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

This study utilizes Nuclear-Electronic Orbital DFT to demonstrate that while nuclear quantum effects significantly stiffen S-H bonds and alter phonon properties in high-pressure H$_3$S, their direct impact on the electronic structure is minimal, indicating that the observed isotope effect on the critical temperature arises primarily from phonon changes rather than electronic modifications.

The Gist

The Big Picture: A Superconductor's Secret Life

Imagine you have a material called H3S (Hydrogen Sulfide). When you squeeze it incredibly hard (like 2 million times the pressure of the atmosphere), it becomes a superconductor. This means electricity can flow through it with zero resistance, but only if it's kept very cold.

For a long time, scientists have known that H3S is a champion superconductor, working at temperatures around -70°C (which is actually "hot" for superconductors!). They knew how it worked: tiny vibrations in the material (called phonons) help electrons pair up and dance together without friction.

But there was a nagging question: Does the fact that the hydrogen atoms are "quantum" matter?

In the quantum world, particles aren't like billiard balls; they are more like fuzzy clouds of probability. They vibrate even when they are supposed to be still (Zero-Point Energy). Because hydrogen is the lightest atom, it is very "fuzzy" and jittery. Scientists wanted to know: Does this quantum fuzziness change the electronic structure (the dance floor for the electrons) enough to explain why the superconducting temperature drops when you swap Hydrogen for its heavier cousin, Deuterium?

The New Tool: The "Two-Player" Simulator

To answer this, the researchers used a new, fancy computer simulation method called NEO-DFT.

• The Old Way (Standard DFT): Imagine a dance floor where the dancers (electrons) are treated as quantum waves, but the furniture (the atomic nuclei) is treated as solid, heavy, classical chairs. The chairs don't move or wiggle; they just sit there.
• The New Way (NEO-DFT): This method treats the hydrogen atoms (the protons) just like the electrons. They are also fuzzy, wiggly quantum clouds. The simulation solves for the electrons and the protons at the same time, as if they are two players in a game who are constantly reacting to each other.

What They Found: Two Different Stories

The researchers ran the simulation and found two very different stories happening at the same time.

1. The Electronic Story: "A Tiny Tweak"

The Analogy: Imagine the electrons are a crowd of people dancing on a floor. The "Fermi energy" is the VIP section of the dance floor. Scientists thought that if the hydrogen atoms (the furniture) started wiggling quantum-mechanically, it might completely reshape the dance floor, changing the music and the number of people who can dance.

The Result: The wiggle didn't change the dance floor much.

• The "fuzziness" of the protons did shift the energy levels slightly, like moving a few chairs an inch to the left.
• It changed the number of available dance spots (Density of States) by only about 5% to 7%.
• The Bottom Line: If this were the only thing happening, the superconducting temperature would go up by a tiny bit (about 5 degrees Kelvin). This is a very small change.

2. The Vibration Story: "The Stiffening Spring"

The Analogy: Now, imagine the hydrogen atoms are connected to the sulfur atoms by springs. In the old "classical" view, these springs are soft. But in the "quantum" view, because the hydrogen is so jittery and energetic, it pulls the springs tight. The springs become stiffer.

The Result: The vibrations (phonons) changed drastically.

• The high-pitched vibrations of the hydrogen atoms became much faster and higher energy because the "springs" stiffened due to quantum effects.
• This stiffening is the real reason why the superconducting temperature changes when you swap Hydrogen for Deuterium. Deuterium is heavier, so it wiggles less, the springs don't stiffen as much, and the superconducting temperature drops.

The Grand Conclusion

The paper solves a mystery about why H3S behaves the way it does.

• The Misconception: Many scientists thought the drop in temperature when using Deuterium was because the electronic structure (the dance floor) changed significantly.
• The Reality: The electronic structure barely changed. The "dance floor" stayed mostly the same.
• The Real Culprit: The change comes from the vibrations. The quantum nature of the protons makes the chemical bonds "stiffer," which changes how the material vibrates. Since superconductivity relies on these vibrations, the change in vibration is what actually controls the temperature.

Why This Matters

This study is like realizing that to understand why a car engine is loud, you don't need to look at the paint job (the electronic structure); you need to look at the pistons (the vibrations).

By using this new "Two-Player" simulation (NEO-DFT), the authors proved that we don't need to worry about the electrons changing their minds when protons get quantum; we just need to focus on how the protons stiffen the bonds. This gives scientists a clearer, more unified way to predict and design future superconductors that might work at room temperature.

In short: The quantum jitter of hydrogen doesn't change the map of the electrons, but it definitely changes the music (vibrations) the electrons dance to. And that music is what makes the superconductor work.

Technical Summary

1. Problem Statement

High-pressure sulfur hydride (H$_3$S) is a benchmark hydrogen-rich superconductor with a critical temperature ($T_c$) exceeding 200 K. While first-principles theories (based on the Born-Oppenheimer approximation) have successfully predicted its existence and properties, a key question remains regarding the role of Nuclear Quantum Effects (NQEs).

• The Gap: Previous studies treated protons classically or applied perturbative zero-point renormalization to the electronic structure. It is unclear to what extent the explicit quantum mechanical nature of protons (delocalization, zero-point energy) alters the electronic band structure and density of states (DOS) near the Fermi level.
• The Specific Question: Do NQEs significantly modify the electronic features (specifically van Hove singularities) that drive the high $T_c$ in H$_3$S, or is the observed isotope effect (reduction in $T_c$ upon deuteration) driven primarily by changes in phonon properties?

2. Methodology

The authors employed Nuclear-Electronic Orbital Density Functional Theory (NEO-DFT), a multicomponent approach that transcends the Born-Oppenheimer approximation.

• Framework: Instead of treating nuclei as classical point charges, selected nuclei (hydrogen/protons) are treated as quantum particles on the same footing as electrons.
• Formalism: The method solves coupled Kohn-Sham equations for both electrons and protons:
  • Electrons: Standard KS equations with effective potentials including electron-proton interactions.
  • Protons: KS-like equations where the proton wavefunction is solved self-consistently. The proton mass ($M_p$) is explicitly included in the kinetic energy term.
• Correlation: The study utilized an electron-proton correlation functional (epc17-2, LDA-type) to account for non-classical electron-proton interactions beyond simple electrostatics.
• Computational Details:
  • System: H$_3$S and D$_3$S in the $Im\bar{3}m$ phase at 200 GPa (lattice constant 2.985 Å).
  • Software: Periodic NEO-DFT implementation within the all-electron FHI-aims code.
  • Basis Sets: Numeric atom-centered orbitals (NAO) for electrons; Gaussian-type basis sets for protons.
  • Phonons: Calculated using the finite-displacement method with a 4$\times$4$\times$4 supercell, treating displaced hydrogen atoms as quantum protons within a constrained NEO-DFT framework.

3. Key Contributions

1. First Self-Consistent Multicomponent Treatment: This is the first application of a fully self-consistent NEO-DFT method to the electronic structure of high-pressure H$_3$S, explicitly including quantum protons in the ground-state energy functional without perturbative corrections.
2. Decoupling Electronic vs. Phononic NQEs: The study rigorously separates the impact of NQEs on the electronic structure (band structure, DOS) from their impact on lattice dynamics (phonons).
3. Validation of Deuteration Mechanism: The work provides a theoretical basis for why deuteration lowers $T_c$, distinguishing between electronic and vibrational contributions.

4. Key Results

A. Quantum Proton Properties

• Localization: Despite the light mass of hydrogen, the protons in H$_3$S remain highly localized between sulfur atoms due to the steep, symmetric single-well potential of the $Im\bar{3}m$ phase. This preserves the high symmetry of the lattice.
• Delocalization: Including electron-proton correlation increases the spatial spread of the proton density (from ~0.009 Å$^2$ to ~0.020 Å$^2$) but has a negligible effect on the overall electronic band structure.
• Zero-Point Energy (ZPE): The calculated ZPE per proton is ~0.389 eV (with correlation), confirming significant quantum effects.

B. Electronic Structure (Band Structure & DOS)

• Subtle Modifications: NQEs induce only minor changes to the electronic band structure near the Fermi energy ($E_F$).
  • The van Hove singularities (vHS) near $E_F$ are shifted slightly upward in energy.
  • The Fermi surface topology remains largely unchanged, though saddle points shift slightly.
• Density of States (DOS): The DOS at the Fermi level, $N(\epsilon_F)$, increases by approximately 6.9% for H$_3$S and 5.0% for D$_3$S due to NQEs.
• Impact on $T_c$ (Electronic): Based on the Allen-Dynes strong-coupling limit ($T_c \propto \sqrt{N(\epsilon_F)}$), this increase in DOS would theoretically raise $T_c$ by only ~2.5% to 3.4% (approx. 5 K). This is considered negligible compared to other theoretical uncertainties.

C. Phonon Dispersion

• Significant Hardening: In contrast to the electronic results, NQEs cause large changes in the phonon dispersion, specifically for the high-frequency hydrogen-dominated modes.
• Mechanism: The zero-point motion of protons effectively "stiffens" the S-H bonds, shifting phonon frequencies to higher energies (hardening). Conversely, S-H bond-bending modes are softened.
• Isotope Effect: The changes are more pronounced for Hydrogen than Deuterium, consistent with the mass difference.
• Impact on $T_c$ (Phononic): Previous studies (cited in the paper) indicate that this phonon renormalization reduces the electron-phonon coupling constant ($\lambda$) by ~30%, leading to a 20–25% reduction in $T_c$.

5. Significance and Conclusion

• Resolution of the Isotope Effect: The study concludes that the experimentally observed reduction in $T_c$ upon deuteration (from H$_3$S to D$_3$S) arises predominantly from changes in phonon properties (specifically the stiffening of S-H bonds and reduction in electron-phonon coupling), rather than from modifications to the electronic structure itself.
• Electronic Structure Stability: The electronic band structure and DOS of H$_3$S are remarkably robust against nuclear quantum effects. The "quantum smearing" of protons does not significantly destroy the van Hove singularities responsible for the high $T_c$.
• Methodological Advancement: The paper demonstrates that NEO-DFT provides a unified, first-principles framework capable of simultaneously capturing both electronic and phonon quantum effects. This validates the use of multicomponent DFT for predicting superconductivity in hydrogen-rich materials, confirming that accurate $T_c$ prediction requires explicit quantum treatment of lattice dynamics, even if the electronic structure can be approximated classically.

In summary, while proton quantum effects are crucial for accurately modeling the phonon spectrum and predicting the correct critical temperature, their direct influence on the electronic density of states in high-pressure H$_3$S is minimal.