Search

The Very Large Area gamma-ray Space Telescope (VLAST) is a mission concept proposed to detect gamma-ray photons through both the Compton scattering and electron-positron pair production mechanisms, enabling the detection of photons with energies ranging from MeV to TeV. This project aims to conduct a comprehensive survey of the gamma-ray sky from a low Earth orbit using an anti-coincidence detector, a tracker detector that also serves as a low energy calorimeter, and a high energy imaging calorimeter. We developed a Monte Carlo simulation application of the detector with the GEANT4 toolkit to evaluate the instrument performance including the effective area, angular resolution and energy resolution, as well as explored specific optimizations of the detector configuration. Our simulation-based analysis indicates that the VLAST's current design is physically feasible, with an acceptance larger than 10~$\rm m^2\ sr$ which is four times larger than Fermi-LAT, an energy resolution better than 2\% at 10~GeV, and an angular resolution better than 0.2 degrees at 10~GeV. The VLAST project is expected to make significant contribution to the field of gamma-ray astronomy and to enhance our understanding of the cosmos., Comment: 15 pages, 16 figures; Nuclear Science and Techniques in press

As the largest gravitationally bound objects in the Universe, galaxy clusters have provided the first piece of evidence for the presence of dark matter and may be suitable targets for indirect dark matter searches. Among various signals, the GeV-TeV $\gamma$-ray line has been taken as the smoking-gun signal of the dark matter annihilation/decay since no known astrophysical/physical process(es) could generate such a peculiar spectrum. With 15.5 years of Fermi-LAT P8R3 publicly available data, we search for the $\gamma$-ray line emission in the directions of a group of 13 nearby massive galaxy clusters with an unbinned likelihood analysis. A $\gamma$-ray line signal at $\sim 43.2$ GeV has a net TS value of $\approx 30$ if we only take into account the data in the directions of Virgo, Fornax and Ophiuchus clusters, three massive clusters with the highest J-factors expected to generate the dark matter annihilation signal. The signal still presents when the data of other 10 nearby massive clusters have also been included, though the TS value decreases to $\approx 21$ likely because of their lower signal-to-noise ratios. The absence of this signal in the inner Galaxy disfavors both the instrumental effect and the canonical dark matter annihilation interpretation, and a more sophisticated dark matter model or very peculiar astrophysical scenario might be needed. This $\gamma$-ray line signal, if intrinsic, could be unambiguously verified by the Very Large Area $\gamma$-ray Space Telescope in its first two years of performance., Comment: 9 pages, 7 figures, 1 table. Comments are welcome!

The extended $\gamma$-ray halos around pulsars are unique probe of transportation of high-energy electrons (and positrons) in vicinities of such pulsars. Observations of morphologies of several such halos indicate that particles diffuse very slowly around pulsars, compared with that in the Milky Way halo. The energy-dependent morphologies are expected to be very important in studying the energy-dependence of the diffusion coefficient. In this work we point out that the spectrum of high-energy electrons takes effect in shaping the $\gamma$-ray morphologies at the ultra-high-energy bands, and thus results in a degeneracy between the electron spectrum and the energy-dependence of the diffusion coefficient. The reasons for such a degeneracy include both the Klein-Nishina effect of the inverse Compton scattering and the curvature (if any) of the electron spectrum. It it thus necessary to take into account the spectral shape of electrons when deriving the energy-dependence of diffusion coefficient using ultra-high-energy $\gamma$-ray measurements of extended pulsar halos.

The recent measurements of cosmic ray deuteron fluxes by AMS-02 show that the rigidity dependence of deuterons is similar with that of protons but flatter than $^3$He, which has been attributed to the existence of primary deuterons with abundance much higher than that from the Big Bang nucleosynthesis. The requirement of highly deuteron-abundant sources imposes a serious challenge on the modern astrophysics since there is no known process to produce a large amount of deuterons without violating other constraints \cite{1976Natur.263..198E}. In this work we demonstrate that the fragmentation of heavy nuclei up to iron plays a crucial role in shaping/enhancing the spectrum/flux of the cosmic ray deuterons. Based on the latest cosmic ray data, the predicted secondary fluxes of deuterons and $^3$He are found to be reasonably consistent with the AMS-02 measurements and a primary deuteron component is not needed. The apparent difference between D/$^4$He (D/p) and $^3$He/$^4$He ($^3$He/p) is probably due to a combined effect of the kinetic-energy-to-rigidity conversion and the solar modulation. More precise measurements of the fragmentation cross sections of various nuclei to produce deuterons, tritons, and $^3$He in a wide energy range will be very helpful in further testing the secondary origin of cosmic ray deuterons., Comment: 5 pages, 5 figures, and 1 table

Strongly scattering media disrupt both the wavefront distribution and the polarization state of the incident light field. Controlling and effectively utilizing depolarization effects are crucial for optical applications in highly scattering environments, such as imaging through dense fog. However, current simulation models have difficulty simulating the evolution of vector light fields within scattering media, posing challenges for studying vector light fields in strongly scattering environments. Here, we propose the Vector Angular Spectrum (VAS) model for simulating the propagation of vector light fields within scattering media. By introducing the angular spectrum distribution of vector light scattering and polarization conversion mechanisms, this model can simulate the depolarization effects of vector light propagating through strongly scattering media. The VAS model has also been used to investigate the focusing of vector scattered light through scattering media. Furthermore, the simulation results of the model have been validated through experiments. The proposed VAS model is expected to play a role in the theoretical research of vector scattered light and optical applications in strongly scattering environments.

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

Precise measurements of energy spectra of different cosmic ray species were obtained in recent years, by particularly the AMS-02 experiment on the International Space Station. It has been shown that apparent differences exist in different groups of the primary cosmic rays. However, it is not straightforward to conclude that the source spectra of different particle groups are different since they will experience different propagation processes (e.g., energy losses and fragmentations) either. In this work, we study the injection spectra of different nuclear species using the measurements from Voyager-1 outside the solar system, and ACR-CRIS and AMS-02 on top of atmosphere, in a physical framework of cosmic ray transportation. Two types of injection spectra are assumed, the broken power-law and the non-parametric spline interpolation form. The non-parametric form fits the data better than the broken power-law form, implying that potential structures beyond the constrained spectral shape of broken power-law may exist. For different nuclei the injection spectra are overall similar in shape but do show some differences among each other. For the non-parametric spectral form, the helium injection spectrum is the softest at low energies and the hardest at high energies. For both spectral shapes, the low-energy injection spectrum of neon is the hardest among all these species, and the carbon and oxygen spectra have more prominent bumps in 1-10 GV in the R2dN/dR presentation. Such differences suggest the existence of differences in the sources or acceleration processes of various nuclei of cosmic rays.

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.

We investigate the impact of inelastic collisions between dark matter (DM) and heavy cosmic ray (CR) nuclei on CR propagation. We approximate the fragmentation cross-sections for DM-CR collisions using collider-measured proton-nuclei scattering cross-sections, allowing us to assess how these collisions affect the spectra of CR Boron and Carbon. We derive new CR spectra from DM-CR collisions by incorporating their cross-sections into the source terms and solving the diffusion equation for the complete network of reactions involved in generating secondary species. In a specific example with a coupling strength of $b_{\chi}=0.1$ and a DM mass of $m_{\chi}=0.1$ GeV, considering a simplified scenario where DM interacts exclusively with Oxygen, a notable modification in the Boron-to-Carbon spectrum due to the DM-CR interaction is observed. Particularly, the peak within the spectrum, spanning from $0.1$ GeV to $10$ GeV, experiences an enhancement of approximately 1.5 times. However, in a more realistic scenario where DM particles interact with all CRs, this peak can be amplified to twice its original value.Utilizing the latest data from AMS-02 and DAMPE on the Boron-to-Carbon ratio, we estimate a 95\% upper limit for the effective inelastic cross-section of DM-proton as a function of DM mass. Our findings reveal that at $m_\chi\simeq 2$ MeV, the effective inelastic cross-section between DM and protons must be less than $\mathcal{O}(10^{-32})~{\rm cm}^2$., Comment: 19 pages, 8 figures

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 axion is a hypothetical elementary particle that could solve the long-standing strong CP problem in particle physics and the dark matter mystery in the cosmos. Due to the stimulation of the ambient photons, the axion dark matter decay into photons is significantly enhanced so that its echo signal could be detected by terrestrial telescopes. As a pathfinder, we study the expected sensitivity of searching for the axion dark matter in the mass range between $0.41$ and $1.6\mu\text{eV}$ with the 21 CentiMeter Array (21CMA). We aim to cover the whole 21CMA frequency range in two years by using a 1MW emitter. We find that the resulting sensitivity on the axion-photon coupling could surpass other existing limits by about one order of magnitude., Comment: 7 pages, 4 figures

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 Giant Radio Array for Neutrino Detection (GRAND) is an envisioned observatory of ultra-high-energy particles of cosmic origin, with energies in excess of 100 PeV. GRAND uses large surface arrays of autonomous radio-detection units to look for the radio emission from extensive air showers that are triggered by the interaction of ultra-high-energy cosmic rays, gamma rays, and neutrinos in the atmosphere or underground. In particular, for ultra-high-energy neutrinos, the future final phase of GRAND aims to be sensitive enough to discover them in spite of their plausibly tiny flux. Presently, three prototype GRAND radio arrays are in operation: GRANDProto300, in China, GRAND@Auger, in Argentina, and GRAND@Nancay, in France. Their goals are to field-test the design of the radio-detection units, understand the radio background to which they are exposed, and develop tools for diagnostic, data gathering, and data analysis. This list of contributions to the 38th International Cosmic Ray Conference (ICRC 2023) presents an overview of GRAND, in its present and future incarnations, and a look at the first data collected by GRANDProto13, the first phase of GRANDProto300., Comment: Note: To access the list of contributions, please follow the "HTML" link that can be found on the arXiv page

The Dark Matter Particle Explorer (DAMPE) is a satellite-borne detector designed to measure high energy cosmic-rays and $\gamma$-rays. As a key sub-detector of DAMPE, the Bismuth Germanium Oxide (BGO) imaging calorimeter is utilized to measure the particle energy with a high resolution. The nonlinear fluorescence response of BGO for large ionization energy deposition, known as the quenching effect, results in an under-estimate of the energy measurement for cosmic-ray nuclei. In this paper, various models are employed to characterize the BGO quenching factors obtained from the experimental data of DAMPE. Applying the proper quenching model in the detector simulation process, we investigate the tuned energy responses for various nuclei and compare the results based on two different simulation softwares, i.e. GEANT4 and FLUKA. The BGO quenching effect results in a decrease of the measured energy by approximately $2.5\%$ ($5.7 \%$) for carbon (iron) at $\sim$10 GeV/n and $<1\%$ above 1 TeV/n, respectively. Accordingly, the correction of the BGO quenching effect leads to an increase of the low-energy flux measurement of cosmic-ray nuclei., Comment: 13 pages, 4 figures, to be published in Nuclear Inst. and Methods in Physics Research, A

The high-energy diffuse gamma-ray emission and neutrino emission are expected from the Galactic plane, generated by hadronuclear interactions between cosmic rays (CR) and interstellar medium (ISM). Therefore, measurements of these diffuse emissions will provide important clues on the origin and nature of Galactic CRs. Comparing the latest observations of LHAASO and IceCube on the diffuse Galactic gamma-ray and neutrino emissions respectively, we suggest that the diffuse gamma-ray emission at multi-TeV energies contains a considerable contribution of a leptonic component. By modelling the gamma-ray halos powered by middle-aged pulsars in our Galaxy with taking into account the magnetic field configuration and the interstellar radiation field in the Galaxy, we demonstrate that the collective contribution of pulsar halos can account for the excess in the measured diffuse gamma-ray emission with respect to the predicted flux from CR-ISM interactions. The resulting one-dimensional profile along the Galactic longitude is also consistent with the observation., Comment: 28 pages, 13 figures; published in Nature Astronomy

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 origin of interstellar and intergalactic magnetic fields remains largely unknown. One possibility is that they are related to the primordial magnetic fields (PMFs) produced by, for instance, the phase transitions of the early Universe. In this paper, we show that the PMF-induced turbulence generated at around the QCD phase transition epoch--the characteristic magnetic field strength $B_{\rm ch}^* \sim \mathcal{O}(1)~\rm{\mu G}$ and coherent length scale $\ell_{\rm ch}^* \sim \mathcal{O}(1)~\rm{pc}$--can naturally accommodate nanohertz gravitational waves reported by pulsar timing array (PTA) collaborations. Moreover, the evolution of the PMFs to the recombination era with the form of $B_{\rm ch}\sim \ell_{\rm ch}^{-\alpha}$ can induce baryon density inhomogeneities, alter the recombination history, and alleviate the tension of the Hubble parameter $H_0$ and the matter clumpiness parameter $S_8$ between early- and late-time measurements for $0.88\leq \alpha \leq 1.17$ (approximate 95\% credible region based on three PTA likelihoods). The further evolved PMFs may account for the $\sim {\cal O}(10^{-16})$ Gauss extragalactic magnetic field inferred with GRB 221009A., Comment: 7 pages, 3 figures, 1 table

The diffuse gamma-ray emission between 10 and 1000 TeV from the Galactic plane was recently measured by the Large High Altitude Air Shower Observatory (LHAASO). These observations will help tremendously in constraining the propagation and interaction of cosmic rays in the Milky Way. Additionally, new measurements of CR spectra reach a very high precision up to 100 TeV energies, revealing multiple spectral structures of various species. In this work, we confront the model prediction of the diffuse gamma-ray emission, based on up-to-date measurements of the local cosmic ray spectra and simplified propagation setup, with the measurements of diffuse gamma-rays. To better constrain the low-energy part of the model, we analyze the 14.6 years of Fermi-LAT data to extract the Galactic diffuse emission between 1 and 500 GeV from the same sky regions of LHAASO, after subtracting the contribution from known sources and the isotropic diffuse gamma-ray background. The joint Fermi-LAT and LHAASO spectra thus cover a very wide energy range from 1 GeV to 1 PeV with small gaps from 0.5 to 10 TeV. Compared with the prediction, we find that clear excesses between several GeV and ~60 TeV of the diffuse emission exist. Possible reasons to explain the excesses may include unresolved sources or more complicated propagation models. We illustrate that an exponential-cutoff-power-law component with an index of -2.40 and cutoff energy of ~30 TeV is able to account for such excesses.

Precise measurements of the spectra of secondary and primary cosmic rays are crucial for understanding the origin and propagation of those energetic particles. The High Energy cosmic-Radiation Detection (HERD) facility onboard China`s Space Station, which is expected to operate in 2027, will push the direct measurements of cosmic ray fluxes precisely up to PeV energies. In this work, we investigate the potential of HERD on studying the propagation of cosmic rays using the measurements of boron, carbon, and oxygen spectra. We find that, compared with the current results, the new HERD measurements can improve the accuracy of the propagation parameters by 8\% to 40\%. The constraints on the injection spectra at high energies will also be improved., Comment: 10 pages, 4 figures

Ultralight dark matter (ULDM) is proposed as a theoretical candidate of dark matter particles with masses of approximately $10^{-22}$ eV. The interactions between ULDM particles and standard model particles would cause variations in pulse arrival times of millisecond pulsars, which means that the pulsar timing array (PTA) can be used to indirectly detect ULDM. In this letter, we use the gamma-ray PTA composed of 29 millisecond pulsars observed by the Fermi Large Area Telescope (Fermi-LAT) to test four ULDM effects, including gravitational effects for generalized ULDM with different Spin-0/1, the fifth-force coupling effect of dark photon, and the modified gravitational effect of the Spin-2 ULDM. The gamma-ray pulsar timing is not affected by the ionized interstellar medium and suffers relatively simple noises, unlike that of the radio band. Our work is the first time that the gamma-ray PTA has been used to search for the ULDM. No significant signals of ULDM are found based on the Fermi-LAT PTA for all four kinds of ULDM models. Constraints on ULDM parameters are set with the 95% confidence level, which provides a complementary check of the nondetection of ULDM for radio PTAs and direct detection experiments., Comment: 11 pages, 4 figures, 2 tables; publish as a Letter in Phys. Rev. D

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.

Precise measurements of the stellar orbits around Sagittarius A* have established the existence of a supermassive black hole (SMBH) at the Galactic centre (GC). Due to the interplay between the SMBH and dark matter (DM), the DM density profile in the innermost region of the Galaxy, which is crucial for the DM indirect detection, is still an open question. Among the most popular models in the literature, the theoretical spike profile proposed by Gondolo and Silk (1999; GS hereafter) is well adopted. In this work, we investigate the DM spike profile using updated data from the Keck and VLT telescopes considering that the presence of such an extended mass component may affect the orbits of the S-stars in the Galactic center. We examine the radius and slope of the generalized NFW spike profile, analyze the Einasto spike, and discuss the influence of DM annihilation on the results. Our findings indicate that an initial slope of $\gamma \gtrsim 0.92$ for the generalized NFW spike profile is ruled out at a 95% confidence level. Additionally, the spike radius $R_{\rm sp}$ larger than 21.5 pc is rejected at 95% probability for the Einasto spike with $\alpha=0.17$, which also contradicts the GS spike model. The constraints with the VLT/GRAVITY upper limits are also projected. Although the GS NFW spike is well constrained by the Keck and VLT observation of S2, an NFW spike with a weak annihilation cusp may still be viable, as long as the DM annihilation cross section satisfies $\left< \sigma v \right> \gtrsim 7.7\times 10^{-27}~{\rm cm^3\,s^{-1}} (m_{\rm DM}/100~{\rm GeV})$ at 95% level., Comment: 12 pages, 10 figures and 1 table. Accepted for publication in MNRAS

Galactic cosmic rays are mostly made up of energetic nuclei, with less than $1\%$ of electrons (and positrons). Precise measurement of the electron and positron component requires a very efficient method to reject the nuclei background, mainly protons. In this work, we develop an unsupervised machine learning method to identify electrons and positrons from cosmic ray protons for the Dark Matter Particle Explorer (DAMPE) experiment. Compared with the supervised learning method used in the DAMPE experiment, this unsupervised method relies solely on real data except for the background estimation process. As a result, it could effectively reduce the uncertainties from simulations. For three energy ranges of electrons and positrons, 80--128 GeV, 350--700 GeV, and 2--5 TeV, the residual background fractions in the electron sample are found to be about (0.45 $\pm$ 0.02)$\%$, (0.52 $\pm$ 0.04)$\%$, and (10.55 $\pm$ 1.80)$\%$, and the background rejection power is about (6.21 $\pm$ 0.03) $\times$ $10^4$, (9.03 $\pm$ 0.05) $\times$ $10^4$, and (3.06 $\pm$ 0.32) $\times$ $10^4$, respectively. This method gives a higher background rejection power in all energy ranges than the traditional morphological parameterization method and reaches comparable background rejection performance compared with supervised machine learning~methods., Comment: 10 pages, 5 figures, 1 table

Large High Altitude Air Shower Observatory has detected $0.2-13$ TeV emission of GRB 221009A within 2000 s since the trigger. Here we report the detection of a 400 GeV photon, without accompanying prominent low-energy emission, by Fermi Large Area Telescope in this direction with a 0.4 days' delay. Given an intergalactic magnetic field strength of about $4 \times 10^{-17}$ G, which is comparable to limits from TeV blazars, the delayed 400 GeV photon can be explained as the cascade emission of about 10 TeV gamma rays. We estimate the probabilities of the cascade emission that can result in one detectable photon beyond 100 GeV by Fermi Large Area Telescope within $0.3-1$ days is about 2$\%$ whereas it is about 20.5$\%$ within $0.3-250$ days. Our results show that Synchrotron Self-Compton explanation is less favored with probabilities lower by a factor of about $3-30$ than the cascade scenario., Comment: 9 pages, 4 figures, Nature Communications, in press

Precise measurements of the boron-to-carbon and boron-to-oxygen ratios by DAMPE show clear hardenings around $100$ GeV/n, which provide important implications on the production, propagation, and interaction of Galactic cosmic rays. In this work we investigate a number of models proposed in literature in light of the DAMPE findings. These models can roughly be classified into two classes, driven by propagation effects or by source ones. Among these models discussed, we find that the re-acceleration of cosmic rays, during their propagation, by random magnetohydrodynamic waves may not reproduce sufficient hardenings of B/C and B/O, and an additional spectral break of the diffusion coefficient is required. The other models can properly explain the hardenings of the ratios. However, depending on simplifications assumed, the models differ in their quality in reproducing the data in a wide energy range. The models with significant re-acceleration effect will under-predict low-energy antiprotons but over-predict low-energy positrons, and the models with secondary production at sources over-predict high-energy antiprotons. For all models high-energy positron excess exists., Comment: 14 pages, 7 figures. Published version

The High Altitude Detection of Astronomical Radiation (HADAR) experiment is a refracting terrestrial telescope array based on the atmospheric Cherenkov imaging technique. It focuses the Cherenkov light emitted by extensive air showers through a large aperture water-lens system for observing very-high-energy-rays and cosmic rays. With the advantages of a large field-of-view (FOV) and low energy threshold, the HADAR experiment operates in a large-scale sky scanning mode to observe galactic sources. This study presents the prospects of using the HADAR experiment for the sky survey of TeV {\gamma}-ray sources from TeVCat and provids a one-year survey of statistical significance. Results from the simulation show that a total of 23 galactic point sources, including five supernova remnant sources and superbubbles, four pulsar wind nebula sources, and 14 unidentified sources, were detected in the HADAR FOV with a significance greater than 5 standard deviations ({\sigma}). The statistical significance for the Crab Nebula during one year of operation reached 346.0 {\sigma} and the one-year integral sensitivity of HADAR above 1TeV was ~1.3%-2.4% of the flux from the Crab Nebula.

Axions provide a natural and well-motivated dark matter candidate, with the capability to convert directly to photons in the presence of an electromagnetic field. A particularly compelling observational target is the conversion of dark matter axions into photons in the magnetospheres of highly magnetised neutron stars, which is expected to produce a narrow spectral peak centred at the frequency of the axion mass. We point the MeerKAT radio telescope towards the isolated neutron star J0806.4$-$4123 for $10$-hours of observation and obtain the radio spectra in the frequency range $769$-$1051$ MHz. By modelling the conversion process of infalling axion dark matter (DM), we then compare these spectra to theoretical expectations for a given choice of axion parameters. Whilst finding no signal above $5\sigma$ in the data, we provide a unique constraint on the Primakoff coupling of axion DM, $g_{{\rm a}\gamma\gamma}\lesssim 9.3 \times 10^{-12}\,{\rm GeV}^{-1}$ at the $95\%$ confidence level, in the mass range $3.18$-$4.35\,\mu$eV. This result serves the strongest constraint in the axion mass range $4.20$-$4.35\,\mu$eV., Comment: 6 pages, 3 figures, accepted by Physical Review D

An axion cloud surrounding a supermassive black hole can be naturally produced through the superradiance process. Its existence can be examined by the axion induced birefringence effect. It predicts an oscillation of the electric vector position angle of linearly polarized radiations. Stringent constraints of the existence of the axion in a particular mass window has been obtained based on the recent Event Horizon Telescope measurement on M87$^\star$. The future Very-Long-Baseline Interferometry (VLBI) observations will be able to measure the vicinity of many supermassive black holes, thus it opens the possibility to search for the existence of axions in a wide mass regime. In this paper, we study how different black hole properties and accretion flows influence the signatures of the axion induced birefringence. We include the impacts of black hole inclination angles, spins, magnetic fields, plasma velocity distributions, the thickness of the accretion flows. We pay special attention to characterize the washout effects induced by the finite thickness of the accretion flows and the lensed photons. Based on this study, we give prospects on how to optimize the axion search using future VLBI observations, such as the next-generation Event Horizon Telescope, to further increase the sensitivity., Comment: 36 pages, 13 figures, published version in JCAP

We search for stochastic gravitational wave background emitted from cosmic strings using the Parkes Pulsar Timing Array data over 15 years. While we find that the common power-law excess revealed by several pulsar timing array experiments might be accounted for by the gravitational wave background from cosmic strings, the lack of the characteristic Hellings-Downs correlation cannot establish its physical origin yet. The constraints on the cosmic string model parameters are thus derived with conservative assumption that the common power-law excess is due to unknown background. Two representative cosmic string models with different loop distribution functions are considered. We obtain constraints on the dimensionless string tension parameter $G\mu<10^{-11}\sim10^{-10}$, which is more stringent by two orders of magnitude than that obtained by the high-frequency LIGO-Virgo experiment for one model, and less stringent for the other. The results provide the chance to test the Grand unified theories, with the spontaneous symmetry breaking scale of $U(1)$ being two-to-three orders of magnitude below $10^{16}$ GeV. The pulsar timing array experiments are thus quite complementary to the LIGO-Virgo experiment in probing the cosmic strings and the underlying beyond standard model physics in the early Universe., Comment: 10 pages, 8 figures, 4 tables. Comments welcome