Cosmic Ray Species Research Articles

Non-parametric determination of H and He interstellar fluxes from cosmic-ray data

Top-of-atmosphere (TOA) cosmic-ray (CR) fluxes from satellites and balloon-borne experiments are snapshots of the solar activity imprinted on the interstellar (IS) fluxes. Given a series of snapshots, the unknown IS flux shape and the level of modulation (for each snapshot) can be recovered. We wish (i) to provide the most accurate determination of the IS H and He fluxes from TOA data only, (ii) to obtain the associated modulation levels (and uncertainties) fully accounting for the correlations with the IS flux uncertainties, and (iii) to inspect whether the minimal Force-Field approximation is sufficient to explain all the data at hand. Using H and He TOA measurements, including the recent high precision AMS, BESS-Polar and PAMELA data, we perform a non-parametric fit of the IS fluxes $J^{\rm IS}_{\rm H,~He}$ and modulation level $\phi_i$ for each data taking period. We rely on a Markov Chain Monte Carlo (MCMC) engine to extract the PDF and correlations (hence the credible intervals) of the sought parameters. Despite H and He being the most abundant and best measured CR species, several datasets must be excluded from the analysis due to inconsistencies with other measurements. From the subset of data passing our consistency cut, we provide ready-to-use best-fit and credible intervals for the H and He IS fluxes from MeV/n to PeV/n energy (with a relative precision in the range $[2-10\%]$ at 1$\sigma$). Given the strong correlation between $J^{\rm IS}$ and $\phi_i$ parameters, the uncertainties on $J^{\rm IS}$ translate into $\Delta\phi\approx \pm 30$ MV (at 1$\sigma$) for all experiments. We also find that the presence of $^3$He in He data biases $\phi$ towards higher $\phi$ values by $\sim 30$ MV. The force-field approximation, despite its limitation, gives an excellent ($\chi^2/$dof$=1.02$) description of the recent high-precision TOA H and He fluxes.

Antideuterons from decaying gravitino dark matter

From the point of view of Dark Matter searches, antideuteron are an interesting species of cosmic rays. In the past eighteen years much progress has been done in the understanding of these cosmic rays but this did not translate into experimental efforts to actually detect these particles. In this work, starting from a specific case of dark matter candidate with relatively generic features, we study some sources of theoretical uncertainties and show that the margin of improvement beyond the limits already set by cosmic-ray antiproton are quite narrow and anyway unachievable for the next generation of experiments.

Cosmic-ray slowing down in molecular clouds: Effects of heavy nuclei

Context. A cosmic ray (CR) spectrum propagated through ISM contains very few low-energy (<100 MeV) particles. Recently, a local CR spectrum, with strong low energy components, has been proposed to be responsible for the over production of H + molecule in some molecular clouds. Aims. We aim to explore the effects of the chemical composition of low-energy cosmic rays (CRs) when they slow down in dense molecular clouds without magnetic fields. We considered both ionization and solid material processing rates. Methods. We used galatic CR chemical composition from proton to iron. We propagated two types of CR spectra through a cloud made of H2: those CR spectra with different contents of low energy CRs and those assumed to be initially identical for all CR species. The stopping and range of ions in matter (SRIM) package provided the necessary stopping powers. The ionization rates were computed with cross sections from recent semi-empirical laws, while effective cross sections were parametrized for solid processing rates using a power law of the stopping power (power 1 to 2). Results. The relative contribution to the cloud ionization of proton and heavy CRs was found identical everywhere in the irradiated cloud, no matter which CR spectrum we used. As compared to classical calculations, using protons and high-energy behaviour of ionization processes (Z 2 scaling), we reduced absolute values of ionization rates by few a tens of percents but only in the case of spectrum with a high content of low-energy CRs. We found, using the same CR spectrum, the solid material processing rates to be reduced between the outer and inner part of thick cloud by a factor 10 (as in case of the ionization rates) or by a factor 100, depending on the type of process.

Cosmic ray transport in astrophysical plasmas

Since the development of satellite space technology about 50 years ago the solar heliosphere is explored almost routinely by several spacecrafts carrying detectors for measuring the properties of the interplanetary medium including energetic charged particles (cosmic rays), solar wind particle densities, and electromagnetic fields. In 2012, the Voyager 1 spacecraft has even left what could be described as the heliospheric modulation region, as indicated by the sudden disappearance of low energy heliospheric cosmic ray particles. With the available in-situ measurements of interplanetary turbulent electromagnetic fields and of the momentum spectra of different cosmic ray species in different interplanetary environments, the heliosphere is the best cosmic laboratory to test our understanding of the transport and acceleration of cosmic rays in space plasmas. I review both the historical development and the current state of various cosmic ray transport equations. Similarities and differences to transport theories for terrestrial fusion plasmas are highlighted. Any progress in cosmic ray transport requires a detailed understanding of the electromagnetic turbulence that is responsible for the scattering and acceleration of these particles.

IMPULSIVE SPOT HEATING AND THERMAL EXPLOSION OF INTERSTELLAR GRAINS REVISITED

The problem of impulsive heating of dust grains in cold, dense interstellar clouds is revisited theoretically, with the aim to better understand leading mechanisms of the explosive desorption of icy mantles. It is rigorously shown that if the heating of a reactive medium occurs within a sufficiently localized spot (e.g., heating of mantles by cosmic rays), then the subsequent thermal evolution is characterized by a single dimensionless number $\lambda$. This number identifies a bifurcation between two distinct regimes: When $\lambda$ exceeds a critical value (threshold), the heat equation exhibits the explosive solution, i.e., the thermal (chemical) explosion is triggered. Otherwise, thermal diffusion causes the deposited heat to spread over the entire grain -- this regime is commonly known as the whole-grain heating. The theory allows us to find a critical combination of the physical parameters that govern the explosion of icy mantles due to impulsive spot heating. In particular, the calculations suggest that heavy cosmic ray species (e.g., iron ions) colliding with dust are able to trigger the explosion. Based on the recently calculated local cosmic-ray spectra, the expected rate of the explosive desorption is estimated. The efficiency of the desorption, which affects all solid species independent of their binding energy, is shown to be comparable with other cosmic-ray desorption mechanisms typically considered in the literature. Also, the theory allows us to estimate maximum abundances of reactive species that may be stored in the mantles, which provides important constraints on available astrochemical models.

Design and performances of prototype muon detectors of LHAASO-KM2A

In the Large High Altitude Air Shower Observatory project, 1221 muon detectors will be built to discriminate cosmic gamma rays from nuclei and between cosmic ray species with energies in 1013–1018eV. The muon detector has been designed to meet the physical requirements, according to which two muon detector prototypes have been built at the Yangbajing Cosmic Ray Observatory at the end of 2012 and in October 2014, respectively. The first one has been stably operated for more than 2years. The performances of the two detectors have been evaluated, especially for position dependence and punch-through effect. In this paper, the design and the performances of the prototype muon detectors are presented.

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 15GeV-1TeV). 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.

Constraining the origin of the rising cosmic ray positron fraction with the boron-to-carbon ratio

The rapid rise in the cosmic ray positron fraction above 10 GeV, as measured by PAMELA and AMS, suggests the existence of nearby primary sources of high energy positrons, such as pulsars or annihilating/decaying dark matter. In contrast, the spectrum of secondary positrons produced through the collisions of cosmic rays in the interstellar medium is predicted to fall rapidly with energy, and thus is unable to account for the observed rise. It has been proposed, however, that secondary positrons could be produced and then accelerated in nearby supernova remnants, potentially explaining the observed rise, without the need of primary positron sources. Yet, if secondary positrons are accelerated in such shocks, other secondary cosmic ray species (such as boron nuclei, and antiprotons) will also be accelerated, leading to rises in the boron-to-carbon and antiproton-to-proton ratios. The measurements of the boron-to-carbon ratio by the PAMELA and AMS collaborations, however, show no sign of such a rise. With this new data in hand, we revisit the secondary acceleration scenario for the rising positron fraction. Assuming that the same supernova remnants accelerate both light nuclei (protons, helium) and heavier cosmic ray species, we find that no more than ~25% of the observed rise in the positron fraction can result from this mechanism (at the 95% confidence level)

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.

Revisiting the hardening of the cosmic ray energy spectrum at TeV energies

Measurements of cosmic rays by experiments such as ATIC, CREAM, and PAMELA indicate a hardening of the cosmic-ray energy spectrum at TeV energies. In our recent work (Thoudam & H\"orandel 2012a), we showed that the hardening can be due to the effect of nearby supernova remnants. We showed it for the case of proton and helium. In this paper, we present an improved and detailed version of our previous work, and extend our study to heavier cosmic-ray species such as boron, carbon, oxygen, and iron nuclei. Unlike our previous study, the present work involves a detailed calculation of the background cosmic rays and follows a consistent treatment of cosmic-ray source parameters between the background and the nearby components. Moreover, we also present a detailed comparison of our results on the secondary-to-primary ratios, secondary spectra, and the diffuse gamma-ray spectrum with the results expected from other existing models, which can be checked by future measurements at high energies.

NaI Detector Network at Aragats

The Aragats Space Environmental Center (ASEC) [1] provides monitoring of different species of secondary cosmic rays and consists of two high altitude research stations on Mt. Aragats in Armenia. Along with solar modulation effects, ASEC detectors register several coherent enhancements associated with thunderstorm activity. The experimental techniques used allowed for the first time to simultaneously measure fluxes of the electrons, muons, gamma rays, and neutrons correlated with thunderstorm activity [2,3]. Ground-based observations by a complex of surface particle detectors, measuring in systematically and repeatable fashion, gamma quanta, electrons, muons and neutrons from atmospheric sources are necessary for proving the theory of particle acceleration and multiplication during thunderstorms. Energy spectra and correlations between fluxes of different particles, measured on Earth's surface address the important issues of research of the solar modulation effects and the atmospheric high-energy phenomena. In May 26 2011, launched 5 NaI(Tl) (thallium-doped sodium iodide) scintillation detectors and 1 plastic one in the new ASEC laboratory on Aragats to detect low energy gamma rays from the thunderclouds and short particle bursts. Including NaI(Tl) detectors in ASEC detectors system is of great importance for investigation thunderstorm phenomena because NaI(Tl) detectors have high efficiency of gamma ray detecting in comparison with plastic ones.

Space Environmental Viewing and Analysis Network (SEVAN) – characteristics and first operation results

Space Environmental Viewing and Analysis Network is a worldwide network of identical particle detectors located at middle and low latitudes aimed to improve fundamental research of space weather conditions and to provide short- and long-term forecasts of the dangerous consequences of space storms. SEVAN detected changing fluxes of different species of secondary cosmic rays at different altitudes and latitudes, thus turning SEVAN into a powerful integrated device used to explore solar modulation effects. Till to now the SEVAN modules are installed at Aragats Space Environmental Centre in Armenia (3 units at altitudes 800, 2000 and 3200 m a.s.l.), Bulgaria (Moussala), Croatia and India (New-Delhi JNU.) and now under installation in Slovakia, LomnitskySchtit). Recently SEVAN detectors were used for research of new high-energy phenomena originated in terrestrial atmosphere – Thunderstorm Ground Enhancements (TGEs). In 2011 first joint measurements of solar modulation effects were detected by SEVAN network, now under analysis.

A study of the performance parameters of the High Altitude Gamma Ray (HAGAR) telescope system at Ladakh in India

We present results of Monte Carlo simulations for the High Altitude Gamma Ray (HAGAR) telescope array which detects very high energy gamma rays from astronomical sources. This telescope array, located at Hanle at an altitude of 4270m in the Ladakh region of the Himalayas in India, is the highest altitude atmospheric Cherenkov experiment in the world. Taking advantage of the high altitude, this experiment could achieve relatively low energy threshold with a modest mirror area coverage. To understand the performance parameters of this telescope system, we have simulated large samples of extensive air showers initiated by gamma rays and various species of cosmic rays, using the CORSIKA package. Cherenkov photons produced in the atmosphere are sampled at ground level. These photons are then passed through the detector simulation program, which takes into account various design details and the data acquisition system of HAGAR. Night sky photons are also considered in the detector simulation program as performance of the telescope depends strongly on the level of night sky background (NSB) at the observation site. We have estimated various performance parameters like energy threshold and effective area for vertically incident showers as well as inclined showers. Details of these parameters, results obtained from simulations and comparison with the observed data are presented. It is shown that the energy threshold of the HAGAR telescope system is about 208GeV, a factor of ∼4 less than for a similar set up at about 1000m altitude, and it is able to detect Crab like sources at 5σ significance in 17h of observation without imposing additional criteria, like gamma-hadron separation, for further rejection of cosmic rays.

Indirect searches for gravitino dark matter

The gravitino in models with a small violation of R-parity is a well-motivated decaying dark matter candidate that leads to a cosmological scenario that is consistent with big bang nucleosynthesis and thermal leptogenesis. The gravitino lifetime is cosmologically long-lived since its decays are suppressed by the Planck-scale as well as the small R-parity violating parameter. We discuss the signals in different cosmic-ray species coming from the decay of gravitino dark matter, namely gamma rays, positrons, antiprotons, antideuterons and neutrinos. Comparison to cosmic-ray data can be used to constrain the parameters of the model.

TESTING THE ORIGIN OF HIGH-ENERGY COSMIC RAYS

Recent accurate measurements of cosmic-ray (CR) protons and nuclei by ATIC-2, CREAM, and PAMELA reveal: a) unexpected spectral hardening in the spectra of CR species above a few hundred GeV per nucleon, b) a harder spectrum of He compared to protons, and c) softening of the CR spectra just below the break energy. These newly-discovered features may offer a clue to the origin of the observed high-energy Galactic CRs. We discuss possible interpretations of these spectral features and make predictions for the secondary CR fluxes and secondary to primary ratios, anisotropy of CRs, and diffuse Galactic {\gamma}-ray emission in different phenomenological scenarios. Our predictions can be tested by currently running or near-future high-energy astrophysics experiments.

Status of the large high altitude air shower observatory project

The Large High Altitude Air Shower Observatory (LHAASO) project is a multipurpose project. The main scientific tasks can be summarized as follows: (1) searching for galactic cosmic ray origins through gamma ray source detection above 30TeV; (2) wide field of view survey for gamma ray sources at energies higher than 100GeV; (3) energy spectrum measurements for individual cosmic ray species from 30TeV to 10PeV. To target above tasks, a complex detector array is designed. This paper describes the progress on the research and development of all kind of detectors. Construction and operation of a prototype detector array at Tibet site with 4300m a.s.l. are also presented.

Cosmic Ray Signatures of Decaying Dark Matter

Astrophysical and cosmological observations do not require the dark matter particles to be absolutely stable. If they are indeed unstable, their decay into Standard Model particles might occur at a sufficiently large rate to allow the indirect detection of dark matter through an anomalous contribution to the high energy cosmic ray fluxes. We analyze the implications of the excess in the total electron plus positron flux and the positron fraction reported by the Fermi and PAMELA collaborations, respectively, for the scenario of decaying dark matter. We also discuss the constraints on this scenario from measurements of other cosmic ray species and the predictions for the diffuse gamma ray flux and the neutrino flux. In particular, we expect a sizable dipole-like anisotropy which may be observed in the near future by the Fermi-LAT.

The GeV-TeV Galactic gamma-ray diffuse emission

The Galactic gamma-ray diffuse emission is currently observed in the GeV-TeV energy range with unprecedented accuracy by the Fermi satellite. Understanding this component is crucial as it provides a background to many different signals such as extragalactic sources or annihilating dark matter. It is timely to reinvestigate how it is calculated and to assess the various uncertainties which are likely to affect the accuracy of the predictions. The Galactic gamma-ray diffuse emission is mostly produced above a few GeV by the interactions of cosmic ray primaries impinging on the interstellar material. The theoretical error on that component is derived by exploring various potential sources of uncertainty. Particular attention is paid to cosmic ray propagation. Nuclear cross sections, the proton and helium fluxes at the Earth, the Galactic radial profile of supernova remnants and the hydrogen distribution can also severely affect the signal. The propagation of cosmic ray species throughout the Galaxy is described in the framework of a semi-analytic two-zone diffusion/convection model. This allows to convert the constraints set by the boron-to-carbon data into a theoretical uncertainty on the diffuse emission. New deconvolutions of the HI and CO sky maps are also used to get the hydrogen distribution within the Galaxy. The thickness of the cosmic ray diffusive halo is found to have a significant effect on the Galactic gamma-ray diffuse emission while the interplay between diffusion and convection has little influence on the signal. The uncertainties related to nuclear cross sections and to the primary cosmic ray fluxes at the Earth are significant. The radial distribution of supernova remnants along the Galactic plane turns out to be a key ingredient. As expected, the predictions are extremely sensitive to the spatial distribution of hydrogen within the Milky Way.

High energy cosmic rays from decaying dark matter

AbstractWe discuss the implications of the excess in the total electron plus positron flux and the positron fraction reported by the Fermi and PAMELA collaborations, respectively, for the scenario of decaying dark matter. We also discuss the constraints on this scenario from measurements of other cosmic ray species and the predictions for the diffuse gamma ray flux and the neutrino flux. In particular, we expect a sizable dipole‐like anisotropy in the diffuse gamma‐ray flux which may be observed by the Fermi‐LAT in the near future.