Showing 1–50 of 171 results for author: Feng, H

Velocity-Resolved Ionization Mapping of Broad Line Region. I. Insights into Diverse Geometry and Kinematics

Authors: Sha-Sha Li, Hai-Cheng Feng, H. T. Liu, J. M. Bai, Xiang Ji, Cheng Cheng, Kai-Xing Lu, Jian-Guo Wang, Rui Li

Abstract: Broad emission lines of active galactic nuclei (AGNs) originate from the broad-line region (BLR), consisting of dense gas clouds in orbit around an accreting supermassive black hole. Understanding the geometry and kinematics of the region is crucial for gaining insights into the physics and evolution of AGNs. Conventional velocity-resolved reverberation mapping may face challenges in disentangling… ▽ More Broad emission lines of active galactic nuclei (AGNs) originate from the broad-line region (BLR), consisting of dense gas clouds in orbit around an accreting supermassive black hole. Understanding the geometry and kinematics of the region is crucial for gaining insights into the physics and evolution of AGNs. Conventional velocity-resolved reverberation mapping may face challenges in disentangling the degeneracy between intricate motion and geometry of this region. To address this challenge, new key constraints are required. Here, we report the discovery of an asymmetric BLR using a novel technique: velocity-resolved ionization mapping, which can map the distance of emitting gas clouds by measuring Hydrogen line ratios at different velocities. By analyzing spectroscopic monitoring data, we find that the Balmer decrement is anticorrelated with the continuum and correlated with the lags across broad emission line velocities. Some line ratio profiles deviate from the expectations for a symmetrically virialized BLR, suggesting that the red-shifted and blue-shifted gas clouds may not be equidistant from the supermassive black hole (SMBH). This asymmetric geometry might represent a formation imprint, provide new perspectives on the evolution of AGNs, and influence SMBH mass measurements. △ Less

Submitted 7 July, 2024; originally announced July 2024.

Comments: 20 pages, 10 figures, Accepted by ApJ

A timing view of the additional high-energy spectral component discovered in the black hole candidate Swift J1727.8-1613

Authors: Zi-Xu Yang, Liang Zhang, Shuang-Nan Zhang, L. Tao, Shu Zhang, Ruican Ma, Qingcui Bu, Yue Huang, He-Xin Liu, Wei Yu, Guang C. Xiao, Peng-Ju Wang, Hua Feng, Li-Ming Song, Xiang Ma, Mingyu Ge, QingChang Zhao, J. L. Qu

Abstract: We present an energy-dependent analysis for the type-C quasi-periodic oscillations (QPOs) observed in the black hole X-ray binary Swift J1727.8-1613 using Insight-HXMT observations. We find that the QPO fractional rms at energies above 40 keV is significantly higher than that below 20 keV. This is the first report of a high energy (HE)-rms excess in the rms spectrum of a black hole X-ray binary. I… ▽ More We present an energy-dependent analysis for the type-C quasi-periodic oscillations (QPOs) observed in the black hole X-ray binary Swift J1727.8-1613 using Insight-HXMT observations. We find that the QPO fractional rms at energies above 40 keV is significantly higher than that below 20 keV. This is the first report of a high energy (HE)-rms excess in the rms spectrum of a black hole X-ray binary. In the high energy band, an extra hard component is observed in additional to the standard thermal Comptonization component at similar energy band. The value of the QPO HE-rms excess is not only correlated with the disk parameters and the photon index of the standard Comptonization component, but also exhibits a moderate positive correlation with the flux of the additional hard spectral component. No features in the QPO phase-lag spectra are seen corresponding to the additional hard component. We propose that the additional hard component in the spectrum may originate from jet emission and the associated QPO HE-rms excess can be explained by the precession of the jet base. △ Less

Submitted 6 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

M17 MIR: A Massive Star is Forming via Episodic Mass Accretion

Authors: Wei Zhou, Zhiwei Chen, Zhibo Jiang, Haoran Feng, Yu Jiang

Abstract: We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) band 6 data for the outbursting massive protostar M17~MIR. The ALMA CO $J=2-1$ data reveal a collimated and bipolar north-south outflow from M17~MIR. The blue-shifted outflow exhibits four CO knots (N1 to N4) along the outflow axis, while the red-shifted outflow appears as a single knot (S1). The extremely high velocity (EHV) emis… ▽ More We analyzed the Atacama Large Millimeter/submillimeter Array (ALMA) band 6 data for the outbursting massive protostar M17~MIR. The ALMA CO $J=2-1$ data reveal a collimated and bipolar north-south outflow from M17~MIR. The blue-shifted outflow exhibits four CO knots (N1 to N4) along the outflow axis, while the red-shifted outflow appears as a single knot (S1). The extremely high velocity (EHV) emissions of N1 and S1 are jet-like and contain sub-knots along the outflow axis. Assuming the nearest EHV sub-knots trace the ejecta from the accretion outbursts in the past decades, a tangential ejection velocity of $\sim421\,\mathrm{km\,s^{-1}}$ is derived for M17~MIR. Assuming the same velocity, the dynamical times of the multiple ejecta, traced by the four blue-shifted CO knots, range from 20 to 364 years. The four blue-shifted CO knots imply four clustered accretion outbursts with a duration of tens of years in the past few hundred years. The intervals between the four clustered accretion outbursts are also about tens of years. These properties of the four clustered accretion outbursts are in line with the disk gravitational instability and fragmentation model. The episodic accretion history of M17~MIR traced by episodic outflow suggests that a massive star can form from a lower-mass protostar via frequent episodic accretion events triggered by disk gravitational instability and fragmentation. The first detection of the knotty outflow from an outbursting massive protostar suggests that mass ejections accompanied with accretion events could serve as an effective diagnostic tool for the episodic accretion histories of massive protostars. △ Less

Submitted 17 June, 2024; v1 submitted 7 June, 2024; originally announced June 2024.

Comments: Accepted for publication in ApJL; typos corrected

Morpho-Photometric Classification of KiDS DR5 Sources Based on Neural Networks: A Comprehensive Star-Quasar-Galaxy Catalog

Authors: Hai-Cheng Feng, Rui Li, Nicola R. Napolitano, Sha-Sha Li, J. M. Bai, Ran Li, H. T. Liu, Kai-Xing Lu, Mario Radovich, Huan-Yuan Shan, Jian-Guo Wang, Wen-Zhe Xi, Ling-Hua Xie, Yang-Wei Zhang

Abstract: We present a novel multimodal neural network for classifying astronomical sources in multiband ground-based observations, from optical to near infrared, to separate sources in stars, galaxies and quasars. Our approach combines a convolutional neural network branch for learning morphological features from $r$-band images with an artificial neural network branch for extracting spectral energy distri… ▽ More We present a novel multimodal neural network for classifying astronomical sources in multiband ground-based observations, from optical to near infrared, to separate sources in stars, galaxies and quasars. Our approach combines a convolutional neural network branch for learning morphological features from $r$-band images with an artificial neural network branch for extracting spectral energy distribution (SED) information. Specifically, we have used 9-band optical ($ugri$) and NIR ($ZYHJK_s$) data from the Kilo-Degree Survey (KiDS) Data Release 5. The two branches of the network are concatenated and feed into fully-connected layers for final classification. We train the network on a spectroscopically confirmed sample from the Sloan Digital Sky Survey cross-matched with KiDS. The trained model achieves 98.76\% overall accuracy on an independent testing dataset, with F1 scores exceeding 95\% for each class. Raising the output probability threshold, we obtain higher purity at the cost of a lower completeness. We have also validated the network using external catalogs cross-matched with KiDS, correctly classifying 99.74\% of a pure star sample selected from Gaia parallaxes and proper motions, and 99.74\% of an external galaxy sample from the Galaxy and Mass Assembly survey, adjusted for low-redshift contamination. We apply the trained network to 27,334,751 KiDS DR5 sources with $r \leqslant 23$ mag to generate a new classification catalog. This multimodal neural network successfully leverages both morphological and SED information to enable efficient and robust classification of stars, quasars, and galaxies in large photometric surveys. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 18 pages, 12 figures, 2 tables, Submitted to ApJS

All-sky Guide Star Catalog for CSST

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

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

Submitted 3 June, 2024; originally announced June 2024.

Comments: published on RAA

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

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

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

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

Comments: 15 pages, 9 figures

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

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

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

Submitted 19 May, 2024; originally announced May 2024.

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

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

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

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

Submitted 13 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

Polarization Perspectives on Hercules X-1: Further Constraining the Geometry

Authors: Qingchang Zhao, Hancheng Li, Lian Tao, Hua Feng, Shuangnan Zhang, Roland Walter, Mingyu Ge, Hao Tong, Long Ji, Liang Zhang, Jinlu Qu, Yue Huang, Xiang Ma, Shu Zhang, Qianqing Yin, Hongxing Yin, Ruican Ma, Shujie Zhao, Panping Li, Zixu Yang, Hexin Liu, Wei Yu, Yiming Huang, Zexi Li, Yajun Li , et al. (2 additional authors not shown)

Abstract: We conduct a comprehensive analysis of the accreting X-ray pulsar, Hercules X-1, utilizing data from IXPE and NuSTAR. IXPE performed five observations of Her X-1, consisting of three in the Main-on state and two in the Short-on state. Our time-resolved analysis uncovers the linear correlations between the flux and polarization degree as well as the pulse fraction and polarization degree. Geometry… ▽ More We conduct a comprehensive analysis of the accreting X-ray pulsar, Hercules X-1, utilizing data from IXPE and NuSTAR. IXPE performed five observations of Her X-1, consisting of three in the Main-on state and two in the Short-on state. Our time-resolved analysis uncovers the linear correlations between the flux and polarization degree as well as the pulse fraction and polarization degree. Geometry parameters are rigorously constrained by fitting the phase-resolved modulations of Cyclotron Resonance Scattering Feature and polarization angle with a simple dipole model and Rotating Vector Model respectively, yielding roughly consistent results. The changes of $χ_{\rm p}$ (the position angle of the pulsar's spin axis on the plane of the sky) between different Main-on observations suggest the possible forced precession of the neutron star crust. Furthermore, a linear association between the energy of Cyclotron Resonance Scattering Feature and polarization angle implies the prevalence of a dominant dipole magnetic field, and their phase-resolved modulations likely arise from viewing angle effects. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: Accepted for MNRAS

Soft X-ray prompt emission from a high-redshift gamma-ray burst EP240315a

Authors: Y. Liu, H. Sun, D. Xu, D. S. Svinkin, J. Delaunay, N. R. Tanvir, H. Gao, C. Zhang, Y. Chen, X. -F. Wu, B. Zhang, W. Yuan, J. An, G. Bruni, D. D. Frederiks, G. Ghirlanda, J. -W. Hu, A. Li, C. -K. Li, J. -D. Li, D. B. Malesani, L. Piro, G. Raman, R. Ricci, E. Troja , et al. (170 additional authors not shown)

Abstract: Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a,… ▽ More Long gamma-ray bursts (GRBs) are believed to originate from core collapse of massive stars. High-redshift GRBs can probe the star formation and reionization history of the early universe, but their detection remains rare. Here we report the detection of a GRB triggered in the 0.5--4 keV band by the Wide-field X-ray Telescope (WXT) on board the Einstein Probe (EP) mission, designated as EP240315a, whose bright peak was also detected by the Swift Burst Alert Telescope and Konus-Wind through off-line analyses. At a redshift of $z=4.859$, EP240315a showed a much longer and more complicated light curve in the soft X-ray band than in gamma-rays. Benefiting from a large field-of-view ($\sim$3600 deg$^2$) and a high sensitivity, EP-WXT captured the earlier engine activation and extended late engine activity through a continuous detection. With a peak X-ray flux at the faint end of previously known high-$z$ GRBs, the detection of EP240315a demonstrates the great potential for EP to study the early universe via GRBs. △ Less

Submitted 25 April, 2024; originally announced April 2024.

Comments: 41 pages, 8 figures, 7 tables

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.

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)

Distribution and Properties of Molecular Gas Toward the Monoceros OB1 Region

Authors: Zi Zhuang, Yang Su, Shiyu Zhang, Xuepeng Chen, Qing-Zeng Yan, Haoran Feng, Li Sun, Xiaoyun Xu, Yan Sun, Xin Zhou, Hongchi Wang, Ji Yang

Abstract: We perform a comprehensive CO study toward the Monoceros OB1 (Mon OB1) region based on the MWISP survey at an angular resolution of about $50''$. The high-sensitivity data, together with the high dynamic range, shows that molecular gas in the $\rm 8^{\circ}\times4^{\circ}$ region displays complicated hierarchical structures and various morphology (e.g., filamentary, cavity-like, shell-like, and ot… ▽ More We perform a comprehensive CO study toward the Monoceros OB1 (Mon OB1) region based on the MWISP survey at an angular resolution of about $50''$. The high-sensitivity data, together with the high dynamic range, shows that molecular gas in the $\rm 8^{\circ}\times4^{\circ}$ region displays complicated hierarchical structures and various morphology (e.g., filamentary, cavity-like, shell-like, and other irregular structures). Based on Gaussian decomposition and clustering for $\mathrm{^{13}CO}$ data, a total of 263 $\mathrm{^{13}CO}$ structures are identified in the whole region, and 88% of raw data flux is recovered. The dense gas with relatively high column density from the integrated CO emission is mainly concentrated in the region where multiple $\rm ^{13}CO$ structures are overlapped. Combining the results of 32 large $\mathrm{^{13}CO}$ structures with distances from Gaia DR3, we estimate an average distance of $\rm 729^{+45}_{-45}~pc$ for the GMC complex. The total mass of the GMC Complex traced by $\mathrm{^{12}CO}$, $\mathrm{^{13}CO}$, and $\mathrm{C^{18}O}$ are $1.1\times10^5~M_\odot$, $4.3\times10^4~M_\odot$, and $8.4\times10^3~M_\odot$, respectively. The dense gas fraction shows a clear difference between Mon OB1 GMC East (12.4%) and Mon OB1 GMC West (3.3%). Our results show that the dense gas environment is closely linked to the nearby star-forming regions. On the other hand, star-forming activities have a great influence on the physical properties of the surrounding molecular gas (e.g., greater velocity dispersion, higher temperatures, and more complex velocity structures, etc.). We also discuss the distribution/kinematics of molecular gas associated with nearby star-forming activities. △ Less

Submitted 14 May, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: 28 pages, 17 figures, match to the version of APJ, 966, 202. The dataset has been released on https://doi.org/10.57760/sciencedb.17451

The Multilayer Nature of Molecular Gas toward the Cygnus Region

Authors: Shiyu Zhang, Yang Su, Xuepeng Chen, Min Fang, Qingzeng Yan, Shaobo Zhang, Yan Sun, Xiaolong Wang, Haoran Feng, Yuehui Ma, Miaomiao Zhang, Zi Zhuang, Xin Zhou, Zhiwei Chen, Ji Yang

Abstract: We study the physical properties and 3D distribution of molecular clouds (MCs) toward the Cygnus region using the MWISP CO survey and Gaia DR3 data. Based on Gaussian decomposition and clustering for $\rm ^{13}CO$ lines, over 70% of the fluxes are recovered. With the identification result of $\rm ^{13}CO$ structures, two models are designed to measure the distances of the molecular gas in velocity… ▽ More We study the physical properties and 3D distribution of molecular clouds (MCs) toward the Cygnus region using the MWISP CO survey and Gaia DR3 data. Based on Gaussian decomposition and clustering for $\rm ^{13}CO$ lines, over 70% of the fluxes are recovered. With the identification result of $\rm ^{13}CO$ structures, two models are designed to measure the distances of the molecular gas in velocity crowding regions. The distances of more than 200 large $\rm ^{13}CO$ structures are obtained toward the 150 square degree region. Additionally, tens of the identified MC structures coincide well with masers and/or intense mid-IR emission. We find multiple gas layers toward the region: (1) the extensive gas structures composing the Cygnus Rift from 700 pc to 1 kpc across the whole region; (2) the $\sim$ 1.3 kpc gas layer mainly in the Cygnus X South region; and (3) the 1.5 kpc dense filament at the Cygnus X North region and many cometary clouds shaped by Cygnus OB2. We also note that the spatial distribution of YSO candidates is generally consistent with the molecular gas structures. The total molecular mass of the Cygnus region is estimated to be $\sim 2.7\times10^{6}~M_{\odot}$ assuming an X-factor ratio $X_{\rm CO} = 2 \times 10^{20} \rm cm^{-2} (K\cdot km\cdot s^{-1})^{-1}$. The foreground Cygnus Rift contributes $\sim$25% of the molecular mass in the whole region. Our work presents a new 3D view of the MCs' distribution toward the Cygnus X region, as well as the exact molecular gas mass distribution in the foreground Cygnus Rift. △ Less

Submitted 23 April, 2024; v1 submitted 29 February, 2024; originally announced March 2024.

Comments: 51 pages, 26 figures, 4 tables, to match the AJ version (2024 AJ 167 220Z). The data can be found at doi: 10.57760/sciencedb.16716

New method for estimating molecular cloud distances based on Gaia, 2MASS, and the TRILEGAL galaxy model

Authors: Juan Mei, Zhiwei Chen, Zhibo Jiang, Sheng Zheng, Haoran Feng

Abstract: We propose a new method for estimating the distances of molecular clouds traced by CO line emission. Stars from 2MASS and Gaia EDR3 are selected as on-cloud stars when they are projected on a cloud. The background on-cloud stars have redder colors on average than the foreground stars. Instead of searching for stars projected away from the cloud, we employed the TRILEGA galaxy model to mimic the st… ▽ More We propose a new method for estimating the distances of molecular clouds traced by CO line emission. Stars from 2MASS and Gaia EDR3 are selected as on-cloud stars when they are projected on a cloud. The background on-cloud stars have redder colors on average than the foreground stars. Instead of searching for stars projected away from the cloud, we employed the TRILEGA galaxy model to mimic the stellar population without cloud extinction along the sightline toward the cloud. Our method does not require an exact boundary of a cloud. The boundaries are highly variable and depend on the sensitivity of the molecular line data. For each cloud, we compared the distributions of on-cloud stars to the TRILEGAL stellar populations in the diagram of $J-K_s$ color versus distance. The intrinsic $J-K_s$ colors of main-sequence and evolved stars from TRILEGAL were considered separately, and they were used as the baseline for subtracting the observed $J-K_s$ colors. The baseline-corrected $J-K_s$ color was deployed with the Bayesian analysis and Markov chain Monte Carlo sampling to determine the distance at which the $J-K_s$ color jump is largest. This method was successfully applied to measure the distances of 27 molecular clouds, which were selected from previously published cloud samples. By replacing TRILEGAL with the GALAXIA galaxy model, we were able to measure the distances for 21 of the 27 clouds. The distances of the 21 clouds based on the GALAXIA model agree well with those based on the TRILEGAL model. The distances of the 27 clouds estimated by this method are consistent with previous estimates. We will apply this new method to a larger region of the gaseous galactic plane, in particular, for the inner galactic region, where a region free of CO emission is hard to separate from the crowded field of clouds. △ Less

Submitted 5 May, 2024; v1 submitted 27 February, 2024; originally announced February 2024.

Comments: 11 pages, 11 figures, accepted by A&A

Journal ref: A&A 685, A39 (2024)

Correcting the Contamination of Second-order Spectra: Improving Hα Measurements in Reverberation Mapping Campaigns

Authors: Wen-Zhe Xi, Kai-Xing Lu, Hai-Cheng Feng, Sha-Sha Li, Jin-Ming Bai, Rui-Lei Zhou, Hong-Tao Liu, Jian-Guo Wang

Abstract: Long-term spectroscopic monitoring campaigns on active galactic nuclei (AGNs) provide a wealth of information about its interior structure and kinematics. However, a number of the observations suffer from the contamination of second-order spectra (SOS) which will introduce some undesirable uncertainties at the red side of the spectra. In this paper, we test the effect of SOS and propose a method t… ▽ More Long-term spectroscopic monitoring campaigns on active galactic nuclei (AGNs) provide a wealth of information about its interior structure and kinematics. However, a number of the observations suffer from the contamination of second-order spectra (SOS) which will introduce some undesirable uncertainties at the red side of the spectra. In this paper, we test the effect of SOS and propose a method to correct it in the time domain spectroscopic data using the simultaneously observed comparison stars. Based on the reverberation mapping (RM) data of NGC 5548 in 2019, one of the most intensively monitored AGNs by the Lijiang 2.4 m telescope, we find that the scientific object, comparison star, and spectrophotometric standard star can jointly introduce up to similar to 30% SOS for Grism 14. This irregular but smooth SOS significantly affects the flux density and profile of the emission line, while having little effect on the light curve. After applying our method to each spectrum, we find that the SOS can be corrected effectively. The deviation between corrected and intrinsic spectra is similar to 2%, and the impact of SOS on time lag is very minor. This method makes it possible to obtain the H alpha RM measurements from archival data provided that the spectral shape of the AGN under investigation does not have a large change. △ Less

Submitted 22 January, 2024; originally announced January 2024.

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

The First Polarimetric View on Quasi-Periodic Oscillations in a Black Hole X-ray Binary

Authors: Qingchang Zhao, Lian Tao, Hancheng Li, Shuangnan Zhang, Hua Feng, Mingyu Ge, Long Ji, Yanan Wang, Yue Huang, Xiang Ma, Liang Zhang, Jinlu Qu, Yanjun Xu, Shu Zhang, Qianqing Yin, Qingcang Shui, Ruican Ma, Shujie Zhao, Panping Li, Zixu Yang, Hexin Liu, WeiYu

Abstract: We present the first polarimetric analysis of Quasi-Periodic Oscillations (QPO) in a black hole binary utilizing \textit{IXPE} data. Our study focuses on Swift J1727.8--1613, which experienced a massive outburst that was observed by various telescopes across different wavelengths. The \textit{IXPE} observation we studied was conducted during the Hard-Intermediate state. The polarization degree (PD… ▽ More We present the first polarimetric analysis of Quasi-Periodic Oscillations (QPO) in a black hole binary utilizing \textit{IXPE} data. Our study focuses on Swift J1727.8--1613, which experienced a massive outburst that was observed by various telescopes across different wavelengths. The \textit{IXPE} observation we studied was conducted during the Hard-Intermediate state. The polarization degree (PD) and polarization angle (PA) were measured at 4.28$\pm$0.20\% and $1.9^{\circ}\pm1.4^{\circ}$, respectively. Remarkably, significant QPO signals were detected during this observation, with a QPO frequency of approximately 1.34 Hz and a fractional root-mean-square (RMS) amplitude of about 12.3\%. Furthermore, we conducted a phase-resolved analysis of the QPO using the Hilbert-Huang transform technique. The photon index showed a strong modulation with respect to the QPO phase. In contrast, the PD and PA exhibit no modulations in relation to the QPO phase, which is inconsistent with the expectation of the Lense-Thirring precession of the inner flow. Further theoretical studies are needed to conform with the observational results. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted for publication in APJL

Measuring the Virial Factor in SDSS DR7 AGNs with Redshifted H$β$ and H$α$ Broad Emission Lines

Authors: H. T. Liu, Hai-Cheng Feng, Sha-Sha Li, J. M. Bai, H. Z. Li

Abstract: Under the hypothesis of gravitational redshift induced by the central supermassive black hole, and based on line widths and shifts of redward shifted H$β$ and H$α$ broad emission lines for more than 8000 SDSS DR7 AGNs, we measure the virial factor in determining supermassive black hole masses. The virial factor had been believed to be independent of accretion radiation pressure on gas clouds in br… ▽ More Under the hypothesis of gravitational redshift induced by the central supermassive black hole, and based on line widths and shifts of redward shifted H$β$ and H$α$ broad emission lines for more than 8000 SDSS DR7 AGNs, we measure the virial factor in determining supermassive black hole masses. The virial factor had been believed to be independent of accretion radiation pressure on gas clouds in broad-line region (BLR), and only dependent on inclination effects of BLR. The virial factor measured spans a very large range. For the vast majority of AGNs ($>$96%) in our samples, the virial factor is larger than $f=1$ usually used in literatures. The $f$ correction makes the percent of high-accreting AGNs decrease by about 100 times. There are positive correlations of $f$ with the dimensionless accretion rate and Eddington ratio. The redward shifts of H$β$ and H$α$ are mainly the gravitational origin, confirmed by a negative correlation between the redward shift and the dimensionless radius of BLR. Our results show that radiation pressure force is a significant contributor to the measured virial factor, containing the inclination effects of BLR. The usually used values of $f$ should be corrected for high-accreting AGNs, especially high redshift quasars. The $f$ correction increases their masses by one--two orders of magnitude, which will make it more challenging to explain the formation and growth of supermassive black holes at high redshifts. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 34 pages with 11 figures and 9 tables, Accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:2202.08949

Journal ref: 2024, ApJ, 963, 30

MeV Astrophysical Spectroscopic Surveyor (MASS): A Compton Telescope Mission Concept

Authors: Jiahuan Zhu, Xutao Zheng, Hua Feng, Ming Zeng, Chien-You Huang, Jr-Yue Hsiang, Hsiang-Kuang Chang, Hong Li, Hao Chang, Xiaofan Pan, Ge Ma, Qiong Wu, Yulan Li, Xuening Bai, Mingyu Ge, Long Ji, Jian Li, Yangping Shen, Wei Wang, Xilu Wang, Binbin Zhang, Jin Zhang

Abstract: We propose a future mission concept, the MeV Astrophysical Spectroscopic Surveyor (MASS), which is a large area Compton telescope using 3D position sensitive cadmium zinc telluride (CZT) detectors optimized for emission line detection. The payload consists of two layers of CZT detectors in a misaligned chessboard layout, with a total geometric area of 4096 cm$^2$ for on-axis observations. The dete… ▽ More We propose a future mission concept, the MeV Astrophysical Spectroscopic Surveyor (MASS), which is a large area Compton telescope using 3D position sensitive cadmium zinc telluride (CZT) detectors optimized for emission line detection. The payload consists of two layers of CZT detectors in a misaligned chessboard layout, with a total geometric area of 4096 cm$^2$ for on-axis observations. The detectors can be operated at room-temperature with an energy resolution of 0.6\% at 0.662 MeV. The in-orbit background is estimated with a mass model. At energies around 1 MeV, a line sensitivity of about $10^{-5}$ photons cm$^{-2}$ s$^{-1}$ can be obtained with a 1 Ms observation. The main science objectives of MASS include nucleosynthesis in astrophysics and high energy astrophysics related to compact objects and transient sources. The payload CZT detectors weigh roughly 40 kg, suggesting that it can be integrated into a micro- or mini-satellite. We have constructed a pathfinder, named as MASS-Cube, to have a direct test of the technique with 4 detector units in space in the near future. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: accepted for publication in Experimental Astronomy

The bright black hole X-ray binary 4U 1543--47 during 2021 outburst: a thick accretion disk inflated by high luminosity

Authors: S. J. Zhao, L. Tao, P. P. Li, R. Soria, H. Feng, Y. X. Zhang, R. C. Ma, W. D. Zhang, E. L. Qiao, Q. Q. Yin, S. N. Zhang, L. Zhang, Q. C. Bu, X. Ma, Y. Huang, M. Y. Ge, X. B. Li, Q. C. Zhao, J. Q. Peng, Y. X. Xiao

Abstract: The black hole X-ray binary source 4U 1543--47 experienced a super-Eddington outburst in 2021, reaching a peak flux of up to $\sim1.96\times10^{-7}\rm erg\ \rm cm^{-2}\ \rm s^{-1}$ ($\sim 8.2$ Crab) in the 2--10\,keV band. Soon after the outburst began, it rapidly transitioned into the soft state. Our goal is to understand how the accretion disk structure deviates from a standard thin disk when th… ▽ More The black hole X-ray binary source 4U 1543--47 experienced a super-Eddington outburst in 2021, reaching a peak flux of up to $\sim1.96\times10^{-7}\rm erg\ \rm cm^{-2}\ \rm s^{-1}$ ($\sim 8.2$ Crab) in the 2--10\,keV band. Soon after the outburst began, it rapidly transitioned into the soft state. Our goal is to understand how the accretion disk structure deviates from a standard thin disk when the accretion rate is near Eddington. To do so, we analyzed spectra obtained from quasi-simultaneous observations conducted by the Hard X-ray Modulation Telescope (Insight-HXMT), the Nuclear Spectroscopic Telescope Array (NuSTAR), and the Neil Gehrels Swift Observatory (Swift). These spectra are well-fitted by a model comprising a disk, a weak corona, and a reflection component. We suggest that the reflection component is caused by disk self-irradiation, that is by photons emitted from the inner disk which return to the accretion disk surface, as their trajectories are bent by the strong gravity field. In this scenario, the best-fitting parameters imply that the reflected flux represents more than half of the total flux. Using general relativistic ray-tracing simulations, we show that this scenario is viable when the disk becomes geometrically thick, with a funnel-like shape, as the accretion rate is near or above the Eddington limit. In the specific case of 4U 1543--47, an angle $\gtrsim$ 45 deg between the disk surface and the equatorial plane can explain the required amount of self-irradiation. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy and Astrophysics. 15 pages, 4 tables, 12 figures

Journal ref: A&A 685, A42 (2024)

Noise discrimination method based on charge distribution of CMOS detectors for soft X-ray

Authors: Xinchao Fang, Jirong Cang, Qiong Wu, Hua Feng, Ming Zeng

Abstract: Complementary metal-oxide semiconductor (CMOS) sensors have been widely used as soft X-ray detectors in several fields owing to their recent developments and unique advantages. The parameters of CMOS detectors have been extensively studied and evaluated. However, the key parameter signal-to-noise ratio in certain fields has not been sufficiently studied. In this study, we analysed the charge distr… ▽ More Complementary metal-oxide semiconductor (CMOS) sensors have been widely used as soft X-ray detectors in several fields owing to their recent developments and unique advantages. The parameters of CMOS detectors have been extensively studied and evaluated. However, the key parameter signal-to-noise ratio in certain fields has not been sufficiently studied. In this study, we analysed the charge distribution of the CMOS detector GSENSE2020BSI and proposed a two-dimensional segmentation method to discriminate signals according to the charge distribution. The effect of the two-dimensional segmentation method on the GSENSE2020BSI dectector was qualitatively evaluated. The optimal feature parameters used in the two-dimensional segmentation method was studied for G2020BSI. However, the two-dimensional segmentation method is insensitive to feature parameters. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 19 pages, 13 figures, submitted to NIM-A

Research on the X-Ray Polarization Deconstruction Method Based on Hexagonal Convolutional Neural Network

Authors: Ya-Nan Li, Jia-Huan Zhu, Huai-Zhong Gao, Hong Li, Ji-Rong Cang, Zhi Zeng, Hua Feng, Ming Zeng

Abstract: Track reconstruction algorithms are critical for polarization measurements. In addition to traditional moment-based track reconstruction approaches, convolutional neural networks (CNN) are a promising alternative. However, hexagonal grid track images in gas pixel detectors (GPD) for better anisotropy do not match the classical rectangle-based CNN, and converting the track images from hexagonal to… ▽ More Track reconstruction algorithms are critical for polarization measurements. In addition to traditional moment-based track reconstruction approaches, convolutional neural networks (CNN) are a promising alternative. However, hexagonal grid track images in gas pixel detectors (GPD) for better anisotropy do not match the classical rectangle-based CNN, and converting the track images from hexagonal to square results in loss of information. We developed a new hexagonal CNN algorithm for track reconstruction and polarization estimation in X-ray polarimeters, which was used to extract emission angles and absorption points from photoelectron track images and predict the uncertainty of the predicted emission angles. The simulated data of PolarLight test were used to train and test the hexagonal CNN models. For individual energies, the hexagonal CNN algorithm produced 15-30% improvements in modulation factor compared to moment analysis method for 100% polarized data, and its performance was comparable to rectangle-based CNN algorithm newly developed by IXPE team, but at a much less computational cost. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 21 pages, 12 figures, submitted to NST

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

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

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

Submitted 25 October, 2023; originally announced October 2023.

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

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

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

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

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

Comments: 49pages, 11figures

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

Systematic Search for Water Fountain Candidates using the Databases of Circumstellar Maser Sources

Authors: Haichen Fan, Jun-ichi Nakashima, D. Engels, Yong Zhang, Jian-Jie Qiu, Huan-Xue Feng, Jia-Yong Xie, Hiroshi Imai, Chih-Hao Hsia

Abstract: Water fountains (WFs) are thought to be objects in the morphological evolution of the circumstellar envelopes of low- and intermediate-mass evolved stars, transitioning from spherically symmetric to asymmetric shapes. We used databases of circumstellar 1612 MHz OH and 22.235 GHz H$_2$O maser sources to search for new WF candidates using the criterion of a larger velocity range of the H$_2$O maser… ▽ More Water fountains (WFs) are thought to be objects in the morphological evolution of the circumstellar envelopes of low- and intermediate-mass evolved stars, transitioning from spherically symmetric to asymmetric shapes. We used databases of circumstellar 1612 MHz OH and 22.235 GHz H$_2$O maser sources to search for new WF candidates using the criterion of a larger velocity range of the H$_2$O maser emission compared to that of the OH maser emission. Thus, it is in principle possible to identify WFs with an H$_2$O velocity ranges smaller than those for the previously known WFs. For the OH maser line, we analyzed database entries of 8,474 observations from 2,195 sources, and 6,085 observations from 3,642 sources for H$_2$O maser line. After a close examination of the velocity ranges and line profiles, we identified 11 sources that meet the criterion mentioned above. We examined the IRAS colors of the selected sources and found that two of them (IRAS 19069+0916 and IRAS 19319+2214) are in the color region for post-AGB stars. We find that the maser velocity criterion can discover other astrophysically interesting objects than just WFs. Such objects may include peculiar planetary nebulae with maser emissions and stellar merger remnants. △ Less

Submitted 8 October, 2023; originally announced October 2023.

Comments: 34 pages, 10 figures, Accepted for publication in ApJS

In-orbit background and sky survey simulation study of POLAR-2/LPD

Authors: Zu-Ke Feng, Hong-Bang Liu, Fei Xie, Huan-Bo Feng, Qian-Nan Mai, Jiang-Chuan Tuo, Qian Zhong, Jian-Chao Sun, Jiang He, Yuan-Hao Wang, Qian Liu, Di-Fan Yi, Rui-Ting Ma, Bin-Long Wang, Zhen-Yu Tang, Shuang-Nan Zhang, En-Wei Liang

Abstract: The Low-Energy X-ray Polarization Detector (LPD) is one of the payloads in the POLAR-2 experiment, designed as an external payload for the China Space Station (CSS) deployment in early 2024. LPD is specifically designed to observe the polarization of Gamma-Ray Bursts (GRBs) prompt emission in the energy range of 2-10 keV, with a wide field of view (FoV) of 90 degrees in preliminary design. This ob… ▽ More The Low-Energy X-ray Polarization Detector (LPD) is one of the payloads in the POLAR-2 experiment, designed as an external payload for the China Space Station (CSS) deployment in early 2024. LPD is specifically designed to observe the polarization of Gamma-Ray Bursts (GRBs) prompt emission in the energy range of 2-10 keV, with a wide field of view (FoV) of 90 degrees in preliminary design. This observation is achieved using an array of X-ray photoelectric polarimeters based on gas pixel detectors. Due to the wide FoV configuration, the in-orbit background count rate in the soft X-ray range is high, while GRBs themselves also exhibit a high flux in this energy band. In order to assess the contribution of various background components to the total count rate, we conducted detailed simulations using the GEANT4 C++ package. Our simulations encompassed the main interactions within the instrument materials and provided insights into various background components within the wide FoV scheme. The simulation results reveal that among the background components, the primary contributors are the cosmic X-ray background (CXB) and bright X-ray sources. The total background count rate of LPD, after applying the charged particle background rejection algorithm, is approximately 0.55 counts/cm^2/s on average, and it varies with the detector's orbit and pointing direction. Furthermore, we performed comprehensive simulations and comparative analyses of the CXB and X-ray bright sources under different FoVs and detector pointings. These analyses provide valuable insights into the background characteristic for soft X-ray polarimeter with wide FoV. △ Less

Submitted 1 October, 2023; originally announced October 2023.

The molecular clouds in a section of the third Galactic quadrant: observational properties and chemical abundance ratio between CO and its isotopologues

Authors: Chen Wang, Haoran Feng, Ji Yang, Xuepeng Chen, Yang Su, Qing-Zeng Yan, Fujun Du, Yuehui Ma, Jiajun Cai

Abstract: We compare the observational properties between $^{12}$CO, $^{13}$CO, and C$^{18}$O and summarize the observational parameters based on 7069 clouds sample from the Milky Way Imaging Scroll Painting (MWISP) CO survey in a section of the third Galactic quadrant. We find that the $^{13}$CO angular area ($A_{\rm ^{13}CO}$) generally increases with that of $^{12}$CO ($A_{\rm ^{12}CO}$), and the ratio o… ▽ More We compare the observational properties between $^{12}$CO, $^{13}$CO, and C$^{18}$O and summarize the observational parameters based on 7069 clouds sample from the Milky Way Imaging Scroll Painting (MWISP) CO survey in a section of the third Galactic quadrant. We find that the $^{13}$CO angular area ($A_{\rm ^{13}CO}$) generally increases with that of $^{12}$CO ($A_{\rm ^{12}CO}$), and the ratio of $A_{\rm ^{13}CO}$ to $A_{\rm ^{12}CO}$ is 0.38 by linear fitting. We find that the $^{12}$CO and $^{13}$CO flux are tightly correlated as $F_{\rm ^{13}CO}~=~0.17~ F_{\rm ^{12}CO}$ with both fluxes calculated within the $^{13}$CO-bright region. This indicates that the abundance $X_{\rm ^{13}CO}$ is a constant to be 6.5$^{+0.1}_{-0.5}$ $\times 10^{-7}$ for all samples under assumption of local thermodynamic equilibrium (LTE). Additionally, we observed that the X-factor is approximately constant in large sample molecular clouds. Similarly, we find $F_{\rm C^{18}O}~=~0.11~F_{\rm ^{13}CO}$ with both fluxes calculated within C$^{18}$O-bright region, which indicates that the abundance ratios ${X_{\rm ^{13}CO}/X_{\rm C^{18}O}}$ stays the same value 9.7$^{+0.6}_{-0.8}$ across the molecular clouds under LTE assumption. The linear relationships of $F_{\rm ^{12}CO}$ vs. $F_{\rm ^{13}CO}$ and $F_{\rm ^{13}CO}$ vs. $F_{\rm C^{18}O}$ hold not only for the $^{13}$CO-bright region or C$^{18}$O-bright region, but also for the entire molecular cloud scale with lower flux ratio. The abundance ratio ${X_{\rm ^{13}CO}/X_{\rm C^{18}O}}$ inside clouds shows a strong correlation with column density and temperature. This indicates that the ${X_{\rm ^{13}CO}/X_{\rm C^{18}O}}$ is dominated by a combination of chemical fractionation, selectively dissociation, and self-shielding effect inside clouds. △ Less

Submitted 21 August, 2023; originally announced August 2023.

Comments: 11 pages, 16 figures, 1 table, accepted by AJ

Black Hole Accretion with Saturated Magnetic Pressure and Disk Wind

Authors: Jiahui Huang, Hua Feng, Wei-Min Gu, Wen-Biao Wu

Abstract: We construct an analytical black hole accretion disk model that incorporates both magnetic pressure and disk wind, which are found to be important from numerical simulations. A saturated magnetic pressure that relates the Alfven velocity with local Keplerian velocity and gas sound speed is assumed in addition to radiation and gas pressures. The mass accretion rate is assumed to have a power-law fo… ▽ More We construct an analytical black hole accretion disk model that incorporates both magnetic pressure and disk wind, which are found to be important from numerical simulations. A saturated magnetic pressure that relates the Alfven velocity with local Keplerian velocity and gas sound speed is assumed in addition to radiation and gas pressures. The mass accretion rate is assumed to have a power-law form in response to mass loss in the wind. We find three sets of self-consistent solutions that are thermally stable and satisfy the model assumptions. At high accretion rates, the disk is geometrically and optically thick, resembling the slim disk solution. At relatively low accretion rates, our model predicts an accretion flow consisting of a geometrically thin and optically thick outer disk (similar to the standard disk), and a geometrically thick and optically thin inner disk (similar to the advection-dominated accretion flow or ADAF). Thus, this is a natural solution for a truncated disk connected with an inner ADAF, which has been proposed to explain some observations. The magnetic pressure plays a more important role than the outflow in shaping the disk structure. The observed disk luminosity tends to saturate around 8 times the Eddington limit, suggesting that supercritical accretion onto black holes can be used for black hole mass estimate, or a standard candle with known black hole masses. △ Less

Submitted 13 July, 2023; originally announced July 2023.

Comments: 10 pages, 9 figures, accepted for publication in ApJ

Toward a stellar population catalog in the Kilo Degree Survey: the impact of stellar recipes on stellar masses and star formation rates

Authors: Linghua Xie, Nicola R. Napolitano, Xiaotong Guo, Crescenzo Tortora, Haicheng Feng, Antonios Katsianis, Rui Li, Sirui Wu, Mario Radovich, Leslie K. Hunt, Yang Wang, Lin Tang, Baitian Tang, Zhiqi Huang

Abstract: The Kilo Degree Survey (KiDS) is currently the only sky survey providing optical ($ugri$) plus near-infrared (NIR, $ZYHJK_S$) seeing matched photometry over an area larger than 1000 $\rm deg^2$. This is obtained by incorporating the NIR data from the VISTA Kilo Degree Infrared Galaxy (VIKING) survey, covering the same KiDS footprint. As such, the KiDS multi-wavelength photometry represents a uniqu… ▽ More The Kilo Degree Survey (KiDS) is currently the only sky survey providing optical ($ugri$) plus near-infrared (NIR, $ZYHJK_S$) seeing matched photometry over an area larger than 1000 $\rm deg^2$. This is obtained by incorporating the NIR data from the VISTA Kilo Degree Infrared Galaxy (VIKING) survey, covering the same KiDS footprint. As such, the KiDS multi-wavelength photometry represents a unique dataset to test the ability of stellar population models to return robust photometric stellar mass ($M_*$) and star-formation rate (SFR) estimates. Here we use a spectroscopic sample of galaxies for which we possess $u g r i Z Y J H K_s$ ``gaussianized'' magnitudes from KiDS data release 4. We fit the spectral energy distribution from the 9-band photometry using: 1) three different popular libraries of stellar {population} templates, 2) single burst, simple and delayed exponential star-formation history models, and 3) a wide range of priors on age and metallicity. As template fitting codes we use two popular softwares: LePhare and CIGALE. We investigate the variance of the stellar masses and the star-formation rates from the different combinations of templates, star formation recipes and codes to assess the stability of these estimates and define some ``robust'' median quantities to be included in the upcoming KiDS data releases. As a science validation test, we derive the mass function, the star formation rate function, and the SFR-$M_*$ relation for a low-redshift ($z<0.5$) sample of galaxies, that result in excellent agreement with previous literature data. The final catalog, containing $\sim290\,000$ galaxies with redshift $0.01<z<0.9$, is made publicly available. △ Less

Submitted 11 December, 2023; v1 submitted 9 July, 2023; originally announced July 2023.

Comments: 39 pages,4 tables, 14 figures, SCIENCE CHINA Physics, Mechanics & Astronomy Accepted

Deep HI Mapping of Stephan's Quintet and Its Neighborhood

Authors: Cheng Cheng, Cong Kevin Xu, P. N. Appleton, P. -A. Duc, N. -Y. Tang, Y. S. Dai, J. -S. Huang, U. Lisenfeld, F. Renaud, Chuan He, Hai-Cheng Feng

Abstract: We carried out deep mapping observations of the atomic hydrogen (HI) 21 cm line emission in a field centered on the famous galaxy group Stephan's Quintet (SQ), using the Five-hundred-meter Aperture Spherical Telescope (FAST) equipped with the 19-Beam Receiver. The final data cube reaches an HI column density sensitivity of $5 σ= 2.1\times 10^{17}$ cm$^{-2}$ per 20 km s$^{-1}$ channel with an angul… ▽ More We carried out deep mapping observations of the atomic hydrogen (HI) 21 cm line emission in a field centered on the famous galaxy group Stephan's Quintet (SQ), using the Five-hundred-meter Aperture Spherical Telescope (FAST) equipped with the 19-Beam Receiver. The final data cube reaches an HI column density sensitivity of $5 σ= 2.1\times 10^{17}$ cm$^{-2}$ per 20 km s$^{-1}$ channel with an angular resolution of $4'.0$. The discovery of a large diffuse feature of the HI emission in the outskirt of the intragroup medium of SQ was reported in a previous paper (Xu et al. 2022). Here we present a new study of the total HI emission of SQ and the detection of several neighboring galaxies, exploiting the high sensitivity and the large sky coverage of the FAST observations. A total HI mass of $M_{\rm HI} = 3.48 \pm 0.35 \times 10^{10}\; M_\odot$ is found for SQ, which is significantly higher than previous measurements in the literature. This indicates that, contrary to earlier claims, SQ is not HI deficient. The excessive HI gas is mainly found in the velocity ranges of 6200 - 6400 km s$^{-1}$ and 6800 - 7000 km s$^{-1}$, which was undetected in previous observations that are less sensitive than ours. Our results suggest that the ``missing HI" in compact groups may be hidden in the low-density diffuse neutral gas instead of in the ionized gas. △ Less

Submitted 19 June, 2023; originally announced June 2023.

Comments: 20 pages, 5 figures, Accepted by ApJ

Identification of Bubble Nebulae around NGC 55 ULX-1 with MUSE Observations

Authors: Changxing Zhou, Hua Feng, Fuyan Bian

Abstract: Using the Multi Unit Spectroscopic Explorer (MUSE) instrument on the Very Large Telescope, we identified three bubble nebulae (denoted as A, B, and C) around an ultraluminous X-ray source (ULX) in NGC 55. Bubble A shows a regular elliptical shape surrounding the ULX, with a morphology similar to the canonical ULX bubble around NGC 1313 X-2. It is most likely inflated by the ULX disk wind with a me… ▽ More Using the Multi Unit Spectroscopic Explorer (MUSE) instrument on the Very Large Telescope, we identified three bubble nebulae (denoted as A, B, and C) around an ultraluminous X-ray source (ULX) in NGC 55. Bubble A shows a regular elliptical shape surrounding the ULX, with a morphology similar to the canonical ULX bubble around NGC 1313 X-2. It is most likely inflated by the ULX disk wind with a mechanical power close to $10^{39}$ erg s$^{-1}$. Bubble B lies 11\arcsec\ away from the ULX on the sky plane and is not contiguous to Bubble A. It displays a bow shock like morphology, and is likely driven by a collimated dark jet from the ULX with a mechanical power of about $3 \times 10^{38}$ erg s$^{-1}$. If this scenario is correct, we predict that Bubble B should present radio emission with a flux of about $1 - 10^2$ $μ$Jy at 5 GHz. Bubble C appears within Bubble A, with a velocity and velocity dispersion distinct from the rest of Bubble A. Its nature is unclear and could be part of Bubble A as a result of low local density. The optical counterpart of ULX-1 exhibits a broad H$α$, consistent with emission from a hot disk wind. △ Less

Submitted 23 August, 2023; v1 submitted 11 June, 2023; originally announced June 2023.

Comments: ApJ to appear

The First LHAASO Catalog of Gamma-Ray Sources

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

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

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

Comments: 40 pages, 13 figures, 4 tables

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

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

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

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

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

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

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

X-ray Polarimetry of the accreting pulsar 1A~0535+262 in the supercritical state with PolarLight

Authors: Xiangyun Long, Hua Feng, Hong Li, Ling-Da Kong, Jeremy Heyl, Long Ji, Lian Tao, Fabio Muleri, Qiong Wu, Jiahuan Zhu, Jiahui Huang, Massimo Minuti, Weichun Jiang, Saverio Citraro, Hikmat Nasimi, Jiandong Yu, Ge Jin, Ming Zeng, Peng An, Luca Baldini, Ronaldo Bellazzini, Alessandro Brez, Luca Latronico, Carmelo Sgrò, Gloria Spandre , et al. (3 additional authors not shown)

Abstract: The X-ray pulsar 1A 0535+262 exhibited a giant outburst in 2020, offering us a unique opportunity for X-ray polarimetry of an accreting pulsar in the supercritical state. Measurement with PolarLight yielded a non-detection in 3-8 keV; the 99% upper limit of the polarization fraction (PF) is found to be 0.34 averaged over spin phases, or 0.51 based on the rotating vector model. No useful constraint… ▽ More The X-ray pulsar 1A 0535+262 exhibited a giant outburst in 2020, offering us a unique opportunity for X-ray polarimetry of an accreting pulsar in the supercritical state. Measurement with PolarLight yielded a non-detection in 3-8 keV; the 99% upper limit of the polarization fraction (PF) is found to be 0.34 averaged over spin phases, or 0.51 based on the rotating vector model. No useful constraint can be placed with phase resolved polarimetry. These upper limits are lower than a previous theoretical prediction of 0.6-0.8, but consistent with those found in other accreting pulsars, like Her X-1, Cen X-3, 4U 1626-67, and GRO J1008-57, which were in the subcritical state, or at least not confidently in the supercritical state, during the polarization measurements. Our results suggest that the relatively low PF seen in accreting pulsars cannot be attributed to the source not being in the supercritical state, but could be a general feature. △ Less

Submitted 26 April, 2023; originally announced April 2023.

Comments: accepted for publication in ApJ

Identification of a Helium Donor Star in NGC 247 ULX-1

Authors: Changxing Zhou, Hua Feng, Fuyan Bian

Abstract: With Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE) observations, we detected highly variable helium emission lines from the optical counterpart of the supersoft ultraluminous X-ray source (ULX) NGC 247 ULX-1. No Balmer lines can be seen in the source spectrum. This is the first evidence for the presence of a helium donor star in ULXs, consistent with a prediction that helium… ▽ More With Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE) observations, we detected highly variable helium emission lines from the optical counterpart of the supersoft ultraluminous X-ray source (ULX) NGC 247 ULX-1. No Balmer lines can be seen in the source spectrum. This is the first evidence for the presence of a helium donor star in ULXs, consistent with a prediction that helium donor stars may be popular in ULXs. The helium lines with a FWHM of about 200 km/s are likely produced on the outer accretion disk. Their strong variation implies that the central X-ray source can be significantly obscured to the outer disk. Also, a ring or a double-ring structure is revealed in the MUSE image. It is unknown whether or not it is related to the progenitor of the ULX binary. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: accepted for publication in ApJ

In search for infalling clumps in molecular clouds -- A catalogue of CO blue-profiles

Authors: Zhibo Jiang, Shaobo Zhang, Zhiwei Chen, Yang Yang, Shuling Yu, Haoran Feng, Ji Yang, the MWISP group

Abstract: We have started a systematic survey of molecular clumps with infall motions to study the very early phase of star formation. Our first step is to utilize the data products by MWISP to make an unbiased survey for blue asymmetric line profiles of CO isotopical molecules. Within a total area of $\sim$ 2400 square degrees nearby the Galactic plane, we have found 3533 candidates showing blue-profiles,… ▽ More We have started a systematic survey of molecular clumps with infall motions to study the very early phase of star formation. Our first step is to utilize the data products by MWISP to make an unbiased survey for blue asymmetric line profiles of CO isotopical molecules. Within a total area of $\sim$ 2400 square degrees nearby the Galactic plane, we have found 3533 candidates showing blue-profiles, in which 3329 are selected from the $^{12}$CO&$^{13}$CO pair and 204 are from the $^{13}$CO&C$^{18}$O pair. Exploration of the parametric spaces suggests our samples are in the cold phase with relatively high column densities ready for star formation. Analysis of the spatial distribution of our samples suggests that they exist virtually in all major components of the Galaxy. The vertical distribution suggest that the sources are located mainly in the thick disk of $\sim$ 85 parsec, but still a small part are located far beyond Galactic midplane. Our follow-up observation indicates that these candidates are a good sample to start a search for infall motions, and to study the condition of very early phase of star formation. △ Less

Submitted 20 March, 2023; originally announced March 2023.

GRB 220408B: A Three-Episode Burst from a Precessing Jet

Authors: Zijian Zhang, Yihan Yin, Chenyu Wang, Xiangyu Ivy Wang, Jun Yang, Yan-Zhi Meng, Zi-Ke Liu, Guo-Yin Chen, Xiaoping Fu, Huaizhong Gao, Sihao Li, Yihui Liu, Xiangyun Long, Yong-Chang Ma, Xiaofan Pan, Yuanze Sun, Wei Wu, Zirui Yang, Zhizhen Ye, Xiaoyu Yu, Shuheng Zhao, Xutao Zheng, Tao Zhou, Qing-Wen Tang, Qiurong Yan , et al. (5 additional authors not shown)

Abstract: Jet precession has previously been proposed to explain the apparently repeating features in the light curves of a few gamma-ray bursts (GRBs). In this {\it Letter}, we further apply the precession model to a bright GRB 220408B by examining both its temporal and spectral consistency with the predictions of the model. As one of the recently confirmed GRBs observed by our GRID CubeSat mission, GRB 22… ▽ More Jet precession has previously been proposed to explain the apparently repeating features in the light curves of a few gamma-ray bursts (GRBs). In this {\it Letter}, we further apply the precession model to a bright GRB 220408B by examining both its temporal and spectral consistency with the predictions of the model. As one of the recently confirmed GRBs observed by our GRID CubeSat mission, GRB 220408B is noteworthy as it exhibits three apparently similar emission episodes. Furthermore, the similarities are reinforced by their strong temporal correlations and similar features in terms of spectral evolution and spectral lags. Our analysis demonstrates that these features can be well explained by the modulated emission of a Fast-Rise-Exponential-Decay (FRED) shape light curve intrinsically produced by a precessing jet with a precession period of $18.4 \pm 0.2$ seconds, a nutation period of $11.1 \pm 0.2$ seconds and viewed off-axis. This study provides a straightforward explanation for the complex yet similar multi-episode GRB light curves. △ Less

Submitted 6 February, 2023; originally announced February 2023.

Comments: 11 pages, 3 tables, 8 figures

Global Three-Dimensional Radiation Magnetohydrodynamic Simulations of Accretion onto a Stellar Mass Black Hole at Sub- and Near-critical Accretion Rates

Authors: Jiahui Huang, Yan-Fei Jiang, Hua Feng, Shane W. Davis, James M. Stone, Matthew J. Middleton

Abstract: We present global 3D radiation magnetohydrodynamical simulations of accretion onto a 6.62 solar mass black hole with quasi-steady state accretion rates reaching 0.016 to 0.9 times the critical accretion rate, which is defined as the accretion rate to power the Eddington luminosity assuming a 10% radiative efficiency, in different runs. The simulations show no sign of thermal instability over hundr… ▽ More We present global 3D radiation magnetohydrodynamical simulations of accretion onto a 6.62 solar mass black hole with quasi-steady state accretion rates reaching 0.016 to 0.9 times the critical accretion rate, which is defined as the accretion rate to power the Eddington luminosity assuming a 10% radiative efficiency, in different runs. The simulations show no sign of thermal instability over hundreds of thermal timescales at 10 $r_{\rm g}$. The energy dissipation happens close to the mid-plane in the near-critical runs and near the disk surface in the low accretion rate run. The total radiative luminosity inside $\sim$20 $r_{\rm g}$ is about 1% to 30% the Eddington limit, with a radiative efficiency of about 6% and 3%, respectively, in the sub- and near-critical accretion regimes. In both cases, self-consistent turbulence generated by the magnetorotational instability (MRI) leads to angular momentum transfer, and the disk is supported by magnetic pressure. Outflows from the central low-density funnel with a terminal velocity of $\sim$0.1$c$ are seen only in the near-critical runs. We conclude that these magnetic pressure dominated disks are thermally stable and thicker than the $α$ disk, and the effective temperature profiles are much flatter than that in the $α$ disks. The magnetic pressure of these disks are comparable within an order of magnitude with the previous analytical magnetic pressure dominated disk model. △ Less

Submitted 30 January, 2023; originally announced January 2023.

Comments: 17 pages, 13 figures, 3 tables, accepted for publication in ApJ

Choked accretion onto Kerr-Sen black holes in Einstein-Maxwell-dilaton-axion gravity

Authors: Haiyuan Feng, Yingdong Wu, Rong-Jia Yang, Leonardo Modesto

Abstract: We investigate the process of an ultrarelativistic fluid accreted onto axisymmetric Kerr-Sen black holes in Einstein-Maxwell-dilaton-axion theory. We obtain the solution describing the velocity potential of a stationary irrotational fluid with a stiff equation of state and the solution for the streamlined diagram of the quadrupolar flow. We also investigate how the solution's coefficients and the… ▽ More We investigate the process of an ultrarelativistic fluid accreted onto axisymmetric Kerr-Sen black holes in Einstein-Maxwell-dilaton-axion theory. We obtain the solution describing the velocity potential of a stationary irrotational fluid with a stiff equation of state and the solution for the streamlined diagram of the quadrupolar flow. We also investigate how the solution's coefficients and the stagnation points are affected by the parameters. The injection rate, the ejection rate, and the critical angle are discussed in detail. We find that with an increasing dilaton parameter the ratio of the ejection rate to the injection rate increases and that the radiative efficiency is larger, while the redshift is lower, compared to the Kerr black hole. △ Less

Submitted 14 March, 2024; v1 submitted 7 January, 2023; originally announced January 2023.

Comments: 25pages,10 figures

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

Bubble in the Whale: Identifying the Optical Counterparts and Extended Nebula for the Ultraluminous X-ray Sources in NGC 4631

Authors: Jing Guo, Jianfeng Wu, Hua Feng, Zheng Cai, Ping Zhou, Changxing Zhou, Shiwu Zhang, Junfeng Wang, Mouyuan Sun, Wei-Min Gu, Shan-Shan Weng, Jifeng Liu

Abstract: We present a deep optical imaging campaign on the starburst galaxy NGC 4631 with CFHT/MegaCam. By supplementing the HST/ACS and Chandra/ACIS archival data, we search for the optical counterpart candidates of the five brightest X-ray sources in this galaxy, four of which are identified as ultraluminous X-ray sources (ULXs). The stellar environments of the X-ray sources are analyzed using the extinc… ▽ More We present a deep optical imaging campaign on the starburst galaxy NGC 4631 with CFHT/MegaCam. By supplementing the HST/ACS and Chandra/ACIS archival data, we search for the optical counterpart candidates of the five brightest X-ray sources in this galaxy, four of which are identified as ultraluminous X-ray sources (ULXs). The stellar environments of the X-ray sources are analyzed using the extinction-corrected color-magnitude diagrams and the isochrone models. We discover a highly asymmetric bubble nebula around X4 which exhibits different morphology in the H$α$ and [O III] images. The [O III]/H$α$ ratio map shows that the H$α$-bright bubble may be formed mainly via the shock ionization by the one-sided jet/outflow, while the more compact [O III] structure is photoionized by the ULX. We constrain the bubble expansion velocity and interstellar medium density with the MAPPINGS V code, and hence estimate the mechanical power injected to the bubble as $P_w \sim 5\times10^{40}$ erg s$^{-1}$ and the corresponding bubble age of $\sim7\times 10^{5}$ yr. Relativistic jets are needed to provide such level of mechanical power with a mass-loss rate of $\sim10^{-7}\ M_{\odot}\ \rm yr^{-1}$. Besides the accretion, the black hole spin is likely an additional energy source for the super-Eddington jet power. △ Less

Submitted 30 December, 2022; originally announced January 2023.

Comments: 17 pages, 10 figures, accepted by ApJ

Molecules in the peculiar age-defying source IRAS 19312+1950

Authors: Jian-Jie Qiu, Yong Zhang, Jun-ichi Nakashima, Jiang-Shui Zhang, Nico Koning, Xin-Di Tang, Yao-Ting Yan, Huan-Xue Feng

Abstract: Context. IRAS 19312+1950 is an isolated infrared source that exhibits a characteristic quasi-point-symmetric morphology in the near- and mid-infrared images and is also very bright in molecular radio lines. Because of its unique observational characteristics, various observational studies have been conducted and several hypotheses have been proposed regarding its origin, which is still unclear. So… ▽ More Context. IRAS 19312+1950 is an isolated infrared source that exhibits a characteristic quasi-point-symmetric morphology in the near- and mid-infrared images and is also very bright in molecular radio lines. Because of its unique observational characteristics, various observational studies have been conducted and several hypotheses have been proposed regarding its origin, which is still unclear. So far, it has been suggested that it could be a peculiar evolved star, a young stellar object, or even a red nova remnant. Regardless of which type of object it is ultimately classified as, IRAS 19312+1950 is exceptionally bright in the infrared and molecular radio lines and therefore will undoubtedly be crucial as a prototype of this kind of object having a peculiar nature or unusual evolutionary phase. Aims. This study aims to reveal the molecular composition of the central part of IRAS 19312+1950 by performing an unbiased molecular radio line survey and discussing the origin of the object from a molecular chemical point of view. Methods. We carried out a spectral line survey with the IRAM 30 m telescope towards the center of IRAS 19312+1950 in the 3 and 1.3 mm windows. Results. In total, 28 transition lines of 22 molecular species and those isotopologues are detected towards IRAS 19312+1950, some of which exhibit a broad and a narrow components. Seventeen thermal lines and 1 maser line are newly detected. The molecular species of C$^{17}$O, $^{30}$SiO, HN$^{13}$C, HC$^{18}$O$^{+}$, H$_{2}$CO, and $c$-C$_{3}$H$_{2}$ are detected for the first time in this object. Conclusions. Our results, in combination with previous studies, favor the hypothesis that IRAS 19312+1950 might be a red nova remnant, in which the progenitors that merged to become a red nova may have contained at least two evolved stars with oxygen-rich and carbon-rich chemistry, respectively. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: 30 pages, 24 figures, accepted for publication in A&A

Journal ref: A&A 669, A121 (2023)

Monte-Carlo simulations on possible collimation effects of outflows to fan-beamed emission of ultraluminous accreting X-ray pulsars

Authors: X. Hou, Y. You, L. Ji, R. Soria, S. N. Zhang, M. Y. Ge, L. Tao, S. Zhang, H. Feng, M. Zhou, Y. L. Tuo, L. M. Song, J. C. Wang

Abstract: Pulsating ultraluminous X-ray sources (PULXs) are accreting pulsars with apparent X-ray luminosity exceeding $10^{39}\, \rm erg\ s^{-1}$. We perform Monte-Carlo simulations to investigate whether high collimation effect (or strong beaming effect) is dominant in the presence of accretion outflows, for the fan beam emission of the accretion column of the neutron stars in PULXs. We show that the thre… ▽ More Pulsating ultraluminous X-ray sources (PULXs) are accreting pulsars with apparent X-ray luminosity exceeding $10^{39}\, \rm erg\ s^{-1}$. We perform Monte-Carlo simulations to investigate whether high collimation effect (or strong beaming effect) is dominant in the presence of accretion outflows, for the fan beam emission of the accretion column of the neutron stars in PULXs. We show that the three nearby PULXs (RX J0209.6$-$7427, Swift J0243.6+6124 and SMC X-3), namely the three musketeers here, have their main pulsed emission not strongly collimated even if strong outflows exist. This conclusion can be extended to the current sample of extragalactic PULXs, if accretion outflows are commonly produced from them. This means that the observed high luminosity of PULXs is indeed intrinsic, which can be used to infer the existence of very strong surface magnetic fields of $\sim10^{13-14}$ G, possibly multipole fields. However, if strong outflows are launched from the accretion disks in PULXs as a consequence of disk spherization by radiation pressure, regular dipole magnetic fields of $\sim10^{12}$ G may be required, comparable to that of the three musketeers, which have experienced large luminosity changes from well below their Eddington limit ($2\times10^{38}\, \rm erg\ s^{-1}$ for a NS) to super-Eddington and their maximum luminosity fills the luminosity gap between Galactic pulsars and extragalactic PULXs. △ Less

Submitted 15 November, 2022; originally announced November 2022.

Comments: 26 pages, 8 figures, 1 table. Accepted by ApJ

A Fast Transient Backend to Detect FRBs with the Tianlai Dish Pathfinder Array

Authors: Zijie Yu, Furen Deng, Shijie Sun, Chenhui Niu, Jixia Li, Fengquan Wu, Wei-Yang Wang, Yougang Wang, Hui Feng, Lin Shu, Jie Hao, Reza Ansari, Albert Stebbins, Xuelei Chen

Abstract: The Tianlai Dish Pathfinder array is a radio interferometer array consisting of 16 six meter dish antennas. The original digital backend integration time is at the seconds level, designed for HI intensity mapping experiment. A new digital backend with millisecond response is added to enable it to search for fast radio burst (FRB) during its observations. The design and calibration of this backend,… ▽ More The Tianlai Dish Pathfinder array is a radio interferometer array consisting of 16 six meter dish antennas. The original digital backend integration time is at the seconds level, designed for HI intensity mapping experiment. A new digital backend with millisecond response is added to enable it to search for fast radio burst (FRB) during its observations. The design and calibration of this backend, and the real time search pipeline for it are described in this paper. It is capable of forming 16 digital beams for each linear polarisation, covering an area of 19.6 square degrees. The search pipeline is capable of searching for, recording and classifying FRBs automatically in real time. In commissioning, we succeeded in capturing the signal pulses from the pulsars PSR B0329+54 and B2021+51. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: 16 pages, 14 figures, RAA accepted

Journal ref: Research in Astronomy and Astrophysics, 22, 125007 (2022)

Supermassive Black Hole and Broad-line Region in NGC~5548: Results from Five-season Reverberation Mapping

Authors: Kai-Xing Lu, Jin-Ming Bai, Jian-Min Wang, Chen Hu, Yan-Rong Li, Pu Du, Ming Xiao, Hai-Cheng Feng, Sha-Sha Li, Jian-Guo Wang, Zhi-Xiang Zhang, Ying-Ke Huang

Abstract: NGC 5548 is one of the active galactic nuclei (AGN) selected for our long-term spectroscopic monitoring with the Lijiang 2.4~m telescope, aiming at investigating the origin and evolution of the broad-line regions (BLRs), accurately measuring the mass of the supermassive black holes (SMBHs), and understanding structure and evolution of the AGN. We have performed five-season observations for NGC~554… ▽ More NGC 5548 is one of the active galactic nuclei (AGN) selected for our long-term spectroscopic monitoring with the Lijiang 2.4~m telescope, aiming at investigating the origin and evolution of the broad-line regions (BLRs), accurately measuring the mass of the supermassive black holes (SMBHs), and understanding structure and evolution of the AGN. We have performed five-season observations for NGC~5548 with the median sampling interval ranging from 1.25 to 3 days. The light curves of the 5100~Å continuum and broad emission lines are measured after subtracting contamination of the host galaxy starlight. The time lags of the broad He~{\sc ii}, He~{\sc i}, H$γ$, and H$β$ lines with respect to the 5100~Å continuum are obtained for each season and their mean time lags over the five seasons are 0.69, 4.66, 4.60, 8.43 days, respectively. The H$γ$ and H$β$ velocity-resolved lag profiles in the seasons of 2015, 2018, 2019, and 2021 are constructed, from which an ``M-shaped'' structure is found in 2015 but disappears after 2018. Our five-season reverberation mapping (RM) yields an averaged virial SMBH mass of $M_\bullet/10^7M_\odot=14.22$, with a small standard deviation of $1.89$. By combining the previous 18 RM campaigns and our five-season campaign for NGC~5548, we find that there exists a time lag of 3.5~years between the changes in the BLR size and optical luminosity. In addition, we also construct the BLR radius$-$luminosity relation and the virial relation for NGC~5548. △ Less

Submitted 23 September, 2022; v1 submitted 22 September, 2022; originally announced September 2022.

Comments: 23 pages, 15 figures, accepted for publication in ApJS

Fan beamed X-ray emission from 1 keV to above 130 keV from the ultraluminous X-ray pulsar RX J0209.6-7427 in the Small Magellanic Cloud

Authors: X. Hou, M. Y. Ge, L. Ji, S. N. Zhang, Y. You, L. Tao, S. Zhang, R. Soria, H. Feng, M. Zhou, Y. L. Tuo, L. M. Song, J. C. Wang

Abstract: We present detailed timing and spectral analyses of the transient X-ray pulsar RX J0209.6$-$7427 in the Small Magellanic Cloud during its 2019 giant outburst. With a better known distance than most galactic X-ray pulsars, its peak luminosity is determined to be $(1.11\pm0.06)\times 10^{39}\, \rm erg\ s^{-1}$; it is thus a {\it bonda fide} pulsating ultraluminous X-ray source (PULX). Owing to the b… ▽ More We present detailed timing and spectral analyses of the transient X-ray pulsar RX J0209.6$-$7427 in the Small Magellanic Cloud during its 2019 giant outburst. With a better known distance than most galactic X-ray pulsars, its peak luminosity is determined to be $(1.11\pm0.06)\times 10^{39}\, \rm erg\ s^{-1}$; it is thus a {\it bonda fide} pulsating ultraluminous X-ray source (PULX). Owing to the broad energy band of \textit{Insight}-HXMT, its pulsed X-ray emission was detected from 1 keV up to the 130$-$180 keV band, which is the highest energy emission detected from any PULXs outside the Milky Way. This allows us to conclude that its main pulsed X-ray emission is from the "fan beam" of the accretion column, and its luminosity is thus intrinsic. We also estimate its magnetic field of (4.8$-$8.6)$\times10^{12}$ G or (1.7$-$2.2)$\times10^{13}$ G, from its spin evolution or transition in the accretion column structure during the outburst; we suggest that the two values of the magnetic field strength correspond to the dipole and multipole magnetic fields of the neutron star, similar to the recent discovery in the Galactic PULX Swift J0243.6+6124. Therefore, the nature of the neutron star and its ULX emission can be understood within the current theoretical frame of accreting neutron stars. This may have implications for understanding the nature of those farther away extragalactic PULXs. △ Less

Submitted 31 August, 2022; originally announced August 2022.

Comments: 31 pages, 14 figures, 5 tables. Accepted by ApJ

Velocity-resolved Reverberation Mapping of Changing-look Active Galactic Nucleus NGC~4151 During Outburst Stage: Evidence for Kinematics Evolution of Broad-line Region

Authors: Sha-Sha Li, Hai-Cheng Feng, H. T. Liu, J. M. Bai, Rui Li, Kai-Xing Lu, Jian-Guo Wang, Ying-Ke Huang, Zhi-Xiang Zhang

Abstract: Changing-look active galactic nucleus NGC~4151, which has attracted a lot of attention, is undergoing the second dramatic outburst stage in its evolutionary history. To investigate the geometry and kinematics of the broad-line region (BLR), and measure the mass of supermassive black hole in NGC~4151, we perform a seven-month photometric and spectroscopic monitoring program in 2020--2021, using the… ▽ More Changing-look active galactic nucleus NGC~4151, which has attracted a lot of attention, is undergoing the second dramatic outburst stage in its evolutionary history. To investigate the geometry and kinematics of the broad-line region (BLR), and measure the mass of supermassive black hole in NGC~4151, we perform a seven-month photometric and spectroscopic monitoring program in 2020--2021, using the 2.4 m telescope at Lijiang Observatory. We successfully measure the time lags of the responses from broad \ha, \hb, \hg, \hei, and \heii\ emission lines to continuum variation, which are $7.63_{-2.62}^{+1.85}$, $6.21_{-1.13}^{+1.41}$, $5.67_{-1.94}^{+1.65}$, $1.59_{-1.11}^{+0.86}$, and $0.46_{-1.06}^{+1.22}$ days, respectively, following radial stratification. The ratios of time lags among these lines are $1.23 : 1.00 : 0.91 : 0.26 : 0.07$. We find that the continuum lag between the ultraviolet and optical bands can significantly affect the lag measurements of \hei\ and \heii. Virial and infalling gas motions coexist in this campaign, which is different from previous results, implying the evolutionary kinematics of BLR. Based on our measurements and previous ones in the literature, we confirm that the BLR of NGC~4151 is basically virialized. Finally, we compute the black hole mass through multiple lines, and the measurement from \hb\ to be $ 3.94_{-0.72}^{+0.90} \times 10^7 M_{\odot}$, which is consistent with previous results. The corresponding accretion rate is $0.02_{-0.01}^{+0.01} L_{\rm Edd} c^{-2}$, implying a sub-Eddington accretor. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: 16 pages, 9 figures, accepted for publication in ApJ

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

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

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

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

Comments: 18 pages, 11 figures

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

VLT MUSE observations of the bubble nebula around NGC 1313 X-2 and evidence for additional photoionization

Authors: Changxing Zhou, Fuyan Bian, Hua Feng, Jiahui Huang

Abstract: The bubble nebula surrounding NGC 1313 X-2 is believed to be powered by high velocity winds from the central ultraluminous X-ray source (ULX) as a result of supercritical accretion. With the Multi-Unit Spectroscopic Explorer (MUSE) observation of the nebula, we find enhanced OIII emission at locations spatially coincident with clusters of stars and the central X-ray source, suggesting that photoio… ▽ More The bubble nebula surrounding NGC 1313 X-2 is believed to be powered by high velocity winds from the central ultraluminous X-ray source (ULX) as a result of supercritical accretion. With the Multi-Unit Spectroscopic Explorer (MUSE) observation of the nebula, we find enhanced OIII emission at locations spatially coincident with clusters of stars and the central X-ray source, suggesting that photoionization in addition to shock-ionization plays an important role in powering the nebula. The X-ray luminosity of the ULX and the number of massive stars in the nebula region can account for the required ionizing luminosity derived with MAPPINGS V, which also confirms that pure shocks cannot explain the observed emission line ratios. △ Less

Submitted 13 July, 2022; originally announced July 2022.

Comments: ApJ to appear

A Catalog of Molecular Clumps and Cores with Infall Signatures

Authors: Shuling Yu, Zhibo Jiang, Yang Yang, Zhiwei Chen, Haoran Feng

Abstract: The research of infall motion is a common means to study molecular cloud dynamics and the early process of star formation. Many works had been done in-depth research on infall. We searched the literature related to infall study of molecular cloud since 1994, summarized the infall sources identified by the authors. A total of 456 infall sources are catalogued. We classify them into high-mass and lo… ▽ More The research of infall motion is a common means to study molecular cloud dynamics and the early process of star formation. Many works had been done in-depth research on infall. We searched the literature related to infall study of molecular cloud since 1994, summarized the infall sources identified by the authors. A total of 456 infall sources are catalogued. We classify them into high-mass and low-mass sources, in which the high-mass sources are divided into three evolutionary stages: prestellar, protostellar and HII region. We divide the sources into clumps and cores according to their sizes. The H$_2$ column density values range from 1.21$\times$ 10$^{21}$ to 9.75 $\times$ 10$^{24}$ cm$^{-2}$, with a median value of 4.17$\times$ 10$^{22}$ cm$^{-2}$. The H$_2$ column densities of high-mass and low-mass sources are significantly separated. The median value of infall velocity for high-mass clumps is 1.12 km s$^{-1}$, and the infall velocities of low-mass cores are virtually all less than 0.5 km s$^{-1}$. There is no obvious difference between different stages of evolution. The mass infall rates of low-mass cores are between 10$^{-7}$ and 10$^{-4}$ M$_{\odot} \text{yr}^{-1}$, and those of high-mass clumps are between 10$^{-4}$ and 10$^{-1}$ M$_{\odot} \text{yr}^{-1}$ with only one exception. We do not find that the mass infall rates vary with evolutionary stages. △ Less

Submitted 29 June, 2022; v1 submitted 28 June, 2022; originally announced June 2022.