Showing 1–50 of 608 results for author: Li, L

Novae: An Important Source of Lithium in the Galaxy

Authors: Jun Gao, Chunhua Zhu, Guoliang Lü, Jinlong Yu, Lin Li, Helei Liu, Sufen Guo

Abstract: The source of the Galactic Lithium (Li) has long been a puzzle. With the discovery of Li in novae, extensive research has been conducted. However, there still exists a significant disparity between the observed abundance of lithium in novae and the existing theoretical predictions. Using the Modules for Experiments in Stellar Astrophysics (MESA), we simulate the evolution of nova with element diff… ▽ More The source of the Galactic Lithium (Li) has long been a puzzle. With the discovery of Li in novae, extensive research has been conducted. However, there still exists a significant disparity between the observed abundance of lithium in novae and the existing theoretical predictions. Using the Modules for Experiments in Stellar Astrophysics (MESA), we simulate the evolution of nova with element diffusion and appropriately increased the amount of 3^He in the mixtures. Element diffusion enhances the transport efficiency between the nuclear reaction zone and the convective region on the surface of the white dwarf during nova eruptions, which results in more 7^Be to be transmitted to the white dwarf surface and ultimately ejected. Compared to the previous predictions, the abundance of 7^Be in novae simulated in our model significantly increases. And the result is able to explain almost all observed novae. Using the method of population synthesis, we calculate Li yield in the Galaxy. We find that the Galactic occurrence rate of nova is about 130 yr^{-1}, and about 110M Li produced by nova eruption is ejected into the interstellar medium (ISM). About 73\% of Li in the Galactic ISM originates from novae, and approximately 15\%-20\% of the entire Galaxy. It means that novae are the important source of Li in the Galactic. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 12 pages, 4 figures. Accepted for publication in Astrophysical Journal

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

Probing the Shock Breakout Signal of SN 2024ggi from the Transformation of Early Flash Spectroscopy

Authors: Jujia Zhang, Luc Dessart, Xiaofeng Wang, Qian Zhai, Yi Yang, Liping Li, Han Lin, Giorgio Valerin, Yongzhi Cai, Zhen Guo, Lingzhi Wang, Zeyi Zhao, Zhenyu Wang, Shengyu Yan

Abstract: We present early-time, hour-to-day cadence spectroscopy of the nearby type II supernova (SN II) 2024ggi, which was discovered at a phase when the SN shock just emerged from the red-supergiant (RSG) progenitor star. Over the first few days after the first light, SN 2024ggi exhibited prominent narrow emission lines formed through intense and persistent photoionization of the nearby circumstellar mat… ▽ More We present early-time, hour-to-day cadence spectroscopy of the nearby type II supernova (SN II) 2024ggi, which was discovered at a phase when the SN shock just emerged from the red-supergiant (RSG) progenitor star. Over the first few days after the first light, SN 2024ggi exhibited prominent narrow emission lines formed through intense and persistent photoionization of the nearby circumstellar material (CSM). In the first 63 hours, spectral lines of He, C, N, and O revealed a rapid rise in ionization, as a result of the progressive sweeping-up of the CSM by the shock. The duration of the IIn-like spectra indicates a dense and relatively confined CSM distribution extending up to $\sim 4 \times 10^{14}$ cm. Spectral modeling reveals a CSM mass loss rate at this region exceeding $5 \times\, 10^{-3}${\rm M}_{\odot} yr$^{-1}$ is required to reproduce low-ionization emissions, which dramatically exceeds that of an RSG. Analyzing H$α$ emission shift implies the velocity of the unshocked outer CSM to be between 20 and 40 km s$^{-1}$, matching the typical wind velocity of an RSG. The differences between the inner and outer layers of the CSM and an RSG progenitor highlight a complex mass loss history before the explosion of SN 2024ggi. △ Less

Submitted 29 June, 2024; v1 submitted 11 June, 2024; originally announced June 2024.

Comments: 10 pages and 5 figures in the main text (16 pages and 9 figures in total). Accepted for publication in ApJL

The destiny of NGC 6530: past and future

Authors: Delong Jia, Heng Yu, Zhengyi Shao, Lu Li

Abstract: Studying the structures of open clusters is crucial for understanding stellar evolution and galactic dynamics. Based on Gaia DR3 data, we apply the hierarchical clustering algorithm to a young open cluster NGC 6530 and group its members into 5 substructures. By linear tracing with the kinematic information of their members, we find that: Sub 1 is the core of the cluster. It is expanding slowly. Su… ▽ More Studying the structures of open clusters is crucial for understanding stellar evolution and galactic dynamics. Based on Gaia DR3 data, we apply the hierarchical clustering algorithm to a young open cluster NGC 6530 and group its members into 5 substructures. By linear tracing with the kinematic information of their members, we find that: Sub 1 is the core of the cluster. It is expanding slowly. Sub 2 consists of less bound members, which began escaping from the core about 0.78 Myr ago. Sub 3 is associated with a young star forming region. It will merge with the core after 0.72 Myr; Sub 4, as an outskirt group, is also moving towards the core, but won't end up falling in. While Sub 5 is composed of less-bound members with field contamination. This work reveals the complex internal structure and evolutionary trends of the cluster NGC 6530. It also shows the potential of the hierarchical clustering algorithm in star cluster structure analysis. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 17 pages, 11 figures, accepted for publication in ApJ

Determination method of binary fractions by the integrated spectrum

Authors: F. Zhang, L. Li, Z. Han, X. Wang

Abstract: We need to resolve the individual stars for binary fraction determinations of stellar systems. Therefore, it is not possible to obtain the binary fractions for dense or distant stellar systems. % We proposed a method to determine the binary fraction of a dense or distant stellar system. The method is to first determine the binary fraction variation for any two adjacent regions and then add up thos… ▽ More We need to resolve the individual stars for binary fraction determinations of stellar systems. Therefore, it is not possible to obtain the binary fractions for dense or distant stellar systems. % We proposed a method to determine the binary fraction of a dense or distant stellar system. The method is to first determine the binary fraction variation for any two adjacent regions and then add up those binary fraction variations along the radial direction to obtain the binary fraction for a stellar system. Binary fraction variation is derived by using ten binary fraction-sensitive spectral absorption feature indices (SAFIs) and the binary fraction variation calibrations in terms of these SAFIs. Using this method, we first presented the binary fraction variations for twenty-one Galactic globular clusters (GCs). By comparisons, we find that they agree well with the binary fractions based on the main-sequence fiducial line method by previous studies. This verifies that the above mentioned method is feasible. Next, we presented the binary fraction variations of thirteen Galactic GCs. We gave the relationships between binary fraction and various parameters, and found that binary fraction is negatively correlated with NHB and NRR, binary fraction of some studies is not strongly correlated with NBS, and the number of GCs with large binary fraction is greater at extreme blue horizontal branch population ratio. At last, if we want to obtain more accurate binary fraction, we suggest that the spectroscopic and photometric observations are conducted at an appropriate area interval for a stellar system. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Comments: 12 pages, 6 figures, accepted by MNRAS

All-sky Guide Star Catalog for CSST

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

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

Submitted 3 June, 2024; originally announced June 2024.

Comments: published on RAA

Kinetic temperature of massive star-forming molecular clumps measured with formaldehyde V. The massive filament DR21

Authors: X. Zhao, X. D. Tang, C. Henkel, Y. Gong, Y. Lin, D. L. Li, Y. X. He, Y. P. Ao, X. Lu, T. Liu, Y. Sun, K. Wang, X. P. Chen, J. Esimbek, J. J. Zhou, J. W. Wu, J. J. Qiu, X. W. Zheng, J. S. Li, C. S. Luo, Q. Zhao

Abstract: The kinetic temperature structure of the massive filament DR21 has been mapped using the IRAM 30 m telescope. This mapping employed the para-H$_2$CO triplet ($J_{\rm K_aK_c}$ = 3$_{03}$--2$_{02}$, 3$_{22}$--2$_{21}$, and 3$_{21}$--2$_{20}$) on a scale of $\sim$0.1 pc. By modeling the averaged line ratios of para-H$_{2}$CO with RADEX under non-LTE assumptions, the kinetic temperature of the dense g… ▽ More The kinetic temperature structure of the massive filament DR21 has been mapped using the IRAM 30 m telescope. This mapping employed the para-H$_2$CO triplet ($J_{\rm K_aK_c}$ = 3$_{03}$--2$_{02}$, 3$_{22}$--2$_{21}$, and 3$_{21}$--2$_{20}$) on a scale of $\sim$0.1 pc. By modeling the averaged line ratios of para-H$_{2}$CO with RADEX under non-LTE assumptions, the kinetic temperature of the dense gas was derived at a density of $n$(H$_{2}$) = 10$^{5}$ cm$^{-3}$. The para-H$_2$CO lines reveal significantly higher temperatures than NH$_3$ (1,1)/(2,2) and FIR wavelengths. The dense clumps appear to correlate with the notable kinetic temperature. Among the four dense cores (N44, N46, N48, and N54), temperature gradients are observed on a scale of $\sim$0.1-0.3 pc. This suggests that the warm dense gas is influenced by internal star formation activity. With the exception of N54, the temperature profiles of these cores were fitted with power-law indices ranging from $-$0.3 to $-$0.5. This indicates that the warm dense gas is heated by radiation emitted from internally embedded protostar(s) and/or clusters. While there is no direct evidence supporting the idea that the dense gas is heated by shocks resulting from a past explosive event in the DR21 region, our measurements toward the DR21W1 region provide compelling evidence that the dense gas is indeed heated by shocks originating from the western DR21 flow. Higher temperatures appear to be associated with turbulence. The physical parameters of the dense gas in the DR21 filament exhibit a remarkable similarity to the results obtained in OMC-1 and N113. This may imply that the physical mechanisms governing the dynamics and thermodynamics of dense gas traced by H$_{2}$CO in diverse star formation regions may be dominated by common underlying principles despite variations in specific environmental conditions. (abbreviated) △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures, 3 tabels. Accepted for publication by Astronomy & Astrophysics

Evidence for a close-in tertiary orbiting around the Algol-type system HZ Dra with tidal splitting and spots activities

Authors: Ping Li, Wen-Ping Liao, Sheng-Bang Qian, Lin-Jia Li, Qi-Bin Sun, Xiang-Dong Shi, Azizbek Matekov, Qi-Huan Zeng, Zhao-Long Deng, Xiao-Hui Fang

Abstract: We reported a cyclic variation of $O-C$ diagram with a semi-amplitude of 0.0033 days and a period of 1.05 years for the pulsating eclipsing binary HZ Dra. The cyclic variation can be explained by the light travel-time effect via the presence of a close-in third body orbiting around HZ Dra in an elliptical orbit with a maximum semi-major axis of 0.92 au. Based on the W-D code, the contribution of t… ▽ More We reported a cyclic variation of $O-C$ diagram with a semi-amplitude of 0.0033 days and a period of 1.05 years for the pulsating eclipsing binary HZ Dra. The cyclic variation can be explained by the light travel-time effect via the presence of a close-in third body orbiting around HZ Dra in an elliptical orbit with a maximum semi-major axis of 0.92 au. Based on the W-D code, the contribution of the third light to the total system is determined to be 29 $\%$, which is in agreement with the estimated value. Our light curve modelling indicates an evolving hot and cool spot on the surface of the primary and secondary components, respectively. Their positions are roughly symmetrical to the inner Lagrangian point L1, which could be used to explain the variation in the O$^{'}$Connell effect. Our frequency analysis detects 1 radial p-mode, 7 non-radial p-modes and 1 non-radial g-mode. In addition, a total of 6 multiplets are identified, spaced by the orbital frequency, which can be explained as a tidally split mode caused by the equilibrium tides of the close binary system with a circular orbit. These pulsating features suggest that the primary of HZ Dra is a $δ$ Scuti star, pulsating in both p- and g-mode and influenced by tidal forces. △ Less

Submitted 30 May, 2024; v1 submitted 25 May, 2024; originally announced May 2024.

Comments: 15 pages, 11 figures. The first version is Under Review in ApJ

SN 2019tua : A Type IIb Supernova with Multiple Bumps in the Light Curves

Authors: Xin-Bo Huang, Xiang-Gao Wang, Long Li, Li-Ping Xin, Jing Wang, Tian-Ci Zheng, Qi Wang, Hui-Ya Liu, Zi-Min Zhou, Xiao-meng Lu, jian-yan Wei, En-Wei Liang

Abstract: We present photometric and spectroscopic observations and analysis of the type IIb supernova (SN) SN 2019tua, which exhibits multiple bumps in its declining light curves between 40 and 65 days after discovery. SN 2019tua shows a time to peak of about 25 days similar to other type IIb SNe. Our observations indicate a decrease in its brightness of about 1 magnitude in the 60 days after the peak. At… ▽ More We present photometric and spectroscopic observations and analysis of the type IIb supernova (SN) SN 2019tua, which exhibits multiple bumps in its declining light curves between 40 and 65 days after discovery. SN 2019tua shows a time to peak of about 25 days similar to other type IIb SNe. Our observations indicate a decrease in its brightness of about 1 magnitude in the 60 days after the peak. At about days 50, and 60, its multiband light curves exhibit bumpy behavior. The complex luminosity evolution of SN 2019tua could not be well modeled with a single currently popular energy source model, e.g., radioactive decay of $^{56}$Ni, magnetar, interaction between the ejecta and a circumstellar shell. Even though the magnetar model has a smaller $χ^2 / \text{dof}$ value, the complex changes in SN 2019tua's brightness suggest that more than one physical process might be involved. We propose a hybrid CSM interaction plus $^{56}$Ni model to explain the bolometric light curve (LC) of SN 2019tua. The fitting results show that the ejecta mass $M_{\rm ej} \approx 2.4~M_\odot$, the total CSM mass $M_{\rm CSM} \approx 1.0~M_\odot$, and the $^{56}$Ni mass $M_{\rm Ni} \approx 0.4~M_\odot$. The total kinetic energy of the ejecta is $E_k\approx 0.5 \times 10^{51}\rm~erg$. Pre-existing multiple shells suggest that the progenitor of SN 2019tua experienced mass ejections within approximately $\sim6 - 44$ years prior to the explosion. △ Less

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

Comments: Accepted for publication in ApJ. 24 pages, 17 figures, 6 tables

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 population synthesis of Wolf-Rayet stars involving binary merger channels

Authors: Zhuowen Li, Chunhua Zhu, Guoliang Lü, Lin Li, Helei Liu, Sufen Guo, Jinlong Yu, Xizhen Lu

Abstract: Wolf-Rayet stars (WRs) are very important massive stars. However, their origin and the observed binary fraction within the entire WR population are still debated. We investigate some possible merger channels for the formation of WRs, including main sequence (MS)/ Hertzsprung Gap (HG) + MS, He + HG/ Giant Branch (GB). We find that many products produced via binary merger can evolve into WRs, the MS… ▽ More Wolf-Rayet stars (WRs) are very important massive stars. However, their origin and the observed binary fraction within the entire WR population are still debated. We investigate some possible merger channels for the formation of WRs, including main sequence (MS)/ Hertzsprung Gap (HG) + MS, He + HG/ Giant Branch (GB). We find that many products produced via binary merger can evolve into WRs, the MS/ HG + MS merger channel can explain WRs with luminosities higher than $\sim 10^{5.4}$\,L$_{\odot}$, while the He + HG/ GB merger channel can explain low-luminosity WRs in the range of $10^{4.7}$\,L$_{\odot}$\,$\sim$\,$10^{5.5}$\,L$_{\odot}$. In the population synthesis analysis of WRs, we assume an initial binary fraction ($f_{\rm ini,bin}$) of 50\% and 100\% for massive stars. We also assume that MS/ HG + MS merger products are non-rotating or rapidly rotating ($ω/ω_{\rm crit}=0.8$). In different cases, the calculated single fractions of WRs range from $22.2\%$ to $60.6\%$ in the Milky Way (MW) and from $8.3\%$ to $70.9\%$ in the Large Magellanic Cloud (LMC). The current observations fall within the range of our calculations. When the merger product of MS/HG+MS rotates rapidly, we estimate that there are approximately 1015 to 1396 WRs in the MW and 128 to 204 WRs in the LMC. Our model also roughly reproduces the observed single-peak luminosity distribution of WRs in the MW. However, the weak bimodal luminosity distribution observed in the LMC is not reproduced in our model. We assess that this may be due to the model underestimating the mass-loss rate in the LMC. In conclusion, we consider that the binary merger is significant formation channel for WR formation, and can explain the observed high fraction of the single WRs in the total population. △ Less

Submitted 19 May, 2024; originally announced May 2024.

Comments: 14 pages, 8 figures , Accepted to APJ

Ten Supernova-rise in Binary Driven Gamma-ray Bursts

Authors: R. Ruffini, C. L. Bianco, Liang Li, M. T. Mirtorabi, R. Moradi, F. Rastegarnia, J. A. Rueda, S. R. Zhang, Y. Wang

Abstract: The observation of a gamma-ray burst (GRB) associated with a supernova (SN) coincides remarkably with the energy output from a binary system comprising a very massive carbon-oxygen (CO) core and an associated binary neutron star (NS) by the Binary-Driven Hypernova (BdHN) model. The dragging effect in the late evolution of such systems leads to co-rotation, with binary periods on the order of minut… ▽ More The observation of a gamma-ray burst (GRB) associated with a supernova (SN) coincides remarkably with the energy output from a binary system comprising a very massive carbon-oxygen (CO) core and an associated binary neutron star (NS) by the Binary-Driven Hypernova (BdHN) model. The dragging effect in the late evolution of such systems leads to co-rotation, with binary periods on the order of minutes, resulting in a very fast rotating core and a binary NS companion at a distance of $\sim 10^5$ km. Such a fast-rotating CO core, stripped of its hydrogen and helium, undergoes gravitational collapse and, within a fraction of seconds, leads to a supernova (SN) and a newly born, fast-spinning neutron star ($ν$NS), we name the emergence of the SN and the $ν$NS as the SN-rise and $ν$NS-rise. Typically, the SN energies range from $10^{51}$ to $10^{53}$ erg. We address this issue by examining 10 cases of Type-I BdHNe, the most energetic ones, in which SN accretion onto the companion NS leads to the formation of a black hole (BH). In all ten cases, the energetics of the SN events are estimated, ranging between $0.18$ and $12 \times 10^{52}$ erg. Additionally, in all 8 sources at redshift $z$ closer than $4.61$, a clear thermal blackbody component has been identified, with temperatures between $6.2$ and $39.99$ keV, as a possible signature of pair-driven SN. The triggering of the X-ray afterglow induced by the $ν$NS-rise are identified in three cases at high redshift where early X-ray observations are achievable, benefits from the interplay of cosmological effects. △ Less

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

Comments: The previous version was uploaded on arXiv before being submitted to any Journal, this new version matches the one actually submitted for publication

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

An Optical Gamma-Ray Burst Catalogue with Measured Redshift PART I: Data Release of 535 Gamma-Ray Bursts and Colour Evolution

Authors: M. G. Dainotti, B. De Simone, R. F. Mohideen Malik, V. Pasumarti, D. Levine, N. Saha, B. Gendre, D. Kido, A. M. Watson, R. L. Becerra, S. Belkin, S. Desai, A. C. C. do E. S. Pedreira, U. Das, L. Li, S. R. Oates, S. B. Cenko, A. Pozanenko, A. Volnova, Y. -D. Hu, A. J. Castro-Tirado, N. B. Orange, T. J. Moriya, N. Fraija, Y. Niino , et al. (27 additional authors not shown)

Abstract: We present the largest optical photometry compilation of Gamma-Ray Bursts (GRBs) with redshifts ($z$). We include 64813 observations of 535 events (including upper limits) from 28 February 1997 up to 18 August 2023. We also present a user-friendly web tool \textit{grbLC} which allows users the visualization of photometry, coordinates, redshift, host galaxy extinction, and spectral indices for each… ▽ More We present the largest optical photometry compilation of Gamma-Ray Bursts (GRBs) with redshifts ($z$). We include 64813 observations of 535 events (including upper limits) from 28 February 1997 up to 18 August 2023. We also present a user-friendly web tool \textit{grbLC} which allows users the visualization of photometry, coordinates, redshift, host galaxy extinction, and spectral indices for each event in our database. Furthermore, we have added a Gamma Ray Coordinate Network (GCN) scraper that can be used to collect data by gathering magnitudes from the GCNs. The web tool also includes a package for uniformly investigating colour evolution. We compute the optical spectral indices for 138 GRBs for which we have at least 4 filters at the same epoch in our sample and craft a procedure to distinguish between GRBs with and without colour evolution. By providing a uniform format and repository for the optical catalogue, this web-based archive is the first step towards unifying several community efforts to gather the photometric information for all GRBs with known redshifts. This catalogue will enable population studies by providing light curves (LCs) with better coverage since we have gathered data from different ground-based locations. Consequently, these LCs can be used to train future LC reconstructions for an extended inference of the redshift. The data gathering also allows us to fill some of the orbital gaps from Swift in crucial points of the LCs, e.g., at the end of the plateau emission or where a jet break is identified. △ Less

Submitted 3 June, 2024; v1 submitted 3 May, 2024; originally announced May 2024.

Comments: 20 pages, 16 figures, 2 tables. Submitted to MNRAS, this version matches the third revision. The Online Materials and data will be available after the publication

Comparison of Ion-Proton Differential Speed between ICMEs and Solar Wind near 1 au

Authors: Xuechao Zhang, Hongqiang Song, Chengxiao Zhang, Hui Fu, Leping Li, Jinrong Li, Xiaoqian Wang, Rui Wang, Yao Chen

Abstract: The elemental abundance of ICMEs and solar wind near 1 au is often adopted to represent the abundance in the corresponding coronal sources. However, the absolute abundance of heavy ions (relative to hydrogen) near 1 au might be different from the coronal abundance due to the ion-proton differential speed ($V_{ip}$). To illustrate the $V_{ip}$ characteristics and explore whether it influences the a… ▽ More The elemental abundance of ICMEs and solar wind near 1 au is often adopted to represent the abundance in the corresponding coronal sources. However, the absolute abundance of heavy ions (relative to hydrogen) near 1 au might be different from the coronal abundance due to the ion-proton differential speed ($V_{ip}$). To illustrate the $V_{ip}$ characteristics and explore whether it influences the absolute abundance analysis for ICMEs and solar wind, we perform a statistical study on the $V_{ip}$ for He$^{2+}$, C$^{5+}$, O$^{6+}$, and Fe$^{10+}$ in both ICMEs and solar wind based on measurements of Advanced Composition Explorer. The results show that the $V_{ip}$ is negligible within ICMEs and slow solar wind ($<$ 400 km s$^{-1}$), while obvious in the intermediate (400 -- 600 km s$^{-1}$) and fast wind ($>$ 600 km s$^{-1}$). Previous studies showed that the $V_{ip}$ in ICMEs keeps negligible during propagation from 0.3 to 5 au, but in solar wind it increases with the decreasing heliocentric distance. Therefore, it might be questionable to infer the absolute abundance of coronal sources through in-situ abundance near 1 au for solar wind. Fortunately, the ion-oxygen (O$^{6+}$) differential speed ($V_{io}$) is negligible for He$^{2+}$, C$^{5+}$, and Fe$^{10+}$ within both ICMEs and solar wind, and previous studies suggested that the $V_{io}$ does not vary significantly with the heliocentric distance. This indicates that various heavy ions always flow at the same bulk speed and their relative abundance (relative to oxygen) near 1 au can represent the coronal abundance for both ICMEs and solar wind. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 14 pages, 8 figures

Gravitational wave radiation from magnetar-driven supernovae

Authors: Lang Xie, Hong-Yu Gong, Long LI, Da-Ming Wei, J. L. Han

Abstract: Rapidly spinning magnetars are potential candidates for the energy source of supernovae (SNe) and gamma-ray bursts and the most promising sources for continuous gravitational waves (GWs) detected by ground-based GW detectors. Continuous GWs can be radiated from magnetars due to magnetic-induced deformation or fluid oscillations, compatible with magnetic dipole (MD) radiation for spin-down energy.… ▽ More Rapidly spinning magnetars are potential candidates for the energy source of supernovae (SNe) and gamma-ray bursts and the most promising sources for continuous gravitational waves (GWs) detected by ground-based GW detectors. Continuous GWs can be radiated from magnetars due to magnetic-induced deformation or fluid oscillations, compatible with magnetic dipole (MD) radiation for spin-down energy. In this paper, we investigate the diverse light curves of magnetar-driven SNe in the scenario that the spin-down is dominated by GW radiation and/or MD radiation. By simulating the light curves of SNe and employing the Markov chain Monte Carlo method, we constrain the parameters of the magnetars and SN explosions and show that the signature of GW radiation may be indicated by the bolometric luminosity curves of SNe Ic-BL 2007ru and 2009bb. We find that the ellipticity of magnetars in the order of $10^{-3}$ can be induced by the magnetic field of $\sim 10^{16} \mbox{ G}$. If such continuous GWs associated with SNe can be detected in the future by the Advanced LIGO and Virgo detectors, this would be a smoking gun for a magnetar engine powering SNe. △ Less

Submitted 28 April, 2024; originally announced April 2024.

Comments: 13pages, 9 figures, Accepted for publication in ApJ

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

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

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

Submitted 25 April, 2024; originally announced April 2024.

Comments: 41 pages, 8 figures, 7 tables

Statistical analysis of pulsar flux density distribution

Authors: H. W. Xu, R. S. Zhao, Erbil Gugercinoglu, H. Liu, D. Li, P. Wang, C. H. Niu, C. Miao, X. Zhu, R. W. Tian, W. L. Li, S. D. Wang, Z. F. Tu, Q. J. Zhi, S. J. Dang, L. H. Shang, S. Xiao

Abstract: This study presents a comprehensive analysis of the spectral properties of 886 pulsars across a wide frequency range from 20MHz to 343.5GHz, including a total of 86 millisecond pulsars. The majority of the pulsars exhibit power-law behavior in their spectra, although some exceptions are observed. Five different spectral models, namely simple power-law, broken power-law, low-frequency turn-over, hi… ▽ More This study presents a comprehensive analysis of the spectral properties of 886 pulsars across a wide frequency range from 20MHz to 343.5GHz, including a total of 86 millisecond pulsars. The majority of the pulsars exhibit power-law behavior in their spectra, although some exceptions are observed. Five different spectral models, namely simple power-law, broken power-law, low-frequency turn-over, high-frequency cut-off, and double turn-over, were employed to explore the spectral behaviors. The average spectral index for pulsars modeled with a simple power-law is found to be -1.64 +/-0.80, consistent with previous studies. Additionally, significant correlations between the spectral index and characteristic parameters are observed particularly in millisecond pulsars, while no strong correlation is observed in normal pulsars. Different models show variations in the most influential characteristic parameters associated with the spectral index, indicating diverse dominant radiation mechanisms in millisecond pulsars.Finally, this study identifies 22 pulsars of the Gigahertz-peaked Spectra (GPS) type for the first time based on the Akaike information criterion. △ Less

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

Comments: 39 papers,17figures

Can AI Understand Our Universe? Test of Fine-Tuning GPT by Astrophysical Data

Authors: Yu Wang, Shu-Rui Zhang, Aidin Momtaz, Rahim Moradi, Fatemeh Rastegarnia, Narek Sahakyan, Soroush Shakeri, Liang Li

Abstract: ChatGPT has been the most talked-about concept in recent months, captivating both professionals and the general public alike, and has sparked discussions about the changes that artificial intelligence (AI) will bring to the world. As physicists and astrophysicists, we are curious about if scientific data can be correctly analyzed by large language models (LLMs) and yield accurate physics. In this… ▽ More ChatGPT has been the most talked-about concept in recent months, captivating both professionals and the general public alike, and has sparked discussions about the changes that artificial intelligence (AI) will bring to the world. As physicists and astrophysicists, we are curious about if scientific data can be correctly analyzed by large language models (LLMs) and yield accurate physics. In this article, we fine-tune the generative pre-trained transformer (GPT) model by the astronomical data from the observations of galaxies, quasars, stars, gamma-ray bursts (GRBs), and the simulations of black holes (BHs), the fine-tuned model demonstrates its capability to classify astrophysical phenomena, distinguish between two types of GRBs, deduce the redshift of quasars, and estimate BH parameters. We regard this as a successful test, marking the LLM's proven efficacy in scientific research. With the ever-growing volume of multidisciplinary data and the advancement of AI technology, we look forward to the emergence of a more fundamental and comprehensive understanding of our universe. This article also shares some interesting thoughts on data collection and AI design. Using the approach of understanding the universe - looking outward at data and inward for fundamental building blocks - as a guideline, we propose a method of series expansion for AI, suggesting ways to train and control AI that is smarter than humans. △ Less

Submitted 14 April, 2024; originally announced April 2024.

Comments: 27 pages, 7 figures. Comments welcome

Eruption of a million-Kelvin warm magnetic flux rope on the Sun

Authors: Leping Li, Hongqiang Song, Hardi Peter, Lakshmi Pradeep Chitta, Xin Cheng, Zhentong Li, Guiping Zhou

Abstract: Solar magnetic flux rope (MFR) plays a central role in the physics of coronal mass ejections (CMEs). It mainly includes a cold filament at typical chromospheric temperatures (10000 K) and a hot channel at high coronal temperatures (10 MK). The warm MFR at quiescent coronal temperatures of a million Kelvin is, however, rarely reported. In this study, using multiwavelength images from Atmospheric Im… ▽ More Solar magnetic flux rope (MFR) plays a central role in the physics of coronal mass ejections (CMEs). It mainly includes a cold filament at typical chromospheric temperatures (10000 K) and a hot channel at high coronal temperatures (10 MK). The warm MFR at quiescent coronal temperatures of a million Kelvin is, however, rarely reported. In this study, using multiwavelength images from Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) and Extreme Ultraviolet Imager (EUVI) on board the Solar Terrestrial Relations Observatory-A (STEREO-A), we present an eruption of a warm channel, that represents an MFR with quiescent coronal temperatures (0.6-2.5 MK). On 2022 May 8, we observed the failed eruption of a hot channel, with the average temperature and emission measure (EM) of 10 MK and 1.1*1028 cm^-5, using AIA high-temperature images in active region (AR) 13007. This failed eruption was associated with a C8.2 flare, with no CME. Subsequently, we observed a warm channel that appeared in AIA and EUVI low-temperature images, rather than AIA high-temperature images. It then erupted, and transformed toward a semi-circular shape. An associated C2.1 flare, along with the signatures of magnetic reconnection in AIA high-temperature images, were identified. Additionally, we observed a CME associated with this event. Compared with the hot channel, the warm channel is cooler and rarer with the average temperature and EM of 1.7 (1.6) MK and 2.0*1026 (2.3*1026) cm^-5. All the results suggest an unambiguous observation of the million-Kelvin warm MFR, that erupted as a CME, and fill a gap in the temperature domain of coronal MFRs. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 18 pages, 7 figures, 2 tables, accepted for publication in ApJ

Multipolar Electromagnetic Emission of Newborn Magnetar

Authors: Yu Wang, Rahim Moradi, Liang Li

Abstract: It is generally recognized that the electromagnetic multipolar emission from magnetars can be used to explain radiation from Soft Gamma Repeaters (SGRs) or Anomalous X-ray Pulsars (AXPs), but they have little impact on the spindown of magnetars. We here present a comprehensive analytical solution for the neutron star multipolar electromagnetic fields and their associated expected luminosities. We… ▽ More It is generally recognized that the electromagnetic multipolar emission from magnetars can be used to explain radiation from Soft Gamma Repeaters (SGRs) or Anomalous X-ray Pulsars (AXPs), but they have little impact on the spindown of magnetars. We here present a comprehensive analytical solution for the neutron star multipolar electromagnetic fields and their associated expected luminosities. We find that for newborn millisecond magnetars, the spin-down luminosity from higher multipolar components can match or even exceed that from the dipole component. Such high-intensity radiation will undoubtedly affect related astrophysical phenomena at the birth of a magnetar. We show that the spin-down luminosity from multipoles can well explain the majority of Gamma-Ray Bursts (GRBs) afterglows, from the plateau starting at several hundred seconds till the normal decay phase lasting for many years. The fitted magnetar parameters for GRB afterglows are all typical values, with spins in the millisecond range and magnetic field strengths in the order of $10^{14} - 10^{15}$ Gauss. Our results in turn, provide support for the hypothesis that GRBs originate from the birth of magnetars with a few millisecond period, thus deepening our understanding of the complex magnetic field structure and the equation of state of magnetars. △ Less

Submitted 14 April, 2024; originally announced April 2024.

Discovery of a new IW And-type dwarf nova with both tilted disk and tidal instability

Authors: Yongkang Sun, Xin Li, Qige Ao, Wenyuan Cui, Bowen Zhang, Yang Huang, Jianrong Shi, Linlin Li, Jifeng Liu

Abstract: IW And-type dwarf novae are anomalous Z Cam stars featured with outbursts happening during standstill states, which are not expected in the standard disk instability model. The physical mechanisms for these variations remain unclear. In this study, we report the discovery of a new candidate IW And-type dwarf nova J0652+2436, identified with its frequent outbursts from the slowly rising standstill… ▽ More IW And-type dwarf novae are anomalous Z Cam stars featured with outbursts happening during standstill states, which are not expected in the standard disk instability model. The physical mechanisms for these variations remain unclear. In this study, we report the discovery of a new candidate IW And-type dwarf nova J0652+2436, identified with its frequent outbursts from the slowly rising standstill states. Luckily, the TESS observations during a long standstill state and the earlier K2 observations give a chance to find the orbital and negative superhump period in the light curve of J0652+2436, allowing the measurement of its mass ratio of 0.366. This mass ratio is marginally possible for the tidal instability to set in according to previous SPH simulations. Thus, we propose that the outbursts in J0652+2436 are likely to be caused by the growing accretion disk during standstills, in favor of the previous hypothesis of the mechanisms lying in all IW And stars. We conclude that J0652+2436 might be the first IW And star with both a precessing tilted disk and tidal instability, which will be an important laboratory for studying the accretion disk dynamics and help understand IW And phenomenon. △ Less

Submitted 13 April, 2024; originally announced April 2024.

Comments: 12 pages, 12 figures, accepted to MNRAS

LHAASO-KM2A detector simulation using Geant4

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

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

Submitted 7 April, 2024; originally announced April 2024.

Clues in Bailey Diagrams Indicate that the Horizontal Branch Evolutionary Effect Is the Direct Reason for Oosterhoff Phenomenon

Authors: L. -J. Li, S. -B. Qian, L. -Y. Zhu, X. -D. Shi, W. -P. Liao

Abstract: Noticing the potential relationship between the pulsation amplitude of fundamental mode RR Lyrae stars and their average effective temperatures, we derived the relation between them from observational data. Using these relationships, we present the zero-age horizontal branch lines and corresponding evolutionary tracks in the Bailey diagrams. We first find that the trends of zero-age horizontal bra… ▽ More Noticing the potential relationship between the pulsation amplitude of fundamental mode RR Lyrae stars and their average effective temperatures, we derived the relation between them from observational data. Using these relationships, we present the zero-age horizontal branch lines and corresponding evolutionary tracks in the Bailey diagrams. We first find that the trends of zero-age horizontal branch lines are consistent with the distributions of Oosterhoff I RR Lyrae stars, and there is a period shift effect in relatively metal-rich samples ([Fe/H] > -1). Main results show that evolutionary effects are the direct cause of the Oosterhoff dichotomy. The vast majority of Oosterhoff I RR Lyrae stars are in the early stage of horizontal branch evolution, while most of the Oosterhoff II stars are near the end of the horizontal branch phase. Our Bailey diagrams can intuitively explain some Oosterhoff phenomena (e.g., long-period RR Lyrae stars in the metal-rich globular clusters NGC 6388 and NGC 6441), and also have the potential to become an important tool for in-depth study of RR Lyrae stars. Further research is needed on the Oosterhoff dichotomy to explore the formations and evolutionary processes of the Milky Way and its substructures. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 4 figures

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)

Suppression of Star Formation in Galaxy Pairs

Authors: Shuai Feng, Shi-Yin Shen, Fang-Ting Yuan, Wen-Xin Zhong, Wen-Yuan Cui, Lin-Lin Li

Abstract: We investigate the suppression of star formation in galaxy pairs based on the isolated galaxy pair sample derived from the SDSS survey. By comparing the star formation rate between late-type galaxies in galaxy pairs and those in the isolated environment, we detect the signal of star formation suppression in galaxy pairs at $d_p < 100$kpc and $200$kpc$ < d_p < 350$kpc. The occurrence of star format… ▽ More We investigate the suppression of star formation in galaxy pairs based on the isolated galaxy pair sample derived from the SDSS survey. By comparing the star formation rate between late-type galaxies in galaxy pairs and those in the isolated environment, we detect the signal of star formation suppression in galaxy pairs at $d_p < 100$kpc and $200$kpc$ < d_p < 350$kpc. The occurrence of star formation suppression in these late-type galaxies requires their companion galaxies to have an early-type morphology ($n_s > 2.5$). Star formation suppression in wide galaxy pairs with $200$kpc$ < d_p < 350$kpc mainly occurs in massive late-type galaxies, while in close galaxy pairs with $d_p < 100$kpc, it only appears in late-type galaxies with a massive companion ( $\log M_\star > 11.0$), nearly independent of their own stellar mass. Based on these findings, we infer that star formation suppression in wide galaxy pairs is actually a result of galaxy conformity, while in close galaxy pairs, it stems from the influence of hot circum-galactic medium surrounding companion galaxies. △ Less

Submitted 14 March, 2024; originally announced March 2024.

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

Super-adiabatic Temperature Gradient at Jupiter's Equatorial Zone and Implications for the Water Abundance

Authors: Cheng Li, Michael Allison, Sushil Atreya, Shawn Brueshaber, Leigh N. Fletcher, Tristan Guillot, Liming Li, Jonathan Lunine, Yamila Miguel, Glenn Orton, Paul Steffes, J. Hunter Waite, Michael H. Wong, Steven Levin, Scott Bolton

Abstract: The temperature structure of a giant planet was traditionally thought to be an adiabat assuming convective mixing homogenizes entropy. The only in-situ measurement made by the Galileo Probe detected a near-adiabatic temperature structure within one of Jupiter's 5$μ$m hot spots with small but definite local departures from adiabaticity. We analyze Juno's microwave observations near Jupiter's equato… ▽ More The temperature structure of a giant planet was traditionally thought to be an adiabat assuming convective mixing homogenizes entropy. The only in-situ measurement made by the Galileo Probe detected a near-adiabatic temperature structure within one of Jupiter's 5$μ$m hot spots with small but definite local departures from adiabaticity. We analyze Juno's microwave observations near Jupiter's equator (0 ~ 5$^o$N) and find that the equatorial temperature structure is best characterized by a stable super-adiabatic temperature profile rather than an adiabatic one. Water is the only substance with sufficient abundance to alter the atmosphere's mean molecular weight and prevent dynamic instability if a super-adiabatic temperature gradient exists. Thus, from the super-adiabaticity, our results indicate a water concentration (or the oxygen to hydrogen ratio) of about 4.9 times solar with a possible range of 1.5 ~ 8.3 times solar in Jupiter's equatorial region. △ Less

Submitted 8 March, 2024; originally announced March 2024.

Comments: 46 pages, 11 figures, accepted by Icraus

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

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

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

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

CSS_J154915.7+375506: A low-mass-ratio marginal contact binary system with a hierarchical third body

Authors: Jin-Feng Wu, Li-Ying Zhu, Azizbek Matekov, Lin-jia Li, Shuhrat Ehgamberdiev, Ildar Asfandiyarov, Jiang-Jiao Wang, Jia Zhang, Fang-Bin Meng

Abstract: We presented the multi-filter light curves of CSS_J154915.7+375506 inaugurally, which were observed by the 1.5 m AZT-22 telescope at Maidanak Astronomical Observatory. A low-resolution spectrum obtained by LAMOST reveals it is an A-type close binary. By analyzing the BVRI total-eclipse light curves, we are able to derive a reliable photometric solution for this system, which indicates that CSS_J15… ▽ More We presented the multi-filter light curves of CSS_J154915.7+375506 inaugurally, which were observed by the 1.5 m AZT-22 telescope at Maidanak Astronomical Observatory. A low-resolution spectrum obtained by LAMOST reveals it is an A-type close binary. By analyzing the BVRI total-eclipse light curves, we are able to derive a reliable photometric solution for this system, which indicates that CSS_J154915.7+375506 is an extremely low-mass-ratio (q=0.138) marginal contact binary system. The location in the HR diagram shows that its secondary component with a much smaller mass is the more evolved one, indicating the mass ratio reversal occurred. The present secondary component had transferred a significant amount of mass to the present primary one. By the combination of a total of 20 times of minimum, we investigated its O-C curve. A periodic oscillation and a possible period decrease have been detected. As the period decreases, the system will evolve towards the contact phase. This makes CSS\_J154915.7+375506 a valuable case to study the formation scenario of contact binaries through mass reversal. The periodic oscillation suggested a third body with a minimal mass of $0.91\,M_{\odot}$, which is larger than that of the less massive component in the central binary. This implies that the secondary body was not replaced by the third body during early stellar interactions, indicating that it is a fossil system and retains its original dynamical information. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 8 pages, 6 figures

AT2023lli: A Tidal Disruption Event with Prominent Optical Early Bump and Delayed Episodic X-ray Emission

Authors: Shifeng Huang, Ning Jiang, Jiazheng Zhu, Yibo Wang, Tinggui Wang, Shan-Qin Wang, Wen-Pei Gan, En-Wei Liang, Yu-Jing Qin, Zheyu Lin, Lin-Na Xu, Min-Xuan Cai, Ji-An Jiang, Xu Kong, Jiaxun Li, Long Li, Jian-Guo Wang, Ze-Lin Xu, Yongquan Xue, Ye-Fei Yuan, Jingquan Cheng, Lulu Fan, Jie Gao, Lei Hu, Weida Hu , et al. (20 additional authors not shown)

Abstract: High-cadence, multiwavelength observations have continuously revealed the diversity of tidal disruption events (TDEs), thus greatly advancing our knowledge and understanding of TDEs. In this work, we conducted an intensive optical-UV and X-ray follow-up campaign of TDE AT2023lli, and found a remarkable month-long bump in its UV/optical light curve nearly two months prior to maximum brightness. The… ▽ More High-cadence, multiwavelength observations have continuously revealed the diversity of tidal disruption events (TDEs), thus greatly advancing our knowledge and understanding of TDEs. In this work, we conducted an intensive optical-UV and X-ray follow-up campaign of TDE AT2023lli, and found a remarkable month-long bump in its UV/optical light curve nearly two months prior to maximum brightness. The bump represents the longest separation time from the main peak among known TDEs to date. The main UV/optical outburst declines as $t^{-4.10}$, making it one of the fastest decaying optically selected TDEs. Furthermore, we detected sporadic X-ray emission 30 days after the UV/optical peak, accompanied by a reduction in the period of inactivity. It is proposed that the UV/optical bump could be caused by the self-intersection of the stream debris, whereas the primary peak is generated by the reprocessed emission of the accretion process. In addition, our results suggest that episodic X-ray radiation during the initial phase of decline may be due to the patched obscurer surrounding the accretion disk, a phenomenon associated with the inhomogeneous reprocessing process. The double TDE scenario, in which two stars are disrupted in sequence, is also a possible explanation for producing the observed early bump and main peak. We anticipate that the multicolor light curves of TDEs, especially in the very early stages, and the underlying physics can be better understood in the near future with the assistance of dedicated surveys such as the deep high-cadence survey of the 2.5-meter Wide Field Survey Telescope (WFST). △ Less

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

Comments: 14 pages, 8 figures,accepted for publication by ApJL

Magnetar as the Central Engine of AT2018cow: Optical, Soft X-Ray, and Hard X-Ray Emission

Authors: Long Li, Shu-Qing Zhong, Di Xiao, Zi-Gao Dai, Shi-Feng Huang, Zhen-Feng Sheng

Abstract: AT2018cow is the most extensively observed and widely studied fast blue optical transient to date; its unique observational properties challenge all existing standard models. In this paper, we model the luminosity evolution of the optical, soft X-ray, and hard X-ray emission, as well as the X-ray spectrum of AT2018cow with a magnetar-centered engine model. We consider a two-zone model with a strip… ▽ More AT2018cow is the most extensively observed and widely studied fast blue optical transient to date; its unique observational properties challenge all existing standard models. In this paper, we model the luminosity evolution of the optical, soft X-ray, and hard X-ray emission, as well as the X-ray spectrum of AT2018cow with a magnetar-centered engine model. We consider a two-zone model with a striped magnetar wind in the interior and an expanding ejecta outside. The soft and hard X-ray emission of AT2018cow can be explained by the leakage of high-energy photons produced by internal gradual magnetic dissipation in the striped magnetar wind, while the luminous thermal UV/optical emission results from the thermalization of the ejecta by the captured photons. The two-component energy spectrum yielded by our model with a quasi-thermal component from the optically thick region of the wind superimposed on an optically thin synchrotron component well reproduces the X-ray spectral shape of AT2018cow. The Markov Chain Monte Carlo fitting results suggest that in order to explain the very short rise time to peak of the thermal optical emission, a low ejecta mass $M_{\rm ej}\approx0.1~M_\odot$ and high ejecta velocity $v_{\rm SN}\approx0.17c$ are required. A millisecond magnetar with $P_0\approx3.7~\rm ms$ and $B_p\approx2.4\times10^{14}~\rm G$ is needed to serve as the central engine of AT2018cow. △ Less

Submitted 22 February, 2024; originally announced February 2024.

Comments: 9 Pages, 5 Figures, 1 Tables, Accepted by ApJL; https://doi.org/10.3847/2041-8213/ad2611

A measurement of the sodium and iodine scintillation quenching factors across multiple NaI(Tl) detectors to identify systematics

Authors: D. Cintas, S. Hedges, W. G. Thompson, P. An, C. Awe, P. S. Barbeau, E. Barbosa de Souza, J. H. Jo, L. Li, M. Martínez, R. H. Maruyama, G. C. Rich, R. Runge, M. L. Sarsa

Abstract: The amount of light produced by nuclear recoils in scintillating targets is strongly quenched compared to that produced by electrons. A precise understanding of the quenching factor is particularly interesting for WIMP searches and CEνNS measurements since both rely on nuclear recoils, whereas energy calibrations are more readily accessible from electron recoils. There is a wide variation among th… ▽ More The amount of light produced by nuclear recoils in scintillating targets is strongly quenched compared to that produced by electrons. A precise understanding of the quenching factor is particularly interesting for WIMP searches and CEνNS measurements since both rely on nuclear recoils, whereas energy calibrations are more readily accessible from electron recoils. There is a wide variation among the current measurements of the quenching factor in sodium iodide (NaI) crystals, especially below 10 keV, the energy region of interest for dark matter and CEνNS studies. A better understanding of the quenching factor in NaI(Tl) is of particular interest for resolving the decades-old puzzle in the field of dark matter between the null results of most WIMP searches and the claim for dark matter detection by the DAMA/LIBRA collaboration. In this work, we measured sodium and iodine quenching factors for five small NaI(Tl) crystals grown with similar thallium concentrations and growth procedures. Unlike previous experiments, multiple crystals were tested, with measurements made in the same experimental setup to control systematic effects. The quenching factors agree in all crystals we investigated, and both sodium and iodine quenching factors are smaller than those reported by DAMA/LIBRA. The dominant systematic effect was due to the electron equivalent energy calibration originating from the non-proportional behavior of the NaI(Tl) light yield at lower energies, potentially the cause for the discrepancies among the previous measurements. △ Less

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

Comments: 22 pages, 27 figures, 7 tables. arXiv admin note: text overlap with arXiv:2310.07339

The Very Early Soft X-ray Plateau of GRB 230307A: Signature of an Evolving Radiative Efficiency in Magnetar Wind Dissipation?

Authors: Shu-Qing Zhong, Long Li, Di Xiao, Hui Sun, Bin-Bin Zhang, Zi-Gao Dai

Abstract: Very recently, a particularly long gamma-ray burst (GRB) 230307A was reported and proposed to originate from a compact binary merger based on its host galaxy property, kilonova, and heavy elements. More intriguingly, a very early plateau followed by a rapid decline in soft X-ray band was detected in its light curve by the Lobster Eye Imager for Astronomy, indicating strong evidence of the existenc… ▽ More Very recently, a particularly long gamma-ray burst (GRB) 230307A was reported and proposed to originate from a compact binary merger based on its host galaxy property, kilonova, and heavy elements. More intriguingly, a very early plateau followed by a rapid decline in soft X-ray band was detected in its light curve by the Lobster Eye Imager for Astronomy, indicating strong evidence of the existence of a magnetar as the merger product. This work explores that the Magnetar Wind Internal Gradual MAgnetic Dissipation (MIGMAD) model, in which the radiative efficiency evolves over time, successfully fits it to the observed data. Our results reinforce the notion that the X-ray plateau serves as a powerful indicator of a magnetar and imply that an evolving efficiency is likely to be a common feature in X-ray plateaus of GRB afterglows. In addition, we also discuss the explanations for the prompt emission, GRB afterglows, as well as kilonova, and predict possible kilonova afterglows in a magnetar central engine. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 10 pages, 4 figures, 2 tables. Accepted for publication in The Astrophysical Journal Letters

Quantifying Systematic Uncertainties in Experimental Physics: An Approximation Method

Authors: Lu Li

Abstract: In the domain of physics experiments, data fitting is a pivotal technique for extracting insights from both experimental and simulated datasets. This article presents an approximation method designed to estimate the systematic errors prevalent in data analyses. By applying our method to the Nab experiment, we compare our findings with simulation-derived results, thereby confirming the concordance… ▽ More In the domain of physics experiments, data fitting is a pivotal technique for extracting insights from both experimental and simulated datasets. This article presents an approximation method designed to estimate the systematic errors prevalent in data analyses. By applying our method to the Nab experiment, we compare our findings with simulation-derived results, thereby confirming the concordance of our approach with established simulation outcomes. This corroboration highlights the versatility of our method as a good tool for validating simulation results across various experimental contexts. △ Less

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

Neutron Star Mergers as the Dominant Contributor to the Production of Heavy $r$-Process Elements

Authors: Meng-Hua Chen, Li-Xin Li, Qiu-Hong Chen, Rui-Chong Hu, En-Wei Liang

Abstract: The discovery of the radioactively powered kilonova AT2017gfo, associated with the short-duration gamma-ray burst GRB 170817A and the gravitational wave source GW170817, has provided the first direct evidence supporting binary neutron star mergers as crucial astrophysical sites for the synthesis of heavy elements beyond iron through $r$-process nucleosysthesis in the universe. However, recent iden… ▽ More The discovery of the radioactively powered kilonova AT2017gfo, associated with the short-duration gamma-ray burst GRB 170817A and the gravitational wave source GW170817, has provided the first direct evidence supporting binary neutron star mergers as crucial astrophysical sites for the synthesis of heavy elements beyond iron through $r$-process nucleosysthesis in the universe. However, recent identifications of kilonovae following long-duration gamma-ray bursts, such as GRB 211211A and GRB 230307A, has sparked discussions about the potential of neutron star-white dwarf mergers to also produce neutron-rich ejecta and contribute to the production of heavy $r$-process elements. In this work, we estimate the contribution of binary neutron star mergers to the total mass of $r$-process elements in the Milky Way and investigate the possibility of neutron star-white dwarf mergers as alternative astrophysical sites for $r$-process nucleosynthesis through an analysis of the total mass of the $r$-process elements in the Milky Way. Our results reveal that binary neutron star mergers can sufficiently account for the Galactic heavy $r$-process elements, suggesting that these events are the dominant contributor to the production of heavy $r$-process elements in the Milky Way. Considering the total mass of $r$-process elements in the Milky Way and the higher occurrence rate of neutron star-white dwarf mergers, it is unlikely that such mergers can produce a significant amount of neutron-rich ejecta, with the generated mass of $r$-process elements being lower than $0.005M_{\odot}$. △ Less

Submitted 12 February, 2024; originally announced February 2024.

Comments: 7 pages, 4 figures, Accepted for publication in MNRAS

An Enduring Rapidly Moving Storm as a Guide to Saturn's Equatorial Jet's Complex Structure

Authors: A. Sánchez-Lavega, E. García-Melendo, S. Perez-Hoyos, R. Hueso, M. H. Wong, A. Simon, J. F. Sanz-Requena, A. Antuñano, N. Barrado-Izagirre, I. Garate-Lopez, J. F. Rojas, T. del Rio Gaztelurrutia, J. M. Gómez-Forrellad, I. de Pater, L. Li, T. Barry, PVOL contributors

Abstract: Saturn has an intense and broad eastward equatorial jet with a complex three-dimensional structure mixed with time variability. The equatorial region experiences strong seasonal insolation variations enhanced by ring shadowing and three of the six known giant planetary-scale storms have developed in it. These factors make Saturn's equator a natural laboratory to test models of jets in giant planet… ▽ More Saturn has an intense and broad eastward equatorial jet with a complex three-dimensional structure mixed with time variability. The equatorial region experiences strong seasonal insolation variations enhanced by ring shadowing and three of the six known giant planetary-scale storms have developed in it. These factors make Saturn's equator a natural laboratory to test models of jets in giant planets. Here we report on a bright equatorial atmospheric feature imaged in 2015 that moved steadily at a high speed of 450 ms-1 not measured since 1980-81 with other equatorial clouds moving within an ample range of velocities. Radiative transfer models show that these motions occur at three altitude levels within the upper haze and clouds. We find that the peak of the jet (latitudes 10\degree N to 10\degree S) suffers intense vertical shears reaching +2.5 ms-1 km-1, two orders of magnitude higher than meridional shears, and temporal variability above 1 bar altitude level. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Journal ref: Nature Communications, vol. 7, article number: 13626 (2016)

Formation Mechanism of Laser-Driven Magnetized "Pillars of Creation"

Authors: Zhu Lei, Lifeng Wang, Jiwei Li, Shiyang Zou, Junfeng Wu, Zhonghai Zhao, Wei Sun, Wenqiang Yuan, Longxing Li, Zheng Yan, Jun Li, Wenhua Ye, Xiantu He, Bin Qiao

Abstract: Pillars of Creation, one of the most recognized objects in the sky, are believed to be associated with the formation of young stars. However, so far, the formation and maintenance mechanism for the pillars are still not fully understood due to the complexity of the nonlinear radiation magneto-hydrodynamics (RMHD). Here, assuming laboratory laser-driven conditions, we studied the self-consistent dy… ▽ More Pillars of Creation, one of the most recognized objects in the sky, are believed to be associated with the formation of young stars. However, so far, the formation and maintenance mechanism for the pillars are still not fully understood due to the complexity of the nonlinear radiation magneto-hydrodynamics (RMHD). Here, assuming laboratory laser-driven conditions, we studied the self-consistent dynamics of pillar structures in magnetic fields by means of two-dimensional (2D) and three-dimensional (3D) RMHD simulations, and these results also support our proposed experimental scheme. We find only when the magnetic pressure and ablation pressure are comparable, the magnetic field can significantly alter the plasma hydrodynamics. For medium magnetized cases ($β_{initial} \approx 3.5$), {the initial magnetic fields undergo compression and amplification. This amplification results in the magnetic pressure inside the pillar becoming large enough to support the sides of the pillar against radial collapse due to pressure from the surrounding hot plasma. This effect is particularly pronounced for the parallel component ($B_y$), which is consistent with observational results.} In contrast, a strong perpendicular ($B_x, B_z$) magnetic field ($β_{initial} < 1$) almost remains its initial distribution and significantly suppresses the expansion of blow-off gas plasma, leading to the inability to form pillar-like structures. The 3D simulations suggest that the bending at the head of `Column \uppercase\expandafter{\romannumeral1}' in pillars of creation may be due to the non-parallel magnetic fields. After similarity scaling transformation, our results can be applied to explain the formation and maintenance mechanism of the pillars, and can also provide useful information for future experimental designs. △ Less

Submitted 30 January, 2024; originally announced January 2024.

The Activity of Black Hole Imprints on the Internal Plateau and the Subsequent Sharp Decay

Authors: Tian-Ci Zheng, Da-Ming Wei, Yun Wang, Hao Zhou, Long Li

Abstract: A stellar-mass black hole (BH) or a millisecond magnetar is believed to be born as the central engine of Gamma-ray bursts (GRBs). The presence of plateaus in the X-ray extended emission or afterglow of GRBs is widely accepted as an indicator of magnetar central engine, particularly those with a sharp decay (faster than $t^{-3}$), so-called internal plateau. However, an alternative model, by taking… ▽ More A stellar-mass black hole (BH) or a millisecond magnetar is believed to be born as the central engine of Gamma-ray bursts (GRBs). The presence of plateaus in the X-ray extended emission or afterglow of GRBs is widely accepted as an indicator of magnetar central engine, particularly those with a sharp decay (faster than $t^{-3}$), so-called internal plateau. However, an alternative model, by taking the evolution of the magnetic flux at the BH horizon into account, suggests that an internal plateau can also arise from a Blandford-Znajek (BZ) mechanism powered jet (hereafter referred to as the BZ jet). In this study, we propose that a precessional BZ jet would manifest a Quasi-Periodic Oscillation (QPO) signature on the internal plateau and the subsequent sharp decay. Such lightcurves cannot be readily explained by the activity of a short-lived, supermassive magnetar, thus favoring a Kerr BH as the central engine. The X-ray afterglow of GRB 050904, comprising nine flares, is characterized by a QPO-modulated plateau and sharp decay, which can be well reproduced by a precessional BZ jet model. Therefore, one potential clue for distinguishing between these two engines lies in whether QPO signature throughout the entire plateau and the subsequent sharp decay, as the magnetar scenario suggests a collapse at the end of the plateau. △ Less

Submitted 30 January, 2024; originally announced January 2024.

Comments: Accepted for publication in ApJ

A spectral data release for 104 Type II Supernovae from the Tsinghua Supernova Group

Authors: Han Lin, Xiaofeng Wang, Jujia Zhang, Danfeng Xiang, Tianmeng Zhang, Xulin Zhao, Xinghan Zhang, Hanna Sai, Liming Rui, Jun Mo, Gaobo Xi, Fang Huang, Xue Li, Yongzhi Cai, Weili Lin, Jie Lin, Chengyuan Wu, Jicheng Zhang, Zhihao Chen, Zhitong Li, Wenxiong Li, Linyi Li, Kaicheng Zhang, Cheng Miao, Juncheng Chen , et al. (11 additional authors not shown)

Abstract: We present 206 unpublished optical spectra of 104 type II supernovae obtained by the Xinglong 2.16m telescope and Lijiang 2.4m telescope during the period from 2011 to 2018, spanning the phases from about 1 to 200 days after the SN explosion. The spectral line identifications, evolution of line velocities and pseudo equivalent widths, as well as correlations between some important spectral paramet… ▽ More We present 206 unpublished optical spectra of 104 type II supernovae obtained by the Xinglong 2.16m telescope and Lijiang 2.4m telescope during the period from 2011 to 2018, spanning the phases from about 1 to 200 days after the SN explosion. The spectral line identifications, evolution of line velocities and pseudo equivalent widths, as well as correlations between some important spectral parameters are presented. Our sample displays a large range in expansion velocities. For instance, the Fe~{\sc ii} $5169$ velocities measured from spectra at $t\sim 50$ days after the explosion vary from ${\rm 2000\ km\ s^{-1}}$ to ${\rm 5500\ km\ s^{-1}}$, with an average value of ${\rm 3872 \pm 949\ km\ s^{-1}}$. Power-law functions can be used to fit the velocity evolution, with the power-law exponent quantifying the velocity decline rate. We found an anticorrelation existing between H$β$ velocity at mid-plateau phase and its velocity decay exponent, SNe II with higher velocities tending to have smaller velocity decay rate. Moreover, we noticed that the velocity decay rate inferred from the Balmer lines (i.e., H$α$ and H$β$) have moderate correlations with the ratio of absorption to emission for H$α$ (a/e). In our sample, two objects show possibly flash-ionized features at early phases. Besides, we noticed that multiple high-velocity components may exist on the blue side of hydrogen lines of SN 2013ab, possibly suggesting that these features arise from complex line forming region. All our spectra can be found in WISeREP and Zenodo. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Formation of millisecond pulsars with wide orbits

Authors: Bo Wang, Dongdong Liu, Yunlang Guo, Hailiang Chen, Wenshi Tang, Luhan Li, Zhanwen Han

Abstract: Millisecond pulsars (MSPs) are a kind of radio pulsars with short spin periods, playing a key role in many aspects of stellar astrophysics. In recent years, some more MSPs with wide orbits ($>30\,\rm d$) have been discovered, but their origin is still highly unclear. In the present work, according to an adiabatic power-law assumption for the mass-transfer process, we carried out a large number of… ▽ More Millisecond pulsars (MSPs) are a kind of radio pulsars with short spin periods, playing a key role in many aspects of stellar astrophysics. In recent years, some more MSPs with wide orbits ($>30\,\rm d$) have been discovered, but their origin is still highly unclear. In the present work, according to an adiabatic power-law assumption for the mass-transfer process, we carried out a large number of complete binary evolution computations for the formation of MSPs with wide orbits through the iron core-collapse supernova (CCSN) channel, in which a neutron star (NS) originating from a CCSN accretes matter from a red-giant (RG) star and spun up to millisecond periods. We found that this channel can form the observed MSPs with wide orbits in the range of $30-1200\,{\rm d}$, in which the WD companions have masses in the range of $0.28-0.55\,\rm M_{\odot}$. We also found that almost all the observed MSPs can be reproduced by this channel in the WD companion mass versus orbital period diagram. We estimate that the Galactic numbers of the resulting MSPs from the CCSN channel are in the range of $\sim 0.9-1.4\times10^{6}$. Compared with the accretion-induced collapse channel, the CCSN channel provides a dominant way to produce MSPs with wide orbits. △ Less

Submitted 9 January, 2024; originally announced January 2024.

Comments: 17 pages, 7 figures, submitted to MNRAS, a revised version after referee's comments

Constraints on photon mass and dark photon from the Jovian magnetic field

Authors: Shi Yan, Lingfeng Li, JiJi Fan

Abstract: The Jovian magnetic field, being the strongest and largest planetary one in the solar system, could offer us new insights into possible microscopic scale new physics, such as a non-zero mass of the Standard Model (SM) photon or a light dark photon kinetically mixing with the SM photon. We employ the immense data set from the latest Juno mission, which provides us unprecedented information about th… ▽ More The Jovian magnetic field, being the strongest and largest planetary one in the solar system, could offer us new insights into possible microscopic scale new physics, such as a non-zero mass of the Standard Model (SM) photon or a light dark photon kinetically mixing with the SM photon. We employ the immense data set from the latest Juno mission, which provides us unprecedented information about the magnetic field of the gas giant, together with a more rigorous statistical approach compared to the literature, to set strong constraints on the dark photon mass and kinetic mixing parameter, as well as the SM photon mass. The constraint on the dark photon parameters is independent of whether dark photon is (part of) dark matter or not, and serves as the most stringent one in a certain regime of the parameter space. △ Less

Submitted 11 December, 2023; originally announced December 2023.

Comments: 24 pages, 8 figures

Identification of Carbon Stars from LAMOST DR7

Authors: Linlin Li, Kecheng Zhang, Wenyuan Cui, Jianrong Shi, Wei Ji, Zhenyan Huo, Yawei Gao, Shuai Zhang, Mingxu Sun

Abstract: Carbon stars are excellent kinematic tracers of galaxies and play important roles in understanding the evolution of the Galaxy. Therefore, it is worthwhile to search for them in a large amount of spectra. In this work, we build a new carbon star catalog based on the LAMOST DR7 spectra. The catalog contains 4542 spectra of 3546 carbon stars, identified through line index and near-infrared color-col… ▽ More Carbon stars are excellent kinematic tracers of galaxies and play important roles in understanding the evolution of the Galaxy. Therefore, it is worthwhile to search for them in a large amount of spectra. In this work, we build a new carbon star catalog based on the LAMOST DR7 spectra. The catalog contains 4542 spectra of 3546 carbon stars, identified through line index and near-infrared color-color diagrams. Through visual inspection of the spectra, we further subclassify them into 925 C--H, 384 C--R, 608 C--N, and 1292 Ba stars. However, 437 stars could not be sub-classified due to their low signal-to-noise. Moreover, by comparing with LAMOST DR7 pipeline we find 567 more carbon stars and visually sub-classify them. We find that on the $J-H$ vs. $H-K_{\rm s}$ two-color diagram, C--N stars can be reliably distinguished from the other three sub-types. Additionally, by utilizing the Gaia distance, we study the distribution of carbon stars in the H-R diagram and identify 258 dwarf carbon stars by the criterion $M_{\rm G}>$5.0\,mag. Finally, we present the spatial distribution in Galactic coordinates of the 3546 carbon stars. The majority of C-N, C-R, and Ba stars are distributed at low Galactic latitudes, while most C--H and dC stars distribute at high Galactic latitudes. △ Less

Submitted 11 December, 2023; originally announced December 2023.

The measurement of masses of OB-type stars from LAMOST DR5

Authors: Zhenyan Huo, Zhicun Liu, Wenyuan Cui, Chao Liu, Jiaming Liu, Mingxu Sun, Shuai Feng, Linlin Li

Abstract: The measurements of masses and luminosities of massive stars play an important role in understanding the formation and evolution of their host galaxies. In this work, we present the measurement of masses and luminosities of 2,946 OB-type stars, including 78 O-type stars and 2,868 B-type stars, based on their stellar parameters (effective temperature, surface gravity, and metallicity) and PARSEC is… ▽ More The measurements of masses and luminosities of massive stars play an important role in understanding the formation and evolution of their host galaxies. In this work, we present the measurement of masses and luminosities of 2,946 OB-type stars, including 78 O-type stars and 2,868 B-type stars, based on their stellar parameters (effective temperature, surface gravity, and metallicity) and PARSEC isochrones model. Our results show that the median mass and luminosity of the 2,946 OB-type stars are 5.4 M$_{\odot}$ and log(L/L$_{\odot}$)=3.2 with the median relative error of 21.4$\%$ and 71.1$\%$, respectively. A good agreement between our results estimated by using our method and those derived by using the orbital motions of binary stars from the literature is found for some B-type stars. In addition, we also fit the mass-luminosity relation of B-type stars by using our derived mass and the luminosity from $Gaia$ DR3. △ Less

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

Comments: 12 pages, 10 figures, 1 table, accepted for publication in ApJS

Identification of Blue Horizontal-Branch Stars From LAMOST DR5

Authors: Jie Ju, Wenyuan Cui, Zhenyan Huo, Chao liu, Xiangxiang Xue, Jiaming Liu, Shuai Feng, Mingxu Sun, Linlin Li

Abstract: We construct a new catalog of the blue horizontal-branch (BHB) stars from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR5 dataset, which contains 5355+81 BHB stars at high Galactic latitude (($|Glat|>20^{\circ}$). We combine the spectral line indices with a set of Balmer line profile selection criteria to identify the BHB stars. During the selection process, we use the l… ▽ More We construct a new catalog of the blue horizontal-branch (BHB) stars from the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) DR5 dataset, which contains 5355+81 BHB stars at high Galactic latitude (($|Glat|>20^{\circ}$). We combine the spectral line indices with a set of Balmer line profile selection criteria to identify the BHB stars. During the selection process, we use the line index of \ion{Ca}{2}\,K to exclude the metal-rich A-type dwarfs. We obtain their atmospheric parameters by cross-matching our BHB stars with the catalog provided by \citet{Xiang2022}. The results show that our sample is consistent with the theoretical $T_{\rm eff}$-log\,$g$ evolutionary tracks of the BHB stars, indicating that our method is robust for identifying BHB stars from the LAMOST spectra. Their spatial distribution indicates that most of our BHB stars are located in the inner halo or the disk of the Milky Way. Combined with other BHB samples from the literature, the BHB stars can cover a large Galactic volume, which makes it a better probe for studying the kinematics, dynamics, and structural characteristics of the Milky Way. △ Less

Submitted 27 November, 2023; originally announced November 2023.

Comments: accepted by ApJS.15 pages, 18 figures

Data downloaded via parachute from a NASA super-pressure balloon

Authors: Ellen L. Sirks, Richard Massey, Ajay S. Gill, Jason Anderson, Steven J. Benton, Anthony M. Brown, Paul Clark, Joshua English, Spencer W. Everett, Aurelien A. Fraisse, Hugo Franco, John W. Hartley, David Harvey, Bradley Holder, Andrew Hunter, Eric M. Huff, Andrew Hynous, Mathilde Jauzac, William C. Jones, Nikky Joyce, Duncan Kennedy, David Lagattuta, Jason S. -Y. Leung, Lun Li, Stephen Lishman , et al. (18 additional authors not shown)

Abstract: In April to May 2023, the superBIT telescope was lifted to the Earth's stratosphere by a helium-filled super-pressure balloon, to acquire astronomical imaging from above (99.5% of) the Earth's atmosphere. It was launched from New Zealand then, for 40 days, circumnavigated the globe five times at a latitude 40 to 50 degrees South. Attached to the telescope were four 'DRS' (Data Recovery System) cap… ▽ More In April to May 2023, the superBIT telescope was lifted to the Earth's stratosphere by a helium-filled super-pressure balloon, to acquire astronomical imaging from above (99.5% of) the Earth's atmosphere. It was launched from New Zealand then, for 40 days, circumnavigated the globe five times at a latitude 40 to 50 degrees South. Attached to the telescope were four 'DRS' (Data Recovery System) capsules containing 5 TB solid state data storage, plus a GNSS receiver, Iridium transmitter, and parachute. Data from the telescope were copied to these, and two were dropped over Argentina. They drifted 61 km horizontally while they descended 32 km, but we predicted their descent vectors within 2.4 km: in this location, the discrepancy appears irreducible below 2 km because of high speed, gusty winds and local topography. The capsules then reported their own locations to within a few metres. We recovered the capsules and successfully retrieved all of superBIT's data - despite the telescope itself being later destroyed on landing. △ Less

Submitted 14 November, 2023; originally announced November 2023.

Comments: 12 pages

Journal ref: Aerospace 2023, 10, 960

Probing Thermal Electrons in GRB Afterglows

Authors: Hao-Xuan Gao, Jin-Jun Geng, Tian-Rui Sun, Liang Li, Yong-Feng Huang, Xue-Feng Wu

Abstract: Particle-in-cell simulations have unveiled that shock-accelerated electrons do not follow a pure power-law distribution, but have an additional low-energy "thermal" part, which owns a considerable portion of the total energy of electrons. Investigating the effects of these thermal electrons on gamma-ray burst (GRB) afterglows may provide valuable insights into the particle acceleration mechanisms.… ▽ More Particle-in-cell simulations have unveiled that shock-accelerated electrons do not follow a pure power-law distribution, but have an additional low-energy "thermal" part, which owns a considerable portion of the total energy of electrons. Investigating the effects of these thermal electrons on gamma-ray burst (GRB) afterglows may provide valuable insights into the particle acceleration mechanisms. We solve the continuity equation of electrons in the energy space, from which multi-wavelength afterglows are derived by incorporating processes including synchrotron radiation, synchrotron self-absorption, synchrotron self-Compton scattering, and gamma-gamma annihilation. First, there is an underlying positive correlation between temporal and spectral indices due to the cooling of electrons. Moreover, thermal electrons would result in the simultaneous non-monotonic variation in both spectral and temporal indices at multi-wavelength, which could be individually recorded by the 2.5-meter Wide Field Survey Telescope and Vera Rubin Observatory Legacy Survey of Space and Time (LSST). The thermal electrons could also be diagnosed from afterglow spectra by synergy observation in the optical (with LSST) and X-ray bands (with the Microchannel X-ray Telescope on board the Space Variable Objects Monitor). Finally, we use Monte Carlo simulations to obtain the distribution of peak flux ratio ($R_{\rm X}$) between soft and hard X-rays, and of the time delay ($Δt$) between peak times of soft X-ray and optical light curves. The thermal electrons significantly raise the upper limits of both $R_{\rm X}$ and $Δt$. Thus the distribution of GRB afterglows with thermal electrons is more dispersive in the $R_{\rm X} - Δt$ plane. △ Less

Submitted 10 June, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: 19 pages, 18 figures, accepted by the Astrophysical Journal

Partial Eruption of Solar Filaments. I. Configuration and Formation of Double-decker Filaments

Authors: Yijun Hou, Chuan Li, Ting Li, Jiangtao Su, Ye Qiu, Shuhong Yang, Liheng Yang, Leping Li, Yilin Guo, Zhengyong Hou, Qiao Song, Xianyong Bai, Guiping Zhou, Mingde Ding, Weiqun Gan, Yuanyong Deng

Abstract: Partial eruptions of solar filaments are the typical representative of solar eruptive behavior diversity. Here we investigate a typical filament partial eruption event and present integrated evidence for configuration of the pre-eruption filament and its formation. The CHASE H$α$ observations reveal structured Doppler velocity distribution within the pre-eruption filament, where distinct redshift… ▽ More Partial eruptions of solar filaments are the typical representative of solar eruptive behavior diversity. Here we investigate a typical filament partial eruption event and present integrated evidence for configuration of the pre-eruption filament and its formation. The CHASE H$α$ observations reveal structured Doppler velocity distribution within the pre-eruption filament, where distinct redshift only appeared in the east narrow part of the south filament region and then disappeared after the partial eruption while the north part dominated by blueshift remained. Combining the SDO, ASO-S observations, and NLFFF modeling results, we verify that there were two independent material flow systems within the pre-flare filament, whose magnetic topology is a special double-decker configuration consisting of two magnetic flux ropes (MFRs) with opposite magnetic twist. During the formation of this filament system, continuous magnetic flux cancellation and footpoint motion were observed around its north end. Therefore, we propose a new double-decker formation scenario that the two MFRs composing such double-decker configuration originated from two magnetic systems with different initial connections and opposite magnetic twist. Subsequent magnetic reconnection with surrounding newly-emerging fields resulted in the motion of footpoint of the upper MFR to the region around footpoint of the lower MFR, thus leading to eventual formation of the double-decker configuration consisting of two MFRs with similar footpoints but opposite signs of magnetic twist. These results provide a potential way to determine unambiguously the progenitor configuration of a partial-eruptive filament and reveal a special type of double-decker MFR configuration and a new double-decker formation scenario. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: 16 pages, 8 figures, 1 table, accepted for publication in ApJ as part of the Focus Issue "Early results from the Chinese Ha Solar Explorer (CHASE)"

Searching for Associations Between Short Gamma-ray Bursts and Fast Radio Burst

Authors: Ming-Xuan Lu, Long Li, Xiang-Gao Wang, Can-Min Deng, Yun-Feng Liang, Da-Bin Lin, En-Wei Liang

Abstract: The physical origin of fast radio bursts (FRBs) is still unclear. However, young magnetars associated with short-duration gamma-ray bursts (SGRBs) have been thought to be possible central engines for some FRBs. In this paper, we perform a systematic search for SGRBs that are associated with FRBs in a sample including 623 FRBs (601 one-off bursts and 22 repeaters) and 168 SGRBs with precise localiz… ▽ More The physical origin of fast radio bursts (FRBs) is still unclear. However, young magnetars associated with short-duration gamma-ray bursts (SGRBs) have been thought to be possible central engines for some FRBs. In this paper, we perform a systematic search for SGRBs that are associated with FRBs in a sample including 623 FRBs (601 one-off bursts and 22 repeaters) and 168 SGRBs with precise localizations. We find that FRB 190309A is spatially associated with GRB 060502B, with a chance probability of 0.05 when temporal and redshift information is taken into account. Considering the high chance probability (the statistical significance is < 3σ), we examine other observational properties such as the host galaxy, the dispersion measure, and the energy budget of the central engine to check the possibility of their association. Although the available observational information is insufficient to determine whether they are physically associated, it does not rule out such a possibility. As the only pair of FRB and GRB that are spatially associated, it remains an interesting case worthy of further attention △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 19 pages, 4 figures, accepted to ApJ

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

Ages and metallicities of stellar clusters using S-PLUS narrow-band integrated photometry: the Small Magellanic Cloud

Authors: Gabriel Fabiano de Souza, Pieter Westera, Felipe Almeida-Fernandes, Guilherme Limberg, Bruno Dias, José A. Hernandez-Jimenez, Fábio R. Herpich, Leandro O. Kerber, Eduardo Machado-Pereira, Hélio D. Perottoni, Rafael Guerço, Liana Li, Laura Sampedro, Antonio Kanaan, Tiago Ribeiro, William Schoenell, Claudia Mendes de Oliveira

Abstract: The Magellanic Clouds are the most massive and closest satellite galaxies of the Milky Way, with stars covering ages from a few Myr up to 13 Gyr. This makes them important for validating integrated light methods to study stellar populations and star-formation processes, which can be applied to more distant galaxies. We characterized a set of stellar clusters in the Small Magellanic Cloud (SMC), us… ▽ More The Magellanic Clouds are the most massive and closest satellite galaxies of the Milky Way, with stars covering ages from a few Myr up to 13 Gyr. This makes them important for validating integrated light methods to study stellar populations and star-formation processes, which can be applied to more distant galaxies. We characterized a set of stellar clusters in the Small Magellanic Cloud (SMC), using the $\textit{Southern Photometric Local Universe Survey}$. This is the first age (metallicity) determination for 11 (65) clusters of this sample. Through its 7 narrow bands, centered on important spectral features, and 5 broad bands, we can retrieve detailed information about stellar populations. We obtained ages and metallicities for all stellar clusters using the Bayesian spectral energy distribution fitting code $\texttt{BAGPIPES}$. With a sample of clusters in the color range $-0.20 < r-z < +0.35$, for which our determined parameters are most reliable, we modeled the age-metallicity relation of SMC. At any given age, the metallicities of SMC clusters are lower than those of both the Gaia Sausage-Enceladus disrupted dwarf galaxy and the Milky Way. In comparison with literature values, differences are $Δ$log(age)$\approx0.31$ and $Δ$[Fe/H]$\approx0.41$, which is comparable to low-resolution spectroscopy of individual stars. Finally, we confirm a previously known gradient, with younger clusters in the center and older ones preferentially located in the outermost regions. On the other hand, we found no evidence of a significant metallicity gradient. △ Less

Submitted 30 November, 2023; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: 12 pages, 11 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 2, January 2024, Pages 1733-1744