⚛️ phenomenology

A Mean-Field Approach to the Dielectric Response of Bulk Superconductors for Light Dark Matter Direct Detection

This paper derives the electronic dielectric function for bulk superconductors within the BCS framework and demonstrates that the standard Lindhard approximation remains a robust model for dark matter direct detection when energy depositions exceed five times the superconducting gap.

Original authors: Chengyao Gan, Hui Li, Wenjing Li, Zheng-Liang Liang, Lin Zhang, Fawei Zheng

Published 2026-09-09

📖 4 min read🧠 Deep dive

CC BY 4.0

Original authors: Chengyao Gan, Hui Li, Wenjing Li, Zheng-Liang Liang, Lin Zhang, Fawei Zheng

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

For decades, the search for dark matter has focused on the heaviest, slowest particles, hoping they would bump into heavy atomic nuclei and leave a detectable thud. But as detectors have become more sensitive, physicists have turned their gaze toward the lighter, faster end of the spectrum. These elusive particles, if they exist, would be too light to shake a heavy nucleus but might be energetic enough to nudge an electron. To catch these faint whispers, scientists are building detectors out of superconductors, materials that conduct electricity with zero resistance when cooled to temperatures near absolute zero. In these exotic states, electrons pair up and move in a synchronized dance, creating a fragile energy barrier that a passing dark matter particle could break, leaving behind a trail of excited particles that the detector can see.

However, to know if a signal is truly dark matter, researchers must first understand exactly how the detector itself reacts. When a particle hits the material, it doesn't just hit a single electron; it triggers a complex response from the entire crowd of electrons, a phenomenon known as screening. For normal metals, scientists have long used a standard mathematical tool to predict this crowd behavior. But superconductors are not normal metals; their electrons are locked in pairs, and it was unclear whether the old tool would still work in this new, paired environment. If the standard tool failed, it could mean that current experiments are either missing dark matter signals or misinterpreting noise as discovery.

A team of researchers has now stepped in to settle this question by deriving a new, more precise description of how electrons behave in a superconductor when struck by a particle. They built a detailed model based on the fundamental theory of superconductivity, accounting for the unique way electrons pair up and the specific factors that govern their movement. Instead of relying on the old, simplified tool, they calculated the response of the electron crowd from first principles, incorporating the full complexity of the superconducting state. They then compared their new, rigorous calculation against the standard tool used in current experiments, specifically looking at materials like aluminum and tungsten silicide, which are the active ingredients in some of the most advanced dark matter detectors today.

The results offer a reassuring validation for the current generation of experiments. The researchers found that for energy deposits greater than five times the superconducting energy gap—a threshold that corresponds to the minimum energy needed to break an electron pair—the standard tool works remarkably well. In the vast majority of the scenarios relevant to detecting dark matter, the simplified model produces results that are nearly identical to the complex, new calculation. This means that the experiments currently running, which rely on the older approximation, are likely interpreting their data correctly. The team confirmed that for the specific materials used in these detectors, the difference between the two methods is negligible in the energy ranges where dark matter is expected to leave a mark.

There are, however, narrow boundaries where the old tool begins to drift from reality. The new model shows that the standard approximation breaks down very close to the minimum energy required to break a pair, and in a specific region where the electron response in a normal metal would vanish but remains active in a superconductor. Yet, these edge cases do not undermine the main conclusion for the bulk of the search. The study suggests that for the typical dark matter particles these detectors are hunting, the simplified approach is robust enough to be trusted.

This work does more than just confirm a calculation; it clears the path for future upgrades. As detectors become larger and more sensitive, aiming to catch even lighter dark matter particles, scientists will need to be certain that their models of the detector's behavior are flawless. By proving that the standard approximation holds up under rigorous scrutiny for the most critical energy ranges, the researchers have provided a solid foundation for the next decade of discovery. They have shown that while the electrons in a superconductor are indeed behaving in a unique, paired way, the simple rules used to describe them are sufficient to guide the search for the universe's most mysterious substance. The only time the rules need to be rewritten is in the very specific, low-energy limits that are currently beyond the reach of most experiments, leaving the main search for dark matter on a firm and verified footing.

Technical Summary: A Mean-Field Approach to the Dielectric Response of Bulk Superconductors for Light Dark Matter Direct Detection

Problem Statement
Direct detection of sub-GeV dark matter (DM) using condensed matter targets, particularly superconducting sensors (e.g., Kinetic Inductance Detectors, Transition-Edge Sensors, and Superconducting Nanowire Single-Photon Detectors), relies on accurately modeling the interaction between DM particles and the detector's electron system. The scattering rate is governed by the dynamic structure factor, which is derived from the imaginary part of the inverse dielectric function. Current experimental analyses for superconducting detectors typically employ the Lindhard dielectric function—originally derived for normal metals within the Random-Phase Approximation (RPA)—to model in-medium screening effects. This approach assumes that the normal-state dielectric response adequately approximates the screening in the superconducting state. However, the validity of this approximation within the superconducting phase, where Cooper pairing and energy gaps fundamentally alter the electronic spectrum, has not been rigorously tested.

Methodology
The authors generalize the RPA framework to the superconducting ground state using the Bardeen-Cooper-Schrieffer (BCS) theory.

Formalism: The electron field is organized using Nambu spinors to naturally couple electron and hole degrees of freedom. The Coulomb interaction is expressed in this basis, and the BCS mean-field Hamiltonian introduces pairing correlations.
Bogoliubov Transformation: The system is diagonalized via a Bogoliubov unitary transformation, yielding quasiparticle (QP) excitations with a dispersion relation $E_p = \sqrt{\epsilon_p^2 + \Delta^2}$ and characteristic coherence factors $(u_p, v_p)$ .
Dielectric Function Derivation: The authors derive the electronic dielectric function $\epsilon_{BCS}(Q, \omega)$ for bulk superconductors. This involves calculating the QP polarization function $\Pi_{QP}(Q, \omega)$ by summing bubble diagrams to all orders in the RPA. Crucially, this derivation incorporates the full Bogoliubov quasiparticle coherence factors, which project the electron-density vertex onto the QP basis.
DM Interaction: A generic model is considered where spin-1/2 DM particles couple to electron density via a Yukawa potential. The tree-level scattering process involves breaking a Cooper pair to excite two QPs. The scattering amplitude is renormalized by the BCS dielectric function.
Numerical Comparison: The derived $\epsilon_{BCS}(Q, \omega)$ is numerically compared against the standard Lindhard function $\epsilon_{Lin}(Q, \omega)$ for two materials relevant to current experiments: Aluminum and Tungsten Silicide (WSi, used in the QROCODILE experiment). Parameters include specific energy gaps ( $\Delta$ ) and Fermi energies ( $E_F$ ) for each material.

Key Contributions

Derivation of $\epsilon_{BCS}$ : The paper provides the first systematic derivation of the electronic dielectric function for bulk superconductors within the BCS framework, explicitly including full coherence factors in the RPA.
Validation of Approximations: It offers a quantitative assessment of the Lindhard function's validity in the superconducting state by comparing it directly with the derived BCS dielectric function.
Kinematic Analysis: The study maps the validity of the approximation against the kinematic regimes relevant for halo DM scattering, considering typical velocities ( $v_\chi \approx 10^{-3}c$ ) and various DM masses.

Results

Agreement Regime: A systematic comparison reveals that $|\epsilon_{BCS}|^2$ and $|\epsilon_{Lin}|^2$ agree closely over most of the parameter space of interest. Specifically, for energy depositions $\omega \gtrsim 5\Delta$ , the Lindhard function provides a robust approximation for the in-medium response of superconducting targets.
Discrepancy Regimes: Significant deviations occur in two specific regions:
1. Near the superconducting gap edge ( $\omega \approx 2\Delta$ ), where QP production is just becoming possible.
2. In the region where the imaginary part of the Lindhard function vanishes (above the red dotted curve $\omega \approx v_F Q$ ), while the BCS function remains finite.
Applicability: For bulk superconductors and momentum transfers $Q \gtrsim 0.2$ keV (corresponding to length scales $\lesssim 1$ nm), the bulk approximation holds. This validates the use of the Lindhard function for DM masses $m_\chi \gtrsim 0.1$ MeV in current superconducting detector configurations.

Significance and Claims
The paper claims that the widely used practice of employing the normal-state Lindhard function to model screening in superconducting DM detectors is justified for energy transfers significantly above the superconducting gap ( $\omega \gtrsim 5\Delta$ ). This finding supports the reliability of current constraints on light DM-electron scattering derived from experiments like QROCODILE.

The authors maintain a modest stance regarding the limits of their mean-field approach. They acknowledge that while the mean-field BCS description is sufficient for direct detection scenarios with energy transfers above $5\Delta$ (where the relative error is at most $\sim 4\%$ ), it may not capture collective dynamics in the long-wavelength limit ( $Q \ll 1/\xi_0$ ). In this regime, the phase degree of freedom of the order parameter becomes prominent, requiring a description beyond mean-field theory to account for phase-field-electron coupling. The authors explicitly state that a detailed analysis of this crossover regime is left for future work.

Technical Summary: A Mean-Field Approach to the Dielectric Response of Bulk Superconductors for Light Dark Matter Direct Detection

More like this