Diffuse Cosmic Rays Shining in the Galactic Center: A Novel Interpretation of H.E.S.S. and Fermi-LAT γ-Ray Data

We present a novel interpretation of the γ-ray diffuse emission measured by Fermi-LAT and H.E.S.S. in the Galactic center (GC) region and the Galactic ridge (GR). In the first part we perform a data-driven analysis based on PASS8 Fermi-LAT data: We extend down to a few GeV the spectra measured by H.E.S.S. and infer the primary cosmic-ray (CR) radial distribution between 0.1 and 3 TeV. In the second part we adopt a CR transport model based on a position-dependent diffusion coefficient. Such behavior reproduces the radial dependence of the CR spectral index recently inferred from the Fermi-LAT observations. We find that the bulk of the GR emission can be naturally explained by the interaction of the diffuse steady-state Galactic CR sea with the gas present in the central molecular zone. Although we confirm the presence of a residual radial-dependent emission associated with a central source, the relevance of the large-scale diffuse component prevents to claim a solid evidence of GC pevatrons.

Possible Dark Matter Annihilation Signal in the AMS-02 Antiproton Data

Using the latest AMS-02 cosmic-ray antiproton flux data, we search for a potential dark matter annihilation signal. The background parameters about the propagation, source injection, and solar modulation are not assumed a priori but based on the results inferred from the recent B/C ratio and proton data measurements instead. The possible dark matter signal is incorporated into the model self-consistently under a Bayesian framework. Compared with the astrophysical background-only hypothesis, we find that a dark matter signal is favored. The rest mass of the dark matter particles is ∼20-80  GeV, and the velocity-averaged hadronic annihilation cross section is about (0.2-5)×10^{-26}  cm^{3} s^{-1}, in agreement with that needed to account for the Galactic center GeV excess and/or the weak GeV emission from dwarf spheroidal galaxies Reticulum 2 and Tucana III. Tight constraints on the dark matter annihilation models are also set in a wide mass region.

Solution of Heliospheric Propagation: Unveiling the Local Interstellar Spectra of Cosmic-ray Species

Local interstellar spectra (LIS) for protons, helium, and antiprotons are built using the most recent experimental results combined with state-of-the-art models for propagation in the Galaxy and heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of cosmic-ray (CR) species at different modulation levels and at both polarities of the solar magnetic field. To do so in a self-consistent way, an iterative procedure was developed, where the GALPROP LIS output is fed into HelMod, providing modulated spectra for specific time periods of selected experiments to compare with the data; the HelMod parameter optimization is performed at this stage and looped back to adjust the LIS using the new GALPROP run. The parameters were tuned with the maximum likelihood procedure using an extensive data set of proton spectra from 1997 to 2015. The proposed LIS accommodate both the low-energy interstellar CR spectra measured by Voyager 1 and the high-energy observations by BESS, Pamela, AMS-01, and AMS-02 made from the balloons and near-Earth payloads; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The found solution is in a good agreement with proton, helium, and antiproton data by AMS-02, BESS, and PAMELA in the whole energy range.

GALACTIC COSMIC RAYS IN THE LOCAL INTERSTELLAR MEDIUM: VOYAGER 1 OBSERVATIONS AND MODEL RESULTS

Since 2012 August Voyager 1 has been observing the local interstellar energy spectra of Galactic cosmic-ray nuclei down to 3 MeV nuc-1 and electrons down to 2.7 MeV. The H and He spectra have the same energy dependence between 3 and 346 MeV nuc-1, with a broad maximum in the 10-50 MeV nuc-1 range and a H/He ratio of 12.2 ± 0.9. The peak H intensity is ~15 times that observed at 1 AU, and the observed local interstellar gradient of 3-346 MeV H is -0.009 ± 0.055% AU-1, consistent with models having no local interstellar gradient. The energy spectrum of electrons (e - + e +) with 2.7-74 MeV is consistent with E -1.30±0.05 and exceeds the H intensity at energies below ~50 MeV. Propagation model fits to the observed spectra indicate that the energy density of cosmic-ray nuclei with >3 MeV nuc-1 and electrons with >3 MeV is 0.83-1.02 eV cm-3 and the ionization rate of atomic H is in the range of 1.51-1.64 × 10-17 s-1. This rate is a factor >10 lower than the ionization rate in diffuse interstellar clouds, suggesting significant spatial inhomogeneity in low-energy cosmic rays or the presence of a suprathermal tail on the energy spectrum at much lower energies. The propagation model fits also provide improved estimates of the elemental abundances in the source of Galactic cosmic rays.

BAYESIAN ANALYSIS OF COSMIC RAY PROPAGATION: EVIDENCE AGAINST HOMOGENEOUS DIFFUSION

We present the results of the most complete scan of the parameter space for cosmic ray (CR) injection and propagation. We perform a Bayesian search of the main GALPROP parameters, using the MultiNest nested sampling algorithm, augmented by the BAMBI neural network machine-learning package. This is the first study to separate out low-mass isotopes (p, p ¯ , and He) from the usual light elements (Be, B, C, N, and O). We find that the propagation parameters that best-fit p, p ¯ , and He data are significantly different from those that fit light elements, including the B/C and 10Be/9Be secondary-to-primary ratios normally used to calibrate propagation parameters. This suggests that each set of species is probing a very different interstellar medium, and that the standard approach of calibrating propagation parameters using B/C can lead to incorrect results. We present posterior distributions and best-fit parameters for propagation of both sets of nuclei, as well as for the injection abundances of elements from H to Si. The input GALDEF files with these new parameters will be included in an upcoming public GALPROP update.

Measurements of cosmic-ray proton and helium spectra from the BESS-Polar long-duration balloon flights over Antarctica

The BESS-Polar Collaboration measured the energy spectra of cosmic-ray protons and helium during two long-duration balloon flights over Antarctica in December 2004 and December 2007, at substantially different levels of solar modulation. Proton and helium spectra probe the origin and propagation history of cosmic rays in the galaxy, and are essential to calculations of the expected spectra of cosmic-ray antiprotons, positrons, and electrons from interactions of primary cosmic-ray nuclei with the interstellar gas, and to calculations of atmospheric muons and neutrinos. We report absolute spectra at the top of the atmosphere for cosmic-ray protons in the kinetic energy range 0.2-160 GeV and helium nuclei 0.15-80 GeV/nucleon. The corresponding magnetic rigidity ranges are 0.6-160 GV for protons and 1.1-160 GV for helium. These spectra are compared to measurements from previous BESS flights and from ATIC-2, PAMELA, and AMS-02. We also report the ratio of the proton and helium fluxes from 1.1 GV to 160 GV and compare to ratios from PAMELA and AMS-02.

Radiation Measurements Performed with Active Detectors Relevant for Human Space Exploration

A reliable radiation risk assessment in space is a mandatory step for the development of countermeasures and long-duration mission planning in human spaceflight. Research in radiobiology provides information about possible risks linked to radiation. In addition, for a meaningful risk evaluation, the radiation exposure has to be assessed to a sufficient level of accuracy. Consequently, both the radiation models predicting the risks and the measurements used to validate such models must have an equivalent precision. Corresponding measurements can be performed both with passive and active devices. The former is easier to handle, cheaper, lighter, and smaller but they measure neither the time dependence of the radiation environment nor some of the details useful for a comprehensive radiation risk assessment. Active detectors provide most of these details and have been extensively used in the International Space Station. To easily access such an amount of data, a single point access is becoming essential. This review presents an ongoing work on the development of a tool that allows obtaining information about all relevant measurements performed with active detectors providing reliable inputs for radiation model validation.

Precision Measurement of the Helium Flux in Primary Cosmic Rays of Rigidities 1.9 GV to 3 TV with the Alpha Magnetic Spectrometer on the International Space Station

Knowledge of the precise rigidity dependence of the helium flux is important in understanding the origin, acceleration, and propagation of cosmic rays. A precise measurement of the helium flux in primary cosmic rays with rigidity (momentum/charge) from 1.9 GV to 3 TV based on 50 million events is presented and compared to the proton flux. The detailed variation with rigidity of the helium flux spectral index is presented for the first time. The spectral index progressively hardens at rigidities larger than 100 GV. The rigidity dependence of the helium flux spectral index is similar to that of the proton spectral index though the magnitudes are different. Remarkably, the spectral index of the proton to helium flux ratio increases with rigidity up to 45 GV and then becomes constant; the flux ratio above 45 GV is well described by a single power law.

Signatures of a Two Million Year Old Supernova in the Spectra of Cosmic Ray Protons, Antiprotons, and Positrons

The locally observed cosmic ray spectrum has several puzzling features, such as the excess of positrons and antiprotons above ~20  GeV and the discrepancy in the slopes of the spectra of cosmic ray protons and heavier nuclei in the TeV-PeV energy range. We show that these features are consistently explained by a nearby source which was active approximately two million years ago and has injected (2-3)×10^{50} erg in cosmic rays. The transient nature of the source and its overall energy budget point to the supernova origin of this local cosmic ray source. The age of the supernova suggests that the local cosmic ray injection was produced by the same supernova that has deposited ^{60}Fe isotopes in the deep ocean crust.

Precision Measurement of the Proton Flux in Primary Cosmic Rays from Rigidity 1 GV to 1.8 TV with the Alpha Magnetic Spectrometer on the International Space Station

A precise measurement of the proton flux in primary cosmic rays with rigidity (momentum/charge) from 1 GV to 1.8 TV is presented based on 300 million events. Knowledge of the rigidity dependence of the proton flux is important in understanding the origin, acceleration, and propagation of cosmic rays. We present the detailed variation with rigidity of the flux spectral index for the first time. The spectral index progressively hardens at high rigidities.

NEW CALCULATION OF ANTIPROTON PRODUCTION BY COSMIC RAY PROTONS AND NUCLEI

A dramatic increase in the accuracy and statistics of space-borne cosmic ray (CR) measurements has yielded several breakthroughs over the last several years. The most puzzling is the rise in the positron fraction above ~10 GeV over the predictions of the propagation models assuming pure secondary production. The accuracy of the antiproton production cross section is critical for astrophysical applications and searches for new physics since antiprotons in CRs seem to hold the keys to many puzzles including the origin of those excess positrons. However, model calculations of antiproton production in CR interactions with interstellar gas are often employing parameterizations that are out of date or are using outdated physical concepts. This may lead to an incorrect interpretation of antiproton data which could have broad consequences for other areas of astrophysics. In this work, we calculate antiproton production in pp-, pA-, and AA-interactions using EPOS-LHC and QGSJET-II-04, two of the most advanced Monte Carlo (MC) generators tuned to numerous accelerator data including those from the Large Hadron Collider (LHC). We show that the antiproton yields obtained with these MC generators differ by up to an order of magnitude from yields of parameterizations commonly used in astrophysics.

Low level irradiation in mice can lead to enhanced trabecular bone morphology

Charged particle radiation such as iron ions and their secondary fragmentation products are of particular concern to the skeleton due to their high charge and energy deposition. However, little is known about the long-term effects of these particles on trabecular and cortical bone morphology when applied at relatively low levels. We hypothesized that even a 4.4 cGy dose of a complex secondary iron ion radiation field will compromise skeletal quantity and architecture in adult mice. One year after radiation exposure and compared to age-matched controls, 4.4 cGy irradiated mice had 51 % more trabecular bone, 56 % greater trabecular bone volume fraction, 16 % greater trabecular number, and 17 % less trabecular separation in the distal metaphysis of the femur. Similar to the metaphysis, trabecular bone of the distal femoral epiphysis in 4.4 cGy mice had 33 % more trabecular bone, 31 % greater trabecular bone volume fraction, and a 33 % smaller structural model index. Cortical bone morphology, whole bone mechanical properties, and lower leg muscle mass were unaffected. When compared to two additional groups, irradiated at either 8.9 or 17.8 cGy, a (negative) dose response relationship was observed for trabecular bone in the metaphysis but not in the epiphysis. In contrast to our original hypothesis, these data indicated that a secondary field of low-level, high-linear energy transfer iron radiation may cause long-term augmentation, rather than deterioration, of trabecular bone in the femoral metaphysis and epiphysis of mice.

NUCLEAR ENHANCEMENT OF THE PHOTON YIELD IN COSMIC RAY INTERACTIONS

The concept of the nuclear enhancement factor has been used since the beginning of γ-ray astronomy. It provides a simple and convenient way to account for the contribution of nuclei (A > 1) in cosmic rays (CRs) and in the interstellar medium (ISM) to the diffuse γ-ray emission. An accurate treatment of the dominant emission process, such as hadronic interactions of CRs with the ISM, enables one to study CR acceleration processes and CR propagation in the ISM, and provides a reliable background model for searches of new phenomena. The Fermi Large Area Telescope launched in 2008 provides excellent quality data in a wide energy range 30 MeV-1 TeV where the diffuse emission accounts for the majority of photons. Exploiting its data to the fullest requires a new study of the processes of γ-ray production in hadronic interactions. In this paper we point out that several commonly used studies of the nuclear enhancement factor fail to account for the spectrally averaged energy loss fraction which ensures that the energy fraction transferred to photons is averaged properly with the spectra of CR species. We present a new calculation of the spectrally averaged energy loss fraction and the nuclear enhancement factor using the QGSJET-II-04 and EPOS-LHC interaction models.

Inferred cosmic-ray spectrum from Fermi large area telescope γ-ray observations of Earth's limb

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly discovered features may offer a clue to the origin of high-energy CRs. We use the Fermi Large Area Telescope observations of the γ-ray emission from Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range ∼90  GeV-6  TeV (derived from a photon energy range 15 GeV-1 TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index 2.68±0.04 and 2.61±0.08 above ∼200  GeV, respectively.

AMS-02 results support the secondary origin of cosmic ray positrons

We show that the recent AMS-02 positron fraction measurement is consistent with a secondary origin for positrons and does not require additional primary sources such as pulsars or dark matter. The measured positron fraction at high energy saturates the previously predicted upper bound for secondary production, obtained by neglecting radiative losses. This coincidence, which will be further tested by upcoming AMS-02 data at higher energy, is a compelling indication for a secondary source. Within the secondary model, the AMS-02 data imply a cosmic ray propagation time in the Galaxy of <10(6) yr and an average traversed interstellar matter density of ~1 cm(-3), comparable to the density of the Milky Way gaseous disk, at a rigidity of 300 GV.

Cosmic Rays in the Milky Way and Beyond

Cosmic rays (CRs) are the probes of the deep space. They allow us to study particle acceleration, chemical composition of the interstellar medium, and global properties of our Galaxy. However, until recently studies of CRs were similar to astronomical observations with blurred lenses that capture only the integral brightness of all stars in the field of view. Thanks to the recent technological developments, our "lenses" are now capable of capturing sharp images and making precise measurements of all CR species. We have a full range of unique instrumentation for direct measurements of CRs in space and for multi-wavelength observations of their emissions and more missions are coming. The unveiling picture is astonishing. This paper gives a short overview of very exciting developments in astrophysics of CRs in the Milky Way and in other star-forming galaxies.

Three-dimensional model of cosmic-ray lepton propagation reproduces data from the Alpha Magnetic Spectrometer on the International Space Station

We study the compatibility of Alpha Magnetic Spectrometer (AMS-02) data on the cosmic-ray (CR) positron fraction with data on the CR electron and positron spectra provided by PAMELA and Fermi LAT. We do that in terms of a novel propagation model in which sources are distributed in spiral arm patterns in agreement with astrophysical observations. While former interpretations assumed an unrealistically steep injection spectrum for astrophysical background electrons, the enhanced energy losses experienced by CR leptons due to the larger average source distance from Earth allow us to reproduce the data with harder injection spectra as expected in a shock acceleration scenario. Moreover, we show that in this approach, and accounting for AMS-02 results, the contribution of nearby accelerators to the fluxes at very high energy can be significantly reduced, thus avoiding any tension with anisotropy upper limits.

The long-term variability of cosmic ray protons in the heliosphere: A modeling approach

Galactic cosmic rays are charged particles created in our galaxy and beyond. They propagate through interstellar space to eventually reach the heliosphere and Earth. Their transport in the heliosphere is subjected to four modulation processes: diffusion, convection, adiabatic energy changes and particle drifts. Time-dependent changes, caused by solar activity which varies from minimum to maximum every ∼11 years, are reflected in cosmic ray observations at and near Earth and along spacecraft trajectories. Using a time-dependent compound numerical model, the time variation of cosmic ray protons in the heliosphere is studied. It is shown that the modeling approach is successful and can be used to study long-term modulation cycles.

Diffuse Galactic gamma rays from shock-accelerated cosmic rays

A shock-accelerated particle flux is proportional to p(-s), where p is the particle momentum, follows from simple theoretical considerations of cosmic-ray acceleration at nonrelativistic shocks followed by rigidity-dependent escape into the Galactic halo. A flux of shock-accelerated cosmic-ray protons with s≈2.8 provides an adequate fit to the Fermi Large Area Telescope γ-ray emission spectra of high-latitude and molecular cloud gas when uncertainties in nuclear production models are considered. A break in the spectrum of cosmic-ray protons claimed by Neronov, Semikoz, and Taylor [Phys. Rev. Lett. 108, 051105 (2012)] when fitting the γ-ray spectra of high-latitude molecular clouds is a consequence of using a cosmic-ray proton flux described by a power law in kinetic energy.