Erratum: "Inference of the Local Interstellar Spectra of Cosmic-Ray Nuclei Z ⩽ 28 with the GALPROP-HELMOD Framework" (2020, ApJS, 250, 27).

[This corrects the article PMC8549769.].

The Discovery of a Low-energy Excess in Cosmic-Ray Iron: Evidence of the Past Supernova Activity in the Local Bubble.

Since its launch, the Alpha Magnetic Spectrometer-02 (AMS-02) has delivered outstanding quality measurements of the spectra of cosmic-ray (CR) species ( p ¯ , e ±, and nuclei, 1H-8O, 10Ne, 12Mg, 14Si) which resulted in a number of breakthroughs. One of the latest long-awaited surprises is the spectrum of 26Fe just published by AMS-02. Because of the large fragmentation cross section and large ionization energy losses, most of CR iron at low energies is local and may harbor some features associated with relatively recent supernova (SN) activity in the solar neighborhood. Our analysis of the new AMS-02 results, together with Voyager 1 and ACE-CRIS data, reveals an unexpected bump in the iron spectrum and in the Fe/He, Fe/O, and Fe/Si ratios at 1-2 GV, while a similar feature in the spectra of He, O, and Si and in their ratios is absent, hinting at a local source of low-energy CRs. The found excess extends the recent discoveries of radioactive 60Fe deposits in terrestrial and lunar samples and in CRs. We provide an updated local interstellar spectrum (LIS) of iron in the energy range from 1 MeV nucleon-1 to ~10 TeV nucleon-1. Our calculations employ the GALPROP-HELMOD framework, which has proved to be a reliable tool in deriving the LIS of CR p ¯ , e -, and nuclei Z ⩽ 28.

Deciphering the Local Interstellar Spectra of Secondary Nuclei with the Galprop/Helmod Framework and a Hint for Primary Lithium in Cosmic Rays.

Local interstellar spectra (LIS) of secondary cosmic-ray (CR) nuclei, lithium, beryllium, boron, and partially secondary nitrogen, are derived in the rigidity range from 10 MV to ~200 TV using the most recent experimental results combined with state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. The lithium spectrum appears somewhat flatter at high energies compared to other secondary species, which may imply a primary lithium component. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for the model-data comparison. The proposed LIS accommodates the low-energy interstellar spectra measured by Voyager 1, the High Energy Astrophysics Observatory-3 (HEAO-3), and the Cosmic Ray Isotope Spectrometer on board of the Advanced Composition Explorer (ACE/CRIS), as well as the high-energy observations by the Payload for Antimatter Matter Exploration and Light-nuclei Astrophysics (PAMELA), Alpha Magnetic Spectrometer-02 (AMS-02), and earlier experiments that are made deep in the heliosphere. The interstellar and heliospheric propagation parameters derived in this study are consistent with our earlier results for propagation of CR protons, helium, carbon, oxygen, antiprotons, and electrons.

Inference of the Local Interstellar Spectra of Cosmic-Ray Nuclei Z ⩽ 28 with the GalProp-HelMod Framework.

Composition and spectra of Galactic cosmic rays (CRs) are vital for studies of high-energy processes in a variety of environments and on different scales, for interpretation of γ-ray and microwave observations, for disentangling possible signatures of new phenomena, and for understanding of our local Galactic neighborhood. Since its launch, AMS-02 has delivered outstanding-quality measurements of the spectra of p ¯ , e ±, and nuclei: 1H-8O, 10Ne, 12Mg, 14Si. These measurements resulted in a number of breakthroughs; however, spectra of heavier nuclei and especially low-abundance nuclei are not expected until later in the mission. Meanwhile, a comparison of published AMS-02 results with earlier data from HEAO-3-C2 indicates that HEAO-3-C2 data may be affected by undocumented systematic errors. Utilizing such data to compensate for the lack of AMS-02 measurements could result in significant errors. In this paper we show that a fraction of HEAO-3-C2 data match available AMS-02 measurements quite well and can be used together with Voyager 1 and ACE-CRIS data to make predictions for the local interstellar spectra (LIS) of nuclei that are not yet released by AMS-02. We are also updating our already-published LIS to provide a complete set from 1H-28Ni in the energy range from 1 MeV nucleon-1 to ~100-500 TeV nucleon-1, thus covering 8-9 orders of magnitude in energy. Our calculations employ the GalProp-HelMod framework, which has proved to be a reliable tool in deriving the LIS of CR p ¯ , e -, and nuclei 1H-8O.

Cosmic-ray antinuclei as messengers of new physics: status and outlook for the new decade.

The precise measurement of cosmic-ray antinuclei serves as an important means for identifying the nature of dark matter and other new astrophysical phenomena, and could be used with other cosmic-ray species to understand cosmic-ray production and propagation in the Galaxy. For instance, low-energy antideuterons would provide a "smoking gun" signature of dark matter annihilation or decay, essentially free of astrophysical background. Studies in recent years have emphasized that models for cosmic-ray antideuterons must be considered together with the abundant cosmic antiprotons and any potential observation of antihelium. Therefore, a second dedicated Antideuteron Workshop was organized at UCLA in March 2019, bringing together a community of theorists and experimentalists to review the status of current observations of cosmic-ray antinuclei, the theoretical work towards understanding these signatures, and the potential of upcoming measurements to illuminate ongoing controversies. This review aims to synthesize this recent work and present implications for the upcoming decade of antinuclei observations and searches. This includes discussion of a possible dark matter signature in the AMS-02 antiproton spectrum, the most recent limits from BESS Polar-II on the cosmic antideuteron flux, and reports of candidate antihelium events by AMS-02; recent collider and cosmic-ray measurements relevant for antinuclei production models; the state of cosmic-ray transport models in light of AMS-02 and Voyager data; and the prospects for upcoming experiments, such as GAPS. This provides a roadmap for progress on cosmic antinuclei signatures of dark matter in the coming years.

Deciphering Residual Emissions: Time-dependent Models for the Nonthermal Interstellar Radiation from the Milky Way.

Cosmic rays (CRs) in the Galaxy are an important dynamical component of the interstellar medium (ISM) that interact with the other major components (interstellar gas and magnetic and radiation fields) to produce broadband interstellar emissions that span the electromagnetic spectrum. The standard modeling of CR propagation and production of the associated emissions is based on a steady-state assumption, where the CR source spatial density is described using a smoothly varying function of position that does not evolve with time. While this is a convenient approximation, reality is otherwise, where primary CRs are produced in and about highly localized regions, e.g., supernova remnants, which have finite lifetimes. In this paper, we use the latest version of the galprop CR propagation code to model time-dependent CR injection and propagation through the ISM from a realistic 3D discretized CR source density distribution, together with full 3D models for the other major ISM components, and make predictions of the associated broadband nonthermal emissions. We compare the predictions for the discretized and equivalent steady-state model, finding that the former predicts novel features in the broadband nonthermal emissions that are absent for the steady-state case. Some of the features predicted by the discretized model may be observable in all-sky observations made by WMAP and Planck, the recently launched eROSITA, the Fermi-LAT, and ground-based observations by HESS, HAWC, and the forthcoming CTA. The nonthermal emissions predicted by the discretized model may also provide explanations of puzzling anomalies in high-energy γ-ray data, such as the Fermi-LAT north/south asymmetry and residuals like the so-called "Fermi bubbles."

AAfrag: Interpolation routines for Monte Carlo results on secondary production in proton-proton, proton-nucleus and nucleus-nucleus interactions.

We provide a compilation of predictions of the QGSJET-II-04m model for the production of secondary species (photons, neutrinos, electrons, positrons, and antinucleons) that are covering a wide range of energies of the beam particles in proton-proton, proton-nucleus, nucleus-proton, and nucleus-nucleus reactions. The current version of QGSJET-II-04m has an improved treatment of the production of secondary particles at low energies: the parameters of the hadronization procedure have been fine-tuned, based on a number of recent benchmark experimental data, notably, from the LHCf, LHCb, and NA61 experiments. Our results for the production spectra are made publicly accessible through the interpolation routines AAfrag which are described below. Besides, we comment on the impact of Feynman scaling violation and isospin symmetry effects on antinucleon production.

AMS-100: The Next Generation Magnetic Spectrometer in Space - An International Science Platform for Physics and Astrophysics at Lagrange Point 2.

The next generation magnetic spectrometer in space, AMS-100, is designed to have a geometrical acceptance of 100 m2 sr and to be operated for at least ten years at the Sun-Earth Lagrange Point 2. Compared to existing experiments, it will improve the sensitivity for the observation of new phenomena in cosmic rays, and in particular in cosmic antimatter, by at least a factor of 1000. The magnet design is based on high temperature superconductor tapes, which allow the construction of a thin solenoid with a homogeneous magnetic field of 1 Tesla inside. The inner volume is instrumented with a silicon tracker reaching a maximum detectable rigidity of 100 TV and a calorimeter system that is 70 radiation lengths deep, equivalent to four nuclear interaction lengths, which extends the energy reach for cosmic-ray nuclei up to the PeV scale, i.e. beyond the cosmic-ray knee. Covering most of the sky continuously, AMS-100 will detect high-energy gamma rays in the calorimeter system and by pair conversion in the thin solenoid, reconstructed with excellent angular resolution in the silicon tracker.

Galactic PeVatrons and helping to find them: Effects of galactic absorption on the observed spectra of very high energy γ-ray sources.

Identification of the cosmic-ray (CR) "PeVatrons," which are sources capable of accelerating particles to ~1015 eV energies and higher, may lead to resolving the long-standing question of the origin of the spectral feature in the all-particle CR spectrum known as the "knee." Because CRs with these energies are deflected by interstellar magnetic fields identification of individual sources and determination of their spectral characteristics is more likely via very high energy γ-ray emissions, which provide the necessary directional information. However, pair production on the interstellar radiation field (ISRF) and cosmic microwave background (CMB) leads to steepening of the high energy tails of γ-ray spectra, and should be corrected for to enable true properties of the spectrum at the source to be recovered. Employing recently developed three-dimensional ISRF models this paper quantifies the pair-absorption effect on spectra for sources in the Galactic center (GC) direction at 8.5 and 23.5 kpc distances, with the latter corresponding to the far side of the Galactic stellar disc where it is expected that discrimination of spectral features >10 TeV is possible by the forthcoming Cherenkov Telescope Array (CTA). The estimates made suggest spectral cutoffs could be underestimated by factors of a few in the energy range so far sampled by TeV γ-ray telescopes. As an example to illustrate this, the recent HESS measurements of diffuse γ-ray emissions possibly associated with injection of CRs nearby Sgr A* are ISRF corrected, and estimates of the spectral cutoff are reevaluated. It is found that it could be higher by up to a factor of ~2, indicating that these emissions may be consistent with a CR accelerator with a spectral cutoff of at least 1 PeV at the 95% confidence level.

HelMod in the Works: From Direct Observations to the Local Interstellar Spectrum of Cosmic-Ray Electrons.

The local interstellar spectrum (LIS) of cosmic-ray (CR) electrons for the energy range 1 MeV to 1 TeV is derived using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The optimized HelMod parameters are then used to adjust GALPROP parameters to predict a refined LIS with the procedure repeated subject to a convergence criterion. The parameter optimization uses an extensive data set of proton spectra from 1997 to 2015. The proposed CR electron LIS accommodates both the low-energy interstellar spectra measured by Voyager 1 as well as the high-energy observations by PAMELA and AMS-02 that are made deep in the heliosphere; it also accounts for Ulysses counting rate features measured out of the ecliptic plane. The interstellar and heliospheric propagation parameters derived in this study agree well with our earlier results for CR protons, helium nuclei, and anti-protons propagation and LIS obtained in the same framework.

Deciphering the Local Interstellar Spectra of Primary Cosmic-Ray Species with helmod.

Local interstellar spectra (LIS) of primary cosmic ray (CR) nuclei, such as helium, oxygen, and mostly primary carbon are derived for the rigidity range from 10 MV to ~200 TV using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. Two propagation packages, GALPROP and helmod, are combined into a single framework that is used to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. The developed iterative maximum-likelihood method uses GALPROP-predicted LIS as input to helmod, which provides the modulated spectra for specific time periods of the selected experiments for model-data comparison. The interstellar and heliospheric propagation parameters derived in this study are consistent with our prior analyses using the same methodology for propagation of CR protons, helium, antiprotons, and electrons. The resulting LIS accommodate a variety of measurements made in the local interstellar space (Voyager 1) and deep inside the heliosphere at low (ACE/CRIS, HEAO-3) and high energies (PAMELA, AMS-02).

High-energy Gamma Rays from the Milky Way: Three-dimensional Spatial Models for the Cosmic-Ray and Radiation Field Densities in the Interstellar Medium.

High-energy γ-rays of interstellar origin are produced by the interaction of cosmic-ray (CR) particles with the diffuse gas and radiation fields in the Galaxy. The main features of this emission are well understood and are reproduced by existing CR propagation models employing 2D galactocentric cylindrically symmetrical geometry. However, the high-quality data from instruments like the Fermi Large Area Telescope reveal significant deviations from the model predictions on few to tens of degrees scales, indicating the need to include the details of the Galactic spiral structure and thus requiring 3D spatial modeling. In this paper, the high-energy interstellar emissions from the Galaxy are calculated using the new release of the GALPROP code employing 3D spatial models for the CR source and interstellar radiation field (ISRF) densities. Three models for the spatial distribution of CR sources are used that are differentiated by their relative proportion of input luminosity attributed to the smooth disk or spiral arms. Two ISRF models are developed based on stellar and dust spatial density distributions taken from the literature that reproduce local near- to far-infrared observations. The interstellar emission models that include arms and bulges for the CR source and ISRF densities provide plausible physical interpretations for features found in the residual maps from high-energy γ-ray data analysis. The 3D models for CR and ISRF densities provide a more realistic basis that can be used for the interpretation of the nonthermal interstellar emissions from the Galaxy.

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.

FERMI-LAT OBSERVATIONS OF HIGH- AND INTERMEDIATE-VELOCITY CLOUDS: TRACING COSMIC RAYS IN THE HALO OF THE MILKY WAY.

It is widely accepted that cosmic rays (CRs) up to at least PeV energies are Galactic in origin. Accelerated particles are injected into the interstellar medium where they propagate to the farthest reaches of the Milky Way, including a surrounding halo. The composition of CRs coming to the solar system can be measured directly and has been used to infer the details of CR propagation that are extrapolated to the whole Galaxy. In contrast, indirect methods, such as observations of γ-ray emission from CR interactions with interstellar gas, have been employed to directly probe the CR densities in distant locations throughout the Galactic plane. In this article we use 73 months of data from the Fermi Large Area Telescope in the energy range between 300 MeV and 10 GeV to search for γ-ray emission produced by CR interactions in several high- and intermediate-velocity clouds (IVCs) located at up to ~7 kpc above the Galactic plane. We achieve the first detection of IVCs in γ rays and set upper limits on the emission from the remaining targets, thereby tracing the distribution of CR nuclei in the halo for the first time. We find that the γ-ray emissivity per H atom decreases with increasing distance from the plane at 97.5% confidence level. This corroborates the notion that CRs at the relevant energies originate in the Galactic disk. The emissivity of the upper intermediate-velocity Arch hints at a 50% decline of CR densities within 2 kpc from the plane. We compare our results to predictions of CR propagation models.

ANALYTIC SOLUTION FOR SELF-REGULATED COLLECTIVE ESCAPE OF COSMIC RAYS FROM THEIR ACCELERATION SITES.

Supernova remnants (SNRs), as the major contributors to the galactic cosmic rays (CRs), are believed to maintain an average CR spectrum by diffusive shock acceleration regardless of the way they release CRs into the interstellar medium (ISM). However, the interaction of the CRs with nearby gas clouds crucially depends on the release mechanism. We call into question two aspects of a popular paradigm of the CR injection into the ISM, according to which they passively and isotropically diffuse in the prescribed magnetic fluctuations as test particles. First, we treat the escaping CR and the Alfvén waves excited by them on an equal footing. Second, we adopt field-aligned CR escape outside the source, where the waves become weak. An exact analytic self-similar solution for a CR "cloud" released by a dimmed accelerator strongly deviates from the test-particle result. The normalized CR partial pressure may be approximated as P ( p , z , t ) = 2 [ | z | 5 / 3 + z dif 5 / 3 ( p , t ) ] - 3 / 5  exp [ - z 2 / 4 D ISM ( p ) t ] , where p is the momentum of CR particle, and z is directed along the field. The core of the cloud expands as z dif ∝ D NL ( p ) t and decays in time as P ∝ 2 z dif - 1 ( t ) . The diffusion coefficient D NL is strongly suppressed compared to its background ISM value D ISM: D NL ~ D ISM exp (-Π) ≪ D ISM for sufficiently high field-line-integrated CR partial pressure, Π. When Π â‰« 1, the CRs drive Alfvén waves efficiently enough to build a transport barrier ( P ≈ 2 / | z | - "  pedestal " ) that strongly reduces the leakage. The solution has a spectral break at p = p br, where p br satisfies the equation D NL ( p br ) ≃ z 2 / t .

Development of laser induced fluorescence diagnostic for measuring the parameters of plasma containing rare gas species.

Laser induced fluorescence (LIF) technique development activity for measurement of plasma parameters in ITER divertor plasma is described. Helium density is the task of priority, but Doppler measurement of ion (atom) temperatures is also the aim of the program. The concept of ITER scenarios includes injection of "extrinsic" impurities (Ne, Ar, and Kr). It is possible to use the species as tracing elements for measurement of T(i), T(a). The program included modeling experiments on PNX-U (a multicusp trap with microwave argon plasma). Helium was added by puffing into discharge. Temperatures T(i)(Ar(1+)) and T(a)(He(0)) have been measured by scanning laser line across absorption line of species. Summarizing of fluorescence signals provided input data for estimation of Ar(1+) and He(0) densities via interpretative collisional-radiative models. Besides, the collisional-radiative model has been used for estimation of electron density using the ratio of fluorescence signals at 388.9 and 706.5 nm helium lines.