Cosmic Rays Research Articles

Evolution of the Secondary Electron Spectrum during Cosmic-Ray Discharge in the Universe

We recently found that streaming cosmic rays (CRs) induce a resistive electric field that can accelerate secondary electrons produced by CR ionization. In this work, we study the evolution of the energy spectrum of secondary electrons by numerically solving the one-dimensional Boltzmann equation and Ohm’s law. We show that the accelerated secondary electrons further ionize a gas, that is, an electron avalanche occurs, resulting in increased ionization and excitation of the gas. Although the resistive electric field becomes weaker than the one before the CR discharge, the weak resistive electric field weakly accelerates the secondary electrons. The quasi-steady state is almost independent of the initial resistive electric field but depends on the electron fraction in the gas. The resistive electric field in the quasi-steady state is larger for the higher electron fraction, which makes the number of secondary electrons that can ionize the gas larger, resulting in a higher ionization rate. The CR discharge could explain the high ionization rate that is observed in some molecular clouds.

The Algorithm of the Two Neutron Monitors for the Analysis of the Rigidity Spectrum Variations of Galactic Cosmic Ray Intensity Flux in Solar Cycle 24

The method of the two neutron monitors was used to analyze the parameters of the rigidity spectrum variations (RSV) of galactic cosmic ray intensity (GCR) flux in solar cycle 24 based on the data from the global network of neutron monitors. This method is an alternative to the least squares method when there are few monitors working stably in a given period, and their use in the least squares method is impossible. Analyses of the changes in exponent γ in the RSV of GCR flux from 2009 to 2019 were studied. The soft RSV (γ = 1.2–1.3) of the GCR flux around the maximum epoch and the hard RSV (γ = 0.6–0.9) around the minimum epoch of solar activity (SA) is the general feature of GCR modulation in the GeV energy scale (5, 50), to which neutron monitors were found to correspond. Therefore, various values of the RSV γ in the considered period show that during the decrease and increase period of SA, the essential changes in the large-scale structure of the heliospheric magnetic field (HMF) fluctuations/turbulence take place. The exponent γ of the RSV of the GCR flux can be considered a significant parameter to investigate the long-period changes in the GCR flux.

X-ray polarization: A view deep inside cosmic ray driven turbulence and particle acceleration in supernova remnants

X-ray polarization: A view deep inside cosmic ray driven turbulence and particle acceleration in supernova remnants

Prospects for joint reconstruction of imaging air Cherenkov telescope array and extensive air shower array

In this paper, we proposed a joint reconstruction of γ-ray events using both extensive air array (EAS) and Imaging air Cherenkov Telescope array (IACT). We considered eight Cherenkov telescopes to be built on the LHAASO (Large High Altitude Air Shower Observatory) site and investigate the improvement in differential sensitivity when combining the information from both IACT and Moun detectors of LHAASO-KM2A. We found that due to the higher cosmic ray background rejection power and higher gamma ray retention ratio provided by muon detectors of LHAASO, such a joint reconstruction can significantly improve the sensitivity of IACTs, especially for extended sources and long exposure time. The sensitivity of the eight-telescopes array can be improved by 25%−60% in different energy ranges using joint reconstructions with muon detectors.

Blind and robust estimation of adaptive optics point spread function and diffuse halo with sharp-edged objects

Context. Initially designed to detect and characterise exoplanets, extreme adaptive optics (AO) systems open a new window onto the Solar System by resolving its small bodies. Nonetheless, their study remains limited by the accuracy of the knowledge of the AO-corrected point spread function (AO-PSF) that degrades their image and produces a bright halo, potentially hiding faint moons in their close vicinity. Aims. To overcome the random nature of AO-PSFs, I aim to develop a method that blindly recovers the PSF and its faint structured extensions directly into the data of interest, without any prior on the instrument or the object’s shape. The objectives are both to deconvolve the object and to properly estimate and remove its surrounding halo to highlight potential faint companions. Methods. My method first estimated the PSF core via a parametric model fit, under the assumption of a sharp-edged flat object. Then, the resolved object and the PSF extensions were alternatively deconvolved with a robust method, insensitive to model outliers, such as cosmic rays or unresolved moons. Finally, the complex halo produced by the AO system was modelled and removed from the data. Results. The method is validated on realistic simulations with an on-sky AO-PSF from the SPHERE/ZIMPOL instrument. On real data, the proposed blind deconvolution algorithm strongly improves the image sharpness and retrieves details on the surface of asteroids. In addition, their moons are visible in all tested epochs despite important variability in turbulence conditions. Conclusions. My method shows the feasibility of retrieving the complex features of AO-PSFs directly from the data of interest. It paves the way towards more precise studies of asteroid surfaces and the discovery and characterisation of Solar System moons in archival data or with future instruments on extremely large telescopes with ever more complex AO-PSFs.

Increasing the period of active use of on-board electronic equipment of spacecraft

For electronic equipment of space systems, and primarily memory devices, the task of protection from the effects of ionizing cosmic radiation and other external factors that distort stored and processed information is relevant. This paper proposes a holographic coding method that allows you to restore information in the event of a large number of errors. The method is based on recording into memory, instead of the original digital hologram data, a virtual digital object corresponding to a data block. The divisibility property of a hologram is used, which makes it possible to reconstruct a recorded data block from its fragment. The achieved level of noise immunity is determined by the size of the hologram. For an 8-bit data block, recording a 256-bit hologram provides information recovery if 75 % of the recorded hologram is lost. The developed decoder corrects a package of dependent (grouping) errors that distort all bits of the hologram. The number of random independent errors that the decoder corrects can be up to 40 % of the recorded information. The information storage system, resistant to ionizing radiation, is a memory array of increased capacity, taking into account the selected redundancy factor, and a memory controller that performs holographic encoding when recording information and decoding with automatic error correction when reading information. The operating algorithm of the controller itself can be implemented in the form of a programmable logic integrated circuit, or stored in a read-only memory device that is not affected by ionizing radiation.

Cosmic-Ray Transport in the Presence of a Fisk-type Heliospheric Magnetic Field: Investigating the Influence of Drift

Drifts due to the gradients in, and curvatures of, the heliospheric magnetic field (HMF) play a significant role in the transport of galactic cosmic rays (GCRs) in the heliosphere. Although this has been well studied for the Parker HMF, the influence of Fisk-type fields, with their unique geometry, has hitherto received less attention. Here, drift velocity profiles computed for a Schwadron–Parker hybrid field are compared with those for a purely Parkerian field. Furthermore, the influence of this field on GCR modulation, as computed with a 3D, ab initio, modulation model, is investigated. Globally, the differences between the computed intensities are small. Nevertheless, local measures, such as azimuthal variations in the GCR intensities, show a significant influence of the Schwadron–Parker hybrid field on GCR transport.

Upper limit on the coronal cosmic ray energy budget in Seyfert galaxies

Abstract The IceCube Collaboration has reported possible detections of high-energy neutrinos from nearby Seyfert galaxies. While central hot coronae are proposed as the primary neutrino production site, the exact coronal cosmic ray energy budget has been loosely constrained. In this study, we propose a new stringent upper bound on the coronal cosmic ray energy budget of Seyfert galaxies, considering both accretion dynamics and observed properties of radio-quiet Seyfert galaxies. Notably, even under the calorimetric condition where cosmic rays lose all their energy, our limit indicates that the coronal neutrino flux of NGC 1068 is about an order of magnitude fainter than the observed levels. This discrepancy suggests the need for further theoretical and observational investigations on the IceCube signals from Seyfert galaxies.

Comparative analysis of leukemia and risk estimation in working age population between provinces of Ecuador.

Leukemia is associated with exposure to radiation, benzene derivatives, and pesticides. Previous research has documented an increase in work-related leukemia in the Latin American Andean region. To date, there are only few studies in Ecuador on the impact of oil exploitation on adjacent indigenous communities. Our study aims to show the impact of leukemia on the working-age population. For the calculation of morbidity and mortality rates, we used hospital discharge and death records from the National Institute of Statistics of Ecuador. These data were collected and adjusted to the corresponding province's population for further analysis. Large differences were observed between provinces in adjusted rates of leukemia mortality and morbidity in the working-age population. The variations in altitude among different areas in Ecuador give the provinces a distinct geographic identity. Likewise, the provinces with the highest morbidity and mortality rankings, such as Azuay, Loja, Imbabura, and Tungurahua, have an average altitude above 2000 meters. As a result, there are variations in the average temperature, exposure to solar and cosmic radiation, and mining and farming methods. The observed differences warrant the future collection of geolocation data for affected individuals. This could help to better understand how leukemia cases have demogrpahic hotspots in the country, identify possible risk factors associated with the disease in each region, and design more effective prevention and control strategies.

Can the Symmetric Fermi and eROSITA Bubbles Be Produced by Tilted Jets?

The Fermi Gamma-Ray Space Telescope reveals two large bubbles in the Galaxy, extending nearly symmetrically ∼50° above and below the Galactic center (GC). Previous simulations of bubble formation invoking active galactic nucleus (AGN) jets have assumed that the jets are vertical to the Galactic disk; however, in general, the jet orientation does not necessarily correlate with the rotational axis of the Galactic disk. Using three-dimensional special relativistic hydrodynamic simulations including cosmic rays (CRs) and thermal gas, we show that the dense clumpy gas within the Galactic disk disrupts jet collimation (“failed jets” hereafter), which causes the failed jets to form hot bubbles. Subsequent buoyancy in the stratified atmosphere renders them vertical to form the symmetric Fermi and eROSITA bubbles (collectively, Galactic bubbles). We find that (1) despite the relativistic jets emanating from the GC at various angles ≤45° with respect to the rotational axis of the Galaxy, the Galactic bubbles nonetheless appear aligned with the axis; (2) the edge of the eROSITA bubbles corresponds to a forward shock driven by the hot bubbles; (3) followed by the forward shock is a tangling contact discontinuity corresponding to the edge of the Fermi bubbles; (4) assuming a leptonic model we find that the observed gamma-ray bubbles and microwave haze can be reproduced with a best-fit CR power-law spectral index of 2.4; The agreements between the simulated and the observed multiwavelength features suggest that forming the Galactic bubbles by oblique AGN failed jets is a plausible scenario.

Extracting the True Cosmic-ray Relative Intensity Sky Map from Normalized Relative Intensity Data Measured by Air Shower Experiments

We use data from the Tibet AS γ experiment for 4 teraelectronvolt (TeV) cosmic rays as an example to perform a nonlinear interstellar distribution model regression according to the way the observed anisotropy is typically presented, from which we extract normalization factors that allow us to obtain a true relative intensity sky map from the measurements. By using various test statistics, we show that the nonlinear fit significantly outperforms the direct linear fit in its ability to model cosmic-ray anisotropy. The procedure also allows us to produce normalization constants that can trace minute latitudinal variations of experimental response to cosmic-ray intensity. Applying the correction of the latitudinal response function to the Tibet ASγ data, we generate a sky map of true relative intensity. As a result, we observe that the measured and corrected sky maps show significant differences in intensity and angular spectral power. Our full anisotropy sky map of true relative intensity contradicts the assumption that the latitudinal variation in longitudinally averaged flux is negligible. The result further confirms that TeV cosmic-ray anisotropy is dominated by a dipole (ℓ = 1) aligned with the interstellar magnetic field’s direction. Our results also confirm the existence of much weaker middle-scale interstellar anisotropy between ℓ = 2 and ℓ = 13.

Interplay between the non-resonant streaming instability and self-generated pressure anisotropies

ABSTRACT The non-thermal particles escaping from collision-less shocks into the surrounding medium can trigger a non-resonant streaming instability that converts parts of their drift kinetic energy into large amplitude magnetic field perturbations, and promote the confinement and acceleration of high energy cosmic rays. We present simulations of the instability using an hybrid-Particle-in-Cell approach including Monte Carlo collisions, and demonstrate that the development of the non-resonant mode is associated with important ion pressure anisotropies in the background plasma. Depending on the initial conditions, the anisotropies may act on the instability by lowering its growth and trigger secondary micro-instabilities. Introducing collisions with neutrals yield a strong reduction of the magnetic field amplification as predicted by linear fluid theory. In contrast, Coulomb collisions in fully ionized plasmas are found to mitigate the self-generated pressure anisotropies and promote the growth of the magnetic field.

Characterization of in situ cosmogenic 14CO production, retention and loss in firn and shallow ice at Summit, Greenland

Abstract. Measurements of carbon-14-containing carbon monoxide (14CO) in glacial ice are useful for studies of the past oxidative capacity of the atmosphere as well as for reconstructing the past cosmic ray flux. The 14CO abundance in glacial ice represents the combination of trapped atmospheric 14CO and in situ cosmogenic 14CO. The systematics of in situ cosmogenic 14CO production and retention in ice are not fully quantified, posing an obstacle to interpretation of ice core 14CO measurements. Here we provide the first comprehensive characterization of 14CO at an ice accumulation site (Summit, Greenland), including measurements in the ice grains of the firn matrix, firn air and bubbly ice below the firn zone. The results are interpreted with the aid of a firn gas transport model into which we implemented in situ cosmogenic 14C. We find that almost all (≈ 99.5 %) of in situ 14CO that is produced in the ice grains in firn is very rapidly (in <1 year) lost to the open porosity and from there mostly vented to the atmosphere. The timescale of this rapid loss is consistent with what is expected from gas diffusion through ice. The small fraction of in situ 14CO that initially stays in the ice grains continues to slowly leak out to the open porosity at a rate of ≈ 0.6 % yr−1. Below the firn zone we observe an increase in 14CO content with depth that is due to in situ 14CO production by deep-penetrating muons, confirming recent estimates of 14CO production rates in ice via the muon mechanisms and allowing for narrowing constraints on these production rates.

Extended gamma-ray emission from particle escape in pulsar wind nebulae. Application to HESS J1809-193 and HESS J1825-137

There is growing evidence from gamma-ray observations at high and very high energies that particle escape is a key aspect shaping the morphological properties of pulsar wind nebulae (PWNe) at various evolutionary stages. We aim to provide a simple model for the gamma-ray emission from these objects including the transport of particles across the different components of the system. We applied it to sources HESS J1809$-$193 and HESS J1825$-$137. We developed a multi-zone framework applicable to dynamically young PWNe, taking into account the diffusive escape of relativistic electron-positron pairs out of the nebula into the parent supernova remnant (SNR) and their confinement downstream of the magnetic barrier of the forward shock until an eventual release into the surrounding interstellar medium (ISM). For a wide range of turbulence properties in the nebula, the GeV-TeV inverse-Compton radiation from pairs that escaped into the remnant can be a significant if not dominant contribution to the emission from the system. It may dominate the pion-decay radiation from cosmic rays accelerated at the forward shock and advected downstream of it. In the TeV-PeV range, the contribution from particles escaped into the ISM can exceed by far that of the SNR+PWN components. Applied to HESS J1809$-$193 and HESS J1825$-$137, we found that spatially extended GeV-TeV emission components can be accounted for mostly from particles escaped into the ISM, while morphologically more compact components above $50-100$ are ascribed to the PWNe. In these two cases, the model suggests high turbulence in the nebula and a forward shock accelerating cosmic rays up to $ at most. The model provides the temporal and spectral properties of the flux of particles originally energized by the pulsar wind and ultimately released in the ISM. It can be used to constrain the transport of particles in the vicinity of pulsar-PWN-SNR systems from broadband gamma-ray observations, or in studies of the contribution of pulsar-related systems to the local positron flux.

Different spectra of cosmic ray H, He, and heavier nuclei escaping compact star clusters

ABSTRACT Cosmic ray acceleration at the termination shock of compact star clusters has recently received much attention, mainly because of the detection of gamma-ray emission from some of such astrophysical sources. Here we focus on the acceleration of nuclei at the termination shock and we investigate the role played by proton energy losses and spallation reactions of nuclei, especially downstream of the shock. We show that for a reasonable choice of the mean gas density in the cavity excavated by the cluster wind, dominated by the presence of dense clouds, the spectrum of He nuclei escaping the bubble is systematically harder than the spectrum of hydrogen, in a manner that appears to be qualitatively consistent with the observed and yet unexplained phenomenon of discrepant hardening. We also find that, in this scenario, the spallation reactions of heavier nuclei are likely to be so severe that their spectra becomes very hard and with a low normalization, meaning that it is unlikely that heavy nuclei escaping star clusters can provide a sizeable contribution to the spectrum of cosmic rays at the Earth. Limitations and implications of this scenario are discussed.

Observational Evidence for Magnetic Field Amplification in SN 1006

We report the first observational evidence for magnetic field amplification in the northeast/southwest (NE/SW) shells of supernova remnant SN 1006, one of the most promising sites of cosmic ray acceleration. In previous studies, the strength of magnetic fields in these shells was estimated to be B SED ≃ 25 μG from the spectral energy distribution, where the synchrotron emission from relativistic electrons accounted for radio to X-rays, along with the inverse Compton emission extending from the GeV to TeV energy bands. However, the analysis of broadband radio data, ranging from 1.37 to 100 GHz, indicated that the radio spectrum steepened from α 1 = 0.52 ± 0.02 to α 2 = 1.34 ± 0.21 by Δα = 0.85 ± 0.21. This is naturally interpreted as a cooling break under a strong magnetic field of B brk ≥ 2 mG. Moreover, the high-resolution MeerKAT image indicated that the width of the radio NE/SW shells was broader than that of the X-ray shell by a factor of only 3−20, as measured by Chandra. Such narrow radio shells can be naturally explained if the magnetic field responsible for the radio emissions is B R ≥ 2 mG. Assuming that the magnetic field is locally enhanced by a factor of approximately a = 100 along the NE/SW shells, we argue that the filling factor, which is the volume ratio of such a magnetically enhanced region to that of the entire shell, must be as low as approximately k = 2.5 × 10−5.

Impact of Cosmic Rays on Radiation Exposures and Scientific Activities at the Atacama Large Millimeter/Submillimeter Array (ALMA) Sites.

This study delves into the investigation of cosmic-ray radiation exposure levels for workers and their impact on the signal correlation subsystems at the Atacama Large Millimeter/submillimeter Array (ALMA) observatory sites. The analysis presents a detailed examination of secondary cosmic ray spectra and flux at the ALMA sites, encompassing the operational period from 2010 to the present day, with a particular focus on the consequences of extreme solar flares. In terms of radiation exposure for ALMA employees, the annual exposure at the highest site (AOS) reaches approximately 4.8 mSv. This value exceeds the exposure level of a typical nuclear fuel cycle worker or those working at high-altitude Antarctica stations. The exposure is approximately 2.7 times lower at the ALMA Operations Support Facility (OSF). Furthermore, the additional ambient dose equivalent resulting from solar events, while low for events similar to those observed since the 1950s, can reach up to approximately 1 mSv when considering more ancient solar events based on environmental archives. Our analysis includes radiation effects measurements in the Baseline Correlator at the AOS and, more generally, underscores the significance of employing accurate modeling and simulation techniques to assess the effects of galactic cosmic rays and extreme solar events on the integrated circuits utilized or planned in the ALMA correlation subsystem.

Mass composition of ultrahigh energy cosmic rays from distribution of their arrival directions with the Telescope Array

Mass composition of ultrahigh energy cosmic rays from distribution of their arrival directions with the Telescope Array

Isotropy of Cosmic Rays beyond 10^{20} eV Favors Their Heavy Mass Composition.

We report an estimation of the injected mass composition of ultrahigh energy cosmic rays (UHECRs) at energies higher than 10EeV. The composition is inferred from an energy-dependent sky distribution of UHECR events observed by the Telescope Array surface detector by comparing it to the Large Scale Structure of the local Universe. In the case of negligible extragalactic magnetic fields (EGMFs), the results are consistent with a relatively heavy injected composition at E∼10 EeV that becomes lighter up to E∼100 EeV, while the composition at E>100 EeV is very heavy. The latter is true even in the presence of highest experimentally allowed extragalactic magnetic fields, while the composition at lower energies can be light if a strong EGMF is present. The effect of the uncertainty in the galactic magnetic field on these results is subdominant.

Design and development of the extended-range Bonner sphere spectrometer at Ciemat

The Neutron Standards Laboratory (LPN) at CIEMAT serves as the national reference for calibrating neutron measurement devices and conducts research on measuring different neutrons fields, arising the necessity of neutron spectrometry. There are different neutron fields which are attracting interest in the recent years, for example ambient neutron spectra, generated basically by cosmic rays in the atmosphere or secondary neutrons on the surroundings of accelerators, for instance in protontherapy centres. These neutrons spectra are characterized by having an important component of high energy neutrons up to hundreds MeV or even GeV. In order to do spectrometry in this kind of neutron fields, a Bonner Sphere Spectrometer extended range was designed, constructed and tested.