Three-dimensional spin susceptibility in Ba$_{0.75}$K$_{0.25}$Fe$_{2}$As$_{2}$: Out-of-plane modulation revealed by neutron spectroscopy and theoretical modeling

This study combines inelastic neutron scattering experiments and DFT-based theoretical modeling to reveal that the spin dynamics in Ba$_{0.75}$K$_{0.25}$Fe$_{2}$As$_{2}$ exhibit a distinct three-dimensional character with out-of-plane modulations at low energies that evolve into a two-dimensional profile at higher energies, a phenomenon accurately captured by a realistic 3D electronic band structure that highlights the crucial role of states away from the Fermi level over simple Fermi surface nesting.

The Gist

The Big Picture: A 3D Dance Floor

Imagine a superconductor (a material that conducts electricity with zero resistance) as a giant, bustling dance floor. In this specific material, Ba₀.₇₅K₀.₂₅Fe₂As₂, the dancers are electrons.

For a long time, scientists thought these electrons mostly danced in a flat, 2D pattern, like people moving only on the surface of a pond. They believed the "magnetic music" (spin fluctuations) that drives the superconductivity was also flat.

This paper proves that view is incomplete. The researchers discovered that at low energies (slow, gentle music), the electrons are actually dancing in 3D. They move up and down, creating a complex, layered structure. But as the music gets faster and more energetic (higher energy), the dancers forget the vertical moves and go back to dancing flatly on the surface.

The Experiment: The "Slow-Motion" Camera

To see this, the team used a massive machine called a neutron spectrometer. Think of this as a super-powered, slow-motion camera that shoots tiny, invisible bullets (neutrons) at the material.

• The Setup: They took a stack of crystal samples and aligned them perfectly, like a deck of cards.
• The Shot: They fired neutrons at different angles and energies.
• The Discovery: When they looked at the "low-energy" shots, they saw a strange pattern. The magnetic signal wasn't uniform; it was modulated. It was strong at certain "floors" (odd layers) and weak at others (even layers).

The Analogy: Imagine a multi-story building. If you shine a flashlight through the whole building, you might expect the light to pass through evenly. But here, the light was only bright on the 1st, 3rd, and 5th floors, and dim on the 2nd and 4th. This proved the magnetic activity had a distinct vertical (out-of-plane) structure.

However, when they turned up the energy (shone a brighter, faster flashlight), that "odd-floor" pattern disappeared. The light became uniform across all floors. The material had switched from a 3D dancer to a 2D dancer.

The Theory: The Digital Twin

To confirm what they saw, the scientists built a digital twin of the material using a supercomputer. They used a method called DFT (Density Functional Theory), which is like creating a hyper-realistic video game simulation of the electrons based on the laws of quantum physics.

• The Result: The computer simulation perfectly matched the real-world experiment. It showed the same "odd-floor" bright spots at low energy and the same fading away at high energy.
• The "Aha!" Moment: This proved that their computer model was accurate. It wasn't just a guess; it was a realistic map of the material's internal structure.

Why This Matters: Breaking the "Nesting" Myth

For years, scientists explained magnetism in these materials using a concept called "Fermi Surface Nesting."

• The Old Idea: Imagine two puzzle pieces (electron pockets and hole pockets). If you slide one over the other, they fit perfectly. This "perfect fit" was thought to be the only reason the magnetic waves happened.
• The New Finding: The researchers found that this "puzzle piece" idea doesn't work for the vertical (3D) part of the magnetism. The puzzle pieces don't actually fit perfectly in 3D space.

Instead, they found that the "magnetic peak" (the strongest dance move) is driven by electrons that are not even on the main dance floor (the Fermi level). It's like the rhythm of the party is being set by the people in the VIP lounge upstairs, not just the people on the main floor. This means we can't just look at the surface to understand the material; we have to look at the whole 3D building.

The Takeaway

1. It's 3D, not 2D: At low energies, the magnetism in this superconductor is a full 3D phenomenon, not just a flat sheet.
2. The Crossover: As energy increases, the material naturally shifts from 3D behavior to 2D behavior.
3. The Model Works: The computer models based on real physics (DFT) can predict this complex 3D behavior without needing to "fudge" the numbers.
4. The Lesson: To understand how these superconductors work, we must stop treating them like flat pancakes and start treating them like layered cakes. The "secret sauce" for superconductivity might be hiding in those vertical layers.

In short: The scientists used a neutron camera and a supercomputer to prove that the magnetic "heartbeat" of this superconductor has a vertical rhythm that changes depending on how fast the energy is moving. This helps us build better models to design future superconductors.

Technical Summary

1. Problem Statement

Iron-based superconductors (FeSCs), particularly the 122 family (e.g., BaFe$_2$As$_2$), exhibit superconductivity in close proximity to antiferromagnetic (AFM) order. While the pairing mechanism is widely believed to be mediated by spin fluctuations, most theoretical models and experimental analyses have historically treated these systems as quasi-two-dimensional (2D).

• The Gap: Although parent compounds like BaFe$_2$As$_2$ exhibit 3D long-range AFM order with a wavevector $q_{AFM} = (0.5, 0.5, 1)$, theoretical studies often assume a 2D electronic structure, neglecting the momentum dependence along the out-of-plane ($L$) direction.
• The Question: Can Density Functional Theory (DFT)-derived models, which inherently capture 3D electronic structures, accurately reproduce the experimentally observed out-of-plane modulation of spin susceptibility? Furthermore, how does this 3D character evolve with energy, and is the "Fermi surface nesting" picture sufficient to explain the 3D magnetic instability?

2. Methodology

The study employs a combined experimental and theoretical approach to investigate the dynamical spin susceptibility of hole-doped Ba$_{0.75}$K$_{0.25}$Fe$_2$As$_2$.

Experimental Approach

• Sample: High-quality single crystals of Ba$_{0.75}$K$_{0.25}$Fe$_2$As$_2$ grown via the FeAs-flux method, co-aligned to a total mass of 5.0 g.
• Technique: Time-of-flight (TOF) inelastic neutron scattering (INS) using the AMATERAS and 4SEASONS spectrometers at J-PARC.
• Data Acquisition: Unlike previous studies that used fixed-geometry scans (assuming $L$-independence), the authors performed multi-orientation scans. By rotating the crystal array over a wide angular range ($\phi \in [-40^\circ, +60^\circ]$) and utilizing continuous rotation techniques, they reconstructed the full 4D scattering function $S(Q, \omega)$ (momentum $Q$ and energy $\omega$).
• Energy Range: Measurements were conducted with incident neutron energies ($E_i$) of 7.74, 15.15, 42.0, 55.6, and 125 meV to cover both low-energy (near $q_{AFM}$) and high-energy regimes.

Theoretical Approach

• Model: A realistic 3D electronic band structure derived from DFT (using Quantum ESPRESSO and Wannier90).
  • The unit cell was unfolded to a 1-Fe basis to match the symmetry of the Fe sublattice.
  • A 5-orbital tight-binding model was constructed.
  • K-doping was modeled via a rigid band shift, and band narrowing was accounted for by rescaling the DFT bands by a factor of 3 to match ARPES data.
• Calculation: The dynamical spin susceptibility $\hat{\chi}_s(q, \omega)$ was calculated using the Random Phase Approximation (RPA).
  • The calculation included the full 3D momentum dependence ($k_z$ dispersion).
  • Interaction parameters (Hubbard $U$, Hund's coupling $J$, etc.) were set to standard values for FeSCs.

3. Key Results

Experimental Findings: 3D-to-2D Crossover

• Low-Energy Regime ($\omega < 15$ meV): The INS data revealed a pronounced out-of-plane modulation of magnetic intensity along the $(0.5, 0.5, L)$ direction.
  • Intensity peaks at odd $L$ positions (e.g., $L=1$) and is suppressed at even $L$.
  • This confirms the 3D nature of low-energy spin fluctuations and the tendency toward AFM order at $q_{AFM} = (0.5, 0.5, 1)$.
• High-Energy Regime ($\omega > 15$ meV): As energy increases, the odd-$L$ modulation gradually fades.
  • At $\omega \approx 20$ meV and above (up to 125 meV), the intensity distribution becomes nearly uniform along $L$, indicating a crossover to a 2D-like magnetic response.
• Temperature Dependence: The out-of-plane modulation persists well above the Néel temperature ($T_N = 90$ K), confirming that these are robust paramagnetic spin correlations rather than just static order.

Theoretical Findings: Validation of DFT-RPA

• Reproduction of Modulation: The DFT-derived RPA model successfully reproduced the experimental features:
  • A strong peak at $L=1$ at low energies ($\omega = 5$ meV).
  • A gradual suppression of the $L$-modulation at higher energies ($\omega = 50$ meV), matching the experimental 3D-to-2D crossover.
• Origin of the Peak:
  • Not just Nesting: Analysis showed that the peak at $L=1$ cannot be explained solely by Fermi surface nesting (geometric overlap of hole and electron pockets). Nesting alone produced multiple peaks and failed to select $L=1$ uniquely.
  • Role of Non-Fermi Surface States: When the calculation was restricted to electronic states within $\pm 5$ meV of the Fermi level, the peak structure was incorrect. Including states far from the Fermi level was crucial for correctly selecting the $L=1$ instability. This highlights the limitation of simple nesting pictures.

4. Key Contributions

1. Experimental Benchmark: Provided the first comprehensive 4D mapping of spin susceptibility in FeSCs, explicitly demonstrating the 3D-to-2D crossover in spin dynamics as a function of energy.
2. Theoretical Validation: Demonstrated that DFT-derived itinerant models (without phenomenological parameter tuning) can accurately capture complex 3D magnetic correlations, validating their use for predicting material-specific electronic structures.
3. Mechanism Insight: Established that electronic states away from the Fermi level play a decisive role in shaping the out-of-plane AFM instability, challenging the conventional view that spin fluctuations in FeSCs are governed purely by Fermi surface nesting.
4. Methodological Advance: Showed that multi-orientation TOF neutron scattering is essential for resolving out-of-plane dynamics, which are often missed in fixed-geometry 2D scans.

5. Significance

• For Superconductivity Theory: The study confirms that the pairing mechanism in 122 FeSCs must account for 3D electronic structures. Since the superconducting gap symmetry (e.g., horizontal nodes) is sensitive to $k_z$ dispersion, accurate 3D models are critical for understanding the pairing glue.
• For Model Development: It serves as a rigorous benchmark for theoretical methods. The success of the DFT-RPA approach suggests it is a reliable tool for describing FeSCs, though the authors note that more strongly correlated systems (like FeSe) may require advanced methods like DMFT.
• General Physics: The finding that non-Fermi-surface states drive magnetic instabilities suggests that a full-band calculation is necessary for understanding magnetic ordering in correlated electron systems, moving beyond simplified low-energy effective models.

In conclusion, the paper bridges the gap between experimental observation and theoretical modeling, proving that the 3D nature of spin fluctuations is a fundamental, energy-dependent feature of iron-based superconductors that cannot be ignored.