Identify signatures of underground muons from atmospheric charm: Simulation study

Muons from the “prompt” decays of charmed mesons in cosmic ray air showers start to show abundance on the atmospheric muon spectrum from few tens of TeV. Study of these prompt muons have broader interest in particle and astroparticle physics. The measurement of prompt muon in air showers is challenging because of their low production rate and the large amount of conventional muons produced in company with them. This paper describes the simulation study of a method that identifies prompt muon signatures based on the pattern of stochastic energy losses by muon bundles in deep under ice. The systematics associated with different hadronic interaction models and cosmic ray primary composition assumptions were estimated. Using IceCube as an example, we briefly discussed the challenge of using this method in experimental data analysis.     

Characterization of the atmospheric muon flux in IceCube

Muons produced in atmospheric cosmic ray showers account for the by far dominant part of the event yield in large-volume underground particle detectors. The IceCube detector, with an instrumented volume of about a cubic kilometer, has the potential to conduct unique investigations on atmospheric muons by exploiting the large collection area and the possibility to track particles over a long distance. Through detailed reconstruction of energy deposition along the tracks, the characteristics of muon bundles can be quantified, and individual particles of exceptionally high energy identified. The data can then be used to constrain the cosmic ray primary flux and the contribution to atmospheric lepton fluxes from prompt decays of short-lived hadrons.In this paper, techniques for the extraction of physical measurements from atmospheric muon events are described and first results are presented. The multiplicity spectrum of TeV muons in cosmic ray air showers for primaries in the energy range from the knee to the ankle is derived and found to be consistent with recent results from surface detectors. The single muon energy spectrum is determined up to PeV energies and shows a clear indication for the emergence of a distinct spectral component from prompt decays of short-lived hadrons. The magnitude of the prompt flux, which should include a substantial contribution from light vector meson di-muon decays, is consistent with current theoretical predictions.The variety of measurements and high event statistics can also be exploited for the evaluation of systematic effects. In the course of this study, internal inconsistencies in the zenith angle distribution of events were found which indicate the presence of an unexplained effect outside the currently applied range of detector systematics. The underlying cause could be related to the hadronic interaction models used to describe muon production in air showers.     

Production of neutrinos and secondary electrons in cosmic sources

We study the individual contribution to secondary lepton production in hadronic interactions of cosmic rays (CRs) including resonances and heavier secondaries. For this purpose we use the same methodology discussed earlier [C.-Y. Huang, S.-E. Park, M. Pohl, C.D. Daniels, Astropart. Phys. 27 (2007) 429], namely the Monte-Carlo particle collision code DPMJET3.04 to determine the multiplicity spectra of various secondary particles with leptons as the final decay states, that result from inelastic collisions of cosmic-ray protons and Helium nuclei with the interstellar medium of standard composition. By combining the simulation results with parametric models for secondary particle (with resonances included) for incident cosmic-ray energies below a few GeV, where DPMJET appears unreliable, we thus derive production matrices for all stable secondary particles in cosmic-ray interactions with energies up to about 10 PeV.

We apply the production matrices to calculate the radio synchrotron radiation of secondary electrons in a young shell-type SNR, RX J1713.7-3946, which is a measure of the age, the spectral index of hadronic cosmic rays, and most importantly the magnetic field strength. We find that the multi-mG fields recently invoked to explain the X-ray flux variations are unlikely to extend over a large fraction of the radio-emitting region, otherwise the spectrum of hadronic cosmic rays in the energy window 0.1–100 GeV must be unusually hard.

We also use the production matrices to calculate the muon event rate in an IceCube-like detector that are induced by muon neutrinos from high-energy γ-ray sources such as RX J1713.7-3946, Vela Jr. and MGRO J2019+37. At muon energies of a few TeV, or in other word, about 10 TeV neutrino energy, an accumulation of data over about 5–10 years would allow testing the hadronic origin of TeV γ-rays.     

Measurement of the cosmic ray composition at the knee with the SPASE-2/AMANDA-B10 detectors

The mass composition of high-energy cosmic rays at energies above 10¹⁵ eV can provide crucial information for the understanding of their origin. Air showers were measured simultaneously with the SPASE-2 air shower array and the AMANDA-B10 Cherenkov telescope at the South Pole. This combination has the advantage to sample almost all high-energy shower muons and is thus a new approach to the determination of the cosmic ray composition. The change in the cosmic ray mass composition was measured versus existing data from direct measurements at low energies. Our data show an increase of the mean log atomic mass lnA by about 0.8 between 500 TeV and 5 PeV. This trend of an increasing mass through the “knee” region is robust against a variety of systematic effects.     

Monte-Carlo simulations of atmospheric muon production: Implication of the past martian environment

In order to evaluate the obliquity-driven atmospheric-density path length effect on nuclide production rate on Mars, we performed a Monte-Carlo simulation to produce the number of secondary particles such as muons, neutrons and protons in the martian atmosphere and to simulate that production of ¹⁰Be and ³⁶Cl in the martian regolith by muons and neutrons depends on how much atmosphere had been present for the past 10 million years. The vertical profile of the present martian atmosphere to generate secondary particles has been determined based on the data provided by the Viking missions. For other thickness profiles, we scaled Linsley's atmospheric model. Atmospheric shower has been generated with the SIBYLL 2.1 for high-energy hadronic interactions and EHSA for low energy photonuclear interactions. With increasing atmospheric thickness, more primary interactions occur in the atmosphere. Consequently the proton flux is reduced and the secondary cosmic ray flux, such as muons or energetic neutrons increases at surface. The result indicates that the muon production is more sensitive to obliquity-driven atmospheric variations than proton reduction. A thicker atmosphere would result in enhanced nuclide production at a place deeper than 5 m below the surface and the nuclides present in detectable concentrations. Application to the polar deposit is described.     

Overview of Direct Measurements of Cosmic Rays

The history of cosmic ray studies can be traced back to the 1910s when Hess and other scientists first discovered them. Cosmic rays are very important laboratories of particle physics, and have led to many important discoveries of fundamental particles, such as the positrons, muons, pions, and a series of strange particles. Cosmic rays are nowadays the key probes of the extremely high-energy physics and dark matter particles. A brief review about the history and recent progresses of direct observations of cosmic rays is presented. In recent years, the new space-borne experiments such as PAMELA and AMS-02, as well as a few of balloon-borne experiments, have measured the energy spectra of cosmic rays very precisely, and revealed several new features/anomalies. Remarkable excesses of positron fraction in the total electron plus positron fluxes have been observed, which may be caused by the annihilation/decay of dark matter particles or by astrophysical pulsars. The cosmic ray antiprotons, which are expected to have the same secondary origin as that of positrons, do not show significant excesses compared with the background prediction. This result also constrains the modeling of the positron excesses. In addition, the spectral hardening above several hundred GeV of cosmic ray nuclei has been revealed. These results have important and interesting implications on our understandings of the origin, acceleration, and propagation of cosmic rays. In particular, China has launched the Dark Matter Particle Explorer (DAMPE) to indirectly search for the dark matter and explore the high-energy universe in the TeV window. Most recently, the DAMPE collaborators reported the new measurements of the cosmic ray electron plus positron fluxes up to about 5 TeV with a very high precision. The DAMPE data revealed clearly a deflection around 0.9 TeV in the electron energy spectrum. Possible fine structures of the electron plus positron spectra can be critically addressed with the accumulation of data in the coming years.     

Comparison of measured and simulated lateral distributions for electrons and muons with KASCADE

Lateral distributions for electrons and muons in extensive air showers measured with the array of the KASCADE experiment are compared to results of simulations based on the high-energy hadronic interaction models QGSJet and SIBYLL. It is shown, that the muon distributions are well described by both models. Deviations are found for the electromagnetic component, where both models predict a steeper lateral shape than observed in the data. For both models the observed lateral shapes of the electron component indicate a transition from a light to a more heavy composition of the cosmic ray spectrum above the knee.     

宇宙线直接探测进展概述

宇宙线从发现起至今已超过百年。在20世纪上半叶,大型粒子加速器技术成熟以前,对宇宙线的研究引领着基本粒子物理的发展,从宇宙线研究中取得的多项成果斩获诺贝尔奖。21世纪,宇宙线因其与极端高能的物理规律和暗物质等新物理现象联系密切而绽放出新的活力,宇宙线起源、加速、传播等相关的天文学及物理学问题也备受关注。简述了近年来在空间直接观测宇宙线实验方面取得的进展,以及其对理解宇宙线物理问题的推动。最后概述了中国在相关领域的研究历程和现状。     

Distributions of secondary muons at sea level from cosmic gamma rays below 10 TeV

The Monte Carlo program is used to predict the distributions of the muons which originate from primary cosmic gamma rays and reach sea level. The main result is the angular distribution of muons produced by vertical gamma rays which is necessary to predict the inherent angular resolution of any instrument utilizing muons to infer properties of gamma ray primaries. Furthermore, various physical effects are discussed which affect these distributions in differing proportions.     

High energy neutrinos from the TeV Blazar 1ES 1959 + 650

The AMANDA neutrino telescope has recently reported the detection of high-energy neutrinos spatially and temporally coincident with the flaring of the TeV blazar 1ES 1959 + 650. If high-energy neutrinos are in fact generated by this blazar, it would be the first strong evidence for the hadronic acceleration of cosmic rays. At present, the statistical significance of this observation cannot be reliably assessed, however. In this letter, we investigate whether circumstances exist where the source can produce the flux implied by the coincident events. We show that if the TeV gamma-ray emission observed from 1ES 1959 + 650 or other nearby TeV blazars is the result of accelerated protons interacting with nucleons, it is reasonable that AMANDA could detect several events during a flaring period. Such rates require that the spectral index of the source be rather high (for instance 2.8 for 1ES 1959 + 650) and that the Lorentz factor of the jet be fairly small (Γ  1).     

Performance of the Tibet hybrid experiment (YAC-II + Tibet-III + MD) to measure the energy spectra of the light primary cosmic rays at energies 50–10,000 TeV

A new hybrid detector system has been constructed by the Tibet ASγ collaboration at Tibet, China, since 2014 to measure the chemical composition of cosmic rays around the knee in the wide energy range. They consist of an air-shower-core detector-grid (YAC-II) to detect high energy electromagnetic component, the Tibet air-shower array (Tibet-III) and a large underground water-Cherenkov muon-detector array (MD). We have carried out a detailed air-shower Monte Carlo (MC) simulation to study the performance of the hybrid detectors by using CORSIKA (version 6.204), which includes QGSJET01c and SIBYLL2.1 hadronic interaction models. Assumed primary cosmic ray models are based on helium poor, helium rich and Gaisser’s fit compositions around the knee. All detector responses are calculated using Geant4 (version 9.5) according to the real detector configurations and the MC events are reconstructed by the same procedure as the experimental data analysis. The energy determination is made by lateral density fitting (LDF) method using modified NKG function and the separation of the light components (proton, helium) is made by means of the artificial neural network (ANN) method and the random forest (RF) method. The systematic errors of the spectra of proton and helium caused by each steps of the analysis procedure are investigated including the dependence of the MC data on the hadronic interaction models and the primary composition models, and the algorithms for the primary mass identification. The systematic errors of the flux to be obtained by the new experiment are summarized as less than 30% in total. Our results show that the new hybrid experiment is powerful enough to study the chemical composition of the cosmic rays, in particular, to obtain the light-component spectra of the primary cosmic rays in 50–10,000 TeV energy range overlapping to the direct observation data at low energy side and ground-based indirect observations at high energy side. It is possible in this energy range to find the break points of the power indices of proton and helium (the knee of individual component spectrum) which are basically important parameter for the study of the cosmic-ray origin.     

The cosmic ray primary composition between 10 and 10 eV from Extensive Air Showers electromagnetic and TeV muon data

The cosmic ray primary composition in the energy range between 10¹⁵ and 10¹⁶ eV, i.e., around the “knee” of the primary spectrum, has been studied through the combined measurements of the EAS-TOP air shower array (2005 m a.s.l., 10⁵ m² collecting area) and the MACRO underground detector (963 m a.s.l., 3100 m w.e. of minimum rock overburden, 920 m² effective area) at the National Gran Sasso Laboratories. The used observables are the air shower size (N_e) measured by EAS-TOP and the muon number (N_μ) recorded by MACRO. The two detectors are separated on average by 1200 m of rock, and located at a respective zenith angle of about 30°. The energy threshold at the surface for muons reaching the MACRO depth is approximately 1.3 TeV. Such muons are produced in the early stages of the shower development and in a kinematic region quite different from the one relevant for the usual N_μ−N_e studies. The measurement leads to a primary composition becoming heavier at the knee of the primary spectrum, the knee itself resulting from the steepening of the spectrum of a primary light component (p, He) of Δγ=0.7±0.4 at E₀4×10¹⁵ eV. The result confirms the ones reported from the observation of the low energy muons at the surface (typically in the GeV energy range), showing that the conclusions do not depend on the production region kinematics. Thus, the hadronic interaction model used (CORSIKA/QGSJET) provides consistent composition results from data related to secondaries produced in a rapidity region exceeding the central one. Such an evolution of the composition in the knee region supports the “standard” galactic acceleration/propagation models that imply rigidity dependent breaks of the different components, and therefore breaks occurring at lower energies in the spectra of the light nuclei.     

Calculation of the TeV prompt muon component in very high energy cosmic ray showers

HEMAS-DPM is a Monte Carlo for the simulation of very high energy cosmic ray showers, which includes the DPMJET-II code based on the two component Dual Parton Model. DPMJET-II provides also charm production in agreement with data and, for p exceeding 5 GeV/c, with perturbative QCD results in proton-proton interactions. In this respect, a new scheme has been considered for the inclusive production of D mesons at large p in hadronic collisions in the framework of perturbative fragmentation functions, allowing an analysis at the Next to Leading Order (NLO) level which goes beyond the fixed I(_s³) perturbative theory of open charm production. We have applied HEMAS-DPM to the calculation of the prompt muon component for E_μ ≥ 1 TeV in air showers considering the two extreme cases of primary protons and Fe nuclei.     

Diffuse TeV emission at the Galactic Centre

The High-Energy Stereoscopic System (HESS) has detected intense diffuse TeV emission correlated with the distribution of molecular gas along the Galactic ridge at the centre of our Galaxy. Earlier HESS observations of this region had already revealed the presence of several point sources at these energies, one of them (HESS J1745−290) coincident with the supermassive black hole Sagittarius A*. It is still not entirely clear what the origin of the TeV emission is, nor even whether it is due to hadronic or leptonic interactions. It is reasonable to suppose, however, that at least for the diffuse emission, the tight correlation of the intensity distribution with the molecular gas indicates a pionic-decay process involving relativistic protons. In this paper, we explore the possible source(s) of energetic hadrons at the Galactic Centre, and their propagation through a turbulent medium. We conclude that though Sagittarius A* itself may be the source of cosmic rays producing the emission in HESS J1745−290, it cannot be responsible for the diffuse emission farther out. A distribution of point sources, such as pulsar wind nebulae dispersed along the Galactic plane, similarly do not produce a TeV emission profile consistent with the HESS map. We conclude that only a relativistic proton distribution accelerated throughout the intercloud medium can account for the TeV emission profile measured with HESS.     

Comparison of AGASA data with CORSIKA simulation

An interpretation of Akeno giant air shower array (AGASA) data by comparing the experimental results with the simulated ones by cosmic ray simulation for KASCADE (CORSIKA) has been made. General features of the electromagnetic component and low energy muons observed by AGASA can be well reproduced by CORSIKA. The form of the lateral distribution of charged particles agrees well with the experimental one between a few hundred metres and 2000 m from the core, irrespective of the hadronic interaction model studied and the primary composition (proton or iron). It does not depend on the primary energy between 10^17.5 and 10²⁰ eV as the experiment shows. If we evaluate the particle density measured by scintillators of 5 cm thickness at 600 m from the core S₀(600), suffix 0 denotes the vertically incident shower) by taking into account the similar conditions as in the experiment, the conversion relation from S₀(600) to the primary energy is expressed as E (eV)=2.15×10¹⁷S₀(600)^1.015 within 10% uncertainty among the models and composition used, which suggests the present AGASA conversion factor is the lower limit. Although the form of the muon lateral distribution fits well to the experiment within 1000 m from the core, the absolute values change with hadronic interaction model and primary composition. The slope of the ρ_μ(600) (muon density above 1 GeV at 600 m from the core) vs. S₀(600) relation in experiment is flatter than that in simulation of any hadronic model and primary composition. As the experimental slope is constant from 10¹⁵ to 10¹⁹ eV, we need to study this relation in a wide primary energy range to infer the rate of change of chemical composition with energy.     

Simulation of extensive air showers at ultra-high energy using the CORSIKA Monte Carlo code

Extensive simulations have been carried out with CORSIKA version 5.62 to investigate the general properties of giant cosmic air shower in the energy range 1–100 EeV. The comparison between protons, heavy nuclei and γ initiated showers exhibits unexpected and interesting features. The apparent muon electron ratio at great distances (1.5 km from axis) tends to be comparable at ultra-high energy in both photon-induced cascades and hadronic cascades (compensation between the enhancement versus energy of photo-production cross-section and of the decrease of both pair production cross-section and the bremsstrahlung cross-section, with Landau–Pomeranchuk–Migdal effect); for proton and nuclei primaries, a correlation with lateral electron profile suggests a new energy estimator, in complement to electrons size or density at 600 m, suitable for the determination of the total primary energy spectrum. Another tendency is the local contrast in the abundance of positive and negative muons (with a possible ellipticity in the lateral muon distribution) induced by the geomagnetic field, especially visible for some azimuthal and zenith angles. These distortions are more intense for heavy primaries; they can be exploited on the most favorable horizontal axis or areas, for the discrimination between nuclei and protons.     

Particle acceleration test with Cas A multiwavelength emission

The investigation of supernova remnants (SNRs) across the electromagnetic spectrum from radio up to very high energy gamma-rays can serve as a test of the particle acceleration and touches on one of the unresolved problems of modern astrophysics, namely the origin of cosmic rays and the Galaxy's contribution to the overall cosmic ray spectrum. The multiwavelength observations of Cas A SNR demonstrated that structure and spectral features have clear signs of young SNRs and its overall properties make this object the best target to test a hypothesis of cosmic ray origin in SNRs. Studies of Cas A at very high energies by SHALON telescope showed the location of TeV gamma-ray emission region relative to the position of reveres shock. Also, the spectral energy distribution was obtained at high and very high energies. To describe the spectral and structural features of this SNR viewed in non-thermal emission, two approaches involving reverse and also both reverse and forward shocks to the mechanism of diffusive shock acceleration of cosmic rays in Cas A were applied. It is demonstrated that the observational properties of Cas A are well reproduced by the hadronic model with significant contribution of both the forward and reverse shocks in the generation of broadband emission. Calculation results suggest that the very high efficiency of particle acceleration in Cas A, which value is up to 25% of the supernova explosion energy with energy of accelerated particles not exceeding of $$ eV. Whereas, the forward shock model predicts the spectral characteristics of the TeV-gamma-emission corresponding to ones detected at 800 GeV–40 TeV that are the evidence of acceleration of the hadronic cosmic rays in shells of SNRs up to $$ eV     

High energy neutrinos from the TeV Blazar 1ES 1959 + 650

The AMANDA neutrino telescope has recently reported the detection of high-energy neutrinos spatially and temporally coincident with the flaring of the TeV blazar 1ES 1959 + 650. If high-energy neutrinos are in fact generated by this blazar, it would be the first strong evidence for the hadronic acceleration of cosmic rays. At present, the statistical significance of this observation cannot be reliably assessed, however. In this letter, we investigate whether circumstances exist where the source can produce the flux implied by the coincident events. We show that if the TeV gamma-ray emission observed from 1ES 1959 + 650 or other nearby TeV blazars is the result of accelerated protons interacting with nucleons, it is reasonable that AMANDA could detect several events during a flaring period. Such rates require that the spectral index of the source be rather high (for instance ∼2.8 for 1ES 1959 + 650) and that the Lorentz factor of the jet be fairly small (Γ ∼ 1).     

Modeling the emission processes in blazars

Blazars are the most violent steady/recurrent sources of high-energy gamma-ray emission in the known Universe. They are prominent emitters of electromagnetic radiation throughout the entire electromagnetic spectrum. The observable radiation most likely originates in a relativistic jet oriented at a small angle with respect to the line of sight. This review starts out with a general overview of the phenomenology of blazars, including results from a recent multiwavelength observing campaign on 3C279. Subsequently, issues of modeling broadband spectra will be discussed. Spectral information alone is not sufficient to distinguish between competing models and to constrain essential parameters, in particular related to the primary particle acceleration and radiation mechanisms in the jet. Short-term spectral variability information may help to break such model degeneracies, which will require snap-shot spectral information on intraday time scales, which may soon be achievable for many blazars even in the gamma-ray regime with the upcoming GLAST mission and current advances in Atmospheric Cherenkov Telescope technology. In addition to pure leptonic and hadronic models of gamma-ray emission from blazars, leptonic/hadronic hybrid models are reviewed, and the recently developed hadronic synchrotron mirror model for TeV γ-ray flares which are not accompanied by simultaneous X-ray flares (“orphan TeV flares”) is revisited. The U.S. Government’s right to retain a non-exclusive, royalty-free license in and to any copyright is acknowledged.