Sommerfeld enhancement from unstable final-state particles in dark matter annihilation

arXiv:2603.02647 · hep-ph, astro-ph.CO · Submitted 2026-03-03 · Read on arXiv

Listen

Radio episode about this paper

Transcript

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: I'm Vera, and with me are Jocelyn and Subrahmanyan, guest researcher.

Jocelyn: Today's paper: "Sommerfeld enhancement from unstable final-state particles in dark matter annihilation".

Vera: This paper investigates the Sommerfeld enhancement (SE) of dark matter annihilation when it proceeds into heavier, unstable final-state particles, a phenomenon that significantly impacts predictions for dark matter relic abundance.

Jocelyn: First, who's behind it and why it matters.

Title and authors: Vera: Well, we're starting today with a paper titled "Sommerfeld enhancement from unstable final-state particles in dark matter annihilation." It looks like it tackles a really specific problem where dark matter annihilates into heavier, unstable particles. It suggests that when these products interact over long distances, their wave functions get distorted, which can boost the annihilation cross section significantly.

Jocelyn: That sounds complicated but interesting from an observational standpoint because it links microscopic particle physics to something we might be able to see in cosmological observations. The authors are Abe, Sato, and Yamanaka from Tokyo University of Science and Osaka University. I wonder how this enhancement factor actually translates into anything measurable in the sky or through relic density calculations.

Subrahmanyan: From a theoretical perspective, the focus on unstable final states is key because it introduces a dependence on the particle's lifetime, which fundamentally alters how we calculate these non-perturbative corrections. This isn't just about standard annihilation; it's about what happens when things decay while they are still interacting.

Vera: Exactly, and I think the title tells us the core issue is how long-range interactions affect particles that aren't stable. It suggests that we need a more complete way to handle these unstable states in our models of dark matter behavior.

Jocelyn: So, when we look at the results, I'm anticipating some constraints on what kind of dark matter particle might be relevant based on how strongly these long-range forces affect the annihilation cross section. It feels like it sets a new benchmark for how we predict relic abundance.

Subrahmanyan: Indeed, this work is important because it moves beyond simpler treatments where you just introduce a cutoff velocity to handle instability. This paper seems to look at the full physics, including those processes with off-shell final state particles that might be missed by approximations like the cutoff method.

Vera: That's what I'm hoping we hear about; a full picture that includes everything near the threshold. It suggests our current predictions for dark matter relic abundance might need serious re-evaluation if these long-range effects are significant.

Jocelyn: It’s exciting because it points toward specific mass ratios and potential structures that will actually lead to observable deviations in cosmological predictions, which is something we can look for in the data.

The paper's summary: Vera: So, digging into the summary of "Sommerfeld enhancement from unstable final-state particles in dark matter annihilation," they are essentially looking at how long-range forces between annihilation products distort their wave functions when they get close to the kinematic threshold. They find that this distortion can cause a substantial enhancement in the annihilation cross section.

Jocelyn: It seems like the core mechanism is that if these products experience long-range interactions, their wave functions aren't just simple plane waves anymore, and this change directly boosts the annihilation rate. The paper shows that this enhancement factor depends heavily on whether those final state particles have a narrow or a large decay width.

Subrahmanyan: That dependence on the decay width is what makes it tricky, because for narrow widths, you get bound states that cause a resonant effect, whereas for large widths, you automatically account for processes involving off-shell particles. This distinction is crucial for understanding the physics in these scenarios.

Vera: So the summary highlights that the lifetime of the final state particle dictates whether we see bound state effects or just general off-shell contributions, which is a vital piece of information for model building. It shows this factor isn't just a minor correction but can be quite substantial when resonances are present.

Jocelyn: And I think what's most telling is that the formulation automatically includes both bound state effects and those contributions from annihilation processes involving off-shell final state particles, which is a big plus for completeness.

Subrahmanyan: Precisely, this unified treatment is superior to previous work that often had to rely on introducing a cutoff velocity, like v cut = p /m two which only includes the enhancement if the dark matter velocity is larger than that threshold.

Vera: That comparison with the cutoff method is really important because it shows where those simpler methods fall short, especially since the full formulation allows for annihilation processes to be kinematically allowed even when energy is below the threshold due to off-shell effects.

Jocelyn: That means we can predict cross sections that are non-zero in regimes where simpler models would say they should be zero, which opens up new avenues for exploring parameter space.

The paper's improvements: Vera: The paper points out several ways this formulation improves upon existing research, specifically by using Schrödinger equations that explicitly include the decay width of the annihilation products. This is a step up from just treating the instability as a simple cutoff.

Jocelyn: They show that this detailed approach allows them to derive an expression for the annihilation cross section, -sigma v rel = two S f(E two), where S f(E two) is the enhancement factor that captures both bound state effects and off-shell contributions.

Subrahmanyan: The paper notes that for a narrow decay width, bound states create a resonant enhancement near the threshold, while for a large decay width, it automatically covers the off-shell particles. This distinction shows how the physics naturally separates itself based on the final state's lifetime.

Vera: And they apply this to specific potentials, like the attractive Coulomb potential and the attractive Hulthén potential, showing that they can identify resonant enhancements at energies corresponding to specific bound states, such as 1s, 2s, and 3s states.

Jocelyn: I'm interested in the condition they set for consistency: the annihilation rate of the chi two two state into chi one one must be much smaller than the total decay width. This consistency check is a rigorous way to ensure their method is physically sound.

Subrahmanyan: And they compare this to previous work by noting that the full formulation is consistent only if two to one where two to one is the partial decay rate of the bound state into chi one one ensuring the relevant annihilation rate is dominated by that specific decay. This tight constraint is important for validating their approach.

Vera: It seems like a really robust set of checks, moving from simple velocity cutoffs to this full quantum mechanical treatment involving coupled Schrödinger equations and careful handling of the decay width.

Conclusion: Jocelyn: So, wrapping up the paper "Sommerfeld enhancement from unstable final-state particles in dark matter annihilation," the main implication is that solving these coupled Schrödinger equations including the decay width is necessary to correctly capture resonant effects. They confirm that this method works consistently when two to one which means we can trust their predictions for dark matter relic abundance in regimes where resonances are expected.

Subrahmanyan: From the cosmic perspective, the impact is that for attractive potentials like the Coulomb potential, they show that the enhancement is indeed present in certain energy ranges and isn't negligible. This means we might be seeing deviations from previous methods for mass ratios above m two/m one one point zero zero one two, which directly alters the predicted mass spectra that reproduce the observed relic density.

Vera: That alteration in the predicted mass spectra is significant because it means that our understanding of what kind of dark matter particle we're looking for to match observations might need adjustment based on these results from this paper. It shows the final-state SE isn't just a small tweak; it can change the fundamental predictions.

Jocelyn: So, in short, this research provides a more rigorous framework for calculating annihilation cross sections when instability is involved, and it offers concrete examples of how to predict resonant effects using both bound states and off-shell particles.

Subrahmanyan: I just want to emphasize that the consistency condition two to one is what makes this method reliable for multi-component models, suggesting it can be extended beyond these two-body scenarios.

Vera: It’s a solid piece of work that solidifies the necessity of this complex treatment when dealing with unstable final states in dark matter annihilation calculations. We'll take a moment to process all this information before we move on to the next topic.

Tomohiro Abe, Ryosuke Sato, Takumu Yamanakab

Department of Physics and Astronomy, Faculty of Science and Technology, Tokyo University of Science · Department of Physics, The University of Osaka

hep-ph, astro-ph.CO

Submitted: 2026-03-03

Updated: 2026-09-29

Comments: 22 pages, 9 figures, v2: published version

Journal ref: JHEP 09 (2026) 275

DOI: 10.1007/JHEP09(2026)275

License: http://creativecommons.org/licenses/by/4.0/

Importance score: 79/100

The gist: This paper investigates the Sommerfeld enhancement (SE) of dark matter annihilation when it proceeds into heavier, unstable final-state particles, a phenomenon that significantly impacts predictions

Key concepts

Sommerfeld enhancement (SE)
This is a phenomenon where long-range interactions between particles in a dark matter annihilation process distort their wave functions near the kinematic threshold, leading to a substantial increase in the annihilation cross section.
Unstable final-state particles
These are the heavier particles produced when dark matter annihilates. Their instability, characterized by a decay width, is crucial because it dictates whether effects are seen as bound states or off-shell contributions during interactions.
Cutoff velocity method
This is a simpler approximation used in previous work to handle particle instability. It only includes the enhancement if the dark matter velocity exceeds a specific threshold, which this paper aims to improve upon.

Terminology

Summary

This paper investigates the Sommerfeld enhancement (SE) of dark matter annihilation when it proceeds into heavier, unstable final-state particles, a phenomenon that significantly impacts predictions for dark matter relic abundance. The research focuses on how long-range interactions between these annihilation products can distort their wave functions to enhance the annihilation cross section, and how this effect is modified by the decay width of the final state. This is crucial because it provides a mechanism where bound states of annihilation products can lead to resonant effects that significantly alter cosmological predictions.

Formulation of Final-State Sommerfeld Enhancement (SE)

The study formulates the DM annihilation process using Schrödinger equations that explicitly include the decay width of the annihilation products. For a two-to-two process involving DM particles with mass 1 and annihilation products with mass 2, the wave functions are governed by coupled Schrödinger equations (Eqs. 2.1 and 2.2). The cross section is then derived from the probability current, leading to the expression:

- vrel = dP/dt / j in 1 = 2 Im(uΨ∗1(0)Ψ2(0)).

The annihilation cross section including the final-state SE is expressed as:

- σvrel = a˜2 Sf (E2, Γ).

The enhancement factor, Sf (E2, Γ), is defined as the ratio of the annihilation cross section with and without the long-range potential. This formulation automatically includes both bound state effects and contributions from annihilation processes with off-shell final state particles.

Impact of Decay Width and Bound States

The paper finds distinct behaviors depending on whether the decay width is narrow or large. For a narrow decay width, bound states of the annihilation products can enhance the annihilation cross section through a resonant effect. Conversely, for a large decay width, the formulation automatically includes processes with off-shell final state particles. The consistency of this formulation requires that the annihilation rate of χ2χ¯2 state into χ1χ¯1 is much smaller than the total decay width Γ.

Comparison with Previous Work and Cutoff Method

Previous studies often treated instability by introducing a cutoff velocity, where SE is included only if the DM velocity is larger than a cutoff velocity, vcut = pΓ/m2. The paper compares its formulation to this cutoff method, noting that the cutoff method neglects annihilation cross sections below the kinematical threshold. The full formulation shows that for negative E2 (below threshold), annihilation processes involving off-shell particles are kinematically allowed, leading to a non-zero cross section, which is missed by the cutoff method.

Examples and Potential Dependence

The formulation is applied to two specific long-range potentials:

  1. Attractive Coulomb potential: V(r) = -α/r. The analysis shows that the full result exhibits resonant enhancement at specific energies corresponding to the binding energies of 1s, 2s, and 3s bound states of χ2χ¯2.

  2. Attractive Hulthén potential: VH(r) = −αm∗e−m∗r/1 − e−m∗r. This potential allows for a finite number of bound states, with the condition for at least one bound state being m∗ ≤ 2µ2α.

Influence on Dark Matter Relic Abundance

The final-state SE significantly modifies the prediction of the dark matter relic abundance, especially when resonances are present. For attractive potentials like the Coulomb potential, the full result shows that it is indeed enhanced in this energy range and the final-state SE is not negligible. This enhancement can lead to a deviation between our formulation and previous methods for specific mass ratios (e.g., m2/m1 ≳ 1.0012), meaning the predicted mass spectra for DM particles that reproduce the observed relic density are altered. The existence of resonances affects the thermally averaged annihilation cross section, leading to a smaller predicted relic abundance in some regimes compared to the cutoff method result.

Conclusion and Consistency

The study concludes that solving Schrödinger equations including the decay width is necessary to correctly evaluate resonant effects. The formulation is consistent only if Γ ≫ Γ2→1, where Γ2→1 is the partial decay rate of the bound state into χ1χ¯1, ensuring that the relevant annihilation rate is dominated by the partial decay rate of the bound state. The method can be extended to multi-component models, such as accidental composite DM.

Improvements for AI systems

Here are specific improvements to AI systems that could be derived from the scientific concepts presented in this paper, categorized by capability:

) 1. Enhanced Dark Matter Phenomenology Simulation (High-Fidelity Physics Modeling)

The core of the paper lies in accurately modeling non-perturbative quantum mechanical effects (Sommerfeld Enhancement, SE) in particle annihilation near kinematic thresholds, which is crucial for predicting dark matter relic abundance.

Improvements for AI systems:

  • AI can be trained to solve the coupled Schrödinger equations (Equations 2.1 and 2.2) with complex potentials (Coulomb, Hulthén) and finite decay widths analytically or numerically at high precision, rather than relying on simplified cutoff methods (Equation 2.34).

  • The system can incorporate the full Green's function formalism (Equations 13-15) to automatically account for both bound state effects and off-shell final states simultaneously, which is superior to sequential approximations.

  • The AI can predict the precise dependence of the annihilation cross section on key physical parameters—specifically mass ratios and coupling strengths—by mapping the complex relationships derived in Section 4 (e.g., how resonance peaks shift based on potential strength or decay width).

Improved System Capability:

  • Predicting dark matter relic abundance with a level of precision that accounts for resonant enhancements (bound states) and off-shell effects, leading to significantly more constrained and accurate models of DM particle masses and interaction strengths.

) 2. Non-Perturbative Cross-Section Estimation (Resonance Detection in High-Energy Collisions)

The paper demonstrates how bound states lead to resonant enhancements in annihilation cross sections near thresholds, which is analogous to processes like top-antitop pair production near threshold (Eqs. 17–21).

Improvements for AI systems:

  • AI can be deployed as a tool for resonance hunting in high-energy collider data simulations. Instead of relying on standard perturbative cross-section calculations, the AI can use the derived analytic forms (like Equation 4.8) to identify regions where non-perturbative enhancements are expected based on the mass ratio and potential structure.

  • The system can distinguish between true resonant enhancements due to bound states and artifacts from simple kinematic cutoffs, which is essential for interpreting experimental data that might otherwise be misinterpreted.

Improved System Capability:

  • Developing advanced signal processing algorithms for experimental dark matter searches that can identify subtle, non-perturbative signatures (like specific resonant peaks in energy spectra) that would be missed by traditional perturbative cross-section models.

) 3. Model Parameter Space Exploration and Constraint Generation (Automated Discovery)

The paper shows how the existence of bound states critically affects the relic abundance prediction, especially as a function of temperature and mass ratios (Figures 10, 6).

Improvements for AI systems:

  • AI can be used for automated parameter scanning. Given experimental constraints on DM density or other cosmological observables, the AI can rapidly explore vast regions of the parameter space (e.g., varying the potential strength parameters like 'α' or 'm∗' in Section 4) to find combinations of these parameters that yield a viable relic abundance prediction.

  • The system can quantify the sensitivity of the relic abundance to specific physical effects, such as determining whether a model is more sensitive to resonance contributions (low temperature/high mass ratio) or off-shell contributions (large decay width).

Improved System Capability:

  • Automated discovery of viable dark matter models by efficiently navigating complex, multi-parameter spaces and identifying the specific physical mechanisms (e.g., bound state resonances vs. off-shell effects) that drive the relic abundance predictions for a given set of experimental constraints.

) 4. Consistency Checking and Model Validation (Theoretical Rigor Enforcement)

The paper explicitly discusses consistency conditions, such as requiring the annihilation rate of the bound state into DM to be much smaller than its total decay width (Equation 2.31).

Improvements for AI systems:

  • AI can be integrated into a theoretical workflow to perform automated consistency checks on proposed dark matter models. Before presenting a final relic abundance prediction, the AI can verify if the underlying assumptions (like the validity of the Born approximation or the consistency condition between decay width and annihilation rate) are met for all relevant bound states.

  • This capability prevents garbage in, garbage out by flagging theoretical inconsistencies arising from overly simplistic approximations.

Improved System Capability:

  • Acting as an automated theoretical validator that ensures complex physical models adhere to necessary consistency conditions (like the Born approximation validity), thereby increasing the scientific rigor and reliability of derived predictions.

Abstract

We study the Sommerfeld enhancement of the annihilation cross section of dark matter into heavier unstable particles. In this process, the annihilation products become non-relativistic near the kinematical threshold. If they experience long-range interactions with each other, their wave function is distorted from a plane wave, and the annihilation cross section can be significantly enhanced. When evaluating the Sommerfeld enhancement from the long-range interactions between the annihilation products, the decay of the products needs to be taken into account. We treat this issue by including the decay width in the Schrödinger equations of the two-body wave function of the annihilation products. We find that bound states of the annihilation products with a narrow decay width enhance the annihilation cross section through a resonant effect. At the same time, this formulation automatically includes the annihilation process with off-shell final state particles, which is relevant for a wide decay width. We show that the resonant effect significantly affects the prediction of the dark matter relic abundance.

Sources

Related papers