Research

Paper

TESTING March 19, 2026

Toward universal coalescence models for antideuteron production

Authors

Mattia di Mauro, Jordan Koechler, Lorenzo Stefanuto, Francesca Bellini, Fiorenza Donato, Nicolao Fornengo

Abstract

Cosmic-ray (CR) antinuclei, especially antideuteron $\overline{\rm D}$ and antihelium-3 nuclei ${}^3\overline{\rm He}$, are among the most promising messengers for indirect dark matter (DM) searches. This is because secondary production in CR interactions with the interstellar medium is strongly suppressed at kinetic energies $K\simeq (0.1 - 1)$ GeV/$n$, typically one to two orders of magnitude below fluxes expected in standard DM scenarios. From the theoretical side, the formation of $\overline{\rm D}$ and ${}^3\overline{\rm He}$ is governed by coalescence, whose dynamics cannot yet be reliably derived from first principles. Phenomenological approaches therefore introduce effective coalescence parameters, possibly dependent on collision energy and production environment (hadronic versus electroweak). In this work we show, for the first time, that a common set of physically motivated coalescence models can simultaneously reproduce collider data in two qualitatively different regimes: ALICE measurements of (anti)deuteron production in $pp$ collisions at $\sqrt{s}=(0.9 - 13)$ TeV and the ALEPH $\overline{\rm D}$ multiplicity in hadronic $Z$ decays at $\sqrt{s}=m_Z$. We test both simple event-by-event prescriptions based on a relative-momentum cutoff, finding a preferred coalescence scale $p_{\rm coal}\simeq 0.2$ GeV, and quantum-mechanical models in the Wigner formalism. In the latter, a Gaussian bound-state wavefunction gives a best-fit momentum width, corresponding to $δ\simeq 1.7$ fm, while a parameter-free implementation using the Argonne $v_{18}$ wavefunction (constrained by proton-neutron scattering data) agrees with ALICE spectra at the $\sim 25\%$ level. Overall, our results support an approximately universal coalescence description across energies and production environments, strengthening the theoretical basis for interpreting upcoming CR antinuclei searches.

Metadata

arXiv ID: 2603.19352
Provider: ARXIV
Primary Category: hep-ph
Published: 2026-03-19
Fetched: 2026-03-23 16:54

Related papers

Raw Data (Debug)
{
  "raw_xml": "<entry>\n    <id>http://arxiv.org/abs/2603.19352v1</id>\n    <title>Toward universal coalescence models for antideuteron production</title>\n    <updated>2026-03-19T18:00:00Z</updated>\n    <link href='https://arxiv.org/abs/2603.19352v1' rel='alternate' type='text/html'/>\n    <link href='https://arxiv.org/pdf/2603.19352v1' rel='related' title='pdf' type='application/pdf'/>\n    <summary>Cosmic-ray (CR) antinuclei, especially antideuteron $\\overline{\\rm D}$ and antihelium-3 nuclei ${}^3\\overline{\\rm He}$, are among the most promising messengers for indirect dark matter (DM) searches. This is because secondary production in CR interactions with the interstellar medium is strongly suppressed at kinetic energies $K\\simeq (0.1 - 1)$ GeV/$n$, typically one to two orders of magnitude below fluxes expected in standard DM scenarios. From the theoretical side, the formation of $\\overline{\\rm D}$ and ${}^3\\overline{\\rm He}$ is governed by coalescence, whose dynamics cannot yet be reliably derived from first principles. Phenomenological approaches therefore introduce effective coalescence parameters, possibly dependent on collision energy and production environment (hadronic versus electroweak). In this work we show, for the first time, that a common set of physically motivated coalescence models can simultaneously reproduce collider data in two qualitatively different regimes: ALICE measurements of (anti)deuteron production in $pp$ collisions at $\\sqrt{s}=(0.9 - 13)$ TeV and the ALEPH $\\overline{\\rm D}$ multiplicity in hadronic $Z$ decays at $\\sqrt{s}=m_Z$. We test both simple event-by-event prescriptions based on a relative-momentum cutoff, finding a preferred coalescence scale $p_{\\rm coal}\\simeq 0.2$ GeV, and quantum-mechanical models in the Wigner formalism. In the latter, a Gaussian bound-state wavefunction gives a best-fit momentum width, corresponding to $δ\\simeq 1.7$ fm, while a parameter-free implementation using the Argonne $v_{18}$ wavefunction (constrained by proton-neutron scattering data) agrees with ALICE spectra at the $\\sim 25\\%$ level. Overall, our results support an approximately universal coalescence description across energies and production environments, strengthening the theoretical basis for interpreting upcoming CR antinuclei searches.</summary>\n    <category scheme='http://arxiv.org/schemas/atom' term='hep-ph'/>\n    <category scheme='http://arxiv.org/schemas/atom' term='astro-ph.HE'/>\n    <published>2026-03-19T18:00:00Z</published>\n    <arxiv:comment>18 pages, 11 figures, 6 tables. The differential cross-section for secondary antiprotons and antideuterons we computed in this work are publicly accessible at the following Zenodo repository: https://doi.org/10.5281/zenodo.19099608</arxiv:comment>\n    <arxiv:primary_category term='hep-ph'/>\n    <author>\n      <name>Mattia di Mauro</name>\n    </author>\n    <author>\n      <name>Jordan Koechler</name>\n    </author>\n    <author>\n      <name>Lorenzo Stefanuto</name>\n    </author>\n    <author>\n      <name>Francesca Bellini</name>\n    </author>\n    <author>\n      <name>Fiorenza Donato</name>\n    </author>\n    <author>\n      <name>Nicolao Fornengo</name>\n    </author>\n  </entry>"
}