Showing 1–50 of 58 results for author: Gu, H

The UV/optical Continuum Reverberation Mapping of Eight Active Galactic Nuclei with Swift: Further Evidence for the Outer Component

Authors: Chi-Zhuo Wang, Xue-Bing Wu, Yuanzhe Jiang, Qinchun Ma, Huapeng Gu, Yuhan Wen

Abstract: In our previous work, we applied the ICCF-Cut method to the continuum reverberation mapping (CRM) of six active galactic nuclei (AGNs) based on the published Swift data. Extending this work, we perform a systematic AGN CRM study utilizing the Swift archive. We enlarge our sample with eight additional AGNs at $z<0.05$ with high-cadence ($<3$ days) and multiband photometric observations. Time series… ▽ More In our previous work, we applied the ICCF-Cut method to the continuum reverberation mapping (CRM) of six active galactic nuclei (AGNs) based on the published Swift data. Extending this work, we perform a systematic AGN CRM study utilizing the Swift archive. We enlarge our sample with eight additional AGNs at $z<0.05$ with high-cadence ($<3$ days) and multiband photometric observations. Time series analysis of these light curves shows two main results: (1) The interband lags are broadly consistent with $τ\propto λ^{4/3}$, while the average interband lags are larger than those predicted by the standard thin accretion disk model. (2) For most targets, there exists a $U$ band lag excess, which is probably due to the diffuse continuum emission from the broad-line region (BLR). We employ the ICCF-Cut method to extract the possible diffuse continuum component from the $U$ band light curves and calculate the diffuse continuum lags ($τ_{\rm cut}$), which are generally consistent with the lags ($τ_{\rm jav}$) derived by the JAVELIN Photometric Reverberation Mapping Model. Further analysis with our sample indicates a positive correlation between the diffuse continuum region size and the BLR size ($R_{\rm DCR}-R_{\rm BLR}$ relation), as well as another correlation with the luminosity ($R_{\rm DCR}-L$ relation). These findings provide further evidence for a significant contribution of diffuse continuum emission from the BLR to the AGN continuum lags. △ Less

Submitted 7 November, 2025; originally announced November 2025.

Comments: Accepted for publication in ApJ, 26 pages, 16 figures, 4 tables. The processed light curves are available on Zenodo at DOI: 10.5281/zenodo.14930797

Evidence of cosmic-ray acceleration up to sub-PeV energies in the supernova remnant IC 443

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

Abstract: Supernova remnants (SNRs) have been considered as the primary contributors to cosmic rays (CRs) in our Galaxy. However, the maximum energy of particles that can be accelerated by shocks of SNRs is uncertain observationally and theoretically, and the role of contribution to CRs around PeV energies by SNRs is unclear. In this study, we present observations of high-energy $γ$-ray emission from the SN… ▽ More Supernova remnants (SNRs) have been considered as the primary contributors to cosmic rays (CRs) in our Galaxy. However, the maximum energy of particles that can be accelerated by shocks of SNRs is uncertain observationally and theoretically, and the role of contribution to CRs around PeV energies by SNRs is unclear. In this study, we present observations of high-energy $γ$-ray emission from the SNR IC 443 using the Large High Altitude Air Shower Observatory (LHAASO). The morphological analysis reveals a pointlike source whose location and spectrum are consistent with those of the Fermi-LAT-detected compact source with $π^0$-decay signature, and a more extended source which is consistent with a newly discovered source, previously unrecognized by Fermi-LAT. The spectrum of the point source can be described by a power-law function with an index of $\sim3.0$, extending beyond $\sim 30$ TeV without apparent cutoff. Assuming a hadronic origin of the $γ$-ray emission, the $95\%$ lower limit of accelerated protons reaches about 300 TeV. The extended source might be coincident with IC 443, SNR G189.6+3.3 or the putative pulsar wind nebula CXOU J061705.3+222127, and can be explained by either a hadronic or leptonic model. The LHAASO results provide compelling evidence that CR protons up to sub-PeV energies can be accelerated by the SNR. △ Less

Submitted 29 October, 2025; originally announced October 2025.

Slitless Spectroscopy Source Detection Using YOLO Deep Neural Network

Authors: Xiaohan Chen, Man I Lam, Yingying Zhou, Hongrui Gu, Jinzhi Lai, Zhou Fan, Jing Li, Xin Zhang, Hao Tian

Abstract: Slitless spectroscopy eliminates the need for slits, allowing light to pass directly through a prism or grism to generate a spectral dispersion image that encompasses all celestial objects within a specified area. This technique enables highly efficient spectral acquisition. However, when processing CSST slitless spectroscopy data, the unique design of its focal plane introduces a challenge: photo… ▽ More Slitless spectroscopy eliminates the need for slits, allowing light to pass directly through a prism or grism to generate a spectral dispersion image that encompasses all celestial objects within a specified area. This technique enables highly efficient spectral acquisition. However, when processing CSST slitless spectroscopy data, the unique design of its focal plane introduces a challenge: photometric and slitless spectroscopic images do not have a one-to-one correspondence. As a result, it becomes essential to first identify and count the sources in the slitless spectroscopic images before extracting spectra. To address this challenge, we employed the You Only Look Once (YOLO) object detection algorithm to develop a model for detecting targets in slitless spectroscopy images. This model was trained on 1,560 simulated CSST slitless spectroscopic images. These simulations were generated from the CSST Cycle 6 and Cycle 9 main survey data products, representing the Galactic and nearby galaxy regions and the high galactic latitude regions, respectively. On the validation set, the model achieved a precision of 88.6% and recall of 90.4% for spectral lines, and 87.0% and 80.8% for zeroth-order images. In testing, it maintained a detection rate >80% for targets brighter than 21 mag (medium-density regions) and 20 mag (low-density regions) in the Galactic and nearby galaxies regions, and >70% for targets brighter than 18 mag in high galactic latitude regions. △ Less

Submitted 12 October, 2025; originally announced October 2025.

Comments: 22 pages, 11 figures, 7 tables. PASP accepted

Photometric Analysis of 30 Contact Binaries in M31

Authors: Xiang Gao, Kai Li, Li-Heng Wang, Hai-bo Yuan, Hong-rui Gu, Ya-Ni Guo

Abstract: M31, as the largest galaxy in the Local Group, is of significant importance for the study of stellar formation and evolution. Based on the data of 5,859 targets observed in M31 by Gu et al (2024), we selected 30 contact binaries by visual inspection for further study. Using the PHOEBE software and employing Bayesian optimization and Markov Chain Monte Carlo sampling, we determined the physical par… ▽ More M31, as the largest galaxy in the Local Group, is of significant importance for the study of stellar formation and evolution. Based on the data of 5,859 targets observed in M31 by Gu et al (2024), we selected 30 contact binaries by visual inspection for further study. Using the PHOEBE software and employing Bayesian optimization and Markov Chain Monte Carlo sampling, we determined the physical parameters of these 30 systems. The results show that 10 systems exhibit the O'Connell effect, which is well explained by introducing a dark spot on the primary star. 11 systems have mass ratios below 0.15, classifying them as extremely low mass ratio contact binaries, making them promising candidates for binary mergers. Six systems have primary star temperatures exceeding 10,000 K, classifying them as early-type contact binaries. The absolute physical parameters reveal that two contact binary systems contain massive stellar components, which may eventually evolve into compact binary star systems. To compare the effects of different galactic environments on the evolution of contact binaries, we constructed evolutionary diagrams for these 30 targets and for contact binaries in the Milky Way. The results show that, at the same mass, our targets have larger radii, higher luminosities, and lower orbital angular momenta than contact binaries in the Milky Way, indicating that they are at more advanced evolutionary stages. This may be attributed to the higher metallicity in M31 compared to the Milky Way. △ Less

Submitted 6 October, 2025; v1 submitted 3 October, 2025; originally announced October 2025.

Comments: 16pages, 5 figures, 3 tables, accepted by ApJ

Artificial Satellite Trails Detection Using U-Net Deep Neural Network and Line Segment Detector Algorithm

Authors: Xiaohan Chen, Hongrui Gu, Cunshi Wang, Haiyang Mu, Jie Zheng, Junju Du, Jing Ren, Zhou Fan, Jing Li

Abstract: With the rapid increase in the number of artificial satellites, astronomical imaging is experiencing growing interference. When these satellites reflect sunlight, they produce streak-like artifacts in photometry images. Such satellite trails can introduce false sources and cause significant photometric errors. As a result, accurately identifying the positions of satellite trails in observational d… ▽ More With the rapid increase in the number of artificial satellites, astronomical imaging is experiencing growing interference. When these satellites reflect sunlight, they produce streak-like artifacts in photometry images. Such satellite trails can introduce false sources and cause significant photometric errors. As a result, accurately identifying the positions of satellite trails in observational data has become essential. In this work, we propose a satellite trail detection model that combines the U-Net deep neural network for image segmentation with the Line Segment Detector (LSD) algorithm. The model is trained on 375 simulated images of satellite trails, generated using data from the Mini-SiTian Array. Experimental results show that for trails with a signal-to-noise ratio (SNR) greater than 3, the detection rate exceeds 99. Additionally, when applied to real observational data from the Mini-SiTian Array, the model achieves a recall of 79.57 and a precision of 74.56. △ Less

Submitted 20 September, 2025; originally announced September 2025.

Comments: 15 pages, 7 figures, 2 tables, PASP accepted

The Stellar Abundances and Galactic Evolution Survey (SAGES). IV. Surface Gravity Estimation and Giant-Dwarf Separation with the DDO51 Filter

Authors: Qiqian Zhang, Zhou Fan, Gang Zhao, Ying Wu, Wei Wang, Kai Xiao, Hongrui Gu, Jie Zheng, Jingkun Zhao, Chun Li, Yuqin Chen, Haibo Yuan, Haining Li, Kefeng Tan, Yihan Song, Ali Luo, Nan Song, Yujuan Liu, Yaqian Wu

Abstract: Reliable estimation of stellar surface gravity (log $g$) for a large sample is crucial for evaluating stellar evolution models and understanding galactic structure; However, it is not easy to accomplish due to the difficulty in gathering a large spectroscopic data set. Photometric sky survey using a specific filter, on the other hand, can play a substantial role in the assessment of log $g$. The S… ▽ More Reliable estimation of stellar surface gravity (log $g$) for a large sample is crucial for evaluating stellar evolution models and understanding galactic structure; However, it is not easy to accomplish due to the difficulty in gathering a large spectroscopic data set. Photometric sky survey using a specific filter, on the other hand, can play a substantial role in the assessment of log $g$. The Stellar Abundances and Galactic Evolution Survey (SAGES) utilizes eight filters to provide accurate stellar parameters for $\sim10^{7}$ stars, with its DDO51 intermediate-band filter specifically designed for robust log $g$ determination. In this work, the observed SAGES $u_{\rm SC}$ and $v_{\rm SAGES}$ photometry, the synthetic photometry in $g$, $r$, $i$, and DDO51 bands derived from \textit{Gaia} XP spectra are employed to investigate the importance of the DDO51 filter in the determination of log $g$. We applied machine-learning-based extinction correction and employed XGBoost models, trained on stellar parameters from LAMOST, to predict log $g$ using photometric data. By comparing model predicted log $g$ with LAMOST values, we find that including DDO51 filter improve the accuracies of log $g$ estimates by 21.0\% (from 0.224\,dex to 0.177\,dex) overall, and by 26.5\% (from 0.302\,dex to 0.222\,dex ) for GK-type stars, as compared to those obtained without DDO51. The DDO51 filter is also validated to be particularly effective for metal-poor stars ([Fe/H]$<$-1.0), where it significantly mitigates systematic biases. Our findings highlight the diagnostic power of the SAGES DDO51 filter, providing enhanced stellar characterization vital for future in-depth studies of the Milky Way. △ Less

Submitted 18 September, 2025; originally announced September 2025.

Comments: 22 pages, 13 figures. Accepted for publication in ApJ

An all-sky 3D dust map Based on Gaia and LAMOST

Authors: Tao Wang, Haibo Yuan, Bingqiu Chen, Maosheng Xiang, Ruoyi Zhang, Bowen Huang, Hongrui Gu, Shuaicong Wang, Jiawei Li

Abstract: We present a comprehensive 3D dust reddening map covering the entire Milky Way, constructed by combining reddening estimates based on LAMOST low-resolution spectra (E(B$-$V)$_{\rm LAMOST}$) with those derived from $Gaia$ XP spectra (E(B$-$V)$_{\rm XP}$), along with revised $Gaia$ distances. E(B$-$V)$_{\rm LAMOST}$ values of $\sim$ 4.6 million unique sources were obtained with the standard-pair ana… ▽ More We present a comprehensive 3D dust reddening map covering the entire Milky Way, constructed by combining reddening estimates based on LAMOST low-resolution spectra (E(B$-$V)$_{\rm LAMOST}$) with those derived from $Gaia$ XP spectra (E(B$-$V)$_{\rm XP}$), along with revised $Gaia$ distances. E(B$-$V)$_{\rm LAMOST}$ values of $\sim$ 4.6 million unique sources were obtained with the standard-pair analysis using LAMOST DR11 stellar parameters and synthesized $B/V$-band photometry from $Gaia$ XP spectra, showing a typical precision of $\sim$ 0.01 mag. The E(B$-$V)$_{\rm XP}$ from the catalog of \citet{zhang2023}, which was derived using forward modeling of $Gaia$ XP spectra, were cross-validated with E(B$-$V)$_{\rm LAMOST}$, leading to the selection of $\sim$ 150 million high-reliability measurements. The combined dataset achieves a median precision of $\sim$ 0.03 mag for E(B$-$V). To model the reddening -- distance relationship along various lines-of-sight, we implemented a parametric approach that accounts for contributions from the local bubble, diffuse interstellar-medium, and multiple potential molecular clouds. The sky was adaptively partitioned based on stellar density, resulting in angular resolutions ranging from 3.4$^{\prime}$ to 58$^{\prime}$, with about half of the sky having a resolution better than 6.9$^{\prime}$. The reddening precision of our 3D map for individual stars reaches $\sim$ 0.01 mag in most regions at $|b| > 20^\circ$, but degrades to 0.01-0.05 mag at $|b| < 20^\circ$. The map reaches a maximum distance of 3-5 kpc in high-extinction regions with $|b| < 5^\circ$, and extends to 10-15 kpc elsewhere. An interactive platform and Python package have been developed for utilization of the 3D dust map. Available online: https://nadc.china-vo.org/data/dustmaps/. △ Less

Submitted 9 September, 2025; originally announced September 2025.

Comments: 25 pages, 20 figures, 3 tables, published in The Astrophysical Journal Supplement Series (ApJS)

Journal ref: ApJS, 280, 15 (2025)

Unlocking New Paths for Science with Extreme-Mass-Ratio Inspirals: Machine Learning-Enhanced MCMC for Accurate Parameter Inversion

Authors: Bo Liang, Chang Liu, Hanlin Song, Zhenwei Lyu, Minghui Du, Peng Xu, Ziren Luo, Sensen He, Haohao Gu, Tianyu Zhao, Manjia Liang Yuxiang Xu, Li-e Qiang, Mingming Sun, Wei-Liang Qian

Abstract: The detection of gravitational waves from extreme-mass-ratio inspirals (EMRIs) in space-borne antennas like Taiji and LISA promises deep insights into strong-field gravity and black hole physics. However, the complex, highly degenerate, and non-convex likelihood landscapes characteristic of EMRI parameter spaces pose severe challenges for conventional Markov chain Monte Carlo (MCMC) methods. Under… ▽ More The detection of gravitational waves from extreme-mass-ratio inspirals (EMRIs) in space-borne antennas like Taiji and LISA promises deep insights into strong-field gravity and black hole physics. However, the complex, highly degenerate, and non-convex likelihood landscapes characteristic of EMRI parameter spaces pose severe challenges for conventional Markov chain Monte Carlo (MCMC) methods. Under realistic instrumental noise and broad priors, these methods demand impractical computational costs but are prone to becoming trapped in local maxima, leading to biased and unreliable parameter estimates. To address this, we introduce Flow-Matching Markov Chain Monte Carlo (FM-MCMC), a novel Bayesian framework that integrates continuous normalizing flows (CNFs) with parallel tempering MCMC (PTMCMC). By generating high-likelihood regions via CNFs and refining them through PTMCMC, FM-MCMC enables robust exploration of the nontrivial parameter spaces, while achieving orders-of-magnitude improvement in computational efficiency and, more importantly, ensuring statistically reliable and unbiased inference. By enabling real-time, unbiased parameter inference, FM-MCMC could unlock the full scientific potential of EMRI observations, and would serve as a scalable pipeline for precision gravitational-wave astronomy. △ Less

Submitted 26 August, 2025; v1 submitted 1 August, 2025; originally announced August 2025.

All-sky search for individual Primordial Black Hole bursts with LHAASO

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

Abstract: Primordial Black Holes~(PBHs) are hypothetical black holes with a wide range of masses that formed in the early universe. As a result, they may play an important cosmological role and provide a unique probe of the early universe. A PBH with an initial mass of approximately $10^{15}$~g is expected to explode today in a final burst of Hawking radiation. In this work, we conduct an all-sky search for… ▽ More Primordial Black Holes~(PBHs) are hypothetical black holes with a wide range of masses that formed in the early universe. As a result, they may play an important cosmological role and provide a unique probe of the early universe. A PBH with an initial mass of approximately $10^{15}$~g is expected to explode today in a final burst of Hawking radiation. In this work, we conduct an all-sky search for individual PBH burst events using the data collected from March 2021 to July 2024 by the Water Cherenkov Detector Array of the Large High Altitude Air Shower Observatory (LHAASO). Three PBH burst durations, 10~s, 20~s, and 100~s, are searched, with no significant PBH bursts observed. The upper limit on the local PBH burst rate density is set to be as low as 181~pc$^{-3}$~yr$^{-1}$ at 99$\%$ confidence level, representing the most stringent limit achieved to date. △ Less

Submitted 2 November, 2025; v1 submitted 30 May, 2025; originally announced May 2025.

Comments: 8 pages, 2 figures

First Identification and Precise Spectral Measurement of the Proton Component in the Cosmic-Ray `Knee'

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

Abstract: We report the first high-purity identification of cosmic-ray (CR) protons and a precise measurement of their energy spectrum from 0.15 to 12 PeV using the Large High Altitude Air Shower Observatory (LHAASO). Abundant event statistics, combined with the simultaneous detection of electrons/photons, muons, and Cherenkov light in air showers, enable spectroscopic measurements with statistical and syst… ▽ More We report the first high-purity identification of cosmic-ray (CR) protons and a precise measurement of their energy spectrum from 0.15 to 12 PeV using the Large High Altitude Air Shower Observatory (LHAASO). Abundant event statistics, combined with the simultaneous detection of electrons/photons, muons, and Cherenkov light in air showers, enable spectroscopic measurements with statistical and systematic accuracy comparable to satellite data at lower energies. The proton spectrum shows significant hardening relative to low-energy extrapolations, culminating at 3 PeV, followed by sharp softening. This distinct spectral structure - closely aligned with the knee in the all-particle spectrum - points to the emergence of a new CR component at PeV energies, likely linked to the dozens of PeVatrons recently discovered by LHAASO, and offers crucial clues to the origin of Galactic cosmic rays. △ Less

Submitted 20 May, 2025; originally announced May 2025.

The Mini-SiTian Array: Evaluation Camera System

Authors: Yu Zhang, Lin Du, Yi Hu, Yang Huang, Ying Wu, Min He, Hongrui Gu, Haiyang Mu, Xunhao Chen, Hong Wu

Abstract: The Mini-SiTian project, which is the pathfinder for the SiTian project, utilizes three 30 cm telescopes equipped with commercial CMOS cameras (ZWO ASI6200MM Pro) to simulate large-area time-domain survey. Due to the avoidance of the traditional mechanical shutter, the CMOS camera is favorable in time-domain survey projects. In the future, the SiTian telescope array will employ a two-by-two scient… ▽ More The Mini-SiTian project, which is the pathfinder for the SiTian project, utilizes three 30 cm telescopes equipped with commercial CMOS cameras (ZWO ASI6200MM Pro) to simulate large-area time-domain survey. Due to the avoidance of the traditional mechanical shutter, the CMOS camera is favorable in time-domain survey projects. In the future, the SiTian telescope array will employ a two-by-two scientific-grade mosaic CMOS camera to survey a 10,000-degree square area every 30 minutes. Therefore, the performance of CMOS directly determines the detection capability of SiTian telescopes for transient sources, and a comprehensive understanding of the performance of CMOS cameras is crucial. In this research, laboratory testing was conducted to thoroughly evaluate three cameras by assessing several critical parameters, including bias stability, dark current, pixel anomalies, linearity, gain, and read noise. We find exceptional short-term bias stability with standard deviations below 0.02 ADU, negligible dark current of approximately 0.002 e$^{-}$ pixel$^{-1}$ s$^{-1}$ at $0^\circ\text{C}$, and excellent linearity with nonlinearity consistently below $\pm$ 0.5\%, and a small proportion (0.06\% to 0.08\%) of pixels with anomalous responses. Furthermore, our analysis demonstrates uniform gain values across all cameras, ranging from 0.252 to 0.255 e$^{-}$ ADU$^{-1}$, with low readout noise, measured to be below 1.6 e$^{-}$ using conventional methods. We also propose a novel method for pixel-level gain and read noise calculation for CMOS sensors, which revealed a narrow gain distribution and a low median read noise of 1.028 e$^-$ for one of the cameras. The laboratory testing of the ZWO ASI6200MM Pro cameras indicates their potential to meet the requirements of time-domain surveys for the Mini-SiTian project. △ Less

Submitted 2 April, 2025; originally announced April 2025.

Comments: 20 pages, 10 figures, accepted for publication in a special issue of Research in Astronomy and Astrophysics on the Mini-SiTian Array

The Mini-SiTian Array: Imaging Processing Pipeline

Authors: Kai Xiao, Zhirui Li, Yang Huang, Jie Zheng, Haibo Yuan, Junju Du, Linying Mi, Hongrui Gu, Yongkang Sun, Bowen Zhang, Shunxuan He, Henggeng Han, Min He, Ruifeng Shi, Yu Zhang, Chuanjie Zheng, Zexi Niu, Guiting Tian, Hu Zou, Yongna Mao, Hong Wu, Jifeng Liu

Abstract: As a pathfinder of the SiTian project, the Mini-SiTian (MST) array, employed three commercial CMOS cameras, represents a next-generation, cost-effective optical time-domain survey project. This paper focuses primarily on the precise data processing pipeline designed for wide-field, CMOS-based devices, including the removal of instrumental effects, astrometry, photometry, and flux calibration. When… ▽ More As a pathfinder of the SiTian project, the Mini-SiTian (MST) array, employed three commercial CMOS cameras, represents a next-generation, cost-effective optical time-domain survey project. This paper focuses primarily on the precise data processing pipeline designed for wide-field, CMOS-based devices, including the removal of instrumental effects, astrometry, photometry, and flux calibration. When applying this pipeline to approximately 3000 observations taken in the Field 02 (f02) region by MST, the results demonstrate a remarkable astrometric precision of approximately 70--80\,mas (about 0.1\,pixel), an impressive calibration accuracy of approximately 1\,mmag in the MST zero points, and a photometric accuracy of about 4\,mmag for bright stars. Our studies demonstrate that MST CMOS can achieve photometric accuracy comparable to that of CCDs, highlighting the feasibility of large-scale CMOS-based optical time-domain surveys and their potential applications for cost optimization in future large-scale time-domain surveys, like the SiTian project. △ Less

Submitted 2 April, 2025; originally announced April 2025.

Comments: 24 pages, 17 figures, 1 table, accepted for publication in a special issue of Research in Astronomy and Astrophysics on the Mini-SiTian Array, see main results in Figures 11, 12, and 15

Searching for kilonova with the SiTian prototype telescope

Authors: Zhirui Li, Hongrui Gu, Yongkang Sun, Yang Huang, Youjun Lu, Jifeng Liu

Abstract: We simulate the optical searching of gravitational-wave electromagnetic counterpart of the binary neutron star (BNS) merger event (i.e., a kilonova) using the ground based {\it SiTian} project prototype telescope with a 5-min limiting magnitude of 22.0, triggered by LIGO and Virgo gravitational wave detectors during the ongoing O4 run. Our simulations show that an average of 0.17-0.25 kilonova eve… ▽ More We simulate the optical searching of gravitational-wave electromagnetic counterpart of the binary neutron star (BNS) merger event (i.e., a kilonova) using the ground based {\it SiTian} project prototype telescope with a 5-min limiting magnitude of 22.0, triggered by LIGO and Virgo gravitational wave detectors during the ongoing O4 run. Our simulations show that an average of 0.17-0.25 kilonova events can be observed over the entire O4 period of $\sim 2$ years in the most optimistic case we set, while no kilonova can be detected in other cases. We note that it is beneficial for {\it SiTian}'s kilonova searching by extending the exposure time to gain deeper limiting magnitude despite the rapid decline of kilonova luminosity. △ Less

Submitted 2 April, 2025; originally announced April 2025.

Comments: 16 pages, 7 figures. Accepted for publication in a special issue of Research in Astronomy and Astrophysics on the Mini-SiTian Array

The Mini-SiTian Array: first-two-year operation

Authors: Min He, Hong Wu, Liang Ge, Jian-feng Tian, Zheng Wang, Hai-yang Mu, Yu Zhang, Yang Huang, Jie Zheng, Zhou Fan, Zheng-yang Li, Hong-hui Gu, Heng-geng Han, Kai Xiao, Zhi-rui Li, Jun-jie Jin, Bei-chuan Wang, Jun Ma, Jin-hang Zou, Ying Wu, Jiu-peng Guo, Li-guo Fang, Zhi-gang Hou, Bo-wen Zhang, Yun-fei Xu , et al. (48 additional authors not shown)

Abstract: The SiTian project, designed to utilize 60 telescopes distributed across multiple sites in China, is a next-generation time-domain survey initiative. As a pathfinder for the SiTian project, the Mini-SiTian (MST) has been proposed and implemented to test the SiTian's brain and data pipeline, and to evaluate the feasibility of its technology and science cases. Mounted at the Xinglong Observatory, th… ▽ More The SiTian project, designed to utilize 60 telescopes distributed across multiple sites in China, is a next-generation time-domain survey initiative. As a pathfinder for the SiTian project, the Mini-SiTian (MST) has been proposed and implemented to test the SiTian's brain and data pipeline, and to evaluate the feasibility of its technology and science cases. Mounted at the Xinglong Observatory, the MST project comprises three 30 cm telescopes and has been operated since Nov. 2022. Each telescope of the MST possesses a large field of view, covering $2.29^{\circ}$ $\times$ $1.53^{\circ}$ FOV, and is equipped with $g'$, $r'$ and $i'$ filters, respectively. Acting as the pioneer of the forthcoming SiTian project, the MST is dedicated to the discovery of variable stars, transients, and outburst events, and has already obtained some interesting scientific results. In this paper, we will summarize the first-two-year operation of the MST project. △ Less

Submitted 2 April, 2025; originally announced April 2025.

Comments: 10 pages, 11 figures, Accepted for publication in a special issue of Research in Astronomy and Astrophysics on the Mini-SiTian Array

The Mini-SiTian Array: Light Curves Analysis of Asteroids

Authors: Zhaoxing Liu, Jian Gao, Hongrui Gu, Yang Huang, Shaoming Hu, Hu Zou, Keyu Xing, Hao Huang, Zehao Zhang

Abstract: The SiTian project, with its vast field of view, will become an ideal platform for asteroid scientific research. In this study, we develop a pipeline to analyze the photometry of asteroids and derive their periods from the data collected by the SiTian pathfinder project Mini-SiTian (MST). The pipeline is applied to the MST f02 region, a MST test region with a sky area of… ▽ More The SiTian project, with its vast field of view, will become an ideal platform for asteroid scientific research. In this study, we develop a pipeline to analyze the photometry of asteroids and derive their periods from the data collected by the SiTian pathfinder project Mini-SiTian (MST). The pipeline is applied to the MST f02 region, a MST test region with a sky area of $2.29^{\circ} \times 1.53^{\circ}$. Rotation periods of 22 asteroids are derived by the obtained light curves analysis. Among them, there are 8 asteroids available in the Asteroid Lightcurve Photometry Database (ALCDEF), and 6 of them with more photometric points ($>$200) have similar period parameters as the ones in ALCDEF. Additionally, the periods for 14 of these asteroids are newly obtained and are not listed in ALCDEF. This study demonstrates the feasibility of asteroid photometric research by the SiTian project. It shows that future observations from the SiTian project will provide even more photometry of asteroids, significantly increasing the number of available light curves. The potential vast photometric data of asteroids will help us to further understand the physics of asteroids, their material composition, and the formation and evolution of the solar system. △ Less

Submitted 2 April, 2025; originally announced April 2025.

Comments: 14 pages, 7 figures, 4 tables. Accepted for publication in a special issue of Research in Astronomy and Astrophysics on the Mini-SiTian Array

The Mini-SiTian Array: White Paper

Authors: Henggeng Han, Yang Huang, Beichuan Wang, Yongkang Sun, Cunshi Wang, Zhirui Li, Junjie Jin, Ningchen Sun, Kai Xiao, Min He, Hongrui Gu, Zexi Niu, Hong Wu, Jifeng Liu

Abstract: This paper outlines the scientific goals and observational strategies of the Mini-SiTian array. Mounted at Xinglong Observatory, the Mini-SiTian array consists of three 30 cm telescopes and has been in operation since 2022. The large field of view, combined with the capability for multi-band photometric observations, enables the Mini-SiTian array to perform rapid follow-up observations to identify… ▽ More This paper outlines the scientific goals and observational strategies of the Mini-SiTian array. Mounted at Xinglong Observatory, the Mini-SiTian array consists of three 30 cm telescopes and has been in operation since 2022. The large field of view, combined with the capability for multi-band photometric observations, enables the Mini-SiTian array to perform rapid follow-up observations to identify optical counterparts of gravitational waves, capture the early light curves of tidal disruption events and supernovae, and monitor stellar flares, Be star outbursts, and cataclysmic variable stars, although its limiting magnitude is not very deep. By collaborating with the Xinglong 2.16-m telescope and leveraging a real-time image processing pipeline, simultaneous photometric and spectroscopic observations could be performed to reveal their underlying physical mechanisms. The observational and research experience provide critical guidance for the implementation of the full-scale SiTian project in the future. △ Less

Submitted 2 April, 2025; originally announced April 2025.

Comments: 11 pages, 2 figures, Accepted for publication in a special issue of Research in Astronomy and Astrophysics on the Mini-SiTian Array

The Mini-SiTian Array: the mini-SiTian Realtime Image Processing pipeline (STRIP)

Authors: Hongrui Gu, Yang Huang, Yongkang Sun, Kai Xiao, Zhirui Li, Beichuan Wang, Zhou Fan, Chuanjie Zheng, Henggeng Han, Hu Zou, Wenxiong Li, Hong Wu, Jifeng Liu

Abstract: This paper provides a comprehensive introduction to the Mini-SiTian Real-Time Image Processing pipeline (STRIP) and evaluates its operational performance. The STRIP pipeline is specifically designed for real-time alert triggering and light curve generation for transient sources. By applying the STRIP pipeline to both simulated and real observational data of the Mini-SiTian survey, it successfully… ▽ More This paper provides a comprehensive introduction to the Mini-SiTian Real-Time Image Processing pipeline (STRIP) and evaluates its operational performance. The STRIP pipeline is specifically designed for real-time alert triggering and light curve generation for transient sources. By applying the STRIP pipeline to both simulated and real observational data of the Mini-SiTian survey, it successfully identified various types of variable sources, including stellar flares, supernovae, variable stars, and asteroids, while meeting requirements of reduction speed within 5 minutes. For the real observational dataset, the pipeline detected 1 flare event, 127 variable stars, and 14 asteroids from three monitored sky regions. Additionally, two datasets were generated: one, a real-bogus training dataset comprising 218,818 training samples, and the other, a variable star light curve dataset with 421 instances. These datasets will be used to train machine learning algorithms, which are planned for future integration into STRIP. △ Less

Submitted 2 April, 2025; originally announced April 2025.

Comments: 17 pages, 9 figures. Accepted for publication in a special issue of Research in Astronomy and Astrophysics on the Mini-SiTian Array

The Mini-SiTian Array: real-bogus classification using deep learning

Authors: Jing-Hang Shi, Hong-Rui Gu, Yang Huang, Yan-Xia Zhang, Peng-Liang Du

Abstract: The Mini-SiTian (MST) project is a pathfinder for China's next-generation large-scale time-domain survey, SiTian, aimed at discovering variable stars, transients, and explosive events. MST generates hundreds of thousands of transient alerts every night, approximately 99\% of which are false alarms, posing a significant challenge to its scientific goals. To mitigate the impact of false positives, w… ▽ More The Mini-SiTian (MST) project is a pathfinder for China's next-generation large-scale time-domain survey, SiTian, aimed at discovering variable stars, transients, and explosive events. MST generates hundreds of thousands of transient alerts every night, approximately 99\% of which are false alarms, posing a significant challenge to its scientific goals. To mitigate the impact of false positives, we propose a deep learning-based solution and systematically evaluate thirteen convolutional neural networks. The results show that ResNet achieves exceptional specificity (99.70\%), EfficientNet achieves the highest recall rate (98.68\%), and DenseNet provides balanced performance with a recall rate of 94.55\% and specificity of 98.66\%. Leveraging these complementary strengths, we developed a bagging-based ensemble classifier that integrates ResNet18, DenseNet121, and EfficientNet\_B0 using a soft voting strategy. This classifier achieved the best AUC value (0.9961) among all models, with a recall rate of 95.37\% and specificity of 99.25\%. It has now been successfully deployed in the MST real-time data processing pipeline. Validation using 5,000 practically processed samples with a classification threshold of 0.798 showed that the classifier achieved 88.31\% accuracy, 91.89\% recall rate, and 99.82\% specificity, confirming its effectiveness and robustness under real application conditions. △ Less

Submitted 2 April, 2025; originally announced April 2025.

Comments: 15 pages, 7 figures, 3 tables. Accepted for publication in a special issue of Research in Astronomy and Astrophysics on the Mini-SiTian Array

Calibration of Complementary Metal-oxide-semiconductor Sensor-based Photometry to a Few-millimagnitude Precision: The Case of the Mini-SiTian Array

Authors: Kai Xiao, Yang Huang, Haibo Yuan, Zhirui Li, Yongkang Sun, Timothy C. Beers, Min He, Jifeng Liu, Hong Wu, Yongna Mao, Bowen Huang, Mingyang Ma, Chuanjie Zheng, Hongrui Gu, Beichuan Wang, Lin Yang, Shuai Xu

Abstract: We present a pioneering achievement in the high-precision photometric calibration of CMOS-based photometry, by application of the Gaia BP/RP (XP) spectra-based synthetic photometry (XPSP) method to the mini-SiTian array (MST) photometry. Through 79 repeated observations of the $\texttt{f02}$ field on the night, we find good internal consistency in the calibrated MST $G_{\rm MST}$-band magnitudes f… ▽ More We present a pioneering achievement in the high-precision photometric calibration of CMOS-based photometry, by application of the Gaia BP/RP (XP) spectra-based synthetic photometry (XPSP) method to the mini-SiTian array (MST) photometry. Through 79 repeated observations of the $\texttt{f02}$ field on the night, we find good internal consistency in the calibrated MST $G_{\rm MST}$-band magnitudes for relatively bright stars, with a precision of about 4\,mmag for $G_{\rm MST}\sim 13$. Results from more than 30 different nights (over 3100 observations) further confirm this internal consistency, indicating that the 4\,mmag precision is stable and achievable over timescales of months. An independent external validation using spectroscopic data from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR10 and high-precision photometric data using CCDs from Gaia DR3 reveals a zero-point consistency better than 1\,mmag. Our results clearly demonstrate that CMOS photometry is on par with CCD photometry for high-precision results, highlighting the significant capabilities of CMOS cameras in astronomical observations, especially for large-scale telescope survey arrays. △ Less

Submitted 16 March, 2025; originally announced March 2025.

Comments: 8 pages, 6 figures, ApJL accepted, see main results in Figures 4

"Frog-eyes" in Astronomy: Monitoring Binary Radial Velocity Variations Through A Pair of Narrow-Band Filters

Authors: Chuanjie Zheng, Yang Huang, Jifeng Liu, Hongrui Gu, Hong Wu, Youjun Lu, Yongkang Sun, Henggeng Han, Song Wang, Timothy C. Beers, Kai Xiao, Zhirui Li, Boweng Zhang, Yongna Mao, Zhengyang Li, Hangxin Ji

Abstract: Spectroscopic observations are a crucial step in driving major discoveries in the era of time-domain surveys. However, the pace of current spectroscopic surveys is increasingly unable to meet the demands of rapidly advancing large-scale time-domain surveys. To address this issue, we propose the ``Frog-eyes" system, which employs a pair of narrow-band filters: one positioned near a strong absorptio… ▽ More Spectroscopic observations are a crucial step in driving major discoveries in the era of time-domain surveys. However, the pace of current spectroscopic surveys is increasingly unable to meet the demands of rapidly advancing large-scale time-domain surveys. To address this issue, we propose the ``Frog-eyes" system, which employs a pair of narrow-band filters: one positioned near a strong absorption line to capture signals from Doppler shifts, and the other placed on the adjacent continuum to monitor intrinsic variations. The combination of observations from the two filters enables the extraction of radial velocity (RV) curves from a large sample of binary stars, and is particularly efficient for single-lined binaries (SB1), using photometric techniques. Comprehensive mock simulations on SB1 demonstrate that the binary orbital parameters can be precisely measured from the extracted RV curves for binary systems where the primary star has an effective temperature greater than 6000 K. With a typical ground-based photometric precision of approximately 0.3%, the uncertainties in the derived semi-amplitude K and eccentricity e are less than 10% and 0.1, respectively, for binary systems with K $\ge$ 30 km/s. These encouraging results are further validated by real observations of the hot subdwarf-white dwarf binary system HD 265435, using a non-specialized ``Frog-eyes" system installed on the Chinese 2.16m telescope. Once this system is properly installed on large-field-of-view survey telescopes, the rate of acquiring RV curves for binaries will approach their detection rate in leading time-domain photometric surveys. △ Less

Submitted 25 February, 2025; originally announced February 2025.

Comments: 17 pages, 11 figures, and 3 tables. Accepted for publication by SCPMA

Ultra-high-energy $γ$-ray emission associated with the tail of a bow-shock pulsar wind nebula

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

Abstract: In this study, we present a comprehensive analysis of an unidentified point-like ultra-high-energy (UHE) $γ$-ray source, designated as 1LHAASO J1740+0948u, situated in the vicinity of the middle-aged pulsar PSR J1740+1000. The detection significance reached 17.1$σ$ (9.4$σ$) above 25$\,$TeV (100$\,$TeV). The source energy spectrum extended up to 300$\,$TeV, which was well fitted by a log-parabola f… ▽ More In this study, we present a comprehensive analysis of an unidentified point-like ultra-high-energy (UHE) $γ$-ray source, designated as 1LHAASO J1740+0948u, situated in the vicinity of the middle-aged pulsar PSR J1740+1000. The detection significance reached 17.1$σ$ (9.4$σ$) above 25$\,$TeV (100$\,$TeV). The source energy spectrum extended up to 300$\,$TeV, which was well fitted by a log-parabola function with $N0 = (1.93\pm0.23) \times 10^{-16} \rm{TeV^{-1}\,cm^{-2}\,s^{-2}}$, $α= 2.14\pm0.27$, and $β= 1.20\pm0.41$ at E0 = 30$\,$TeV. The associated pulsar, PSR J1740+1000, resides at a high galactic latitude and powers a bow-shock pulsar wind nebula (BSPWN) with an extended X-ray tail. The best-fit position of the gamma-ray source appeared to be shifted by $0.2^{\circ}$ with respect to the pulsar position. As the (i) currently identified pulsar halos do not demonstrate such offsets, and (ii) centroid of the gamma-ray emission is approximately located at the extension of the X-ray tail, we speculate that the UHE $γ$-ray emission may originate from re-accelerated electron/positron pairs that are advected away in the bow-shock tail. △ Less

Submitted 24 February, 2025; v1 submitted 21 February, 2025; originally announced February 2025.

Comments: Corrected spelling errors in several author names

Journal ref: The Innovation (2025), 100802

Broadband $γ$-ray spectrum of supernova remnant Cassiopeia A

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

Abstract: The core-collapse supernova remnant (SNR) Cassiopeia A (Cas A) is one of the brightest galactic radio sources with an angular radius of $\sim$ 2.5 $\arcmin$. Although no extension of this source has been detected in the $γ$-ray band, using more than 1000 days of LHAASO data above $\sim 0.8$ TeV, we find that its spectrum is significantly softer than those obtained with Imaging Air Cherenkov Telesc… ▽ More The core-collapse supernova remnant (SNR) Cassiopeia A (Cas A) is one of the brightest galactic radio sources with an angular radius of $\sim$ 2.5 $\arcmin$. Although no extension of this source has been detected in the $γ$-ray band, using more than 1000 days of LHAASO data above $\sim 0.8$ TeV, we find that its spectrum is significantly softer than those obtained with Imaging Air Cherenkov Telescopes (IACTs) and its flux near $\sim 1$ TeV is about two times higher. In combination with analyses of more than 16 years of \textit{Fermi}-LAT data covering $0.1 \, \mathrm{GeV} - 1 \, \mathrm{TeV}$, we find that the spectrum above 30 GeV deviates significantly from a single power-law, and is best described by a smoothly broken power-law with a spectral index of $1.90 \pm 0.15_\mathrm{stat}$ ($3.41 \pm 0.19_\mathrm{stat}$) below (above) a break energy of $0.63 \pm 0.21_\mathrm{stat} \, \mathrm{TeV}$. Given differences in the angular resolution of LHAASO-WCDA and IACTs, TeV $γ$-ray emission detected with LHAASO may have a significant contribution from regions surrounding the SNR illuminated by particles accelerated earlier, which, however, are treated as background by IACTs. Detailed modelling can be used to constrain acceleration processes of TeV particles in the early stage of SNR evolution. △ Less

Submitted 7 February, 2025; originally announced February 2025.

The Stellar Abundances and Galactic Evolution Survey (SAGES). II. Machine Learning-Based Stellar parameters for 21 million stars from the First Data Release

Authors: Hongrui Gu, Zhou Fan, Gang Zhao, Yang Huang, Timothy C. Beers, Wei Wang, Jie Zheng, Jingkun Zhao, Chun Li, Yuqin Chen, Haibo Yuan, Haining Li, Kefeng Tan, Yihan Song, Ali Luo, Nan Song, Yujuan Liu

Abstract: Stellar parameters for large samples of stars play a crucial role in constraining the nature of stars and stellar populations in the Galaxy. An increasing number of medium-band photometric surveys are presently used in estimating stellar parameters. In this study, we present a machine-learning approach to derive estimates of stellar parameters, including [Fe/H], logg, and Teff, based on a combinat… ▽ More Stellar parameters for large samples of stars play a crucial role in constraining the nature of stars and stellar populations in the Galaxy. An increasing number of medium-band photometric surveys are presently used in estimating stellar parameters. In this study, we present a machine-learning approach to derive estimates of stellar parameters, including [Fe/H], logg, and Teff, based on a combination of medium-band and broad-band photometric observations. Our analysis employs data primarily sourced from the SAGE Survey , which aims to observe much of the Northern Hemisphere. We combine the $uv$-band data from SAGES DR1 with photometric and astrometric data from Gaia EDR3, and apply the random forest method to estimate stellar parameters for approximately 21 million stars. We are able to obtain precisions of 0.09 dex for [Fe/H], 0.12 dex for logg, and 70 K for Teff. Furthermore, by incorporating 2MASS and WISE infrared photometric and GALEX ultraviolet data, we are able to achieve even higher precision estimates for over 2.2 million stars. These results are applicable to both giant and dwarf stars. Building upon this mapping, we construct a foundational dataset for research on metal-poor stars, the structure of the Milky Way, and beyond. With the forthcoming release of additional bands from SAGE Survey such DDO51 and H-alpha, this versatile machine learning approach is poised to play an important role in upcoming surveys featuring expanded filter sets △ Less

Submitted 5 February, 2025; originally announced February 2025.

Comments: Accepted by ApJS.12 pages, 12 figures, 3 tables

The Mini-SiTian Array: Prospects of Searching for Tidal Disruption Events

Authors: Bei-chuan Wang, Jun-jie Jin, Yu Zhang, Yanan Wang, Song Wang, Hong-rui Gu, Min He, Hai-yang Mu, Kai Xiao, Zhi-rui Li, Zhou Fan, Liang Ge, Jian-feng Tian, Yang Huang, Jie Zheng, Hong Wu, Jifeng Liu

Abstract: We assess the detectability of tidal disruption events (TDEs) using mock observations from the Mini-SiTian array. We select 100 host galaxy samples from a simulated galaxy catalog based on specific criteria such as redshift, BH mass, and event rate. Taking into account the site conditions and survey strategy, we simulate observations over a 440 deg$^2$ field. The results indicate that… ▽ More We assess the detectability of tidal disruption events (TDEs) using mock observations from the Mini-SiTian array. We select 100 host galaxy samples from a simulated galaxy catalog based on specific criteria such as redshift, BH mass, and event rate. Taking into account the site conditions and survey strategy, we simulate observations over a 440 deg$^2$ field. The results indicate that $0.53\pm 0.73$ TDEs can be detected per year when observing in both $g$ and $r$ bands with 300-second exposures every 3 days. Applying this method to the SiTian project, we expect to discover approximately 204 TDEs annually, heralding a new era in TDE science. △ Less

Submitted 2 April, 2025; v1 submitted 16 January, 2025; originally announced January 2025.

Comments: 19 pages, 9 figures, accepted for publication in a special issue of Research in Astronomy and Astrophysics on the Mini-SiTian Array

The S-PLUS Ultra-Short Survey: Photometric Re-calibration with the BEst STar Database

Authors: Xiaolu Li, Kai Xiao, Yang Huang, Haibo Yuan, Yanke Tang, Timothy C. Beers, Bowen Huang, Mingyang Ma, Pedro K. Humire, Alvaro Alvarez-Candal, Federico Sestito, Ning Gai, Yongna Mao, Hongrui Gu, Zhenzhao Tao, Lin Yang, Shuai Xu, Rong Hu

Abstract: We present an independent validation and comprehensive re-calibration of S-PLUS Ultra-Short Survey (USS) DR1 12-band photometry using about 30,000--70,000 standard stars from the BEst STar (BEST) database. We identify spatial variation of zero-point offsets, up to 30--40\,mmag for blue filters ($u$, $J0378$, $J0395$) and 10\,mmag for others, predominantly due to the higher uncertainties of the tec… ▽ More We present an independent validation and comprehensive re-calibration of S-PLUS Ultra-Short Survey (USS) DR1 12-band photometry using about 30,000--70,000 standard stars from the BEst STar (BEST) database. We identify spatial variation of zero-point offsets, up to 30--40\,mmag for blue filters ($u$, $J0378$, $J0395$) and 10\,mmag for others, predominantly due to the higher uncertainties of the technique employed in the original USS calibration. Moreover, we detect large- and medium-scale CCD position-dependent systematic errors, up to 50\,mmag, primarily caused by different aperture and flat-field corrections. We then re-calibrate the USS DR1 photometry by correcting the systematic shifts for each tile using second-order two-dimensional polynomial fitting combined with a numerical stellar flat-field correction method. The re-calibrated results from the XPSP and the SCR standards are consistent within 6\,mmag in the USS zero-points, demonstrating both the typical precision of re-calibrated USS photometry and a sixfold improvement in USS zero-point precision. Further validation using SDSS and Pan-STARRS1, as well as LAMOST DR10 and Gaia photometry, also confirms this precision for the re-calibrated USS photometry. Our results clearly demonstrate the capability and the efficiency of the BEST database in improving calibration precision to the milli-magnitude level for wide-field photometric surveys. The re-calibrated USS DR1 photometry is publicly available: doi: 10.12149/101503 (https://nadc.china-vo.org/res/r101504/). △ Less

Submitted 8 January, 2025; v1 submitted 6 January, 2025; originally announced January 2025.

Comments: 28 pages, 25 figures, ApJS accepted, see main results in Figures 4 and 7

Photometric Stellar Parameters for 195,478 Kepler Input Catalog (KIC) Stars

Authors: Bowen Zhang, Yang Huang, Timothy C. Beers, Kai Xiao, Jifeng Liu, Lei Jia, Henggeng Han, Zhirui Li, Chuanjie Zheng, Yongkang Sun, Ruifeng Shi, Hongrui Gu

Abstract: The stellar atmospheric parameters and physical properties of stars in the Kepler Input Catalog (KIC) are of great significance for the study of exoplanets, stellar activity, and asteroseismology. However, despite extensive effort over the past decades, accurate spectroscopic estimates of these parameters are available for only about half of the stars in the full KIC catalog. In our work, by train… ▽ More The stellar atmospheric parameters and physical properties of stars in the Kepler Input Catalog (KIC) are of great significance for the study of exoplanets, stellar activity, and asteroseismology. However, despite extensive effort over the past decades, accurate spectroscopic estimates of these parameters are available for only about half of the stars in the full KIC catalog. In our work, by training relationships between photometric colors and spectroscopic stellar parameters from Gaia DR3, the Kepler Issac-Newton Survey, LAMOST DR10, and APOGEE DR17, we have obtained atmospheric-parameter estimates for over 195,000 stars, accounting for 97$\%$ of the total sample of KIC stars. We obtain 1$σ$ uncertainties of 0.1 dex on metallicity [Fe/H], 100 K on effective temperature $T_{\rm eff}$, and 0.2 dex on surface gravity log $g$. In addition, based on these atmospheric parameters, we estimated the ages, masses, radii, and surface gravities of these stars using the commonly adopted isochrone-fitting approach. External comparisons indicate that the resulting precision for turn-off stars is 20$\%$ in age; for dwarf stars, it is 0.07 $M_{\odot}$ in mass, 0.05 $R_{\odot}$ in radius, and 0.12 dex in surface gravity; and for giant stars, it is 0.14 $M_{\odot}$ in mass, 0.73 $R_{\odot}$ in radius, and 0.11 dex in surface gravity. △ Less

Submitted 27 December, 2024; v1 submitted 22 December, 2024; originally announced December 2024.

Comments: 26 pages, 22 figures. To be published in ApJS. For associated catalogs of the result, see https://doi.org/10.5281/zenodo.14546166

Detection of two TeV gamma-ray outbursts from NGC 1275 by LHAASO

Authors: Zhen Cao, F. Aharonian, 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, T. L. Chen , et al. (254 additional authors not shown)

Abstract: The Water Cherenkov Detector Array (WCDA) is one of the components of Large High Altitude Air Shower Observatory (LHAASO) and can monitor any sources over two-thirds of the sky for up to 7 hours per day with >98\% duty cycle. In this work, we report the detection of two outbursts of the Fanaroff-Riley I radio galaxy NGC 1275 that were detected by LHAASO-WCDA between November 2022 and January 2023… ▽ More The Water Cherenkov Detector Array (WCDA) is one of the components of Large High Altitude Air Shower Observatory (LHAASO) and can monitor any sources over two-thirds of the sky for up to 7 hours per day with >98\% duty cycle. In this work, we report the detection of two outbursts of the Fanaroff-Riley I radio galaxy NGC 1275 that were detected by LHAASO-WCDA between November 2022 and January 2023 with statistical significance of 5.2~$σ$ and 8.3~$σ$. The observed spectral energy distribution in the range from 500 GeV to 3 TeV is fitted by a power-law with a best-fit spectral index of $α=-3.37\pm0.52$ and $-3.35\pm0.29$, respectively. The outburst flux above 0.5~TeV was ($4.55\pm 4.21)\times~10^{-11}~\rm cm^{-2}~s^{-1}$ and ($3.45\pm 1.78)\times~10^{-11}~\rm cm^{-2}~s^{-1}$, corresponding to 60\%, 45\% of Crab Nebula flux. Variation analysis reveals the variability time-scale of days at the TeV energy band. A simple test by one-zone synchrotron self-Compton model reproduces the data in the gamma-ray band well. △ Less

Submitted 18 April, 2025; v1 submitted 2 November, 2024; originally announced November 2024.

Comments: 11 pages, 8 figures, 3 tables

LHAASO detection of very-high-energy gamma-ray emission surrounding PSR J0248+6021

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 report the detection of an extended very-high-energy (VHE) gamma-ray source coincident with the location of middle-aged (62.4~\rm kyr) pulsar PSR J0248+6021, by using the LHAASO-WCDA data of live 796 days and LHAASO-KM2A data of live 1216 days. A significant excess of \gray induced showers is observed both by WCDA in energy bands of 1-25~\rm TeV and KM2A in energy bands of $>$ 25~\rm TeV with 7… ▽ More We report the detection of an extended very-high-energy (VHE) gamma-ray source coincident with the location of middle-aged (62.4~\rm kyr) pulsar PSR J0248+6021, by using the LHAASO-WCDA data of live 796 days and LHAASO-KM2A data of live 1216 days. A significant excess of \gray induced showers is observed both by WCDA in energy bands of 1-25~\rm TeV and KM2A in energy bands of $>$ 25~\rm TeV with 7.3 $σ$ and 13.5 $σ$, respectively. The best-fit position derived through WCDA data is R.A. = 42.06$^\circ \pm$ 0.12$^\circ$ and Dec. = 60.24$^\circ \pm $ 0.13$^\circ$ with an extension of 0.69$^\circ\pm$0.15$^\circ$ and that of the KM2A data is R.A.= 42.29$^\circ \pm $ 0.13$^\circ$ and Dec. = 60.38$^\circ \pm$ 0.07$^\circ$ with an extension of 0.37$^\circ\pm$0.07$^\circ$. No clear extended multiwavelength counterpart of this LHAASO source has been found from the radio band to the GeV band. The most plausible explanation of the VHE \gray emission is the inverse Compton process of highly relativistic electrons and positrons injected by the pulsar. These electrons/positrons are hypothesized to be either confined within the pulsar wind nebula or to have already escaped into the interstellar medium, forming a pulsar halo. △ Less

Submitted 3 December, 2024; v1 submitted 6 October, 2024; originally announced October 2024.

Comments: 12 pages, 10 figures, Accepted by Sci. China-Phys. Mech. Astron

J-PLUS: Beyond Spectroscopy III. Stellar Parameters and Elemental-abundance Ratios for Five Million Stars from DR3

Authors: Yang Huang, Timothy C. Beers, Kai Xiao, Haibo Yuan, Young Sun Lee, Hongrui Gu, Jihye Hong, Jifeng Liu, Zhou Fan, Paula Coelho, Patricia Cruz, F. J. Galindo-Guil, Simone Daflon, Fran Jiménez-Esteban, Javier Cenarro, David Cristóbal-Hornillos, Carlos Hernández-Monteagudo, Carlos López-Sanjuan, Antonio Marín-Franch, Mariano Moles, Jesús Varela, Héctor Vázquez Ramírez, Jailson Alcaniz, Renato Dupke, Alessandro Ederoclite , et al. (2 additional authors not shown)

Abstract: We present a catalog of stellar parameters (effective temperature $T_{\rm eff}$, surface gravity $\log g$, age, and metallicity [Fe/H]) and elemental-abundance ratios ([C/Fe], [Mg/Fe], and [$α$/Fe]) for some five million stars (4.5 million dwarfs and 0.5 million giants stars) in the Milky Way, based on stellar colors from the Javalambre Photometric Local Universe Survey (J-PLUS) DR3 and \textit{Ga… ▽ More We present a catalog of stellar parameters (effective temperature $T_{\rm eff}$, surface gravity $\log g$, age, and metallicity [Fe/H]) and elemental-abundance ratios ([C/Fe], [Mg/Fe], and [$α$/Fe]) for some five million stars (4.5 million dwarfs and 0.5 million giants stars) in the Milky Way, based on stellar colors from the Javalambre Photometric Local Universe Survey (J-PLUS) DR3 and \textit{Gaia} EDR3. These estimates are obtained through the construction of a large spectroscopic training set with parameters and abundances adjusted to uniform scales, and trained with a Kernel Principal Component Analysis. Owing to the seven narrow/medium-band filters employed by J-PLUS, we obtain precisions in the abundance estimates that are as good or better than derived from medium-resolution spectroscopy for stars covering a wide range of the parameter space: 0.10-0.20 dex for [Fe/H] and [C/Fe], and 0.05 dex for [Mg/Fe] and [$α$/Fe]. Moreover, systematic errors due to the influence of molecular carbon bands on previous photometric-metallicity estimates (which only included two narrow/medium-band blue filters) have now been removed, resulting in photometric-metallicity estimates down to [Fe/H] $\sim -4.0$, with typical uncertainties of 0.25 dex and 0.40 dex for dwarfs and giants, respectively. This large photometric sample should prove useful for the exploration of the assembly and chemical-evolution history of our Galaxy. △ Less

Submitted 4 August, 2024; originally announced August 2024.

Comments: 31 pages, 20 figures, 3 tables, accepeted by ApJ

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

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

H$α$ Reverberation Mapping from Broad-band Photometry of Dwarf Seyfert 1 Galaxy NGC 4395

Authors: Huapeng Gu, Xue-Bing Wu, Yuhan Wen, Qinchun Ma, Hengxiao Guo

Abstract: NGC 4395 is a dwarf Seyfert 1 galaxy with a possible intermediate-mass black hole of several $\rm{10^4}$ solar masses in its center. As a well-studied object, its broad line region size has been measured via H$\rmα$ time lag in numerous spectroscopic reverberation mapping (SRM) and narrow-band photometric reverberation mapping (PRM) campaigns. Here we present its H$\rmα$ time lag measurement using… ▽ More NGC 4395 is a dwarf Seyfert 1 galaxy with a possible intermediate-mass black hole of several $\rm{10^4}$ solar masses in its center. As a well-studied object, its broad line region size has been measured via H$\rmα$ time lag in numerous spectroscopic reverberation mapping (SRM) and narrow-band photometric reverberation mapping (PRM) campaigns. Here we present its H$\rmα$ time lag measurement using broad-band photometric data, with the application of our newly-developed ICCF-Cut method as well as the JAVELIN and $χ^2$ methods. utilizing the minute-cadence multi-band light curves obtained from the $\rm{2}$m FTN and $\rm{10.4}$m GTC telescopes in recent works, we measured its H$\rmα$ lag as approximately $40 \sim 90$ minutes from broad-band PRM. With the H$\rmα$ emission line velocity dispersion, we calculated its central black hole mass as $\rm M_{\rm BH} = (8\pm4) \times 10^3\, M_{\rm \odot}$. These results are comparable with previous results obtained by narrow-band PRM and SRM, providing further support to an intermediate-mass black hole in NGC 4395. In addition, our study also validates the ICCF-Cut as an effective method for broad-band PRM, which holds the potential for widespread application in the era of large multi-epoch, high-cadence photometric surveys. △ Less

Submitted 17 April, 2024; v1 submitted 17 April, 2024; originally announced April 2024.

Comments: 11 pages, 9 figures, accepted for publication in MNRAS

Searching for short-period variables in M31: method and catalogs

Authors: Hongrui Gu, Haibo Yuan, Subo Dong, Chenfa Zheng, Shenzhe Cui, Yi Ren, Haozhu Fu, Yang Huang, Zhou Fan

Abstract: Utilizing high-cadence and continuous g- and r-band data over three nights acquired from the 3.6-meter Canada France Hawaii Telescope (CFHT) aimed to find short-duration microlensing events, we conduct a systematic search for variables, transients, and asteroids across a $\sim1^\circ$ field of view of the Andromeda Galaxy (M 31). We present a catalog of 5859 variable stars, yielding the most exten… ▽ More Utilizing high-cadence and continuous g- and r-band data over three nights acquired from the 3.6-meter Canada France Hawaii Telescope (CFHT) aimed to find short-duration microlensing events, we conduct a systematic search for variables, transients, and asteroids across a $\sim1^\circ$ field of view of the Andromeda Galaxy (M 31). We present a catalog of 5859 variable stars, yielding the most extensive compilation of short-period variable sources of M 31. We also detected 19 flares, predominantly associated with foreground M dwarfs in the Milky Way. In addition, we discovered 17 previously unknown asteroid candidates, and we subsequently reported them to the Minor Planet Center. Lastly, we report a microlensing event candidate C-ML-1 and present a preliminary analysis. △ Less

Submitted 10 April, 2024; originally announced April 2024.

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.

Evidence for an Outer Component in the Continuum Reverberation Mapping of Active Galactic Nuclei

Authors: Yuanzhe Jiang, Xue-Bing Wu, Qinchun Ma, Huapeng Gu, Yuhan Wen

Abstract: The continuum reverberation mapping is widely used in studying accretion disk of active galactic nuclei (AGN). While some indirect evidence and simulations indicated that the diffuse continuum, especially the strong Balmer continuum from the broad line region (BLR), may contribute to the continuum in the u/U band. Here, we present direct evidence for this contribution. In this work, we apply the I… ▽ More The continuum reverberation mapping is widely used in studying accretion disk of active galactic nuclei (AGN). While some indirect evidence and simulations indicated that the diffuse continuum, especially the strong Balmer continuum from the broad line region (BLR), may contribute to the continuum in the u/U band. Here, we present direct evidence for this contribution. In this work, we apply the ICCF-Cut method to continuum reverberation mapping to extract the possible diffuse continuum light curves of 6 AGNs with high cadence, high quality and multi-band observations. We find the existence of an outer component out of the accretion disk for each of 6 AGNs in the Swift U band. Meanwhile, similar results can be derived by JAVELIN Photometric Reverberation Mapping Model for 4 of them. The lags of the outer components are consistent with the predicted Balmer continuum lags, which are about half of the H$β$ lag values. Our result directly reinforces that an outer component, especially the Balmer continuum in the rest-frame u/U band, can contribute significantly to the continuum reverberation lags of AGNs. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 17 pages, 10 figures, accepted for publication in the Astrophysical Journal

H$α$ Time Delays of AGNs from the Zwicky Transcient Facility Broadband Photometry

Authors: Qinchun Ma, Yuhan Wen, Xue-Bing Wu, Huapeng Gu, Yuming Fu

Abstract: In our previous work on broadband photometric reverberation mapping (PRM), we proposed the ICCF-Cut process to obtain the time lags of H$α$ emission line from two broadband lightcurves via subtracting the continuum emission from the line band. Extending the work, we enlarge our sample to the Zwicky Transient Facility (ZTF) database. We adopt two criteria to select 123 type 1 AGNs with sufficient v… ▽ More In our previous work on broadband photometric reverberation mapping (PRM), we proposed the ICCF-Cut process to obtain the time lags of H$α$ emission line from two broadband lightcurves via subtracting the continuum emission from the line band. Extending the work, we enlarge our sample to the Zwicky Transient Facility (ZTF) database. We adopt two criteria to select 123 type 1 AGNs with sufficient variability and smooth lightcurves from 3537 AGNs at $z<0.09$ with more than 100 epoch observations in the $g$ and $r$ bands from the ZTF database. We calculate the H$α$ time lags for 23 of them which have previous spectroscopic reverberation mapping (SRM) results using ICCF-Cut, JAVELIN and $χ^2$ methods. Our obtained H$α$ time lags are slightly larger than the H$β$ time lags, which is consistent with the previous SRM results and the theoretical model of the AGN broad line region. The comparisons between SRM and PRM lag distributions and between the subtracted emission line lightcurves indicate that after selecting AGNs with the two criteria, combining the ICCF-Cut, JAVELIN and $χ^2$ methods provides an efficient way to get the reliable H$α$ lags from the broadband PRM. Such techniques can be used to estimate the black hole masses of a large sample of AGNs in the large multi-epoch photometric sky surveys such as the Legacy Survey of Space and Time (LSST) and the survey from the Wide Field Survey Telescope (WFST) in the near future. △ Less

Submitted 15 March, 2024; originally announced March 2024.

Comments: 23 pages, 34 figures, accepted for publication in the Astrophysical Journal

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)

A comprehensive correction of the Gaia DR3 XP spectra

Authors: Bowen Huang, Haibo Yuan, Maosheng Xiang, Yang Huang, Kai Xiao, Shuai Xu, Ruoyi Zhang, Lin Yang, Zexi Niu, Hongrui Gu

Abstract: By combining spectra from the CALSPEC and NGSL, as well as spectroscopic data from the LAMOST Data Release 7 (DR7), we have analyzed and corrected the systematic errors of the Gaia DR3 BP/RP (XP) spectra. The errors depend on the normalized spectral energy distribution (simplified by two independent ``colors'') and $G$ magnitude. Our corrections are applicable in the range of approximately… ▽ More By combining spectra from the CALSPEC and NGSL, as well as spectroscopic data from the LAMOST Data Release 7 (DR7), we have analyzed and corrected the systematic errors of the Gaia DR3 BP/RP (XP) spectra. The errors depend on the normalized spectral energy distribution (simplified by two independent ``colors'') and $G$ magnitude. Our corrections are applicable in the range of approximately $-0.5<BP-RP<2$, $3<G<17.5$ and $E(B-V)<0.8$. To validate our correction, we conduct independent tests by comparing with the MILES and LEMONY spectra. The results demonstrate that the systematic errors of $BP-RP$ and $G$ have been effectively corrected, especially in the near ultraviolet. The consistency between the corrected Gaia XP spectra and the MILES and LEMONY is better than 2 per cent in the wavelength range of $336-400$\,nm and 1 per cent in redder wavelengths. A global absolute calibration is also carried out by comparing the synthetic Gaia photometry from the corrected XP spectra with the corrected Gaia DR3 photometry. Our study opens up new possibilities for using XP spectra in many fields. A Python package is publicly available to do the corrections (https://doi.org/10.12149/101375 or https://github.com/HiromonGON/GaiaXPcorrection). △ Less

Submitted 22 January, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: 20 pages, 17 figures, ApJS in Press

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

Constraints on charged black holes from merger-ringdown signals in GWTC-3 and prospects for the Einstein Telescope

Authors: Hua-Peng Gu, Hai-Tian Wang, Lijing Shao

Abstract: Whether astrophysical black holes (BHs) can have charge is a question to be addressed by observations. In the era of gravitational wave (GW) astronomy, one can constrain the charge of a merged BH remnant using the merger-ringdown signal of the GW data. Extending earlier studies, we analyze five GW events in GWTC-3, assuming Kerr-Newman BHs. Our results show no strong evidence for a charged BH, and… ▽ More Whether astrophysical black holes (BHs) can have charge is a question to be addressed by observations. In the era of gravitational wave (GW) astronomy, one can constrain the charge of a merged BH remnant using the merger-ringdown signal of the GW data. Extending earlier studies, we analyze five GW events in GWTC-3, assuming Kerr-Newman BHs. Our results show no strong evidence for a charged BH, and give a limit on the charge-to-mass-ratio $Q<0.37$ at $90\%$ credible level (CL). Due to the charge-spin degeneracy in the waveform and the limited signal-to-noise ratios (SNRs), it is challenging for LIGO/Virgo/KAGRA observations to provide better constraints. We further simulate data for the Einstein Telescope (ET), where SNRs can be as large as $\sim270$ in the ringdown signal. These simulated events allow us to consider the 220, 221, and 330 ringdown modes altogether, which can help break the charge-spin degeneracy. The analysis of a simulated GW150914-like signal shows that ET can improve the constraints on the charge-to-mass-ratio to $Q \lesssim 0.2$ at $90\%$ CL with one ringdown signal. △ Less

Submitted 10 January, 2024; v1 submitted 16 October, 2023; originally announced October 2023.

Comments: 13 pages, 7 figures; accepted by PRD

Journal ref: Phys. Rev. D 109 (2024) 024058

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

Beyond spectroscopy. II. Stellar parameters for over twenty million stars in the northern sky from SAGES DR1 and Gaia DR3

Authors: Yang Huang, Timothy C. Beers, Hai-Bo Yuan, Ke-Feng Tan, Wei Wang, Jie Zheng, Chun Li, Young Sun Lee, Hai-Ning Li, Jing-Kun Zhao, Xiang-Xiang Xue, Yu-Juan Liu, Hua-Wei Zhang, Xue-Ang Sun, Ji Li, Hong-Rui Gu, Christian Wolf, Christopher A. Onken, Ji-Feng Liu, Zhou Fan, Gang Zhao

Abstract: We present precise photometric estimates of stellar parameters, including effective temperature, metallicity, luminosity classification, distance, and stellar age, for nearly 26 million stars using the methodology developed in the first paper of this series, based on the stellar colors from the Stellar Abundances and Galactic Evolution Survey (SAGES) DR1 and Gaia EDR3. The optimal design of stella… ▽ More We present precise photometric estimates of stellar parameters, including effective temperature, metallicity, luminosity classification, distance, and stellar age, for nearly 26 million stars using the methodology developed in the first paper of this series, based on the stellar colors from the Stellar Abundances and Galactic Evolution Survey (SAGES) DR1 and Gaia EDR3. The optimal design of stellar-parameter sensitive $uv$ filters by SAGES has enabled us to determine photometric-metallicity estimates down to $-3.5$, similar to our previous results with the SkyMapper Southern Survey (SMSS), yielding a large sample of over five million metal-poor (MP; [Fe/H]$\le -1.0$) stars and nearly one million very metal-poor (VMP; [Fe/H]$\le -2.0$) stars. The typical precision is around $0.1$ dex for both dwarf and giant stars with [Fe/H]$>-1.0$, and 0.15-0.25/0.3-0.4 dex for dwarf/giant stars with [Fe/H]$<-1.0$. Using the precise parallax measurements and stellar colors from Gaia, effective temperature, luminosity classification, distance and stellar age are further derived for our sample stars. This huge data set in the Northern sky from SAGES, together with similar data in the Southern sky from SMSS, will greatly advance our understanding of the Milky Way, in particular its formation and evolution. △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: 14 pages, 14 figures, 3 tables, accepted by ApJ. arXiv admin note: text overlap with arXiv:2104.14154

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

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)

The H$α$ broadband photometric reverberation mapping of four Seyfert 1 galaxies

Authors: Qinchun Ma, Xue-Bing Wu, Huapeng Gu, Yuhan Wen, Yuming Fu

Abstract: Broadband photometric reverberation mapping (PRM) have been investigated for AGNs in recent years, but mostly on accretion disk continuum RM. Due to the small fraction of broad emission lines in the broadband, PRM for emission lines is very challenging. Here we present an ICCF-Cut method for broadband PRM to obtain the H$α$ broad line lag and apply it to four Seyfert 1 galaxies, MCG+08-11-011, NGC… ▽ More Broadband photometric reverberation mapping (PRM) have been investigated for AGNs in recent years, but mostly on accretion disk continuum RM. Due to the small fraction of broad emission lines in the broadband, PRM for emission lines is very challenging. Here we present an ICCF-Cut method for broadband PRM to obtain the H$α$ broad line lag and apply it to four Seyfert 1 galaxies, MCG+08-11-011, NGC 2617, 3C 120 and NGC 5548. All of them have high quality broadband lightcurves with daily/sub-daily cadence, which enables us to extract H$α$ lightcurves from the line band by subtracting the contributions from the continuum and host galaxy. Their extracted H$α$ lightcurves are compared with the lagged continuum band lightcurves, as well as the lagged H$β$ lightcurves obtained by spectroscopic RM (SRM) at the same epochs. The consistency of these lightcurves and the comparison with the SRM H$β$ lags provide supports to the H$α$ lags of these AGNs, in a range from 9 to 19 days, obtained by the ICCF-Cut, JAVELIN and $χ^2$ methods. The simulations to evaluate the reliability of H$α$ lags and the comparisons between SRM H$β$ and PRM H$α$ lags indicate that the consistency of the ICCF-Cut, JAVELIN and $χ^2$ results can ensure the reliability of the derived H$α$ lags. These methods may be used to estimate the broad line region sizes and black hole masses of a large sample of AGNs in the large multi-epoch high cadence photometric surveys such as LSST in the future. △ Less

Submitted 18 March, 2023; originally announced March 2023.

Comments: 22 pages, 19 figures, accepted for publication in ApJ

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)

Peta-electron volt gamma-ray emission from the Crab Nebula

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

Abstract: The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ pet… ▽ More The Crab pulsar and the surrounding nebula powered by the pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind is a bright source of gamma-rays carrying crucial information about this complex conglomerate. We report the detection of $γ$-rays with a spectrum showing gradual steepening over three energy decades, from $5\times 10^{-4}$ to $1.1$ petaelectronvolt (PeV). The ultra-high-energy photons exhibit the presence of a PeV electron accelerator (a pevatron) with an acceleration rate exceeding 15% of the absolute theoretical limit. Assuming that unpulsed $γ$-rays are produced at the termination of the pulsar's wind, we constrain the pevatron's size, between $0.025$ and $0.1$ pc, and the magnetic field $\approx 110 μ$G. The production rate of PeV electrons, $2.5 \times 10^{36}$ erg $\rm s^{-1}$, constitutes 0.5% of the pulsar's spin-down luminosity, although we do not exclude a non-negligible contribution of PeV protons to the production of the highest energy $γ$-rays. △ Less

Submitted 11 November, 2021; originally announced November 2021.

Comments: 43 pages, 13 figures, 2 tables; Published in Science

Journal ref: Science, 2021, Vol 373, Issue 6553, pp. 425-430

Calibration of the Air Shower Energy Scale of the Water and Air Cherenkov Techniques in the LHAASO experiment

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (233 additional authors not shown)

Abstract: The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galact… ▽ More The Wide Field-of-View Cherenkov Telescope Array (WFCTA) and the Water Cherenkov Detector Arrays (WCDA) of LHAASO are designed to work in combination for measuring the energy spectra of various cosmic ray species over a very wide energy range from a few TeV to 10 PeV. The energy calibration of WCDA can be achieved with a proven technique of measuring the westward shift of the Moon shadow of galactic cosmic rays due to the geomagnetic field. This deflection angle $Δ$ is inversely proportional to the energy of the cosmic rays. The precise measurements of the shifts by WCDA allows us to calibrate its energy scale for energies as high as 35 TeV. The energy scale measured by WCDA can be used to cross calibrate the energy reconstructed by WFCTA, which spans the whole energy range up to 10 PeV. In this work, we will demonstrate the feasibility of the method using the data collected from April 2019 to January 2020 by the WFCTA array and WCDA-1 detector, the first of the three water Cherenkov ponds, already commissioned at LHAASO site. △ Less

Submitted 13 April, 2021; v1 submitted 11 April, 2021; originally announced April 2021.

Journal ref: Phys. Rev. D 104, 062007 (2021)

Construction and On-site Performance of the LHAASO WFCTA Camera

Authors: F. Aharonian, Q. An, Axikegu, L. X. Bai, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, H. Cai, J. T. Cai, Z. Cao, Z. Cao, J. Chang, J. F. Chang, X. C. Chang, B. M. Chen, J. Chen, L. Chen, L. Chen, L. Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (234 additional authors not shown)

Abstract: The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this… ▽ More The focal plane camera is the core component of the Wide Field-of-view Cherenkov/fluorescence Telescope Array (WFCTA) of the Large High-Altitude Air Shower Observatory (LHAASO). Because of the capability of working under moonlight without aging, silicon photomultipliers (SiPM) have been proven to be not only an alternative but also an improvement to conventional photomultiplier tubes (PMT) in this application. Eighteen SiPM-based cameras with square light funnels have been built for WFCTA. The telescopes have collected more than 100 million cosmic ray events and preliminary results indicate that these cameras are capable of working under moonlight. The characteristics of the light funnels and SiPMs pose challenges (e.g. dynamic range, dark count rate, assembly techniques). In this paper, we present the design features, manufacturing techniques and performances of these cameras. Finally, the test facilities, the test methods and results of SiPMs in the cameras are reported here. △ Less

Submitted 4 July, 2021; v1 submitted 29 December, 2020; originally announced December 2020.

Comments: 45 pages, 21 figures, article

Journal ref: Eur. Phys. J. C 81, 657 (2021)