Search

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 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.

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

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

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

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