Search

Since the discovery of cosmic rays (CRs) over a century ago, their origin has remained a mystery and a key research question. Recently, the LHAASO experiment identified the first CR super-acceleration source, the Cygnus bubble, which can accelerate CRs to energies exceeding 10 PeV. A pertinent question is: how much does the Cygnus bubble contribute to the CR spectrum observed on Earth? With the aim of answering that question, a 3D propagation analysis was conducted on CRs in this study. The Cygnus bubble was incorporated into our propagation model in order to determine its contributions to the observed spectra. First, we calculated the spectrum and spatial morphology of the Cygnus bubble to reproduce the observed LHAASO data. Subsequently, we calculated the diffuse $\gamma$-ray emissions produced by the CRs from the Cygnus bubble and the energy spectrum of the cosmic ray particles near Earth after propagation. Finally, we utilized a CR spatial-dependent propagation model to calculate the large-scale CR energy spectrum and the resulting diffuse $\gamma$-ray emissions. Our results indicate that: (1) the Cygnus bubble contributes minimally to the CR spectrum observed on Earth, (2) the emissions produced by the CR particles from the Cygnus bubble dominates the diffuse $\gamma$-ray emissions in that region, (3) the structural fluctuations of the diffuse $\gamma$-ray emissions observed by LHAASO are likely due to the local CR halo. We anticipate that LHAASO will identify more CR halo sources to validate our model.

The Sun stands out as the closest and clearest astrophysical accelerator of cosmic rays, while other objects within and beyond the galaxy remain enigmatic. It is probable that the cosmic ray spectrum and mass components from these celestial sources share similarities, offering a novel approach to study their origin. In this study, we analyze of spectra and mass in the energy range from MeV to 10~EeV. We find: (1) the mean-logarithmic mass $\rm\left\langle lnA \right\rangle$ distribution with energy exhibits much clearer feature structures than the spectra; (2) a 100~TeV bump is presented in the $\rm\left\langle lnA \right\rangle$ distribution; (3) for protons, the knee is located at $\sim2$ PeV, the boundary between the galaxy and extra-galaxy occurs at $\sim30$ PeV, marked by a sharp dip; (4) the all-particle spectrum exhibits hardening at $\sim30$~PeV due to the contribution of nearby galaxies, and the extra-galactic dominate $\sim0.7$~EeV. We hope the LHAASO experiment can perform spectral measurements of individual species to validate our results., Comment: 6 pages, 3 figures

Recently LHAASO Collaboration gives precise measurements of cosmic rays (CR) all particle energy spectrum and mean logarithmic mass $\left \langle \ln A \right \rangle$ from 0.3 PeV to 30 PeV. Combining the CR measurements by AMS-02 and DAMPE in space and that by LHAASO and Auger on the ground we construct a model to recover all these measurements from tens of GeV to tens of EeV. We find the LHAASO measurement is crucial in the model construction by connecting the Galactic component to the extragalactic component. The precise measurements of CR spectra for individual species by AMS-02 and DAMPE together with the newest LHAASO results clearly indicates three Galactic CR components, that is, a soft low energy background, a hard high energy component, and a local source contribution. However, the LHAASO data show that above $\sim 10^{16}$ eV a nonnegligible extragalactic component must be included. Combining the Auger results and the LHAASO results we figure out the extragalactic CRs which need at least two components at lower and higher energies. Thanks to the precise measurements by LHAASO the constraints on the model parameters are quite stringent. The spectra features and mass measurements in all energy range are all well reproduced in the model., Comment: 11 pages, 2 figures, 4 tables

Cosmic rays travel throughout the Galaxy, leaving traces from radio to ultra-high-energy gamma rays due to interactions with the interstellar gas, radiation field and magnetic field. Therefore, it is necessary to utilize multi-wavelength investigations on the Galactic diffuse emission to shed light on the physics of CR production and propagation. In this work, we present a spatially dependent propagation scenario, taking account of a local source contribution, while making allowances for an additional CR component freshly accelerated near their sources. In this picture, after reproducing the particle measurements at the Solar system, we calculated the intensity and compared the spectral energy distribution to observations from Fermi-LAT and LHAASO-KM2A in the gamma-ray band, and from WMAP and Planck among other radio surveys at lower energies. Multi-band data considered in conjunction, the former comparison exhibits sufficiently good consistency in favor of our model, while the latter calls for improvement in data subtraction and processing. From this standpoint, there remains potential for advanced observations at energies from milli-eVs to MeVs towards the Galactic plane, in order to evaluate our model further and more comprehensively in the future., Comment: Accepted by ApJ

The detection of GW170817/GRB170817A implied the strong association between short gamma-ray bursts (SGRBs) and binary neutron star (BNS) mergers which produce gravitational waves (GWs). More evidence is needed to confirm the association and reveal the physical processes of BNS mergers. The upcoming High Altitude Detection of Astronomical Radiation (HADAR) experiment, excelling in a wide field of view (FOV) and a large effective area above tens of GeV, is a hope for the prompt detection of very-high-energy (VHE; > 10 GeV) SGRBs. The aim of this paper is to simulate and analyse GW/SGRB joint detections by future GW detector networks in synergy with HADAR, including the second generation LIGO, Virgo and KAGRA and the third generation ET and CE. We provide a brief introduction of the HADAR experiment for SGRB simulations and its expected SGRB detections. For GW simulations, we adopt a phenomenological model to describe GWs produced by BNS mergers and introduce the signal-noise ratios (SNRs) as detector responses. Following a theoretical analysis we compute the redshift-dependent efficiency functions of GW detector networks. We then construct the simulation of GW detection by Monte Carlo sampling. We compare the simulated results of LIGO-Virgo O2 and O3 runs with their actual detections as a check. The combination of GW and SGRB models is then discussed for joint detection, including parameter correlations, triggered SNRs and efficiency skymaps. The estimated joint detection rates are 0.09-2.52 per year for LHVK network with HADAR under different possible configurations, and approximately 0.27-7.89 per year for ET+CE network with HADAR.

The AMS-02 experiment has observed new properties of primary cosmic rays (CRs) categorized into two groups: He-C-O-Fe and Ne-Mg-Si-S, which are independent of CR propagation. In this study, we investigate the unexpected properties of these nuclei using a spatial propagation model. All nuclei spectra are accurately reproduced and separated into primary and secondary contributions. Our findings include: 1. Primary CR spectra are identical. 2. Our calculations align with AMS-02 results for primary-dominated nuclei within a 10\% difference, but show significant discrepancies for the secondary-dominated nuclei. 3. The primary element abundance is presented for the first time. We anticipate that the DAMPE and future HERD experiments will provide observations of nuclei spectra above TeV energy., Comment: 8 pages, 8 figures. accepted by PRD

The HADAR experiment, which will be constructed in Tibet, China, combines the wide-angle advantages of traditional EAS array detectors with the high sensitivity advantages of focused Cherenkov detectors. Its physics objective is to observe transient sources such as gamma-ray bursts and counterparts of gravitational waves. The aim of this study is to utilize the latest AI technology to enhance the sensitivity of the HADAR experiment. We have built training datasets and models with distinctive creativity by incorporating relevant physical theories for various applications. They are able to determine the kind, energy, and direction of incident particles after careful design. We have obtained a background identification accuracy of 98.6%, a relative energy reconstruction error of 10.0%, and an angular resolution of 0.22-degrees in a test dataset at 10 TeV. These findings demonstrate the enormous potential for enhancing the precision and dependability of detector data analysis in astrophysical research. Thanks to deep learning techniques, the HADAR experiment's observational sensitivity to the Crab Nebula has surpassed that of MAGIC and H.E.S.S. at energies below 0.5 TeV and remains competitive with conventional narrow-field Cherenkov telescopes at higher energies. Additionally, our experiment offers a fresh approach to dealing with strongly connected scattered data.

Recent observations of cosmic rays (CRs) have revealed a two-component anomaly in the spectra of primary and secondary particles, as well as their ratios, prompting investigation into their common origin. In this study, we incorporate the identification of slow diffusion zones around sources as a common phenomenon into our calculations, which successfully reproduces all previously described anomalies except for the positron spectrum. Crucially, our research offers a clear physical picture of the origin of CR: while high-energy ($\textrm{>200~GV}$, including the knee) particles are primarily produced by fresh accelerators and are confined to local regions, low energy ($\textrm{<200~GV}$) components come from distant sources and travel through the outer diffusive zone outside of the galactic disk. This scenario can be universally applied in the galactic disk, as evidenced by ultra-high energy diffuse $\rm\gamma$-ray emissions detected by the AS$\rm\gamma$ experiment. Furthermore, our results predict that the spectrum of diffuse $\rm\gamma$-ray is spatial-dependent, resting with local sources, which can be tested by LHAASO experiment., Comment: 9 pages, 7 figures, accepted by PhysRevD

The observations of the energy spectra of cosmic-ray have revealed complicated structures. Especially, spectral hardenings in the boron-to-carbon and boron-to-oxygen ratios above $\sim 200$ GV has been revealed by AMS-02 and DAMPE experiments. One scenario to account for the hardenings of secondary-to-primary ratios is the nuclear fragmentation of freshly accelerated particles around sources. In this work, we further study this scenario based on new observations of Galactic diffuse gamma rays by LHAASO and neutrinos by IceCube. We find that the spectra of cosmic ray nuclei, the diffuse ultra-high-energy gamma rays, and the Galactic component of neutrinos can be simultaneously explained, given an average confinement and interaction time of $\sim 0.25$ Myr around sources. These multi-messenger data thus provide evidence of non-negligible grammage of Galactic cosmic rays surrounding sources besides the traditional one during the propagation., Comment: 10 papes, 5 pages

The observation of short gamma ray bursts (SGRBs) in the TeV energy range plays an important role in understanding the radiation mechanism and probing new areas of physics such as Lorentz invariance violation. However, no SGRB has been observed in this energy range due to the short duration of SGRBs and the weakness of current experiments. New experiments with new technology are required to detect sub-TeV SGRBs. In this work, we observe the very high energy (VHE) $\gamma$-ray emissions from SGRBs and calculate the annual detection rate with the High Altitude Detection of Astronomical Radiation HADAR (HADAR) experiment. First, a set of pseudo-SGRB samples is generated and checked using the observations of Fermi-GBM, Fermi-LAT, and SWIFT measurements. The annual detection rate is calculated from these SGRB samples based on the performance of the HADAR instrument. As a result, the HADAR experiment can detect 0.5 SGRB per year if the spectral break-off of $\gamma$-rays caused by the internal absorption is larger than 100 GeV. For a GRB09010-like GRB in HADAR's view, it should be possible to detect approximately 2000 photons considering the internal absorption. With a time delay assumption due to the Lorentz invariance violation effects, a simulated light curve of GRB090510 has evident energy dependence. We hope that the HADAR experiment can perform the SGRB observations and test our calculations in the future.

The multi-messenger anomalies, including spectral hardening or excess for nuclei, leptons, ratios of $\bar p/p$ and B/C, and anisotropic reversal, were observed in past years. AMS-02 experiment also revealed different spectral break for positron and electron at 284 GeV and beyond TeV respectively. It is natural to ask whether all those anomalies originate from one unified physical scenario. In this work, the spatially-dependent propagation (SDP) with a nearby SNR source is adopted to reproduce above mentioned anomalies. There possibly exists dense molecular cloud(DMC) around SNRs and the secondary particles can be produced by pp-collision or fragmentation between the accelerated primary cosmic rays and DMC. As a result, the spectral hardening for primary, secondary particles and ratios of $B/C$ and $\bar p/p$ can be well reproduced. Due to the energy loss at source age of 330 kyrs, the characteristic spectral break-off for primary electron is at about 1 TeV hinted from the measurements. The secondary positron and electron from charged pion take up $5\%$ energy from their mother particles, so the positron spectrum has a cut-off at $\sim$250 GeV. Therefore, the different spectral break for positron and electron together with other anomalies can be fulfilled in this unified physical scenario. More interesting is that we also obtain the featured structures as spectral break-off at 5 TV for secondary particles of Li, Be, B, which can be served to verify our model. We hope that those tagged structures can be observed by the new generation of space-borne experiment HERD in future., Comment: 16 pages, 12 figures

More and more experiments have identified that the energy spectra of both primary and secondary cosmic-rays exhibit a hardening above $\sim 200$ GV. Most recently, the DAMPE experiment has reported a hardening of boron-to-carbon ratio at $200$ GV. These signs call for modifications of the conventional cosmic-ray (CR) picture. In this work, we propose that the plethoric secondary cosmic rays, for example, boron, antiprotons, originate from the hadronic interactions of freshly accelerated cosmic rays with the interstellar gas near the sources. We find that secondary-to-primary ratios, for example, boron-to-carbon, boron-to-oxygen and antiproton-to-proton ratios, could be well described. The measurements of electrons and positrons could also be accounted for., Comment: 13 pages, 4 figures, published in JCAP

The energy spectra and anisotropies are very important probes of the origin of cosmic rays. Recent measurements show that complicated but very interesting structures exist, at similar energies, in both the spectra and energy-dependent anisotropies, indicating a common origin of these structures. Particularly interesting phenomenon is that there is a reversal of the phase of the dipole anisotropies, which challenges a theoretical modeling. In this work, for the first time, we identify that there might be an additional phase reversal at $\sim 100$ GeV energies of the dipole anisotropies as indicated by a few underground muon detectors and the first direct measurement by the Fermi satellite, coincident with the hundreds of GV hardenings of the spectra. We propose that these two phase reversals, together with the energy-evolution of the amplitudes and spectra, can be naturally explained with a nearby source overlapping onto the diffuse background. As a consequence, the spectra and anisotropies can be understood as the scalar and vector components of this model, and the two reversals of the phases characterize just the competition of the cosmic ray streamings between the nearby source and the background. The alignment of the cosmic ray streamings along the local large-scale magnetic field may play an important but sub-dominant role in regulating the cosmic ray propagation. More precise measurements of the anisotropy evolution at both low energies by space detectors and high energies by air shower experiments for individual species will be essential to further test this scenario., Comment: 9 pages, 2 figures, accepted by APJ

The diffuse $\gamma$-ray emission was measured up to $957$ TeV by the Tibet-AS$\gamma$ experiment recently. Assuming that it is produced by the hadronic interaction between cosmic ray nuclei and the interstellar medium, it requires that the cosmic ray nuclei should be accelerated well beyond PeV energies. Measurements of the cosmic ray spectra for different species show diverse results at present. The Tibet experiments showed that the spectrum of proton plus helium has an early knee below PeV. If this is correct, the diffuse $\gamma$-ray emission would suggest an additional component of Galactic cosmic rays above PeV energies. This second component may originate from a source population of so-called PeVatrons revealed by recent ultra-high energy $\gamma$-ray observations, and could contribute to the cosmic ray fluxes up to the energy of the second knee. On the other hand, the KASCADE measurement showed that the knee of protons is higher than PeV. In this case, the diffuse $\gamma$-rays observed by Tibet-AS$\gamma$ can be well accounted for by only one cosmic ray component. These two scenarious (ie. the Tibet and KASCADE knees) could be distinguished by the spectral structures of diffuse $\gamma$-rays and cosmic ray nuclei. Future measurements of spectra of individual nuclei by HERD and LHAASO experiments and diffuse $\gamma$-rays by LHAASO can jointly constrain these two scenarios., Comment: 9 pages,4 figures. accepted by ApJ

The local sources, such as Geminga SNR, may play important role for the anomaly of proton, electron and anisotropy in the past works. In fact, there exists twelve SNRs around solar system within $1$ kpc. One question is that can other SNRs also possibly contribute the spectra of nuclei and electron and explain the special structure of anisotropy? In this work, under the spatial-dependent propagation, we systematically study the contribution of all local SNRs within 1 kpc around solar to the spectra of nuclei and electron, as well as the energy dependence of anisotropy. As a result, only Geminga, Monogem, and Vela SNRs have quantitive contribution to the nuclei and electron spectra and anisotropy. Here, Geminga SNR is the sole optimal candidate and Monogem SNR is controversial due to the tension of anisotropy between model calculation and observations. The Vela SNR contributes a new spectral structure beyond TeV energy, hinted by HESS, VERITAS, DAMPE and CALET measurements. More interesting is that the electron anisotropy satisfies the Fermi-LAT limit below TeV energy, but rises greatly and reaches $10\%$ at several TeV. This new structure will shed new light to check our model. We hope that the new structure of electron spectrum and anisotropy can be observed by space-borne DAMPE and HERD and ground-based HAWC and LHAASO experiments in the near future., Comment: 11 pages, 8 figures, accepted by APJ

The diffuse $\gamma$-ray spectrum at sub-PeV energy region has been measured for the first time by the Tibet-AS$\gamma$ experiment. It will shed new light on the understanding of origin and propagation of Galactic cosmic rays at very high energies. It has been pointed out that the traditional cosmic ray propagation model based on low energy measurements undershoot the new data, and modifications of the model with new ingredients or alternative propagation framework is required. In this work, we propose that the hadronic interactions between freshly accelerated cosmic rays and the medium surrounding the sources, which was neglected in the traditional model, can naturally account for the Tibet-AS$\gamma$ diffuse emission. We show that this scenario gives a consistent description of other secondary species such as the positron spectrum, the boron-to-carbon ratio, and the antiproton-to-proton ratio. As a result, the electron spectrum above 10 TeV will have a hardening due to this secondary component, which may be tested by future measurements., Comment: 9 pages, 4 figures

Recent studies have shown that the anisotropy is of great value to decipher cosmic rays' origin and propagation. We have built an unified scenario to describe the observations of the energy spectra and the large-scale anisotropy and called attention to their synchronously evolution with energy. In this work, the impact of of the local regular magnetic field (LRMF) and corresponding anisotropic diffusion on large-scale anisotropy have been investigated. When the perpendicular diffusion coefficient is much smaller than the parallel one, the dipole anisotropy points to the LRMF and the observational phase below $100$ TeV could be reproduced. Moreover we find that the dipole phase above $100$ TeV strongly depends on the evolution of local diffusion. But the current measurements at that energy are still scarce. We suggest that more precise measurements at that energy could be carried out to unveil the local diffusion and further the local turbulence., Comment: 10 pages, 5 figures

The precise measurements of energy spectra and anisotropy could help us uncover the local cosmic-ray accelerators. Our recent works have shown that spectral hardening above $200$ GeV in the energy spectra and transition of large-scale anisotropy at $\sim 100$ TeV are of local source origin. Less than $100$ TeV, both spectral hardening and anisotropy explicitly indicate the dominant contribution from nearby sources. In this work, we further investigate the parameter space of sources allowed by the observational energy spectra and anisotropy amplitude. To obtain the best-fit source parameters, a numerical package to compute the parameter posterior distributions based on Bayesian inference, which is applied to perform an elaborate scan of parameter space. We find that by combining the energy spectra and anisotropy data, the permissible range of location and age of local source is considerably reduced. When comparing with the current local SNR catalog, only Geminga SNR could be the proper candidate of the local cosmic-ray source., Comment: 9 pages, 9 figures, 2 tables

The Tibet ASgamma experiment just reported their measurement of sub-PeV diffuse gamma ray emission from the Galactic disk, with the highest energy up to 957 TeV. These gamma-rays are most likely the hadronic origin by cosmic ray interaction with interstellar gas in the Galaxy. This measurement provides direct evidence to the hypothesis that the Galactic cosmic rays can be accelerated beyond PeV energies. In this work, we try to explain the sub-PeV diffuse gamma-ray spectrum within cosmic rays diffusive propagation model. We find there is a tension between the sub-PeV diffuse gamma rays and the local cosmic ray spectrum. To describe the sub-PeV diffuse gamma-ray flux, it generally requires larger local cosmic-ray flux than measurement in the knee region. We further calculate the PeV neutrino flux from the cosmic ray propagation model. Even all of these sub-PeV diffuse gamma rays originate from the propagation, the Galactic neutrinos only account for less than ~15% of observed flux, most of which are still from extragalactic sources., Comment: 8 pages, 5 figures

The observation of very-high-energy (VHE, $\rm >10~GeV$) $\gamma$-ray emission from $\rm \gamma$-ray bursts (GRBs), especially in the prompt phase, will provide critical information for understanding many aspects of their nature including the physical environment, the relativistic bulk motion, the mechanisms of particle acceleration of GRBs and for studying Lorentz invariance violation, etc. For the afterglow phase, the highest energy photons detected to date by the imaging atmospheric Cherenkov telescopes extend to the TeV regime. However, for the prompt phase, years of efforts in searching for the VHE emission has yielded no statistically significant detections. A wide field-of-view (FOV) and large effective area above tens of GeV are essential for detecting the VHE emissions from GRBs in the prompt phase. The High Altitude Detection of Astronomical Radiation (HADAR) experiment has such merits. In this paper, we report the estimates of its expected annual GRB detection rate, which are obtained by combining the performance of the HADAR instrument with the theoretical calculations based on a phenomenological model to generate the pseudo-GRB population. The expected detectable gamma-ray signal from GRBs above the background is then obtained to give the detection rate. In the spectral model, an extra component is assigned to every GRB event in addition to the Band function. The results indicate that if the energy of the cutoff due to internal absorption is higher than 50 GeV, the detection rate for GRBs for the HADAR experiment is approximately two or three GRBs per year, which varies slightly depending upon the characteristics of the extra component., Comment: accepted by ApJ

As a basic characteristic of cosmic ray (CR) propagation, the diffusive halo can advance our understanding of many CR-related studies and indirect dark matter. The method to derive the halo size usually has degeneracy problems thus affected by large uncertainties. The diffusion gamma ray from high-latitude clouds might shed light on the halo size independently. Since the spatially dependent propagation (SDP) model has a better agreement with the observed CRs, compared with conventional propagation model, in this work, we investigate the halo thickness based on SDP model with Fermi-LAT $\rm\gamma$-ray observation on the high- and intermediate-velocity clouds. As a result, in order not to exceed the relative $\gamma$-ray emissivity in the high-latitude clouds, halo thickness should be in the range of $\rm ~3.3\sim9~ kpc$. Moreover, the spatial morphology of $\rm\gamma$-rays estimated based on SDP model under different values of halo thickness are distinctive, which provides us a tool to determine the halo size. We hope that our model could be tested and tuned by multi-wavelength observations in the future., Comment: 5 pages, 3 figures

The spatial-dependent propagation (SDP) model with a nearby source works well to reproduce the co-evolving features of both cosmic ray (CR) nuclei spectra and anisotropy. However, it is well known that the Sun is offset from the Galactic plane. This will lead to a dominating anisotropy in perpendicular direction, which is discrepant with observations. Thus it is necessary to reboot further investigation into the effect of the solar offset. In this work, for the first time the combined studies of the solar offset, nuclei spectra and anisotropy are performed based on the SDP model. As a result, to reproduce CR spectra and anisotropy, the thickness of the halo $\rm (\xi z_h)$ increases linearly with the displacement of the Sun. When the offset is $\rm \sim8~pc$ as estimated from the matter-borne methods, $\rm \xi z_h$ is about 0.9 kpc, which is a much thicker halo than usually. Moreover, the PeV anisotropy could estimate the value of diffusion coefficient, thus breaking the degeneracy of diffusion coefficient and halo thickness. Therefore it is a good messenger to constrain the halo thickness. On the other hand, the anisotropy in PeV energy region, as a new probe, might also shed new light to constrain the solar offset. We hope that the anisotropy at the energies of $\rm \sim TeV$ to $\rm PeV$ can be finely measured by LHAASO experiment, leading to a better understanding about the thick halo., Comment: 7 pages, 7 figures

The local pulsar and its progenitor, SNR, can together accelerate the positron, electron and nuclei to very high energy. The famous excesses of positron(nuclei) above $20$($200$) GeV possibly come from such kind of local source. This hints that the primary electron should also hold "excess" above $200$ GeV, synchronously accelerated along with the nuclei. The recent precise measurement of sharp dropoff at 284 GeV of positron by AMS-02 experiment takes chance to study this expected electron excess. In this work, the spatially-dependent propagation with a local source is used to reproduce the spectrum of positron, electron and proton. When considering the dropoff at 284 GeV of positron, a sharp bump structure for primary electron above 284 GeV is required to fit the total spectrum of positron and electron. Then we systematically study the common origin of the excesses of positron, electron and nuclei from Geminga pulsar and SNR. Those excesses can be reproduced under this unified single-source model. Lastly, we hope that the fine bump structure can be observed to support our model by AMS-02 experiment in future., Comment: 15 pages, 4 figures published by JCAP

Recently the AMS-02 reported the precise measurements of the energy spectra of medium-mass compositions (Neon, Magnesium, Silicon) of primary cosmic rays, which reveal different properties from those of light compositions (Helium, Carbon, Oxygen). Here we propose a nearby source scenario, together with the background source contribution, to explain the newly measured spectra of cosmic ray Ne, Mg, Si, and particularly their differences from that of He, C, O. Their differences at high energies can be naturally accounted for by the element abundance of the nearby source. Specifically, the abundance ratio of the nearby source to the background of the Ne, Mg, Si elements is lower by a factor of $\sim1.7$ than that of the He, C, O elements. Such a difference could be due to the abundance difference of the stellar evolution of the progenitor star or the acceleration process/environment, of the nearby source. This scenario can simultaneously explain the high-energy spectral softening features of cosmic ray spectra revealed recently by CREAM/NUCLEON/DAMPE, as well as the energy-dependent behaviors of the large-scale anisotropies. It is predicted that the dipole anisotropy amplitudes below PeV energies of the Ne, Mg, Si group are smaller than that of the He, C, O group, which can be tested with future measurements., Comment: 6 pages, 3 figures; accepted for publication in Frontiers of Physics

Conventional cosmic-ray propagation models usually assume an isotropic diffusion coefficient to account for the random deflection of cosmic rays by the turbulent interstellar magnetic field. Such a picture is very successful in explaining many observational phenomena related to the propagation of Galactic cosmic rays, such as broken power-law energy spectra, secondary-to-primary ratios, etc. However, the isotropic diffusion presupposition is facing severe challenges from recent observations. In particular, such observations on the large-scale anisotropy of TeV cosmic rays show that the dipole direction differs from the prediction of the conventional model. One possible reason is that the large-scale regular magnetic field, which leads to an anisotropic diffusion of cosmic rays, has not been included in the model provided by the public numerical packages. In this work, we propose two numerical schemes to solve the $3$-dimensional anisotropic transport equation: the pseudo source method and Hundsdorfer-Verwer scheme. Both methods are verified by reproducing the measured B/C and proton spectrum and the radial variation of spectral index expected by former 2D simulation. As a demonstration of the prediction capability, dipole anisotropy is also calculated by a toy simulation with a rough magnetic field., Comment: 11 pages, 5 figures

The spectral hardenings of cosmic ray nuclei above $\sim 200$ GV followed by softenings around 10 TV, the knee of the all-particle spectrum around PeV energies, as well as the pattern change of the amplitude and phase of the large-scale anisotropies around 100 TeV indicate the complexities of the origin and transportation of Galactic cosmic rays. It has been shown that nearby source(s) are most likely to be the cause of such spectral features of both the spectra and the anisotropies. In this work, we study the anisotropy features of different mass composition (or mass groups) of cosmic rays in this nearby source model. We show that even if the spectral features from the nearby source component is less distinctive compared with the background component from e.g., the population of distant sources, the anisotropy features are more remarkable to be identified. Measurements of the anisotropies of each mass composition (group) of cosmic rays by the space experiments such as DAMPE and HERD and the ground-based experiments such as LHAASO in the near future are expected to be able to critically test this scenario., Comment: 9 pages, 5 figures

Since the science white paper of the Large High Altitude Air Shower Observatory (LHAASO) published on arXiv in 2019 [e-Print: 1905.02773 (astro-ph.HE)], LHAASO has completed the transition from a project to an operational gamma-ray astronomical observatory LHAASO is a new generation multi-component facility located in Daocheng, Sichuan province of China, at an altitude of 4410 meters. It aims at measuring with unprecedented sensitivity the spectrum, composition, and anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$~eV, and acting simultaneously as a wide aperture (one stereoradiant) continuously operating gamma-ray telescope in the energy range between 10$^{11}$ and $10^{15}$~eV with the designed sensitivity of 1.3\% of the Crab Unit (CU) above 100 TeV. LHAASO's capability of measuring simultaneously different shower components (electrons, muons, and Cherenkov/fluorescence light), will allow it to investigate the origin, acceleration, and propagation of CR through measurement of the energy spectrum, elemental composition, and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO will play a key role in CR physics and gamma-ray astronomy for a general and comprehensive exploration of the high energy universe and will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. The LHAASO Collaboration organized an editorial working group and finished all editorial work of this science book, to summarize the instrumental features and outline the prospects of scientific researches with the LHAASO experiment., Comment: This document is a collaborative effort, 185 pages, 110 figures

The spatial-dependent propagation model has been successfully used to explain diverse observational phenomena, including the spectral hardening of cosmic-ray nuclei above $200$ GV, the large-scale dipole anisotropy and the diffusive gamma distribution. In this work, we further apply the spatial-dependent propagation model to both electrons and positrons. To account for the excess of positrons above $10$ GeV, an additional local source is introduced. And we also consider a more realistic spiral distribution of background sources. We find that due to the gradual hardening above $10$ GeV, the hardening of electron spectrum above tens of GeV can be explained in the SDP model and both positron and electron spectra less than TeV energies could be naturally described. The spatial-dependent propagation with spiral-distributed sources could conforms with the total electron spectrum in the whole energy. Meanwhile compared with the conventional model, the spatial-dependent propagation with spiral-distributed sources could produce larger background positron flux, so that the multiplier of background positron flux is $1.42$, which is much smaller than the required value by the conventional model. Thus the shortage of background positron flux could be solved. Furthermore we compute the anisotropy of electron under spatial-dependent propagation model, which is well below the observational limit of Fermi-LAT experiment.

The origin of Galactic cosmic rays (GCRs) remains a mystery after more than one century of their discovery. The diffusive propagation of charged particles in the turbulent Galactic magnetic field makes us unable to trace back to their acceleration sites. Nevertheless, nearby GCR source(s) may leave imprints on the locally measured energy spectra and the anisotropies of the arrival direction. In this work we propose a simple but natural description of the GCR production and propagation, within a two-zone disk-halo diffusion scenario together with a nearby source, to understand the up-to-date precise measurements of the energy spectra and anisotropies of GCRs. We find that a common energy scale of $\sim100$ TeV appears in both energy spectra of protons and Helium nuclei measured recently by CREAM and large-scale anisotropies detected by various experiments. These results indicate that one or more local sources are very likely important contributors to GCRs below $100$ TeV. This study provides a probe to identify source(s) of GCRs by means of joint efforts of spectral and anisotropy measurements., Comment: 17 pages, 5 figures

Recently AMS-02 collaboration publish their measurements of light cosmic-ray nuclei, including lithium, beryllium, boron, carbon and oxygen. All of them reveal a prominent excess above $\sim$ 200 GV, coinciding with proton and helium. Particularly, the secondary cosmic rays even harden than the primary components above that break. One of the viable interpretations for above anomalies is the spatial-dependent diffusion process. Such model has been successfully applied to multiple observational phenomena, for example primary cosmic ray nuclei, diffuse gamma ray and anisotropy. In this work, we investigate the spatial-dependent propagation model in light of the new observational data. We find that such model is able to explain the upturn of secondary spectrum as well as the primary's. All the spectra can be well reproduced and the calculated ratios are also in good agreement with the observations., Comment: 8 pages, 5 figures

Recent direct measurements of Galactic cosmic ray spectra by balloon/space-borne detectors reveal spectral hardenings of all major nucleus species at rigidities of a few hundred GV. The all-sky diffuse gamma-ray emissions measured by the Fermi Large Area Telescope also show spatial variations of the intensities and spectral indices of cosmic rays. These new observations challenge the traditional simple acceleration and/or propagation scenario of Galactic cosmic rays. In this work we propose a spatially-dependent diffusion scenario to explain all these phenomena. The diffusion coefficient is assumed to be anti-correlated with the source distribution, which is a natural expectation from the charged particle transportation in turbulent magnetic field. The spatially-dependent diffusion model also gives a lower level of anisotropies of cosmic rays, which are consistent with observations by underground muons and air shower experiments. The spectral variations of cosmic rays across the Galaxy can be properly reproduced by this model., Comment: accepted by PRD

The DArk Matter Particle Explorer (DAMPE), a high energy cosmic ray and $\gamma$-ray detector in space, has recently reported the new measurement of the total electron plus positron flux between 25 GeV and 4.6 TeV. A spectral softening at $\sim0.9$ TeV and a tentative peak at $\sim1.4$ TeV have been reported. We study the physical implications of the DAMPE data in this work. The presence of the spectral break significantly tightens the constraints on the model parameters to explain the electron/positron excesses. The spectral softening can either be explained by the maximum acceleration limits of electrons by astrophysical sources, or a breakdown of the common assumption of continuous distribution of electron sources at TeV energies in space and time. The tentive peak at $\sim1.4$ TeV implies local sources of electrons/positrons with quasi-monochromatic injection spectrum. We find that the cold, ultra-relativistic $e^+e^-$ winds from pulsars may give rise to such a structure. The pulsar is requird to be middle-aged, relatively slowly-rotated, mildly magnetized, and isolated in a density cavity. The annihilation of DM particles ($m_{\chi}\sim1.5$ TeV) into $e^+e^-$ pairs in a nearby clump or an over-density region may also explain the data. In the DM scenario, the inferred clump mass (or density enhancement) is about $10^7-10^8$ M$_\odot$ (or $17-35$ times of the canonical local density) assuming a thermal production cross section, which is relatively extreme compared with the expectation from numerical simulations. A moderate enhancement of the annihilation cross section via, e.g., the Sommerfeld mechanism or non-thermal production, is thus needed., Comment: 16 pages, 11 figures, and 3 tables

More than fifty years after the discovery of the knee in the cosmic ray spectra, its physical origin remains a mystery. This is partly due to the ambiguity of the energy spectrum of individual composition. Recently, direct measurements from several balloon/satellite-borne detectors found significant spectral hardenings of cosmic ray nuclei at a few hundred GV rigidities. A joint modeling of the direct and indirect measurements of the cosmic ray spectra may help to understand the experimental systematics and probably the physics of the knee. In this work, we update the phenomenological "poly-gonato" model of the knee to include the spectral hardenings, with a changing spectral index of $\gamma + \beta \cdot \log E$ with respect to energy. We find that there is some tension between the direct and indirect measurements. Specifically, the hard spectra of CREAM suggest a relatively low energy cutoff of the individual cosmic ray spectrum around the knee, which is inconsistent with the light component spectra of the air shower measurements. Furthermore, the light component data from AS$\gamma$/ARGO-YBJ/WFCTA/KASCADE are consistent with the all-particle spectra. We expect that more precise measurements of the cosmic ray spectrum of individual composition at TeV energies and beyond by the operating space detectors such as CALET and DAMPE, and the ground based experiments such as LHAASO, will eventually solve these discrepancies. Finally, as an illustration, we show that a spatial-dependent diffusion under two-halo model is able to reproduce the results of this "modified poly-gonato" model.

Spectra of Cosmic Rays (CRs), particularly their features, may hold essential information about astroparticle and fundamental physics. Realizing that both $4$ PeV CR knee and electron's TeV break have nearly identical Lorentz factor, i.e. $\gamma \sim 10^6$, we propose that both CR nuclei and electron experience a threshold interaction with unknown light particle $\rm X$, widespread in the diffusive halo. The interaction occurs only when the Lorentz factor of CRs is larger than $10^6$. Via this process, the energy spectrum of each component has deformation above the threshold Lorentz factor, so both knee and spectral break of electrons can be naturally reproduced. Given the large uncertainties of the available measurements and simplified assumption, our model accommodates a wide mass range of the particle X, which is allowable from ultralight value to $\sim 1$ eV., Comment: 7 pages, 7 figures

The multi-TeV $\gamma$-rays from the Galactic Center (GC) have a cutoff at tens of TeV, whereas the diffuse emission has no such cutoff, which is regarded as an indication of PeV proton acceleration by the HESS experiment. It is important to understand the inconsistency and study the possibility that PeV cosmic-ray acceleration could account for the apparently contradictory point and diffuse $\gamma$-ray spectra. In this work, we propose that the cosmic rays are accelerated up to $>$PeV in GC. The interaction between cosmic rays and molecular clouds is responsible for the multi-TeV $\gamma$-ray emissions from both the point source and diffuse sources today. Enhanced by the small volume filling factor (VFF) of the clumpy structure, the absorption of the $\gamma$-rays leads to a sharp cutoff spectrum at tens of TeV produced in the GC. Away from galactic center, the VFF grows and the absorption enhancement becomes negligible. As a result, the spectra of $\gamma$-ray emissions for both point source and diffuse sources can be successfully reproduced under such self-consistent picture. In addition, a "surviving-tail" at $\sim$100 TeV is expected from the point source, which can be observed by future projects CTA and LHAASO. Neutrinos are simultaneously produced during proton-proton (PP) collision. With 5-10 years observations, the KM3Net experiment will be able to detect the PeV source according to our calculation.

Recent precise measurements of cosmic ray spectral revealed an anomalous hardening at ~200 GV for nuclei from PAMELA, CREAM, ATIC, AMS02 experiments and at tens of GeV for primary electron derived from AMS02 experiment. Particularly, the latest observation of pbar/p ratio by AMS02 demonstrated a flat distribution, which further validated the spectrum anomalies of secondary particles. All those new phenomena indicated that the conventional propagation model of cosmic rays meet challenge. In this work, the spatial-dependent propagation coefficient D(r,z,\rho) is employed by tracing the source distribution under the physical picture of two-halo model in DRAGON package. Under such scenario, the model calculation will result in a two-component spectral for primary nuclei and electron. Simultaneously, due to the smaller rigidity dependence of D(r,z,\rho) in galactic disk, the ratio of secondary-to-primary will be inevitablly flatter than the calculation in the conventional propagation model. As a result, we can reproduce the spectral hardening of proton, electron and the flat ratio of pbar/p and B/C by only adopting the spatial-dependent propagation coefficient D(r,z,\rho) in galactic disk.