Search

Extensive air showers (EAS), cascades of electromagnetic radiation and ionized particles produced when primary cosmic rays (CRs) interact with atmospheric nuclei, generate a plethora of secondary particles, including electrons, muons, alpha particles, X-rays, and neutrons. These phenomena, pivotal in the exploration of high-energy CR interactions, have prompted advancements in simulating the lateral distribution function (LDF) of secondary particles. Utilizing AIR-shower extended simulations (AIRES) software (version 19.04.10), this study presents simulations of EAS to elucidate the impact on the LDF at the high energy thresholds of 1016, 1018 and 1019 eV. The focus is on secondary muons and electrons, considering the influence of primary protons and iron nuclei across zenith angles of 0°, 20°, and 40°. Through the application of an exponential function model, novel coefficients were derived, reflecting the variation of the LDF relative to the distance from the shower core. These coefficients are instrumental in astrophysics and particle physics, offering a nuanced understanding of CR interactions with earth's atmosphere. By delineating the spatial particle distribution during interactions, insights into the shower's lateral extent and configuration are gleaned, facilitating the deduction of the initiating CR's energy. Such knowledge is crucial for CR spectrum analysis. The investigation of electron and muon LDFs not only advances fundamental particle physics but also holds implications for space exploration, radiation protection, and other domains. [ABSTRACT FROM AUTHOR]

To catalyze transformative advancements in High-energy Physics in the Atmosphere (HEPA), a comprehensive understanding of particle fluxes, electric fields, and lightning occurrences across atmospheric layers is imperative. This paper elucidates the instrumentation and capabilities of the Aragats Space-Environmental Center (ASEC), which encompasses measurement tools for various cosmic ray species, near-surface electric fields, and lightning events integrated across high-mountain research station at slopes of Mt. Aragats and the highest mountains of Eastern Europe and Germany. Through these measurements, we aim to elucidate models of particle acceleration mechanisms and the charge distribution within the lower atmosphere. We introduce an Advanced Data Extraction Infrastructure (ADEI) integrated with sophisticated statistical analysis tools to facilitate rapid access to this wealth of data. Despite the significance of these atmospheric processes, the intricate interplay between thundercloud electrification, lightning activity, wideband radio emissions, and particle fluxes remains poorly understood. A particularly compelling avenue of inquiry lies in exploring the relationship between high-energy atmospheric phenomena, intracloud electric fields, and lightning initiation. Furthermore, investigations into accelerated structures within geospace plasmas hold promise for shedding light on particle acceleration processes, potentially extending to higher energies within analogous structures in cosmic plasmas. This paper also examines practical methodologies for extracting meaningful physical insights from temporal datasets, such as correlating surges in particle flux intensity with variations in near-surface electric field strength and precipitation patterns. Additionally, we explore the utility of wideband field and interferometer antenna signals in this context, offering valuable avenues for further research and analysis. Through these endeavors, we aim to deepen our understanding of high-energy atmospheric processes and their broader implications for terrestrial and cosmic phenomena. [ABSTRACT FROM AUTHOR]

The outflow structures of gamma-ray bursts (GRBs) can provide insights into the origins and radiation mechanisms of these cosmic explosions. We systematically study the GRB outflow structures by modelling their afterglow light curves and check if the structures of long gamma-ray bursts (LGRBs) and short gamma-ray bursts (SGRBs) are different. The sample consists of Swift -XRT afterglows with sufficient coverage and known redshift, which includes 195 well-fit LGRBs and 13 well-fit SGRBs. The model we use is a two-parameter 'boosted fireball' model, which consists of a family of outflows, with shapes varying smoothly from a quasi-spherical outflow to a highly collimated jet. We use the jetfit package to fit afterglow light curves and obtain the jet parameters. We find that there are no statistical differences in the distributions of jet parameters between LGRBs and SGRBs by performing K–S test and 74 per cent of the ratios of the observer angle to jet opening angle are in the range of 0.2 to 1. Our analysis indicates that the majority of GRB afterglows are viewed off-axis and there has no statistical difference between LGRBs and SGRBs. We also find that both the LGRBs and SGRBs exhibit two similar correlations: the jet opening angle is positively correlated with the observer angle, with the correlation coefficient 0.61 for LGRBs and 0.63 for SGRBs; the circumburst density is inversely correlated with the explosion energy with the correlation coefficient −0.89 for LGRBs and −0.69 for SGRBs. Our results suggest that the outflow structures are similar for the LGRBs and SGRBs. [ABSTRACT FROM AUTHOR]

Introduction: In this article, we revisit the pitch-angle scattering equation describing the propagation of energetic particles through magnetized plasma. In this case, solar energetic particles and cosmic rays interact with magnetohydrodynamic turbulence and experience stochastic changes in the pitch-angle. Since this happens over an extended period of time, a pitch-angle isotropization process occurs, leading to parallel spatial diffusion. This process is described well by the pitch-angle scattering equation. However, the latter equation is difficult to solve analytically even when considering special cases for the scattering coefficient. Methods: In the past, a so-called subspace approximation was proposed, which has important applications in the theory of perpendicular diffusion. Alternatively, an approach based on the telegraph equation (also known as telegrapher's equation) has been developed. We show that two-dimensional subspace approximation and the description based on the telegraph equation are equivalent. However, it is also shown that the obtained distribution functions contain artifacts and inaccuracies that cannot be found in the numerical solution to the problem. Therefore, an N-dimensional subspace approximation is proposed corresponding to a semi-analytical/semi-numerical approach. This is a useful alternative compared to standard numerical solvers. Results and Discussion: Depending on the application, the N-dimensional subspace approximation can be orders of magnitude faster. Furthermore, the method can easily be modified so that it can be used for any pitch-angle scattering equation. [ABSTRACT FROM AUTHOR]

Electron densities in different locations of our galaxy are obtained in pulsar astronomy by dividing the dispersion measure (DM) by the distance of the pulsar to Earth. The properties of the interstellar plasma are related to its heating by cosmic rays. Following the present analysis, DM's measurements are obtained with different properties of the corresponding plasmas at different temperatures. (1) For relatively low temperatures (around and below 10,000 (K)), the state of molecular, atomic and ionized hydrogen is analyzed by the interstellar medium model with partially ionized plasma. In this region, various spectroscopic effects are obtained. (2) For temperatures above this limit, the interstellar gas is found to be a completely ionized medium, and this plasma is defined as warm ionized medium (WIM) plasma. We show in the first part of the paper that this plasma is transparent, as obtained from the solution of Saha's equation. The index of refraction of the WIM plasma is real, but the plasma can be observed by dispersion measurements. (3) For very high temperatures (around 1 million (K)), the plasma is defined as hot ionized medium, where X-rays are obtained. We concentrate in the present work on the analysis of WIM plasma. We calculate the mass densities of this plasma and compare it with dark matter mass densities, which are found to be larger. But some factors which may reduce this difference are discussed. [ABSTRACT FROM AUTHOR]

The accuracy of cosmic ray observations by the Large High Altitude Air Shower Observatory Wide Field-of-View Cherenkov/Fluorescence Telescope Array (LHAASO-WFCTA) is influenced by variations in aerosols in the atmosphere. The solar photometer (CE318-T) is extensively utilized within the Aerosol Robotic Network as a highly precise and reliable instrument for aerosol measurements. With this CE318-T 23, 254 sets of valid data samples over 394 days from October 2020 to October 2022 at the LHAASO site were obtained. Data analysis revealed that the baseline Aerosol Optical Depth (AOD) and Ångström Exponent (AE) at 440–870 nm (AE440–870nm) of the aerosols were calculated to be 0.03 and 1.07, respectively, suggesting that the LHAASO site is among the most pristine regions on Earth. The seasonality of the mean AOD is in the order of spring > summer > autumn = winter. The monthly average maximum of AOD440nm occurred in April (0.11 ± 0.05) and the minimum was in December (0.03 ± 0.01). The monthly average of AE440–870nm exhibited slight variations. The seasonal characterization of aerosol types indicated that background aerosol predominated in autumn and winter, which is the optimal period for the absolute calibration of the WFCTA. Additionally, the diurnal daytime variations of AOD and AE across the four seasons are presented. Our analysis also indicates that the potential origins of aerosol over the LHAASO in four seasons were different and the atmospheric aerosols with higher AOD probably originate mainly from Northern Myanmar and Northeast India regions. These results are presented for the first time, providing a detailed analysis of aerosol seasonality and origins, which have not been thoroughly documented before in this region, also enriching the valuable materials on aerosol observation in the Hengduan Mountains and Tibetan Plateau. [ABSTRACT FROM AUTHOR]

Galaxies obey a number of empirical correlations between their radio, γ-ray, and infrared emission, but the physical origins of these correlations remain uncertain. Here, we use the CONGRuENTS model for broad-band non-thermal emission from star-forming galaxies, which self-consistently calculates energy-dependent transport and non-thermal emission from cosmic ray hadrons and leptons, to predict radio and γ-ray emission for a synthetic galaxy population with properties drawn from a large deep-field survey. We show that our synthetic galaxies reproduce observed relations such as the far infrared (FIR)–radio correlation, the FIR–γ correlation, and the distribution of radio spectral indices, and we use the model to explain the physical origins of these relations. Our results show that the FIR–radio correlation arises because the amount of cosmic ray electron power ultimately radiated as synchrotron emission varies only weakly with galaxy star formation rate as a result of the constraints imposed on gas properties by hydrostatic balance and turbulent dynamo action; the same physics dictates the extent of proton calorimetry in different galaxies, and thus sets the FIR–γ–ray correlation. We further show that galactic radio spectral indices result primarily from competition between thermal free–free emission and energy-dependent loss of cosmic ray electrons to bremsstrahlung and escape into galactic haloes, with shaping of the spectrum by inverse Compton, synchrotron, and ionization processes typically playing a subdominant role. In addition to explaining existing observations, we use our analysis to predict a heretofore unseen correlation between the curvature of galaxies' radio spectra and their pion-driven γ-ray emission, a prediction that will be testable with upcoming facilities. [ABSTRACT FROM AUTHOR]

Radiolysis of α-pinene by 61.3 MeV 84Kr15 + ions was analysed with the scope to simulate the effects of heavy ion cosmic ray bombardment on chiral molecules in the interstellar medium. The α-pinene ice samples were irradiated at 10 K and their chemical evolution was monitored by mid-infrared Fourier transform (FTIR) spectroscopy to characterize the reaction products and to determine the extent of racemization. The integrated band strengths have been obtained for all the neutral α-pinene vibrational bands using the experimental band integrated absorbances and the theoretical absolute intensities calculated along the column densities. In the current heavy ion bombardment experiments, small molecules were formed and the precursor, α-pinene, was destroyed instead of being racemized. Twelve hydrocarbons were produced (final fluence of 2.0 × 1012 ions cm−2): methane (CH4), acetylene (C2H2), ethylene (C2H4), propylene (C3H6), propane (C3H8), n-butane (C4H10), butene (C4H8), propyne (C3H4), benzene (C6H6), ethane (C2H6), vinylacetylene (C4H4), and 2-methyl-1,3-butadiene or isoprene (C5H8). The highest formation cross-section (∼ 40 × 10−15 cm2) was observed for the C3H4 and the lowest was for C3H8 (∼ 3 × 10−15 cm2). The radiochemical yields for these molecules follow the same trends as those of their cross-sections. The atom budget calculation confirms that all the expected products have been generated during the radiolysis and supports the conclusion that the proposed A values are accurate. The α-pinene sputtering yield for this ion beam was found to be Y0  = 1.84 × 106 molecules per impact. [ABSTRACT FROM AUTHOR]

Meter-scale, submicron-thick lightsail spacecraft, propelled to relativistic velocities via photon pressure using high-power density laser radiation, offer a potentially new route to space exploration within and beyond the solar system, posing substantial challenges for materials science and engineering. We analyze the structural and photonic design of flexible lightsails by developing a mesh-based multiphysics simulator based on linear elastic theory. We observe spin-stabilized flexible lightsail shapes and designs that are immune to shape collapse during acceleration and exhibit beam-riding stability despite deformations caused by photon pressure and thermal expansion. Excitingly, nanophotonic lightsails based on planar silicon nitride membranes patterned with suitable optical metagratings exhibit both mechanically and dynamically stable propulsion along the pump laser axis. These advances suggest that laser-driven acceleration of membrane-like lightsails to the relativistic speeds needed to access interstellar distances is conceptually feasible, and that their fabrication could be achieved by scaling up modern microfabrication technology. Ultrathin laser-driven lightsails represent a unique vision for interstellar space exploration. Here, the authors show how spinning flexible membranes can be both shape- and trajectory-stable with multiphysics structural and nanophotonic engineering. [ABSTRACT FROM AUTHOR]

The Moon is a unique natural laboratory for the study of the deep space plasma and energetic particles environment. During more than 3/4 of its orbit around the Earth it is exposed to the solar wind. Being an unmagnetized body and lacking a substantial atmosphere, solar wind and solar energetic particles bombard the Moon's surface, interacting with the lunar regolith and the tenuous lunar exosphere. Energetic particles arriving at the Moon's surface can be absorbed, or scattered, or can remove another particle from the lunar regolith by sputtering or desorption. A similar phenomenon occurs also with the galactic cosmic rays, which have fluxes and energy spectra representative of interplanetary space. During the remaining part of its orbit the Moon crosses the tail of the terrestrial magnetosphere. It then provides the opportunity to study in-situ the terrestrial magnetotail plasma environment as well as atmospheric escape from the Earth's ionosphere, in the form of heavy ions accelerated and streaming downtail. The lunar environment is thus a unique natural laboratory for analysing the interaction of the solar wind, the cosmic rays and the Earth's magnetosphere with the surface, the immediate subsurface, and the surface-bounded exosphere of an unmagnetized planetary body. This article is part of a discussion meeting issue 'Astronomy from the Moon: the next decades (part 2)'. [ABSTRACT FROM AUTHOR]

Odd radio circles (ORCs) are a newly discovered class of extended faint radio sources of unknown origin. We report the first detection of diffuse X-ray gas at the location of a low-redshift ORC (z = 0.046) known as Cloverleaf ORC. This observation was performed with the XMM-Newton X-ray telescope. The physical extent of the diffuse X-ray emission corresponds to a region of approximately 230 kpc by 160 kpc, lying perpendicular to the radio emission detected by ASKAP. The X-ray spectrum shows characteristics of thermal multiphase gas with temperatures of 1.10 ± 0.08 keV and 0.22 ± 0.01 keV and a central density of (4.9 ± 0.6)×10−4 cm−3, indicating that the Cloverleaf ORC resides in a low-mass galaxy group. Using X-ray observations, with hydrostatic equilibrium and isothermal assumptions, we measure the galaxy group to have a gas mass and a total mass of (7.7 ± 0.8)×1011 M⊙ and 2.6 ± 0.3 × 1013 M⊙ within the overdensity radius R500. The presence of a high-velocity subgroup identified in optical data, the orientation of the brightest cluster galaxy, the disturbed morphologies of galaxies toward the east of the Cloverleaf ORC, and the irregular morphology of the X-ray emission suggest that this system is undergoing a galaxy group merger. The radio power of the ORC could be explained by the shock reacceleration of fossil cosmic rays generated by a previous episode of black hole activity in the central active galactic nucleus. [ABSTRACT FROM AUTHOR]

A thorough inspection of known Galactic supernova remnants (SNRs) along the Galactic plane with SRG/eROSITA yielded the detection of the X-ray counterpart of the SNR G279.0+01.1. The SNR is located just 1.5° above the Galactic plane. Its X-ray emission emerges as an incomplete, partial shell of ~3° angular extension. It is strongly correlated to the fragmented shell-type morphology of its radio continuum emission. The X-ray spatial morphology of the SNR seems to be linked to the presence of dust clouds in the surroundings. The majority of its X-ray emission is soft (exhibiting strong O, Ne, and Mg lines), and it occurs in a narrow range of energies between 0.3 and 1.5 keV. Above 2.0 keV the remnant remains undetected. The remnant's X-ray spectrum is purely of a thermal nature. Constraining the X-ray absorption column to values which are consistent with optical extinction data from the remnant's location favors nonequilibrium over equilibrium models. A nonequilibrium two-temperature plasma model of kT ~ 0.3 keV and kT ~ 0.6 keV, as well as an absorption column density of NH ~ 0.3 cm−2 describe the spectrum of the entire remnant well. Significant temperature variations across the remnant have been detected. Employing 14.5 yr of Fermi-LAT data, we carried out a comprehensive study of the extended giga-electronvolt source 4FGL J1000.0-5312e. By refining and properly modeling the giga-electronvolt excess originating from the location of the remnant, we conclude that the emission is likely related to the remnant itself rather than being colocated by chance. The remnant's properties as determined by the X-ray spectra are consistent with the ~2.5 kpc distance estimates from the literature, which implies a source diameter of ~140 pc and an old age of >7 × 105 yr. However, if the source is associated with any of the pulsars previously considered to be associated with the SNR, then the updated nearby pulsar distance estimates from the YMW16 electron density model rather place the SNR at a distance of ~0.4 kpc. This would correspond to a ~20 pc linear size and a younger age of 104− < 7 × 105 yr, which would be more in line with the nonequilibrium state of the plasma. [ABSTRACT FROM AUTHOR]

Recent observations have revealed detailed structures of radio relics across a wide range of frequencies. In this work, we performed three-dimensional magnetohydrodynamical (3D MHD) simulations of merger shocks propagating through a turbulent magnetized intracluster medium. We then employed on-the-fly Lagrangian particles to explore the physical processes behind the origination of radio substructures and their appearance in high and low-frequency observations. We employed two cosmic-ray (CR) electron acceleration models, with a fresh injection of electrons from the thermal pool and the re-acceleration of mildly relativistic electrons. We used the relative surface brightness fluctuations, δSν, to define a "degree of patchiness." First, we found that patchiness is produced if the shock's surface has a distribution of Mach numbers, rather than a single Mach number. Second, radio relics appear patchier if the Mach number distribution consists of a large percentage of low Mach numbers (ℳ ≲ 2.5). Furthermore, as the frequency increases, the patchiness also becomes larger. Nevertheless, if radio relics are patchy at high frequencies (e.g., 18.6 GHz), they necessarily will also be patchy at low frequencies (e.g., 150 MHz). Then, to produce noticeable differences in the patchiness at low and high frequencies, the shock front should have a Mach number spread of σℳ ≳ 0.3 − 0.4. Finally, the extent of the patchiness depends on the Mach number distribution as well as the CR acceleration model. We propose δSν as a potential tool for extracting merger shock properties and information on particle acceleration processes at shocks in radio observations. [ABSTRACT FROM AUTHOR]

Context. Hadronic interactions between cosmic rays (CRs) and interstellar gas have been probed in γ rays across the Galaxy. A fairly uniform CR distribution is observed up to a few hundred parsecs from the Sun, except in the Eridu cloud, which shows an unexplained 30–50% deficit in GeV to TeV CR flux. Aims. To explore the origin of this deficit, we studied the Reticulum cloud, which shares notable traits with Eridu: a comparable distance in the low-density region of the Local Valley and a filamentary structure of atomic hydrogen extending along a bundle of ordered magnetic-field lines that are steeply inclined to the Galactic plane. Methods. We measured the γ-ray emissivity per gas nucleon in the Reticulum cloud in the 0.16–63 GeV energy band using 14 years of Fermi-LAT data. We also derived interstellar properties that are important for CR propagation in both the Eridu and Reticulum clouds, at the same parsec scale. Results. The γ-ray emissivity in the Reticulum cloud is fully consistent with the average spectrum measured in the solar neighbourhood, but this emissivity, and therefore the CR flux, is 1.57 ± 0.09 times larger than in Eridu across the whole energy band. The difference cannot be attributed to uncertainties in gas mass. Nevertheless, we find that the two clouds are similar in many respects: both have magnetic-field strengths of a few micro-Gauss in the plane of the sky; both are in approximate equilibrium between magnetic and thermal pressures; they have similar turbulent velocities and sonic Mach numbers; and both show magnetic-field regularity with a dispersion in orientation lower than 10°–15° over large zones. The gas in Reticulum is colder and denser than in Eridu, but we find similar parallel diffusion coefficients around a few times 1028 cm2 s−1 in both clouds if CRs above 1 GV in rigidity diffuse on resonant, self-excited Alfvén waves that are damped by ion-neutral interactions. Conclusions. The loss of CRs in Eridu remains unexplained, but these two clouds provide important test cases to further study how magnetic turbulence, line tangling, and ion-neutral damping regulate CR diffusion in the dominant gas phase of the interstellar medium. [ABSTRACT FROM AUTHOR]

Context. Starburst galaxies are undergoing intense episodes of star formation. In these galaxies, gas is ejected into the surrounding environment through winds created by the effect of hot stars and supernova explosions. When interacting with the intergalactic medium, these winds can produce strong shocks capable of accelerating cosmic rays. The radiation from these cosmic rays mainly occurs in radio and gamma rays. The radio halo can be characterized using the scale height, which is an important parameter for understanding cosmic ray acceleration and transport. Aims. We searched for the presence of radio halos in a sample of edge-on starburst galaxies gathered from the MeerKAT 1.28 GHz Atlas of Southern Sources in the IRAS Revised Bright Galaxy Sample. The investigation of how the radio halos relate to the global properties of the galaxies can shed light on the understanding of the halo origin and the underlying cosmic ray population. Methods. We selected a sample of 25 galaxies with inclinations i > 80° from the original sample and modeled their disk and halo contributions. We determined the scale heights and the radio luminosity of the halos when detected. Results. We have detected and characterized 11 radio halos from a sample of 25 edge-on galaxies. Seven of them are reported here for the first time. The average radio scale height is ∼1 kpc. We found that the halo scale heights increase linearly with the radio diameters and this relation does not depend on the star formation rate. All galaxies in our sample follow the radio-infrared relation with a q parameter value of 2.5 ± 0.1. The halo luminosity linearly increases with the infrared luminosity and star formation rate. Conclusions. The dependence of the halo luminosity on the star formation rate and the infrared luminosity supports the hypothesis that the radio halos are the result of synchrotron radiation produced by relativistic electrons and points toward the fact that the star formation activity plays a crucial role in halo creation. The average scale height of 1 kpc implies a dynamical range of 4 Myr, several orders of magnitude greater than the synchrotron losses for electrons of 10 TeV. This suggests that some process must exist to reaccelerate cosmic rays in the halo if gamma-ray emission of a leptonic origin is detected from the halo. According to the relation between the radio and gamma-ray luminosities, we found that NGC 4666 is a potential gamma-ray source for future observations. [ABSTRACT FROM AUTHOR]

Context. Supernovae remnants (SNRs) represent a major source of feedback from stars on the interstellar medium of galaxies. During the latest stage of supernova explosions (which lasts 10–100 kyr), shock waves produced by the initial blast modify the chemistry of gas and dust, inject kinetic energy in the surroundings, and may alter star formation characteristics. Simultaneously, γ-ray emission is generated by the interaction between the ambient medium and cosmic rays, in particular those locally accelerated in the early stages of the explosion. Aims. We aim to estimate the total molecular mass, local density, and total column density of H2 and the temperature structure in a shocked clump interacting with the supernova remnant IC443 located in a region where cosmic rays interact with the interstellar medium. Measuring the mass of the dense and neutral component of the medium is a prerequisite to understanding the chemistry, energetics, and GeV to TeV γ-ray emission. Methods. Assuming that the emission of H2 pure rotational lines is produced by a collection of gas layers with variable temperature, we compared Spitzer/IRS emission maps for the ν = 0–0 S(0) to S(7) lines with a thermal admixture model. Our description is based on a power-law distribution of thermalized components with temperatures varying between Tmin = 25 K and Tmax = 1500 K. Results. Our thermal admixture model allows the level populations of H2 to be described by a power-law distribution dN = ΛT−ΓdT, with Γ ~ 2.2−4.7. We measured a total mass MH2 = 220−80+110M⊙ across the Spitzer/IRS field of observations. Conclusions. Our analysis shows that an estimate of the cold molecular gas temperature is paramount to accurately constraining the H2 mass, although the mass remains affected by significant uncertainties due to the assumptions on the gas temperature distribution. [ABSTRACT FROM AUTHOR]

Context. Cosmic rays (CRs) heavily impact the chemistry and physics of cold and dense star-forming regions. However, the characterisation of their ionisation rate continues to pose a challenge from the observational point of view. Aims. In the past, a few analytical formulas have been proposed to infer the cosmic-ray ionisation rate, ζ2, from molecular line observations. These have been derived from the chemical kinetics of the involved species, but they have not yet been validated using synthetic data processed with a standard observative pipeline. In this work, we aim to bridge this gap. Methods. We performed a radiative transfer on a set of three-dimensional magneto-hydrodynamical simulations of prestellar cores, exploring different initial ζ2, evolutionary stages, types of radiative transfer (for instance assuming local-thermodynamic-equilibrium conditions), and telescope responses. We then computed the column densities of the involved tracers to determine ζ2, employing a recently proposed method based on the detection of H2D+. We compared this approach with a previous method, based on more common tracers. Both approaches are commonly used. Results. Our results confirm that the equation based on the detection of H2D+ accurately retrieves the actual ζ2 within a factor of two to three in the physical conditions explored in our tests. Since we have also explored a non-local thermodynamic equilibrium (non-LTE) radiative transfer, this work indirectly offers insights into the excitation temperatures of common transitions at moderate volume densities (n ≈ 105 cm−3). We also performed a few tests using a previous methodology that is independent of H2D+, which overestimates the actual ζ2 by at least two orders of magnitude. We considered a new derivation of this method, however, we found that it still leads to high over-estimations. Conclusions. The method based on H2D+ is further validated in this work and demonstrates a reliable method for estimating ζ2 in cold and dense gas. On the contrary, the former analytical equation, as already pointed out by its authors, has no global domain of application. Thus, we find that it ought to be employed with caution. [ABSTRACT FROM AUTHOR]

We study here a novel application of Kim and Lee charged wormholes assuming them to be dark halo objects playing the role of lenses in the Galactic microlensing with source stars belonging to the Galactic Bulge and the Large Magellanic Cloud. First, we observe that both the backreacted scalar (α ) and electrically (Q) charged wormholes have the same zero ADM mass as has the background Ellis–Bronnikov wormhole having a special equation of state parameter γ = - 1 . In particular, we argue that, for α ≠ 0 , the solution formally resembles, but can at best be sourcewise different from, that of the background wormhole. The charge ( Q ≠ 0 ) thus provides an extra degree of freedom that introduces a non-trivial redshift function Φ to the background, alters its throat radius to r th , yet keeps the wormhole massless. Second, we focus on this electrically charged case and calculate the light deflection angle up to 4th PPN order, analyze the effect of Q on the lensing observables such as the image positions, magnification, centroid and time delay of images of the source stars. Third, we analyze the probabilistic features such as optical depth and event rate estimated on the basis of the hypothesis that the wormhole lens could be bound or unbound to our Galaxy. Finally, we report an intriguing qualitative prediction that, compared to the Schwarzschild black hole, the Paczyński light curves of the electrically charged wormhole are much dimmer that also show characteristic gutters at the times the source enters and exits the Einstein ring. The gutters gradually come together as Q approaches the extreme limit r th 2 , at which the Einstein radius R E vanishes so that the source crosses it instantly. It is speculated that re-analyzing past data on Galactic microlensing may betray the presence of charged wormholes. [ABSTRACT FROM AUTHOR]

We report the detection of γ-ray emission by the Fermi Large Area Telescope towards the young massive star cluster RCW 38 in the 1–500 GeV photon energy range. We found spatially extended GeV emission towards the direction of RCW 38, which is best modelled by a Gaussian disc of 0.23° radius with a significance of the extension of ∼11.4σ. Furthermore, the spatial correlation with the ionized and molecular gas content favours the hadronic origin of the γ-ray emission. The γ-ray spectrum of RCW 38 has a relatively hard photon index of 2.44 ± 0.03, which is similar to other young massive star clusters. We argue that the diffuse GeV γ-ray emission in this region likely originates from the interaction of accelerated protons in the stellar cluster with the ambient gas. [ABSTRACT FROM AUTHOR]

Ultra-high energy cosmic rays are the most extreme energetic particles detected on Earth, however, their acceleration sites are still mysterious. We explore the contribution of low-luminosity gamma-ray bursts to the ultra-high energy cosmic ray flux, since they form the bulk of the nearby population. We analyse a representative sample of these bursts detected by BeppoSAX, INTEGRAL , and Swift between 1998–2016, and found that in order to reconcile our theoretical flux with the observed flux, these bursts should accelerate at most 10−13 M⊙ of ultra-high energy cosmic rays. [ABSTRACT FROM AUTHOR]

Cosmic rays are highly energetic messengers propagating in magnetized plasma, which are, possibly but not exclusively, accelerated at astrophysical shocks. Amongst the variety of astrophysical objects presenting shocks, the huge circumstellar stellar wind bubbles forming around very massive stars, are potential non-thermal emitters. We present the 1D magneto-hydrodynamical simulation of the evolving magnetized surroundings of a single, OB-type main-sequence |$60\, \rm {\rm M}_{\odot }$| star, which is post-processed to calculate the re-acceleration of pre-existing non-thermal particles of the Galactic cosmic ray background. It is found that the forward shock of such circumstellar bubble can, during the early phase (⁠|$1\, \rm Myr$|⁠) of its expansion, act as a substantial re-accelerator of pre-existing interstellar cosmic rays. This results in an increasing excess emission flux by a factor of 5, the hadronic component producing γ-rays by π0 decay being more important than those by synchrotron and inverse Compton radiation mechanisms. We propose that this effect is at work in the circumstellar environments of massive stars in general and we conjecture that other nebulae such as the stellar wind bow shocks of runaway massive stars also act as Galactic cosmic ray re-accelerators. Particularly, this study supports the interpretation of the enhanced hadronic emission flux measured from the surroundings of κ Ori as originating from the acceleration of pre-existing particles at the forward shock of its wind bubble. [ABSTRACT FROM AUTHOR]

The circumgalactic medium (CGM) contains information on the cumulative effect of galactic outflows over time, generally thought to be caused by feedback from star formation and active galactic nuclei. Observations of such outflows via absorption in CGM gas of quasar sightlines show a significant amount of cold (⁠|${\lesssim}10^4\,{\rm K}$|⁠) gas, which cosmological simulations struggle to reproduce. Here, we use the adaptive mesh refinement hydrodynamical code Ramses to investigate the effect of cosmic rays (CR) on the cold gas content of the CGM using three zoom realizations of a z  = 1 star-forming galaxy with supernova mechanical feedback: one with no CR feedback (referred to as no-CR), one with a medium CR diffusion coefficient |$\kappa = 10^{28} \, \rm {cm^{2}\, s^{-1}}$| (CR−κmed), and one with a high rate of diffusion of |$\kappa = 3\times 10^{29} \, \rm {cm^{2}\,\, s^{-1}}$| (CR−κhigh). We find that, for CR−κmed, the effects of CRs are largely confined to the galaxy itself as CRs do not extend far into the CGM. However, for CR−κhigh, the CGM temperature is lowered and the amount of outflowing gas is boosted. Our CR simulations fall short of the observed Mg  ii covering fraction, a tracer of gas at temperatures |${\lesssim}10^4\,{\rm K}$|⁠ , but the CR−κhigh simulation is more in agreement with covering fractions of C  iv and O  vi , which trace higher temperature gas. [ABSTRACT FROM AUTHOR]

The Snake is a remarkable Galactic Centre radio filament with a morphology characterized by two kinks along its ∼20 arcmin extent. The major and minor kinks are located where the filament is most distorted from a linear magnetized structure running perpendicular to the Galactic plane. We present Chandra , VLA, and MeerKAT data and report the detection of an X-ray and radio source at the location of the major kink. High-resolution radio images of the major kink reveal a compact source with a steep spectrum with spectral index α ∼ −2.7 surrounded by extended emission. The radio luminosity and steep spectrum of the compact source are consistent with a pulsar. We also show flattening of the spectrum and enhanced synchrotron emissivity away from the position of the major kink along the Snake, which suggests injection of relativistic particles along the Snake. We argue that the major kink is created by a fast-moving (∼500–1000 km s−1) object punching into the Snake, distorting its magnetic structure, and producing X-ray emission. X-ray emission pinpoints an active acceleration site where the interaction is taking place. A secondary kink is argued to be induced by the impact of the high-velocity object producing the major kink. [ABSTRACT FROM AUTHOR]

We investigate the nature of a Galactic Centre source, G0.17+0.15, lying along the northern extension of the radio arc near l ∼ 0.2°. G0.17+0.15 is an H  ii region located towards the eastern edge of the radio bubble, embedded within the highly polarized Galactic Centre eastern lobe where a number of radio filaments appear to cross through the H  ii region. We report the detection of hydrogen and helium recombination lines with a radial velocity exceeding 140 km s−1 based on Green Bank Telescope and Very Large Array observations. The morphology of G0.17+0.15, aided by kinematics, and spectral index characteristics, suggests the presence of an external pressure dragging and shredding the ionized gas. We argue that this ionized cloud is interacting with a bundle of radio filaments and is entrained by the ram pressure of the radio bubble, which itself is thought to be produced by cosmic ray driven outflows at the Galactic Centre. In this interpretation, the gas streamers on the western side of G0.17+0.15 are stripped, accelerated from 0 to |$\delta v\sim \, 35$| km s−1 over a time-scale roughly 8 × 104 yr, implying that ablating ram pressure is |$\sim 700\, \mathrm{eV\, cm^{-3}}$|⁠ , comparable to the |$\sim 10^3\, \mathrm{eV \, cm^{-3}}$| cosmic ray driven wind pressure in the Galactic Centre region. [ABSTRACT FROM AUTHOR]

The presence of live 60Fe nuclei (lifetime of 3.8 Myr) in cosmic rays detected by the Advanced Composition Explorer /Cosmic Ray Isotope Spectrometer instrument suggests a nearby nucleosynthesis source. 60Fe is primarily produced in core-collapse supernovae, and we aim to clarify whether the detected 60Fe nuclei can be associated with a particular local supernova. We consider 25 OB associations and subgroups located within 1 kpc of the Solar system based on recent Gaia census. A model is developed that combines stellar population synthesis within these OB associations, cosmic ray acceleration within associated superbubbles, and cosmic ray transport to the Solar system. The most critical model parameter impacting 60Fe cosmic ray production is the explodability criterion, which determines if a massive star ends its life as a supernova. Our study points to the Scorpius–Centaurus (Sco–Cen) OB association as the most probable origin of the observed 60Fe nuclei, particularly suggesting they were accelerated in the Sco–Cen superbubble by a young supernova aged ≤500 kyr with a progenitor mass of approximately 13–20 M⊙. A less likely source is the supernova at the origin of the Geminga pulsar 342 kyr ago, if the progenitor originated in the Orion OB1 association. The contribution of local OB associations to the cosmic ray density of stable 56Fe is estimated to be around 20 per cent, with some sensitivity to cosmic ray acceleration efficiency and diffusion coefficient. These findings shed light on the origins of cosmic ray nuclei, connecting them to nucleosynthesis events within our local cosmic neighbourhood. [ABSTRACT FROM AUTHOR]

In the 1960s, the remnants of supernova explosions (SNRs) were indicated as a possible source of galactic cosmic rays through the Diffusive Shock Acceleration (DSA) mechanism. Since then, the observation of gamma-ray emission from relativistic ions in these objects has been one of the main goals of high-energy astrophysics. A few dozen SNRs have been detected at GeV and TeV photon energies in the last two decades. However, these observations have shown a complex phenomenology that is not easy to reduce to the standard paradigm based on DSA acceleration. Although the understanding of these objects has greatly increased, and their nature as efficient electron and proton accelerators has been observed, it remains to be clarified whether these objects are the main contributors to galactic cosmic rays. Here, we review the observations of γ -ray emission from SNRs and the perspectives for the future. [ABSTRACT FROM AUTHOR]

Cosmic rays (CRs) may drive outflows and alter the phase structure of the circumgalactic medium, with potentially important implications on galaxy formation. However, these effects ultimately depend on the dominant mode of transport of CRs within and around galaxies, which remains highly uncertain. To explore potential observable constraints on CR transport, we investigate a set of cosmological fire -2 CR-magnetohydrodynamic simulations of L* galaxies which evolve CRs with transport models motivated by self-confinement (SC) and extrinsic turbulence (ET) paradigms. To first order, the synchrotron properties diverge between SC and ET models due to a CR physics-driven hysteresis. SC models show a higher tendency to undergo 'ejective' feedback events due to a runaway buildup of CR pressure in dense gas due to the behaviour of SC transport scalings at extremal CR energy densities. The corresponding CR wind-driven hysteresis results in brighter, smoother, and more extended synchrotron emission in SC runs relative to ET and constant diffusion runs. The differences in synchrotron arise from different morphology, interstellar medium gas, and B properties, potentially ruling out SC as the dominant mode of CR transport in typical star-forming L* galaxies, and indicating the prospect for non-thermal radio continuum observations to constrain CR transport physics. [ABSTRACT FROM AUTHOR]

• ORCA is an instrument for the observation of cosmic rays in the Antarctic Peninsula. • Capable of detecting neutrons, muons and muon incidence directions. • Changes observed in the cosmic ray spectrum during a ICME passage in Nov. 2021. • Secondary muons show changes in their arrival directions along the ICME passage. • Different behavior, both spectrum and arrival directions, in shock, sheath and MC. The Antarctic Cosmic Ray Observatory (ORCA) has been in continuous operation since January 2020 and in nominal phase since March 2021 at the Spanish Antarctic Base Juan Carlos I (Livingston Island, Antarctica 62 ° 39 ′ 46 ′ ′ S , 60 ° 23 ′ 20 ′ ′ W at 12 m above sea level) at an effective vertical cutoff rigidity of R = 2.37 GV. ORCA consists of two neutron monitors, ORC-A with three BF 3 counter tubes following the NM64 standard, and ORC-B, with three bare BF 3 counter tubes; and a muon telescope, ORC-M consisting of two 1 m2 scintillators with the neutron monitors between them. ORCA provides neutron count rates at two different energy thresholds (provided by ORC-A and ORC-B) and muon count rates and muon incidence directions throughout the ORCA volume. The neutron-producing cosmic rays observed in ORC-A have a lower rigidity threshold than the muon-producing ones observed in ORC-M. The relationship between them gives an idea about the change in the slope of the cosmic ray spectrum at the corresponding thresholds, 2.37 GV for ORC-A and ≈ 12 GV for ORC-M. This relation allows to investigate the change in the cosmic ray spectrum along the different regions in an Interplanetary Coronal Mass Ejection, i.e. the shock, the sheath and the magnetic cloud. ORCA first observations are presented, focusing on the analysis of the Forbush Decrease (FD) that happened at the beginning of November 2021. We present these observations and show the promising capabilities of the Antarctic Cosmic Ray Observatory, which represents a new concept of cosmic ray and solar energetic particle detector. The cosmic ray spectrum changes during the Forbush decrease detected in early November 2021 as seen in the change in the ratio of ORC-A to ORC-M count rates. This change is clear after the arrival of the shock and during the magnetic cloud being smoother or even flat in the sheath. In addition, the secondary muons show changes in their arrival directions along the ICME, being clear the different behavior in shock, sheath and magnetic cloud. [ABSTRACT FROM AUTHOR]

According to UNSCEAR, cosmic radiation contributes to ~16% (0.39 mSv/y) of the total dose received by the public at sea level. The exposure to cosmic rays at a specific location is therefore a non-negligible parameter that contributes to the assessment of the overall public exposure to radiation. In this study, simulations were conducted with the Particle and Heavy Ion Transport code System, a Monte Carlo code, to determine the fluxes and effective dose due to cosmic rays received by the population of Douala. In minimum solar activity, the total effective dose considering the contribution of neutron, muon+, muon-, electron, positron and photon, was found to be 0.31 ± 0.02 mSv/y at the ground level. For maximum solar activity, it was found to be 0.27 ± 0.02 mSv/y at ground level. During maximum solar activity, galactic cosmic rays are reduced by solar flares and winds, resulting in an increase in the solar cosmic-ray component and a decrease in the galactic cosmic-ray component on Earth. This ultimately leads to a decrease in the total cosmic radiation on Earth. These results were found to be smaller than the UNSCEAR values, thus suggesting a good estimation for the population of Douala city located near the equatorial line. In fact, the cosmic radiation is more deflected at the equator than near the pole. Muons+ were found to be the main contributors to human exposure to cosmic radiation at ground level, with ~38% of the total effective dose due to cosmic exposure. However, electrons and positrons were found to be the less contributors to cosmic radiation exposure. As regards the obtained results, the population of Douala is not significantly exposed to cosmic radiation. [ABSTRACT FROM AUTHOR]

Reanalysis soil moisture products are valuable for diverse applications, but their quality assessment is limited due to scale discrepancies when compared to traditional in situ point-scale measurements. The emergence of cosmic ray neutron sensors (CRNSs) with field-scale soil moisture estimates (∼ 250 m radius, up to 0.7 m deep) is more suitable for the product evaluation owing to their larger footprint. In this study, we perform a comprehensive evaluation of eight widely used reanalysis soil moisture products (ERA5-Land, CFSv2, MERRA2, JRA55, GLDAS-Noah, CRA40, GLEAM and SMAP L4 datasets) against 135 CRNS sites from the COSMOS-UK, COSMOS-Europe, COSMOS USA and CosmOz Australia networks. We evaluate the products using six metrics capturing different aspects of soil moisture dynamics. Results show that all reanalysis products generally exhibit good temporal correlation with the measurements, with the median temporal correlation coefficient (R) values spanning 0.69 to 0.79, though large deviations are found at sites with seasonally varying vegetation cover. Poor performance is observed across products for soil moisture anomalies time series, with R values varying from 0.46 to 0.66. The performance of reanalysis products differs greatly across regions, climate, land covers and topographic conditions. In general, all products tend to overestimate data in arid climates and underestimate data in humid regions as well as grassland. Most reanalysis products perform poorly in steep terrain. Relatively low temporal correlation and high bias are detected in some sites from the west of the UK, which might be associated with relatively low bulk density and high soil organic carbon. Overall, ERA5-Land, CRA40, CFSv2, SMAP L4 and GLEAM exhibit superior performance compared to MERRA2, GLDAS-Noah and JRA55. We recommend that ERA5-Land and CFSv2 could be used in humid climates, whereas SMAP L4 and CRA40 perform better in arid regions. SMAP L4 has good performance for cropland, while GLEAM is more effective in shrubland regions. Our findings also provide insights into directions for improvement of soil moisture products for product developers. [ABSTRACT FROM AUTHOR]

We present six case studies from a comprehensive mass range (1-109 M⊙) of astrophysical objects, each of which possess jets, emit high-energy gamma radiation and in some instances spatially identifiable ultra-high-energy cosmic rays (UHECRs). All sources are strong candidates for UHECR emission, if not already known to emit them. We surmise that wakefield acceleration in conjunction with the magnetorotational instability of the accretion disc explains both structural properties of the jets and details in their emission signals, such as correlations in neutrino and gamma-ray bursts, and in the case of blazars, anticorrelations in flux and spectral index. Furthermore, our model predicts an upper bound for the energy of UHECRs emitted from a source given the mass of its central compact object and total jet luminosity. To provide context for our model predictions, we quantitatively compare them with observational data, however, we have not accounted for the GZK limit and assumed universal values for several model parameters (e.g. jet-spreading index, p) that likely differ between sources. Since the accretion and acceleration mechanisms are independent of mass, aside from determining maximum values, blazars (∼109 M⊙), radio galaxies (∼ 108 M⊙), Seyfert galaxies (∼ 106 M⊙ ), starburst galaxies (∼ 103 M⊙ ), even microquasars (1-10 M⊙) interestingly exhibit the same physics. Other radiation bands, such as X-ray, ultraviolet, or radio, may harbour additional information, but we chose not to focus on them for brevity. However, such an endeavour may open the door to a new multimessenger approach for understanding these objects.

Starting from the perspective of reports published in Brazilian newspapers at the time, as well as letters exchanged between scientists who worked in Brazil and North American colleagues, and documents from the symposium on cosmic rays, a chronological sequence of how the so-called Compton mission in Brazil took place and was perceived by the literate public will be presented. For a better understanding of the article, research on cosmic rays in Brazil and around the world is contextualized and the reasons for the creation of the mission and its coming to Brazil are discussed.

Muography is a rapidly developing and non-destructive tomographic technology that uses cosmic ray muons. Due to the natural presence and deeper penetration of cosmic ray muons, scientists have performed various pioneer studies in fields, such as customs security, the internal imaging of volcanoes, scientific archaeology, and others. With unique advantages, muography has gained increasing attention from archaeologists as a novel and innovative tool to investigate large-scale archaeological sites. This approach may be especially helpful for identifying endangered cultural relics and monuments. In the work, we employ a compact, rugged, and portable muon imaging system, CORMIS (COsmic Ray Muon Imaging System), deployed at up to six measurement locations to perform a case study of three-dimensional muography in Xi'an city, China. Cultural cities, such as Xi'an, have long histories and could benefit from innovative techniques used to investigate, conserve, and protect large historical sites. In this paper, we present in detail a high resolution survey on a rampart of a Xi'an defensive wall in demand of urgent protection. The survey data are carefully processed with advanced statistical methods newly introduced in muography, and the results indicate density anomalies inside the rampart with unprecedented levels of precision. The density anomalies are potential safety hazards and need to be eliminated as soon as possible. The successful implementation of this survey significantly encourages more engagement on the tangible application of high-precision 3D muography in archaeological investigations and protection projects around the world. [ABSTRACT FROM AUTHOR]