Showing 1–50 of 231 results for author: Cao, Z

The Velocity Aberration Effect of the CSST Main Survey Camera

Authors: Hui-Mei Feng, Zi-Huang Cao, Man I Lam, Ran Li, Hao Tian, Xin Zhang, Peng Wei, Xin-Feng Li, Wei Wang, Hugh R. A. Jones, Mao-Yuan Liu, Chao Liu

Abstract: In this study, we conducted simulations to find the geometric aberrations expected for images taken by the Main Survey Camera (MSC) of the Chinese Space Station Telescope (CSST) due to its motion. As anticipated by previous work, our findings indicate that the geometric distortion of light impacts the focal plane's apparent scale, with a more pronounced influence as the size of the focal plane inc… ▽ More In this study, we conducted simulations to find the geometric aberrations expected for images taken by the Main Survey Camera (MSC) of the Chinese Space Station Telescope (CSST) due to its motion. As anticipated by previous work, our findings indicate that the geometric distortion of light impacts the focal plane's apparent scale, with a more pronounced influence as the size of the focal plane increases. Our models suggest that the effect consistently influences the pixel scale in both the vertical and parallel directions. The apparent scale variation follows a sinusoidal distribution throughout one orbit period. Simulations reveal that the effect is particularly pronounced in the center of the Galaxy and gradually diminishes along the direction of ecliptic latitude. At low ecliptic latitudes, the total aberration leads to about 0.94 pixels offset (a 20-minute exposure) and 0.26 pixels offset (a 300-second exposure) at the edge of the field of view, respectively. Appropriate processings for the geometric effect during the CSST pre- and post-observation phases are presented. △ Less

Submitted 23 August, 2024; originally announced August 2024.

Comments: 13 pages, 8 figures; accepted by RAA

Interference detection in radio astronomy applying Shapiro-Wilks normality test, spectral entropy, and spectral relative entropy

Authors: Zhicheng Cao, Natalia A. Schmid, Kevin Bandura, Duncan R. Lorimer, Morgan Dameron, Katelyn Crockett, Clayton Grubick, Andreas Schmid, Shaonan Zheng

Abstract: Radio-frequency interference (RFI) is becoming an increasingly significant problem for most radio telescopes. Working with Green Bank Telescope data from PSR J1730+0747 in the form of complex-valued channelized voltages and their respective high-resolution power spectral densities, we evaluate a variety of statistical measures to characterize RFI. As a baseline for performance comparison, we use m… ▽ More Radio-frequency interference (RFI) is becoming an increasingly significant problem for most radio telescopes. Working with Green Bank Telescope data from PSR J1730+0747 in the form of complex-valued channelized voltages and their respective high-resolution power spectral densities, we evaluate a variety of statistical measures to characterize RFI. As a baseline for performance comparison, we use median absolute deviation (MAD) in complex channelized voltage data and spectral kurtosis (SK) in power spectral density data to characterize and filter out RFI. From a new perspective, we implement the Shapiro-Wilks (SW) test for normality and two information theoretical measures, spectral entropy (SE) and spectral relative entropy (SRE), and apply them to mitigate RFI. The baseline RFI mitigation algorithms are compared against our novel RFI detection algorithms to determine how effective and robust the performance is. Except for MAD, we find significant improvements in signal-to-noise ratio through the application of SE, symmetrical SRE, asymmetrical SRE, SK, and SW. These algorithms also do a good job of characterizing broadband RFI. Time- and frequency-variable RFI signals are best detected by SK and SW tests. △ Less

Submitted 12 August, 2024; originally announced August 2024.

Comments: 13 pages, 26 figures, accepted for publication in RAS Techniques and Instruments

Anisotropy of Nanohertz Gravitational Wave Background and Individual Sources from Supermassive Binary Black Holes: Probe of Cosmic Large Scale Structure

Authors: Qing Yang, Xiao Guo, Zhoujian Cao, Xiaoyun Shao, Xi Yuan

Abstract: Several pulsar timing array (PTA) groups have recently claimed the detection of nanohertz gravitational wave (GW) background, but the origin of this GW signal remains unclear. Nanohertz GWs generated by supermassive binary black holes (SMBBHs) are one of the most important GW sources in the PTA band. Utilizing data from numerical cosmology simulation, we generate mock SMBBHs within the observable… ▽ More Several pulsar timing array (PTA) groups have recently claimed the detection of nanohertz gravitational wave (GW) background, but the origin of this GW signal remains unclear. Nanohertz GWs generated by supermassive binary black holes (SMBBHs) are one of the most important GW sources in the PTA band. Utilizing data from numerical cosmology simulation, we generate mock SMBBHs within the observable universe and treat them as PTA band GW sources. We present their statistical properties, and analyze the isotropic and anisotropic characteristics of the gravitational wave background (GWB) signal they produce. Specifically, we derive the characteristic amplitude and spectrum of the GWB signal, and calculate the angular power spectrum for both GW strains/energy density and the position distribution of GW sources. We predict that the angular power spectrum of GWB energy density has $C_1/C_0\approx0.40\pm0.32$, and $C_l/C_0\simeq \frac{1}{2(2l+1)}$ (for $l>1$). Furthermore, for the upcoming Chinese Pulsar Timing Array (CPTA) and Square Kilometre Array (SKA) PTA, we predict the spatial distribution, numbers and signal-to-noise ratio (SNR) distribution of individual GW sources that may be detected with SNR>8, and study the anisotropy property in the spatial distribution of these individual GW sources. △ Less

Submitted 9 August, 2024; originally announced August 2024.

Comments: 26 pages,18 figures,1 table

Gravitational Wave Mixture Separation for Future Gravitational Wave Observatories Utilizing Deep Learning

Authors: Cunliang Ma, Weiguang Zhou, Zhoujian Cao

Abstract: Future GW observatories, such as the Einstein Telescope (ET), are expected to detect gravitational wave signals, some of which are likely to overlap with each other. This overlap may lead to misidentification as a single GW event, potentially biasing the estimated parameters of mixture GWs. In this paper, we adapt the concept of speech separation to address this issue by applying it to signal sepa… ▽ More Future GW observatories, such as the Einstein Telescope (ET), are expected to detect gravitational wave signals, some of which are likely to overlap with each other. This overlap may lead to misidentification as a single GW event, potentially biasing the estimated parameters of mixture GWs. In this paper, we adapt the concept of speech separation to address this issue by applying it to signal separation of overlapping GWs. We show that deep learning models can effectively separate overlapping GW signals. The proposed method may aid in eliminating biases in parameter estimation for such signals. △ Less

Submitted 18 July, 2024; originally announced July 2024.

Comments: 15 pages, 8 figures

Swift-BAT GUANO follow-up of gravitational-wave triggers in the third LIGO-Virgo-KAGRA observing run

Authors: Gayathri Raman, Samuele Ronchini, James Delaunay, Aaron Tohuvavohu, Jamie A. Kennea, Tyler Parsotan, Elena Ambrosi, Maria Grazia Bernardini, Sergio Campana, Giancarlo Cusumano, Antonino D'Ai, Paolo D'Avanzo, Valerio D'Elia, Massimiliano De Pasquale, Simone Dichiara, Phil Evans, Dieter Hartmann, Paul Kuin, Andrea Melandri, Paul O'Brien, Julian P. Osborne, Kim Page, David M. Palmer, Boris Sbarufatti, Gianpiero Tagliaferri , et al. (1797 additional authors not shown)

Abstract: We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wav… ▽ More We present results from a search for X-ray/gamma-ray counterparts of gravitational-wave (GW) candidates from the third observing run (O3) of the LIGO-Virgo-KAGRA (LVK) network using the Swift Burst Alert Telescope (Swift-BAT). The search includes 636 GW candidates received in low latency, 86 of which have been confirmed by the offline analysis and included in the third cumulative Gravitational-Wave Transient Catalogs (GWTC-3). Targeted searches were carried out on the entire GW sample using the maximum--likelihood NITRATES pipeline on the BAT data made available via the GUANO infrastructure. We do not detect any significant electromagnetic emission that is temporally and spatially coincident with any of the GW candidates. We report flux upper limits in the 15-350 keV band as a function of sky position for all the catalog candidates. For GW candidates where the Swift-BAT false alarm rate is less than 10$^{-3}$ Hz, we compute the GW--BAT joint false alarm rate. Finally, the derived Swift-BAT upper limits are used to infer constraints on the putative electromagnetic emission associated with binary black hole mergers. △ Less

Submitted 13 July, 2024; originally announced July 2024.

Comments: 50 pages, 10 figures, 4 tables

Extraction of binary neutron star gravitational wave waveforms from Einstein Telescope using deep learning

Authors: Cunliang Ma, Xinyao Yu, Zhoujian Cao, Mingzhen Jia

Abstract: In the future, the third generation (3G) gravitational wave (GW) detectors, exemplified by the Einstein Telescope (ET), will be operational. The detection rate of GW from binary neutron star (BNS) is expected to reach approximately $10^4$ per year. To address the challenges posed by BNS GW data processing for 3G GW detectors, this paper explores the extraction of BNS waveforms from ET. Drawing ins… ▽ More In the future, the third generation (3G) gravitational wave (GW) detectors, exemplified by the Einstein Telescope (ET), will be operational. The detection rate of GW from binary neutron star (BNS) is expected to reach approximately $10^4$ per year. To address the challenges posed by BNS GW data processing for 3G GW detectors, this paper explores the extraction of BNS waveforms from ET. Drawing inspiration from SPIIR's matched filtering approach, we introduce a novel framework leveraging deep learning for BNS waveform extraction. By integrating denoised outputs of time-delayed strain, we can reconstruct the embedded BNS waveform. We have established three distinct BNS GW denoising models, each tailored to address the early inspiral, later inspiral, and merger phases of BNS GW, respectively. To further regulate the waveform shape, we propose the Amplitude Regularity Model that takes denoised output as input and regulated waveform as output. The experiments conducted on test data demonstrate the efficacy of the denoising models, the Amplitude Regularity Models, as well as the overall waveform construction method. To the best of our knowledge, this marks the first instance of deep learning being utilized for the task of BNS waveform extraction. We believe that the proposed method holds promise for early warning, searching, and localization of BNS GWs. △ Less

Submitted 22 June, 2024; originally announced June 2024.

Comments: 14 pages, 16 figures

The projected sensitivity of SCEP experiment to Magnetic Monopole

Authors: Changqing Ye, Beige Liu, Zhe Cao, Lingzhi Han, Xinming Huang, Min Jiang, Dong Liu, Qing Lin, Shitian Wan, Yusheng Wu, Lei Zhao, Yue Zhang, Xinhua Peng, Zhengguo Zhao

Abstract: The investigation of beyond-Standard-Model particles is a compelling direction in the pursuit of new physics. One such hypothetical particle, the magnetic monopole, has garnered considerable attention due to its strong theoretical motivation and potential to unveil profound physical phenomena. The magnetic monopole is intricately linked to the long-standing enigma surrounding the quantization of e… ▽ More The investigation of beyond-Standard-Model particles is a compelling direction in the pursuit of new physics. One such hypothetical particle, the magnetic monopole, has garnered considerable attention due to its strong theoretical motivation and potential to unveil profound physical phenomena. The magnetic monopole is intricately linked to the long-standing enigma surrounding the quantization of electric charge. In this manuscript, we propose a novel detection scenario for magnetic monopoles by employing a coincidence measurement technique that combines a room-temperature magnetometer with plastic scintillators. This setup allows for the collection of both the induction and scintillation signals generated by the passage of a monopole. The estimation of the sensitivity using a simple benchmark setup is given. △ Less

Submitted 19 June, 2024; v1 submitted 18 June, 2024; originally announced June 2024.

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

All-sky Guide Star Catalog for CSST

Authors: Hui-Mei Feng, Zi-Huang Cao, Man I Lam, Ran Li, Hao Tian, Da-Yi Yin, Yuan-Yu Yang, Xin Zhang, Dong-Wei Fan, Yi-Qiao Dong, Xin-Feng Li, Wei Wang, Long Li, Hugh R. A. Jones, Yi-Han Tao, Jia-Lu Nie, Pei-Pei Wang, Mao-Yuan Liu, He-jun Yang, Chao Liu

Abstract: The China Space Station Telescope (CSST) is a two-meter space telescope with multiple back-end instruments. The Fine Guidance Sensor (FGS) is an essential subsystem of the CSST Precision Image Stability System to ensure the required absolute pointing accuracy and line-of-sight stabilization. In this study, we construct the Main Guide Star Catalog for FGS. To accomplish this, we utilize the informa… ▽ More The China Space Station Telescope (CSST) is a two-meter space telescope with multiple back-end instruments. The Fine Guidance Sensor (FGS) is an essential subsystem of the CSST Precision Image Stability System to ensure the required absolute pointing accuracy and line-of-sight stabilization. In this study, we construct the Main Guide Star Catalog for FGS. To accomplish this, we utilize the information about the FGS and object information from the Gaia Data Release 3. We provide an FGS instrument magnitude and exclude variables, binaries, and high proper motion stars from the catalog to ensure uniform FGS guidance capabilities. Subsequently, we generate a HEALPix index, which provides a hierarchical tessellation of the celestial sphere, and employ the Voronoi algorithm to achieve a homogeneous distribution of stars across the catalog. This distribution ensures adequate coverage and sampling of the sky. The performance of the CSST guide star catalog was assessed by simulating the field of view of the FGS according to the CSST mock survey strategy catalog. The analysis of the results indicates that this catalog provides adequate coverage and accuracy. The catalog's performance meets the FGS requirements, ensuring the functioning of the FGS and its guidance capabilities. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: published on RAA

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

Tidal disruption event AT2020ocn: early-time X-ray flares caused by a possible disc alignment process

Authors: Z. Cao, P. G. Jonker, D. R. Pasham, S. Wen, N. C. Stone, A. I. Zabludoff

Abstract: A tidal disruption event (TDE) may occur when a star is torn apart by the tidal force of a black hole (BH). Eventually, an accretion disc is thought to form out of stellar debris falling back towards the BH. If the star's orbital angular momentum vector prior to disruption is not aligned with the BH spin angular momentum vector, the disc will be tilted with respect to the BH equatorial plane. The… ▽ More A tidal disruption event (TDE) may occur when a star is torn apart by the tidal force of a black hole (BH). Eventually, an accretion disc is thought to form out of stellar debris falling back towards the BH. If the star's orbital angular momentum vector prior to disruption is not aligned with the BH spin angular momentum vector, the disc will be tilted with respect to the BH equatorial plane. The disc will eventually be drawn into the BH equatorial plane due to a combination of the Bardeen-Petterson effect and internal torques. Here, we analyse the X-ray and UV observations of the TDE AT2020ocn obtained by Swift, XMM-Newton, and NICER. The X-ray light curve shows strong flares during the first $\approx100$ days, while, over the same period, the UV emission decays gradually. We find that the X-ray flares can be explained by a model that also explains the spectral evolution. This model includes a slim disc viewed under a variable inclination plus an inverse-Comptonisation component processing the slim disc emission. A scenario where the ongoing Lense-Thirring precession during the disc alignment process is responsible for the observed inclination variations is consistent with the data. In later observations, we find that the X-ray spectrum of AT2020ocn becomes harder, while the mass accretion rate remains at super-Eddington levels, suggesting the formation of a corona in line with accretion onto other compact objects. We constrain the BH mass to be $(7^{+13}_{-3})\times10^{5}$ M$_\odot$ at the 1$σ$ (68%) confidence level. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 14 pages, 8 figures, 5 tables, with supplementary materials. Accepted for publication in ApJ

AT2018fyk: Candidate Tidal Disruption Event by a (Super)massive Black Hole Binary

Authors: S. Wen, P. G. Jonker, A. J. Levan, D. Li, N. C. Stone, A. I. Zabludoff, Z. Cao, T. Wevers, D. R. Pasham, C. Lewin, E. Kara

Abstract: The tidal disruption event (TDE) AT2018fyk has unusual X-ray, UV, and optical light curves that decay over the first $\sim$600d, rebrighten, and decay again around 1200d. We explain this behavior as a one-off TDE associated with a massive black hole (BH) \emph{binary}. The sharp drop-offs from $t^{-5/3}$ power laws at around 600d naturally arise when one BH interrupts the debris fallback onto the… ▽ More The tidal disruption event (TDE) AT2018fyk has unusual X-ray, UV, and optical light curves that decay over the first $\sim$600d, rebrighten, and decay again around 1200d. We explain this behavior as a one-off TDE associated with a massive black hole (BH) \emph{binary}. The sharp drop-offs from $t^{-5/3}$ power laws at around 600d naturally arise when one BH interrupts the debris fallback onto the other BH. The BH mass $M_\bullet$ derived from fitting X-ray spectra with a slim disk accretion model and, independently, from fitting the early UV/optical light curves, is smaller by two orders of magnitude than predicted from the $M_\bullet$--$σ_*$ host galaxy relation, suggesting that the debris is accreted onto the secondary, with fallback cut off by the primary. Furthermore, if the rebrightening were associated with the primary, it should occur around 5000d, not the observed 1200d. The secondary's mass and dimensionless spin is $M_{\bullet,{\rm s}}=2.7^{+0.5}_{-1.5} \times 10^5 M_\odot$ and $a_{\bullet,{\rm s}}>0.3$ (X-ray spectral fitting), while the primary's mass is $M_{\bullet,{\rm p}}=10^{7.7\pm0.4}M_\odot$ ($M_\bullet$-$σ_*$ relation). An intermediate mass BH secondary is consistent with the observed UV/optical light curve decay, i.e., the secondary's outer accretion disk is too faint to produce a detectable emission floor. The time of the first accretion cutoff constrains the binary separation to be $(6.7\pm 1.2) \times 10^{-3}~{\rm pc}$. X-ray spectral fitting and timing analysis indicate that the hard X-rays arise from a corona above the secondary's disk. The early UV/optical emission, suggesting a super-Eddington phase for the secondary, possibly originates from shocks arising from debris circularization. △ Less

Submitted 29 July, 2024; v1 submitted 1 May, 2024; originally announced May 2024.

Comments: 5 tables, 7 figures

The hunt of PeVatrons as the origin of the most energetic photons observed in our Galaxy

Authors: Emma de Oña Wilhelmi, Rubén López-Coto, Felix Aharonian, Elena Amato, Zhen Cao, Stefano Gabici, Jim Hinton

Abstract: Ultrarelativistic particles called cosmic rays permeate the Milky Way, propagating through the Galactic turbulent magnetic fields. The mechanisms under which these particles increase their energy can be reasonably described by current theories of acceleration and propagation of cosmic rays. There are, however, still many open questions as to how to reach petaelectronvolt (PeV) energies, the maximu… ▽ More Ultrarelativistic particles called cosmic rays permeate the Milky Way, propagating through the Galactic turbulent magnetic fields. The mechanisms under which these particles increase their energy can be reasonably described by current theories of acceleration and propagation of cosmic rays. There are, however, still many open questions as to how to reach petaelectronvolt (PeV) energies, the maximum energy believed to be attained in our Galaxy, and in which astrophysical sources (dubbed {\it PeVatrons}) this ultra-high energy acceleration happens. In this article, we describe the theoretical conditions for plasma acceleration to these energies, and the Galactic sources in which these conditions are possible. These theoretical predictions are then confronted with the latest experimental results, summarising the state-of-the-art of our current knowledge of PeVatrons. We finally describe the prospects to keep advancing the understanding of these elusive objects, still unidentified more than one hundred years after the discovery of cosmic rays. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: Published in Nature Astronomy, Volume 8, pages 425-431 (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.

Polarization Holes as an Indicator of Magnetic Field-Angular Momentum Alignment I. Initial Tests

Authors: Lijun Wang, Zhuo Cao, Xiaodan Fan, Hua-bai Li

Abstract: The formation of protostellar disks is still a mystery, largely due to the difficulties in observations that can constrain theories. For example, the 3D alignment between the rotation of the disk and the magnetic fields (B-fields) in the formation environment is critical in some models, but so far impossible to observe. Here, we study the possibility of probing the alignment between B-field and di… ▽ More The formation of protostellar disks is still a mystery, largely due to the difficulties in observations that can constrain theories. For example, the 3D alignment between the rotation of the disk and the magnetic fields (B-fields) in the formation environment is critical in some models, but so far impossible to observe. Here, we study the possibility of probing the alignment between B-field and disk rotation using ``polarization holes'' (PHs). PHs are widely observed and are caused by unresolved B-field structures. With ideal magnetohydrodynamic (MHD) simulations, we demonstrate that different initial alignments between B-field and angular momentum (AM) can result in B-field structures that are distinct enough to produce distinguishable PHs. Thus PHs can potentially serve as probes for alignments between B-field and AM in disk formation. △ Less

Submitted 22 March, 2024; originally announced March 2024.

Comments: accepted by 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)

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

Prospects for joint reconstruction of imaging air Cherenkov Telescope array and extensive air shower array

Authors: Zhipeng Zhang, Ruizhi Yang, Shoushan Zhang, Liqiao Yin, Jiali Liu, Yudong Wang, Lingling Ma, Zhen Cao

Abstract: In this paper we proposed a joint reconstruction of \gray events using both extensive air array (EAS) and Imaging air Cherenkov Telescope array (IACT). We considered eight Cherenkov telescopes to be built on the LHAASO (Large High Altitude Air Shower Observatory) site and investigate the improvement in differential sensitivity when combining the information from both IACT and Moun detectors of LHA… ▽ More In this paper we proposed a joint reconstruction of \gray events using both extensive air array (EAS) and Imaging air Cherenkov Telescope array (IACT). We considered eight Cherenkov telescopes to be built on the LHAASO (Large High Altitude Air Shower Observatory) site and investigate the improvement in differential sensitivity when combining the information from both IACT and Moun detectors of LHAASO-KM2A. We found that due to the higher cosmic ray background rejection power and higher gamma ray retention ratio provided by muon detectors of LHAASO, such a joint reconstruction can significantly improve the sensitivity of IACTs, especially for extended sources and long exposure time. In this article, we have shown the performance of an eight-telescopes mini array, and our results indicate that above $10~\rm TeV$, the sensitivity can be improved by muon detector from $25\% - 60\%$ in different energy ranges. △ Less

Submitted 28 June, 2024; v1 submitted 17 February, 2024; originally announced February 2024.

Comments: 6 pages, 6 figures

GWAI: Harnessing Artificial Intelligence for Enhancing Gravitational Wave Data Analysis

Authors: Tianyu Zhao, Yue Zhou, Ruijun Shi, Zhoujian Cao, Zhixiang Ren

Abstract: Gravitational wave (GW) astronomy has opened new frontiers in understanding the cosmos, while the integration of artificial intelligence (AI) in science promises to revolutionize data analysis methodologies. However, a significant gap exists, as there is currently no dedicated platform that enables scientists to develop, test, and evaluate AI algorithms efficiently. To address this gap, we introdu… ▽ More Gravitational wave (GW) astronomy has opened new frontiers in understanding the cosmos, while the integration of artificial intelligence (AI) in science promises to revolutionize data analysis methodologies. However, a significant gap exists, as there is currently no dedicated platform that enables scientists to develop, test, and evaluate AI algorithms efficiently. To address this gap, we introduce GWAI, a pioneering AI-centered software platform designed for gravitational wave data analysis. GWAI contains a three-layered architecture that emphasizes simplicity, modularity, and flexibility, covering the entire analysis pipeline. GWAI aims to accelerate scientific discoveries, bridging the gap between advanced AI techniques and astrophysical research. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 10 pages, 5 figures

Accuracy of numerical relativity waveforms with respect to space-based gravitational wave detectors

Authors: Zun Wang, Junjie Zhao, Zhoujian Cao

Abstract: As with the laser interferometer gravitational-wave observatory (LIGO), the matched filtering technique will be critical to the data analysis of gravitational wave detection by space-based detectors, including LISA, Taiji and Tianqin. Waveform templates are the basis for such matched filtering techniques. To construct ready-to-use waveform templates, numerical relativity waveforms are a starting p… ▽ More As with the laser interferometer gravitational-wave observatory (LIGO), the matched filtering technique will be critical to the data analysis of gravitational wave detection by space-based detectors, including LISA, Taiji and Tianqin. Waveform templates are the basis for such matched filtering techniques. To construct ready-to-use waveform templates, numerical relativity waveforms are a starting point. Therefore, the accuracy issue of numerical relativity waveforms is critically important. There are many investigations regarding this issue with respect to LIGO. But unfortunately there are few results on this issue with respect to space-based detectors. The current paper investigates this problem. Our results indicate that the existing numerical relativity waveforms are as accurate as 99% with respect to space-based detectors, including LISA, Taiji and Tianqin. Such an accuracy level is comparable to that with respect to LIGO. △ Less

Submitted 27 January, 2024; originally announced January 2024.

Simultaneously search for multi-target Galactic binary gravitational waves in reduced parameter space with LMPSO-CV

Authors: Pin Gao, Xilong Fan, Zhoujian Cao

Abstract: We propose an innovative approach to the concurrent exploration of gravitational waves originating from Galactic binaries through the development of a new Local Maxima Particle Swarm Optimization (LMPSO) algorithm. Our methodology employs strategic Create Voids (CV) to streamline parameter space, maximizing the identification of local maxima for the $\mathcal{F}$-statistic even in the overlapped s… ▽ More We propose an innovative approach to the concurrent exploration of gravitational waves originating from Galactic binaries through the development of a new Local Maxima Particle Swarm Optimization (LMPSO) algorithm. Our methodology employs strategic Create Voids (CV) to streamline parameter space, maximizing the identification of local maxima for the $\mathcal{F}$-statistic even in the overlapped signals case. Subsequently, a "find-real-$\mathcal{F}$-statistic-analysis", which implements the astrophysical models and properties of $\mathcal{F}$-statistic in parameter space, is conducted to reveal Galactic binary gravitational wave signals within the dataset. Our new approach eliminates inaccuracies associated with signal subtraction contamination, which is a challenge for traditional iterative-subtraction method when addressing low signal-to-noise ratio signals (e.g., SNR $<$ 15). We also consider the effect of overlapping signals and degeneracy noise. To demonstrate the efficacy of our approach, we utilize the residuals from the LISA mock data challenge (LDC1-4), where 10982 injection sources with SNR $\ge$ 15 have been eliminated. For the remaining low SNR sources (SNR $<$ 15), the LMPSO-CV method efficiently identifies 8995 signals with a 52.3\% detection rate or 3463 signals with a 73.1\% detection rate when the correlation coefficient threshold $R_t$ is set to 0.8. △ Less

Submitted 10 May, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Testing an exact diffraction formula with gravitational wave source lensed by a supermassive black hole in binary systems

Authors: Xiao Guo, Zhoujian Cao

Abstract: When it comes to long-wavelength gravitational waves (GWs), diffraction effect becomes significant when these waves are lensed by celestial bodies. Typically, the traditional diffraction integral formula neglects large-angle diffraction, which is often adequate for most of cases. Nonetheless, there are specific scenarios, such as when a GW source is lensed by a supermassive black hole in a binary… ▽ More When it comes to long-wavelength gravitational waves (GWs), diffraction effect becomes significant when these waves are lensed by celestial bodies. Typically, the traditional diffraction integral formula neglects large-angle diffraction, which is often adequate for most of cases. Nonetheless, there are specific scenarios, such as when a GW source is lensed by a supermassive black hole in a binary system, where the lens and source are in close proximity, where large-angle diffraction can play a crucial role. In our prior research, we have introduced an exact, general diffraction integral formula that accounts for large-angle diffraction as well. This paper explores the disparities between this exact diffraction formula and the traditional, approximate one under various special conditions. Our findings indicate that, under specific parameters, such as a lens-source distance of $D_{\rm LS}=0.1$\,AU and a lens mass of $M_{\rm L}=4\times10^{6}M_\odot$, the amplification factor for the exact diffraction formula is notably smaller than that of the approximate formula, differing by a factor of approximately $r_F\simeq0.806$. This difference is substantial enough to be detectable. Furthermore, our study reveals that the proportionality factor $r_F$ gradually increases from 0.5 to 1 as $D_{\rm LS}$ increases, and decreases as $M_{\rm L}$ increases. Significant differences between the exact and approximate formulas are observable when $D_{\rm LS}\lesssim0.2$ AU (assuming $M_{\rm L}=4\times10^{6}M_\odot$) or when $M_{\rm L}\gtrsim2\times10^6M_\odot$ (assuming $D_{\rm LS}=0.1$ AU). These findings suggest that there is potential to validate our general diffraction formula through future GW detections. △ Less

Submitted 15 April, 2024; v1 submitted 3 January, 2024; originally announced January 2024.

Comments: 17 pages, 7 figures, comments are welcome

Satellite Impact on Astronomical Observations Based on Elliptical Orbit Model

Authors: Tianzhu Hu, Yong Zhang, Xiangqun Cui, Zihuang Cao, Kang Huang, Jingyi Cai, Jun Li, Tong Zhou

Abstract: Space-based and ground-based telescopes have extensively documented the impact of satellites on astronomical observations. With the proliferation of satellite mega-constellation programs, their influence on astronomical observations has become undeniable. It is crucial to quantify the impact of satellites on telescopes. To address this need, we have enhanced the circular orbit model for satellites… ▽ More Space-based and ground-based telescopes have extensively documented the impact of satellites on astronomical observations. With the proliferation of satellite mega-constellation programs, their influence on astronomical observations has become undeniable. It is crucial to quantify the impact of satellites on telescopes. To address this need, we have enhanced the circular orbit model for satellites and introduced a methodology based on two-line element (TLE) orbit data. This involves constructing a satellite probability distribution model to evaluate the impact of satellites on telescopes. Using our method, we assessed the satellite impact on global observatories. The results indicate that the regions most severely affected by satellite interference currently are those near the equator, with latitudes around 50 and 80 degrees experiencing the most significant impact from low Earth orbit satellites. Furthermore, we validated the reliability of our method using imaging data obtained from the focal surface acquisition camera of the LAMOST telescope. △ Less

Submitted 21 December, 2023; originally announced December 2023.

Comments: 9 pages, 7 figures, submitted to A&A, any comments are welcome

Journal ref: A&A 687, A135 (2024)

Scientific Preparation for CSST: Classification of Galaxy and Nebula/Star Cluster Based on Deep Learning

Authors: Yuquan Zhang, Zhong Cao, Feng Wang, Lam, Man I, Hui Deng, Ying Mei, Lei Tan

Abstract: The Chinese Space Station Telescope (abbreviated as CSST) is a future advanced space telescope. Real-time identification of galaxy and nebula/star cluster (abbreviated as NSC) images is of great value during CSST survey. While recent research on celestial object recognition has progressed, the rapid and efficient identification of high-resolution local celestial images remains challenging. In this… ▽ More The Chinese Space Station Telescope (abbreviated as CSST) is a future advanced space telescope. Real-time identification of galaxy and nebula/star cluster (abbreviated as NSC) images is of great value during CSST survey. While recent research on celestial object recognition has progressed, the rapid and efficient identification of high-resolution local celestial images remains challenging. In this study, we conducted galaxy and NSC image classification research using deep learning methods based on data from the Hubble Space Telescope. We built a Local Celestial Image Dataset and designed a deep learning model named HR-CelestialNet for classifying images of the galaxy and NSC. HR-CelestialNet achieved an accuracy of 89.09% on the testing set, outperforming models such as AlexNet, VGGNet and ResNet, while demonstrating faster recognition speeds. Furthermore, we investigated the factors influencing CSST image quality and evaluated the generalization ability of HR-CelestialNet on the blurry image dataset, demonstrating its robustness to low image quality. The proposed method can enable real-time identification of celestial images during CSST survey mission. △ Less

Submitted 8 December, 2023; originally announced December 2023.

Dawning of a New Era in Gravitational Wave Data Analysis: Unveiling Cosmic Mysteries via Artificial Intelligence -- A Systematic Review

Authors: Tianyu Zhao, Ruijun Shi, Yue Zhou, Zhoujian Cao, Zhixiang Ren

Abstract: Background: Artificial intelligence (AI), with its vast capabilities, has become an integral part of our daily interactions, particularly with the rise of sophisticated models like Large Language Models. These advancements have not only transformed human-machine interactions but have also paved the way for significant breakthroughs in various scientific domains. Aim of review: This review is cente… ▽ More Background: Artificial intelligence (AI), with its vast capabilities, has become an integral part of our daily interactions, particularly with the rise of sophisticated models like Large Language Models. These advancements have not only transformed human-machine interactions but have also paved the way for significant breakthroughs in various scientific domains. Aim of review: This review is centered on elucidating the profound impact of AI, especially deep learning, in the field of gravitational wave data analysis (GWDA). We aim to highlight the challenges faced by traditional GWDA methodologies and how AI emerges as a beacon of hope, promising enhanced accuracy, real-time processing, and adaptability. Key scientific concepts of review: Gravitational wave (GW) waveform modeling stands as a cornerstone in the realm of GW research, serving as a sophisticated method to simulate and interpret the intricate patterns and signatures of these cosmic phenomena. This modeling provides a deep understanding of the astrophysical events that produce gravitational waves. Next in line is GW signal detection, a refined technique that meticulously combs through extensive datasets, distinguishing genuine gravitational wave signals from the cacophony of background noise. This detection process is pivotal in ensuring the authenticity of observed events. Complementing this is the GW parameter estimation, a method intricately designed to decode the detected signals, extracting crucial parameters that offer insights into the properties and origins of the waves. Lastly, the integration of AI for GW science has emerged as a transformative force. AI methodologies harness vast computational power and advanced algorithms to enhance the efficiency, accuracy, and adaptability of data analysis in GW research, heralding a new era of innovation and discovery in the field. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Compact Binary Systems Waveform Generation with Generative Pre-trained Transformer

Authors: Ruijun Shi, Yue Zhou, Tianyu Zhao, Zhoujian Cao, Zhixiang Ren

Abstract: Space-based gravitational wave (GW) detection is one of the most anticipated GW detection projects in the next decade, which promises to detect abundant compact binary systems. At present, deep learning methods have not been widely explored for GW waveform generation and extrapolation. To solve the data processing difficulty and the increasing waveform complexity caused by the detector's response… ▽ More Space-based gravitational wave (GW) detection is one of the most anticipated GW detection projects in the next decade, which promises to detect abundant compact binary systems. At present, deep learning methods have not been widely explored for GW waveform generation and extrapolation. To solve the data processing difficulty and the increasing waveform complexity caused by the detector's response and second-generation time-delay interferometry (TDI 2.0), an interpretable pre-trained large model named CBS-GPT (Compact Binary Systems Waveform Generation with Generative Pre-trained Transformer) is proposed. For compact binary system waveforms, three models were trained to predict the waveforms of massive black hole binaries (MBHB), extreme mass-ratio inspirals (EMRIs), and galactic binaries (GB), achieving prediction accuracies of at most 99%, 91%, and 99%, respectively. The CBS-GPT model exhibits notable generalization and interpretability, with its hidden parameters effectively capturing the intricate information of waveforms, even with the complex instrument response and a wide parameter range. Our research demonstrates the potential of large models in the GW realm, opening up new opportunities and guidance for future researches such as complex waveforms generation, gap completion, and deep learning model design for GW science. △ Less

Submitted 5 March, 2024; v1 submitted 31 October, 2023; originally announced October 2023.

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

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

Eccentricity enables the earliest warning and localization of gravitational waves with ground-based detectors

Authors: Tao Yang, Rong-Gen Cai, Zhoujian Cao, Hyung Mok Lee

Abstract: The early and precise localization of gravitational waves (GWs) is pivotal in detecting their electromagnetic (EM) counterparts, especially for binary neutron stars (BNS) and neutron star-black hole binaries (NSBH). In this letter, we pioneer the exploration of utilizing the higher harmonic modes induced by the eccentricity of compact binaries to localize GWs with ground-based detectors even befor… ▽ More The early and precise localization of gravitational waves (GWs) is pivotal in detecting their electromagnetic (EM) counterparts, especially for binary neutron stars (BNS) and neutron star-black hole binaries (NSBH). In this letter, we pioneer the exploration of utilizing the higher harmonic modes induced by the eccentricity of compact binaries to localize GWs with ground-based detectors even before the quadrupole baseline $\ell=2$ mode enters the detector band. Our theoretical analysis marks a first in proposing a strategy for gaining the earliest possible warning and maximizing preparation time for observing pre- and/or post-merger EM counterparts. We simulate three typical binaries from GWTC-3 with eccentricities ranging from 0.05 to 0.4. Our results reveal that the third-generation (3G) detectors (low frequency cut-off $f_0=5$ Hz) can accumulate sufficient signal-to-noise ratios through higher modes before the onset of the baseline $\ell=2$ mode entry into the band. Notably, relying solely on the higher modes, the 3G detector network ET+2CE achieves an average localization on the order of $1-10^2~\rm deg^2$ around 1-1.8 hours before the merger of a GW170817-like BNS, and $10-10^3~\rm deg^2$ approximately 18-30 minutes prior to the merger of a GW200115-like NSBH. A $100~\rm deg^2$ localization is attainable even 2-4 hours prior to a BNS merger. Moreover, in the near face-on orientations which are generally more favorable for EM counterpart detection, the localization can be further improved. △ Less

Submitted 13 May, 2024; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: 8 pages, 2 figures, published in PRD

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

Ultra-High-Energy Gamma-Ray Astronomy

Authors: Zhen Cao, Songzhan Chen, Ruoyu Liu, Ruizhi Yang

Abstract: Ultra-High Energy (UHE, $>$0.1\,PeV) $γ$-ray Astronomy is rapidly evolving into an expanding branch of the $γ$-ray astronomy with the surprising discovery of 12 PeVatrons and the detection of a handful of photons above 1 PeV. Nearly all known celestial object types that have emissions in the TeV band are found also emitting UHE photons. UHE $γ$-rays have a well-defined horizon inside our galaxy du… ▽ More Ultra-High Energy (UHE, $>$0.1\,PeV) $γ$-ray Astronomy is rapidly evolving into an expanding branch of the $γ$-ray astronomy with the surprising discovery of 12 PeVatrons and the detection of a handful of photons above 1 PeV. Nearly all known celestial object types that have emissions in the TeV band are found also emitting UHE photons. UHE $γ$-rays have a well-defined horizon inside our galaxy due to the absorption of infrared and cosmic microwave backgrounds in the universe. With the last 30 years, traditional cosmic ray (CR) detection techniques allow the detection of UHE $γ$-rays, and opened up the last observation window. For leptonic sources, UHE radiation is in the deep Klein-Nishina regime which is largely suppressed. Therefore UHE $γ$-ray detection will help to locate and identify hadronic radiation sources, tracing the historic pursuit for the origin of CRs around the knee of the spectrum. The Crab Nebula is again the focus of attention with measured photon emissions above 1\,PeV. In the absence of hadronic processes, this may indicate the existence of an extreme accelerator of e$^+$/e$^-$. Utilization of the CR extensive air shower detection techniques broadens the field of view of the source observations, enabling the measurement of UHE radiation surrounding the sources. These observations can probe the particle propagation inside and outside the accelerators and the subsequent injection/escape into the interstellar medium. △ Less

Submitted 2 October, 2023; originally announced October 2023.

Comments: 23 pages, 6 figures, to appear in Annual Review of Nuclear and Particle Science

Journal ref: Annual Review of Nuclear and Particle Science, Volume 73 (2023), page 341-363

Impacts of Gravitational-Wave Background from Supermassive Black Hole Binaries on the Detection of Compact Binaries by LISA

Authors: Fan Huang, Yan-Chen Bi, Zhoujian Cao, Qing-Guo Huang

Abstract: In the frequency band of Laser Interferometer Space Antenna (LISA), extensive research has been conducted on the impact of foreground confusion noise generated by galactic binaries within the Milky Way galaxy. Additionally, the recent evidence for a stochastic signal, announced by the NANOGrav, EPTA, PPTA, CPTA and InPTA, indicates that the stochastic gravitational-wave background generated by sup… ▽ More In the frequency band of Laser Interferometer Space Antenna (LISA), extensive research has been conducted on the impact of foreground confusion noise generated by galactic binaries within the Milky Way galaxy. Additionally, the recent evidence for a stochastic signal, announced by the NANOGrav, EPTA, PPTA, CPTA and InPTA, indicates that the stochastic gravitational-wave background generated by supermassive black hole binaries (SMBHBs) can contribute a strong background noise within in LISA band. Given the presence of such strong noise, it is expected to have a considerable impacts on LISA's scientific missions. In this work, we investigate the impacts of the SGWB generated by SMBHBs on the detection of massive black hole binaries (MBHBs), verified galactic binaries (VGBs) and extreme mass ratio inspirals (EMRIs) in the context of LISA, and find it crucial to resolve and eliminate the exceed noise from the SGWB to ensure the success of LISA's missions. △ Less

Submitted 25 September, 2023; originally announced September 2023.

Comments: 6 pages, 3 figures

Fast resolving Galactic binaries in LISA data and its ability to study the Milky Way

Authors: Pin Gao, Xi-Long Fan, Zhou-Jian Cao, Xue-Hao Zhang

Abstract: Resolving individual gravitational waves from tens of millions of double white dwarf (DWD) binaries in the Milky Way is a challenge for future space-based gravitational wave detection programs. By using previous data to define the priors for the next search, we propose an accelerated approach of searching the DWD binaries and demonstrate its efficiency based on the GBSIEVER detection pipeline. Com… ▽ More Resolving individual gravitational waves from tens of millions of double white dwarf (DWD) binaries in the Milky Way is a challenge for future space-based gravitational wave detection programs. By using previous data to define the priors for the next search, we propose an accelerated approach of searching the DWD binaries and demonstrate its efficiency based on the GBSIEVER detection pipeline. Compared to the traditional GBSIEVER method, our method can obtain $\sim 50\%$ of sources with 2.5\% of the searching time for LDC1-4 data. In addition, we find that both methods have a similar ability to detect the Milky Way structure by their confirmed sources. The relative error of distance and chirp mass is about 20\% for DWD binaries whose gravitational wave frequency is higher than $4\times10^{-3}$ Hz, even if they are close to the Galactic center. Finally, we propose a signal-to-noise ratio (SNR) threshold for LISA to confirm the detection of DWD binaries. The threshold should be 16 when the gravitational wave frequency is lower than $4\times10^{-3}$ Hz and 9 when the frequency range is from $4\times10^{-3}$ Hz to $1.5\times10^{-2}$ Hz. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: 16 pages, 19 figures

Journal ref: Phys. Rev. D 107, 123029, 2023

Matched filtering for gravitational wave detection without template bank driven by deep learning template prediction model bank

Authors: CunLiang Ma, Sen Wang, Wei Wang, Zhoujiang Cao

Abstract: The existing matched filtering method for gravitational wave (GW) search relies on a template bank. The computational efficiency of this method scales with the size of the templates within the bank. Higher-order modes and eccentricity will play an important role when third-generation detectors operate in the future. In this case, traditional GW search methods will hit computational limits. To spee… ▽ More The existing matched filtering method for gravitational wave (GW) search relies on a template bank. The computational efficiency of this method scales with the size of the templates within the bank. Higher-order modes and eccentricity will play an important role when third-generation detectors operate in the future. In this case, traditional GW search methods will hit computational limits. To speed up the computational efficiency of GW search, we propose the utilization of a deep learning (DL) model bank as a substitute for the template bank. This model bank predicts the latent templates embedded in the strain data. Combining an envelope extraction network and an astrophysical origin discrimination network, we realize a novel GW search framework. The framework can predict the GW signal's matched filtering signal-to-noise ratio (SNR). Unlike the end-to-end DL-based GW search method, our statistical SNR holds greater physical interpretability than the $p_{score}$ metric. Moreover, the intermediate results generated by our approach, including the predicted template, offer valuable assistance in subsequent GW data processing tasks such as parameter estimation and source localization. Compared to the traditional matched filtering method, the proposed method can realize real-time analysis. The minor improvements in the future, the proposed method may expand to other scopes of GW search, such as GW emitted by the supernova explosion. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 20 pages,12 figures

Journal ref: Physical Review D, 2024, 109(4): 043009

Neutron Star vs Quark Star in the Multimessenger Era

Authors: Zheng Cao, Lie-Wen Chen

Abstract: Neutron stars (NSs) which could contain exotic degrees of freedom in the core and the self-bound quark stars (QSs) made purely of absolutely stable deconfined quark matter are still two main candidates for the compact objects observed in pulsars and gravitational wave (GW) events in binary star mergers. We perform a Bayesian model-agnostic inference of the properties of NSs and QSs by combining mu… ▽ More Neutron stars (NSs) which could contain exotic degrees of freedom in the core and the self-bound quark stars (QSs) made purely of absolutely stable deconfined quark matter are still two main candidates for the compact objects observed in pulsars and gravitational wave (GW) events in binary star mergers. We perform a Bayesian model-agnostic inference of the properties of NSs and QSs by combining multi-messenger data of GW170817, GW190425, PSR J0030+0451, PSR J0740+6620, PSR J1614-2230, PSR J0348+0432 as well as ab initio calculations from perturbative quantum chromodynamics and chiral effective field theory. We find the NS scenario is strongly favored against the QS scenario with a Bayes factor of NS over QS $\mathcal{B}^\text{NS}_\text{QS} = 11.5$. In addition, the peak of the squared sound velocity $c_s^2 \sim 0.5c^2$ around $3.5$ times nuclear saturation density $n_0$ observed in the NS case disappears in the QS case which suggests that the $c_s^2$ first increases and then saturates at $c_s^2 \sim 0.5c^2$ above $\sim 4n_0$. The sound velocity and trace anomaly are found to approach the conformal limit in the core of heavy NSs with mass $M \gtrsim 2M_{\odot}$, but not in the core of QSs. △ Less

Submitted 31 August, 2023; originally announced August 2023.

Comments: 6 pages, 2 figures

Dilated convolutional neural network for detecting extreme-mass-ratio inspirals

Authors: Tianyu Zhao, Yue Zhou, Ruijun Shi, Zhoujian Cao, Zhixiang Ren

Abstract: The detection of Extreme Mass Ratio Inspirals (EMRIs) is intricate due to their complex waveforms, extended duration, and low signal-to-noise ratio (SNR), making them more challenging to be identified compared to compact binary coalescences. While matched filtering-based techniques are known for their computational demands, existing deep learning-based methods primarily handle time-domain data and… ▽ More The detection of Extreme Mass Ratio Inspirals (EMRIs) is intricate due to their complex waveforms, extended duration, and low signal-to-noise ratio (SNR), making them more challenging to be identified compared to compact binary coalescences. While matched filtering-based techniques are known for their computational demands, existing deep learning-based methods primarily handle time-domain data and are often constrained by data duration and SNR. In addition, most existing work ignores time-delay interferometry (TDI) and applies the long-wavelength approximation in detector response calculations, thus limiting their ability to handle laser frequency noise. In this study, we introduce DECODE, an end-to-end model focusing on EMRI signal detection by sequence modeling in the frequency domain. Centered around a dilated causal convolutional neural network, trained on synthetic data considering TDI-1.5 detector response, DECODE can efficiently process a year's worth of multichannel TDI data with an SNR of around 50. We evaluate our model on 1-year data with accumulated SNR ranging from 50 to 120 and achieve a true positive rate of 96.3% at a false positive rate of 1%, keeping an inference time of less than 0.01 seconds. With the visualization of three showcased EMRI signals for interpretability and generalization, DECODE exhibits strong potential for future space-based gravitational wave data analyses. △ Less

Submitted 14 May, 2024; v1 submitted 30 August, 2023; originally announced August 2023.

Comments: 11 pages, 5 figures, and 2 tables

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

Including higher harmonics in gravitational-wave parameter estimation and cosmological implications for LISA

Authors: Yi Gong, Zhoujian Cao, Junjie Zhao, Lijing Shao

Abstract: Massive black holes (MBHs) are crucial in shaping their host galaxies. How the MBH co-evolves with its host galaxy is a pressing problem in astrophysics and cosmology. The valuable information carried by the binary MBH is encoded in the gravitational waves (GWs), which will be detectable by the space-borne GW detector LISA. In the GW data analysis, usually, only the dominant $(2,2)$ mode of the GW… ▽ More Massive black holes (MBHs) are crucial in shaping their host galaxies. How the MBH co-evolves with its host galaxy is a pressing problem in astrophysics and cosmology. The valuable information carried by the binary MBH is encoded in the gravitational waves (GWs), which will be detectable by the space-borne GW detector LISA. In the GW data analysis, usually, only the dominant $(2,2)$ mode of the GW signal is considered in the parameter estimation for LISA. However, including the higher harmonics in parameter estimation can break the degeneracy between the parameters, especially for the inclination angle and luminosity distance. This may enable the identification of GW signals without electromagnetic counterparts, known as ''dark sirens''. Thus, incorporating higher harmonics will be beneficial to resolve the Hubble tension and constrain the cosmological model. In this paper, we investigate the role of higher harmonics in the parameter estimation for GWs emitted by binary MBHs. We demonstrate that including $(3,3)$ mode can lead to a $10^3$-times improvement in angular resolution and a $10^4$-times improvement in luminosity distance. Meanwhile, our results indicate that considering higher harmonics increases the probability of identifying over 70% host galaxies from $10^{-2}\,\rm{Gpc}^3$ cosmological volume threshold (corresponding $10^5$ host galaxies), while the probability less than 8% for only the $(2,2)$ mode. Thus, our results underscore the importance of including higher modes in the GW signal from binary MBHs, for LISA at least $(3,3)$ mode. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Comments: 14 pages, 11 figures, 4 tables; accepted by Physical Review D

Journal ref: Phys. Rev. D 108 (2023) 064046

A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run

Authors: C. Fletcher, J. Wood, R. Hamburg, P. Veres, C. M. Hui, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, A. Goldstein, B. A. Hristov, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, S. Poolakkil, A. von Kienlin, C. A. Wilson-Hodge, The Fermi Gamma-ray Burst Monitor Team, M. Crnogorčević, J. DeLaunay, A. Tohuvavohu, R. Caputo, S. B. Cenko , et al. (1674 additional authors not shown)

Abstract: We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,… ▽ More We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Efficient production of nuclear isomer $^{93m}$Mo with laser-accelerated proton beam and an astrophysical implication on $^{92m}$Mo production

Authors: Wenru Fan, Wei Qi, Jingli Zhang, Zongwei Cao, Haoyang Lan, Xinxiang Li, Yi Xu, Yuqiu Gu, Zhigang Deng, Zhimeng Zhang, Changxiang Tan, Wen Luo, Yun Yuan, Weimin Zhou

Abstract: Nuclear isomers play a key role in the creation of the elements in the universe and have a number of fascinating potential applications related to the controlled release of nuclear energy on demand. Particularly, $^{93m}$Mo isomer is a good candidate for studying the depletion of nuclear isomer via nuclear excitation by electron capture. For such purposes, efficient approach for $^{93m}$Mo product… ▽ More Nuclear isomers play a key role in the creation of the elements in the universe and have a number of fascinating potential applications related to the controlled release of nuclear energy on demand. Particularly, $^{93m}$Mo isomer is a good candidate for studying the depletion of nuclear isomer via nuclear excitation by electron capture. For such purposes, efficient approach for $^{93m}$Mo production needs to be explored. In the present work, we demonstrate experimentally an efficient production of $^{93m}$Mo through $^{93}$Nb(p, n) reaction induced by intense laser pulse. When a ps-duration, 100-J laser pulse is employed, the $^{93m}$Mo isomer at 2425 keV (21/2$^+$, $T_{1/2}$ = 6.85 h) are generated with a high yield of $1.8\times10^6$ particles/shot. The resulting peak efficiency is expected to be $10^{17}$ particles/s, which is at least five orders of magnitudes higher than using classical proton accelerator. The effects of production and destruction of $^{93m}$Mo on the controversial astrophysical p-isotope $^{92}$Mo are studied. It is found that the $^{93}$Nb(p, n)-$^{93m}$Mo reaction is an important production path for ^{93m}Mo seed nucleus, and the influence of ^{93m}Mo-^{92}Mo reaction flow on ^{92}Mo production cannot be ignored. In addition, we propose to directly measure the astrophysical rate of (p, n) reaction using laser-induced proton beam since the latter one fits the Maxwell-Boltzmann distribution well. We conclude that laser-induced proton beam opens a new path to produce nuclear isomers with high peak efficiency towards the understanding of p-nuclei nucleosythesis. △ Less

Submitted 5 August, 2023; originally announced August 2023.

Comments: 21 pages, 9 figure

Probe the gravitational constant variation via the propagation of gravitational waves

Authors: Bing Sun, Jiachen An, Zhoujian Cao

Abstract: The gravitational constant variation means the breakdown of the strong equivalence principle. As the cornerstone of general relativity, the validity of general relativity can be examined by studying the gravitational constant variation. Such variations have the potential to affect both the generation and propagation of gravitational waves. In this paper, our focus lies on the effect of gravitation… ▽ More The gravitational constant variation means the breakdown of the strong equivalence principle. As the cornerstone of general relativity, the validity of general relativity can be examined by studying the gravitational constant variation. Such variations have the potential to affect both the generation and propagation of gravitational waves. In this paper, our focus lies on the effect of gravitational constant variation specifically on the propagation of gravitational waves. We employ two analytical methods, namely based on the Fierz-Pauli action and the perturbation of Einstein-Hilbert action around Minkowski spacetime, both leading to the the same gravitational wave equation. By solving this equation, we find the effects of gravitational constant variation on gravitational wave propagation. The result is consistent with previous investigations based on Maxwell-like equations for gravitational waves. Notably, we find that small variations in the gravitational constant result in an amplitude correction at the leading order and a phase correction at the sub-leading order for gravitational waves. These results provide valuable insights for probing gravitational constant variation and can be directly applied to gravitational wave data analysis. △ Less

Submitted 31 July, 2023; originally announced August 2023.

Comments: 9 pages, 1 figure

Journal ref: Physics Letters B, Volume 848, January 2024, 138350

Post-Keplerian waveform model for binary compact object as sources of space-based gravitational wave detector and its implications

Authors: Li-Fang Li, Zhoujian Cao

Abstract: Binary compact objects will be among the important sources for the future space-based gravitational wave detectors. Such binary compact objects include stellar massive binary black hole, binary neutron star, binary white dwarf and mixture of these compact objects. Regarding to the relatively low frequency, the gravitational interaction between the two objects of the binary is weak. Post-Newtonian… ▽ More Binary compact objects will be among the important sources for the future space-based gravitational wave detectors. Such binary compact objects include stellar massive binary black hole, binary neutron star, binary white dwarf and mixture of these compact objects. Regarding to the relatively low frequency, the gravitational interaction between the two objects of the binary is weak. Post-Newtonian approximation of general relativity is valid. Previous works about the waveform model for such binaries in the literature consider the dynamics for specific situations which involve detailed complicated matter dynamics between the two objects. We here take a different idea. We adopt the trick used in pulsar timing detection. For any gravity theories and any detailed complicated matter dynamics, the motion of the binary can always be described as a post-Keplerian expansion. And a post-Keplerian gravitational waveform model will be reduced. Instead of object masses, spins, matter's equation of state parameters and dynamical parameters beyond general relativity, the involved parameters in our post-Keplerian waveform model are the Keplerian orbit elements and their adiabatic variations. Respect to current planning space-based gravitational wave detectors including LISA, Taiji and Tianqin, we find that the involved waveform model parameters can be well determined. And consequently the detail matter dynamics of the binary can be studied then. For binary with purely gravitational interactions, gravity theory can be constrained well. △ Less

Submitted 17 June, 2023; originally announced June 2023.

Comments: 14 pages, 3 figures

Journal ref: General Relativity and Gravitation, 55, 76 (2023)

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)

Site-testing at the Muztagh-ata Site.V. Nighttime Cloud Amount during the Last Five Years

Authors: Jing Xu, Guo-jie Feng, Guang-xin Pu, Le-tian Wang, Zi-Huang Cao, Li-Qing Ren, Xuan Zhang, Shu-guo Ma, Chun-hai Bai, Ali Esamdin, Jian Li, Yuan Tian, Zheng Wang, Yong-heng Zhao, Jian-rong Shi

Abstract: The clarity of nights is the major factor that should be carefully considered for optical/infrared astronomical observatories in site-testing campaigns. Cloud coverage is directly related to the amount of time available for scientific observations at observatories. In this article, we report on the results of detailed night-time cloud statistics and continuous observing me derived from ground-base… ▽ More The clarity of nights is the major factor that should be carefully considered for optical/infrared astronomical observatories in site-testing campaigns. Cloud coverage is directly related to the amount of time available for scientific observations at observatories. In this article, we report on the results of detailed night-time cloud statistics and continuous observing me derived from ground-based all-sky cameras at the Muztagh-ata site from 2017 to 2021. Results obtained from acquisition data show that the proportion of the annual observing me at the Muztagh-ata site is 65%, and the best period with the least cloud coverage and longer continuous observing time is from September to February. We made a comparison of the monthly mean observing nights obtained from our all-sky cameras and CLARA dataset, results show that the discrepancy between them may depend on the cloud top heights. On average, this site can provide 175 clear nights and 169 nights with at least 4 hours of continuous observing time per year. △ Less

Submitted 8 May, 2023; v1 submitted 4 May, 2023; originally announced May 2023.

Comments: 11 pages, 9 figures

Journal ref: [J]. Research in Astronomy and Astrophysics, 2023, 23(4): 045015

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated… ▽ More Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 28 pages, 11 figures

Report number: LIGO-P2200031

Turbulence in Zeeman Measurements from Molecular Clouds

Authors: Zhuo Cao, Hua-bai Li

Abstract: Magnetic fields (B-fields) play an important role in molecular cloud fragmentation and star formation, but are very difficult to detect. The temporal correlation between the field strength (B) and gas density (n) of an isolated cloud has been suggested as an indication of the dynamical importance of B-fields relative to self-gravity. This temporal B-n relation is, however, unobservable. What can b… ▽ More Magnetic fields (B-fields) play an important role in molecular cloud fragmentation and star formation, but are very difficult to detect. The temporal correlation between the field strength (B) and gas density (n) of an isolated cloud has been suggested as an indication of the dynamical importance of B-fields relative to self-gravity. This temporal B-n relation is, however, unobservable. What can be observed using Zeeman measurements are the "spatial B-n relations" from the current plane of the sky. Nevertheless, the temporal B-n relation argument has still been widely used to interpret observations. Here we present the first numerical test of the legitimacy of this interpretation. From a simulation that can reproduce the observed Zeeman spatial B~n^2/3 relation, we found that temporal B-n relations of individual cores bear no resemblance to the spatial B-n relations. This result inspired us to discover that the true mechanism behind the 2/3 index is random turbulence compression instead of symmetrical gravitational contraction. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: accepted by ApJL

L dwarfs detection from SDSS images using improved Faster R-CNN

Authors: Zhi Cao, Zhenping Yi, Jingchang Pan, Hao Su, Yude Bu, Xiao Kong, Ali Luo

Abstract: We present a data-driven approach to automatically detect L dwarfs from Sloan Digital Sky Survey(SDSS) images using an improved Faster R-CNN framework based on deep learning. The established L dwarf automatic detection (LDAD) model distinguishes L dwarfs from other celestial objects and backgrounds in SDSS field images by learning the features of 387 SDSS images containing L dwarfs. Applying the L… ▽ More We present a data-driven approach to automatically detect L dwarfs from Sloan Digital Sky Survey(SDSS) images using an improved Faster R-CNN framework based on deep learning. The established L dwarf automatic detection (LDAD) model distinguishes L dwarfs from other celestial objects and backgrounds in SDSS field images by learning the features of 387 SDSS images containing L dwarfs. Applying the LDAD model to the SDSS images containing 93 labeled L dwarfs in the test set, we successfully detected 83 known L dwarfs with a recall rate of 89.25% for known L dwarfs. Several techniques are implemented in the LDAD model to improve its detection performance for L dwarfs,including the deep residual network and the feature pyramid network. As a result, the LDAD model outperforms the model of the original Faster R-CNN, whose recall rate of known L dwarfs is 80.65% for the same test set. The LDAD model was applied to detect L dwarfs from a larger validation set including 843 labeled L dwarfs, resulting in a recall rate of 94.42% for known L dwarfs. The newly identified candidates include L dwarfs, late M and T dwarfs, which were estimated from color (i-z) and spectral type relation. The contamination rates for the test candidates and validation candidates are 8.60% and 9.27%, respectively. The detection results indicate that our model is effective to search for L dwarfs from astronomical images. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: 12 pages, 10 figures, accepted to be published in AJ

Measuring mass transfer of AM CVn binaries with a space-based gravitational wave detector

Authors: Zijian Wang, Zhoujian Cao, Xian-Fei Zhang

Abstract: The formation mechanism of AM CVn binary has not been well understood yet. Accurate measurements of the mass transfer rate can help to determine the formation mechanism. But unfortunately such observation by electromagnetic means is quite challenging. One possible formation channel of AM CVn binary is a semi-detached white dwarf binary. Such system emits strong gravitational wave radiation which c… ▽ More The formation mechanism of AM CVn binary has not been well understood yet. Accurate measurements of the mass transfer rate can help to determine the formation mechanism. But unfortunately such observation by electromagnetic means is quite challenging. One possible formation channel of AM CVn binary is a semi-detached white dwarf binary. Such system emits strong gravitational wave radiation which could be measured by the future space-based detectors. We can simultaneously extract the mass transfer rate and the orbital period from the gravitational wave signal. We employ a post-Keplerian waveform model of gravitational wave and carry out a Fisher analysis to estimate the measurement accuracy of mass transfer rate through gravitational wave detection. Special attention is paid to the observed sources in Gaia Data Release 2. We found that we can accurately measure the mass transfer rate for those systems. Comparing to electromagnetic observations, gravitational wave detection improves the accuracy more than two orders of magnitude. Our results imply that the gravitational wave detection will help much in understanding the formation mechanism of AM CVn binaries. △ Less

Submitted 17 August, 2023; v1 submitted 17 February, 2023; originally announced February 2023.

Lorentz transformation of three dimensional gravitational wave tensor

Authors: Xiaokai He, Xiaolin Liu, Zhoujian Cao

Abstract: Recently there are more and more interest on the gravitational wave of moving sources. This introduces a Lorentz transformation problem of gravitational wave. Although Bondi-Metzner-Sachs (BMS) theory has in principle already included the Lorentz transformation of gravitational wave, the transformation of the three dimensional gravitational wave tensor has not been explicitly calculated before. Wi… ▽ More Recently there are more and more interest on the gravitational wave of moving sources. This introduces a Lorentz transformation problem of gravitational wave. Although Bondi-Metzner-Sachs (BMS) theory has in principle already included the Lorentz transformation of gravitational wave, the transformation of the three dimensional gravitational wave tensor has not been explicitly calculated before. Within four dimensional spacetime, gravitational wave have property of `boost weight zero' and `spin weight 2'. This fact makes the Lorentz transformation of gravitational wave difficult to understand. In the current paper we adopt the traditional three dimensional tensor description of gravitational wave. Such a transverse-traceless tensor describes the gravitational wave freedom directly. We derive the explicit Lorentz transformation of the gravitational wave tensor. The transformation is similar to the Lorentz transformation for electric field vector and magnetic field vector which are three dimensional vectors. Based on the deduced Lorentz transformation of the gravitational wave three dimensional tensor, we can construct the gravitational waveform of moving source with any speed if only the waveform of the corresponding rest waveform is given. As an example, we apply our method to the effect of kick velocity of binary black hole. The adjusted waveform by the kick velocity is presented. △ Less

Submitted 15 February, 2023; originally announced February 2023.

Comments: 17 pages, 8 figures

Journal ref: CQG 39, 195020 (2022)