1 Introduction

# The curious case of high-energy deuterons in Galactic cosmic rays

## Abstract

A new analysis of cosmic ray (CR) data collected by the SOKOL experiment in space found that the deuteron-to-helium ratio at energies between 500 and 2000 GeV/nucleon takes the value d/He 1.5. As we will show, this result cannot be explained by standard models of secondary CR production in the interstellar medium and points to the existence of a high-energy source of CR deuterons. To account for the deuteron excess in CRs, we argue that the only viable solution is hadronic interaction processes of accelerated particles inside old supernova remnants (SNRs). From this mechanism, however, the B/C ratio is also expected to increase at energy above  50 of GeV/nucleon, in conflict with new precision data just released by the AMS-02 experiment. Hence, if this phenomenon is a real physical effect, hadronic production of CR deuterons must occur in SNRs characterized by low metal abundance. In such a scenario, the sources accelerating C-N-O nuclei are not the same as those accelerating helium or protons, so that the connection between d/He ratio and B/C ratio is broken, and the latter cannot be used to place constraints on the production of light isotopes or antiparticles.

cosmic rays — acceleration of particles — ISM: supernova remnants

Università degli Studi di Perugia & INFN-Perugia, I-06100 Perugia, Italy; \addressSchool of Physics, Sun Yat-Sen University, Guangzhou 510275, China; \addressInstitute of Physics, Academia Sinica, Nankang, Taipei 11529, Taiwan;

## 1. Introduction

Deuteron isotopes H are rare particles in Galactic cosmic rays (CRs). They are destroyed rather than formed in thermonuclear reactions in stellar interior so that, from supernova remnants (SNRs) as sources of CRs, no significant amount of deuteron is expected to be released by diffusive shock acceleration (DSA). An important process by which high-energy deuterons can be created in the Galaxy is nuclear fragmentation of CR nuclei with the gas of the interstellar medium (ISM). The main source of deuteron production is fragmentation of He and He isotopes, along with C-N-O nuclei, in collisions with interstellar hydrogen and helium. An important contribution comes from the reaction p+pd+. From nuclear fragmentation, interactions of CRs with the ISM are also known to generate Li-Be-B nuclei and antiparticles, that are otherwise rare from stellar nucleosynthesis processes. The measured abundances of these elements in the cosmic radiation enables us to pose tight constraints on the astrophysical models of CR propagation in the Galaxy and, in the context of indirect searches of dark matter, to asses the astrophysical background of CRs antiparticles (Grenier et al., 2015).

Very recently, the Alpha Magnetic Spectrometer (AMS-02) experiment has reported a new precise measurement of the B/C ratio in CRs at energy between  0.5 and 1000 GeV/n (Aguilar et al., 2016). At kinetic energies above  10 GeV/n, the ratio is found to decrease steadily with increasing energy. Since boron nuclei are mainly produced by C-N-O fragmentation in the ISM, the observed trend of the B/C ratio reflects the conception that CRs diffusively propagate in the Galactic magnetic fields with an average diffusion coefficient that increases with energy. According to this picture, the ratio between H and its main progenitor He must follow the same behavior in the GeV-TeV energy range. Recent calculations have shown, indeed, that the observed boron and deuteron abundances at GeV/n energies can both be self-consistently described by diffusion models of CR propagation (Coste et al., 2012; Tomassetti, 2012a). In the GeV-TeV energy region, the d/He ratio is expected to decrease rapidly, as fast as the B/C ratio does, but CR deuterons have never been detected at these energies.

Quite unexpectedly, a new analysis of the data collected by the satellite mission SOKOL has determined the d/He ratio at 0.5-2 TeV/n energy (Turundaevskiy & Podorozhnyi, 2016). In this measurement, the discrimination between H and He isotopes was performed by means of neural network analyses of the topology of hadronic showers developed by these particles. By means of two different Monte Carlo simulations, consistent results have been obtained: 0.1140.023 and 0.0990.021 for the H/(HH) ratio, 1.640.30 and 1.430.27 for the H/He ratio.

In this Letter, we show that the above results represent a striking anomaly in CR physics that cannot be explained by standard models of CR propagation and secondary production in the ISM. Then, we argue that a deuteron excess can be explained in terms of hadronic production occurring inside SNRs, which was proposed in Blasi (2009) to explain the positron excess in CRs (see also Kachelrießet al., 2011; Serpico, 2012). Using new evaluations of fragmentation cross-sections, we calculate for the first time the high-energy production of CR deuterons in SNRs and in the ISM, demonstrating that this mechanism can account for the new d/He data. Along with the d/He ratio, we also compute the B/C ratio under the same framework, showing that there are conflicting results in the model predictions for the two observables. We therefore conclude that this tension can be resolved if the deuterons progenitors are not accelerated in the same sources of boron progenitors. We discuss our results and their implications for the interpretation of antiproton data.

## 2. Calculations

We compute the spectrum of CR nuclei accelerated in SNRs within the linear DSA theory and including the production of secondary fragments. Similar calculations are done in earlier works (Blasi, 2009; Blasi & Serpico, 2009; Mertsch & Sarkar, 2009, 2014; Tomassetti & Donato, 2015; Herms et al., 2016). We follow closely the derivation of Tomassetti & Donato (2012). In the shock rest-frame (), the upstream plasma flows in from with speed (density ) and the downstream plasma flows out to with speed (density ). The compression ratio is . For a nucleus with charge and mass number , the equation describing diffusion and convection at the shock reads

 u∂f∂x=D∂2f∂x2+13dudxp∂f∂p−Γtotf+Q, (1)

where is the phase space density, is the diffusion coefficient near the shock, is the fluid speed, and is the total fragmentation rate for cross-sections and background density , which is assumed to be composed of H and He like the average ISM. The source term includes particle injection at the shock, , with and GV for all nuclei. The -constants set the normalization of each species. We assume strong shocks () and a diffusion coefficient , where parameterizes the deviation of from the Bohm value due to magnetic damping. The resulting acceleration rate at momentum is . For an SNR of age , the condition defines the maximum momentum scale attainable by DSA. In the presence of hadronic interactions, the additional requirement must be fulfilled. The downstream solution reads

 f2(x,p)=f0(p)(1−Γtot2u2x)+q2u2x, (2)

where is the distribution function at the shock. As found in Mertsch & Sarkar (2009), is given by

 f0(p)=α∫p0(p′p)αYδ(p′−pinj)e−χ(p,p′)dp′p′ (3) +α∫p0(p′p)αq1Du21(1+r2)e−χ(p,p′)dp′p′,

with , , and the subscript () indicates the upstream (downstream) region. The first term of Eq. 3 describes primary particles injected at the shock and it is of the form . The second term of Eq. 3 describes the production and acceleration of CR fragments from heavier progenitors. For each process, the –term of Eq. 3 is given by , where is the secondary production rate, and is the corresponding cross-section. The momentum inelasticity factor expresses the conservation of kinetic energy per nucleon between progenitor and fragment.

Equation 1 is solved for all relevant species. We considered primary nuclei (with ) p, He, C, N, O tuned to recent data as in Mertsch & Sarkar (2014), and secondary B production from C-N-O collisions with hydrogen and helium gas. The adopted fragmentation cross-sections are those re-evaluated in Tomassetti (2015b). We account for deuteron production from collisions of p, He, and He off hydrogen and helium. Measurements and calculations for these reactions are available only below 10 GeV/n of energy. Thus, we have performed new calculations at 10-1000 GeV/n using the hadronic Monte Carlo generator QGSJET-II-04 (Ostapchenko, 2011). The cross-sections for the dominant channels are shown in Fig. 1. The solid lines are our parameterizations adapted from Tomassetti (2012a) and extended to higher energy. Cross-sections for collisions off He target are larger by nearly a factor of two and contribute to 20 %. The contribution from fragmentation of heavier nuclei is negligible. In all reactions, kinetic energy per nucleon is approximately conserved except for the pp reaction, which peaks at a proton energy of  600 MeV/n and produces deuterons within a narrow energy range around 150 MeV/n (Meyer, 1972). Hence, these deuterons experience DSA acceleration similarly to primary components with non-zero -factor, i.e., to a power-law spectrum .

In contrast, deuterons ejected by helium fragmentation have an initial spectrum so that, for those produced within a distance from the shock and undergo DSA, their final spectrum reaches the form (see Eq. 3). This means that, for Bohm-type diffusion, the secondary deuteron spectrum is one power harder than that of its progenitors. The total CR flux produced by SNRs is computed as

 Ssnr(p)=4πp2RSN∫τsnru204πx2f2(x,p)dx (4)

where  25 Myr kpc is the explosion rate per unit volume and  40 kyr is the age of the SNR.

To model the subsequent propagation of CRs in the ISM, we adopt a two-halo model of CR diffusion and nuclear interactions (Tomassetti, 2012b, 2015a). The Galaxy is modeled as a disk of half-thickness  100 pc containing SNRs and gas with number density  1 cm. The disk is surrounded by a diffusive halo of half-thickness and zero matter density. We give a one-dimensional description in the thin disk limit (). For each CR nucleus, the transport equation reads

 ∂N∂t=∂∂z(K(z)∂N∂z)−2h^δ(z)~ΓtotN+2h^δ(z)Stot, (5)

where is its density, is a Dirac function of the -coordinate, is the diffusion coefficient of CRs in the Galaxy, and is the destruction rate in the ISM at velocity and cross-section . The source term is split into a primary term , obtained from Eq. 4 as solution of the DSA equation, and a secondary production term , from fragmentation of -type nuclei in the ISM with rate . To compute the interaction rates in the ISM we adopt the same fragmentation network (and same cross-sections) as occurring inside SNRs. Equation 5 is solved in steady-state conditions . The derivation of the full solution is in Tomassetti (2012b). The diffusion coefficient is taken of the type where expresses its normalization. For the scaling index, , we adopt in the region of (inner halo) and for (outer halo). Our default parameters are set as  1/3,  0.55, and 0.01 kpc Myr. The differential energy fluxes of each species are given by . Solar modulation is described in force-field approximation using the parameter  MV for a medium-level modulation strength. The d/He and B/C ratios as function of kinetic energy per nucleon are eventually calculated as and , respectively. In the following, we consider two model implementations representing two alternative scenarios:

• Scenario # 1 (B/C-driven, conservative) — standard model without interactions in sources, which is the case of a CR flux released by young SNRs with amplified magnetic fields ( 100 ) and/or low background density ( cm). In this model, secondary production of CRs deuterons or Li-Be-B occurs only in the ISM. This model is tuned to match the new B/C data from AMS-02 at GeV/n-TeV/n energies.

• Scenario # 2 (d/He-driven, speculative) — model with copious production and acceleration of secondary particles in SNR shockwaves, which is the case for a GeV-TeV flux provided by old SNRs with damped magnetic fields, slow shock speed, or dense ambient medium, i.e., with the combination  400. This model is tuned against the d/He data including the new SOKOL measurement at TeV/n energies.

## 3. Results and Discussion

Model calculations are shown in Fig. 2 for the B/C ratio the d/He ratio at energies between  0.5 GeV and 2 TeV per nucleon. Scenario # 1 is plotted as green solid lines. In this model, secondary nuclei such as d or B are entirely generated in the ISM, i.e., without SNR components, and thus secondary/primary ratios decrease steadily as where . In the high-energy limit one has . It can be seen that this model fits remarkably well the new AMS-02 data on the B/C ratio and dictates a similar trend for the d/He ratio, which is predicted to reach the level of in the TeV/n energy scale. We therefore conclude that the SOKOL measurement of the d/He ratio is at least two orders of magnitude higher than that expected from standard models where CR deuterons are produced by fragmentation in the ISM.

In Scenario # 2, shown as red dashed lines, hadronic interaction processes inside SNRs generate a source component of secondary nuclei, which is harder than that arising from CR collisions with the ISM and, as discussed, even harder than that of primary p-He spectra. It is then possible, with fragmentation inside SNRs, having secondary/primary ratios that increase with energy. Figure 2 shows that Scenario # 2 matches fairly well the d/He ratio measurements at GeV/n and at TeV/n energies, therefore providing an explanation for the new SOKOL data. Under this model, however, the B/C ratio is also predicted to increase, at energies above  50 GeV/n, in remarkable contrast with the new AMS-02 data. While interactions inside SNRs seem to be the only mechanism capable of explaining a rise in the d/He ratio, it is apparent that the observed decrease of the B/C ratio conflicts with this mechanism. We also note that, at the  1 GeV/n energy region where secondary CR production in the ISM dominates, the two ratios are consistently described by both models # 1 and # 2, at least within the precision of the current data.

As we see it, the only solution to this tension is a situation where the connection between d/He ratio and B/C ratio is broken. This situation is realized if the sources accelerating helium and protons (and producing deuterons) are not the same as those accelerating heavier C-N-O nuclei and, in particular, deuterons must be accelerated by a low-metallicity source, which may be the case for SNRs expanding over H-dominated molecular clouds. Such a possibility was also discussed in Cholis & Hooper (2014) (see Sect. VI), and proposed in other works (Fujita et al., 2009; Kohri et al., 2016), all focused on antiparticle excesses in CRs. In such a scenario, the B/C ratio can no longer be used to place constraints on antimatter spectra. In this respect, it is important to note that the connection between d/He ratio and antimatter/matter ratios would still be preserved because, in contrast to Li-Be-B nuclei, secondary deuterons share their progenitors with positrons and antiprotons.

It is then interesting to calculate the antiproton/proton ratio in light of new data released by AMS-02 (Aguilar et al., 2016). Calculations were performed in Blasi & Serpico (2009) and subsequent works (e.g., Mertsch & Sarkar, 2014). In this work, the antiproton distribution at the shock is calculated numerically, as done in Herms et al. (2016), in order to drop the “inelasticity approximation” that links the antiproton momentum to the primary proton momentum through an assumedly constant factor . It was noted that such an approximation leads to an overestimate of the high-energy antiproton production in SNRs (Kachelrießet al., 2011). Along with antiproton production from p-He collisions with hydrogen and helium gas, we also account for tertiary reactions (such as ) and for destruction processes. All these processes are implemented in both acceleration and propagation. The corresponding cross-sections are taken from Feng et al. (2016). The resulting / ratio is shown in Fig. 3. It can be seen that Scenario #2 is preferred by the AMS-02 data. Reducing nuclear uncertainties in antiproton production is clearly essential for a complete discrimination between the different models (Feng et al., 2016).

## 4. Experimental Challenges

Given the implications of these new data on the phenomenology of CR propagation, we believe that the situation deserves more clarification on the experimental side. The SOKOL analysis relies on unconventional techniques of deuteron/proton mass separation, which is always a very challenging task. For instance the d/He measurement might be overestimated due to undetected background arising, e.g., from the mass distribution tails of CR protons or from He nuclei fragmenting in the top of the instrument. On the other hand, it is very unlikely for such a background to affect the results by two orders of magnitude. The deuteron spectrum and the d/He ratio are being precisely measured by the AMS-02 experiment at  0.1-10 GeV/n with standard spectrometric techniques. In addition, AMS-02 is also equipped with a Transition Radiation Detector (TRD), designed for lepton/hadron mass separation, which can provide direct measurements of the Lorentz factor at TeV/n energies (Obermeier & Korsmeyer, 2015). In standard magnetic spectrometers, the CR mass is derived from velocity and momentum measurements, , so that its corresponding resolution is given by

 (δMM)2=(δpp)2+γ4(δββ)2, (6)

showing that the mass resolution degradates rapidly, at relativistic energies, due to the -factor. In contrast, with the opportunity of performing direct TRD-based -measurements, AMS-02 may have the capability to detect CR deuterons at the energy scale.

## 5. Conclusions

This work is aimed at interpreting new data, registered by the SOKOL experiment in space, that revealed a surprisingly high abundance of CR deuterons in the TeV/n region. In contrast to antiparticle excesses that can be explained, e.g., by pulsar models or dark matter annihilation (Serpico, 2012; Feng & Zhang, 2016), the SOKOL data demand an enhanced high-energy production of CR deuterons from hadronic interactions. We found that no explanation for this measurement can be provided in terms of standard collisions of CRs with the gas of the ISM. As we have shown, a viable solution for this puzzle is the occurrence of nuclear fragmentation inside SNRs, but this mechanism conflicts with the new AMS-02 data on the B/C ratio.

Thus, if the SOKOL measurement is taken as face value, we conclude that the sources accelerating helium and protons (thereby producing deuteron) may not be the same as those accelerating C-N-O nuclei (otherwise producing Li-Be-B nuclei), and that the former are more efficient in the production and acceleration of secondary particles. Under such a scenario, the connection between B/C ratio and antiparticle/particle ratios would also be broken, and thus the B/C ratio should not be used to place constraints on the astrophysical antimatter background. On the other hand, the d/He ratio would still represent a direct diagnostic tool for assessing this background.
We are grateful to Tanguy Pierog, Colin Baus, and Ralf Ulrich for the CRMC interface to MC generators. N.T. acknowledges support from the H2020-MSCA-IF-2015 action, grant No. 707543-MAtISSE. J.F. acknowledges support from China Scholarship Council and the Taiwanese Ministry of Science and Technology under grants No. 104-2112-M-001-024 and No. 105-2112-M-001-003.

### Footnotes

1. slugcomment: December 2016

### References

1. Adriani, O., et al., 2010, PRL, 105, 121101
2. Aguilar, M., et al., 2016, PRL, 117, 231102
3. Aguilar, M., et al., 2016, PRL, 117, 091103
4. Aguilar, M., et al., 2011, ApJ, 736, 105
5. Blasi, P., 2009, PRL, 103, 051104
6. Blasi, P., & Serpico, P. D., 2009, PRL, 103, 081103
7. Cholis, I., & Hooper, D., 2014, PRD, 89, 043013
8. Coste, B., et al., 2012, A&A, 539, A88
9. de Nolfo, G., et al., 2000, AIP, 528, 425
10. Fujita, Y., et al., 2009, PRD, 80, 063003
11. Feng, J., Tomassetti, N., & Oliva, A., 2016, PRD, 94, 123007
12. Feng, J., & Zhang, H. H., 2016, Eur.Phys.J. C 76, 229
13. Grenier, I. A., Black, J. H., & Strong, A. W., 2015, Annu. Rev. Astron. Astrophys., 53, 199-246
14. Herms, J., et al., 2016, arXiv:1610.00699
15. Kachelrieß, M., Ostapchenko, & S., Tomàs, R., 2011, ApJ, 733, 11
16. Kim, K. C., et al., 2013, ASR, 51, 234
17. Kohri, K., et al., 2016, Prog. Theor. Exp. Phys., 021E01
18. Mertsch, P., & Sarkar, S., 2009, PRL, 103, 081104
19. Mertsch, P., & Sarkar, S., 2014, PRD, 90, 061301
20. Meyer, J. P., 1972, A&A Suppl., 7, 417–467
21. Myers, Z. D., et al., 2005, ASR, 35, 151
22. Obermeier, A., & Korsmeyer, M., 2015, ASR, 55, 716
23. Ostapchenko, S., 2011, PRD, 83, 014018
24. Papini, R., et al., 2004, ApJ, 615, 259
25. Serpico, P., 2012, Astropart. Phy. 39-40, 2-11
26. Tomassetti, N., 2012, Astrophys. Space Sci., 342, 131
27. Tomassetti, N., 2012, ApJ, 752, L13
28. Tomassetti, N., 2015, PRD, 92, 081301(R)
29. Tomassetti, N., 2015, PRC, 92, 045808
30. Tomassetti, N., & Donato, F., 2012, A&A, 544, A16
31. Tomassetti, N., & Donato, F., 2015, ApJ, 803, L15
32. Turundaevskiy, A., & Podorozhnyi, D., 2017, ASR, 59, 496
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minumum 40 characters