🔬 materials science

Many-body interactions in the dielectric theory of stopping power of solids for classical and quantum projectiles

This paper incorporates wave-vector and frequency-dependent exchange-correlation effects into the dielectric theory of stopping power for both classical and quantum projectiles in crystals, demonstrating that these many-body interactions significantly improve agreement with experimental data at low velocities while proving negligible at high velocities.

Original authors: Vladimir U. Nazarov, Vyacheslav M. Silkin

Published 2026-07-15

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Vladimir U. Nazarov, Vyacheslav M. Silkin

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 a high-speed projectile, like a tiny bullet made of a proton or an electron, zooming through a solid crystal. As it flies, it drags a wake behind it, knocking electrons around and losing energy. This energy loss is called "stopping power." For decades, scientists tried to predict exactly how fast these particles would slow down, but their old maps were missing some crucial terrain.

In this study, Vladimir Nazarov and Vyacheslav Silkin decided to redraw the map by adding two missing ingredients: the crystalline structure of the solid (how the atoms are neatly stacked) and many-body interactions (the complex, messy dance where electrons don't just react to the bullet, but also to each other).

The Old Map vs. The New Map

Think of the old way of calculating stopping power as trying to predict how a boat slows down in a perfectly smooth, featureless ocean (a model called "jellium"). It's a nice, simple idea, but real solids aren't smooth oceans; they are like coral reefs with intricate patterns. Furthermore, the old models often ignored how the water molecules (electrons) interact with each other, treating them as if they were swimming solo.

The authors' new approach is like switching to a high-definition satellite view. They used a sophisticated framework called "dielectric theory" but upgraded it with a dynamic "exchange-correlation kernel" (a fancy math tool, labeled $f_{xc}$ , that accounts for those electron-to-electron interactions). They tested this on both heavy, classical projectiles (like protons) and light, quantum ones (like electrons).

What They Found: The Sweet Spot

The results were a bit like finding the perfect gear in a car.

1. The Low-Speed Zone (The "Slow Lane"):
When the projectile is moving slower than its maximum stopping speed, the old models were missing the mark.

The Problem: The simple "jellium" model (the smooth ocean) consistently underestimated how much the particle slows down. It was too optimistic.
The Fix: When the authors added the "many-body" interactions (the electron dance) and the crystal structure, the predictions got much better. For targets like aluminum and silicon, the new theory lined up much closer with real-world experiments. It's as if they finally accounted for the friction caused by the coral reef's shape and the water's internal turbulence.
The Catch: However, this new map hit a wall with alkali metals like Lithium and Rubidium. At low speeds, the theory broke down completely. The authors suggest this is because, in these specific metals, the projectile gets "stuck" in a temporary trap (a bound state) that the standard linear theory simply cannot see. It's like trying to predict a boat's speed while ignoring that it just ran aground on a hidden sandbar.

2. The High-Speed Zone (The "Fast Lane"):
When the projectile is zooming very fast, things change again.

The Surprise: In this high-velocity regime, the complex "many-body" interactions (the $f_{xc}$ kernel) turned out to be negligible. The authors proved analytically that for very fast particles, you don't need to worry about the complicated electron dance; the simple model works just fine.
The Reality Check: Even with the new theory, the calculations still underestimated the stopping power for aluminum and silicon at high speeds. Why? Because at these speeds, the projectile is hitting the "core" electrons deep inside the atoms—electrons that are tightly bound and not part of the general "sea" of electrons the model usually tracks. The authors note that their current tools couldn't fully simulate these deep-core interactions because it would require too much computing power.

The Quantum Twist

The paper also looked at electrons as projectiles. Here, the rules of quantum mechanics (where particles act like waves) come into play. The authors included a special factor to account for the "exchange" between the incoming electron and the target electrons. The results were generally good, showing that the theory holds up well for electrons, too, though the many-body effects weren't as critical here as they were for protons.

The Bottom Line

This study didn't just tweak the numbers; it clarified where and why our current theories work and where they fail.

It suggests that for most solid targets (like aluminum and silicon) at moderate speeds, you absolutely need to include both the crystal structure and the complex electron interactions to get the right answer.
It argues against the idea that these complex interactions matter at very high speeds; there, they fade into the background.
It highlights a failure for alkali metals at low speeds, indicating that the standard linear theory isn't enough to explain what happens when a projectile gets trapped in a bound state.

The authors didn't solve the whole mystery of stopping power, but they successfully mapped out the terrain where the old models were wrong and showed us exactly where the new, more complex models shine—and where they still hit a wall.

Technical Summary: Many-body interactions in the dielectric theory of stopping power of solids for classical and quantum projectiles

Problem Statement
The stopping power (SP) of solids—the energy loss of a moving charge per unit path length—is a fundamental problem in condensed matter physics. While the dielectric theory of SP, rooted in the random phase approximation (RPA) and the homogeneous electron gas (HEG) model, has been established for decades, it often neglects many-body dynamic interactions, specifically the exchange and correlation (xc) effects. Furthermore, standard approaches frequently rely on the adiabatic local density approximation (ALDA), which fails to capture spatial and temporal non-locality in the electronic response. Additionally, the applicability of the dielectric approach to real crystalline solids, as opposed to idealized jellium models, and its validity for both classical and quantum projectiles require rigorous verification against experimental data, particularly in regimes where the first Born approximation may break down.

Methodology
The authors implement a comprehensive dielectric theory of stopping power that integrates several advanced theoretical components:

Framework: The study utilizes the linear-response time-dependent density functional theory (TDDFT) within the dielectric formalism. The projectile-target interaction is treated perturbatively (first Born approximation).
Many-Body Interactions: The core innovation is the inclusion of the wave-vector and frequency-dependent exchange-correlation kernel, $f_{xc}(q, \omega)$ , to account for dynamic many-body effects. The authors employ the constraint-based revised modified Constantin–Pitarke (rMCP07) parametrization, considered accurate across fluid phase densities, alongside the ALDA for comparison.
Projectile Treatment: The theory distinguishes between classical and quantum projectiles. For finite-mass projectiles, the kinematic constraints are explicitly handled via quantum mechanical integration limits involving $\hbar$ , distinguishing them from the classical limit.
Crystalline Structure: To address real solids, the theory is generalized from the jellium model to periodic crystals. This involves calculating the microscopic dielectric matrix $\epsilon_{GG'}(q, \omega)$ using the all-electron full-potential linearized augmented-plane wave (LAPW) code ELK. The spatial non-locality of $f_{xc}$ is approximated using the average valence electron density of the unit cell.
Scope: Calculations are performed for various targets (Al, Si, Diamond, Be, Li, Rb) and projectiles (protons, positrons, electrons) across a range of velocities, comparing results against available experimental datasets.

Key Contributions and Results

Enhanced Dielectric Theory: The paper presents the dielectric theory of SP in its "unabridged form," combining the accurate evaluation of the dielectric matrix for periodic solids with non-adiabatic linear-response TDDFT to include dynamic xc effects.
Role of xc Effects in Crystals: The authors find that xc effects play a more pronounced role in slowing charges in crystalline solids than in the jellium model. Specifically, for velocities below the stopping maximum, the inclusion of xc (particularly via rMCP07) significantly improves the agreement between theory and experiment for most solid targets (e.g., Al, Si, Diamond).
Velocity Dependence:
- Low to Intermediate Velocities: The inclusion of xc effects corrects the systematic underestimation of SP seen in RPA and the overestimation seen in ALDA, bringing theoretical predictions closer to experimental data.
- High Velocities: In the high-velocity regime, the role of dynamic xc is found to be negligible. The authors substantiate this analytically within the jellium model, noting that discrepancies in this regime are likely due to the excitation of core electrons, which the jellium model cannot capture, rather than a failure of the xc treatment.
Breakdown of the Dielectric Approach: The study identifies specific limitations. For alkali metals (Lithium and Rubidium) at low velocities, the dielectric theory breaks down. This is attributed to the formation of bound states of the pseudo-atom (occurring at $r_s \gtrsim 2$ a.u. for protons), a process the first Born approximation cannot capture. Consequently, the theory fails to describe SP for Li and Rb in this regime, whereas it remains marginally valid for Al and Si where the bound state is shallow.
Electron Projectiles: For electron projectiles, the theory shows good agreement with experiment, though xc effects in the target response play a less dominant role compared to proton projectiles. The rapid disappearance of SP at low velocities for electrons aligns with nonlinear theory predictions for the friction coefficient.

Significance
The paper establishes the bounds of applicability for the dielectric approach to stopping power, demonstrating that while the inclusion of dynamic many-body interactions and crystallinity significantly refines theoretical predictions for a wide range of solids, the approach is not universally applicable. The work highlights that the failure of the dielectric theory in low-velocity alkali metals is an inherent limitation of the perturbative treatment rather than a deficiency in the xc kernel approximation. By validating the theory against experimental data and analytically proving the negligible role of xc at high velocities, the authors provide a robust framework for calculating stopping power that balances computational feasibility with physical accuracy, while clearly delineating where more complex, non-linear theories are required.

The Old Map vs. The New Map

What They Found: The Sweet Spot

The Quantum Twist

The Bottom Line

More like this