Cosmic Ray Flux Research Articles

Cosmic-Ray Acceleration of Galactic Outflows in Multiphase Gas

We investigate the dynamical interaction between cosmic rays (CRs) and the multiphase interstellar medium (ISM) using numerical magnetohydrodynamic (MHD) simulations with a two-moment CR solver and TIGRESS simulations of star-forming galactic disks. We previously studied the transport of CRs within TIGRESS outputs using a “postprocessing” approach, and we now assess the effects of the MHD backreaction to CR pressure. We confirm our previous conclusion that there are three quite different regimes of CR transport in multiphase ISM gas, while also finding that simulations with “live MHD” predict a smoother CR pressure distribution. The CR pressure near the midplane is comparable to other pressure components in the gas, but the scale height of CRs is far larger. Next, with a goal of understanding the role of CRs in driving galactic outflows, we conduct a set of controlled simulations of the extraplanar region above z = 500 pc, with imposed boundary conditions flowing from the midplane into this region. We explore a range of thermal and kinematic properties for the injected thermal gas, encompassing both hot, fast-moving outflows, and cooler, slower-moving outflows. The boundary conditions for CR energy density and flux are scaled from the supernova rate in the underlying TIGRESS model. Our simulations reveal that CRs efficiently accelerate extraplanar material if the latter is mostly warm/warm-hot gas, in which CRs stream at the Alfvén speed, and the effective sound speed increases as density decreases. In contrast, CRs have very little effect on fast, hot outflows where the Alfvén speed is small, even when the injected CR momentum flux exceeds the injected MHD momentum flux.

CosmiXs: cosmic messenger spectra for indirect dark matter searches

The energy spectra of particles produced from dark matter (DM) annihilation or decay are one of the fundamental ingredients to calculate the predicted fluxes of cosmic rays and radiation searched for in indirect DM detection. We revisit the calculation of the source spectra for annihilating and decaying DM using the Vincia shower algorithm in Pythia to include QED and QCD final state radiation and diagrams for the EW corrections with massive bosons, not present in the default Pythia shower model. We take into account the spin information of the particles during the entire EW shower and the off-shell contributions from massive gauge bosons. Furthermore, we perform a dedicated tuning of the Vincia and Pythia parameters to LEP data on the production of pions, photons, and hyperons at the Z resonance and discuss the underlying uncertainties. To enable the use of our results in DM studies, we provide the tabulated source spectra for the most relevant cosmic messenger particles, namely antiprotons, positrons, γ rays and the three neutrino flavors, for all the fermionic and bosonic channels and DM masses between 5 GeV and 100 TeV, on github.

Intelligent methods for natural data analysis: application to space weather

The paper describes methods for detecting anomalies in geophysical monitoring data. This work studies a highly relevant class of problems in this area, aimed at creating methods for space weather forecasting. The negative impact of space weather anomalies on human health and practically all modern infrastructure objects requires the development of methods and the creation of effective means of detecting anomalies. Threshold wavelet filtering methods widely used for data analysis and anomaly detection allow one to obtain fairly accurate estimates using a “greed” strategy, even in the case of incomplete noise data. With this approach, the signal is estimated by isolating coherent structures. But these methods have high computational complexity, failing to provide accurate estimates when the signal-to-noise ratio is low. For such signals, we propose using adaptive probabilistic thresholds. Threshold function parameters are introduced that make it possible to estimate the variability of the process, suppress noise, and detect nonstationary features of different time-frequency structures. The paper also considers ways to combine threshold wavelet filtering with neural networks of the NARX and Autoencoder architectures. Schemes for the implementation of such approaches in the problems of detecting space weather anomalies are proposed. Using the problem of detecting ionospheric anomalies as an example, the efficiency of combining threshold wavelet filtering with the NARX network is shown. The efficiency of network sharing an Autoencoder with adaptive wavelet threshold filtering is shown in the problem of anomaly detection in cosmic ray flux intensity data.

Relationship between magnetic rigidity cutoff and chaotic behavior in cosmic ray time series using visibility graph and network analysis techniques.

Cosmic rays are highly energetic particles originating from astrophysical events outside the Solar System. In this study, we analyze the time series of cosmic ray flux measured by neutron detectors at 16 monitoring stations distributed worldwide. By applying visibility graph analysis, we explore the relationship between the magnetic rigidity cutoff (Rc) and the fractality exhibited from topology of the cosmic ray time series. Our results reveal a significant association between the magnetic rigidity cutoff and the fractality of the cosmic ray time series. Specifically, the analysis of visibility graphs and network properties demonstrates that the magnetic rigidity is inversely related to the magnetic rigidity cutoff. The identified relationship between magnetic rigidity and fractality provides insights into the chaotic nature of cosmic ray variations and their potential uses for predictability.

Cutoff Rigidities, Galactic Cosmic Ray Flux, and Heavy Ion Detections at Jupiter

AbstractA map of vertical cutoff rigidities has been calculated for galactic cosmic ray (GCR) entry into the atmosphere of Jupiter at the 1 bar pressure radius (1 RJ = 71,492 km) using the JRM33 comprehensive model of Jupiter's magnetic field (based on 32 close flybys of Jupiter by the Juno satellite) along with a particle trajectory code (Geomagnetic Cutoff Rigidity Computer Program). The map was combined with measurements of the GCR proton flux at Earth, from the BESS‐Polar ii campaign, to calculate a corresponding proton flux map at Jupiter. Additional cutoff rigidity maps were calculated for 1,000 km above the 1 bar level, and for 1.41 RJ. Furthermore, detections of heavy particles from Juno's Stellar Reference Unit were analyzed for their cutoff rigidities in multiple directions. Cutoff rigidities of 3.5–7.5 GV were found for five of the detections furthest out making them possible GCR candidates. The majority of points, located below 1.6 RJ are not likely to be GCR. Assuming instead that they are trapped particles we have calculated upper and lower limits on their equatorial pitch angles, resulting in a range from 10.1° to 27.1°, which can help constraining Jupiter's energetic radiation.

Origins of cosmic positrons and electrons

The precision measurements of the cosmic-ray positron and electron fluxes collected by the Alpha Magnetic Spectrometer on the International Space Station during first ten years of operations are presented. The positron flux exhibits complex energy dependence. It is described by the sum of a term associated with the positrons produced in the collision of cosmic rays, which dominates at low energies, and a new source term, which dominates at high energies and is associated with either dark matter or astrophysical origin. The positron source term also manifest itself in the measured electron spectrum. This is the first indication of the existence of an identical charge symmetric source term both in the positron and in the electron spectra. If confirmed with a higher level of confidence, these observations might suggest the existence of new physics.

Does or Did the Supernova Remnant Cassiopeia A Operate as a PeVatron?

For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE; E γ ≥ 100 TeV) γ-rays. In this context, the historical SNR Cassiopeia A (Cas A) is considered one of the most promising targets for UHE observations. This paper presents the observation of Cas A and its vicinity by the LHAASO KM2A detector. The exceptional sensitivity of LHAASO KM2A in the UHE band, combined with the young age of Cas A, enabled us to derive stringent model-independent limits on the energy budget of UHE protons and nuclei accelerated by Cas A at any epoch after the explosion. The results challenge the prevailing paradigm that Cas A–type SNRs are major suppliers of PeV CRs in the Milky Way.

Galactic Gamma-Ray Diffuse Emission at TeV Energies with HAWC Data

Galactic gamma-ray diffuse emission (GDE) is emitted by cosmic rays (CRs), ultra-relativistic protons, and electrons, interacting with gas and electromagnetic radiation fields in the interstellar medium. Here we present the analysis of teraelectronvolt diffuse emission from a region of the Galactic plane over the range in longitude of l ∈ [43°, 73°], using data collected with the High Altitude Water Cherenkov (HAWC) detector. Spectral, longitudinal, and latitudinal distributions of the teraelectronvolt diffuse emission are shown. The radiation spectrum is compatible with the spectrum of the emission arising from a CR population with an index similar to that of the observed CRs. When comparing with the DRAGON base model, the HAWC GDE flux is higher by about a factor of 2. Unresolved sources such as pulsar wind nebulae and teraelectronvolt halos could explain the excess emission. Finally, deviations of the Galactic CR flux from the locally measured CR flux may additionally explain the difference between the predicted and measured diffuse fluxes.

Time calibration and synchronization of the scintillation light detection system in ICARUS-T600

The ICARUS-T600 Liquid Argon (LAr) Time Projection Chamber (TPC) is presently taking data in the Short Baseline Neutrino (SBN) program at Fermilab (U.S.A.) to search for a possible LSND-like sterile neutrino signal at Δm 2 ≈ 1 eV2 with the Booster Neutrino Beam (BNB). A light detection system, based on 360 large area Photo-Multiplier Tubes (PMTs), has been realized for ICARUS-T600 to detect VUV photons produced by the passage of ionizing particles in LAr. This system is fundamental for the TPC operation, providing an efficient trigger and contributing to the 3D reconstruction of events. Moreover, since the detector is exposed to a huge flux of cosmic rays due to its shallow depths installation, the light detection system allows for the time reconstruction of events, contributing to the identification and to the selection of genuine neutrino interactions. The correct time reconstruction of events requires the precise knowledge of the delay of each PMT channel and a good synchronization of recording electronics, this last based on fast sampling digitizers. To achieve a time resolution better than 1 ns, we perform three consecutive timing corrections deployed at different stages of the optical data flow. Results demonstrate the capability of the ICARUS-T600 light detection system to allow a precise reconstruction of the temporal evolution of each event occurring in the detector and the association of neutrino events with the bunched structure of BNB.

Constraining (anti)nuclei measurements relevant for astrophysics with ALICE

Antinuclei can be produced in space either by collisions of highenergy cosmic rays with the interstellar medium or from the annihilation of dark matter particles stemming into standard model particles. High-energy hadronic collisions at accelerators create a suitable environment for producing light (anti)nuclei. Hence, studying the production of (anti)nuclei in pp collisions at the LHC can provide crucial insights into the production mechanisms of nuclear states in our Universe. Recent measurements of the production of (anti)nuclei in and out of jets, and as a function of rapidity in pp collisions at √S = 13 TeV have been carried out with ALICE. The latter allow for the extrapolation of the nuclear production models at forward rapidity, region of interest for indirect searches of dark matter. Recent results on the annihilation crosssection of antinuclei are also discussed in the context of astrophysical measurements of cosmic ray flux. Such information is essential to study the different sources of antinuclei in our Universe and to interpret any future measurement of antinuclei in space.

Solar influences on the Earth’s atmosphere: solved and unsolved questions

Solar influences on the Earth’s atmosphere: solved and unsolved questions

H i in Molecular Clouds: Irradiation by FUV Plus Cosmic Rays

We extend the analytic theory presented by Sternberg et al. and Bialy & Sternberg for the production of atomic hydrogen (H i) via far-ultraviolet (FUV) photodissociation at the boundaries of dense interstellar molecular (H2) clouds to also include the effects of penetrating (low-energy) cosmic rays (CRs) for the growth of the total H i column densities. We compute the steady-state abundances of the H i and H2 in one-dimensional gas slabs in which the FUV photodissociation rates are reduced by depth-dependent H2 self-shielding and dust absorption and the CR ionization rates are either constant or reduced by transport effects. The solutions for the H i and H2 density profiles and the integrated H i columns depend primarily on the ratios I UV/Rn and ζ/Rn, where I UV is the intensity of the photodissociating FUV field, ζ is the H2 CR ionization rate, n is the hydrogen gas density, and R is the dust surface H2 formation rate coefficient. We present computations for a wide range of FUV field strengths, CR ionization rates, and dust-to-gas ratios. We develop analytic expressions for the growth of the H i column densities. For Galactic giant molecular clouds (GMCs) with multiphased (warm/cold) H i envelopes, the interior CR zones will dominate the production of the H i only if s−1, where M GMC is the GMC mass, and including attenuation of the CR fluxes. For most Galactic GMCs and conditions, FUV photodissociation dominates over CR ionization for the production of the H i column densities. Furthermore, the CRs do not affect the H i-to-H2 transition points.

Solar Activity, Weather, and Climate: The Elusive Connection

Abstract The search for an explanation of correlations of weather and climate with solar activity has uncovered some subtle effects of the varying solar inputs, which include irradiance, energetic particles, electric fields induced by the solar wind, and the modulation of the galactic cosmic ray flux. Relevant phenomena have been identified in stratospheric chemistry and dynamics and in aerosols, atmospheric electricity, and clouds. The effects are small and appear intermittently in statistical analyses of observations. Correlations and theory permit responses to several solar inputs effective on the decadal time scales, making it difficult to distinguish them observationally. No mechanism for tropospheric decadal change can be considered to be definitely established, but for the day-to-day time scale there is clear observational evidence for clouds responding to solar wind–induced current flow (JZ) in the global electric circuit, due to changes in ionospheric potential, and independently from Forbush decreases of the cosmic ray flux. The responses to JZ also correlate separately with ionospheric potential changes due to changes in day-to-day thunderstorm activity. The identified mechanism is for electric charge effects on in-cloud scavenging, and this implies a decadal response, in view of decadal changes in JZ. However, a complete and quantitative model of the inferred electrically induced changes in cloud microphysics is not yet available. The failures, successes, and controversies of the search illustrate the somewhat chaotic process of scientific discovery.

The Effect of Magnetic Field Line Topology on ICME-related GCR Modulation

The large-scale magnetic structure of interplanetary coronal mass ejections (ICMEs) has been shown to affect the galactic cosmic-ray (GCR) flux measured in situ by spacecraft, causing temporary decreases known as Forbush decreases (FDs). In some ICMEs, the magnetic ejecta exhibits a magnetic flux rope (FR) structure; the strong magnetic field strength and closed field line geometry of such ICME FRs have been proposed to act as a shield to GCR transport. In this study, we identify four ICMEs near Earth that drove FDs with similar mean magnetic field strengths (20–25 nT), two ICMEs with more typical mean speeds (∼400 km s−1), and two fast (∼750 km s−1) ICMEs. Within each speed pairing, we identify an ICME that exhibited an open magnetic field line topology and compare its effect on the GCR flux to that which exhibited a mostly closed topology. We investigate the different mechanisms that contribute to the resulting ICME-related FDs and their recovery and determine which properties, if any, play a more important role than others in driving FDs. We find that much of the GCR response to the ICME events in this study is independent of the open or closed magnetic field line topology of the FR and that features such as the fluctuations in speed, magnetic field structure, and expansion within the FR may play more of a role in determining the smaller-scale structure of the FD profile.

Monte Carlo Simulation and measurement of Calibration Neutron Monitor count rate dependence on proximity to water

Due to their global availability, neutron monitors play a crucial role in measuring time variations in the Galactic cosmic ray flux. A portable calibration neutron monitor (CalMon) is useful for intercalibrating various neutron monitors to ensure accurate measurements. A common technique to ensure that the calibration is done in a consistent environment is to place the CalMon at some height above a wide container (such as a portable swimming pool) filled with water. This study investigates the impact of CalMon height and water depth on the count rate ratio relative to a standard 18NM64 count rate recorded nearby (CalMon/18NM64). We compare simulated data from the FLUKA Monte Carlo package to experimental data from [1] to demonstrate the statistical accuracy of our simulation. Using the simulation results, we then extend the study of the proximity-to-water effect on the counting rate. In this work, we present a preliminary empirical model by analyzing the CalMon/18NM64 as a function of CalMon to water distance. Overall, our study enhances understanding of the response of calibration monitors (now often called “mini-neutron monitors”) operated in various locations worldwide, and validates the Monte Carlo techniques used to model the response of the global neutron monitor network.

Analysis of Cosmic Ray Fluxes at Different Stations during Geomagnetic Storms using Wavelet Based Approaches: Continuous Wavelet Transform and Multi-Resolution Analysis

Analysis of Cosmic Ray Fluxes at Different Stations during Geomagnetic Storms using Wavelet Based Approaches: Continuous Wavelet Transform and Multi-Resolution Analysis

Modification of the dipole in arrival directions of ultra-high-energy cosmic rays due to the Galactic magnetic field

The direction and magnitude of the dipole anisotropy of ultra-high-energy cosmic rays with energies above 8 EeV observed by the Pierre Auger Observatory indicate their extragalactic origin. The observed dipole on Earth does not necessarily need to correspond to the anisotropy of the extragalactic cosmic-ray flux due to the effects of propagation in the Galactic magnetic field. We estimate the size of these effects via numerical simulations using the CRPropa 3 package. The Jansson-Farrar and Terral-Ferrière models of the Galactic magnetic field are used to propagate particles within the Galaxy. We identify allowed directions and amplitudes of the dipole outside the Galaxy that are compatible with the measured features of the dipole on Earth for various mass composition scenarios at the 68% and 95% confidence level.

Research on rapid imaging with cosmic ray muon scattering tomography

Cosmic ray muons tomography is a non-destructive imaging technique that uses the natural radiation of cosmic ray muons to create tomographic images of objects. This article presents a novel imaging algorithm that effectively utilizes experimental measurement data to achieve rapid and clear imaging of cosmic ray muons. A clear image can be obtained with only 20 min of measurement time and approximately 200 effective muons. However, the current detection flux is only about 0.044 cm^{-2} min^{-1}, which is significantly lower than the natural cosmic ray flux of about 1 cm^{-2} min^{-1}.

Application of Geant4 simulation in measurement of cosmic-ray muon flux and studies of muon-induced background

Application of Geant4 simulation in measurement of cosmic-ray muon flux and studies of muon-induced background

The Maximum Energy of Shock-accelerated Cosmic Rays

Identifying the accelerators of Galactic cosmic ray (CR) protons with energies up to a few PeV (1015 eV) remains a theoretical and observational challenge. Supernova remnants (SNRs) represent strong candidates because they provide sufficient energetics to reproduce the CR flux observed at Earth. However, it remains unclear whether they can accelerate particles to PeV energies, particularly after the very early stages of their evolution. This uncertainty has prompted searches for other source classes and necessitates comprehensive theoretical modeling of the maximum proton energy, , accelerated by an arbitrary shock. While analytic estimates of have been put forward in the literature, they do not fully account for the complex interplay between particle acceleration, magnetic field amplification, and shock evolution. This paper uses a multizone, semianalytic model of particle acceleration based on kinetic simulations to place constraints on for a wide range of astrophysical shocks. In particular, we develop relationships between , shock velocity, size, and ambient medium. We find that SNRs can only accelerate PeV particles under a select set of circumstances, namely, if the shock velocity exceeds ∼104 km s−1 and escaping particles drive magnetic field amplification. However, older and slower SNRs may still produce observational signatures of PeV particles due to populations accelerated when the shock was younger. Our results serve as a reference for modelers seeking to quickly produce a self-consistent estimate of the maximum energy accelerated by an arbitrary astrophysical shock. 1 1 Presented as a thesis to the Department of Astronomy and Astrophysics, The University of Chicago, in partial fulfillment of the requirements for a Ph.D. degree.