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

Cancer and circulatory disease risks for the largest solar particle events in the space age

In this paper we use the NASA Space Cancer Risk (NSCR version 2022) model to predict cancer and circulatory disease risks using energy spectra representing the largest SPE's observed in the space age. Because tissue dose-rates behind shielding for large SPE's lead to low dose-rates (<0.2 Gy/h) we consider the integrated risk for several historical periods of high solar activity, including July-November, 1960 events and August-October 1989 events along with the February 1956 and August 1972 events. The galactic cosmic ray (GCR) contribution to risks is considered in predictions. Results for these largest historical events show risk of exposure induced death (REID) are mitigated to < 1.2 % with a 95 % confidence interval with passive radiation shielding of 20 g/cm2 aluminum, while larger amounts would support the application of the ALARA principle. Annual GCR risks are predicted to surpass the risks from large SPEs by ∼30 g/cm2 of aluminum shielding.

Role of Low-Energy (<20 eV) Secondary Electrons in the Extraterrestrial Synthesis of Prebiotic Molecules

We demonstrate for the first time that Galactic cosmic rays with energies as high as ∼1010 eV can trigger a cascade of low-energy (<20 eV) secondary electrons that could be a significant contributor to the interstellar synthesis of prebiotic molecules whose delivery by comets, meteorites, and interplanetary dust particles may have kick-started life on Earth. For the energetic processing of interstellar ice mantles inside dark, dense molecular clouds, we explore the relative importance of low-energy (<20 eV) secondary electrons-agents of radiation chemistry-and low-energy (<10 eV), nonionizing photons-instigators of photochemistry. Our calculations indicate fluxes of ∼102 electrons cm-2 s-1 for low-energy secondary electrons produced within interstellar ices due to attenuated Galactic cosmic-ray protons. Consequently, in certain star-forming regions where internal high-energy radiation sources produce ionization rates that are observed to be a thousand times greater than the typical interstellar Galactic ionization rate, the flux of low-energy secondary electrons should far exceed that of nonionizing photons. Because reaction cross sections can be several orders of magnitude larger for electrons than for photons, even in the absence of such enhancement, our calculations indicate that secondary low-energy (<20 eV) electrons are at least as significant as low-energy (<10 eV) nonionizing photons in the interstellar synthesis of prebiotic molecules. Most importantly, our results demonstrate the pressing need for explicitly incorporating low-energy electrons in current and future astrochemical simulations of cosmic ices. Such models are critically important for interpreting James Webb Space Telescope infrared measurements, which are currently being used to probe the origins of life by studying complex organic molecules found in ices near star-forming regions.

Cosmic ray feedback in galaxies and galaxy clusters: A pedagogical introduction and a topical review of the acceleration, transport, observables, and dynamical impact of cosmic rays

Understanding the physical mechanisms that control galaxy formation is a fundamental challenge in contemporary astrophysics. Recent advances in the field of astrophysical feedback strongly suggest that cosmic rays (CRs) may be crucially important for our understanding of cosmological galaxy formation and evolution. The appealing features of CRs are their relatively long cooling times and relatively strong dynamical coupling to the gas. In galaxies, CRs can be close to equipartition with the thermal, magnetic, and turbulent energy density in the interstellar medium, and can be dynamically very important in driving large-scale galactic winds. Similarly, CRs may provide a significant contribution to the pressure in the circumgalactic medium. In galaxy clusters, CRs may play a key role in addressing the classic cooling flow problem by facilitating efficient heating of the intracluster medium and preventing excessive star formation. Overall, the underlying physics of CR interactions with plasmas exhibit broad parallels across the entire range of scales characteristic of the interstellar, circumgalactic, and intracluster media. Here we present a review of the state-of-the-art of this field and provide a pedagogical introduction to cosmic ray plasma physics, including the physics of wave-particle interactions, acceleration processes, CR spatial and spectral transport, and important cooling processes. The field is ripe for discovery and will remain the subject of intense theoretical, computational, and observational research over the next decade with profound implications for the interpretation of the observations of stellar and supermassive black hole feedback spanning the entire width of the electromagnetic spectrum and multi-messenger data.

Effect of He2+ ion irradiation on the mechanical properties of automated fibre placement (AFP) CF-PEEK thermoplastics composites

Carbon fibre-reinforced polyetheretherketone (CF-PEEK) composites have gained significant usage across diverse industries like automotive and aerospace due to their desirable characteristics. These properties encompass recyclability, low density, high strength, wear resistance and thermal stability. The components made from CF-PEEK composites for space applications will be subjected to a high radiation environment due to the incoming cosmic rays, comprising protons, α particles, electrons, γ rays, etc., once they escape the Earth's atmosphere. The ion irradiation of CF-PEEK is accompanied by radiation-induced effects, which drastically change the structure and properties of irradiated material. Since the resistance of CF-PEEK to radiation damage has not been studied extensively, this study aims to understand the effect of high-energy He2+ ions on the microstructure and properties of CF-PEEK composites manufactured using automated fibre placement (AFP) under different processing conditions. The samples have been radiated with 5 MeV He2+ ions using an energy degrader wheel to create a layer with relatively uniform damage. Then, were characterized using optical and scanning electron microscopy and their hardness was evaluated using nanoindentation. It was observed that, irradiation increases the hardness of the fibres in all cases. Also, fibre orientation affects the hardness in a statistically significant manner in both unirradiated and irradiated conditions.

Properties of Cosmic-Ray Sulfur and Determination of the Composition of Primary Cosmic-Ray Carbon, Neon, Magnesium, and Sulfur: Ten-Year Results from the Alpha Magnetic Spectrometer

We report the properties of primary cosmic-ray sulfur (S) in the rigidity range 2.15 GV to 3.0 TV based on 0.38×10^{6} sulfur nuclei collected by the Alpha Magnetic Spectrometer experiment (AMS). We observed that above 90 GV the rigidity dependence of the S flux is identical to the rigidity dependence of Ne-Mg-Si fluxes, which is different from the rigidity dependence of the He-C-O-Fe fluxes. We found that, similar to N, Na, and Al cosmic rays, over the entire rigidity range, the traditional primary cosmic rays S, Ne, Mg, and C all have sizeable secondary components, and the S, Ne, and Mg fluxes are well described by the weighted sum of the primary silicon flux and the secondary fluorine flux, and the C flux is well described by the weighted sum of the primary oxygen flux and the secondary boron flux. The primary and secondary contributions of the traditional primary cosmic-ray fluxes of C, Ne, Mg, and S (even Z elements) are distinctly different from the primary and secondary contributions of the N, Na, and Al (odd Z elements) fluxes. The abundance ratio at the source for S/Si is 0.167±0.006, for Ne/Si is 0.833±0.025, for Mg/Si is 0.994±0.029, and for C/O is 0.836±0.025. These values are determined independent of cosmic-ray propagation.

The cosmic ray ionization and γ-ray budgets of star-forming galaxies

Cosmic rays in star-forming galaxies are a dominant source of both diffuse γ-ray emission and ionization in gas too deeply shielded for photons to penetrate. Though the cosmic rays responsible for γ-rays and ionization are of different energies, they are produced by the same star formation-driven sources, and thus galaxies' star formation rates, γ-ray luminosities, and ionization rates should all be linked. In this paper, we use up-to-date cross-section data to determine this relationship, finding that cosmic rays in a galaxy of star formation rate [Formula: see text] and gas depletion time t dep produce a maximum primary ionization rate ζ ≈ 1 × 10-16(t dep/Gyr)-1 s-1 and a maximum γ-ray luminosity [Formula: see text] erg s-1 in the 0.1-100 GeV band. These budgets imply either that the ionization rates measured in Milky Way molecular clouds include a significant contribution from local sources that elevate them above the Galactic mean, or that CR-driven ionization in the Milky Way is enhanced by sources not linked directly to star formation. Our results also imply that ionization rates in starburst systems are only moderately enhanced compared to those in the Milky Way. Finally, we point out that measurements of γ-ray luminosities can be used to place constraints on galactic ionization budgets in starburst galaxies that are nearly free of systematic uncertainties on the details of cosmic ray acceleration.

Data driven analysis of cosmic rays in the heliosphere: diffusion of cosmic protons

Understanding the time-dependent relationship between the Sun’s variability and cosmic rays (GCR) is essential for developing predictive models of energetic radiation in space. When traveling inside the heliosphere, GCRs are affected by magnetic turbulence and solar wind disturbances which result in the so-called solar modulation effect. To investigate this phenomenon, we have performed a data-driven analysis of the temporal dependence of the GCR flux over the solar cycle. With a global statistical inference of GCR data collected in space by AMS-02 and PAMELA on monthly basis, we have determined the rigidity and time dependence of the GCR diffusion mean free path. Here we present our results for GCR protons, we discuss their interpretation in terms of basic processes of particle transport and their relations with the dynamics of the heliospheric plasma.

Future Exploration of the Outer Heliosphere and Very Local Interstellar Medium by Interstellar Probe

A detailed overview of the knowledge gaps in our understanding of the heliospheric interaction with the largely unexplored Very Local Interstellar Medium (VLISM) are provided along with predictions of with the scientific discoveries that await. The new measurements required to make progress in this expanding frontier of space physics are discussed and include in-situ plasma and pick-up ion measurements throughout the heliosheath, direct sampling of the VLISM properties such as elemental and isotopic composition, densities, flows, and temperatures of neutral gas, dust and plasma, and remote energetic neutral atom (ENA) and Lyman-alpha (LYA) imaging from vantage points that can uniquely discern the heliospheric shape and bring new information on the interaction with interstellar hydrogen. The implementation of a pragmatic Interstellar Probe mission with a nominal design life to reach 375 Astronomical Units (au) with likely operation out to 550 au are reported as a result of a 4-year NASA funded mission study.

Search for Cosmic-Ray Boosted Sub-GeV Dark Matter Using Recoil Protons at Super-Kamiokande

We report a search for cosmic-ray boosted dark matter with protons using the 0.37  megaton×years data collected at Super-Kamiokande experiment during the 1996-2018 period (SKI-IV phase). We searched for an excess of proton recoils above the atmospheric neutrino background from the vicinity of the Galactic Center. No such excess is observed, and limits are calculated for two reference models of dark matter with either a constant interaction cross section or through a scalar mediator. This is the first experimental search for boosted dark matter with hadrons using directional information. The results present the most stringent limits on cosmic-ray boosted dark matter and exclude the dark matter-nucleon elastic scattering cross section between 10^{-33}cm^{2} and 10^{-27}cm^{2} for dark matter mass from 1  MeV/c^{2} to 300  MeV/c^{2}.

Explaining the GeV Antiproton Excess, GeV γ-Ray Excess, and W-Boson Mass Anomaly in an Inert Two Higgs Doublet Model

For the newly discovered W-boson mass anomaly, one of the simplest dark matter (DM) models that can account for the anomaly without violating other astrophysical and experimental constraints is the inert two Higgs doublet model, in which the DM mass (m_{S}) is found to be within ∼54-74  GeV. In this model, the annihilation of DM via SS→bb[over ¯] and SS→WW^{*} would produce antiprotons and gamma rays, and may account for the excesses identified previously in both particles. Motivated by this, we reanalyze the AMS-02 antiproton and Fermi-LAT Galactic center γ-ray data. For the antiproton analysis, the novel treatment is the inclusion of the charge-sign-dependent three-dimensional solar modulation model as constrained by the time-dependent proton data. We find that the excess of antiprotons is more distinct than previous results based on the force-field solar modulation model. The interpretation of this excess as the annihilation of SS→WW^{*} (SS→bb[over ¯]) requires a DM mass of ∼40-80 (40-60) GeV and a velocity-averaged cross section of O(10^{-26})  cm^{3} s^{-1}. As for the γ-ray data analysis, besides adopting the widely used spatial template fitting, we employ an orthogonal approach with a data-driven spectral template analysis. The fitting to the GeV γ-ray excess yields DM model parameters overlapped with those to fit the antiproton excess via the WW^{*} channel. The consistency of the DM particle properties required to account for the W-boson mass anomaly, the GeV antiproton excess, and the GeV γ-ray excess suggests a common origin of them.

Galactic Cosmic Rays Throughout the Heliosphere and in the Very Local Interstellar Medium

We review recent observations and modeling developments on the subject of galactic cosmic rays through the heliosphere and in the Very Local Interstellar Medium, emphasizing knowledge that has accumulated over the past decade. We begin by highlighting key measurements of cosmic-ray spectra by Voyager, PAMELA, and AMS and discuss advances in global models of solar modulation. Next, we survey recent works related to large-scale, long-term spatial and temporal variations of cosmic rays in different regimes of the solar wind. Then we highlight new discoveries from beyond the heliopause and link these to the short-term evolution of transients caused by solar activity. Lastly, we visit new results that yield interesting insights from a broader astrophysical perspective.

The Structure of the Large-Scale Heliosphere as Seen by Current Models

This review summarizes the current state of research aiming at a description of the global heliosphere using both analytical and numerical modeling efforts, particularly in view of the overall plasma/neutral flow and magnetic field structure, and its relation to energetic neutral atoms. Being part of a larger volume on current heliospheric research, it also lays out a number of key concepts and describes several classic, though still relevant early works on the topic. Regarding numerical simulations, emphasis is put on magnetohydrodynamic (MHD), multi-fluid, kinetic-MHD, and hybrid modeling frameworks. Finally, open issues relating to the physical relevance of so-called "croissant" models of the heliosphere, as well as the general (dis)agreement of model predictions with observations are highlighted and critically discussed.

Inhomogeneity in the Local ISM and Its Relation to the Heliosphere

This paper reviews past research and new studies underway of the local interstellar environment and its changing influence on the heliosphere. The size, shape, and physical properties of the heliosphere outside of the heliopause are determined by the surrounding environment - now the outer region of the Local Interstellar Cloud (LIC). The temperature, turbulence, and velocity vector of neutral atoms and ions in the LIC and other partially ionized interstellar clouds are measured from high-resolution spectra of interstellar absorption lines observed with the STIS instrument on the HST. Analysis of such spectra led to a kinematic model with many interstellar clouds defined by velocity vectors derived from radial velocity measurements. This analysis identified fifteen clouds located within about 10 pc of the Sun and their mean temperatures, turbulence, and velocity vectors. With the increasing number of sight lines now being analyzed, we find that temperatures and turbulent velocities have spatial variations within the LIC and other nearby clouds much larger than measurement uncertainties, and that these spatial variations appear to be randomly distributed and can be fit by Gaussians. The inhomogeneous length scale is less than 4,000 AU, a distance that the heliosphere will traverse in less than 600 years. The temperatures and turbulent velocities do not show significant trends with stellar distance or angle from the LIC center. If/when the Sun enters an inter-cloud medium, the physical properties of the future heliosphere will be very different from the present. For the heliosheath and the very local interstellar medium (VLISM) just outside of the heliopause, the total pressures are approximately equal to the gravitational pressure of overlying material in the Galaxy. The internal pressure in the LIC is far below that in the VLISM, but there is an uncertain ram pressure term produced by the flow of the LIC with respect to its environment.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11214-022-00884-5.

Direct Measurement of the Nickel Spectrum in Cosmic Rays in the Energy Range from 8.8 GeV/n to 240 GeV/n with CALET on the International Space Station

The relative abundance of cosmic ray nickel nuclei with respect to iron is by far larger than for all other transiron elements; therefore it provides a favorable opportunity for a low background measurement of its spectrum. Since nickel, as well as iron, is one of the most stable nuclei, the nickel energy spectrum and its relative abundance with respect to iron provide important information to estimate the abundances at the cosmic ray source and to model the Galactic propagation of heavy nuclei. However, only a few direct measurements of cosmic-ray nickel at energy larger than ∼3  GeV/n are available at present in the literature, and they are affected by strong limitations in both energy reach and statistics. In this Letter, we present a measurement of the differential energy spectrum of nickel in the energy range from 8.8  to 240  GeV/n, carried out with unprecedented precision by the Calorimetric Electron Telescope (CALET) in operation on the International Space Station since 2015. The CALET instrument can identify individual nuclear species via a measurement of their electric charge with a dynamic range extending far beyond iron (up to atomic number Z=40). The particle's energy is measured by a homogeneous calorimeter (1.2 proton interaction lengths, 27 radiation lengths) preceded by a thin imaging section (3 radiation lengths) providing tracking and energy sampling. This Letter follows our previous measurement of the iron spectrum [1O. Adriani et al. (CALET Collaboration), Phys. Rev. Lett. 126, 241101 (2021).PRLTAO0031-900710.1103/PhysRevLett.126.241101], and it extends our investigation on the energy dependence of the spectral index of heavy elements. It reports the analysis of nickel data collected from November 2015 to May 2021 and a detailed assessment of the systematic uncertainties. In the region from 20 to 240  GeV/n our present data are compatible within the errors with a single power law with spectral index -2.51±0.07.

Time-resolved hadronic particle acceleration in the recurrent nova RS Ophiuchi

Recurrent novae are repeating thermonuclear explosions in the outer layers of white dwarfs, due to the accretion of fresh material from a binary companion. The shock generated when ejected material slams into the companion star's wind can accelerate particles. We report very-high-energy (VHE, [Formula: see text]) gamma rays from the recurrent nova RS Ophiuchi, up to a month after its 2021 outburst, observed using the High Energy Stereoscopic System. The VHE emission has a similar temporal profile to lower-energy GeV emission, indicating a common origin, with a two-day delay in peak flux. These observations constrain models of time-dependent particle energization, favoring a hadronic emission scenario over the leptonic alternative. Shocks in dense winds provide favorable environments for efficient acceleration of cosmic-rays to very high energies.

A Perspective on the Solar Modulation of Cosmic Anti-Matter

Global modulation studies with comprehensive numerical models contribute meaningfully to the refinement of very local interstellar spectra (VLISs) for cosmic rays. Modulation of positrons and anti-protons are investigated to establish how the ratio of their intensity, and with respect to electrons and protons, are changing with solar activity. This includes the polarity reversal of the solar magnetic field which creates a 22-year modulation cycle. Modeling illustrates how they are modulated over time and the particle drift they experience which is significant at lower kinetic energy. The VLIS for anti-protons has a peculiar spectral shape in contrast to protons so that the total modulation of anti-protons is awkwardly different to that for protons. We find that the proton-to-anti-proton ratio between 1–2 GeV may change by a factor of 1.5 over a solar cycle and that the intensity for anti-protons may decrease by a factor of ~2 at 100 MeV during this cycle. A composition is presented of VLIS for protons, deuteron, helium isotopes, electrons, and particularly for positrons and anti-protons. Gaining knowledge of their respective 11 and 22 year modulation is useful to interpret observations of low-energy anti-nuclei at the Earth as tests of dark matter annihilation.

Simulations of cosmic ray propagation

We review numerical methods for simulations of cosmic ray (CR) propagation on galactic and larger scales. We present the development of algorithms designed for phenomenological and self-consistent models of CR propagation in kinetic description based on numerical solutions of the Fokker–Planck equation. The phenomenological models assume a stationary structure of the galactic interstellar medium and incorporate diffusion of particles in physical and momentum space together with advection, spallation, production of secondaries and various radiation mechanisms. The self-consistent propagation models of CRs include the dynamical coupling of the CR population to the thermal plasma. The CR transport equation is discretized and solved numerically together with the set of MHD equations in various approaches treating the CR population as a separate relativistic fluid within the two-fluid approach or as a spectrally resolved population of particles evolving in physical and momentum space. The relevant processes incorporated in self-consistent models include advection, diffusion and streaming propagation as well as adiabatic compression and several radiative loss mechanisms. We discuss, applications of the numerical models for the interpretation of CR data collected by various instruments. We present example models of astrophysical processes influencing galactic evolution such as galactic winds, the amplification of large-scale magnetic fields and instabilities of the interstellar medium.

Erratum: Properties of a New Group of Cosmic Nuclei: Results from the Alpha Magnetic Spectrometer on Sodium, Aluminum, and Nitrogen [Phys. Rev. Lett. 127, 021101 (2021)]

This corrects the article DOI: 10.1103/PhysRevLett.127.021101.

Stochastic Fluctuations of Low-Energy Cosmic Rays and the Interpretation of Voyager Data

Data from the Voyager probes have provided us with the first measurement of cosmic ray intensities at MeV energies, an energy range that had previously not been explored. Simple extrapolations of models that fit data at GeV energies, e.g., from AMS-02, however, fail to reproduce the Voyager data in that the predicted intensities are too high. Oftentimes, this discrepancy is addressed by adding a break to the source spectrum or the diffusion coefficient in an ad hoc fashion, with a convincing physical explanation yet to be provided. Here, we argue that the discrete nature of cosmic ray sources, which is usually ignored, is instead a more likely explanation. We model the distribution of intensities expected from a statistical model of discrete sources and show that its expectation value is not representative but has a spectral shape different from that for a typical configuration of sources. The Voyager proton and electron data are however compatible with the median of the intensity distribution.

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.

Gamma Rays as Probes of Cosmic-Ray Propagation and Interactions in Galaxies

Continuum gamma-ray emission produced by interactions of cosmic rays with interstellar matter and radiation fields is a probe of non-thermal particle populations in galaxies. After decades of continuous improvements in experimental techniques and an ever-increasing sky and energy coverage, gamma-ray observations reveal in unprecedented detail the properties of galactic cosmic rays. A variety of scales and environments are now accessible to us, from the local interstellar medium near the Sun and the vicinity of cosmic-ray accelerators, out to the Milky Way at large and beyond, with a growing number of gamma-ray emitting star-forming galaxies. Gamma-ray observations have been pushing forward our understanding of the life cycle of cosmic rays in galaxies and, combined with advances in related domains, they have been challenging standard assumptions in the field and have spurred new developments in modelling approaches and data analysis methods. We provide a review of the status of the subject and discuss perspectives on future progress.

The TeV Cosmic-Ray Bump: A Message from the Epsilon Indi or Epsilon Eridani Star?

A recently observed bump in the cosmic-ray (CR) spectrum from 0.3 to 30 TV is likely caused by a stellar bow shock that reaccelerates preexisting CRs, which further propagate to the Sun along the magnetic field lines. Along their way, these particles generate an Iroshnikov-Kraichnan (I-K) turbulence that controls their propagation and sustains the bump. Ad hoc fitting of the bump shape requires six adjustable parameters. Our model requires none, merely depending on three physical unknowns that we constrain using the fit. These are the shock Mach number, M, its size, l ⊥, and the distance to it, ζ obs. Altogether, they define the bump rigidity R 0. With M ≈ 1.5-1.6 and R 0 ≈ 4.4 TV, the model fits the data with ≈0.08% accuracy. The fit critically requires the I-K spectrum predicted by the model and rules out the alternatives. These attributes of the fit make an accidental agreement highly unlikely. In turn, the R 0 and M derived from the fit impose the distance-size relation on the shock:ζ obs  ( pc ) ~ 10 2 l ⊥ ( pc ) . For sufficiently large bow shocks, l ⊥ = 10-3-10-2 pc, we find the distance of ζ obs = 3-10 pc. Three promising stars in this range are the Scholz's Star at 6.8 pc, Epsilon Indi at 3.6 pc, and Epsilon Eridani at 3.2 pc. Based on their current positions and velocities, we propose that Epsilon Indi and Epsilon Eridani can produce the observed spectral bump. Moreover, Epsilon Eridani's position is only ~6°.7 off of the magnetic field direction in the solar neighborhood, which also changes the CR arrival direction distribution. Given the proximity of these stars, the bump appearance may change in a relatively short time.

Diurnal Effect of Sub-GeV Dark Matter Boosted by Cosmic Rays

We point out a new type of diurnal effect for the cosmic ray boosted dark matter (DM). The DM-nucleon interactions not only allow the direct detection of DM with nuclear recoils but also allow cosmic rays to scatter with and boost the nonrelativistic DM to higher energies. If the DM-nuclei scattering cross sections are sufficiently large, the DM flux is attenuated as it propagates through the Earth, leading to a strong diurnal modulation. This diurnal modulation provides another prominent signature for the direct detection of boosted sub-GeV DM, in addition to signals with higher recoil energy.

A global environmental crisis 42,000 years ago

Geological archives record multiple reversals of Earth's magnetic poles, but the global impacts of these events, if any, remain unclear. Uncertain radiocarbon calibration has limited investigation of the potential effects of the last major magnetic inversion, known as the Laschamps Excursion [41 to 42 thousand years ago (ka)]. We use ancient New Zealand kauri trees (Agathis australis) to develop a detailed record of atmospheric radiocarbon levels across the Laschamps Excursion. We precisely characterize the geomagnetic reversal and perform global chemistry-climate modeling and detailed radiocarbon dating of paleoenvironmental records to investigate impacts. We find that geomagnetic field minima ~42 ka, in combination with Grand Solar Minima, caused substantial changes in atmospheric ozone concentration and circulation, driving synchronous global climate shifts that caused major environmental changes, extinction events, and transformations in the archaeological record.

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.

New Insights From Long-Term Measurements of Inner Belt Protons (10s of MeV) by SAMPEX, POES, Van Allen Probes, and Simulation Results

The Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX) mission provided long-term measurements of 10s of megaelectron volt (MeV) inner belt (L < 2) protons (1992-2009) as did the Polar-orbiting Operational Environmental Satellite-18 (POES-18, 2005 to present). These long-term measurements at low-Earth orbit (LEO) showed clear solar cycle variations which anticorrelate with sunspot number. However, the magnitude of the variation is much greater than the solar cycle variation of galactic cosmic rays (>GeV) that are regarded as a source of these trapped protons. Furthermore, the proton fluxes and their variations sensitively depend on the altitude above the South Atlantic Anomaly (SAA) region. With respect to protons (>36 MeV) mirroring near the magnetic equator, both POES measurements and simulations show no obvious solar cycle variations at L > 1.2. This is also confirmed by recent measurements from the Van Allen Probes (2012-2019), but there are clear solar cycle variations and a strong spatial gradient of the proton flux below L = 1.2. A direct comparison between measurements and simulations leads to the conclusion that energy loss of trapped protons due to collisions with free and bound electrons in the ionosphere and atmosphere is the dominant mechanism for the strong spatial gradient and solar cycle variation of the inner belt protons. This fact is also key of importance for spacecraft and instrument design and operation in near-Earth space.

Cosmic-Ray Database Update: Ultra-High Energy, Ultra-Heavy, and Antinuclei Cosmic-Ray Data (CRDB v4.0)

We present an update on CRDB, the cosmic-ray database for charged species. CRDB is based on MySQL, queried and sorted by jquery and table-sorter libraries, and displayed via PHP web pages through the AJAX protocol. We review the modifications made on the structure and outputs of the database since the first release (Maurin et al., 2014). For this update, the most important feature is the inclusion of ultra-heavy nuclei (Z>30), ultra-high energy nuclei (from 1015 to 1020 eV), and limits on antinuclei fluxes (Z≤−1 for A>1); more than 100 experiments, 350 publications, and 40,000 data points are now available in CRDB. We also revisited and simplified how users can retrieve data and submit new ones. For questions and requests, please contact crdb@lpsc.in2p3.fr.

Predictions of cognitive detriments from galactic cosmic ray exposures to astronauts on exploration missions

For the first-time we report on predictions on cognitive detriments from galactic cosmic ray (GCR) exposures on long-duration space missions outside the protection of the Earth's magnetosphere and solid body shielding. Estimates are based on a relative risk (RR) model of the fluence response for proton and heavy ion in rodent studies using the widely used novel object recognition (NOR) test, which estimates detriments in recognition or object memory. Our recent meta-analysis showed that linear and linear-quadratic dose response models were not accurate, while exponential increasing fluence response models based on particle track structure provided good descriptions of rodent data for doses up to 1 Gy. Using detailed models of the GCR environment and particle transport in shielding and tissue, we predict the excess relative risk (ERR) for NOR detriments for several long-term space mission scenarios. Predictions suggest ERR < 0.15 for most space mission scenarios with ERR<0.1 for 1-year lunar surface missions, and about ERR~0.1 for a 1000 day Mars mission for average solar cycle conditions. We discuss possible implications of these ERR levels of cognitive performance detriments relative to other neurological challenges such as rodent models of Alzheimer's disease (AD), Parkinson's disease (PD) and traumatic brain injury (TBI). Comparisons suggest a small but potentially clinically significant risk for possible space mission scenarios.

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.

Primordial Black Holes as a Dark Matter Candidate Are Severely Constrained by the Galactic Center 511 keV γ-Ray Line

We derive the strongest constraint on the fraction of dark matter that can be composed of low mass primordial black holes by using the observation of the Galactic Center 511 keV γ-ray line. Primordial black holes of masses ≲10^{15}  kg will evaporate to produce e^{±} pairs. The positrons will lose energy in the Galactic Center, become nonrelativistic, then annihilate with the ambient electrons. We derive robust and conservative bounds by assuming that the rate of positron injection via primordial black hole evaporation is less than what is required to explain the SPI/INTEGRAL observation of the Galactic Center 511 keV γ-ray line. Depending on the primordial black hole mass function and other astrophysical uncertainties, these constraints are the most stringent in the literature and show that primordial black holes contribute to less than 1% of the dark matter density. Our technique also probes part of the mass range which was completely unconstrained by previous studies.

H3+, the ideal probe for in situ measurement of the Galactic cosmic rays

The reason why [Formula: see text] is an ideal probe for in situ measurement of cosmic ray energy density in the Galaxy is discussed. The variation of measured values of cosmic ray ionization rate ζ of H₂ since the theoretical prediction by Spitzer & Tomasko (Spitzer & Tomasko 1968 Astrophys. J. 152, 971-986) is reviewed. How the measurements by [Formula: see text] established ζ as of the order of 10^-17 s^-1 for dense clouds, 10^-16 s^-1 for diffuse clouds in the Galactic disc, and 10^-14 s^-1 for warm and diffuse gas in the Galactic centre is discussed. This article is part of a discussion meeting issue 'Advances in hydrogen molecular ions: H₃⁺, H₅⁺ and beyond'.

Fermi-LAT Observations of γ-Ray Emission toward the Outer Halo of M31

The Andromeda galaxy is the closest spiral galaxy to us and has been the subject of numerous studies. It harbors a massive dark matter halo, which may span up to ~600 kpc across and comprises ~90% of the galaxy's total mass. This halo size translates into a large diameter of 42° on the sky, for an M31-Milky Way (MW) distance of 785 kpc, but its presumably low surface brightness makes it challenging to detect with γ-ray telescopes. Using 7.6 yr of Fermi Large Area Telescope (Fermi-LAT) observations, we make a detailed study of the γ-ray emission between 1-100 GeV toward M31's outer halo, with a total field radius of 60° centered at M31, and perform an in-depth analysis of the systematic uncertainties related to the observations. We use the cosmic-ray propagation code GALPROP to construct specialized interstellar emission models to characterize the foreground γ-ray emission from the MW, including a self-consistent determination of the isotropic component. We find evidence for an extended excess that appears to be distinct from the conventional MW foreground, having a total radial extension upward of ~120-200 kpc from the center of M31. We discuss plausible interpretations of the excess emission, but emphasize that uncertainties in the MW foreground-and in particular, modeling of the H i-related components-have not been fully explored and may impact the results.

Cosmic-Ray Propagation in Light of the Recent Observation of Geminga

The High Altitude Water Cherenkov (HAWC) telescope recently observed extended emission around the Geminga and PSR B0656+14 pulsar wind nebulae (PWNe). These observations have been used to estimate cosmic-ray (CR) diffusion coefficients near the PWNe that appear to be more than two orders of magnitude smaller than the diffusion coefficients typically derived for the interstellar medium from measured abundances of secondary species in CRs. Two-zone diffusion models have been proposed as a solution to this discrepancy, where the slower diffusion zone (SDZ) is confined to a small region around the PWN. Such models are shown to successfully reproduce the HAWC observations of the Geminga PWN while retaining consistency with other CR data. It is found that the size of the SDZ influences the predicted positron flux and the spectral shape of the extended γ-ray emission at lower energies that can be observed with the Fermi Large Area Telescope. If the two observed PWNe are not unique, then it is likely that there are similar pockets of slow diffusion around many CR sources elsewhere in the Milky Way. The consequences of such a picture for Galactic CR propagation is explored.

Risks of cognitive detriments after low dose heavy ion and proton exposures

Purpose: Heavy ion and proton brain irradiations occur during space travel and in Hadron therapy for cancer. Heavy ions produce distinct patterns of energy deposition in neuron cells and brain tissues compared to X-rays leading to large uncertainties in risk estimates. We make a critical review of findings from research studies over the last 25 years for understanding risks at low dose. Conclusions: A large number of mouse and rat cognitive testing measures have been reported for a variety of particle species and energies for acute doses. However, tissue reactions occur above dose thresholds and very few studies were performed at the heavy ion doses to be encountered on space missions (<0.04 Gy/y) or considered dose-rate effects, such that threshold doses are not known in rodent models. Investigations of possible mechanisms for cognitive changes have been limited by experimental design with largely group specific and not subject specific findings reported. Persistent oxidative stress and activated microglia cells are common mechanisms studied, while impairment of neurogenesis, detriments in neuron morphology, and changes to gene and protein expression were each found to be important in specific studies. Future research should focus on estimating threshold doses carried out with experimental designs aimed at understating causative mechanisms, which will be essential for extrapolating rodent findings to humans and chronic radiation scenarios, while establishing if mitigation are needed.

Progress in the Global Modeling of the Galactic Magnetic Field

We discuss the global modeling of the properties of the Galactic Magnetic Field (GMF). Several improvements and variations of the model of the GMF from Jansson & Farrar (2012) (JF12) are investigated in an analysis constrained by all-sky rotation measures of extragalactic sources and polarized and unpolarized synchrotron emission data from WMAP and Planck. We present the impact of the investigated model variations on the propagation of ultrahigh-energy cosmic rays in the Galaxy

Practical Modeling of Large-Scale Galactic Magnetic Fields: Status and Prospects

This is a review of the status of efforts to model the large-scale Galactic magnetic field (GMF). Though important for a variety of astrophysical processes, the GMF remains poorly understood despite some interesting new tracers being used in the field. Though we still have too many models that might fit the data, this is not to say that the field has not developed in the last few years. In particular, surveys of polarized dust have given us a new observable that is complementary to the more traditional radio tracers, and a variety of other new tracers and related measurements are becoming available to improve current modeling. This paper reviews: the tracers available; the models that have been studied; what has been learned so far; what the caveats and outstanding issues are; and one opinion of where the most promising future avenues of exploration lie.

Voyager 1 e^{±} Further Constrain Primordial Black Holes as Dark Matter

Primordial black holes (PBHs) with a mass M≲10^{17}  g are expected to inject sub-GeV electrons and positrons in the Galaxy via Hawking radiation. These cosmic rays are shielded by the solar magnetic field for Earth-bound detectors, but not for Voyager 1, which is now beyond the heliopause. We use its data to constrain the fraction of PBHs to the dark matter in the Galaxy, finding that PBHs with M<10^{16}  g cannot contribute more than 0.1% (or less for a log-normal mass distribution). Our limits are based on local Galactic measurements and are thus complementary to those derived from cosmological observations.

Testing Diffusion of Cosmic Rays in the Heliosphere with Proton and Helium Data from AMS

After six years of continuous observations in space, the Alpha Magnetic Spectrometer experiment has released new data on the temporal evolution of the proton and helium fluxes in cosmic rays. These data revealed that the ratio between proton and helium fluxes at the same value of rigidity R=p/Z (momentum/charge ratio) is not constant at R≲3  GV. In particular, the ratio is found to decrease steadily during the descending phase of Solar Cycle 24 toward the next minimum. We show that such a behavior is a remarkable signature of the β×λ(R) dependence in the diffusion of cosmic rays in heliosphere, where β is their adimensional speed and λ(R) is their mean free path, a universal function of rigidity for all nuclei. This dependence is responsible for distinctive charge or mass dependent effects in the time-dependent modulation of low-rigidity particles.

Molecular Ionization Rates and Ultracompact Dark Matter Minihalos

Molecular ionization in the central molecular zone of our Galaxy is enhanced over the typical galactic value by an order of magnitude or more. This cannot be easily explained for dense Galactic center molecular complexes in the absence of embedded sources of low energy cosmic rays. We provide such a source in the form of ultracompact minihalos of self-annihilating dark matter for a variety of annihilation channels that depend on the particle mass and model. Such sources are motivated for plausible inflationary power spectrum parameters, and while possibly subdominant in terms of the total dark matter mass within the Galactic bulge, might also account for, or at least not be in tension with, the Fermi Galactic center γ-ray excess.

Observation of Fine Time Structures in the Cosmic Proton and Helium Fluxes with the Alpha Magnetic Spectrometer on the International Space Station

We present the precision measurement from May 2011 to May 2017 (79 Bartels rotations) of the proton fluxes at rigidities from 1 to 60 GV and the helium fluxes from 1.9 to 60 GV based on a total of 1×10^{9} events collected with the Alpha Magnetic Spectrometer aboard the International Space Station. This measurement is in solar cycle 24, which has the solar maximum in April 2014. We observed that, below 40 GV, the proton flux and the helium flux show nearly identical fine structures in both time and relative amplitude. The amplitudes of the flux structures decrease with increasing rigidity and vanish above 40 GV. The amplitudes of the structures are reduced during the time period, which started one year after solar maximum, when the proton and helium fluxes steadily increase. Above ∼3  GV the p/He flux ratio is time independent. We observed that below ∼3  GV the ratio has a long-term decrease coinciding with the period during which the fluxes start to rise.