Electron-phonon coupling in magnetic materials using the local spin density approximation

This paper extends the EPW package to calculate electron-phonon coupling in magnetic materials using the local spin density approximation, revealing that while electron-phonon scattering dominates resistivity in ferromagnetic iron, it plays a minor role in nickel, and confirming that phonon-driven superconductivity is intrinsically suppressed in both systems.

The Gist

The Big Picture: A City in Motion

Imagine a metal like Iron (Fe) or Nickel (Ni) as a busy, bustling city.

• The Electrons are the commuters trying to get from point A to point B (conducting electricity).
• The Atoms in the metal are the buildings that make up the city streets.
• The Phonons are the vibrations of those buildings (like an earthquake or people dancing on the floor).
• The Magnetism is the mood of the city. In these metals, the commuters have a specific "spin" (like wearing a red shirt or a blue shirt), and they all agree to wear the same color, creating a magnetic field.

The goal of this paper is to understand how traffic jams (resistance) happen in this city when the buildings start shaking (vibrating).

---

1. The Problem: The Old Map Was Wrong

For a long time, scientists had a map (a computer code called EPW) to predict how these commuters move. But this map had a blind spot: it assumed the city was neutral (no magnetic mood). It treated everyone the same.

However, in magnetic materials like Iron and Nickel, the "mood" (magnetism) changes everything.

• The Analogy: Imagine trying to navigate a city where the traffic lights only work if you know which side of the street you are on. If you ignore the "magnetic mood," your map predicts the buildings are stable. But in reality, if you ignore the mood, the buildings might start to collapse (become unstable).

The authors fixed the map. They upgraded the EPW software to understand magnetism. Now, it can tell the difference between "Red Shirt" commuters and "Blue Shirt" commuters and how they interact with the shaking buildings.

2. The Test Drive: Iron vs. Nickel

The team tested their new map on two famous cities: Iron (Fe) and Nickel (Ni). They wanted to see how much the shaking buildings slowed down the traffic.

City Iron (Fe): The "Shaking" City

• What happened: In Iron, the "Red Shirt" commuters (the majority) are the main drivers of traffic.
• The Finding: The traffic jams in Iron are caused almost entirely by the buildings shaking (electron-phonon scattering).
• The Analogy: Imagine a city where the only reason you are late is because the sidewalks are vibrating. If you stop the shaking, the traffic flows perfectly.
• The Danger: If you tried to model Iron without magnetism (ignoring the mood), the map would predict the buildings are shaking so violently they would fall apart (unstable phonons). This proves you must include magnetism to get the physics right.

City Nickel (Ni): The "Chaos" City

• What happened: Nickel is different. It has a mix of "Red" and "Blue" commuters, but the "Blue" ones are much more numerous near the exit.
• The Finding: Surprisingly, the shaking buildings are not the main cause of traffic jams in Nickel. In fact, the shaking buildings only account for less than one-third of the delay.
• The Mystery: What causes the rest of the delay? The paper suggests it's magnetic fluctuations (like the commuters arguing or changing their minds).
• The Analogy: In Nickel, the sidewalks are actually quite stable. The traffic jam happens because the commuters are constantly switching teams or getting distracted by the magnetic mood. If you only look at the shaking sidewalks (the old way), you would think the traffic is fine, but in reality, it's a mess.

3. The "Superconductivity" Dream

Scientists often hope that if they can make the buildings vibrate in a specific way, the commuters can pair up and glide through the city without any friction at all. This is called Superconductivity.

• The Result: The authors checked if Iron or Nickel could do this.
• The Verdict: No. Even with their new, perfect map, they found that the magnetic mood of these cities is too strong. It acts like a wall that prevents the commuters from pairing up.
• The Metaphor: It's like trying to get two people to dance perfectly in sync while a loud, chaotic band is playing right next to them. The magnetic noise is just too loud; the dancers can't hear each other. So, Iron and Nickel will never be superconductors under normal conditions.

4. Why This Matters

Why should you care about traffic jams in a metal city?

1. Better Electronics: We are running out of space for traditional computer chips (Moore's Law is ending). The future is Spintronics, where we use the "spin" (the shirt color) to carry information. To build these devices, we need to know exactly how the "spin" interacts with the "shaking buildings."
2. Energy Efficiency: Every time electricity hits a traffic jam, it turns into heat (energy loss). By understanding that Iron and Nickel lose energy for different reasons, engineers can design better materials to carry electricity with less waste.
3. The "Don't Ignore the Mood" Lesson: The biggest takeaway is that you cannot ignore magnetism. If you try to calculate how Iron works without considering its magnetic nature, you get a result that is not just slightly wrong, but completely broken (predicting the metal falls apart).

Summary

The authors built a new, high-tech GPS (the updated EPW code) that understands magnetic cities. They drove it through Iron and Nickel and discovered:

• Iron is slowed down mostly by vibrations.
• Nickel is slowed down mostly by magnetic chaos, not vibrations.
• Neither city can become a friction-free super-highway (superconductor) because the magnetic noise is too loud.

This helps scientists design the next generation of energy-efficient computers and magnetic storage devices.

Technical Summary

1. Problem Statement

Electron-phonon (e-ph) interactions are fundamental to understanding transport properties (resistivity, mobility) and superconductivity in materials. While first-principles methods like Density Functional Perturbation Theory (DFPT) combined with Maximally Localized Wannier Functions (MLWF) have successfully modeled these interactions in non-magnetic systems, a significant gap exists for magnetic materials.

• The Challenge: Standard implementations of e-ph interpolation (e.g., in the EPW code) assume time-reversal symmetry, which is broken in magnetic systems. Neglecting spin degrees of freedom in magnetic metals (like Iron and Nickel) leads to:
  • Unstable phonon modes (imaginary frequencies) due to incorrect lattice dynamics.
  • Inaccurate predictions of carrier resistivity.
  • An inability to distinguish between scattering mechanisms (e.g., electron-phonon vs. electron-magnon).
• The Need: A robust, first-principles framework capable of handling collinear magnetism to accurately interpolate electron-phonon matrix elements (EPME) and compute transport properties without the prohibitive cost of non-collinear (spin-orbit coupling) calculations.

2. Methodology

The authors extended the EPW (Electron-Phonon Wannier) code to support collinear magnetism, integrating it with the Abinit and Quantum ESPRESSO (QE) DFT codes.

• Theoretical Framework:
  • Spin-Resolved EPME: The electron-phonon matrix element $g_{mn\nu}^{\sigma}(k, q)$ was redefined to include spin indices ($\sigma, \sigma'$). In the collinear approximation, spin-up ($\uparrow$) and spin-down ($\downarrow$) channels are treated separately but self-consistently.
  • Deformation Potential: Defined as $D_{mn\nu}^{\sigma}(k, q)$, measuring the change in the self-consistent potential due to lattice vibrations for specific spin channels.
  • Transport Theory: The authors solved the linearized Boltzmann Transport Equation (BTE) for each spin channel. The total resistivity is calculated using the parallel resistor model: $\rho^{-1} = (\rho^{\uparrow})^{-1} + (\rho^{\downarrow})^{-1}$, contrasting with older additive approximations (LOVA/Ziman) that were shown to be less accurate for anisotropic magnetic systems.
• Computational Setup:
  • Materials: Ferromagnetic Body-Centered Cubic (BCC) Iron (Fe) and Face-Centered Cubic (FCC) Nickel (Ni).
  • Functionals: Local Density Approximation (LDA) and Generalized Gradient Approximation (PBE).
  • Validation: The implementation was validated by comparing direct DFPT calculations in QE and Abinit, and by cross-verifying transport results with a custom ElectronPhonon.jl package.
  • Grids: Used dense momentum grids (up to $112^3$) via Wannier interpolation to ensure convergence of transport properties.

3. Key Contributions

1. Code Development: First implementation of collinear magnetic support in the EPW package, enabling the interpolation of spin-resolved electron-phonon matrix elements.
2. Validation: Rigorous verification against Abinit and QE, demonstrating agreement in total energies, magnetization, phonon frequencies, and deformation potentials (errors < 0.02 meV/bohr).
3. Physical Insight:
   • Demonstrated that neglecting magnetism in Fe leads to dynamically unstable phonons (imaginary modes), whereas Ni remains stable but yields incorrect transport physics.
   • Established that the density of states (DOS) at the Fermi level and Fermi surface nesting are the primary drivers for spin-resolved coupling strength differences, rather than just the deformation potential.
4. Transport Mechanism Differentiation: Quantified the distinct roles of electron-phonon scattering versus magnon scattering in Fe and Ni, resolving long-standing discrepancies in resistivity modeling.

4. Key Results

A. Electron-Phonon Coupling Strength ($\lambda$)

• Fe: The coupling is dominated by the spin-up (majority) channel ($\lambda_{\uparrow} \approx 0.22$, $\lambda_{\downarrow} \approx 0.08$). This is driven by a high DOS in the majority channel at the Fermi level.
• Ni: The coupling is dominated by the spin-down (minority) channel ($\lambda_{\downarrow} \approx 0.27$, $\lambda_{\uparrow} \approx 0.004$). This is due to the high DOS of minority states near the Fermi level (a result of d-band filling).
• Superconductivity: Calculations of the phonon-mediated critical temperature ($T_c$) using the Allen-Dynes formula yielded vanishingly small values for both materials in their ferromagnetic state, confirming that neither is a phonon-mediated superconductor under ambient conditions.

B. Electrical Resistivity

• Iron (Fe):
  • Electron-phonon scattering is the dominant mechanism up to room temperature.
  • Including magnetism reduces the calculated resistivity by a factor of 3 compared to non-magnetic calculations (which overestimate due to soft phonon modes).
  • Theoretical results match experimental data well (accounting for ~75% of resistivity at 300K using PBE).
• Nickel (Ni):
  • Electron-phonon scattering contributes less than one-third of the total resistivity.
  • The majority of the resistivity arises from magnon (spin fluctuation) scattering, which is not captured by e-ph calculations alone.
  • Non-magnetic calculations for Ni yield a curve that accidentally matches experiments at low temperatures due to error cancellation, but fails to capture the correct physics.
  • Including magnetism reduces the calculated resistivity by a factor of 7.

C. Stability of Non-Magnetic Approximations

• Fe: Non-magnetic calculations result in imaginary phonon frequencies, indicating the BCC structure is unstable without magnetic ordering.
• Ni: Non-magnetic phonons are stable, but the resulting transport properties are physically incorrect regarding the source of scattering.

5. Significance and Impact

• Methodological Advancement: This work removes a major barrier in computational materials science, allowing for accurate, high-throughput screening of transport properties in magnetic materials (spintronics, magnetic memory).
• Physical Understanding: It clarifies that the transport properties of magnetic metals cannot be understood by simply "turning off" magnetism. The magnetic phase fundamentally alters lattice stability and scattering mechanisms.
• Design Guidelines:
  • For Fe-like materials, optimizing electron-phonon coupling is key to reducing resistive losses.
  • For Ni-like materials, controlling spin fluctuations (magnons) is the critical factor for improving conductivity.
• Future Directions: The framework paves the way for studying complex magnetic superconductors (e.g., iron-based superconductors, cuprates) and non-collinear magnetism (with spin-orbit coupling), which are essential for next-generation energy-efficient devices and quantum computing applications.

In summary, the paper provides a validated, first-principles toolset that reveals the critical and distinct roles of magnetism in determining the lattice stability and transport behavior of ferromagnetic metals, correcting previous misconceptions derived from non-magnetic approximations.