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.

Anisotropic diffusion is one of the potential interpretations for the morphology of the Geminga pulsar halo. It interprets the observed slow-diffusion phenomenon through a geometric effect, assuming the mean magnetic field direction around Geminga is closely aligned with the line of sight toward it. However, this direction should not extend further than the correlation length of the turbulent magnetic field $L_c$, which could be $100$ pc or less. We first revisit the $L_c=\infty$ scenario and show that the halo asymmetry predicted by this scenario is mainly contributed by the electrons located beyond the ``core" section around Geminga, which has a length of $100$ pc. Then, considering the directional variation of the magnetic field beyond the core section, we take one magnetic field configuration as an example to investigate the possible halo morphology. The predicted morphology has some different features compared to the $L_c=\infty$ scenario. The current experiments may already be able to test these features. In addition, we use a semi-analytical method to solve the anisotropic propagation equation, which offers significant convenience compared to numerical approaches., Comment: 15 pages, 7 figures

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

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 discrepancy between the theoretical prediction of primordial lithium abundances and astronomical observations is called the Lithium Problem. We find that extra contributions from non-thermal hydrogen and helium during Big Bang nucleosynthesis can explain the discrepancy, for both Li-7 and Li-6, and will change the deuterium abundance only little. The allowed parameter space of such an amount of non-thermal particles and the energy range is shown. The hypothesis is stable regardless of the cross-section uncertainty of relevant reactions and the explicit shape of the energy spectrum., Comment: 20 pages, 9 figures

The remnant of supernova explosion is widely believed to be the acceleration site of high-energy cosmic ray particles. The acceleration timescale is, however, typically very long. Here we report the detection of a variable $\gamma$-ray source with the Fermi Large Area Telescope, which is positionally and temporally consistent with a peculiar supernova, iPTF14hls. A quasi-stellar object SDSS J092054.04+504251.5, which is probably a blazar according to the infrared data, is found in the error circle of the $\gamma$-ray source. More data about the $\gamma$-ray source and SDSS J092054.04+504251.5 are needed to confirm their association. On the other hand, if the association between the $\gamma$-ray source and the supernova is confirmed, this would be the first time to detect high-energy $\gamma$-ray emission from a supernova, suggesting very fast particle acceleration by supernova explosions., Comment: 12 pages, 5 figures; accepted for publication in ApJL

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

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

There may be non-thermal cosmic rays during big-bang nucleosynthesis (BBN) epoch (dubbed as BBNCRs). This paper investigated whether such BBNCRs can be the origin of Lithium problem or not. It can be expected that BBNCRs flux will be small in order to keep the success of standard BBN (SBBN). With favorable assumptions on the BBNCR spectrum between 0.09 -- 4 MeV, our numerical calculation showed that extra contributions from BBNCRs can account for the $^7$Li abundance successfully. However $^6$Li abundance is only lifted an order of magnitude, which is still much lower than the observed value. As the deuteron abundance is very sensitive to the spectrum choice of BBNCRs, the allowed parameter space for the spectrum is strictly constrained. We should emphasize that the acceleration mechanism for BBNCRs in the early universe is still an open question. For example, strong turbulent magnetic field is probably the solution to the problem. Whether such a mechanism can provide the required spectrum deserves further studies., Comment: 34 pages, 21 figures, published version

The anisotropy of cosmic rays (CRs) in the solar vicinity is generally at- tributed to the CR streaming due to the discrete distribution of CR sources or local magnetic field modulation. Recently, the two dimensional large scale CR anisotropy has been measured by many experiments in TeV-PeV energy range in both hemispheres. The tail-in excess along the tangential direction of the local spiral arm and the loss cone deficit pointing to the north Galactic pole direction agree with what have been obtained in tens to hundreds of GeV. The persistence of the two large scale anisotropy structures in such a wide range of energy suggests that the anisotropy might be due to a global streaming of the Galactic CRs (GCRs). This work tries to extend the observed CR anisotropy picture from solar system to the whole galaxy. In such a case, we can find a new interesting signature, a loop of GCR streaming, of the GCR propagation. We further calculate the overall GCR streaming induced magnetic field, and find a qualitative consistence with the observed structure of the halo magnetic field.

Supernova remnants have long been regarded as sources of the Galactic cosmic rays up to petaelectronvolts, but convincing evidence is still lacking. In this work we explore the common origin of the subtle features of the cosmic ray spectra, such as the knee of cosmic ray spectra and the excesses of electron/positron fluxes recently observed by ATIC, H.E.S.S., Fermi-LAT and PAMELA. Numerical calculation shows that those features of cosmic ray spectra can be well reproduced in a scenario with e$^+$e$^-$ pair production by interactions between high energy cosmic rays and background photons in an environment similar to the young supernova remnant. The success of such a coherent explanation serves in turn as an evidence that at least a portion of cosmic rays might be accelerated at young supernova remnants., Comment: 11 pages, 3 figures and 1 table. Accepted for publication in ApJL

The CACTUS experiment recently observed a gamma ray excess above 50 GeV from the direction of the Draco dwarf spheroidal galaxy. Considering that Draco is dark matter dominated the gamma rays may be generated through dark matter annihilation in the Draco halo. In the framework of the minimal supersymmetric extension of the standard model we explore the parameter space to account for the gamma ray signals at CACTUS. We find that the neutralino mass is constrained to be approximately in the range between 100 GeV ~ 400 GeV and a sharp central cuspy of the dark halo profile in Draco is necessary to explain the CACTUS results. We then discuss further constraints on the supersymmetric parameter space by observations at the ground based ARGO detector. It is found that the parameter space can be strongly constrained by ARGO if no excess from Draco is observed above 100 GeV., Comment: 15 pages, 4 figures

Adopting a recent calculation of the Galactic interstellar radiation field, we calculate the attenuation of the very high energy gamma rays from the Galactic sources. The infra-red radiation background near the Galactic Center is very intense due to the new calculation and our result shows that a cutoff of high energy gamma ray spectrum begins at about 20 TeV and reaches about 10% for 50 TeV gamma rays., Comment: 6 pages, 1 figure, figure is changed, conclusion not changed