Showing 1–50 of 479 results for author: Chen, L

Diffuse X-ray Explorer: a high-resolution X-ray spectroscopic sky surveyor on the China Space Station

Authors: Hai Jin, Junjie Mao, Liubiao Chen, Naihui Chen, Wei Cui, Bo Gao, Jinjin Li, Xinfeng Li, Jiejia Liu, Jia Quan, Chunyang Jiang, Guole Wang, Le Wang, Qian Wang, Sifan Wang, Aimin Xiao, Shuo Zhang

Abstract: DIffuse X-ray Explorer (DIXE) is a proposed high-resolution X-ray spectroscopic sky surveyor on the China Space Station (CSS). DIXE will focus on studying hot baryons in the Milky Way. Galactic hot baryons like the X-ray emitting Milky Way halo and eROSITA bubbles are best observed in the sky survey mode with a large field of view. DIXE will take advantage of the orbital motion of the CSS to scan… ▽ More DIffuse X-ray Explorer (DIXE) is a proposed high-resolution X-ray spectroscopic sky surveyor on the China Space Station (CSS). DIXE will focus on studying hot baryons in the Milky Way. Galactic hot baryons like the X-ray emitting Milky Way halo and eROSITA bubbles are best observed in the sky survey mode with a large field of view. DIXE will take advantage of the orbital motion of the CSS to scan a large fraction of the sky. High-resolution X-ray spectroscopy, enabled by superconducting microcalorimeters based on the transition-edge sensor (TES) technology, will probe the physical properties (e.g., temperature, density, elemental abundances, kinematics) of the Galactic hot baryons. This will complement the high-resolution imaging data obtained with the eROSITA mission. Here we present the preliminary design of DIXE. The payload consists mainly of a detector assembly and a cryogenic cooling system. The key components of the detector assembly are a microcalorimeter array and frequency-domain multiplexing readout electronics. To provide a working temperature for the detector assembly, the cooling system consists of an adiabatic demagnetization refrigerator and a mechanical cryocooler system. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Comments: 12 pages, 6 figures, the full version is published by Journal of Low Temperature Physics

Interstellar Nitrogen Isotope Ratios: Measurements on tracers of C$^{14}$N and C$^{15}$N

Authors: J. L. Chen, J. S. Zhang, C. Henkel, Y. T. Yan, H. Z. Yu, Y. X. Wang, Y. P. Zou, J. Y. Zhao, X. Y. Wang

Abstract: The nitrogen isotope ratio 14N/15N is a powerful tool to trace Galactic stellar nucleosynthesis and constraining Galactic chemical evolution. Previous observations have found lower 14N/15N ratios in the Galactic center and higher values in the Galactic disk. This is consistent with the inside-out formation scenario of our Milky Way. However, previous studies mostly utilized double isotope ratios a… ▽ More The nitrogen isotope ratio 14N/15N is a powerful tool to trace Galactic stellar nucleosynthesis and constraining Galactic chemical evolution. Previous observations have found lower 14N/15N ratios in the Galactic center and higher values in the Galactic disk. This is consistent with the inside-out formation scenario of our Milky Way. However, previous studies mostly utilized double isotope ratios also including 12C/13C, which introduces additional uncertainties. Here we therefore present observations of C14N and its rare isotopologue, C15N, toward a sample of star forming regions, measured by the IRAM 30 m and/or the ARO 12 m telescope at $λ$ ~3 mm wavelength. For those 35 sources detected in both isotopologues, physical parameters are determined. Furthermore we have obtained nitrogen isotope ratios using the strongest hyperfine components of CN and C15N. For those sources showing small deviations from Local Thermodynamical Equilibrium and/or self-absorption, the weakest hyperfine component, likely free of the latter effect, was used to obtain reliable 14N/15N values. Our measured 14N/15N isotope ratios from C14N and C15N measurements are compatible with those from our earlier measurements of NH3 and 15NH3 (Paper I), i.e., increasing ratios to a Galacticentric distance of ~9 kpc. The unweighted second order polynomial fit yields $\frac{{\rm C^{14}N}}{{\rm C^{15}N}} = (-4.85 \pm 1.89)\;{\rm kpc^{-2}} \times R_{\rm GC}^{2} + (82.11 \pm 31.93) \;{\rm kpc^{-1}} \times R_{\rm GC} - (28.12 \pm 126.62)$. Toward the outer galaxy, the isotope ratio tends to decrease, supporting an earlier finding by H13CN/HC15N. Galactic chemical evolution models are consistent with our measurements of the 14N/15N isotope ratio, i.e. a rising trend from the Galactic center region to approximately 9 kpc, followed by a decreasing trend with increasing $R_{\rm GC}$ toward the outer Galaxy. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 34 pages, 9 figures, 6 tables

Journal ref: The Astrophysical Journal (2004)

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

Simulation Models for Exploring Magnetic Reconnection

Authors: Michael Shay, Subash Adhikari, Naoki Beesho, Joachim Birn, Jorg Buechner, Paul Cassak, Li-Jen Chen, Yuxi Chen, Giulia Cozzani, Jim Drake, Fan Guo, Michael Hesse, Neeraj Jain, Yann Pfau-Kempf, Yu Lin, Yi-Hsin Liu, Mitsuo Oka, Yuri A. Omelchenko, Minna Palmroth, Oreste Pezzi, Patricia H. Reiff, Marc Swisdak, Frank Toffoletto, Gabor Toth, Richard A. Wolf

Abstract: Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This chapter reviews many of these… ▽ More Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This chapter reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM). △ Less

Submitted 9 June, 2024; originally announced June 2024.

Comments: Chapter 5.2 of ISSI Book on Magnetic Reconnection, submitted to Space Science Reviews

PREX and CREX: Evidence for Strong Isovector Spin-Orbit Interaction

Authors: Tong-Gang Yue, Zhen Zhang, Lie-Wen Chen

Abstract: The recent PREX-2 and CREX data on the model-independent extraction of the charge-weak form factor difference $ΔF_{\rm CW}$ in $^{208}$Pb and $^{48}$Ca challenge modern nuclear energy density functionals (EDFs) as well as our present understanding on the neutron skin and nuclear symmetry energy. Within the Skyrme-like EDFs, we demonstrate that the isovector spin-orbit interaction can strongly chan… ▽ More The recent PREX-2 and CREX data on the model-independent extraction of the charge-weak form factor difference $ΔF_{\rm CW}$ in $^{208}$Pb and $^{48}$Ca challenge modern nuclear energy density functionals (EDFs) as well as our present understanding on the neutron skin and nuclear symmetry energy. Within the Skyrme-like EDFs, we demonstrate that the isovector spin-orbit interaction can strongly change the $ΔF_{\rm CW}$ in $^{48}$Ca while it has essentially no influence on the $ΔF_{\rm CW}$ in $^{208}$Pb, mainly due to the eight spin-orbit unpaired $1f_{7/2}$ neutrons in $^{48}$Ca. To simultaneously describe PREX-2 and CREX data in $1σ$ error, we find the strength of isovector spin-orbit interaction should be larger than about four times of that in the conventional Skyrme-like EDFs, implying the neutrons and protons have significantly different spin-orbit interaction. To further reconcile the data on electric dipole polarizability in $^{208}$Pb and $^{48}$Ca, we obtain $L \approx 55$ MeV for the slope parameter of the symmetry energy, $Δr_{\rm np}(^{208}\rm{Pb}) \approx 0.19$ fm and $Δr_{\rm np}(^{48}\rm{Ca}) \approx 0.12$ fm for the neutron skin thickness. The implications of the strong isovector spin-orbit interaction are discussed. △ Less

Submitted 6 June, 2024; originally announced June 2024.

Comments: 14 pages, 6 figures (including Supplemental Material)

A $γ$-Ray Emitting Blazar at Redshift 3.64: Fermi-LAT and OVRO Observations of PKS 0201+113

Authors: Hai Lei, Ying-Kang Zhang, Xiong Jiang, Sebastian Kiehlmann, Anthony C. S. Readhead, Liang Chen, Neng-Hui Liao, Tao An

Abstract: High-redshift ($z >3$) $γ$-ray blazars are rare, but they are crucial for our understanding of jet evolution, $γ$-ray production and propagation, and the growth of supermassive black holes in the early universe. A new analysis of Fermi-LAT data reveals a significant (5$σ$), spectrally soft ($Γ\simeq$ 3.0) $γ$-ray source in a specific 4-month epoch, cospatial with PKS 0201+113 ($z$ = 3.64). Monitor… ▽ More High-redshift ($z >3$) $γ$-ray blazars are rare, but they are crucial for our understanding of jet evolution, $γ$-ray production and propagation, and the growth of supermassive black holes in the early universe. A new analysis of Fermi-LAT data reveals a significant (5$σ$), spectrally soft ($Γ\simeq$ 3.0) $γ$-ray source in a specific 4-month epoch, cospatial with PKS 0201+113 ($z$ = 3.64). Monitoring of PKS 0201+113 at 15 GHz by the Owens Valley Radio Observatory 40 m Telescope from 2008 to 2023 shows a prominent flare that dominates the radio light curve. The maximum of the radio flare coincides with the $γ$-ray flare, strongly suggesting an association ($\textrm{p-value}=0.023$) between the $γ$-ray and the radio sources. PKS 0201+113 is only the third $γ$-ray blazar to be identified with $z> 3.5$, and it is the first such object to be identified by the detection of quasi-simultaneous $γ$-ray and radio flares. The jet properties of this peculiar blazar have been investigated. A detailed study of a two-zone leptonic model is presented that fits the broadband spectral energy distribution. An alternative scenario is also briefly discussed. △ Less

Submitted 5 June, 2024; originally announced June 2024.

Comments: 22 pages, 6 figures, 2 Tables, accepted for publication in ApJ

On the Determination of Stellar Mass and Binary Fraction of Open Clusters within 500 pc from the Sun

Authors: Yueyue Jiang, Jing Zhong, Songmei Qin, Tong Tang, Li Chen, Jinliang Hou

Abstract: We investigated the stellar mass function and the binary fraction of 114 nearby open clusters (OCs) using the high-precision photometric data from Gaia Data Release 3 (Gaia DR3). We estimated the mass of member stars by using a ridge line (RL) that is better in line with the observed color-magnitude diagram (CMD), thus obtaining more accurate stellar mass and binary mass ratio ($q$) at the low-mas… ▽ More We investigated the stellar mass function and the binary fraction of 114 nearby open clusters (OCs) using the high-precision photometric data from Gaia Data Release 3 (Gaia DR3). We estimated the mass of member stars by using a ridge line (RL) that is better in line with the observed color-magnitude diagram (CMD), thus obtaining more accurate stellar mass and binary mass ratio ($q$) at the low-mass region. By analyzing the present-day mass function (PDMF) of star clusters, we found that 108 OCs follow a two-stage power-law distribution, whereas 6 OCs present a single power-law PDMF. Subsequently, we fitted the high(low)-mass index of PDMF ($dN/dm \propto m^{-α}$), denoted as $α_{\rm h}$($α_{\rm l}$), and segmentation point $m_{\rm c}$. For our cluster sample, the median values of $α_{\rm h}$ and $α_{\rm l}$ are 2.65 and 0.95, respectively, which are approximately consistent with the initial mass function (IMF) results provided by Kroupa (2001). We utilized the cumulative radial number distribution of stars with different masses to quantify the degree of mass segregation. We found a significant positive correlation between the state of dynamical evolution and mass segregation in OCs. We also estimated the fraction of binary stars with $q \geq 0.5$, ranging from 6% to 34% with a median of 17%. Finally, we provided a catalog of 114 nearby cluster properties, including the total mass, the binary fraction, the PDMF, and the dynamical state. △ Less

Submitted 21 May, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 20 pages, 15 figures, accepted for publication in ApJ

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 Red Supergiant Progenitor of Type II Supernova 2024ggi

Authors: Danfeng Xiang, Jun Mo, Xiaofeng Wang, Lingzhi Wang, Jujia Zhang, Han Lin, Liyang Chen, Cuiying Song, Liang-Duan Liu, Zhenyu Wang, Gaici Li

Abstract: We present a detailed analysis of the progenitor and its local environment for the recently discovered type II supernova (SN) 2024ggi at a distance of about 6.7~Mpc, by utilizing the pre-explosion images from the Hubble Space Telescope (HST) and \textit{Spitzer} Space Telescope. The progenitor is identified as a red, bright variable star, with absolute $F814W$-band magnitudes being $-$6.2 mag in 1… ▽ More We present a detailed analysis of the progenitor and its local environment for the recently discovered type II supernova (SN) 2024ggi at a distance of about 6.7~Mpc, by utilizing the pre-explosion images from the Hubble Space Telescope (HST) and \textit{Spitzer} Space Telescope. The progenitor is identified as a red, bright variable star, with absolute $F814W$-band magnitudes being $-$6.2 mag in 1995 to $-$7.2 mag in 2003, respectively, consistent with that of a normal red supergiant (RSG) star. Combining with the historical mid-infrared light curves, a pulsational period of about 379~days can be inferred for the progenitor star. Fitting its spectral energy distribution with stellar spectral models yields the stellar parameters of temperature, radius and bolometric luminosity as $T_*=3290_{-27}^{+19}$~K, $R_*=887_{-51}^{+60}$~R$_{\odot}$, and log($L$/L$_{\odot}$)$=4.92_{-0.04}^{+0.05}$, respectively. The above parameters indicate that the progenitor of SN 2024ggi is consistent with the stellar evolutionary track of a solar-metallicity massive star with an initial mass of $13_{-1}^{+1}$~M$_{\odot}$. Moreover, our analysis indicates a relatively low mass loss rate (i.e., $< 3\times10^{-6}$~M$_{\odot}$~yr$^{-1}$) for the progenitor compared to that inferred from the flashed spectra and X-ray detection (i.e., $10^{-2}$$-$$ 10$$^{-5}$~M$_{\odot}$~yr$^{-1}$), implying a significant enhancement in mass loss within a few years prior to the explosion. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 4 figures, 1 table

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

Exploration of morphological coherence in open clusters with a "core-shell'' structure

Authors: Qingshun Hu, Yu Zhang, Songmei Qin, Jing Zhong, Li Chen, Yangping Luo

Abstract: The study of their morphological coherence allows us to obtain a better understanding of the morphological evolution of open clusters. We employed the ellipsoid fitting method to delineate the 3D spatial structure of the sample clusters while using the morphological dislocation (MD) defined in our previous work and the ellipticity ratio (ER) of the clusters' inner and outer structures to character… ▽ More The study of their morphological coherence allows us to obtain a better understanding of the morphological evolution of open clusters. We employed the ellipsoid fitting method to delineate the 3D spatial structure of the sample clusters while using the morphological dislocation (MD) defined in our previous work and the ellipticity ratio (ER) of the clusters' inner and outer structures to characterize the morphological coherence of the sample clusters. The results show an inverse correlation between the ER of the sample clusters and the number of their members, indicating that sample clusters with a much more elliptical external morphology than internal shape generally tend to host a large number of members. Meanwhile, a slight shrinking of the MD of the sample clusters with their members' number may shed light on the significant role of the gravitational binding of the sample clusters in maintaining their morphological stability. Moreover, there are no correlations between the MD and ER of the sample clusters and their age. They are also not significantly correlated with the X-axis, the Y-axis, their orbital eccentricities, and the radial and vertical forces on them. However, the ER of the sample clusters displays some fluctuations in the distributions between it and the above covariates, implying that the morphologies of the sample clusters are sensitive to the external environment if sample effects are not taken into account. Finally, the analysis of the 3D spatial shapes of sample clusters with a small ER or a large ER demonstrates that the number of members lays an important foundation for forming a dense internal system for sample clusters. At the same time, the MD of the sample clusters can serve well as an indicator of their morphological stability, which is built upon a certain amount of member stars. △ Less

Submitted 11 May, 2024; v1 submitted 9 May, 2024; originally announced May 2024.

Comments: 17 pages, 13 figures, 1 table. Accepted by Astronomy & Astrophysics

Magnetically Driven Relativistic Jet in the High-Redshift Blazar OH~471

Authors: S. Guo, T. An, Y. Liu, Y. Sotnikova, A. Volvach, T. Mufakharov, L. Chen, L. Cui, A. Wang, Z. Xu, Y. Zhang, W. Xu, Y. A. Kovalev, Y. Y. Kovalev, M. Kharinov, A. Erkenov, T. Semenova, L. Volvach

Abstract: Context : Understanding the mechanisms that launch and shape powerful relativistic jets from supermassive black holes (SMBHs) in high-redshift active galactic nuclei (AGN) is crucial for probing the co-evolution of SMBHs and galaxies over cosmic time. Aims :We study the high-redshift ($z=3.396$) blazar OH~471 to explore the jet launching mechanism in the early Universe. Methods : Using multi-f… ▽ More Context : Understanding the mechanisms that launch and shape powerful relativistic jets from supermassive black holes (SMBHs) in high-redshift active galactic nuclei (AGN) is crucial for probing the co-evolution of SMBHs and galaxies over cosmic time. Aims :We study the high-redshift ($z=3.396$) blazar OH~471 to explore the jet launching mechanism in the early Universe. Methods : Using multi-frequency radio monitoring observations and high-resolution Very Long Baseline Interferometry imaging over three decades, we study the milliarcsecond structure and long-term variability of OH~471. Results : Spectral modelling of the radio flux densities reveals a synchrotron self-absorbed spectrum indicating strong magnetic fields within the compact core. By applying the flux freezing approximation, we estimate the magnetic flux carried by the jet and find that it reaches or exceeds theoretical predictions for jets powered by black hole spin energy via the Blandford-Znajek mechanism. This implies that OH~471 was in a magnetically arrested disk (MAD) state where the magnetic flux accumulated near the horizon regulates the accretion flow, allowing efficient extraction of black hole rotational energy. Conclusions : Our study demonstrates the dominance of MAD accretion in powering the prominent radio flares and relativistic jets observed in the radio-loud AGN OH~471 and statistical studies of large samples of high-redshift AGN will shed light on the role of MAD accretion in launching and accelerating the earliest relativistic jets. △ Less

Submitted 20 May, 2024; v1 submitted 25 April, 2024; originally announced April 2024.

Comments: 16 pages, 7 figures, 3 tables

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

LHAASO-KM2A detector simulation using Geant4

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

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

Submitted 7 April, 2024; originally announced April 2024.

A fast-rotating blue straggler star in the tidal tail of the open cluster NGC 752

Authors: Chunyan Li, Jing Zhong, Songmei Qin, Dengkai Jiang, Xingmei Shan, Li Chen

Abstract: NGC 752 is a famous Galactic open cluster of intermediate age. In recent works, a very long and asymmetric tail was newly revealed. A blue straggler star (BSS) at the periphery of the tidal tail of the cluster has been identified subsequently. We aim to perform a detailed analysis of the newly detected BSS based on the available comprehensive spectroscopic and photometric data. We also explored th… ▽ More NGC 752 is a famous Galactic open cluster of intermediate age. In recent works, a very long and asymmetric tail was newly revealed. A blue straggler star (BSS) at the periphery of the tidal tail of the cluster has been identified subsequently. We aim to perform a detailed analysis of the newly detected BSS based on the available comprehensive spectroscopic and photometric data. We also explored this BSS's possible formation pathway and age limitation based on the collected spectroscopic and photometric data. We estimated the projected rotational velocity $v\ \mathrm{sin}i$ and the mass of the BSS from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope low-resolution spectra and multiband photometric data from various catalogs, respectively. The newly discovered BSS is confirmed as a genuine member of NGC 752. The lack of ultraviolet excess in the SED and no significant variations in the light curve imply that this BSS is likely a single star ($mass=1.86^{+3.62}_{-0.94}\ M_{\odot}$) formed through stellar mergers. The fast rotation velocity ($v\ \mathrm{sin}i=206.9\pm4.9$~km $\rm s^{-1}$) of the BSS may provide constraints on its age (less than a hundred million years), but more formation details require further investigation. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 7 pages, 6 figures. Accepted by A&A

Journal ref: A&A 686, A215 (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)

Type IV-like Solar Radio Burst Consisting of a Series of Spikes Observed by PSP

Authors: Bing Ma, Ling Chen, De-Jin Wu, Marc Pulupa, Stuart D. Bale

Abstract: Solar and interplanetary radio bursts can reflect the existence and motion of energetic electrons and are therefore a kind of vital phenomenon in solar activities. The present study reported a solar radio burst (SRB) event observed by Parker Solar Probe (PSP) in its 8th orbital encounter phase, and it lasted about 20 hours in a frequency range of 0.5-15 MHz, called the type IV-like SRB. This type… ▽ More Solar and interplanetary radio bursts can reflect the existence and motion of energetic electrons and are therefore a kind of vital phenomenon in solar activities. The present study reported a solar radio burst (SRB) event observed by Parker Solar Probe (PSP) in its 8th orbital encounter phase, and it lasted about 20 hours in a frequency range of 0.5-15 MHz, called the type IV-like SRB. This type IV-like SRB consists of a series of numerous spikes with the center-frequency drifting slowly from ~5 MHz to ~1 MHz, and each individual spike appears a much faster frequency drifting and has a narrow frequency range of a few MHz and short duration of a few minutes. Based on the empirical models of the solar atmosphere adopted commonly, combining the in-situ measurement by PSP, we propose that these small-scale spikes were generated by a group of solitary kinetic Alfvén waves (SKAWs) in a magnetic loop accompanying coronal mass ejection (CME) and moving outwards, in which the frequency drifting of individual spike is caused by the SKAW's propagation and the center-frequency drifting may be attributed to the motion of the magnetic loop. △ Less

Submitted 16 June, 2024; v1 submitted 12 March, 2024; originally announced March 2024.

Comments: There are some questions about models and the emission mechanisms to be discussed more carefully. We need to revise this manuscript

Reconstructions of Jupiter's magnetic field using physics informed neural networks

Authors: Philip W. Livermore, Leyuan Wu, Longwei Chen, Sjoerd A. L. de Ridder

Abstract: Magnetic sounding using data collected from the Juno mission can be used to provide constraints on Jupiter's interior. However, inwards continuation of reconstructions assuming zero electrical conductivity and a representation in spherical harmonics are limited by the enhancement of noise at small scales. Here we describe new reconstructions of Jupiter's internal magnetic field based on physics-in… ▽ More Magnetic sounding using data collected from the Juno mission can be used to provide constraints on Jupiter's interior. However, inwards continuation of reconstructions assuming zero electrical conductivity and a representation in spherical harmonics are limited by the enhancement of noise at small scales. Here we describe new reconstructions of Jupiter's internal magnetic field based on physics-informed neural networks and either the first 33 (PINN33) or the first 50 (PINN50) of Juno's orbits. The method can resolve local structures, and allows for weak ambient electrical currents. Our models are not hampered by noise amplification at depth, and offer a much clearer picture of the interior structure. We estimate that the dynamo boundary is at a fractional radius of 0.8. At this depth, the magnetic field is arranged into longitudinal bands, and strong local features such as the great blue spot appear to be rooted in neighbouring structures of oppositely signed flux. △ Less

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

Astrochemical effect of the fundamental grain surface processes I. The diffusion of grain surface species and the pre-exponential factor

Authors: Long-Fei Chen, Donghui Quan, Jiao He, Yao Wang, Di Li, Thomas Henning

Abstract: Abbreviation. Thermal diffusion is one of the basic processes for the mobility and formation of species on cosmic dust grains. Recent laboratory measurements have found that the diffusion pre-exponential factor can differ from that for desorption by several orders of magnitude. We aim to evaluate the effect of the newly experimentally measured diffusion pre-exponential factor on the chemistry unde… ▽ More Abbreviation. Thermal diffusion is one of the basic processes for the mobility and formation of species on cosmic dust grains. Recent laboratory measurements have found that the diffusion pre-exponential factor can differ from that for desorption by several orders of magnitude. We aim to evaluate the effect of the newly experimentally measured diffusion pre-exponential factor on the chemistry under cold molecular cloud conditions. We found that statistically, more than half of the total gas-phase and grain surface species are not affected by the new pre-exponential factor after a chemical evolution of 10$^5$ yr. The most abundant gas-phase CO and grain surface water ice are not affected by the new pre-exponential factor. For the grain surface species that are affected, compared to the commonly adopted value of the pre-exponential factor for diffusion used in the chemical models, they could be either overproduced or underproduced with the lower diffusion pre-factor used in this work. The former case applies to radicals and the species that serve as reactants, while the latter case applies to complex organic molecules (COMs) on the grain and the species that rarely react with other species. Gas-phase species could also be affected due to the desorption of the grain surface species. The abundance of some gas-phase COMs could be varied by over one order of magnitude depending on the adopted grain surface temperature and/or the ratio of diffusion to desorption energy in the model. Key species whose diffusion pre-exponential factor significantly affects the model predictions were also evaluated, and these specie include CH3OH, H2CO, and NO. The results presented in this study show that the pre-exponential factor is one of the basic and important parameters in astrochemical models. △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 14 pages, 15 figures, 2 tables, accepted for publication in A&A

Journal ref: A&A, 2024, 685, A55

Field-Aligned Current Structures during the Terrestrial Magnetosphere's Transformation into Alfven Wings and Recovery

Authors: Jason M. H. Beedle, Li-Jen Chen, Jason R. Shuster, Harsha Gurram, Dan J. Gershman, Yuxi Chen, Rachel C. Rice, Brandon L. Burkholder, Akhtar S. Ardakani, Kevin J. Genestreti, Roy B. Torbert

Abstract: On April 24th, 2023, a CME event caused the solar wind to become sub-Alfvenic, leading to the development of an Alfven Wing configuration in the Earth's Magnetosphere. Alfven Wings have previously been observed as cavities of low flow in Jupiter's magnetosphere, but the observing satellites did not have the ability to directly measure the Alfven Wings' current structures. Through in situ measureme… ▽ More On April 24th, 2023, a CME event caused the solar wind to become sub-Alfvenic, leading to the development of an Alfven Wing configuration in the Earth's Magnetosphere. Alfven Wings have previously been observed as cavities of low flow in Jupiter's magnetosphere, but the observing satellites did not have the ability to directly measure the Alfven Wings' current structures. Through in situ measurements made by the Magnetospheric Multiscale (MMS) spacecraft, the April 24th event provides us with the first direct measurements of current structures during an Alfven Wing configuration. We have found two distinct types of current structures associated with the Alfven Wing transformation as well as the magnetosphere recovery. These structures are observed to be significantly more anti-field-aligned and electron-driven than typical magnetopause currents, indicating the disruptions caused to the magnetosphere current system by the Alfven Wing formation. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Abundant sub-micron grains revealed in newly discovered extreme debris discs

Authors: A. Moór, P. Ábrahám, K. Y. L. Su, T. Henning, S. Marino, L. Chen, Á. Kóspál, N. Pawellek, J. Varga, K. Vida

Abstract: Extreme debris discs (EDDs) are bright and warm circumstellar dusty structures around main sequence stars. They may represent the outcome of giant collisions occuring in the terrestrial region between large planetesimals or planetary bodies, and thus provide a rare opportunity to peer into the aftermaths of these events. Here, we report on results of a mini-survey we conducted with the aim to incr… ▽ More Extreme debris discs (EDDs) are bright and warm circumstellar dusty structures around main sequence stars. They may represent the outcome of giant collisions occuring in the terrestrial region between large planetesimals or planetary bodies, and thus provide a rare opportunity to peer into the aftermaths of these events. Here, we report on results of a mini-survey we conducted with the aim to increase the number of known EDDs, investigate the presence of solid-state features around 10μm in eight EDDs, and classify them into the silica or silicate dominated groups. We identify four new EDDs and derive their fundamental properties. For these, and for four other previously known discs, we study the spectral energy distribution around 10μm by means of VLT/VISIR photometry in three narrow-band filters and conclude that all eight objects likely exhibit solid-state emission features from sub-micron grains. We find that four discs probably belong to the silicate dominated subgroup. Considering the age distribution of the entire EDD sample, we find that their incidence begins to decrease only after 300 Myr, suggesting that the earlier common picture that these objects are related to the formation of rocky planets may not be exclusive, and that other processes may be involved for older objects (>100 Myr). Because most of the older EDD systems have wide, eccentric companions, we suggest that binarity may play a role in triggering late giant collisions. △ Less

Submitted 23 February, 2024; originally announced February 2024.

Comments: 20 pages, 9 figures, published in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 528, Issue 3, March 2024, pp.4528-4546

Is it possible to detect coronal mass ejections on solar-type stars through extreme-ultraviolet spectral observations?

Authors: Zihao Yang, Hui Tian, Yingjie Zhu, Yu Xu, Linyi Chen, Zheng Sun

Abstract: Stellar coronal mass ejections (CMEs) from host stars are an important factor that affects the habitability of exoplanets. Although their solar counterparts have been well observed for decades, it is still very difficult to find solid evidence of stellar CMEs. Using the spectral line profile asymmetry caused by the Doppler shift of erupting plasma, several stellar CME candidates have been identifi… ▽ More Stellar coronal mass ejections (CMEs) from host stars are an important factor that affects the habitability of exoplanets. Although their solar counterparts have been well observed for decades, it is still very difficult to find solid evidence of stellar CMEs. Using the spectral line profile asymmetry caused by the Doppler shift of erupting plasma, several stellar CME candidates have been identified from spectral lines formed at chromospheric or transition region temperatures of the stars. However, a successful detection of stellar CME signals based on the profile asymmetries of coronal lines is still lacking. It is unclear whether we can detect such signals. Here we construct an analytical model for CMEs on solar-type stars, and derive an expression of stellar extreme-ultraviolet (EUV) line profiles during CME eruptions. For different instrumental parameters, exposure times, CME propagation directions and stellar activity levels, we synthesized the corresponding line profiles of Fe IX 171.07 Å and Fe XV 284.16 Å. Further investigations provide constraints on the instrumental requirements for successful detection and characterization of stellar CMEs. Our results show that it is possible to detect stellar CME signals and infer their velocities based on spectral profile asymmetries using an EUV spectrometer with a moderate spectral resolution and signal-to-noise ratio. Our work provides important references for the design of future EUV spectrometers for stellar CME detection and the development of observation strategies. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: 20 pages, 9 figures. Accepted for publication in ApJ

Multiwavelength Polarization Observations of Mrk 501

Authors: Xin-Ke Hu, Yu-Wei Yu, Jin Zhang, Xiang-Gao Wang, Kishore C. Patra, Thomas G. Brink, Wei-Kang Zheng, Qi Wang, De-Feng Kong, Liang-Jun Chen, Ji-Wang Zhou, Jia-Xin Cao, Ming-Xuan Lu, Zi-Min Zhou, Yi-Ning Wei, Xin-Bo Huang, Xing-Lin Li, Hao Lou, Ji-Rong Mao, En-Wei Liang, Alexei V. Filippenko

Abstract: Mrk 501 is a prototypical high-synchrotron-peaked blazar (HBL) and serves as one of the primary targets for the {\it Imaging X-ray Polarimetry Explorer} ({\it IXPE}). In this study, we report X-ray polarization measurements of Mrk 501 based on six {\it IXPE} observations. The detection of X-ray polarization at a confidence level exceeding 99\% is achieved in four out of the six observations conduc… ▽ More Mrk 501 is a prototypical high-synchrotron-peaked blazar (HBL) and serves as one of the primary targets for the {\it Imaging X-ray Polarimetry Explorer} ({\it IXPE}). In this study, we report X-ray polarization measurements of Mrk 501 based on six {\it IXPE} observations. The detection of X-ray polarization at a confidence level exceeding 99\% is achieved in four out of the six observations conducted across the entire energy range (2--8 keV) of {\it IXPE}. The maximum polarization degree ($Π_{\rm X}$) is measured to be $15.8\%\pm2.8\%$, accompanied by a polarization angle ($ψ_{\rm X}$) of $98.0°\pm5.1°$ at a confidence level of $5.6 σ$. During the remaining two observations, only an upper limit of $Π_{\rm X}<$12\% could be derived at the 99\% confidence level. No temporal variability in polarization is observed throughout all six {\it IXPE} observations for Mrk 501. A discernible trend of energy-dependent variation in the polarization degree is detected in optical spectropolarimetry; however, no analogous indication is observed in $Π_{\rm X}$. The chromatic behavior of $Π$ and the consistent values of $ψ$ across different frequencies from X-rays to radio waves, along with the agreement between $ψ$ and jet position angle, strongly support the interpretation of the energy-stratified model with shock-accelerated particles in the jet of Mrk 501. Additionally, the possibility of the presence of a global helical magnetic field in the jet of Mrk 501 is discussed. △ Less

Submitted 3 July, 2024; v1 submitted 19 February, 2024; originally announced February 2024.

Comments: 20 pages, 8 figures, 4 tables, accepted for publication in ApJL

Earth's Alfvén wings driven by the April 2023 Coronal Mass Ejection

Authors: Li-Jen Chen, Daniel Gershman, Brandon Burkholder, Yuxi Chen, Menelaos Sarantos, Lan Jian, James Drake, Chuanfei Dong, Harsha Gurram, Jason Shuster, Daniel Graham, Olivier Le Contel, Steven Schwartz, Stephen Fuselier, Hadi Madanian, Craig Pollock, Haoming Liang, Matthew Argall, Richard Denton, Rachel Rice, Jason Beedle, Kevin Genestreti, Akhtar Ardakani, Adam Stanier, Ari Le , et al. (11 additional authors not shown)

Abstract: We report a rare regime of Earth's magnetosphere interaction with sub-Alfvénic solar wind in which the windsock-like magnetosphere transforms into one with Alfvén wings. In the magnetic cloud of a Coronal Mass Ejection (CME) on April 24, 2023, NASA's Magnetospheric Multiscale mission distinguishes the following features: (1) unshocked and accelerated cold CME plasma coming directly against Earth's… ▽ More We report a rare regime of Earth's magnetosphere interaction with sub-Alfvénic solar wind in which the windsock-like magnetosphere transforms into one with Alfvén wings. In the magnetic cloud of a Coronal Mass Ejection (CME) on April 24, 2023, NASA's Magnetospheric Multiscale mission distinguishes the following features: (1) unshocked and accelerated cold CME plasma coming directly against Earth's dayside magnetosphere; (2) dynamical wing filaments representing new channels of magnetic connection between the magnetosphere and foot points of the Sun's erupted flux rope; (3) cold CME ions observed with energized counter-streaming electrons, evidence of CME plasma captured due to reconnection between magnetic-cloud and Alfvén-wing field lines. The reported measurements advance our knowledge of CME interaction with planetary magnetospheres, and open new opportunities to understand how sub-Alfvénic plasma flows impact astrophysical bodies such as Mercury, moons of Jupiter, and exoplanets close to their host stars. △ Less

Submitted 3 May, 2024; v1 submitted 12 February, 2024; originally announced February 2024.

Comments: 14 pages, including 4 figures, Under review in Geophys. Res. Lett

Doppler Tracking Data of Martian Mission Tianwen-I and Upper Limit of Stochastic Gravitational Wave Background

Authors: Xiaoming Bi, Zhongkai Guo, Xiaobo Zou, Yong Huang, Peijia Li, Jianfeng Cao, Lue Chen, Wenlin Tang, Yun Kau Lau

Abstract: Two way ranging data for spacecraft tracking of China's first Martian mission Tianwen-I is analysed. Shortly before the spacecraft entered the Mars parking orbit, the two way coherent microwave link between the spacecraft and the Earth resembles a long arm gravitational wave interferometer, with both the spacecraft and the Earth regarded as in an approximate free falling state. By carefully select… ▽ More Two way ranging data for spacecraft tracking of China's first Martian mission Tianwen-I is analysed. Shortly before the spacecraft entered the Mars parking orbit, the two way coherent microwave link between the spacecraft and the Earth resembles a long arm gravitational wave interferometer, with both the spacecraft and the Earth regarded as in an approximate free falling state. By carefully selecting and analysing data segments of the time series of the two way ranging data during this time span, a parametric statistical model is built for the data segments and an upper limit for the stochastic gravitational waves background (SGWB) is then estimated within the frequency window 0.1Hz to 0.1 mHz. The upper bound improves considerably on those obtained before. In particular, around the deci-Hz band, there is a three orders improvement on the bound obtained previously by the two way ranging data of the Chang e 3 mission. Scientific applications of the upper bound is then considered and a weak upper bound is worked out for axions which is a promising candidate for ultra light dark matter. △ Less

Submitted 8 February, 2024; originally announced February 2024.

Comments: 10 pages, 8 figures

Minute-Cadence Observations of the LAMOST Fields with the TMTS V. Machine Learning Classification of TMTS Catalogues of Periodic Variable Stars

Authors: Fangzhou Guo, Jie Lin, Xiaofeng Wang, Xiaodian Chen, Tanda Li, Liyang Chen, Qiqi Xia, Jun Mo, Gaobo Xi, Jicheng Zhang, Qichun Liu, Xiaojun Jiang, Shengyu Yan, Haowei Peng, Jialian Liu, Wenxiong Li, Weili Lin, Danfeng Xiang, Xiaoran Ma, Yongzhi Cai

Abstract: Periodic variables are always of great scientific interest in astrophysics. Thanks to the rapid advancement of modern large-scale time-domain surveys, the number of reported variable stars has experienced substantial growth for several decades, which significantly deepened our comprehension of stellar structure and binary evolution. The Tsinghua University-Ma Huateng Telescopes for Survey (TMTS) h… ▽ More Periodic variables are always of great scientific interest in astrophysics. Thanks to the rapid advancement of modern large-scale time-domain surveys, the number of reported variable stars has experienced substantial growth for several decades, which significantly deepened our comprehension of stellar structure and binary evolution. The Tsinghua University-Ma Huateng Telescopes for Survey (TMTS) has started to monitor the LAMOST sky areas since 2020, with a cadence of 1 minute. During the period from 2020 to 2022, this survey has resulted in densely sampled light curves for ~ 30,000 variables of the maximum powers in the Lomb-Scargle periodogram above the 5sigma threshold. In this paper, we classified 11,638 variable stars into 6 main types using XGBoost and Random Forest classifiers with accuracies of 98.83% and 98.73%, respectively. Among them, 5301 (45.55%) variables are newly discovered, primarily consisting of Delta Scuti stars, demonstrating the capability of TMTS in searching for short-period variables. We cross-matched the catalogue with Gaia's second Data Release (DR2) and LAMOST's seventh Data Release (DR7) to obtain important physical parameters of the variables. We identified 5504 Delta Scuti stars (including 4876 typical Delta Scuti stars and 628 high-amplitude Delta Scuti stars), 5899 eclipsing binaries (including EA-, EB- and EW-type) and 226 candidates of RS Canum Venaticorum. Leveraging the metal abundance data provided by LAMOST and the Galactic latitude, we discovered 8 candidates of SX Phe stars within the class of "Delta Scuti stars". Moreover, with the help of Gaia color-magnitude diagram, we identified 9 ZZ ceti stars. △ Less

Submitted 4 February, 2024; originally announced February 2024.

Comments: 20 pages, 14 figures, accepted to MNRAS

Evolution of magnetic field of the Quasar 1604+159 at pc scale

Authors: Xu-Zhi Hu, Xiaoyu Hong, Wei Zhao, Liang Chen, Wei-Yang Wang, Linhui Wu

Abstract: We have analyzed the total intensity, spectral index, linear polarization, and RM distributions at pc scale for the quasar 1604+159. The source was observed in 2002 and 2020 with the VLBA. Combining the MOJAVE results, we studied the evolution of the magnetic field. We detected a core-jet structure. The jet extends to a distance of ~25 mas. The jet shape varies slightly with time. We divided the s… ▽ More We have analyzed the total intensity, spectral index, linear polarization, and RM distributions at pc scale for the quasar 1604+159. The source was observed in 2002 and 2020 with the VLBA. Combining the MOJAVE results, we studied the evolution of the magnetic field. We detected a core-jet structure. The jet extends to a distance of ~25 mas. The jet shape varies slightly with time. We divided the source structure into the central region and the jet region. In the jet region, we find the polarized emission varies with time. The flatter spectral index values and EVPA direction indicate the possible existence of shocks, contributing to the variation. In the central region, the derived core shift index k_r values indicate that the core in 2002 is close to the equipartition case while deviating from it in 2020. The measured magnetic field strength in 2020 is two orders of magnitude lower than that in 2002. We detected transverse RM gradients, evidence of a helical magnetic field, in the core. At 15 GHz, in the place close to the jet base, the polarization direction changes significantly with time from perpendicular to parallel to the jet direction. The evolution of RM and magnetic field structure are potential reasons for the observed polarization change. The core |RM| in 2020 increases with frequency following a power law with index a = 2.7, suggesting a fast electron density fall-off in the medium with distance from the jet base. △ Less

Submitted 1 February, 2024; v1 submitted 30 January, 2024; originally announced January 2024.

Comments: 24 pages, 14 figures, accepted for publication in ApJ

Resonant amplification of curvature perturbations in inflation model with periodical derivative coupling

Authors: Li-Yang Chen, Hongwei Yu, Puxun Wu

Abstract: In this paper, we introduce a weak, transient and periodical derivative coupling between the inflaton field and gravity, and find that the square of the sound speed of the curvature perturbations becomes a periodic function, which results in that the equation of the curvature perturbations can be transformed into the form of the Mathieu equation in the sub-horizon limit. Thus, the parametric reson… ▽ More In this paper, we introduce a weak, transient and periodical derivative coupling between the inflaton field and gravity, and find that the square of the sound speed of the curvature perturbations becomes a periodic function, which results in that the equation of the curvature perturbations can be transformed into the form of the Mathieu equation in the sub-horizon limit. Thus, the parametric resonance will amplify the curvature perturbations so as to generate a formation of abundant primordial black holes (PBHs). We show that the generated PBHs can make up most of dark matter. Associated with the generation of PBHs, the large scalar perturbations will give rise to the scalar induced gravitational waves which may be detected by future gravitational wave projects. △ Less

Submitted 27 January, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Comments: 19 pages, 4 figures, refs added

Journal ref: Phys. Lett. B 849 (2024) 138457

Mid-infrared evidence for iron-rich dust in the multi-ringed inner disk of HD 144432

Authors: J. Varga, L. B. F. M. Waters, M. Hogerheijde, R. van Boekel, A. Matter, B. Lopez, K. Perraut, L. Chen, D. Nadella, S. Wolf, C. Dominik, Á. Kóspál, P. Ábrahám, J. -C. Augereau, P. Boley, G. Bourdarot, A. Caratti o Garatti, F. Cruz-Sáenz de Miera, W. C. Danchi, V. Gámez Rosas, Th. Henning, K. -H. Hofmann, M. Houllé, J. W. Isbell, W. Jaffe , et al. (18 additional authors not shown)

Abstract: Context. Rocky planets form by the concentration of solid particles in the inner few au regions of planet-forming disks. Their chemical composition reflects the materials in the disk available in the solid phase at the time the planets were forming. Aims. We aim to constrain the structure and dust composition of the inner disk of the young star HD 144432, using an extensive set of infrared interfe… ▽ More Context. Rocky planets form by the concentration of solid particles in the inner few au regions of planet-forming disks. Their chemical composition reflects the materials in the disk available in the solid phase at the time the planets were forming. Aims. We aim to constrain the structure and dust composition of the inner disk of the young star HD 144432, using an extensive set of infrared interferometric data taken by the Very Large Telescope Interferometer (VLTI), combining PIONIER, GRAVITY, and MATISSE observations. Methods. We introduced a new physical disk model, TGMdust, to image the interferometric data, and to fit the disk structure and dust composition. We also performed equilibrium condensation calculations with GGchem. Results. Our best-fit model has three disk zones with ring-like structures at 0.15, 1.3, and 4.1 au. Assuming that the dark regions in the disk at ~0.9 au and at ~3 au are gaps opened by planets, we estimate the masses of the putative gap-opening planets to be around a Jupiter mass. We find evidence for an optically thin emission ($τ<0.4$) from the inner two disk zones ($r<4$ au) at $λ>3\ μ$m. Our silicate compositional fits confirm radial mineralogy gradients. To identify the dust component responsible for the infrared continuum emission, we explore two cases for the dust composition, one with a silicate+iron mixture and the other with a silicate+carbon one. We find that the iron-rich model provides a better fit to the spectral energy distribution. Conclusions. We propose that in the warm inner regions ($r<5$ au) of typical planet-forming disks, most if not all carbon is in the gas phase, while iron and iron sulfide grains are major constituents of the solid mixture along with forsterite and enstatite. Our analysis demonstrates the need for detailed studies of the dust in inner disks with new mid-infrared instruments such as MATISSE and JWST/MIRI. △ Less

Submitted 7 January, 2024; originally announced January 2024.

Comments: 29 pages, 24 figures

Journal ref: A&A, 681, A47 (2024)

Extended Skyrme effective interactions for transport model and neutron stars

Authors: Si-Pei Wang, Rui Wang, Jun-Ting Ye, Lie-Wen Chen

Abstract: It is important to develop a unified theoretical framework to describe the nuclear experiments and astrophysical observations based on the same effective nuclear interactions. Based on the so-called Skyrme pseudopotential up to next-to-next-to-next-to-leading order, we construct a series of extended Skyrme interactions by modifying the density-dependent term and fitting the empirical nucleon optic… ▽ More It is important to develop a unified theoretical framework to describe the nuclear experiments and astrophysical observations based on the same effective nuclear interactions. Based on the so-called Skyrme pseudopotential up to next-to-next-to-next-to-leading order, we construct a series of extended Skyrme interactions by modifying the density-dependent term and fitting the empirical nucleon optical potential up to above $1$ GeV, the empirical properties of isospin symmetric nuclear matter, the microscopic calculations of pure neutron matter and the properties of neutron stars from astrophysical observations. The modification of the density-dependent term in the extended Skyrme interactions follows the idea of Fermi momentum expansion and this leads to a highly flexible density behavior of the symmetry energy. In particular, the values of the density slope parameter $L$ of the symmetry energy for the new extended Skyrme interactions range from $L = -5$ MeV to $L = 125$ MeV by construction, to cover the large uncertainty of the density dependence of the symmetry energy. Furthermore, in order to consider the effects of isoscalar and isovector nucleon effective masses, we adjust the momentum dependency of the single-nucleon optical potential and the symmetry potential of these new extended Skyrme interactions and construct a parameter set family, by which we systematically study the impacts of the symmetry energy and the nucleon effective masses on the properties of nuclear matter and neutron stars. The new extended Skyrme interactions constructed in the present work will be useful to determine the equation of state of isospin asymmetric nuclear matter, especially the symmetry energy, as well as the nucleon effective masses and their isospin splitting, in transport model simulations for heavy-ion collisions, nuclear structure calculations and neutron star studies. △ Less

Submitted 22 April, 2024; v1 submitted 28 December, 2023; originally announced December 2023.

Comments: 38 pages, 16 figures, 8 tables. Discussions added and typos fixed. Accepted version to appear in PRC

A seven-Earth-radius helium-burning star inside a 20.5-min detached binary

Authors: Jie Lin, Chengyuan Wu, Heran Xiong, Xiaofeng Wang, Peter Nemeth, Zhanwen Han, Jiangdan Li, Nancy Elias-Rosa, Irene Salmaso, Alexei V. Filippenko, Thomas G. Brink, Yi Yang, Xuefei Chen, Shengyu Yan, Jujia Zhang, Sufen Guo, Yongzhi Cai, Jun Mo, Gaobo Xi, Jialian Liu, Jincheng Guo, Qiqi Xia, Danfeng Xiang, Gaici Li, Zhenwei Li , et al. (6 additional authors not shown)

Abstract: Binary evolution theory predicts that the second common envelope (CE) ejection can produce low-mass (0.32-0.36 Msun) subdwarf B (sdB) stars inside ultrashort-orbital-period binary systems, as their helium cores are ignited under nondegenerate conditions. With the orbital decay driven by gravitational-wave (GW) radiation, the minimum orbital periods of detached sdB binaries could be as short as ~20… ▽ More Binary evolution theory predicts that the second common envelope (CE) ejection can produce low-mass (0.32-0.36 Msun) subdwarf B (sdB) stars inside ultrashort-orbital-period binary systems, as their helium cores are ignited under nondegenerate conditions. With the orbital decay driven by gravitational-wave (GW) radiation, the minimum orbital periods of detached sdB binaries could be as short as ~20 minutes. However, only four sdB binaries with orbital periods below an hour have been reported so far, while none of them has an orbital period approaching the above theoretical limit. Here we report the discovery of a 20.5-minute-orbital-period ellipsoidal binary, TMTS J052610.43+593445.1, in which the visible star is being tidally deformed by an invisible carbon-oxygen white dwarf (WD) companion. The visible component is inferred to be an sdB star with a mass of ~0.33 Msun, approaching that of helium-ignition limit, although a He-core WD cannot be completely ruled out. In particular, the radius of this low-mass sdB star is only 0.066 Rsun, about seven Earth radii, possibly representing the most compact nondegenerate star ever known. Such a system provides a key clue to map the binary evolution scheme from the second CE ejection to the formation of AM CVn stars having a helium-star donor, and it will also serve as a crucial verification binary of space-borne GW detectors in the future. △ Less

Submitted 10 February, 2024; v1 submitted 21 December, 2023; originally announced December 2023.

Comments: 24 pages, 11 figures, 1 table, published on Nature Astronomy, URL: https://www.nature.com/articles/s41550-023-02188-2

Weak Solar Radio Bursts from the Solar Wind Acceleration Region Observed by Parker Solar Probe and Its Probable Emission Mechanism

Authors: Ling Chen, Bing Ma, Dejin Wu, Xiaowei Zhou, Marc Pulupa, PeiJin Zhang, Pietro Zucca, Stuart D. Bale, Justin C. Kasper, SuPing Duan

Abstract: The Parker Solar Probe (PSP) provides us the unprecedentedly close approach observation to the Sun, and hence the possibility of directly understanding the "elementary process" which occurs in the kinetic scale of particles collective interactioin in solar coronal plasmas. We reported a kind of weak solar radio bursts (SRBs), which are detected by PSP when it passed a low-density magnetic channel… ▽ More The Parker Solar Probe (PSP) provides us the unprecedentedly close approach observation to the Sun, and hence the possibility of directly understanding the "elementary process" which occurs in the kinetic scale of particles collective interactioin in solar coronal plasmas. We reported a kind of weak solar radio bursts (SRBs), which are detected by PSP when it passed a low-density magnetic channel during its second encounter phase. These weak SRBs have low starting frequecny $\sim 20$ MHz and narrow frequency range from a few tens MHz to a few hundres kHz. Their dynamic spectra display a strongly evolving feature of the intermediate relative drift rate decreasing rapidly from above 0.01/s to below 0.01/s. Analyses based on common empirical models of solar coronal plasmas indicate that these weak SRBs originate from the heliocentric distance $\sim 1.1-6.1~R_S$ (the solar radius), a typical solar wind acceleration region with a low-$β$ plasma, and indicate that their soruces have a typic motion velociy $\sim v_A$ (Alfvén velocity) obviously lower than that of fast electrons required by effectively exciting SRBs. We propose that solitary kinetic Alfvén waves with kinetic scales can be responsible for the generation of these small-scalevweak SRBs, called solitary wave radiation (SWR). △ Less

Submitted 29 November, 2023; originally announced November 2023.

ALMA High-resolution Spectral Survey of Thioformaldehyde (H2CS) Towards Massive Protoclusters

Authors: Li Chen, Sheng-Li Qin, Tie Liu, Hong-Li Liu, Sheng-Yuan Liu, Meizhu Liu, Hongqiong Shi, Chuanshou Li, Mengyao Tang, Tianwei Zhang, Ken'ichi Tatematsu, Xiaohu Li, Fengwei Xu, Yuefang Wu, Dongting Yang

Abstract: Investigating the temperature and density structures of gas in massive protoclusters is crucial for understanding the chemical properties therein. In this study, we present observations of the continuum and thioformaldehyde (H2CS) lines at 345 GHz of 11 massive protoclusters using the Atacama Large Millimeter/submillimeter Array (ALMA) telescope. High spatial resolution and sensitivity observation… ▽ More Investigating the temperature and density structures of gas in massive protoclusters is crucial for understanding the chemical properties therein. In this study, we present observations of the continuum and thioformaldehyde (H2CS) lines at 345 GHz of 11 massive protoclusters using the Atacama Large Millimeter/submillimeter Array (ALMA) telescope. High spatial resolution and sensitivity observations have detected 145 continuum cores from the 11 sources. H2CS line transitions are observed in 72 out of 145 cores, including line-rich cores, warm cores and cold cores. The H2 column densities of the 72 cores are estimated from the continuum emission which are larger than the density threshold value for star formation, suggesting that H2CS can be widely distributed in star-forming cores with different physical environments. Rotation temperature and column density of H2CS are derived by use of the XCLASS software. The results show the H2CS abundances increase as temperature rises and higher gas temperatures are usually associated with higher H2CS column densities. The abundances of H2CS are positively correlated with its column density, suggesting that the H2CS abundances are enhanced from cold cores, warm cores to line-rich cores in star forming regions. △ Less

Submitted 29 November, 2023; originally announced November 2023.

Comments: 23 pages, 7 figures, 6 tables, accepted by ApJ

A Review of Axion Lasing in Astrophysics

Authors: Liang Chen, Thomas W. Kephart

Abstract: Axions can be stimulated to decay to photons by ambient photons of the right frequency or by photons from decay of neighboring axions. If the axion density is high enough the photon intensity can be amplified, which is a type of lasing or an axion maser. Here we review the astrophysical situations where axion lasing can appear and possibly be detected. Axions can be stimulated to decay to photons by ambient photons of the right frequency or by photons from decay of neighboring axions. If the axion density is high enough the photon intensity can be amplified, which is a type of lasing or an axion maser. Here we review the astrophysical situations where axion lasing can appear and possibly be detected. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: 19 pages

The first Ka-band (26.1-35 GHz) blind line survey towards Orion KL

Authors: Xunchuan Liu, Tie Liu, Zhiqiang Shen, Sheng-Li Qin, Qiuyi Luo, Yan Gong, Yu Cheng, Christian Henkel, Qilao Gu, Fengyao Zhu, Tianwei Zhang, Rongbing Zhao, Yajun Wu, Bin Li, Juan Li, Zhang Zhao, Jinqing Wang, Weiye Zhong, Qinghui Liu, Bo Xia, Li Fu, Zhen Yan, Chao Zhang, Lingling Wang, Qian Ye , et al. (9 additional authors not shown)

Abstract: We conducted a Ka-band (26.1--35 GHz) line survey towards Orion KL using the TianMa 65-m Radio Telescope (TMRT). It is the first blind line survey in the Ka band, and achieves a sensitivity of mK level (1--3 mK at a spectral resolution of $\sim$1 km s$^{-1}$). In total, 592 Gaussian features are extracted. Among them, 257 radio recombination lines (RRLs) are identified. The maximum $Δn$ of RRLs of… ▽ More We conducted a Ka-band (26.1--35 GHz) line survey towards Orion KL using the TianMa 65-m Radio Telescope (TMRT). It is the first blind line survey in the Ka band, and achieves a sensitivity of mK level (1--3 mK at a spectral resolution of $\sim$1 km s$^{-1}$). In total, 592 Gaussian features are extracted. Among them, 257 radio recombination lines (RRLs) are identified. The maximum $Δn$ of RRLs of H, He and C are 20, 15, and 5, respectively. Through stacking, we have detected the $β$ lines of ion RRLs (RRLs of C$^+$ with possible contribution of other ions like O$^+$) for the first time, and tentative signal of the $γ$ lines of ion RRLs can also be seen on the stacked spectrum. Besides, 318 other line features were assigned to 37 molecular species, and ten of these species were not detected in the Q-band survey of TMRT. The vibrationally excited states of nine species were also detected. Emission of most species can be modeled under LTE. A number of transitions of E-CH3OH ($J_2-J_1$) display maser effects, which are confirmed by our modeling, and besides the bumping peak at $J\sim 6$ there is another peak at $J\sim 13$. Methylcyanoacetylene (CH$_3$C$_3$N) is detected in Orion KL for the first time. This work emphasizes that the Ka band, which was long-ignored for spectral line surveys, is very useful for surveying RRLs and molecular lines simultaneously. △ Less

Submitted 20 November, 2023; originally announced November 2023.

Comments: accepted by ApJS

The disk of the eruptive protostar V900 Mon; a MATISSE/VLTI and MUSE/VLT perspective

Authors: F. Lykou, P. Ábrahám, F. Cruz-Sáenz de Miera, J. Varga, Á. Kóspál, J. Bouwman, L. Chen, S. Kraus, M. L. Sitko, R. W. Russell, M. Pikhartova

Abstract: In this work, we study the silicate dust content in the disk of one of the youngest eruptive stars, V900 Mon, at the highest angular resolution probing down to the inner 10 au of said disk, and study the historical evolution of the system traced in part by a newly discovered emission clump. We performed high-angular resolution mid-infrared interferometric observations of V900 Mon with MATISSE/VLTI… ▽ More In this work, we study the silicate dust content in the disk of one of the youngest eruptive stars, V900 Mon, at the highest angular resolution probing down to the inner 10 au of said disk, and study the historical evolution of the system traced in part by a newly discovered emission clump. We performed high-angular resolution mid-infrared interferometric observations of V900 Mon with MATISSE/VLTI with a spatial coverage ranging from 38-m to 130-m baselines, and compared them to archival MIDI/VLTI data. We also mined and re-analyzed archival optical and infrared photometry of the star to study its long-term evolution since its eruption in the 1990s. We complemented our findings with integral field spectroscopy data from MUSE/VLT. The MATISSE/VLTI data suggest a radial variation of the silicate feature in the dusty disk, whereby at large spatial scales ($\geq10$ au) the protostellar disk's emission is dominated by large-sized ($\geq1\,μm$) silicate grains, while at smaller spatial scales and closer to the star ($\leq5$ au), silicate emission is absent suggesting self-shielding. We propose that the self-shielding may be the result of small dust grains at the base of the collimated CO outflow previously detected by ALMA. A newly discovered knot in the MUSE/VLT data, located at a projected distance approximately 27,000 au from the star, is co-aligned with the molecular gas outflow at a P.A. of $250^o$ ($\pm5^o$) consistent with the position angle and inclination of the disk. The knot is seen in emission in H$α$, [N II], and the [S II] doublet and its kinematic age is about 5150 years. This ejected material could originate from a previous eruption. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: 18 pages, 16 figures, accepted for publication by Astronomy & Astrophysics

Effects of stimulated emission and superradiant growth of non-spherical axion cluster

Authors: Liang Chen, Da Huang, Chao-Qiang Geng

Abstract: We explore the stimulated emission of photons in non-spherical axion clusters with or without the axion source from the superradiance of a rotating black hole (BH). In particular, we focus on the cluster with the initial axion distribution in the $(l,m)=(1,1)$ mode which mimics the shape of an axion cloud induced by the BH superradiance. After establishing the hierarchy of Boltzmann equations gove… ▽ More We explore the stimulated emission of photons in non-spherical axion clusters with or without the axion source from the superradiance of a rotating black hole (BH). In particular, we focus on the cluster with the initial axion distribution in the $(l,m)=(1,1)$ mode which mimics the shape of an axion cloud induced by the BH superradiance. After establishing the hierarchy of Boltzmann equations governing a general non-spherical axion-photon system, we examine the evolution of photon and axion distributions in the cluster and possible stimulated emission signals. In the case without the axion source, the resultant signal would be a single photon pulse. As for the system with the BH superradiance as the axion source, multiple pulses are predicted. We also show that, for the latter case, the combined effects of stimulated emissions and the axion production from the BH superradiance could reach a balance where the axion cluster becomes uniformly and spherically distributed. Due to the energy and temporal characteristics of the obtained pulses, we demonstrate that the stimulated emissions from the axion cluster with axions sourced by the BH superradiance provide a candidate explanation to the observed fast radio bursts. △ Less

Submitted 30 May, 2024; v1 submitted 3 November, 2023; originally announced November 2023.

Comments: 27 pages, 10 figures, 1 table; to appear in JCAP

A Systematic Observational Study on Galactic Interstellar Ratio 18O/17O. II. C18O and C17O J=2-1 Data Analysis

Authors: Y. P. Zou, J. S. Zhang, C. Henkel, D. Romano, W. Liu, Y. H. Zheng, Y. T. Yan, J. L. Chen, Y. X. Wang, J. Y. Zhao

Abstract: To investigate the relative amount of ejecta from high-mass versus intermediate-mass stars and to trace the chemical evolution of the Galaxy, we have performed with the IRAM 30m and the SMT 10m telescopes a systematic study of Galactic interstellar 18O/17O ratios toward a sample of 421 molecular clouds, covering a galactocentric distance range of 1-22 kpc. The results presented in this paper are b… ▽ More To investigate the relative amount of ejecta from high-mass versus intermediate-mass stars and to trace the chemical evolution of the Galaxy, we have performed with the IRAM 30m and the SMT 10m telescopes a systematic study of Galactic interstellar 18O/17O ratios toward a sample of 421 molecular clouds, covering a galactocentric distance range of 1-22 kpc. The results presented in this paper are based on the J=2-1 transition and encompass 364 sources showing both C18O and C17O detections. The previously suggested 18O/17O gradient is confirmed. For the 41 sources detected with both facilities, good agreement is obtained. A correlation of 18O/17O ratios with heliocentric distance is not found, indicating that beam dilution and linear beam sizes are not relevant. For the subsample of IRAM 30 m high-mass star-forming regions with accurate parallax distances, an unweighted fit gives 18O/17O = (0.12+-0.02)R_GC+(2.38+-0.13) with a correlation coefficient of R = 0.67. While the slope is consistent with our J=1-0 measurement, ratios are systematically lower. This should be caused by larger optical depths of C18O 2-1 lines, w.r.t the corresponding 1-0 transitions, which is supported by RADEX calculations and the fact that C18O/C17O is positively correlated with 13CO/C18O. After considering optical depth effects with C18O J=2-1 reaching typically an optical depth of 0.5, corrected 18O/17O ratios from the J=1-0 and J=2-1 lines become consistent. A good numerical fit to the data is provided by the MWG-12 model, including both rotating stars and novae. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: 17 pages, 11 figures, published in ApJS

Astronomical Knowledge Entity Extraction in Astrophysics Journal Articles via Large Language Models

Authors: Wujun Shao, Pengli Ji, Dongwei Fan, Yaohua Hu, Xiaoran Yan, Chenzhou Cui, Linying Mi, Lang Chen, Rui Zhang

Abstract: Astronomical knowledge entities, such as celestial object identifiers, are crucial for literature retrieval and knowledge graph construction, and other research and applications in the field of astronomy. Traditional methods of extracting knowledge entities from texts face challenges like high manual effort, poor generalization, and costly maintenance. Consequently, there is a pressing need for im… ▽ More Astronomical knowledge entities, such as celestial object identifiers, are crucial for literature retrieval and knowledge graph construction, and other research and applications in the field of astronomy. Traditional methods of extracting knowledge entities from texts face challenges like high manual effort, poor generalization, and costly maintenance. Consequently, there is a pressing need for improved methods to efficiently extract them. This study explores the potential of pre-trained Large Language Models (LLMs) to perform astronomical knowledge entity extraction (KEE) task from astrophysical journal articles using prompts. We propose a prompting strategy called Prompt-KEE, which includes five prompt elements, and design eight combination prompts based on them. Celestial object identifier and telescope name, two most typical astronomical knowledge entities, are selected to be experimental object. And we introduce four currently representative LLMs, namely Llama-2-70B, GPT-3.5, GPT-4, and Claude 2. To accommodate their token limitations, we construct two datasets: the full texts and paragraph collections of 30 articles. Leveraging the eight prompts, we test on full texts with GPT-4 and Claude 2, on paragraph collections with all LLMs. The experimental results demonstrated that pre-trained LLMs have the significant potential to perform KEE tasks in astrophysics journal articles, but there are differences in their performance. Furthermore, we analyze some important factors that influence the performance of LLMs in entity extraction and provide insights for future KEE tasks in astrophysical articles using LLMs. △ Less

Submitted 17 January, 2024; v1 submitted 27 October, 2023; originally announced October 2023.

Baryon density dependence of viscosities of the quark-gluon plasma at hadronization

Authors: Zhidong Yang, Yifeng Sun, Lie-Wen Chen

Abstract: The $φ$ meson and $Ω$ baryon provide unique probes of the properties of the quark-gluon plasma (QGP) at hadronization in relativistic heavy-ion collisions. Using the quark recombination model with the quark phase-space information parameterized in a viscous blastwave, we perform Bayesian inference of the shear and bulk viscosities of the QGP at hadronization with a temperature of $T\sim 160$ MeV b… ▽ More The $φ$ meson and $Ω$ baryon provide unique probes of the properties of the quark-gluon plasma (QGP) at hadronization in relativistic heavy-ion collisions. Using the quark recombination model with the quark phase-space information parameterized in a viscous blastwave, we perform Bayesian inference of the shear and bulk viscosities of the QGP at hadronization with a temperature of $T\sim 160$ MeV by analyzing the $φ$ and $Ω$ data in Au+Au collisions at $\sqrt{s_{\rm NN}}=$ 19.6-200 GeV and Pb+Pb collisions at $\sqrt{s_{\rm NN}}=$ 2.76 TeV, corresponding to a baryon chemical potential variation from $μ_B\approx 0$ (at $\sqrt{s_{\rm NN}}= 2.76$ TeV) to $200$ MeV (at $\sqrt{s_{\rm NN}}= 19.6$ GeV). We find that the shear viscosity to enthalpy ratio $ηT/(ε+P)$ of the QGP at hadronization decreases as $μ_B$ increases, with $ηT/(ε+P)\approx 0.18$ at $μ_B=0$ and $ηT/(ε+P)\approx 0.08$ at $μ_B=200$ MeV, while the corresponding specific bulk viscosity is essentially constant with $ζT/(ε+ P)=0.02\sim 0.04$ for $μ_B<200$ MeV. Our results suggest that the QGP at hadronization ($T\sim 160$ MeV) with finite baryon density is more close to perfect fluid than that with zero baryon density. △ Less

Submitted 17 April, 2024; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 9 pages, 3 figures, 2 tables. Some details and discussions added. Accepted version to appear in PRC

Age Determination of LAMOST Red Giant Branch stars based on the Gradient Boosting Decision Tree method

Authors: Hai-Feng Wang, Giovanni Carraro, Xin Li, Qi-Da Li, Lorenzo Spina, Li Chen, Guan-Yu Wang, Li-Cai Deng

Abstract: In this study we estimate the stellar ages of LAMOST DR8 Red Giant Branch (RGB) stars based on the Gradient Boosting Decision Tree algorithm (GBDT). We used 2,643 RGB stars extracted from the APOKASC-2 astero-seismological catalog as training data-set. After selecting the parameterses ([$α$/Fe], [C/Fe], T$_{eff}$, [N/Fe], [C/H], log g) highly correlated with age using GBDT, we apply the same GBDT… ▽ More In this study we estimate the stellar ages of LAMOST DR8 Red Giant Branch (RGB) stars based on the Gradient Boosting Decision Tree algorithm (GBDT). We used 2,643 RGB stars extracted from the APOKASC-2 astero-seismological catalog as training data-set. After selecting the parameterses ([$α$/Fe], [C/Fe], T$_{eff}$, [N/Fe], [C/H], log g) highly correlated with age using GBDT, we apply the same GBDT method to the new catalog of more than 590,000 stars classified as RGB stars. The test data-set shows that the median relative error is around 11.6$\%$ for the method. We also compare the predicted ages of RGB stars with other studies (e.g., based on APOGEE), and find systematic differences. The final uncertainty is about 15 to 30$\%$ compared to open clusters' ages. Then we present the spatial distribution of the RGB sample having an age determination, which could recreate the expected result, and discuss systematic biases. All these diagnostics show that one can apply the GBDT method to other stellar samples to estimate atmospheric parameters and age. △ Less

Submitted 7 April, 2024; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 17 pages, 17 figures, Accepted for publication in The Astrophysical Journal

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

A High-Mass Young Star-forming Core Escaping from Its Parental Filament

Authors: Zhiyuan Ren, Xi Chen, Tie Liu, Emma Mannfors, Leonardo Bronfman, Fengwei Xu, Siyi Feng, Hongli Liu, Fanyi Meng, Amelia. M. Stutz, Shanghuo Li, Chang Won Lee, Ke Wang, Jianwen Zhou, Di Li, Chen Wang, Chakali Eswaraiah, Anandmayee Tej, Long-Fei Chen, Hui Shi

Abstract: We studied the unique kinematic properties in massive filament G352.63-1.07 at $10^3$-AU spatial scale with the dense molecular tracers observed with the Atacama Large Millimeter/submillimeter Array (ALMA). We find the central massive core M1 (12 $M_\odot$) being separated from the surrounding filament with a velocity difference of $v- {v}_{sys}=-2$ km/s and a transverse separation within 3 arcsec… ▽ More We studied the unique kinematic properties in massive filament G352.63-1.07 at $10^3$-AU spatial scale with the dense molecular tracers observed with the Atacama Large Millimeter/submillimeter Array (ALMA). We find the central massive core M1 (12 $M_\odot$) being separated from the surrounding filament with a velocity difference of $v- {v}_{sys}=-2$ km/s and a transverse separation within 3 arcsec. Meanwhile, as shown in multiple dense-gas tracers, M1 has a spatial extension closely aligned with the main filament and is connected to the filament towards its both ends. M1 thus represents a very beginning state for a massive young star-forming core escaping from the parental filament, within a time scale of $\sim 4000$ years. Based on its kinetic energy ($3.5\times10^{44}$ erg), the core escape is unlikely solely due to the original filament motion or magnetic field, but requires more energetic events such as a rapid intense anisotropic collapse. The released energy also seems to noticeably increase the environmental turbulence. This may help the filament to become stabilized again. △ Less

Submitted 12 October, 2023; originally announced October 2023.

Comments: 15 pages, 6 figures, accepted for publication in the Astrophysical Journal

A low-mass line-rich core found in Massive Star-forming Region IRAS 16351-4722

Authors: Meizhu Liu, Sheng-Li Qin, Tie Liu, Mengyao Tang, Sheng-Yuan Liu, Li Chen, ChuanShou Li, HongQiong Shi, Xiaohu Li, Tianwei Zhang, Ken'ichi Tatematsu, Fengwei Xu, Yuefang Wu

Abstract: We present ALMA sub-arcsecond-resolution observations of both continuum and molecular lines at 345 GHz towards the massive star-forming region IRAS 16351-4722 (hereafter I16351). A total of 12 dust cores were detected based on high spatial resolution observations of the continuum. Among them, a high-mass core (11.6 Msun) and a low-mass core (1.7 Msun) show abundant molecular line emissions. 164 mo… ▽ More We present ALMA sub-arcsecond-resolution observations of both continuum and molecular lines at 345 GHz towards the massive star-forming region IRAS 16351-4722 (hereafter I16351). A total of 12 dust cores were detected based on high spatial resolution observations of the continuum. Among them, a high-mass core (11.6 Msun) and a low-mass core (1.7 Msun) show abundant molecular line emissions. 164 molecular transitions from 29 species and 104 molecular transitions from 25 species are identified in the high-mass and low-mass cores, respectively. Complex organic molecules (COMs) such as CH3OH, CH3OCHO, CH3OCH3, C2H5OH, and C2H5CN are detected in the two cores. Under the assumption of local thermodynamic equilibrium (LTE), rotational temperatures and column densities of the COMs are derived with the XCLASS software. The maximum rotation temperature values in the low-mass core and the high-mass core were found to be approximately 130 K and 198 K, respectively. Additionally, the line widths in the high-mass core are larger than those in the low-mass one. Abundant complex organic molecular line transitions, high gas temperatures, and smaller line widths indicate the presence of a low-mass line-rich core in the massive star formation region for the first time, while the high-mass line-rich core shows hot core property. When comparing the molecular abundances of CH3OH, CH3OCHO, CH3OCH3 and C2H5OH of the two cores with other hot cores and hot corinos reported in the literature, we further confirm that both a hot core and a low-mass line-rich core are simultaneously detected in I16351. △ Less

Submitted 6 October, 2023; originally announced October 2023.

Comments: 22 pages, 5 figures, 5 tables, 70 references, accepted by ApJ

The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE): Core Growth, Cluster Contraction, and Primordial Mass Segregation

Authors: Fengwei Xu, Ke Wang, Tie Liu, Mengyao Tang, Neal J. Evans II, Aina Palau, Kaho Morii, Jinhua He, Patricio Sanhueza, Hong-Li Liu, Amelia Stutz, Qizhou Zhang, Xi Chen, Pak Shing Li, Gilberto C. Gómez, Enrique Vázquez-Semadeni, Shanghuo Li, Xiaofeng Mai, Xing Lu, Meizhu Liu, Li Chen, Chuanshou Li, Hongqiong Shi, Zhiyuan Ren, Di Li , et al. (18 additional authors not shown)

Abstract: The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (Mclump>1000 Msun), luminous (Lbol>10,000 Lsun), and blue-profile (infall signature) clumps by ALMA with resolution of 2200-5500 au a… ▽ More The ALMA Survey of Star Formation and Evolution in Massive Protoclusters with Blue Profiles (ASSEMBLE) aims to investigate the process of mass assembly and its connection to high-mass star formation theories in protoclusters in a dynamic view. We observed 11 massive (Mclump>1000 Msun), luminous (Lbol>10,000 Lsun), and blue-profile (infall signature) clumps by ALMA with resolution of 2200-5500 au at 350 GHz (870 um) in continuum and line emission. 248 dense cores were identified, including 106 cores showing protostellar signatures and 142 prestellar core candidates. Compared to early-stage infrared dark clouds (IRDCs) by ASHES, the core mass and surface density within the ASSEMBLE clumps exhibited significant increment, suggesting concurrent core accretion during the evolution of the clumps. The maximum mass of prestellar cores was found to be 2 times larger than that in IRDCs, indicating evolved protoclusters have the potential to harbor massive prestellar cores. The mass relation between clumps and their most massive core (MMCs) is observed in ASSEMBLE but not in IRDCs, which is suggested to be regulated by multiscale mass accretion. The mass correlation between the core clusters and their MMCs has a steeper slope compared to that observed in stellar clusters, which can be due to fragmentation of the MMC and stellar multiplicity. We observe a decrease in core separation and an increase in central concentration as protoclusters evolve. We confirm primordial mass segregation in the ASSEMBLE protoclusters, possibly resulting from gravitational concentration and/or gas accretion. △ Less

Submitted 26 September, 2023; originally announced September 2023.

Comments: 37 pages, 13 figures, 5 tables; accepted for publication in ApJS

Dwarf galaxies with the highest concentration are not thicker than ordinary dwarf galaxies

Authors: Lijun Chen, Hong-Xin Zhang, Zesen Lin, Guangwen Chen, Bojun Tao, Zhixiong Liang, Zheyu Lin, Xu Kong

Abstract: The formation mechanism of high-concentration dwarf galaxies is still a mystery. We perform a comparative study of the intrinsic shape of nearby low-mass galaxies with different stellar concentration. The intrinsic shape is parameterized by the intermediate-to-major axis ratios B/A and the minor-to-major axis ratios C/A of triaxial ellipsoidal models. Our galaxies ($10^{7.5} M_\odot$ < $M_\star$ <… ▽ More The formation mechanism of high-concentration dwarf galaxies is still a mystery. We perform a comparative study of the intrinsic shape of nearby low-mass galaxies with different stellar concentration. The intrinsic shape is parameterized by the intermediate-to-major axis ratios B/A and the minor-to-major axis ratios C/A of triaxial ellipsoidal models. Our galaxies ($10^{7.5} M_\odot$ < $M_\star$ < $10^{10.0} M_\odot$) are selected to have spectroscopic redshift from SDSS or GAMA, and have broadband optical images from the HSC-SSP Wide layer survey. The deep HSC-SSP images allow to measure the apparent axis ratios $q$ at galactic radii beyond the central star-forming area of our galaxies. We infer the intrinsic axis ratios based on the $q$ distributions. We find that 1) our galaxies have typical intrinsic shape similarly close to be oblate ($μ_{B/A}$ $\sim$ 0.9--1), regardless of the concentration, stellar mass, star formation activity, and local environment (being central or satellite); 2) galaxies with the highest concentration tend to have intrinsic thickness similar to or (in virtually all cases) slightly thinner (i.e. smaller mean $μ_{C/A}$ or equivalently lower triaxiality) than ordinary galaxies, regardless of other properties explored here. This appears to be in contrast with the expectation of the classic merger scenario for high-concentration galaxies. Given the lack of a complete understanding of dwarf-dwarf merger, we cannot draw a definite conclusion about the relevance of mergers in the formation of high-concentration dwarfs. Other mechanisms such as halo spin may also play important roles in the formation of high-concentration dwarf galaxies. △ Less

Submitted 10 September, 2023; originally announced September 2023.

Comments: 12 pages, 8 figures, 2 tables, accepted for publication in ApJ

The role of tidal interactions in the formation of slowly rotating early-type stars in young star clusters

Authors: Chenyu He, Chengyuan Li, Weijia Sun, Richard de Grijs, Lu Li, Jing Zhong, Songmei Qin, Li Chen, Li Wang, Baitian Tang, Zhengyi Shao, Cheng Xu

Abstract: The split main sequences found in the colour-magnitude diagrams of star clusters younger than ~600 Myr are suggested to be caused by the dichotomy of stellar rotation rates of upper main-sequence stars. Tidal interactions have been suggested as a possible explanation of the dichotomy of the stellar rotation rates. This hypothesis proposes that the slow rotation rates of stars along the split main… ▽ More The split main sequences found in the colour-magnitude diagrams of star clusters younger than ~600 Myr are suggested to be caused by the dichotomy of stellar rotation rates of upper main-sequence stars. Tidal interactions have been suggested as a possible explanation of the dichotomy of the stellar rotation rates. This hypothesis proposes that the slow rotation rates of stars along the split main sequences are caused by tidal interactions in binaries. To test this scenario, we measured the variations in the radial velocities of slowly rotating stars along the split main sequence of the young Galactic cluster NGC 2422 (~90 Myr) using spectra obtained at multiple epochs with the Canada-France-Hawai'i Telescope. Our results show that most slowly rotating stars are not radial-velocity variables. Using the theory of dynamical tides, we find that the binary separations necessary to fully or partially synchronise our spectroscopic targets, on time-scales shorter than the cluster age, predict much larger radial velocity variations across multiple-epoch observations, or a much larger radial velocity dispersion at a single epoch, than the observed values. This indicates that tidal interactions are not the dominant mechanism to form slowly rotating stars along the split main sequences. As the observations of the rotation velocity distribution among B- and A-type stars in binaries of larger separations hint at a much stronger effect of braking with age, we discuss the consequences of relaxing the constraints of the dynamical tides theory. △ Less

Submitted 5 September, 2023; originally announced September 2023.

Comments: 14 pages, 10 figures, 2 tables, accepted for publication in MNRAS

Neutron Star vs Quark Star in the Multimessenger Era

Authors: Zheng Cao, Lie-Wen Chen

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

Submitted 31 August, 2023; originally announced August 2023.

Comments: 6 pages, 2 figures

Observation of gamma rays up to 320 TeV from the middle-aged TeV pulsar wind nebula HESS J1849$-$000

Authors: M. Amenomori, S. Asano, Y. W. Bao, X. J. Bi, D. Chen, T. L. Chen, W. Y. Chen, Xu Chen, Y. Chen, Cirennima, S. W. Cui, Danzengluobu, L. K. Ding, J. H. Fang, K. Fang, C. F. Feng, Zhaoyang Feng, Z. Y. Feng, Qi Gao, A. Gomi, Q. B. Gou, Y. Q. Guo, Y. Y. Guo, Y. Hayashi, H. H. He , et al. (93 additional authors not shown)

Abstract: Gamma rays from HESS J1849$-$000, a middle-aged TeV pulsar wind nebula (PWN), are observed by the Tibet air shower array and the muon detector array. The detection significance of gamma rays reaches $4.0\, σ$ and $4.4\, σ$ levels above 25 TeV and 100 TeV, respectively, in units of Gaussian standard deviation $σ$. The energy spectrum measured between $40\, {\rm TeV} < E < 320\, {\rm TeV}$ for the f… ▽ More Gamma rays from HESS J1849$-$000, a middle-aged TeV pulsar wind nebula (PWN), are observed by the Tibet air shower array and the muon detector array. The detection significance of gamma rays reaches $4.0\, σ$ and $4.4\, σ$ levels above 25 TeV and 100 TeV, respectively, in units of Gaussian standard deviation $σ$. The energy spectrum measured between $40\, {\rm TeV} < E < 320\, {\rm TeV}$ for the first time is described with a simple power-law function of ${\rm d}N/{\rm d}E = (2.86 \pm 1.44) \times 10^{-16}(E/40\, {\rm TeV})^{-2.24 \pm 0.41}\, {\rm TeV}^{-1}\, {\rm cm}^{-2}\, {\rm s}^{-1}$. The gamma-ray energy spectrum from the sub-TeV ($E < 1\, {\rm TeV}$) to sub-PeV ($100\, {\rm TeV} < E < 1\, {\rm PeV}$) ranges including the results of previous studies can be modeled with the leptonic scenario, inverse Compton scattering by high-energy electrons accelerated by the PWN of PSR J1849$-$0001. On the other hand, the gamma-ray energy spectrum can also be modeled with the hadronic scenario in which gamma rays are generated from the decay of neutral pions produced by collisions between accelerated cosmic-ray protons and the ambient molecular cloud found in the gamma-ray emitting region. The cutoff energy of cosmic-ray protons $E_{\rm p\, cut}$, cut is estimated at ${\rm log}_{10}(E_{\rm p,\, cut}/{\rm TeV}) = 3.73^{+2.98}_{-0.66}$, suggesting that protons are accelerated up to the PeV energy range. Our study thus proposes that HESS J1849$-$000 should be further investigated as a new candidate for a Galactic PeV cosmic-ray accelerator, PeVatron. △ Less

Submitted 26 August, 2023; originally announced August 2023.

Comments: 10 pages, 2 figures, Accepted for publication from the Astrophysical Journal

Measurement of the Gamma-Ray Energy Spectrum beyond 100 TeV from the HESS J1843$-$033 Region

Authors: M. Amenomori, S. Asano, Y. W. Bao, X. J. Bi, D. Chen, T. L. Chen, W. Y. Chen, Xu Chen, Y. Chen, Cirennima, S. W. Cui, Danzengluobu, L. K. Ding, J. H. Fang, K. Fang, C. F. Feng, Zhaoyang Feng, Z. Y. Feng, Qi Gao, A. Gomi, Q. B. Gou, Y. Q. Guo, Y. Y. Guo, H. H. He, Z. T. He , et al. (91 additional authors not shown)

Abstract: HESS J1843$-$033 is a very-high-energy gamma-ray source whose origin remains unidentified. This work presents, for the first time, the energy spectrum of gamma rays beyond $100\, {\rm TeV}$ from the HESS J1843$-$033 region using the data recorded by the Tibet air shower array and its underground muon detector array. A gamma-ray source with an extension of $0.34^{\circ} \pm 0.12^{\circ}$ is success… ▽ More HESS J1843$-$033 is a very-high-energy gamma-ray source whose origin remains unidentified. This work presents, for the first time, the energy spectrum of gamma rays beyond $100\, {\rm TeV}$ from the HESS J1843$-$033 region using the data recorded by the Tibet air shower array and its underground muon detector array. A gamma-ray source with an extension of $0.34^{\circ} \pm 0.12^{\circ}$ is successfully detected above $25\, {\rm TeV}$ at $(α,\, δ) = (281.09^{\circ}\pm 0.10^{\circ},\, -3.76^{\circ}\pm 0.09^{\circ})$ near HESS J1843$-$033 with a statistical significance of $6.2\, σ$, and the source is named TASG J1844$-$038. The position of TASG J1844$-$038 is consistent with those of HESS J1843$-$033, eHWC J1842$-$035, and LHAASO J1843$-$0338. The measured gamma-ray energy spectrum in $25\, {\rm TeV} < E < 130\, {\rm TeV}$ is described with ${\rm d}N/{\rm d}E = (9.70\pm 1.89)\times 10^{-16} (E/40\, {\rm TeV})^{-3.26\pm 0.30}\, {\rm TeV}^{-1} {\rm cm}^{-2} {\rm s}^{-1}$, and the spectral fit to the combined spectra of HESS J1843$-$033, LHAASO J1843$-$0338, and TASG J1844$-$038 implies the existence of a cutoff at $49.5\pm 9.0\, {\rm TeV}$. Associations of TASG J1844-038 with SNR G28.6$-$0.1 and PSR J1844-0346 are also discussed in detail for the first time. △ Less

Submitted 26 August, 2023; originally announced August 2023.

Comments: 11 pages, 4 figures, 1 table