Showing 1–50 of 94 results for author: Ma, Q

Inferring IGM parameters from the redshifted 21-cm Power Spectrum using Artificial Neural Networks

Authors: Madhurima Choudhury, Raghunath Ghara, Saleem Zaroubi, Benedetta Ciardi, Leon V. E. Koopmans, Garrelt Mellema, Abinash Kumar Shaw, Anshuman Acharya, I. T. Iliev, Qing-Bo Ma, Sambit K. Giri

Abstract: The high redshift 21-cm signal promises to be a crucial probe of the state of the intergalactic medium (IGM). Understanding the connection between the observed 21-cm power spectrum and the physical quantities intricately associated with the IGM is crucial to fully understand the evolution of our Universe. In this study, we develop an emulator using artificial neural network (ANN) to predict the 21… ▽ More The high redshift 21-cm signal promises to be a crucial probe of the state of the intergalactic medium (IGM). Understanding the connection between the observed 21-cm power spectrum and the physical quantities intricately associated with the IGM is crucial to fully understand the evolution of our Universe. In this study, we develop an emulator using artificial neural network (ANN) to predict the 21-cm power spectrum from a given set of IGM properties, namely, the bubble size distribution and the volume averaged ionization fraction. This emulator is implemented within a standard Bayesian framework to constrain the IGM parameters from a given 21-cm power spectrum. We compare the performance of the Bayesian method to an alternate method using ANN to predict the IGM parameters from a given input power spectrum, and find that both methods yield similar levels of accuracy, while the ANN is significantly faster. We also use this ANN method of parameter estimation to predict the IGM parameters from a test set contaminated with noise levels expected from the SKA-LOW instrument after 1000 hours of observation. Finally, we train a separate ANN to predict the source parameters from the IGM parameters directly, at a redshift of $z=9.1$, demonstrating the possibility of a non-analytic inference of the source parameters from the IGM parameters for the first time. We achieve high accuracies, with R2-scores ranging between $0.898-0.978$ for the ANN emulator and between $0.966-0.986$ and $0.817-0.981$ for the predictions of IGM parameters from 21-cm power spectrum and source parameters from IGM parameters, respectively. The predictions of the IGM parameters from the Bayesian method incorporating the ANN emulator leads to tight constraints with error bars around $\pm{0.14}$ on the IGM parameters. △ Less

Submitted 3 July, 2024; originally announced July 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

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

A quantification of the effects using different stellar population synthesis models for epoch of reionization

Authors: Peiai Liu, Qingbo Ma, Yunkun Han, Rongxin Luo

Abstract: The luminosity and spectral energy distribution (SED) of high-$z$ galaxies are sensitive to the stellar population synthesis (SPS) models. In this paper, we study the effects of different SPS models on the measurements of high-$z$ galaxies and the budget of ionizing photons during the epoch of reionization, by employing each of them in the semi-analytical galaxy formation model {\sc L-Galaxies 202… ▽ More The luminosity and spectral energy distribution (SED) of high-$z$ galaxies are sensitive to the stellar population synthesis (SPS) models. In this paper, we study the effects of different SPS models on the measurements of high-$z$ galaxies and the budget of ionizing photons during the epoch of reionization, by employing each of them in the semi-analytical galaxy formation model {\sc L-Galaxies 2020}. We find that the different SPS models lead to $\lesssim 0.5$ dex differences on the amplitudes of UV luminosity functions, while the two modes of the same SPS model with and without the inclusion of binary stars leads to similar UV luminosity functions at $z \ge 6$. Instead, the binary stars produce $\sim 40\%$ more ionizing photons than the single stars, while such differences are smaller than those caused by different SPS models, e.g. the BPASS model produces $\sim 100\%$ more ionizing photons than other models. △ Less

Submitted 2 May, 2024; originally announced May 2024.

Comments: 12 pages, 7 figures, APJ accepted

Probing the intergalactic medium during the Epoch of Reionization using 21-cm signal power spectra

Authors: Raghunath Ghara, Abinash Kumar Shaw, Saleem Zaroubi, Benedetta Ciardi, Garrelt Mellema, Léon V. E. Koopmans, Anshuman Acharya, Madhurima Choudhury, Sambit K. Giri, Ilian T. Iliev, Qing-Bo Ma, Florent Mertens

Abstract: The redshifted 21-cm signal from the epoch of reionization (EoR) directly probes the ionization and thermal states of the intergalactic medium during that period. In particular, the distribution of the ionized regions around the radiating sources during EoR introduces scale-dependent features in the spherically-averaged EoR 21-cm signal power spectrum. The goal is to study these scale-dependent fe… ▽ More The redshifted 21-cm signal from the epoch of reionization (EoR) directly probes the ionization and thermal states of the intergalactic medium during that period. In particular, the distribution of the ionized regions around the radiating sources during EoR introduces scale-dependent features in the spherically-averaged EoR 21-cm signal power spectrum. The goal is to study these scale-dependent features at different stages of reionization using numerical simulations and build a source model-independent framework to probe the properties of the intergalactic medium using EoR 21-cm signal power spectrum measurements. Under the assumption of high spin temperature, we modelled the redshift evolution of the ratio of EoR 21-cm brightness temperature power spectrum and the corresponding density power spectrum using an ansatz consisting of a set of redshift and scale-independent parameters. This set of eight parameters probes the redshift evolution of the average ionization fraction and the quantities related to the morphology of the ionized regions. We have tested this ansatz on different reionization scenarios generated using different simulation algorithms and found that it is able to recover the redshift evolution of the average neutral fraction within an absolute deviation $\lesssim 0.1$. Our framework allows us to interpret 21-cm signal power spectra in terms of parameters related to the state of the IGM. This source model-independent framework can efficiently constrain reionization scenarios using multi-redshift power spectrum measurements with ongoing and future radio telescopes such as LOFAR, MWA, HERA, and SKA. This will add independent information regarding the EoR IGM properties. △ Less

Submitted 17 April, 2024; originally announced April 2024.

Comments: 16 pages, 13 figures, 2 tables, Accepted for publication in Astronomy and Astrophysics

Report number: NORDITA 2024-009

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

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

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

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

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

LHAASO-KM2A detector simulation using Geant4

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

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

Submitted 7 April, 2024; originally announced April 2024.

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

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

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

Submitted 20 March, 2024; originally announced March 2024.

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

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

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

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

Submitted 15 March, 2024; originally announced March 2024.

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

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

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

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

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

Comments: 8 pages, 3 figures

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

21-cm Signal from the Epoch of Reionization: A Machine Learning upgrade to Foreground Removal with Gaussian Process Regression

Authors: Anshuman Acharya, Florent Mertens, Benedetta Ciardi, Raghunath Ghara, Léon V. E. Koopmans, Sambit K. Giri, Ian Hothi, Qing-Bo Ma, Garrelt Mellema, Satyapan Munshi

Abstract: In recent years, a Gaussian Process Regression (GPR) based framework has been developed for foreground mitigation from data collected by the LOw-Frequency ARray (LOFAR), to measure the 21-cm signal power spectrum from the Epoch of Reionization (EoR) and Cosmic Dawn. However, it has been noted that through this method there can be a significant amount of signal loss if the EoR signal covariance is… ▽ More In recent years, a Gaussian Process Regression (GPR) based framework has been developed for foreground mitigation from data collected by the LOw-Frequency ARray (LOFAR), to measure the 21-cm signal power spectrum from the Epoch of Reionization (EoR) and Cosmic Dawn. However, it has been noted that through this method there can be a significant amount of signal loss if the EoR signal covariance is misestimated. To obtain better covariance models, we propose to use a kernel trained on the {\tt GRIZZLY} simulations using a Variational Auto-Encoder (VAE) based algorithm. In this work, we explore the abilities of this Machine Learning based kernel (VAE kernel) used with GPR, by testing it on mock signals from a variety of simulations, exploring noise levels corresponding to $\approx$10 nights ($\approx$141 hours) and $\approx$100 nights ($\approx$1410 hours) of observations with LOFAR. Our work suggests the possibility of successful extraction of the 21-cm signal within 2$σ$ uncertainty in most cases using the VAE kernel, with better recovery of both shape and power than with previously used covariance models. We also explore the role of the excess noise component identified in past applications of GPR and additionally analyse the possibility of redshift dependence on the performance of the VAE kernel. The latter allows us to prepare for future LOFAR observations at a range of redshifts, as well as compare with results from other telescopes. △ Less

Submitted 29 November, 2023; v1 submitted 28 November, 2023; originally announced November 2023.

Comments: 13 pages, 7 figures, 3 tables. Accepted for publication in the Monthly Notices of the Royal Astronomical Society

Report number: NORDITA 2023-074

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

Cosmic variance suppression in radiation-hydrodynamic modeling of the reionization-era 21-cm signal

Authors: Anshuman Acharya, Enrico Garaldi, Benedetta Ciardi, Qing-bo Ma

Abstract: The 21-cm line emitted by neutral hydrogen is the most promising probe of the Epoch of Reionization (EoR). Multiple radio interferometric instruments are on the cusp of detecting its power spectrum. It is therefore essential to deliver robust theoretical predictions, enabling sound inference of the coeval Universe properties. The nature of this signal traditionally required the modelling of… ▽ More The 21-cm line emitted by neutral hydrogen is the most promising probe of the Epoch of Reionization (EoR). Multiple radio interferometric instruments are on the cusp of detecting its power spectrum. It is therefore essential to deliver robust theoretical predictions, enabling sound inference of the coeval Universe properties. The nature of this signal traditionally required the modelling of $\mathcal{O}(10^{7-8} \, {\rm Mpc}^3)$ volumes to suppress the impact of cosmic variance. However, the recently-proposed Fixed & Paired (F&P) approach uses carefully-crafted simulation pairs to achieve equal results in smaller volumes. In this work, we thoroughly test the applicability of and improvement granted by this technique to different observables of the 21-cm signal from the EoR. We employ radiation-magneto-hydrodynamics simulations to ensure the most realistic physical description of this epoch, greatly improving over previous studies using a semi-numerical approach without accurate galaxy formation physics and radiative transfer. We estimate the statistical improvement granted by the F&P technique on predictions of the skewness, power spectrum, bispectrum and ionized regions size distribution of the 21-cm signal at redshift $7 \leq z \leq 10$ (corresponding to $\geq80\%$ of the gas being neutral). We find that the effective volume of F&P simulations is at least 3.5 times larger than traditional simulations. This directly translates into an equal improvement in the computational cost (in terms of time and memory). Finally, we confirm that a combination of different observables like skewness, power spectrum and bispectrum across different redshifts can be utilised to maximise the improvement. △ Less

Submitted 14 March, 2024; v1 submitted 20 October, 2023; originally announced October 2023.

Comments: 13 pages, 11 figures, 2 tables. Accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS)

CatNorth: An Improved Gaia DR3 Quasar Candidate Catalog with Pan-STARRS1 and CatWISE

Authors: Yuming Fu, Xue-Bing Wu, Yifan Li, Yuxuan Pang, Ravi Joshi, Shuo Zhang, Qiyue Wang, Jing Yang, FanLam Ng, Xingjian Liu, Yu Qiu, Rui Zhu, Huimei Wang, Christian Wolf, Yanxia Zhang, Zhi-Ying Huo, Y. L. Ai, Qinchun Ma, Xiaotong Feng, R. J. Bouwens

Abstract: A complete and pure sample of quasars with accurate redshifts is crucial for quasar studies and cosmology. In this paper, we present CatNorth, an improved Gaia DR3 quasar candidate catalog with more than 1.5 million sources in the 3$π$ sky built with data from Gaia, Pan-STARRS1, and CatWISE2020. The XGBoost algorithm is used to reclassify the original Gaia DR3 quasar candidates as stars, galaxies,… ▽ More A complete and pure sample of quasars with accurate redshifts is crucial for quasar studies and cosmology. In this paper, we present CatNorth, an improved Gaia DR3 quasar candidate catalog with more than 1.5 million sources in the 3$π$ sky built with data from Gaia, Pan-STARRS1, and CatWISE2020. The XGBoost algorithm is used to reclassify the original Gaia DR3 quasar candidates as stars, galaxies, and quasars. To construct training/validation datasets for the classification, we carefully built two different master stellar samples in addition to the spectroscopic galaxy and quasar samples. An ensemble classification model is obtained by averaging two XGBoost classifiers trained with different master stellar samples. Using a probability threshold of $p_{\mathrm{QSO\_mean}}>0.95$ in our ensemble classification model and an additional cut on the logarithmic probability density of zero proper motion, we retrieved 1,545,514 reliable quasar candidates from the parent Gaia DR3 quasar candidate catalog. We provide photometric redshifts for all candidates with an ensemble regression model. For a subset of 89,100 candidates, accurate spectroscopic redshifts are estimated with the Convolutional Neural Network from the Gaia BP/RP spectra. The CatNorth catalog has a high purity of ~ 90% while maintaining high completeness, which is an ideal sample to understand the quasar population and its statistical properties. The CatNorth catalog is used as the main source of input catalog for the LAMOST phase III quasar survey, which is expected to build a highly complete sample of bright quasars with $i < 19.5$. △ Less

Submitted 13 February, 2024; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 27 pages, 15 figures, accepted by ApJS. Table 4 (The CatNorth quasar candidate catalog) is available at https://nadc.china-vo.org/res/r101313/

Journal ref: Fu et al 2024 ApJS 271 54

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

Sulfur isotope ratios in the Large Magellanic Cloud

Authors: Y. Gong, C. Henkel, K. M. Menten, C. -H. R. Chen, Z. Y. Zhang, Y. T. Yan, A. Weiss, N. Langer, J. Z. Wang, R. Q. Mao, X. D. Tang, W. Yang, Y. P. Ao, M. Wang

Abstract: Sulfur isotope ratios have emerged as a promising tool for tracing stellar nucleosynthesis, quantifying stellar populations, and investigating the chemical evolution of galaxies. While extensively studied in the Milky Way, in extragalactic environments they remain largely unexplored. We focus on investigating the sulfur isotope ratios in the Large Magellanic Cloud (LMC) to gain insights into sulfu… ▽ More Sulfur isotope ratios have emerged as a promising tool for tracing stellar nucleosynthesis, quantifying stellar populations, and investigating the chemical evolution of galaxies. While extensively studied in the Milky Way, in extragalactic environments they remain largely unexplored. We focus on investigating the sulfur isotope ratios in the Large Magellanic Cloud (LMC) to gain insights into sulfur enrichment in this nearby system and to establish benchmarks for such ratios in metal-poor galaxies. We conducted pointed observations of CS and its isotopologues toward N113, one of the most prominent star-formation regions in the LMC, utilizing the Atacama Pathfinder EXperiment 12~m telescope. We present the first robust detection of C$^{33}$S in the LMC by successfully identifying two C$^{33}$S transitions on a large scale of $\sim$5 pc. Our measurements result in an accurate determination of the $^{34}$S/$^{33}$S isotope ratio, which is 2.0$\pm$0.2. Our comparative analysis indicates that the $^{32}$S/$^{33}$S and $^{34}$S/$^{33}$S isotope ratios are about a factor of 2 lower in the LMC than in the Milky Way. Our findings suggest that the low $^{34}$S/$^{33}$S isotope ratio in the LMC can be attributed to a combination of the age effect, low metallicity, and star formation history. △ Less

Submitted 18 October, 2023; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: 6 pages, 1 figures, 2 tables, accepted for publication in A&A, adjusted to the final version

Journal ref: A&A 679, L6 (2023)

Real-time Monitoring for the Next Core-Collapse Supernova in JUNO

Authors: Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Weidong Bai, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato, Marco Beretta, Antonio Bergnoli , et al. (606 additional authors not shown)

Abstract: The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neu… ▽ More The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neutrinos at the Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton liquid scintillator detector currently under construction in South China. The real-time monitoring system is designed to ensure both prompt alert speed and comprehensive coverage of progenitor stars. It incorporates prompt monitors on the electronic board as well as online monitors at the data acquisition stage. Assuming a false alert rate of 1 per year, this monitoring system exhibits sensitivity to pre-SN neutrinos up to a distance of approximately 1.6 (0.9) kiloparsecs and SN neutrinos up to about 370 (360) kiloparsecs for a progenitor mass of 30 solar masses, considering both normal and inverted mass ordering scenarios. The pointing ability of the CCSN is evaluated by analyzing the accumulated event anisotropy of inverse beta decay interactions from pre-SN or SN neutrinos. This, along with the early alert, can play a crucial role in facilitating follow-up multi-messenger observations of the next galactic or nearby extragalactic CCSN. △ Less

Submitted 4 December, 2023; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: 24 pages, 9 figures, accepted for the publication at JCAP

Individual and Averaged Power Density Spectra of X-ray bursts from SGR J1935+2154: Quasiperiodic Oscillation Search and Slopes

Authors: Shuo Xiao, Xiao-Bo Li, Wang-Chen Xue, Shao-Lin Xiong, Shuang-Nan Zhang, Wen-Xi Peng, Ai-Jun Dong, You-Li Tuo, Ce Cai, Xi-Hong Luo, Jiao-Jiao Yang, Yue Wang, Chao Zheng, Yan-Qiu Zhang, Jia-Cong Liu, Wen-Jun Tan, Chen-Wei Wang, Ping Wang, Cheng-Kui Li, Shu-Xu Yi, Shi-Jun Dang, Lun-Hua Shang, Ru-Shuang Zhao, Qing-Bo Ma, Wei Xie , et al. (7 additional authors not shown)

Abstract: The study of quasi-periodic oscillations (QPOs) and power density spectra (PDS) continuum properties can help shed light on the still illusive emission physics of magnetars and as a window into the interiors of neutron stars using asteroseismology. In this work, we employ a Bayesian method to search for the QPOs in the hundreds of X-ray bursts from SGR J1935+2154 observed by {\it Insight}-HXMT, GE… ▽ More The study of quasi-periodic oscillations (QPOs) and power density spectra (PDS) continuum properties can help shed light on the still illusive emission physics of magnetars and as a window into the interiors of neutron stars using asteroseismology. In this work, we employ a Bayesian method to search for the QPOs in the hundreds of X-ray bursts from SGR J1935+2154 observed by {\it Insight}-HXMT, GECAM and Fermi/GBM from July 2014 to January 2022. Although no definitive QPO signal (significance $>3σ$) is detected in individual bursts or the averaged periodogram of the bursts grouped by duration, we identify several bursts exhibiting possible QPO at $\sim$ 40 Hz, which is consistent with that reported in the X-ray burst associated with FRB 200428. We investigate the PDS continuum properties and find that the distribution of the PDS slope in the simple power-law model peaks $\sim$ 2.5, which is consistent with other magnetars but higher than 5/3 commonly seen in gamma-ray bursts. Besides, the distribution of the break frequency in the broken power-law model peaks at $\sim$ 60 Hz. Finally, we report that the power-law index of PDS has an anti-correlation and power-law dependence on the burst duration as well as the minimum variation timescale. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: comments welcome

The Minimum Variation Timescales of X-ray bursts from SGR J1935+2154

Authors: Shuo Xiao, Jiao-Jiao Yang, Xi-Hong Luo, Shao-Lin Xiong, Yuan-Hong Qu, Shuang-Nan Zhang, Wang-Chen Xue, Xiao-Bo Li, You-Li Tuo, Ai-Jun Dong, Ru-Shuang Zhao, Shi-Jun Dang, Lun-Hua Shang, Qing-Bo Ma, Ce Cai, Jin Wang, Ping Wang, Cheng-Kui Li, Shu-Xu Yi, Zhen Zhang, Ming-Yu Ge, Shi-Jie Zheng, Li-Ming Song, Wen-Xi Peng, Xiang-Yang Wen , et al. (12 additional authors not shown)

Abstract: The minimum variation timescale (MVT) of soft gamma-ray repeaters can be an important probe to estimate the emission region in pulsar-like models, as well as the Lorentz factor and radius of the possible relativistic jet in gamma-ray burst (GRB)-like models, thus revealing their progenitors and physical mechanisms. In this work, we systematically study the MVTs of hundreds of X-ray bursts (XRBs) f… ▽ More The minimum variation timescale (MVT) of soft gamma-ray repeaters can be an important probe to estimate the emission region in pulsar-like models, as well as the Lorentz factor and radius of the possible relativistic jet in gamma-ray burst (GRB)-like models, thus revealing their progenitors and physical mechanisms. In this work, we systematically study the MVTs of hundreds of X-ray bursts (XRBs) from SGR J1935+2154 observed by {\it Insight}-HXMT, GECAM and Fermi/GBM from July 2014 to Jan 2022 through the Bayesian Block algorithm. We find that the MVTs peak at $\sim$ 2 ms, corresponding to a light travel time size of about 600 km, which supports the magnetospheric origin in pulsar-like models. The shock radius and the Lorentz factor of the jet are also constrained in GRB-like models. Interestingly, the MVT of the XRB associated with FRB 200428 is $\sim$ 70 ms, which is longer than that of most bursts and implies its special radiation mechanism. Besides, the median of MVTs is 7 ms, shorter than the median MVTs of 40 ms and 480 ms for short GRBs or long GRBs, respectively. However, the MVT is independent of duration, similar to GRBs. Finally, we investigate the energy dependence of MVT and suggest that there is a marginal evidence for a power-law relationship like GRBs but the rate of variation is at least about an order of magnitude smaller. These features may provide an approach to identify bursts with a magnetar origin. △ Less

Submitted 13 July, 2023; originally announced July 2023.

Comments: accepted for publication in ApJS

JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo

Authors: JUNO Collaboration, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Tsagkarakis Alexandros, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Weidong Bai, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato , et al. (581 additional authors not shown)

Abstract: We discuss JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo via detecting inverse beta decay reactions of electron anti-neutrinos resulting from the annihilation. We study possible backgrounds to the signature, including the reactor neutrinos, diffuse supernova neutrino background, charged- and neutral-current interactions of atmospheric neutrinos, backgrounds from muon… ▽ More We discuss JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo via detecting inverse beta decay reactions of electron anti-neutrinos resulting from the annihilation. We study possible backgrounds to the signature, including the reactor neutrinos, diffuse supernova neutrino background, charged- and neutral-current interactions of atmospheric neutrinos, backgrounds from muon-induced fast neutrons and cosmogenic isotopes. A fiducial volume cut, as well as the pulse shape discrimination and the muon veto are applied to suppress the above backgrounds. It is shown that JUNO sensitivity to the thermally averaged dark matter annihilation rate in 10 years of exposure would be significantly better than the present-day best limit set by Super-Kamiokande and would be comparable to that expected by Hyper-Kamiokande. △ Less

Submitted 13 September, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 25 pages, 9 figures, matches the publised version

Journal ref: JCAP 09 (2023) 001

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

Skew spectrum and smoothed skewness of 21-cm signals from epoch of reionization

Authors: Qing-Bo Ma, Ling Peng

Abstract: Due to the non-linear ionizing and heating processes, the 21-cm signals from epoch of reionization (EoR) are expected to have strong non-Gaussian fluctuations. In this paper, we use the semi-numerical simulations to study the non-Gaussian statistics i.e. skew spectrum and smoothed skewness of the 21-cm signals from EoR. We find the 21-cm skew spectrum and smoothed skewness have similar evolution f… ▽ More Due to the non-linear ionizing and heating processes, the 21-cm signals from epoch of reionization (EoR) are expected to have strong non-Gaussian fluctuations. In this paper, we use the semi-numerical simulations to study the non-Gaussian statistics i.e. skew spectrum and smoothed skewness of the 21-cm signals from EoR. We find the 21-cm skew spectrum and smoothed skewness have similar evolution features with the 21-cm bispectrum. All of them are sensitive to the EoR models, while not too much to the cosmic volume applied. With the SKA1-low telescope as reference, we find both the skew spectrum and smoothed skewness have much higher S/N ratios than the 21-cm bispectrum. △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: 7 pages, 5 figures, mnras accepted

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)

POLAR -- I: linking the 21-cm signal from the epoch of reionization to galaxy formation

Authors: Qing-Bo Ma, Raghunath Ghara, Benedetta Ciardi, Ilian T. Iliev, Léon V. E. Koopmans, Garrelt Mellema, Rajesh Mondal, Saleem Zaroubi

Abstract: To self-consistently model galactic properties, reionization of the intergalactic medium, and the associated 21-cm signal, we have developed the algorithm polar by integrating the one-dimensional radiative transfer code grizzly with the semi-analytical galaxy formation code L-Galaxies 2020. Our proof-of-concept results are consistent with observations of the star formation rate history, UV luminos… ▽ More To self-consistently model galactic properties, reionization of the intergalactic medium, and the associated 21-cm signal, we have developed the algorithm polar by integrating the one-dimensional radiative transfer code grizzly with the semi-analytical galaxy formation code L-Galaxies 2020. Our proof-of-concept results are consistent with observations of the star formation rate history, UV luminosity function and the CMB Thomson scattering optical depth. We then investigate how different galaxy formation models affect UV luminosity functions and 21-cm power spectra, and find that while the former are most sensitive to the parameters describing the merger of halos, the latter have a stronger dependence on the supernovae feedback parameters, and both are affected by the escape fraction model. △ Less

Submitted 19 April, 2023; originally announced April 2023.

Comments: MNRAS accepted

Pulsar and Magnetar Navigation with Fermi/GBM and GECAM

Authors: Xi-Hong Luo, Shuo Xiao, Shi-Jie Zheng, Ming-Yu Ge, You-Li Tuo, Shao-Lin Xiong, Shuang-Nan Zhang, Fang-Jun Lu, Yue Huang, Cheng Yang, Qi-Jun Zhi, Li-Ming Song, Wen-Xi Peng, Xiang-Yang Wen, Xin-Qiao Li, Zheng-Hua An, Jin Wang, Ping Wang, Ce Cai, Cheng-Kui Li, Xiao-Bo Li, Fan Zhang, Ai-Jun Dong, Wei Xie, Jian-Chao Feng , et al. (7 additional authors not shown)

Abstract: The determination of the absolute and relative position of a spacecraft is critical for its operation, observations, data analysis, scientific studies, as well as deep space exploration in general. A spacecraft that can determine its own absolute position autonomously may perform more than that must rely on transmission solutions. In this work, we report an absolute navigation accuracy of $\sim$ 2… ▽ More The determination of the absolute and relative position of a spacecraft is critical for its operation, observations, data analysis, scientific studies, as well as deep space exploration in general. A spacecraft that can determine its own absolute position autonomously may perform more than that must rely on transmission solutions. In this work, we report an absolute navigation accuracy of $\sim$ 20 km using 16-day Crab pulsar data observed with $Fermi$ Gamma ray Burst Monitor (GBM). In addition, we propose a new method with the inverse process of the triangulation for joint navigation using repeated bursts like that from the magnetar SGR J1935+2154 observed by the Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM) and GBM. △ Less

Submitted 25 March, 2023; originally announced March 2023.

Comments: accepted for publication in ApJS

Discovery of the linear energy-dependence of the spectral lag of X-ray bursts from SGR J1935+2154

Authors: Shuo Xiao, You-Li Tuo, Shuang-Nan Zhang, Shao-Lin Xiong, Lin Lin, Yan-Qiu Zhang, Yue Wang, Wang-Chen Xue, Ce Cai, He Gao, Cheng-Kui Li, Xiao-Bo Li, Chao Zheng, Jia-Cong Liu, Ping Wang, Jin Wang, Wen-Xi Peng, Cong Zhan Liu, Xin-Qiao Li, Xiang-Yang Wen, Zheng-Hua An, Li-Ming Song, Shi-Jie Zheng, Fan Zhang, Ai-Jun Dong , et al. (9 additional authors not shown)

Abstract: Spectral lag of the low-energy photons with respect to the high-energy ones is a common astrophysical phenomenon (such as Gamma-ray bursts and the Crab pulsar) and may serve as a key probe to the underlying radiation mechanism. However, spectral lag in keV range of the magnetar bursts has not been systematically studied yet. In this work, we perform a detailed spectral lag analysis with the Li-CCF… ▽ More Spectral lag of the low-energy photons with respect to the high-energy ones is a common astrophysical phenomenon (such as Gamma-ray bursts and the Crab pulsar) and may serve as a key probe to the underlying radiation mechanism. However, spectral lag in keV range of the magnetar bursts has not been systematically studied yet. In this work, we perform a detailed spectral lag analysis with the Li-CCF method for SGR J1935+2154 bursts observed by {\it Insight}-HXMT, GECAM and Fermi/GBM from July 2014 to Jan 2022. We discover that the spectral lags of about 61\% (non-zero significance >1$σ$) bursts from SGR J1935+2154 are linearly dependent on the photon energy ($E$) with $t_{\rm lag}(E)=α(E/{\rm keV})+C$, which may be explained by a linear change of the temperature of the blackbody-emitting plasma with time. The distribution of the slope ($α$) approximately follows a Gaussian function with mean and standard deviation of 0.02 ms/keV (i.e. high-energy photons arrive earlier) and 0.02 ms/keV, respectively. We also find that the distribution can be well fitted with three Gaussians with mean values of $\sim$ -0.009, 0.013 and 0.039 ms/keV, which may correspond to different origins of the bursts. These spectral lag features may have important implications on the magnetar bursts. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: accepted for publication in MNRAS

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

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

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

Submitted 18 March, 2023; originally announced March 2023.

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

The Large Sky Area Multi-object Fiber Spectroscopic Telescope (LAMOST) Quasar Survey: Quasar Properties from Data Release Six to Nine

Authors: Jun-Jie Jin, Xue-Bing Wu, Yuming Fu, Su Yao, Yan-Li Ai, Xiao-Tong Feng, Zi-Qi He, Qin-Chun Ma, Yu-Xuan Pang, Rui Zhu, Yan-xia Zhang, Hai-long Yuan, Zhi-ying Huo

Abstract: We report the fourth installment in the series of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Quasar Survey, which includes quasars observed between September, 2017 and June, 2021. There are in total 13,066 quasars reliably identified, of which 6,685 are newly discovered that are not reported in the SDSS DR14 quasar catalog or Million Quasars catalog. Because LAMOST does… ▽ More We report the fourth installment in the series of the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) Quasar Survey, which includes quasars observed between September, 2017 and June, 2021. There are in total 13,066 quasars reliably identified, of which 6,685 are newly discovered that are not reported in the SDSS DR14 quasar catalog or Million Quasars catalog. Because LAMOST does not provide accurate absolute flux calibration, we re-calibrate the spectra with the SDSS/Pan-STARRS1 multi-band photometric data. The emission line properties of H$α$, H$β$, Mg\,{\sc ii} and C\,{\sc iv}, and the continuum luminosities are measured by fitting the re-calibrated spectra. We also estimate the single-epoch virial black hole masses ($\rm M_{BH}$) using the derived emission line and continuum parameters. This is the first time that the emission line and continuum fluxes were estimated based on LAMOST re-calibrated quasar spectra. The catalog and spectra for these quasars are available online. After the nine-year LAMOST quasar survey, there are in total 56,175 identified quasars, of which 24,127 are newly discovered. The LAMOST quasar survey not only discovers a great number of new quasars, but also provides a database for investigating the spectral variability of the quasars observed by both LAMOST and SDSS, and finding rare quasars including changing-look quasars and broad absorption line quasars. △ Less

Submitted 29 December, 2022; v1 submitted 25 December, 2022; originally announced December 2022.

Comments: arXiv admin note: text overlap with arXiv:1811.01570, Accepted by ApJS

Search for Quasi-Periodical Oscillations in Precursors of Short and Long Gamma Ray Bursts

Authors: Shuo Xiao, Wen-Xi Peng, Shuang-Nan Zhang, Shao-Lin Xiong, Xiao-Bo Li, You-Li Tuo, He Gao, Yue Wang, Wang-Chen Xue, Chao Zheng, Yan-Qiu Zhang, Jia-Cong Liu, Cheng-Kui Li, Shu-Xu Yi, Xi-Lu Wang, Zhen Zhang, Ce Cai, Ai-Jun Dong, Wei Xie, Jian-Chao Feng, Qing-Bo Ma, De-Hua Wang, Xi-Hong Luo, Qi-Jun Zhi, Li-Ming Song , et al. (1 additional authors not shown)

Abstract: The precursors of short and long Gamma Ray Bursts (SGRBs and LGRBs) can serve as probes of their progenitors, as well as shedding light on the physical processes of mergers or core-collapse supernovae. Some models predict the possible existence of Quasi-Periodically Oscillations (QPO) in the precursors of SGRBs. Although many previous studies have performed QPO search in the main emission of SGRBs… ▽ More The precursors of short and long Gamma Ray Bursts (SGRBs and LGRBs) can serve as probes of their progenitors, as well as shedding light on the physical processes of mergers or core-collapse supernovae. Some models predict the possible existence of Quasi-Periodically Oscillations (QPO) in the precursors of SGRBs. Although many previous studies have performed QPO search in the main emission of SGRBs and LGRBs, so far there was no systematic QPO search in their precursors. In this work, we perform a detailed QPO search in the precursors of SGRBs and LGRBs detected by Fermi/GBM from 2008 to 2019 using the power density spectrum (PDS) in frequency domain and Gaussian processes (GP) in time domain. We do not find any convinced QPO signal with significance above 3 $σ$, possibly due to the low fluxes of precursors. Finally, the PDS continuum properties of both the precursors and main emissions are also studied for the first time, and no significant difference is found in the distributions of the PDS slope for precursors and main emissions in both SGRBs and LGRBs. △ Less

Submitted 16 October, 2022; originally announced October 2022.

Comments: submitted

Model Independent Approach of the JUNO $^8$B Solar Neutrino Program

Authors: JUNO Collaboration, Jie Zhao, Baobiao Yue, Haoqi Lu, Yufeng Li, Jiajie Ling, Zeyuan Yu, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Tsagkarakis Alexandros, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Weidong Bai , et al. (579 additional authors not shown)

Abstract: The physics potential of detecting $^8$B solar neutrinos will be exploited at the Jiangmen Underground Neutrino Observatory (JUNO), in a model independent manner by using three distinct channels of the charged-current (CC), neutral-current (NC) and elastic scattering (ES) interactions. Due to the largest-ever mass of $^{13}$C nuclei in the liquid-scintillator detectors and the {expected} low backg… ▽ More The physics potential of detecting $^8$B solar neutrinos will be exploited at the Jiangmen Underground Neutrino Observatory (JUNO), in a model independent manner by using three distinct channels of the charged-current (CC), neutral-current (NC) and elastic scattering (ES) interactions. Due to the largest-ever mass of $^{13}$C nuclei in the liquid-scintillator detectors and the {expected} low background level, $^8$B solar neutrinos would be observable in the CC and NC interactions on $^{13}$C for the first time. By virtue of optimized event selections and muon veto strategies, backgrounds from the accidental coincidence, muon-induced isotopes, and external backgrounds can be greatly suppressed. Excellent signal-to-background ratios can be achieved in the CC, NC and ES channels to guarantee the $^8$B solar neutrino observation. From the sensitivity studies performed in this work, we show that JUNO, with ten years of data, can reach the {1$σ$} precision levels of 5%, 8% and 20% for the $^8$B neutrino flux, $\sin^2θ_{12}$, and $Δm^2_{21}$, respectively. It would be unique and helpful to probe the details of both solar physics and neutrino physics. In addition, when combined with SNO, the world-best precision of 3% is expected for the $^8$B neutrino flux measurement. △ Less

Submitted 6 March, 2024; v1 submitted 15 October, 2022; originally announced October 2022.

Comments: 19 pages, 7 figures, accepted version to appear in The Astrophysical Journal. Yufeng Li and Jiajie Ling are corresponding authors

Journal ref: Astrophysical Journal 965 (2024) 122

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

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

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

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

Comments: 18 pages, 11 figures

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

Finding Quasars behind the Galactic Plane. II. Spectroscopic Identifications of 204 Quasars at $|b|< 20°$

Authors: Yuming Fu, Xue-Bing Wu, Linhua Jiang, Yanxia Zhang, Zhi-Ying Huo, Y. L. Ai, Qian Yang, Qinchun Ma, Xiaotong Feng, Ravi Joshi, Wei Jeat Hon, Christian Wolf, Jiang-Tao Li, Junjie Jin, Su Yao, Yuxuan Pang, Jian-Guo Wang, Kai-Xing Lu, Chuan-Jun Wang, Jie Zheng, Liang Xu, Xiao-Guang Yu, Bao-Li Lun, Pei Zuo

Abstract: Quasars behind the Galactic plane (GPQs) are important astrometric references and valuable probes of Galactic gas, yet the search for GPQs is difficult due to severe extinction and source crowding in the Galactic plane. In this paper, we present a sample of 204 spectroscopically confirmed GPQs at |b|<20°, 191 of which are new discoveries. This GPQ sample covers a wide redshift range from 0.069 to… ▽ More Quasars behind the Galactic plane (GPQs) are important astrometric references and valuable probes of Galactic gas, yet the search for GPQs is difficult due to severe extinction and source crowding in the Galactic plane. In this paper, we present a sample of 204 spectroscopically confirmed GPQs at |b|<20°, 191 of which are new discoveries. This GPQ sample covers a wide redshift range from 0.069 to 4.487. For the subset of 230 observed GPQ candidates, the lower limit of the purity of quasars is 85.2%, and the lower limit of the fraction of stellar contaminants is 6.1%. Using a multicomponent spectral fitting, we measure the emission line and continuum flux of the GPQs, and estimate their single-epoch virial black hole masses. Due to selection effects raised from Galactic extinction and target magnitude, these GPQs have higher black hole masses and continuum luminosities in comparison to the SDSS DR7 quasar sample. The spectral-fitting results and black hole mass estimates are compiled into a main spectral catalog, and an extended spectral catalog of GPQs. The successful identifications prove the reliability of both our GPQ selection methods and the GPQ candidate catalog, shedding light on the astrometric and astrophysical programs that make use of a large sample of GPQs in the future. △ Less

Submitted 29 July, 2022; v1 submitted 12 June, 2022; originally announced June 2022.

Comments: 17 pages, 10 figures, published in ApJS. Tables and spectra in this paper are available at https://doi.org/10.12149/101095

Journal ref: 2022 ApJS 261 32

HI Vertical Structure of Nearby Edge-on Galaxies from CHANG-ES

Authors: Yun Zheng, Jing Wang, Judith Irwin, Q. Daniel Wang, Jiangtao Li, Jayanne English, Qingchuan Ma, Ran Wang, Ke Wang, Marita Krause, Toky H. Randriamampandry, Rainer Beck

Abstract: We study the vertical distribution of the highly inclined galaxies from the Continuum Halos in Nearby Galaxies - an EVLA Survey (CHANG-ES). We explore the feasibility of photometrically deriving the HI disk scale-heights from the moment-0 images of the relatively edge-on galaxies with inclination >80 deg, by quantifying the systematic broadening effects and thus deriving correction equations for d… ▽ More We study the vertical distribution of the highly inclined galaxies from the Continuum Halos in Nearby Galaxies - an EVLA Survey (CHANG-ES). We explore the feasibility of photometrically deriving the HI disk scale-heights from the moment-0 images of the relatively edge-on galaxies with inclination >80 deg, by quantifying the systematic broadening effects and thus deriving correction equations for direct measurements. The corrected HI disk scale-heights of the relatively edge-on galaxies from the CHANG-ES sample show trends consistent with the quasi-equilibrium model of the vertical structure of gas disks. The procedure provide a convenient way to derive the scale-heights and can easily be applied to statistical samples in the future. △ Less

Submitted 20 May, 2022; originally announced May 2022.

Comments: Accepted for publication in Research in Astronomy nd Astrophysics (RAA) journal

Prospects for Detecting the Diffuse Supernova Neutrino Background with JUNO

Authors: JUNO Collaboration, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato, Antonio Bergnoli, Thilo Birkenfeld, Sylvie Blin , et al. (577 additional authors not shown)

Abstract: We present the detection potential for the diffuse supernova neutrino background (DSNB) at the Jiangmen Underground Neutrino Observatory (JUNO), using the inverse-beta-decay (IBD) detection channel on free protons. We employ the latest information on the DSNB flux predictions, and investigate in detail the background and its reduction for the DSNB search at JUNO. The atmospheric neutrino induced n… ▽ More We present the detection potential for the diffuse supernova neutrino background (DSNB) at the Jiangmen Underground Neutrino Observatory (JUNO), using the inverse-beta-decay (IBD) detection channel on free protons. We employ the latest information on the DSNB flux predictions, and investigate in detail the background and its reduction for the DSNB search at JUNO. The atmospheric neutrino induced neutral current (NC) background turns out to be the most critical background, whose uncertainty is carefully evaluated from both the spread of model predictions and an envisaged \textit{in situ} measurement. We also make a careful study on the background suppression with the pulse shape discrimination (PSD) and triple coincidence (TC) cuts. With latest DSNB signal predictions, more realistic background evaluation and PSD efficiency optimization, and additional TC cut, JUNO can reach the significance of 3$σ$ for 3 years of data taking, and achieve better than 5$σ$ after 10 years for a reference DSNB model. In the pessimistic scenario of non-observation, JUNO would strongly improve the limits and exclude a significant region of the model parameter space. △ Less

Submitted 13 October, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 29 pages, 11 figures, final published version in JCAP

Journal ref: JCAP 10 (2022) 033

A CRASH simulation of the contribution of binary stars to the epoch of reionization

Authors: Qing-Bo Ma, Silvia Fiaschi, Benedetta Ciardi, Philipp Busch, Marius B. Eide

Abstract: We use a set of 3D radiative transfer simulations to study the effect that a large fraction of binary stars in galaxies during the epoch of reionization has on the physical properties of the intergalactic medium (i.e. the gas temperature and the ionization state of hydrogen and helium), on the topology of the ionized bubbles and on the 21 cm power spectra. Consistently to previous literature, we f… ▽ More We use a set of 3D radiative transfer simulations to study the effect that a large fraction of binary stars in galaxies during the epoch of reionization has on the physical properties of the intergalactic medium (i.e. the gas temperature and the ionization state of hydrogen and helium), on the topology of the ionized bubbles and on the 21 cm power spectra. Consistently to previous literature, we find that the inclusion of binary stars can speed up the reionization process of HI and HeI, while HeII reionization is still dominated by more energetic sources, especially accreting black holes. The earlier ionization attained with binary stars allows for more time for cooling and recombination, so that gas fully ionized by binary stars is typically colder than that ionized by single stars at any given redshift. With the same volume averaged ionization fraction, the inclusion of binary stars results in fewer small ionized bubbles and more large ones, with visible effects also on the large scales of the 21 cm power spectrum. △ Less

Submitted 8 April, 2022; originally announced April 2022.

Comments: 14 pages, 11 figures, MNRAS accepted

CHANG-ES XXV: HI Imaging of Nearby Edge-on Galaxies -- Data Release 4

Authors: Yun Zheng, Jing Wang, Judith Irwin, Jayanne English, Qingchuan Ma, Ran Wang, Ke Wang, Q. Daniel Wang, Marita Krause, Toky H. Randriamampandry, Jiangtao Li, Rainer Beck

Abstract: We present the HI distribution of galaxies from the Continuum Halos in Nearby Galaxies - an EVLA Survey (CHANG-ES). Though the observational mode was not optimized for detecting HI, we successfully produce HI cubes for 19 galaxies. The moment-0 maps from this work are available on CHANG-ES data release website, i.e., https://www.queensu.ca/changes. Our sample is dominated by star-forming, HI-rich… ▽ More We present the HI distribution of galaxies from the Continuum Halos in Nearby Galaxies - an EVLA Survey (CHANG-ES). Though the observational mode was not optimized for detecting HI, we successfully produce HI cubes for 19 galaxies. The moment-0 maps from this work are available on CHANG-ES data release website, i.e., https://www.queensu.ca/changes. Our sample is dominated by star-forming, HI-rich galaxies at distances from 6.27 to 34.1 Mpc. HI interferometric images on two of these galaxies (NGC 5792 and UGC 10288) are presented here for the first time, while 12 of our remaining sample galaxies now have better HI spatial resolutions and/or sensitivities of intensity maps than those in existing publications. We characterize the average scale heights of the HI distributions for a subset of most inclined galaxies (inclination > 80 deg), and compare them to the radio continuum intensity scale heights, which have been derived in a similar way. The two types of scale heights are well correlated, with similar dependence on disk radial extension and star formation rate surface density but different dependence on mass surface density. This result indicates that the vertical distribution of the two components may be governed by similar fundamental physics but with subtle differences. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: 27 pages, 10 figures, accepted for publication in MNRAS

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

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

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

Submitted 11 November, 2021; originally announced November 2021.

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

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

Hint of a truncated primordial spectrum from the CMB large-scale anomalies

Authors: Fulvio Melia, Qingbo Ma, Jun-Jie Wei, Bo Yu

Abstract: Several satellite missions have uncovered a series of potential anomalies in the fluctuation spectrum of the cosmic microwave background temperature, including: (1) an unexpectedly low level of correlation at large angles, manifested via the angular correlation function, C(theta); and (2) missing power in the low multipole moments of the angular power spectrum, C_ell. Their origin is still debated… ▽ More Several satellite missions have uncovered a series of potential anomalies in the fluctuation spectrum of the cosmic microwave background temperature, including: (1) an unexpectedly low level of correlation at large angles, manifested via the angular correlation function, C(theta); and (2) missing power in the low multipole moments of the angular power spectrum, C_ell. Their origin is still debated, however, due to a persistent lack of clarity concerning the seeding of quantum fluctuations in the early Universe. A likely explanation for the first of these appears to be a cutoff, k_min=(3.14 +/- 0.36) x 10^{-4} Mpc^{-1}, in the primordial power spectrum, P(k). Our goal in this paper is twofold: (1) we examine whether the same k_min can also self-consistently explain the missing power at large angles, and (2) we confirm that the of this cutoff in P(k) does not adversely affect the remarkable consistency between the prediction of Planck-LCDM and the Planck measurements at ell > 30. We use the publicly available code CAMB to calculate the angular power spectrum, based on a line-of-sight approach. The code is modified slightly to include the additional parameter (i.e., k_min) characterizing the primordial power spectrum. In addition to this cutoff, the code optimizes all of the usual standard-model parameters. In fitting the angular power spectrum, we find an optimized cutoff, k_min = 2.04^{+1.4}_{-0.79} x 10^{-4} Mpc^{-1}, when using the whole range of ell's, and k_min=3.3^{+1.7}_{-1.3} x 10^{-4} Mpc^{-1}, when fitting only the range ell < 30, where the Sachs-Wolfe effect is dominant. These are fully consistent with the value inferred from C(theta), suggesting that both of these large-angle anomalies may be due to the same truncation in P(k). △ Less

Submitted 12 September, 2021; originally announced September 2021.

Comments: 9 pages, 4 figures, 1 table. Accepted for publication in A&A

Journal ref: A&A 655, A70 (2021)

Accretion Disk Sizes from Continuum Reverberation Mapping of AGN Selected from the ZTF Survey

Authors: Vivek Kumar Jha, Ravi Joshi, Hum Chand, Xue-Bing Wu, Luis C. Ho, Shantanu Rastogi, Qinchun Ma

Abstract: We present the accretion disk size estimates for a sample of 19 active galactic nuclei (AGN) using the optical $g$, $r$, and $i$ band light curves obtained from the Zwicky Transient Facility (ZTF) survey. All the AGN have reliable supermassive black hole (SMBH) mass estimates based on previous reverberation mapping measurements. The multi-band light curves are cross-correlated, and the reverberati… ▽ More We present the accretion disk size estimates for a sample of 19 active galactic nuclei (AGN) using the optical $g$, $r$, and $i$ band light curves obtained from the Zwicky Transient Facility (ZTF) survey. All the AGN have reliable supermassive black hole (SMBH) mass estimates based on previous reverberation mapping measurements. The multi-band light curves are cross-correlated, and the reverberation lag is estimated using the Interpolated Cross-Correlation Function (ICCF) method and the Bayesian method using the {\sc javelin} code. As expected from the disk reprocessing arguments, the $g-r$ band lags are shorter than the $g-i$ band lags for this sample. The interband lags for all, but 5 sources, are larger than the sizes predicted from the standard Shakura Sunyaev (SS) analytical model. We fit the light curves directly using a thin disk model implemented through the {\sc javelin} code to get the accretion disk sizes. The disk sizes obtained using this model are on an average 3.9 times larger than the prediction based on the SS disk model. We find a weak correlation between the disk sizes and the known physical parameters, namely, the luminosity and the SMBH mass. In the near future, a large sample of AGN covering broader ranges of luminosity and SMBH mass from large photometric surveys would be helpful in a better understanding of the structure and physics of the accretion disk. △ Less

Submitted 12 January, 2022; v1 submitted 10 September, 2021; originally announced September 2021.

Comments: 12 pages, 9 Figures. Accepted for publication in MNRAS. This version matches with the accepted version

Sensitivity tests of cosmic velocity fields to massive neutrinos

Authors: Shuren Zhou, Zhenjie Liu, Qinglin Ma, Yu Liu, Le Zhang, Xiao-Dong Li, Yang Wang, Xin Wang, Yu Yu, Hao-Ran Yu, Yi Zheng

Abstract: We investigate impacts of massive neutrinos on the cosmic velocity fields, employing high-resolution cosmological $N$-body simulations provided by the information-optimized CUBE code, where cosmic neutrinos are evolved using collisionless hydrodynamics and their perturbations can be accurately resolved. In this study we focus, for the first time, on the analysis of massive-neutrino induced suppres… ▽ More We investigate impacts of massive neutrinos on the cosmic velocity fields, employing high-resolution cosmological $N$-body simulations provided by the information-optimized CUBE code, where cosmic neutrinos are evolved using collisionless hydrodynamics and their perturbations can be accurately resolved. In this study we focus, for the first time, on the analysis of massive-neutrino induced suppression effects in various cosmic velocity field components of velocity magnitude, divergence, vorticity and dispersion. By varying the neutrino mass sum $M_ν$ from 0 -- 0.4 eV, the simulations show that, the power spectra of vorticity -- exclusively sourced by non-linear structure formation that is affected by massive neutrinos significantly -- is very sensitive to the mass sum, which potentially provide novel signatures in detecting massive neutrinos. Furthermore, using the chi-square statistic, we quantitatively test the sensitivity of the density and velocity power spectra to the neutrino mass sum. Indeed, we find that, the vorticity spectrum has the highest sensitivity, and the null hypothesis of massless neutrinos is incompatible with both vorticity and divergence spectra from $M_ν=0.1$ eV at high significance ($p$-value $= 0.03$ and $0.07$, respectively). These results demonstrate clearly the importance of peculiar velocity field measurements, in particular of vorticity and divergence components, in determination of neutrino mass and mass hierarchy. △ Less

Submitted 5 March, 2022; v1 submitted 28 August, 2021; originally announced August 2021.

Comments: 13 pages, 9 figures. Monthly Notices of the Royal Astronomical Society, 2022

Predicting exoplanet mass from radius and incident flux: A Bayesian mixture model

Authors: Qi Ma, Sujit K Ghosh

Abstract: The relationship between mass and radius (M-R relation) is the key for inferring the planetary compositions and thus valuable for the studies of formation and migration models. However, the M-R relation alone is not enough for planetary characterization due to the dependence of it on other confounding variables. This paper provides a non-trivial extension of the M-R relation by including the incid… ▽ More The relationship between mass and radius (M-R relation) is the key for inferring the planetary compositions and thus valuable for the studies of formation and migration models. However, the M-R relation alone is not enough for planetary characterization due to the dependence of it on other confounding variables. This paper provides a non-trivial extension of the M-R relation by including the incident flux as an additional variable. By using Bayesian hierarchical modeling (BHM) that leverages the flexibility of finite mixture models, a probabilistic mass-radius-flux relationship (M-R-F relation) is obtained based on a sample of 319 exoplanets. We find that the flux has nonnegligible impact on the M-R relation, while such impact is strongest for hot-Jupiters. On the population level, the planets with higher level of flux tend to be denser, and high flux could trigger significant mass loss for plants with radii larger than $13R_{\oplus}$. As a result, failing to account for the flux in mass prediction would cause systematic over or under-estimation. With the recent advent of computing power, although a lot of complex statistical models can be fitted using Monte Carlo methods, it has largely remain illusive how to validate these complex models when the data are observed with large measurement errors. We present two novel methods to examine model assumptions, which can be used not only for the models we present in this paper but can also be adapted for other statistical models. △ Less

Submitted 2 June, 2021; originally announced June 2021.

Comments: Accepted for publication in MNRAS

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

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

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

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

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

Investigating X-ray sources during the epoch of reionization with the 21 cm signal

Authors: Qing-Bo Ma, Benedetta Ciardi, Marius B. Eide, Philipp Busch, Yi Mao, Qi-Jun Zhi

Abstract: Heating of neutral gas by energetic sources is crucial for the prediction of the 21 cm signal during the epoch of reionization (EoR). To investigate differences induced on statistics of the 21 cm signal by various source types, we use five radiative transfer simulations which have the same stellar UV emission model and varying combinations of more energetic sources, such as X-ray binaries (XRBs),… ▽ More Heating of neutral gas by energetic sources is crucial for the prediction of the 21 cm signal during the epoch of reionization (EoR). To investigate differences induced on statistics of the 21 cm signal by various source types, we use five radiative transfer simulations which have the same stellar UV emission model and varying combinations of more energetic sources, such as X-ray binaries (XRBs), accreting nuclear black holes (BHs) and hot interstellar medium emission (ISM). We find that the efficient heating from the ISM increases the average global 21~cm signal, while reducing its fluctuations and thus power spectrum. A clear impact is also observed in the bispectrum in terms of scale and timing of the transition between a positive and a negative value. The impact of XRBs is similar to that of the ISM, although it is delayed in time and reduced in intensity because of the less efficient heating. Due to the paucity of nuclear BHs, the behaviour of the 21~cm statistics in their presence is very similar to that of a case when only stars are considered, with the exception of the latest stages of reionization, when the effect of BHs is clearly visible. We find that differences between the source scenarios investigated here are larger than the instrumental noise of SKA1-low at $z \gtrsim 7-8$, suggesting that in the future it might be possible to constrain the spectral energy distribution of the sources contributing to the reionization process. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 22pages,10 figures, apj accepted

Finding Quasars behind the Galactic Plane. I. Candidate Selections with Transfer Learning

Authors: Yuming Fu, Xue-Bing Wu, Qian Yang, Anthony G. A. Brown, Xiaotong Feng, Qinchun Ma, Shuyan Li

Abstract: Quasars behind the Galactic plane (GPQs) are important astrometric references and useful probes of Milky Way gas. However, the search for GPQs is difficult due to large extinctions and high source densities in the Galactic plane. Existing selection methods for quasars developed using high Galactic latitude (high-$b$) data cannot be applied to the Galactic plane directly because the photometric dat… ▽ More Quasars behind the Galactic plane (GPQs) are important astrometric references and useful probes of Milky Way gas. However, the search for GPQs is difficult due to large extinctions and high source densities in the Galactic plane. Existing selection methods for quasars developed using high Galactic latitude (high-$b$) data cannot be applied to the Galactic plane directly because the photometric data obtained from high-$b$ regions and the Galactic plane follow different probability distributions. To alleviate this dataset shift problem for quasar candidate selection, we adopt a Transfer Learning Framework at both data and algorithm levels. At the data level, to make a training set in which dataset shift is modeled, we synthesize quasars and galaxies behind the Galactic plane based on SDSS sources and Galactic dust map. At the algorithm level, to reduce the effect of class imbalance, we transform the three-class classification problem for stars, galaxies, and quasars to two binary classification tasks. We apply XGBoost algorithm on Pan-STARRS1 (PS1) and AllWISE photometry for classification, and additional cut on Gaia proper motion to remove stellar contaminants. We obtain a reliable GPQ candidate catalog with 160,946 sources located at $|b|\leq 20^{\circ}$ in PS1-AllWISE footprint. Photometric redshifts of GPQ candidates achieved with XGBoost regression algorithm show that our selection method can identify quasars in a wide redshift range ($0<z\lesssim5$). This study extends the systematic searches for quasars to the dense stellar fields and shows the feasibility of using astronomical knowledge to improve data mining under complex conditions in the Big Data era. △ Less

Submitted 19 February, 2021; originally announced February 2021.

Comments: 28 pages, 15 figures, accepted for publication in ApJS. Table 5 can be downloaded at http://paperdata.china-vo.org/FuYuming/2021/ApJS/GPQ_One/table-05.fits

Journal ref: 2021 ApJS 254 6

Construction and On-site Performance of the LHAASO WFCTA Camera

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

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

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

Comments: 45 pages, 21 figures, article

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

Physical and chemical structure of the Serpens filament -- fast formation and gravity-driven accretion

Authors: Y. Gong, A. Belloche, F. J. Du, K. M. Menten, C. Henkel, G. X. Li, F. Wyrowski, R. Q. Mao

Abstract: The Serpens filament, a prominent elongated structure in a relatively nearby molecular cloud, is believed to be at an early evolutionary stage, so studying its physical and chemical properties can shed light on filament formation and early evolution. The main goal is to address the physical and chemical properties as well as the dynamical state of the Serpens filament at a spatial resolution of… ▽ More The Serpens filament, a prominent elongated structure in a relatively nearby molecular cloud, is believed to be at an early evolutionary stage, so studying its physical and chemical properties can shed light on filament formation and early evolution. The main goal is to address the physical and chemical properties as well as the dynamical state of the Serpens filament at a spatial resolution of $\sim$0.07 pc and a spectral resolution of $\lesssim$0.1~km~s$^{-1}$. We performed $^{13}$CO (1--0), C$^{18}$O (1--0), C$^{17}$O (1--0), $^{13}$CO (2--1), C$^{18}$O (2--1), and C$^{17}$O (2--1) imaging observations toward the Serpens filament with the Institut de Radioastronomie Millim{é}trique 30-m (IRAM-30 m) and Atacama Pathfinder EXperiment (APEX) telescopes. Widespread narrow $^{13}$CO (2--1) self-absorption is observed in this filament, causing the $^{13}$CO morphology to be different from the filamentary structure traced by C$^{18}$O and C$^{17}$O. Our excitation analysis suggests that the opacities of C$^{18}$O transitions become higher than unity in most regions, and this analysis confirms the presence of widespread CO depletion. Further we show that the local velocity gradients have a tendency to be perpendicular to the filament's long axis in the outskirts and parallel to the large-scale magnetic field direction. The magnitudes of the local velocity gradients decrease toward the filament's crest. The observed velocity structure can be a result of gravity-driven accretion flows. The isochronic evolutionary track of the C$^{18}$O freeze-out process indicates the filament is young with an age of $\lesssim$2 Myr. We propose that the Serpens filament is a newly-formed slightly-supercritical structure which appears to be actively accreting material from its ambient gas. △ Less

Submitted 20 January, 2021; v1 submitted 22 December, 2020; originally announced December 2020.

Comments: 26 pages, 24 figures, 2 tables, accepted for publication in A&A; A 3D view of the C18O(2-1) cube can be found in https://gongyan2444.github.io/gallery.html

Journal ref: A&A 646, A170 (2021)

The observation of the Crab Nebula with LHAASO-KM2A for the performance study

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

Abstract: As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the… ▽ More As a sub-array of the Large High Altitude Air Shower Observatory (LHAASO), KM2A is mainly designed to cover a large fraction of the northern sky to hunt for gamma-ray sources at energies above 10 TeV. Even though the detector construction is still underway, a half of the KM2A array has been operating stably since the end of 2019. In this paper, we present the pipeline of KM2A data analysis and the first observation on the Crab Nebula, a standard candle in very high energy gamma-ray astronomy. We detect gamma-ray signals from the Crab Nebula in both energy ranges of 10$-$100 TeV and $>$100 TeV with high significance, by analyzing the KM2A data of 136 live days between December 2019 and May 2020. With the observations, we test the detector performance including angular resolution, pointing accuracy and cosmic ray background rejection power. The energy spectrum of the Crab Nebula in the energy range 10-250 TeV fits well with a single power-law function dN/dE =(1.13$\pm$0.05$_{stat}$$\pm$0.08$_{sys}$)$\times$10$^{-14}$$\cdot$(E/20TeV)$^{-3.09\pm0.06_{stat}\pm0.02_{sys}}$ cm$^{-2}$ s$^{-1}$ TeV$^{-1}$. It is consistent with previous measurements by other experiments. This opens a new window of gamma-ray astronomy above 0.1 PeV through which ultrahigh-energy gamma-ray new phenomena, such as cosmic PeVatrons, might be discovered. △ Less

Submitted 13 October, 2020; originally announced October 2020.

Comments: 13 pages, 15 figures,submitted to CPC

Journal ref: Chinese Phys. C 45 025002 (2021)

Using the Mark Weighted Correlation Functions to Improve the Constraints on Cosmological Parameters

Authors: Yizhao Yang, Haitao Miao, Qinglin Ma, Miaoxin Liu, Cristiano G. Sabiu, Jaime Forero-Romero, Yuanzhu Huang, Limin Lai, Qiyue Qian, Yi Zheng, Xiao-Dong Li

Abstract: We used the mark weighted correlation functions (MCFs), $W(s)$, to study the large scale structure of the Universe. We studied five types of MCFs with the weighting scheme $ρ^α$, where $ρ$ is the local density, and $α$ is taken as $-1,\ -0.5,\ 0,\ 0.5$, and 1. We found that different MCFs have very different amplitudes and scale-dependence. Some of the MCFs exhibit distinctive peaks and valleys th… ▽ More We used the mark weighted correlation functions (MCFs), $W(s)$, to study the large scale structure of the Universe. We studied five types of MCFs with the weighting scheme $ρ^α$, where $ρ$ is the local density, and $α$ is taken as $-1,\ -0.5,\ 0,\ 0.5$, and 1. We found that different MCFs have very different amplitudes and scale-dependence. Some of the MCFs exhibit distinctive peaks and valleys that do not exist in the standard correlation functions. Their locations are robust against the redshifts and the background geometry, however it is unlikely that they can be used as ``standard rulers'' to probe the cosmic expansion history. Nonetheless we find that these features may be used to probe parameters related with the structure formation history, such as the values of $σ_8$ and the galaxy bias. Finally, after conducting a comprehensive analysis using the full shapes of the $W(s)$s and $W_{Δs}(μ)$s, we found that, combining different types of MCFs can significantly improve the cosmological parameter constraints. Compared with using only the standard correlation function, the combinations of MCFs with $α=0,\ 0.5,\ 1$ and $α=0,\ -1,\ -0.5,\ 0.5,\ 1$ can improve the constraints on $Ω_m$ and $w$ by $\approx30\%$ and $50\%$, respectively. We find highly significant evidence that MCFs can improve cosmological parameter constraints. △ Less

Submitted 1 August, 2020; v1 submitted 6 July, 2020; originally announced July 2020.

Comments: 15pages, 17figures, APJ accepted

Probing the high-z IGM with the hyperfine transition of $^3$He$^+$

Authors: Shivan Khullar, Qingbo Ma, Philipp Busch, Benedetta Ciardi, Marius B. Eide, Koki Kakiichi

Abstract: The hyperfine transition of $^3$He$^+$ at 3.5cm has been thought as a probe of the high-z IGM since it offers a unique insight into the evolution of the helium component of the gas, as well as potentially give an independent constraint on the 21cm signal from neutral hydrogen. In this paper, we use radiative transfer simulations of reionization driven by sources such as stars, X-ray binaries, accr… ▽ More The hyperfine transition of $^3$He$^+$ at 3.5cm has been thought as a probe of the high-z IGM since it offers a unique insight into the evolution of the helium component of the gas, as well as potentially give an independent constraint on the 21cm signal from neutral hydrogen. In this paper, we use radiative transfer simulations of reionization driven by sources such as stars, X-ray binaries, accreting black holes and shock heated interstellar medium, and simulations of a high-z quasar to characterize the signal and analyze its prospects of detection. We find that the peak of the signal lies in the range 1-50 $μ$K for both environments, but while around the quasar it is always in emission, in the case of cosmic reionization a brief period of absorption is expected. As the evolution of HeII is determined by stars, we find that it is not possible to distinguish reionization histories driven by more energetic sources. On the other hand, while a bright QSO produces a signal in 21cm that is very similar to the one from a large collection of galaxies, its signature in 3.5cm is very peculiar and could be a powerful probe to identify the presence of the QSO. We analyze the prospects of the signal's detectability using SKA1-mid as our reference telescope. We find that the noise power spectrum dominates over the power spectrum of the signal, although a modest S/N ratio can be obtained when the wavenumber bin width and the survey volume are sufficiently large. △ Less

Submitted 2 July, 2020; originally announced July 2020.

Comments: 10 pages, 13 figures, accepted for publication in MNRAS