Production and propagation of secondary antideuteron in the Galaxy

Luis Fernando Galicia Cruztitla,Diego Mauricio Gomez Coral

2024-11-06

Abstract:This work reviews the current state of the antideuteron ($\bar{d}$) production cross-sections in cosmic ray interactions and its uncertainties, considering the coalescence model and measurements in accelerator experiments. These cross-sections have been included in a simulation of cosmic rays propagation in the Galaxy using GALPROP v.57, with updated parameters of the diffusive reacceleration model. An estimation of the expected antideuteron flux at Earth is presented.

High Energy Astrophysical Phenomena

What problem does this paper attempt to address?

The problem that this paper attempts to solve is about the production and propagation mechanisms of antideuterons (\(\bar{d}\)) in the Milky Way galaxy, especially the production cross - section of secondary antideuterons and its uncertainty. The research mainly focuses on the following aspects:

1. **Production of antideuterons**:
- The paper reviews the current situation of the cross - section of antideuterons produced in cosmic - ray interactions and its uncertainty, taking into account the coalescence model and the measurement results of accelerator experiments.
- Through the data generated by Monte Carlo (MC) simulations, the distributions of antiprotons (\(\bar{p}\)) and antineutrons (\(\bar{n}\)) in phase space were studied, and how they combine to form antideuterons under specific conditions was explored.

2. **Propagation of antideuterons**:
- The GALPROP v.57 software package was used to simulate the propagation process of cosmic rays in the Milky Way galaxy, taking into account the latest parameters of the diffusion - reacceleration model.
- Through the simulation, the antideuteron flux received by the Earth was estimated.

3. **Comparison between theory and experiment**:
- The antideuteron flux obtained from the simulation was compared with the experimental data of the Alpha Magnetic Spectrometer (AMS - 02) to verify the accuracy of the model.
- Special attention was paid to the 7 antideuteron candidate events reported by AMS - 02 in the energy region above 1 GeV, and whether these events may come from secondary sources or dark matter annihilation was explored.

### Main formulas

1. **Transport equation**:
\[
\frac{\partial \psi(\bar{r}, p, t)}{\partial t}=Q(\bar{r}, p)+\nabla \cdot [D_{xx}\nabla \psi-\bar{V}\psi]+\frac{\partial}{\partial p}\left(p^{2}D_{pp}\frac{\partial}{\partial p}\frac{\psi}{p^{2}}\right)-\frac{\partial}{\partial p}\left(\frac{dp}{dt}\psi-\frac{1}{3}p\nabla \cdot \bar{V}\psi\right)-\frac{\psi}{\tau_f}-\frac{\psi}{\tau_d}
\]
where \(\psi\) represents the particle density per unit momentum, \(Q(\bar{r}, p)\) represents the cosmic - ray source term, \(D_{xx}\) is the diffusion coefficient, \(\bar{V}\) is the convection velocity, \(D_{pp}\) is the diffusion coefficient in momentum space, and \(\tau_f\) and \(\tau_d\) represent the time scales of fragmentation loss and radioactive decay respectively.

2. **Antideuteron production spectrum**:
\[
\gamma_{\bar{d}}\frac{\partial^{3}N_{\bar{d}}}{\partial p^{3}_{\bar{d}}}=\frac{4\pi}{3}p_{0}^{3}\left(\gamma_{\bar{p}}\frac{\partial^{3}N_{\bar{p}}}{\partial p^{3}_{\bar{p}}}\right)\left(\gamma_{\bar{n}}\frac{\partial^{3}N_{\bar{n}}}{\partial p^{3}_{\bar{n}}}\right)
\]
where \(\gamma\) is the Lorentz factor of each particle and \(p_{0}\) is the coalescence momentum.

3. **Secondary source term**:
\[
Q_{\text{sec}}^{\bar{d}}(E_{\text{kin}}^{\bar{d}}, r)=\sum_{