🔬 materials science

Opposite impact of thermal expansion and phonon anharmonicity on the phonon-limited resistivity of elemental metals from first principles

This study demonstrates that incorporating the opposing effects of thermal expansion, which enhances electron-phonon coupling and overestimates resistivity, and phonon anharmonicity, which reduces it, provides a more accurate first-principles description of electrical resistivity in elemental metals like Pb, Nb, and Al.

Original authors: Ao Wang, Junwen Yin, Félix Antoine Goudreault, Michel Côté, Olle Hellman, Samuel Poncé

Published 2026-02-04

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Ao Wang, Junwen Yin, Félix Antoine Goudreault, Michel Côté, Olle Hellman, Samuel Poncé

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict how easily electricity flows through a metal, like a crowd of people trying to walk through a busy hallway. The "resistance" they feel is the electrical resistivity. For a long time, scientists have used computer models to predict this, but they often missed two crucial factors that happen when the metal gets hot: the hallway getting wider (thermal expansion) and the walls starting to wiggle unpredictably (phonon anharmonicity).

This paper, by Wang and colleagues, reveals that these two missing factors are actually opposites that cancel each other out. If you ignore both, you get a lucky guess that happens to be right. If you include only one, you get a very wrong answer. You have to include both to get the true picture.

Here is the breakdown using simple analogies:

1. The Two Opposing Forces

Factor A: Thermal Expansion (The Hallway Gets Wider)
When a metal heats up, it physically expands, like a balloon inflating. In the world of electrons, this means the "hallway" they walk through gets stretched out.

The Paper's Finding: This stretching actually makes it harder for electrons to move. It's like stretching a rubber band; the atoms get further apart, and the electrons bump into things more often.
The Result: If you only calculate this effect, your computer predicts that the metal will become a terrible conductor (high resistance). In fact, for Lead (Pb), this alone made the predicted resistance nearly double what was actually measured at high temperatures.

Factor B: Phonon Anharmonicity (The Walls Start Wiggling)
"Phonons" are vibrations of the atoms. Usually, scientists pretend these atoms vibrate like perfect springs (back and forth in a straight line). But in reality, as things get hot, the atoms get "anharmonic"—they start wiggling in a messy, non-linear way, almost like a jelly shaking.

The Paper's Finding: This messy wiggling actually stiffens the vibrations (a phenomenon called "phonon hardening"). It's as if the chaotic movement of the atoms somehow organizes the path for the electrons, making it easier for them to slip through.
The Result: If you only calculate this effect, your computer predicts the metal will be too conductive (low resistance).

2. The "Perfect Cancellation" (Lead and Aluminum)

The authors tested this on Lead (Pb) and Aluminum (Al). They found a fascinating "tug-of-war":

Thermal Expansion tries to increase resistance.
Anharmonicity tries to decrease resistance.
The Magic: These two forces are almost equal in strength but point in opposite directions. They cancel each other out perfectly.

The Analogy: Imagine you are trying to walk down a hallway.

Thermal Expansion is like someone stretching the hallway so the floor tiles are further apart, making you trip more often.
Anharmonicity is like the walls suddenly vibrating in a way that creates a smooth, sliding path for you.
The Reality: The stretching makes you trip, but the sliding walls help you recover. The net result is that you walk at your normal speed.

If you only looked at the stretching, you'd think you'd fall down. If you only looked at the sliding walls, you'd think you'd fly. But when you look at the whole picture, you just walk normally. This is why previous models that ignored both factors accidentally got the right answer—they missed two errors that canceled each other out.

3. The Exception: Niobium (The Complex Dance)

The team also tested Niobium (Nb), and the story was slightly different.

In Niobium, the "walls" (the electron energy levels) have a very complex shape (a "nesting Fermi surface").
When the metal heats up, the stretching and the wiggling don't happen in the same places. The stretching affects one part of the hallway, while the wiggling affects another.
The Result: They don't cancel out perfectly. The "wiggling" (anharmonicity) is stronger, so the metal ends up conducting slightly better than the "stretching" alone would suggest, but not quite as perfectly as in Lead or Aluminum.

The Bottom Line

For a long time, scientists calculated electrical resistance by ignoring how metals expand and how atoms wiggle messily when hot. They got lucky because the errors canceled out.

This paper proves that to truly understand how metals conduct electricity at high temperatures, you must include both the expansion and the messy wiggling. When you do, the computer models finally match real-world experiments perfectly, showing us that nature often balances opposing forces to create stability.

Technical Summary: Opposite Impact of Thermal Expansion and Phonon Anharmonicity on Phonon-Limited Resistivity

Problem Statement
Understanding electrical resistivity in metals at finite temperatures remains a central challenge in quantifying charge transport. While current first-principles calculations based on the Boltzmann transport equation (BTE) often agree with experiments, they frequently neglect two critical physical effects: thermal expansion (TE) and phonon anharmonicity. Recent studies indicated that neglecting TE leads to a significant overestimation of resistivity in lead (Pb) at high temperatures (e.g., +92% at 610 K), a discrepancy that could not be resolved by adjusting spin-orbit coupling or exchange-correlation functionals alone. This work addresses the need to explicitly incorporate both TE and phonon anharmonicity to achieve a complete description of resistivity in elemental metals.

Methodology
The authors employ a first-principles approach combining density functional theory (DFT) with the iterative Boltzmann transport equation (IBTE).

Thermal Expansion (TE): The temperature-dependent lattice parameters are determined by calculating the Helmholtz free energy $F(V, T)$ , which includes both electronic contributions (clamped-ion energy and thermal electronic excitation) and phonon contributions within the quasi-harmonic approximation (QHA). This is performed using Quantum ESPRESSO with a third-order Birch-Murnaghan isothermal fit.
Phonon Anharmonicity: To capture anharmonic effects, the study utilizes the Temperature Dependent Effective Potential (TDEP) method. This involves stochastic sampling of atomic configurations generated from the canonical ensemble to determine effective interatomic force constants (IFCs) via linear least-squares solutions. The Hamiltonian is expanded to the second order with respect to atomic displacements.
Transport Calculations: The electrical resistivity is calculated by solving the linearized BTE using the EPW code, interfaced with TDEP for anharmonic IFCs. Electron-phonon coupling (EPC) matrix elements are interpolated using Wannier90. Spin-orbit coupling (SOC) is included for Pb. The scattering rates are computed considering electron-phonon interactions, neglecting phonons with frequencies below 5 cm $^{-1}$ .

Key Contributions and Results
The study investigates the elemental metals Pb, Nb, and Al, revealing a distinct antagonistic relationship between TE and phonon anharmonicity:

Opposing Mechanisms in Lead (Pb):
- Thermal Expansion Effect: Including TE alone causes significant phonon softening, particularly at the $q=X$ point. This softening enhances electron-phonon scattering rates, leading to a drastic overestimation of resistivity (up to 92% at 610 K) and even predicting thermodynamic instability near the melting point.
- Anharmonicity Effect: Conversely, including phonon anharmonicity induces phonon hardening with increasing temperature. This reduces the electron-phonon scattering rates, resulting in an underestimation of resistivity when TE is ignored.
- Cancellation and Agreement: When both effects are included simultaneously, the opposing impacts on phonon frequencies largely cancel out. The resulting phonon dispersion matches experimental neutron scattering data perfectly. Consequently, the calculated electrical resistivity aligns closely with experimental values across the entire temperature range (10–610 K), resolving the previous overestimation issue.
Microscopic Origin:
Analysis of electron scattering rates ( $\tau^{-1}_{nk}$ ) at 310 K and 610 K shows that TE increases scattering by approximately 20% at 310 K, while the inclusion of anharmonicity reduces this rate back to levels comparable to calculations ignoring both effects. Spectral decomposition reveals that when both effects are considered, transport electrons in Pb are primarily scattered by phonons with energies between 3 and 5 meV.
Generality and Exceptions (Al and Nb):
- Aluminum (Al): Similar to Pb, Al exhibits opposing effects from TE and anharmonicity that largely cancel each other, yielding results that match experiments well below 500 K.
- Niobium (Nb): Nb presents a different case due to its Kohn anomaly related to Fermi surface nesting. While TE and anharmonicity still exert opposite influences, they do not fully cancel. Here, the anharmonic effect is stronger, bringing results closer to experiment, though a slight overestimation remains, attributed to saturation effects. The lack of perfect cancellation is attributed to the highly dispersive nature of phonon shifts in Nb, where volume changes (rigid bond stretching) and temperature changes (Fermi surface smearing) affect different regions of the Brillouin zone.

Significance
The paper establishes that thermal expansion and phonon anharmonicity exert antagonistic impacts on the electrical resistivity of simple metals. The authors demonstrate that neglecting either effect leads to significant errors: TE alone causes overestimation, while anharmonicity alone causes underestimation. The primary claim is that a complete and accurate theoretical description of phonon-limited resistivity requires the simultaneous inclusion of both effects. This approach successfully reproduces experimental data for Pb, Nb, and Al, correcting previous discrepancies and providing a more robust framework for quantifying charge transport in metals at finite temperatures.

1. The Two Opposing Forces

2. The "Perfect Cancellation" (Lead and Aluminum)

3. The Exception: Niobium (The Complex Dance)

The Bottom Line

Technical Summary: Opposite Impact of Thermal Expansion and Phonon Anharmonicity on Phonon-Limited Resistivity

More like this