We revise the bound from the supernova SN1987A on the coupling of ultralight axion-like particles (ALPs) to photons. In a core-collapse supernova, ALPs would be emitted via the Primakoff process, and eventually convert into gamma rays in the magnetic field of the Milky Way. The lack of a gamma-ray signal in the GRS instrument of the SMM satellite in coincidence with the observation of the neutrinos emitted from SN1987A therefore provides a strong bound on their coupling to photons. Due to the large uncertainty associated with the current bound, we revise this argument, based on state-of-the-art physical inputs both for the supernova models and for the Milky-Way magnetic field. Furthermore, we provide major amendments, such as the consistent treatment of nucleon-degeneracy effects and of the reduction of the nuclear masses in the hot and dense nuclear medium of the supernova. With these improvements, we obtain a new upper limit on the photon-ALP coupling: g aγ 5.3 × 10 −12 GeV −1 , for m a 4.4 × 10 −10 eV , and we also give its dependence at larger ALP masses m a . Moreover, we discuss how much the Fermi-LAT satellite experiment could improve this bound, should a close-enough supernova explode in the near future.

We model the transport of cosmic ray nuclei in the Galaxy by means of a new numerical code. Differently from previous numerical models we account for a generic spatial distribution of the diffusion coefficient. We found that in the case of radially uniform diffusion, the main secondary/primary ratios (B/C, N/O and sub-Fe/Fe) and the modulated antiproton spectrum match consistently the available observations. Convection and re-acceleration do not seem to be required in the energy range we consider: 1 ≤ E ≤ 10 3 GeV/nucleon. We generalize these results accounting for radial dependence of the diffusion coefficient, which is assumed to trace that of the cosmic ray sources. While this does not affect the prediction of secondary/primary ratios, the simulated longitude profile of the diffuse γ-ray emission is significantly different from the uniform case and may agree with EGRET measurements without invoking ad hoc assumptions on the galactic gas density distribution. ‡ This unit of measure corresponds to the energy E competing on average to each nucleon within a nucleus with A nucleons and energy A × E.

Abstract. We present version 2 of the DRAGON code designed for computing realistic predictions of the CR densities in the Galaxy. The code numerically solves the interstellar CR transport equation (including inhomogeneous and anisotropic diffusion, either in space and momentum, advective transport and energy losses), under realistic conditions.The new version includes an updated numerical solver and several models for the astrophysical ingredients involved in the transport equation. Improvements in the accuracy of the numerical solution are proved against analytical solutions and in reference diffusion scenarios.The novel features implemented in the code allow to simulate the diverse scenarios proposed to reproduce the most recent measurements of local and diffuse CR fluxes, going beyond the limitations of the homogeneous galactic transport paradigm. To this end, several applications using DRAGON2 are presented as well.This new version facilitates the users to include their own physical models by means of a modular C++ structure.

The unprecedented quality of the data collected by the AMS-02 experiment onboard the International Space Station allowed us to address subtle questions concerning the origin and propagation of cosmic rays. Here we discuss the implications of these data for the injection spectrum of elements with different masses and for the diffusion coefficient probed by cosmic rays through their journey from the sources to the Earth. We find that the best fit to the spectra of primary and secondary nuclei requires (1) a break in the energy dependence of the diffusion coefficient at energies ∼ 300 GV; (2) an injection spectrum that is the same for all nuclei heavier than helium, and different injections for both protons and helium. Moreover, if to force the injection spectrum of helium to be the same as for heavier nuclei, the fit to oxygen substantially worsens. Accounting for a small, Xs ∼ 0.4 g cm −2 , grammage accumulated inside the sources leads to a somewhat better fit to the B/C ratio but makes the difference between He and other elements even more evident. The statistic and systematic error bars claimed by the AMS collaboration exceed the error that is expected from calculations once the uncertainties in the cross sections of production of secondary nuclei are taken into account. In order to make this point more quantitative, we present a novel parametrization of a large set of cross sections, relevant for cosmic ray physics, and we introduce the uncertainty in the branching ratios in a way that its effect can be easily grasped.