Showing 1–50 of 267 results for author: Yuan, Q

The AMS-02 cosmic ray deuteron flux is consistent with a secondary origin

Authors: Qiang Yuan, Yi-Zhong Fan

Abstract: 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 as… ▽ More 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. △ Less

Submitted 27 June, 2024; originally announced June 2024.

Comments: 5 pages, 5 figures, and 1 table

Constraints on Ultra Heavy Dark Matter Properties from Dwarf Spheroidal Galaxies with LHAASO Observations

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes… ▽ More In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes of astrophysical $γ$-ray background while large amount of dark matter. By analyzing more than 700 days observational data at LHAASO, no significant dark matter signal from 1 TeV to 1 EeV is detected. Accordingly we derive the most stringent constraints on the ultra-heavy dark matter annihilation cross-section up to EeV. The constraints on the lifetime of dark matter in decay mode are also derived. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, accepted by PRL

Data quality control system and long-term performance monitor of the LHAASO-KM2A

Authors: Zhen Cao, F. Aharonian, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, H. X. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen , et al. (263 additional authors not shown)

Abstract: The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To… ▽ More The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To ensure the reliability of the LHAASO-KM2A data, a three-level quality control system has been established. It is used to monitor the status of detector units, stability of reconstructed parameters and the performance of the array based on observations of the Crab Nebula and Moon shadow. This paper will introduce the control system and its application on the LHAASO-KM2A data collected from August 2021 to July 2023. During this period, the pointing and angular resolution of the array were stable. From the observations of the Moon shadow and Crab Nebula, the results achieved using the two methods are consistent with each other. According to the observation of the Crab Nebula at energies from 25 TeV to 100 TeV, the time averaged pointing errors are estimated to be $-0.003^{\circ} \pm 0.005^{\circ}$ and $0.001^{\circ} \pm 0.006^{\circ}$ in the R.A. and Dec directions, respectively. △ Less

Submitted 13 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures

The Intermediate-Mass Black Hole Reverberation Mapping Project: Initial Results for a candidate IMBH in a nearby Seyfert 1 Galaxy

Authors: Wenwen Zuo, Hengxiao Guo, Jingbo Sun, Qi Yuan, Paulina Lira, Minfeng Gu, Philip G. Edwards, Alok C. Gupta, Shubham Kishore, Jamie Stevens, Tao An, Zhen-Yi Cai, Haicheng Feng, Luis C. Ho, Dragana Ilić, Andjelka B. Kovačević, ShaSha Li, Mar Mezcua, Luka Č. Popović, Mouyuan Sun, Tushar Tripathi, Vivian U., Oliver Vince, Jianguo Wang, Junxian Wang , et al. (3 additional authors not shown)

Abstract: To investigate the short-term variability and determine the size of the optical continuum emitting size of intermediate-mass black holes (IMBHs), we carried out high-cadence, multi-band photometric monitoring of a Seyfert 1 galaxy J0249-0815 across two nights, together with a one-night single-band preliminary test. The presence of the broad Ha component in our target was confirmed by recent Paloma… ▽ More To investigate the short-term variability and determine the size of the optical continuum emitting size of intermediate-mass black holes (IMBHs), we carried out high-cadence, multi-band photometric monitoring of a Seyfert 1 galaxy J0249-0815 across two nights, together with a one-night single-band preliminary test. The presence of the broad Ha component in our target was confirmed by recent Palomar/P200 spectroscopic observations, 23 years after Sloan Digital Sky Survey, ruling out the supernovae origin of the broad Ha line. The photometric experiment was primarily conducted utilizing four-channel imagers MuSCAT 3 & 4 mounted on 2-meter telescopes within the Las Cumbres Observatory Global Telescope Network. Despite the expectation of variability, we observed no significant variation (<1.4%) on timescales of 6-10 hours. This non-detection is likely due to substantial host galaxy light diluting the subtle AGN variability. Dual-band preliminary tests and tailored simulations may enhance the possibility of detecting variability and lag in future IMBH reverberation campaigns. △ Less

Submitted 19 May, 2024; originally announced May 2024.

Comments: 14 pages, 6 figures, submitted to ApJ, comments welcome

Discovery of Very-high-energy Gamma-ray Emissions from the Low Luminosity AGN NGC 4278 by LHAASO

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) i… ▽ More The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) is compatible with NGC 4278 within $\sim0.03$ degree. Variation analysis shows an indication of the variability at a few months level in the TeV band, which is consistent with low frequency observations. Based on these observations, we report the detection of TeV $γ$-ray emissions from this low-luminosity AGN NGC 4278. The observations by LHAASO-WCDA during active period has a significance level of 8.8\,$σ$ with best-fit photon spectral index $\varGamma=2.56\pm0.14$ and a flux $f_{1-10\,\rm{TeV}}=(7.0\pm1.1_{\rm{sta}}\pm0.35_{\rm{syst}})\times10^{-13}\,\rm{photons\,cm^{-2}\,s^{-1}}$, or approximately $5\%$ of the Crab Nebula. The discovery of VHE from NGC 4278 indicates that the compact, weak radio jet can efficiently accelerate particles and emit TeV photons. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

SDSS J222428.53+261423.2: unique emission lines properties unveil a sub-kiloparsec dual AGN candidate

Authors: Qi Zheng, XueGuang Zhang, QiRong Yuan, Paola Severgnini, Cristian Vignali

Abstract: In this paper, we presented a detailed analysis of the Sloan Digital Sky Survey optical spectrum of a new sub-kpc scale dual AGN candidate SDSS J222428.53+261423.2 (=SDSS J2224). The target is one of the few AGNs with all the optical narrow emission lines characterized by double-peaked profiles and with peak separations in velocity units of about 930 km/s. If the double-peaked narrow emission line… ▽ More In this paper, we presented a detailed analysis of the Sloan Digital Sky Survey optical spectrum of a new sub-kpc scale dual AGN candidate SDSS J222428.53+261423.2 (=SDSS J2224). The target is one of the few AGNs with all the optical narrow emission lines characterized by double-peaked profiles and with peak separations in velocity units of about 930 km/s. If the double-peaked narrow emission lines (DPNELs) are due to a dual AGN in \obj, the estimated physical separation between the two cores is about 500 pc. Meanwhile, three alternative explanations are also discussed in this paper, however, we can not find solid evidence to completely rule them out. Our results support the presence of a sub-kpc dual AGN with DPNELs in all lines, indicating a key episode of galaxy merging evolution at sub-kpc scale. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: 9 pages, 5 figures, accepted to be published in MNRAS Letter

LHAASO-KM2A detector simulation using Geant4

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (254 additional authors not shown)

Abstract: KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with… ▽ More KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with large altitude difference (30 m) and huge coverage (1.3 km^2). In this paper, the design of the KM2A simulation code G4KM2A based on Geant4 is introduced. The process of G4KM2A is optimized mainly in memory consumption to avoid memory overffow. Some simpliffcations are used to signiffcantly speed up the execution of G4KM2A. The running time is reduced by at least 30 times compared to full detector simulation. The particle distributions and the core/angle resolution comparison between simulation and experimental data of the full KM2A array are also presented, which show good agreement. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Precise measurement of the cosmic-ray spectrum and $\left \langle \ln A \right \rangle$ by LHAASO -- connecting the Galactic to the extragalactic components

Authors: Xing-Jian Lv, Xiao-Jun Bi, Kun Fang, Yi-Qing Guo, Hui-Hai He, Ling-Ling Ma, Peng-Fei Yin, Qiang Yuan, Meng-Jie Zhao

Abstract: 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 LHAAS… ▽ More 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. △ Less

Submitted 18 March, 2024; originally announced March 2024.

Comments: 11 pages, 2 figures, 4 tables

Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (256 additional authors not shown)

Abstract: We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at… ▽ More We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at $3.67 \pm 0.05 \pm 0.15$ PeV. Below the knee, the spectral index is found to be -$2.7413 \pm 0.0004 \pm 0.0050$, while above the knee, it is -$3.128 \pm 0.005 \pm 0.027$, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -$0.1200 \pm 0.0003 \pm 0.0341$. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components. △ Less

Submitted 26 March, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comments: 8 pages, 3 figures

Journal ref: Physical Review Letters 132, 131002 (2024)

Injection spectra of different species of cosmic rays from AMS-02, ACE-CRIS and Voyager-1

Authors: Xu Pan, Qiang Yuan

Abstract: 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 w… ▽ More 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. △ Less

Submitted 11 March, 2024; originally announced March 2024.

Journal ref: Res. Astron. Astrophys. 23 115002 (2023)

Galactic diffuse emission from radio to ultra-high-energy gamma rays in light of up-to-date cosmic ray measurements

Authors: Xin-Yu He, Pei-Pei Zhang, Qiang Yuan, Yi-Qing Guo

Abstract: 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 prop… ▽ More 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. △ Less

Submitted 20 February, 2024; originally announced February 2024.

Comments: Accepted by ApJ

Journal ref: Xin-Yu He et al 2024 ApJ 964 28

Prospects for Joint Detection of Gravitational Waves with Counterpart Gamma-Ray Bursts Detected by the HADAR Experiment

Authors: Pei-Jin Hu, Qi-Ling Chen, Tian-Lu Chen, Ming-Ming Kang, Yi-Qing Guo, Dan-Zeng Luo-Bu, You-Liang Feng, Qi Gao, Quan-Bu Gou, Hong-Bo Hu, Hai-Jin Li, Cheng Liu, Mao-Yuan Liu, Wei Liu, Xiang-Li Qian, Bing-Qiang Qiao, Jing-Jing Su, Hui-Ying Sun, Xu Wang, Zhen Wang, Guang-Guang Xin, Chao-Wen Yang, Yu-Hua Yao, Qiang Yuan, Yi Zhang

Abstract: 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 v… ▽ More 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. △ Less

Submitted 20 January, 2024; originally announced January 2024.

Galaxy Clusters from the DESI Legacy Imaging Surveys.II. Environment effects on size-mass relation

Authors: Zhaobin Chen, Yizhou Gu, Hu Zou, Qirong Yuan

Abstract: To investigate the environment effects on size growth of galaxies, we study the size-mass relation across a broad range of environment with a vast sample of approximately 32 million galaxies at z < 0.5 from the DESI Legacy Imaging Surveys. This sample is divided into 3 subsamples representing galaxies within three different environments: brightest cluster galaxies (BCGs), other cluster galaxies, a… ▽ More To investigate the environment effects on size growth of galaxies, we study the size-mass relation across a broad range of environment with a vast sample of approximately 32 million galaxies at z < 0.5 from the DESI Legacy Imaging Surveys. This sample is divided into 3 subsamples representing galaxies within three different environments: brightest cluster galaxies (BCGs), other cluster galaxies, and field galaxies. The BCGs in our large sample are dominated by quiescent galaxies (QGs), while only a minority (~13%) of BCGs are star-forming galaxies (SFGs). To demonstrate the influence of environment on size growth, we attempt to observe the difference in size-mass relation for these three subsamples. In general, the slope differences between QGs and SFGs within various environments are significant, and tend to be greater at higher redshifts. For the mass-complete subsamples at z < 0.5, BCGs are found to have the highest slope of size-mass relation, and no difference in size-mass relation is found between cluster members and field galaxies. To assess whether the observed slope differences stem from the variations in environment or mass distribution, we construct the mass-matched subsamples for QGs and SFGs. As a result, both QGs and SFGs show negligible difference in slope of size-mass relation among the galaxies within three distinct environments, indicating that stellar mass is the most fundamental factor driving the size evolution at z < 0.5, though the mass growth mode for QGs and SFGs may have been affected by galaxy environment. △ Less

Submitted 28 December, 2023; originally announced December 2023.

Comments: 18 pages, 12 figures, Accepted by the ApJ

Does or did the supernova remnant Cassiopeia A operate as a PeVatron?

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE;… ▽ More For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE; $E_γ\geq 100$~TeV) $γ$-rays. In this context, the historical SNR Cassiopeia A (Cas A) is considered one of the most promising target for UHE observations. This paper presents the observation of Cas A and its vicinity by the LHAASO KM2A detector. The exceptional sensitivity of LHAASO KM2A in the UHE band, combined with the young age of Cas A, enabled us to derive stringent model-independent limits on the energy budget of UHE protons and nuclei accelerated by Cas A at any epoch after the explosion. The results challenge the prevailing paradigm that Cas A-type SNRs are major suppliers of PeV CRs in the Milky Way. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 11 pages, 3 figures, Accepted by the APJL

Inelastic Scattering of Dark Matter with Heavy Cosmic Rays

Authors: Keyu Lu, Yue-Lin Sming Tsai, Qiang Yuan, Le Zhang

Abstract: 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 inc… ▽ More 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_χ=0.1$ and a DM mass of $m_χ=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_χ\simeq 2$ MeV, the effective inelastic cross-section between DM and protons must be less than $\mathcal{O}(10^{-32})~{\rm cm}^2$. △ Less

Submitted 7 June, 2024; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 19 pages, 8 figures

Journal ref: Research in Astronomy and Astrophysics 24, 065007 (2024)

Very high energy gamma-ray emission beyond 10 TeV from GRB 221009A

Authors: Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the t… ▽ More The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the trigger. The intrinsic energy spectrum of gamma-rays can be described by a power-law after correcting for extragalactic background light (EBL) absorption. Such a hard spectrum challenges the synchrotron self-Compton (SSC) scenario of relativistic electrons for the afterglow emission above several TeV. Observations of gamma-rays up to 13 TeV from a source with a measured redshift of z=0.151 hints more transparency in intergalactic space than previously expected. Alternatively, one may invoke new physics such as Lorentz Invariance Violation (LIV) or an axion origin of very high energy (VHE) signals. △ Less

Submitted 22 November, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 49pages, 11figures

Journal ref: Science Advances, 9, eadj2778 (2023) 15 November 2023

Application of Deep Learning Methods Combined with Physical Background in Wide Field of View Imaging Atmospheric Cherenkov Telescopes

Authors: Ao-Yan Cheng, Hao Cai, Shi Chen, Tian-Lu Chen, Xiang Dong, You-Liang Feng, Qi Gao, Quan-Bu Gou, Yi-Qing Guo, Hong-Bo Hu, Ming-Ming Kang, Hai-Jin Li, Chen Liu, Mao-Yuan Liu, Wei Liu, Fang-Sheng Min, Chu-Cheng Pan, Bing-Qiang Qiao, Xiang-Li Qian, Hui-Ying Sun, Yu-Chang Sun, Ao-Bo Wang, Xu Wang, Zhen Wang, Guang-Guang Xin , et al. (3 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 11 October, 2023; v1 submitted 30 September, 2023; originally announced October 2023.

Listening for the Axion Echo with the 21 CentiMeter Array

Authors: Ariel Arza, Quan Guo, Lei Wu, Qiaoli Yang, Xiaolong Yang, Qiang Yuan, Bin Zhu

Abstract: 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… ▽ More 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μ\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. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Comments: 7 pages, 4 figures

A common origin of multi-messenger spectral anomaly of galactic cosmic rays

Authors: Yu-Hua Yao, Xu-Lin Dong, Yi-Qing Guo, Qiang Yuan

Abstract: 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 anomal… ▽ More 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γ$-ray emissions detected by the AS$\rmγ$ experiment. Furthermore, our results predict that the spectrum of diffuse $\rmγ$-ray is spatial-dependent, resting with local sources, which can be tested by LHAASO experiment. △ Less

Submitted 7 March, 2024; v1 submitted 30 August, 2023; originally announced August 2023.

Comments: 9 pages, 7 figures, accepted by PhysRevD

The Giant Radio Array for Neutrino Detection (GRAND) Collaboration -- Contributions to the 38th International Cosmic Ray Conference (ICRC 2023)

Authors: GRAND Collaboration, Rafael Alves Batista, Aurélien Benoit-Lévy, Teresa Bister, Mauricio Bustamante, Yiren Chen, LingMei Cheng, Simon Chiche, Jean-Marc Colley, Pablo Correa, Nicoleta Cucu Laurenciu, Zigao Dai, Beatriz de Errico, Sijbrand de Jong, João R. T. de Mello Neto, Krijn D. de Vries, Peter B. Denton, Valentin Deocoene, Kaikai Duan, Bohao Duan, Ralph Engel, Yizhong Fan, Arsène Ferrière, QuanBu Gou, Junhua Gu , et al. (74 additional authors not shown)

Abstract: 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 at… ▽ More 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. △ Less

Submitted 27 July, 2023; originally announced August 2023.

Comments: Note: To access the list of contributions, please follow the "HTML" link that can be found on the arXiv page

BGO quenching effect on spectral measurements of cosmic-ray nuclei in DAMPE experiment

Authors: Zhan-Fang Chen, Chuan Yue, Wei Jiang, Ming-Yang Cui, Qiang Yuan, Ying Wang, Cong Zhao, Yi-Feng Wei

Abstract: The Dark Matter Particle Explorer (DAMPE) is a satellite-borne detector designed to measure high energy cosmic-rays and $γ$-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, res… ▽ More The Dark Matter Particle Explorer (DAMPE) is a satellite-borne detector designed to measure high energy cosmic-rays and $γ$-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. △ Less

Submitted 24 July, 2023; originally announced July 2023.

Comments: 13 pages, 4 figures, to be published in Nuclear Inst. and Methods in Physics Research, A

Insights from LHAASO and IceCube into the origin of the Galactic diffuse TeV--PeV emission

Authors: Kai Yan, Ruo-Yu Liu, Rui Zhang, Chao-Ming Li, Qiang Yuan, Xiang-Yu Wang

Abstract: 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… ▽ More 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. △ Less

Submitted 15 March, 2024; v1 submitted 23 July, 2023; originally announced July 2023.

Comments: 28 pages, 13 figures; published in Nature Astronomy

Multi-messenger observations support cosmic ray interactions surrounding acceleration sources

Authors: Dong-Xu Sun, Pei-Pei Zhang, Yi-Qing Guo, Wei Liu, Qiang Yuan

Abstract: 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 t… ▽ More 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. △ Less

Submitted 5 July, 2023; originally announced July 2023.

Comments: 10 papes, 5 pages

Primordial magnetic field as a common solution of nanohertz gravitational waves and the Hubble tension

Authors: Yao-Yu Li, Chi Zhang, Ziwei Wang, Ming-Yang Cui, Yue-Lin Sming Tsai, Qiang Yuan, Yi-Zhong Fan

Abstract: 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… ▽ More 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{μ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}^{-α}$ 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 α\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. △ Less

Submitted 4 March, 2024; v1 submitted 29 June, 2023; originally announced June 2023.

Comments: 7 pages, 3 figures, 1 table

Journal ref: Phys. Rev. D 109, 043538 (2024)

Multi-wavelength temporal variability of the blazar PKS 1510-089

Authors: Q. Yuan, Pankaj Kushwaha, Alok C. Gupta, Ashutosh Tripathi, Paul J. Wiita, M. Zhang, X. Liu, Anne Lahteenmaki, Merja Tornikoski, Joni Tammi, Venkatessh Ramakrishnan, L. Cui, X. Wang, M. F. Gu, Cosimo Bambi, A. E. Volvach

Abstract: We perform correlation and periodicity search analyses on long-term multi-band light curves of the FSRQ 1510-089 observed by the space-based Fermi--Large Area Telescope in gamma-rays, the SMARTS and Steward Observatory telescopes in optical and near-infrared (NIR) and the 13.7 m radio telescope in Metsahovi Radio Observatory between 2008 and 2018. The z-transform discrete correlation function meth… ▽ More We perform correlation and periodicity search analyses on long-term multi-band light curves of the FSRQ 1510-089 observed by the space-based Fermi--Large Area Telescope in gamma-rays, the SMARTS and Steward Observatory telescopes in optical and near-infrared (NIR) and the 13.7 m radio telescope in Metsahovi Radio Observatory between 2008 and 2018. The z-transform discrete correlation function method is applied to study the correlation and possible time lags among these multi band light curves. Among all pairs of wavelengths, the gamma-ray vs. optical/NIR and optical vs. NIR correlations show zero time lags; however, both the gamma-ray and optical/NIR emissions precede the radio radiation. The Generalized Lomb-Scargle periodogram, Weighted Wavelet Z-transform, and REDFIT techniques are employed to investigate the unresolved-core-emission dominated 37 GHz light curve and yield evidence for a quasi-period around 1540 days, although given the length of the whole data set it cannot be claimed to be significant. We also investigate the optical/NIR color variability and find that this source shows a simple redder-when-brighter behavior over time, even in the low flux state. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ; 20 pages, 9 figures, 4 tables

The First LHAASO Catalog of Gamma-Ray Sources

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022.… ▽ More We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022. This catalog represents the main result from the most sensitive large coverage gamma-ray survey of the sky above 1 TeV, covering declination from $-$20$^{\circ}$ to 80$^{\circ}$. In total, the catalog contains 90 sources with an extended size smaller than $2^\circ$ and a significance of detection at $> 5σ$. Based on our source association criteria, 32 new TeV sources are proposed in this study. Among the 90 sources, 43 sources are detected with ultra-high energy ($E > 100$ TeV) emission at $> 4σ$ significance level. We provide the position, extension, and spectral characteristics of all the sources in this catalog. △ Less

Submitted 27 November, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 40 pages, 13 figures, 4 tables

Journal ref: The Astrophysical Journal Supplement Series, 271 (2024) 25

What boost galaxy mergers in two massive galaxy protoclusters at z = 2.24

Authors: Shuang Liu, Xianzhong Zheng, Dongdong Shi, Zheng Cai, Xiaohui Fan, Xin Wang, Qirong Yuan, Haiguang Xu, Zhizheng Pan, Wenhao Liu, Jianbo Qin, Yuheng Zhang, Run Wen

Abstract: Characterizing the structural properties of galaxies in high-redshift protoclusters is key to our understanding of the environmental effects on galaxy evolution in the early stages of galaxy and structure formation. In this study, we assess the structural properties of 85 and 87 Halpha emission-line candidates (HAEs) in the densest regions of two massive protoclusters, BOSS1244 and BOSS1542, respe… ▽ More Characterizing the structural properties of galaxies in high-redshift protoclusters is key to our understanding of the environmental effects on galaxy evolution in the early stages of galaxy and structure formation. In this study, we assess the structural properties of 85 and 87 Halpha emission-line candidates (HAEs) in the densest regions of two massive protoclusters, BOSS1244 and BOSS1542, respectively, using HST H-band imaging data. Our results show a true pair fraction of 22+-5 (33+-6) percent in BOSS1244 (BOSS1542), which yields a merger rate of 0.41+-0.09 (0.52+-0.04) per Gyr for massive HAEs with log (M_*/M_sun) > 10.3. This rate is 1.8 (2.8) times higher than that of the general fields at the same epoch. Our sample of HAEs exhibits half-light radii and Sersic indices that cover a broader range than field star-forming galaxies. Additionally, about 15 percent of the HAEs are as compact as the most massive (log(M_*/M_sun) > 11) spheroid-dominated population. These results suggest that the high galaxy density and cold dynamical state (i.e., velocity dispersion of <400 km/s) are key factors that drive galaxy mergers and promote structural evolution in the two protoclusters. Our findings also indicate that both the local environment (on group scales) and the global environment play essential roles in shaping galaxy morphologies in protoclusters. This is evident in the systematic differences observed in the structural properties of galaxies between BOSS1244 and BOSS1542. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 20 pages, 5 figures, 1 table. Accepted for publication in MNRAS

Galactic diffuse gamma-ray emission from GeV to PeV energies in light of up-to-date cosmic ray measurements

Authors: Rui Zhang, Xiaoyuan Huang, Zhi-Hui Xu, Shiping Zhao, Qiang Yuan

Abstract: 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… ▽ More 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. △ Less

Submitted 8 September, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Journal ref: The Astrophysical Journal, 2023, 957(1): 43

Measurement of ultra-high-energy diffuse gamma-ray emission of the Galactic plane from 10 TeV to 1 PeV with LHAASO-KM2A

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer ar… ▽ More The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer array of the Large High Altitude Air Shower Observatory (LHAASO). Diffuse emissions from the inner ($15^{\circ}<l<125^{\circ}$, $|b|<5^{\circ}$) and outer ($125^{\circ}<l<235^{\circ}$, $|b|<5^{\circ}$) Galactic plane are detected with $29.1σ$ and $12.7σ$ significance, respectively. The outer Galactic plane diffuse emission is detected for the first time in the very- to ultra-high-energy domain ($E>10$~TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of $-2.99\pm0.04$, which is different from the curved spectrum as expected from hadronic interactions between locally measured cosmic rays and the line-of-sight integrated gas content. Furthermore, the measured flux is higher by a factor of $\sim3$ than the prediction. A similar spectrum with an index of $-2.99\pm0.07$ is found in the outer Galaxy region, and the absolute flux for $10\lesssim E\lesssim60$ TeV is again higher than the prediction for hadronic cosmic ray interactions. The latitude distributions of the diffuse emission are consistent with the gas distribution, while the longitude distributions show clear deviation from the gas distribution. The LHAASO measurements imply that either additional emission sources exist or cosmic ray intensities have spatial variations. △ Less

Submitted 19 August, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: 12 pages, 8 figures, 5 tables; accepted for publication in Physical Review Letters; source mask file provided as ancillary file

Journal ref: Phys. Rev. Lett. 131, 151001 (2023)

Potential of Constraining Propagation Parameters of Galactic Cosmic Rays with the High Energy cosmic-Radiation Detection facility onboard China's Space Station

Authors: Zhi-Hui Xu, Qiang Yuan, Zhi-Cheng Tang, Xiao-Jun Bi

Abstract: 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 investigat… ▽ More 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. △ Less

Submitted 21 April, 2023; originally announced April 2023.

Comments: 10 pages, 4 figures

Constraining ultralight dark matter using the Fermi-LAT pulsar timing array

Authors: Zi-Qing Xia, Tian-Peng Tang, Xiaoyuan Huang, Qiang Yuan, Yi-Zhong Fan

Abstract: 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… ▽ More 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. △ Less

Submitted 29 June, 2023; v1 submitted 30 March, 2023; originally announced March 2023.

Comments: 11 pages, 4 figures, 2 tables; publish as a Letter in Phys. Rev. D

Journal ref: Phys. Rev. D 107, L121302 (Published 28 June 2023)

Prospects for detection rate of very-high-energy γ-ray emissions from short γ-ray bursts with the HADAR experiment

Authors: Qi-Ling Chen, Pei-Jin Hu, Jing-Jing Su, Ming-Ming Kang, Yi-Qing Guo, Tian-Lu Chen, Dan-Zeng Luo-Bu, Yu-fan Fan, You-Liang Feng, Qi Gao, Quan-Bu Gou, Hong-Bo Hu, Hai-Jin Li, Cheng Liu, Mao-Yuan Liu, Wei Liu, Xiang-Li Qian, Bing-Qiang Qiao, Hui-Ying Sun, Xu Wang, Zhen Wang, Guang-Guang Xin, Yu-Hua Yao, Qiang Yuan, Yi Zhang

Abstract: 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… ▽ More 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) $γ$-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 $γ$-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. △ Less

Submitted 3 April, 2023; v1 submitted 27 March, 2023; originally announced March 2023.

Journal ref: Chin. Phys. C 47, 095001 (2023)

Exploring dark matter spike distribution around the Galactic centre with stellar orbits

Authors: Zhao-Qiang Shen, Guan-Wen Yuan, Cheng-Zi Jiang, Yue-Lin Sming Tsai, Qiang Yuan, Yi-Zhong Fan

Abstract: 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 liter… ▽ More 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 $γ\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 $α=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< σv \right> \gtrsim 7.7\times 10^{-27}~{\rm cm^3\,s^{-1}} (m_{\rm DM}/100~{\rm GeV})$ at 95% level. △ Less

Submitted 24 October, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 12 pages, 10 figures and 1 table. Accepted for publication in MNRAS

Journal ref: MNRAS 527, 3196 (2024)

Correlation analysis between several observables of OJ 287 radio jet

Authors: Qi Yuan, Ming Zhang, Xiang Liu, Pengfei Jiang, Gulchehra I. Kokhirova

Abstract: We collected the archival data of blazar OJ~287 from heterogeneous very long baseline interferometry (VLBI) monitoring programs at 2.3 GHz, 8.6 GHz, 15 GHz and 43 GHz. The data reduction and observable extraction of those multi-band multi-epoch observations are batch-processed consistently with our automated pipeline. We present the multivariate correlation analysis on the observables at each band… ▽ More We collected the archival data of blazar OJ~287 from heterogeneous very long baseline interferometry (VLBI) monitoring programs at 2.3 GHz, 8.6 GHz, 15 GHz and 43 GHz. The data reduction and observable extraction of those multi-band multi-epoch observations are batch-processed consistently with our automated pipeline. We present the multivariate correlation analysis on the observables at each band. We employ the cross-correlation function to search the correlations and the Monte Carlo (MC) technique to verify the certainty of correlations. Several correlations are found. The foremost findings are the correlations between the core flux density and the jet position angles on different scales, which validated the plausible predictions of the jet with precession characteristics. Meanwhile, there is a variation in the offset between the core EVPA and the inner-jet position angle over time at 15~GHz and 43~GHz. △ Less

Submitted 3 April, 2023; v1 submitted 16 March, 2023; originally announced March 2023.

Comments: 18 pages, 10 figures, 3 tables. Accepted for publication in The Astrophysical Journal

Velocity offset between emission and absorption lines might be an effective indicator of dual core system

Authors: Qi Zheng, Shuang Liu, Xueguang Zhang, Qirong Yuan

Abstract: This paper presents a detection of significant velocity offset between emission and absorption lines for a dual core system in SDSS~J155708.82+273518.74 (= SDSS~J1557). The photometric image of SDSS~J1557 exhibits clear two cores with a projected separation of $\sim$2.2 arcseconds (4.9 kpc) determined by GALFIT. Based on the applications of the commonly accepted pPXF code with 636 theoretical SSP… ▽ More This paper presents a detection of significant velocity offset between emission and absorption lines for a dual core system in SDSS~J155708.82+273518.74 (= SDSS~J1557). The photometric image of SDSS~J1557 exhibits clear two cores with a projected separation of $\sim$2.2 arcseconds (4.9 kpc) determined by GALFIT. Based on the applications of the commonly accepted pPXF code with 636 theoretical SSP templates, the host galaxy contribution can be well determined. Then, the emission line features of SDSS~J1557 can be well measured after subtraction of host starlight. It is found that the velocity offset of emission lines with respect to absorption lines reaches $458 \pm 13$ km/s. According to the Baldwin-Phillips-Terlevich (BPT) diagram, SDSS J1557 is a composite galaxy. In addition, SDSS J1557 can well fit the $M_{\rm BH}-σ_{\ast}$ relation of bulges and the galaxy merger would not change this relation. Two reasonable models (say, AGN-driven outflow vs. dual core system) have been discussed to explain this velocity offset. The model of AGN-driven outflow fails to interpret the systematic redshift of emission lines and similar velocity offsets for various emission lines in SDSS~J1557. A significant velocity offset between emission and absorption lines might be an effective indicator of dual core system. △ Less

Submitted 12 December, 2022; originally announced December 2022.

Comments: 13 pages, 7 figures, accepted to be published in ApJ

An Unsupervised Machine Learning Method for Electron--Proton Discrimination of the DAMPE Experiment

Authors: Zhihui Xu, Xiang Li, Mingyang Cui, Chuan Yue, Wei Jiang, Wenhao Li, Qiang Yuan

Abstract: 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 Pa… ▽ More 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. △ Less

Submitted 4 December, 2022; originally announced December 2022.

Comments: 10 pages, 5 figures, 1 table

Journal ref: Universe 2022, 8(11), 570

A delayed 400 GeV photon from GRB 221009A and implication on the intergalactic magnetic field

Authors: Zi-Qing Xia, Yun Wang, Qiang Yuan, Yi-Zhong Fan

Abstract: 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 t… ▽ More 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. △ Less

Submitted 9 May, 2024; v1 submitted 24 October, 2022; originally announced October 2022.

Comments: 9 pages, 4 figures, Nature Communications, in press

Journal ref: Nature Communications 15, 4280 (2024)

Interpretations of the cosmic ray secondary-to-primary ratios measured by DAMPE

Authors: Peng-Xiong Ma, Zhi-Hui Xu, Qiang Yuan, Xiao-Jun Bi, Yi-Zhong Fan, Igor V. Moskalenko, Chuan Yue

Abstract: 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… ▽ More 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. △ Less

Submitted 6 March, 2023; v1 submitted 17 October, 2022; originally announced October 2022.

Comments: 14 pages, 7 figures. Published version

Journal ref: Front. Phys. 18, 44301 (2023)

Prospective Study on Observations of γ-Ray Sources in the Galaxy Using the HADAR Experiment

Authors: Xiangli Qian, Huiying Sun, Tianlu Chen, Danzengluobu, Youliang Feng, Qi Gao, Quanbu Gou, Yiqing Guo, Hongbo Hu, Mingming Kang, Haijin Li, Cheng Liu, Maoyuan Liu, Wei Liu, Bingqiang Qiao, Xu Wang, Zhen Wang, Guangguang Xin, Yuhua Yao, Qiang Yuan, Yi Zhang

Abstract: 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 thr… ▽ More 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 γ-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 (σ). The statistical significance for the Crab Nebula during one year of operation reached 346.0 σ and the one-year integral sensitivity of HADAR above 1TeV was ~1.3%-2.4% of the flux from the Crab Nebula. △ Less

Submitted 23 September, 2022; originally announced September 2022.

Searching for axion dark matter with MeerKAT Radio Telescope

Authors: Yun-Fan Zhou, Nick Houston, Gyula I. G. Jozsa, Hao Chen, Yin-Zhe Ma, Qiang Yuan, Tao An, Yogesh Chandola, Ran Ding, Fujun Du, Shao-Guang Guo, Xiaoyuan Huang, Mengtian Li, Chandreyee Sengupta

Abstract: 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… ▽ More 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σ$ in the data, we provide a unique constraint on the Primakoff coupling of axion DM, $g_{{\rm a}γγ}\lesssim 9.3 \times 10^{-12}\,{\rm GeV}^{-1}$ at the $95\%$ confidence level, in the mass range $3.18$-$4.35\,μ$eV. This result serves the strongest constraint in the axion mass range $4.20$-$4.35\,μ$eV. △ Less

Submitted 27 October, 2023; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 6 pages, 3 figures, accepted by Physical Review D

Journal ref: Physical Review D 106, 083006 (2022)

Search for relativistic fractionally charged particles in space

Authors: DAMPE Collaboration, F. Alemanno, C. Altomare, Q. An, P. Azzarello, F. C. T. Barbato, P. Bernardini, X. J. Bi, M. S. Cai, E. Casilli, E. Catanzani, J. Chang, D. Y. Chen, J. L. Chen, Z. F. Chen, M. Y. Cui, T. S. Cui, Y. X. Cui, H. T. Dai, A. De-Benedittis, I. De Mitri, F. de Palma, M. Deliyergiyev, A. Di Giovanni, M. Di Santo , et al. (126 additional authors not shown)

Abstract: More than a century after the performance of the oil drop experiment, the possible existence of fractionally charged particles FCP still remains unsettled. The search for FCPs is crucial for some extensions of the Standard Model in particle physics. Most of the previously conducted searches for FCPs in cosmic rays were based on experiments underground or at high altitudes. However, there have been… ▽ More More than a century after the performance of the oil drop experiment, the possible existence of fractionally charged particles FCP still remains unsettled. The search for FCPs is crucial for some extensions of the Standard Model in particle physics. Most of the previously conducted searches for FCPs in cosmic rays were based on experiments underground or at high altitudes. However, there have been few searches for FCPs in cosmic rays carried out in orbit other than AMS-01 flown by a space shuttle and BESS by a balloon at the top of the atmosphere. In this study, we conduct an FCP search in space based on on-orbit data obtained using the DArk Matter Particle Explorer (DAMPE) satellite over a period of five years. Unlike underground experiments, which require an FCP energy of the order of hundreds of GeV, our FCP search starts at only a few GeV. An upper limit of $6.2\times 10^{-10}~~\mathrm{cm^{-2}sr^{-1} s^{-1}}$ is obtained for the flux. Our results demonstrate that DAMPE exhibits higher sensitivity than experiments of similar types by three orders of magnitude that more stringently restricts the conditions for the existence of FCP in primary cosmic rays. △ Less

Submitted 9 September, 2022; originally announced September 2022.

Comments: 19 pages, 6 figures, accepted by PRD

Report number: 106, 063026

Journal ref: Physical Review D 106.6 (2022): 063026

Birefringence Tomography for Axion Cloud

Authors: Yifan Chen, Chunlong Li, Yosuke Mizuno, Jing Shu, Xiao Xue, Qiang Yuan, Yue Zhao, Zihan Zhou

Abstract: 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 rec… ▽ More 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. △ Less

Submitted 28 September, 2022; v1 submitted 11 August, 2022; originally announced August 2022.

Comments: 36 pages, 13 figures, published version in JCAP

Journal ref: JCAP09(2022)073

Flux Variations of Cosmic Ray Air Showers Detected by LHAASO-KM2A During a Thunderstorm on 10 June 2021

Authors: LHAASO Collaboration, F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Zhe Cao, Zhen Cao, J. Chang, J. F. Chang, E. S. Chen, Liang Chen, Liang Chen, Long Chen, M. J. Chen, M. L. Chen, S. H. Chen, S. Z. Chen, T. L. Chen, X. J. Chen , et al. (248 additional authors not shown)

Abstract: The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations… ▽ More The Large High Altitude Air Shower Observatory (LHAASO) has three sub-arrays, KM2A, WCDA and WFCTA. The flux variations of cosmic ray air showers were studied by analyzing the KM2A data during the thunderstorm on 10 June 2021. The number of shower events that meet the trigger conditions increases significantly in atmospheric electric fields, with maximum fractional increase of 20%. The variations of trigger rates (increases or decreases) are found to be strongly dependent on the primary zenith angle. The flux of secondary particles increases significantly, following a similar trend with that of the shower events. To better understand the observed behavior, Monte Carlo simulations are performed with CORSIKA and G4KM2A (a code based on GEANT4). We find that the experimental data (in saturated negative fields) are in good agreement with simulations, assuming the presence of a uniform upward electric field of 700 V/cm with a thickness of 1500 m in the atmosphere above the observation level. Due to the acceleration/deceleration and deflection by the atmospheric electric field, the number of secondary particles with energy above the detector threshold is modified, resulting in the changes in shower detection rate. △ Less

Submitted 6 December, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 18 pages, 11 figures

Journal ref: Chinese Phys. C 47 015001 (2023)

About 300 days optical quasi-periodic oscillations in the long-term light curves of the blazar PKS 2155-304

Authors: Qi Zheng, Xueguang Zhang, Qirong Yuan

Abstract: Based on the long-term light curves collected from the Catalina Sky Survey (CSS) (from 2005 to 2013) and the All-Sky Automated Survey for Supernovae (ASAS-SN) (from 2014 to 2018), optical quasi-periodic oscillations (QPOs) about 300 days can be well determined in the well-known blazar PKS~2155-304 through four different methods: the generalized Lomb-Scargle periodogram (GLSP) method, the weighted… ▽ More Based on the long-term light curves collected from the Catalina Sky Survey (CSS) (from 2005 to 2013) and the All-Sky Automated Survey for Supernovae (ASAS-SN) (from 2014 to 2018), optical quasi-periodic oscillations (QPOs) about 300 days can be well determined in the well-known blazar PKS~2155-304 through four different methods: the generalized Lomb-Scargle periodogram (GLSP) method, the weighted wavelet Z-transform (WWZ) technique, the epoch-folded method and redfit method. The GLSP determined significance level for the periodicity is higher than 99.9999\% based on a false alarm probability. The redfit provided confidence level for the periodicity is higher than 99\% in ASAS-SN light curve, after considering the effects of red noise. Based on continuous autoregressive (CAR) process created artificial light curves, the probability of detecting fake QPOs is lower than 0.8\%. The determined optical periodicity of 300 days from CSS and ASAS-SN light curves is well consistent with the reported optical periodicity in the literature. Moreover, three possible models are discussed to explain the optical QPOs in PKS 2155-304: the relativistic Frame-dragging effect, the binary black hole (BBH) model and the jet precession model. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 13 pages, 9 figures, accepted to be published in RAA

Searching for cosmic string induced stochastic gravitational wave background with the Parkes Pulsar Timing Array

Authors: Ligong Bian, Jing Shu, Bo Wang, Qiang Yuan, Junchao Zong

Abstract: 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 p… ▽ More 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μ<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. △ Less

Submitted 15 May, 2022; originally announced May 2022.

Comments: 10 pages, 8 figures, 4 tables. Comments welcome

Geminga contribution to the cosmic-ray positron excess according to the gamma-ray observations

Authors: Guang-Yao Zhou, Zhao-Huan Yu, Qiang Yuan, Hong-Hao Zhang

Abstract: We attempt to interpret the cosmic-ray positron excess by injection from the nearby pulsar Geminga, assuming a two-zone diffusion scenario and an injection spectrum with a low energy cutoff. Since the high energy positrons and electrons from Geminga can induce $γ$ rays via inverse Compton scattering, we take into account the extended $γ$-ray observations around Geminga from HAWC for $\sim 10$ TeV… ▽ More We attempt to interpret the cosmic-ray positron excess by injection from the nearby pulsar Geminga, assuming a two-zone diffusion scenario and an injection spectrum with a low energy cutoff. Since the high energy positrons and electrons from Geminga can induce $γ$ rays via inverse Compton scattering, we take into account the extended $γ$-ray observations around Geminga from HAWC for $\sim 10$ TeV and from Fermi-LAT for $\mathcal{O}(10)$ GeV. According to the extended $γ$-ray observation claimed by an analysis of Fermi-LAT data, we find that Geminga could explain the positron excess for a $30\%$ energy conversion efficiency into positrons and electrons. However, based on the constraint on the extended $γ$ rays given by another Fermi-LAT analysis, positrons from Geminga would be insufficient to account for the positron excess. A further robust analysis of Fermi-LAT data for the extended $γ$ rays would be crucial to determine whether Geminga can explain the positron excess or not. △ Less

Submitted 26 September, 2022; v1 submitted 14 May, 2022; originally announced May 2022.

Comments: 14 pages, 2 figures; revisions to match the published version

Journal ref: Commun. Theor. Phys. 74 (2022) 105403

Constraining ultralight bosonic dark matter with Keck observations of S2's orbit and kinematics

Authors: Guan-Wen Yuan, Zhao-Qiang Shen, Yue-Lin Sming Tsai, Qiang Yuan, Yi-Zhong Fan

Abstract: Ultralight bosonic dark matter is expected to be able to form a cloud surrounding the supermassive black hole (SMBH) in the Galactic center. With increasing precision of the observations of the stellar kinematics around the SMBH, tiny effects from such a dark matter cloud, including its gravitational perturbation and the direct coupling with the ordinary matter may be detectable. In this work, we… ▽ More Ultralight bosonic dark matter is expected to be able to form a cloud surrounding the supermassive black hole (SMBH) in the Galactic center. With increasing precision of the observations of the stellar kinematics around the SMBH, tiny effects from such a dark matter cloud, including its gravitational perturbation and the direct coupling with the ordinary matter may be detectable. In this work, we search for possible evidence of the scalar cloud using accurate orbital measurements of the S2 star around Sgr~A*. We solve the first order Post-Newtonian equation, considering simultaneously the extended mass distribution of the scalar cloud and the frequency shift induced by the additional coupling via Higgs portal or photon portal interaction. Furthermore, we also investigate the impact of an astrophysical power-law component from the gas and stellar remnants. We find that the astrometric and spectroscopic data of the S2 star are well consistent with the scenario of a point-like mass of Sgr~A*. We thus derive upper limits of the coupling of the new interaction and the extended mass, with and without the contribution from the astrophysical component. The limits of the Higgs/photon coupling and the extended mass of the scalar cloud are the most stringent ones for the scalar mass window between $3.2\times 10^{-19}$~eV and $1.6\times 10^{-18}$~eV. △ Less

Submitted 21 November, 2022; v1 submitted 10 May, 2022; originally announced May 2022.

Comments: 12 pages, 6 figures, published at PRD

Journal ref: Phys.Rev.D 106 (2022) 10, 103024

Explaining the GeV antiproton/$γ-$ray excesses and W-boson mass anomaly in an inert two Higgs doublet model

Authors: Cheng-Rui Zhu, Ming-Yang Cui, Zi-Qing Xia, Zhao-Huan Yu, Xiaoyuan Huang, Qiang Yuan, Yi-Zhong Fan

Abstract: For the newly discovered $W$-boson mass anomaly, one of the simplest dark matter (DM) models that can account for the anomaly without violating other astrophysical/experimental constraints is the inert two Higgs doublet model, in which the DM mass ($m_{S}$) is found to be within $\sim 54-74$ GeV. In this model, the annihilation of DM via $SS\to b\bar{b}$ and $SS\to WW^{*}$ would produce antiproton… ▽ More For the newly discovered $W$-boson mass anomaly, one of the simplest dark matter (DM) models that can account for the anomaly without violating other astrophysical/experimental constraints is the inert two Higgs doublet model, in which the DM mass ($m_{S}$) is found to be within $\sim 54-74$ GeV. In this model, the annihilation of DM via $SS\to b\bar{b}$ and $SS\to WW^{*}$ would produce antiprotons and gamma rays, and may account for the excesses identified previously in both particles. Motivated by this, we re-analyze the AMS-02 antiproton and Fermi-LAT Galactic center gamma-ray data. For the antiproton analysis, the novel treatment is the inclusion of the charge-sign-dependent three-dimensional solar modulation model as constrained by the time-dependent proton data. We find that the excess of antiprotons is more distinct than previous results based on the force-field solar modulation model. The interpretation of this excess as the annihilation of $SS\to WW^{*}$ ($SS\to b\bar{b}$) requires a DM mass of $\sim 40-80$ ($40-60$) GeV and a velocity-averaged cross section of $O(10^{-26})~{\rm cm^3~s^{-1}}$. As for the $γ$-ray data analysis, besides adopting the widely-used spatial template fitting, we employ an orthogonal approach with a data-driven spectral template analysis. The fitting to the GeV $γ$-ray excess yields DM model parameters overlapped with those to fit the antiproton excess via the $WW^{*}$ channel. The consistency of the DM particle properties required to account for the $W$-boson mass anomaly, the GeV antiproton excess, and the GeV $γ$-ray excess suggest a common origin of them. △ Less

Submitted 6 November, 2022; v1 submitted 7 April, 2022; originally announced April 2022.

Comments: 10 pages, 10 figures, accepted by PRL

Journal ref: Phys.Rev.Lett. 129 (2022) 23, 231101

Radio pulsations from a neutron star within the gamma-ray binary LS I +61$^{\circ}$ 303

Authors: Shan-Shan Weng, Lei Qian, Bo-Jun Wang, D. F. Torres, A. Papitto, Peng Jiang, Renxin Xu, Jian Li, Jing-Zhi Yan, Qing-Zhong Liu, Ming-Yu Ge, Qi-Rong Yuan

Abstract: LS I +61$^{\circ}$ 303 is one of the rare gamma-ray binaries, emitting most of their luminosity in photons with energies beyond 100 MeV. The $\sim$26.5 d orbital period is clearly detected at many wavelengths. Additional aspects of its multi-frequency behavior make it the most interesting example of the class. The morphology of high-resolution radio images changes with orbital phase displaying a c… ▽ More LS I +61$^{\circ}$ 303 is one of the rare gamma-ray binaries, emitting most of their luminosity in photons with energies beyond 100 MeV. The $\sim$26.5 d orbital period is clearly detected at many wavelengths. Additional aspects of its multi-frequency behavior make it the most interesting example of the class. The morphology of high-resolution radio images changes with orbital phase displaying a cometary tail pointing away from the high-mass star. LS I +61$^{\circ}$ 303 also shows superorbital variability. A couple of energetic ($\sim 10^{37}$ erg s$^{-1}$), short, magnetar-like bursts have been plausibly ascribed to it. LS I +61$^{\circ}$ 303's phenomenology has been put under theoretical scrutiny for decades, but the lack of certainty regarding the nature of the compact object in the binary has prevented advancing our understanding of the source. Here, using observations done with the Five-hundred-meter Aperture Spherical radio Telescope (FAST), we report on the existence of transient radio pulsations from the direction of LS I +61$^{\circ}$ 303. We find a period $P=269.15508 \pm 0.00016$ ms at a significance of $> 20σ$. This is the first evidence for pulsations from this source at any frequency, and strongly argues for the existence of a rotating neutron star in LS I +61$^{\circ}$ 303. △ Less

Submitted 17 March, 2022; originally announced March 2022.

Comments: Published in Nature Astronomy; authors' version

Role of the companion lensing galaxy in the CLASS gravitational lens B1152+199

Authors: M. Zhang, Q. Yuan, J. -Y. Liu, L. Zhang

Abstract: We reinvestigate the Cosmic Lens All-Sky Survey (CLASS) gravitational lens B1152+199 using archived Hubble Space Telescope (HST) data and Very Long Baseline Interferometry (VLBI) data. A consistent luminosity ratio within effective radius between the host galaxy and the X-galaxy is measured from HST tri-band images, which leads to a mass ratio between the two galaxies as $r_b\sim 2$. To determin… ▽ More We reinvestigate the Cosmic Lens All-Sky Survey (CLASS) gravitational lens B1152+199 using archived Hubble Space Telescope (HST) data and Very Long Baseline Interferometry (VLBI) data. A consistent luminosity ratio within effective radius between the host galaxy and the X-galaxy is measured from HST tri-band images, which leads to a mass ratio between the two galaxies as $r_b\sim 2$. To determine the role of the X-galaxy in the lens system, we modelled the dual-lens system with constraints from the VLBI-resolved jet components and the HST images. The 8.4-GHz global-VLBI data currently provide the most stringent constraints to the mass model, especially to the radial power-law slope. The optimized models for this two-image three-component radio lens favour a steeper-than-isothermal inner slope. The jet bending in image B was also investigated and it turned out to be rather a misalignment than a curvature. The goodness of fit indicates that the role of the X-galaxy is crucial in the lens system if three pairs of resolved jet components are to be fitted. When we imported the optimal model from radio constraints to optical modelling with the HST tri-band data, the optimization kept the consistency of the optimal model and successfully reproduced the features observed in the HST images. This implies that the diffuse emission discovered in the HST images is actually a detection of the secondary lensing effects from the companion lens. △ Less

Submitted 14 February, 2022; originally announced February 2022.

Comments: Published in MNRAS; 20 pages, 14 figures