Showing 1–50 of 204 results for author: Lu, H

The Preliminary Mauve Science Programme: Science themes identified for the first year of operations

Authors: Mauve Science Collaboration, Marcel Agueros, Don Dixon, Chuanfei Dong, Girish M. Duvvuri, Patrick Flanagan, Christopher Johns-Krull, Hongpeng Lu, Hiroyuki Maehara, Kosuke Namekata, Alejandro Nunez, Elena Pancino, Sharmila Rani, Anusha Ravikumar, T. A. A. Sigut, Keivan Stassun, Jamie Stewart, Krisztián Vida, Emma Whelan, Benjamin Wilcock, Sharafina Razin, Arianna Saba, Giovanna Tinetti, Marcell Tessenyi, Jonathan Tennyson

Abstract: Mauve is a low-cost small satellite developed and operated by Blue Skies Space Ltd. The payload features a 13 cm telescope connected with a fibre that feeds into a UV-Vis spectrometer. The detector covers the 200-700 nm range in a single shot, obtaining low resolution spectra at R~20-65. Mauve has launched on 28th November 2025, reaching a 510 km Low-Earth Sun-synchronous orbit. The satellite will… ▽ More Mauve is a low-cost small satellite developed and operated by Blue Skies Space Ltd. The payload features a 13 cm telescope connected with a fibre that feeds into a UV-Vis spectrometer. The detector covers the 200-700 nm range in a single shot, obtaining low resolution spectra at R~20-65. Mauve has launched on 28th November 2025, reaching a 510 km Low-Earth Sun-synchronous orbit. The satellite will enable UV and visible observations of a variety of stellar objects in our Galaxy, filling the gaps in the ultraviolet space-based data. The researchers that have already joined the mission have defined the science themes, observational strategy and targets that Mauve will observe in the first year of operations. To date, 10 science themes have been developed by the Mauve science collaboration for year 1, with observational strategies that include both long duration monitoring and short cadence snapshots. Here, we describe these themes and the science that Mauve will undertake in its first year of operations. △ Less

Submitted 18 December, 2025; originally announced December 2025.

Comments: 16 pages, 15 figures, 2 tables

A Machine-Learning Approach for Identifying CME-Associated Stellar Flares in TESS Observations

Authors: Yu Shi, Hong-Peng Lu, Li-Yun Zhang, Tian-Hao Su, Chao Tan

Abstract: Coronal mass ejections (CMEs) are major drivers of stellar space weather and can strongly influence the habitability of exoplanets. However, compared to the frequent occurrence of white-light flares, confirmed stellar CMEs remain extremely rare. Whether such flares are commonly accompanied by CMEs is a key question for solar-stellar comparative studies. Using Sun-as-a-star soft X-ray flare light c… ▽ More Coronal mass ejections (CMEs) are major drivers of stellar space weather and can strongly influence the habitability of exoplanets. However, compared to the frequent occurrence of white-light flares, confirmed stellar CMEs remain extremely rare. Whether such flares are commonly accompanied by CMEs is a key question for solar-stellar comparative studies. Using Sun-as-a-star soft X-ray flare light curves observed by the GOES XRS 1--8~Å channel, we compiled a sample of 1,766 M-class and larger solar flares and extracted features with both deep convolutional neural networks and manual methods. Five machine-learning classifiers were trained to distinguish eruptive from confined flares, with the random forest model achieving the best performance (true skill statistic; TSS = 0.31). This TSS value indicates that the model possesses a moderate ability to discriminate between eruptive and confined flares. Normalized white-light and GOES XRS flare light curves show broadly consistent temporal evolution, reflecting their shared energy-release history and supporting a probabilistic transfer of the model to white-light flare data. We applied the best-performing RF model to 41,405 TESS-detected flares on FGKM-type main-sequence stars, predicting that approximately 47% of events show CME-like morphological characteristics, with the model-implied intrinsic association fraction lying in the range 35%--60%. Intriguingly, the CME occurrence rate decreases with increasing flare energy, indicating that the most energetic flares may be more strongly confined by overlying magnetic fields. These results provide new insight into flare-CME connections in diverse stellar environments and have important implications for assessing the impact of stellar eruptive activity on exoplanetary atmospheres. △ Less

Submitted 17 December, 2025; originally announced December 2025.

Comments: Accepted by ApJS

Estimating stellar atmospheric parameters and elemental abundances using fully connected residual network

Authors: Shuo Li, Yin-Bi Li, A-Li Luo, Jun-Chao Liang, Hai-Ling Lu, Hugh R. A. Jones

Abstract: Stellar atmospheric parameters and elemental abundances are traditionally determined using template matching techniques based on high-resolution spectra. However, these methods are sensitive to noise and unsuitable for ultra-low-resolution data. Given that the Chinese Space Station Telescope (CSST) will acquire large volumes of ultra-low-resolution spectra, developing effective methods for ultra-l… ▽ More Stellar atmospheric parameters and elemental abundances are traditionally determined using template matching techniques based on high-resolution spectra. However, these methods are sensitive to noise and unsuitable for ultra-low-resolution data. Given that the Chinese Space Station Telescope (CSST) will acquire large volumes of ultra-low-resolution spectra, developing effective methods for ultra-low-resolution spectral analysis is crucial. In this work, we investigated the Fully Connected Residual Network (FCResNet) for simultaneously estimating atmospheric parameters ($T_\text{eff}$, $\log g$, [Fe/H]) and elemental abundances ([C/Fe], [N/Fe], [Mg/Fe]). We trained and evaluated FCResNet using CSST-like spectra (\textit{R} $\sim$ 200) generated by degrading LAMOST spectra (\textit{R} $\sim$ 1,800), with reference labels from APOGEE. FCResNet significantly outperforms traditional machine learning methods (KNN, XGBoost, SVR) and CNN in prediction precision. For spectra with g-band signal-to-noise ratio greater than 20, FCResNet achieves precisions of 78 K, 0.15 dex, 0.08 dex, 0.05 dex, 0.10 dex, and 0.05 dex for $T_\text{eff}$, $\log g$, [Fe/H], [C/Fe], [N/Fe] and [Mg/Fe], respectively, on the test set. FCResNet processes one million spectra in only 42 seconds while maintaining a simple architecture with just 348 KB model size. These results suggest that FCResNet is a practical and promising tool for processing the large volume of ultra-low-resolution spectra that will be obtained by CSST in the future. △ Less

Submitted 11 December, 2025; originally announced December 2025.

Journal ref: Res. Astron. Astrophys. 25 125009 (2025)

Refined M-type Star Catalog from LAMOST DR10: Measurements of Radial Velocities, $T_\text{eff}$, log $g$, [M/H] and [$α$/M]

Authors: Shuo Li, Yin-Bi Li, A-Li Luo, Jun-Chao Liang, You-Fen Wang, Jing Chen, Shuo Zhang, Mao-Sheng Xiang, Hugh R. A. Jones, Zhong-Rui Bai, Xiao-Xiao Ma, Yun-Jin Zhang, Hai-Ling Lu

Abstract: Precise stellar parameters for M-type stars, the Galaxy's most common stellar type, are crucial for numerous studies. In this work, we refined the LAMOST DR10 M-type star catalog through a two-stage process. First, we purified the catalog using techniques including deep learning and color-magnitude diagrams to remove 22,496 non-M spectra, correct 2,078 dwarf/giant classifications, and update 12,90… ▽ More Precise stellar parameters for M-type stars, the Galaxy's most common stellar type, are crucial for numerous studies. In this work, we refined the LAMOST DR10 M-type star catalog through a two-stage process. First, we purified the catalog using techniques including deep learning and color-magnitude diagrams to remove 22,496 non-M spectra, correct 2,078 dwarf/giant classifications, and update 12,900 radial velocities. This resulted in a cleaner catalog containing 870,518 M-type spectra (820,493 dwarfs, 50,025 giants). Second, applying a label transfer strategy using values from APOGEE DR16 for parameter prediction with a ten-fold cross-validated CNN ensemble architecture, we predicted $T_\text{eff}$, $\log g$, [M/H], and [$α$/M] separately for M dwarfs and giants. The average internal errors for M dwarfs/giants are respectively: $T_\text{eff}$ 30/17 K, log $g$ 0.07/0.07 dex, [M/H] 0.07/0.05 dex, and [$α$/M] 0.02/0.02 dex. Comparison with APOGEE demonstrates external precisions of 34/14 K, 0.12/0.07 dex, 0.09/0.04 dex, and 0.03/0.02 dex for M dwarfs/giants, which represents precision improvements of over 20\% for M dwarfs and over 50\% for M giants compared to previous literature results. The catalog is available at https://nadc.china-vo.org/res/r101668/. △ Less

Submitted 11 December, 2025; originally announced December 2025.

Journal ref: ApJS 281 58 (2025)

Probing the Nature of High-Redshift Long GRB 250114A and Its Magnetar Central Engine

Authors: Wen-Yuan Yu, Hou-Jun Lü, Xiao Tian, Liang-Jun Chen, En-Wei Liang

Abstract: GRB 250114A is a long-duration gamma-ray burst (GRB) which triggered the Swift/BAT with a spectroscopic high-redshift at $z = 4.732$. The light curve of the prompt emission is composed of three distinct emission episodes, which are separated by quiescent gaps ranging from tens to hundreds of seconds. While the X-ray light curve exhibits the canonical X-ray emission which is composed of several pow… ▽ More GRB 250114A is a long-duration gamma-ray burst (GRB) which triggered the Swift/BAT with a spectroscopic high-redshift at $z = 4.732$. The light curve of the prompt emission is composed of three distinct emission episodes, which are separated by quiescent gaps ranging from tens to hundreds of seconds. While the X-ray light curve exhibits the canonical X-ray emission which is composed of several power-law segments superposition of a giant X-ray flare. More interestingly, there is still significant X-ray emission during the quiescent time in the prompt emission, suggesting a continuously active central engine whose power fluctuates across the $γ$-ray detectability threshold. In this paper, we propose a magnetar as the central engine of GRB 250114A by fitting the X-ray light curve, and infer a magnetic field strength $B_{\rm p}=13.24^{+1.73}_{-5.84} \, \times10^{15}\ \mathrm{G}$ and an initial spin period $P_{0}=14.31^{+0.93}_{-3.16} \, \mathrm{ms}$ of magnetar, with a jet correction, fall within a reasonable range. Furthermore, we also compare the prompt emission, X-ray afterglow, $E_{\mathrm p}$-$E_{γ,\mathrm{iso}}$, and $\varepsilon-$distribution of GBR 250114A with those of other high-$z$ sample-GRBs, and find no significant statistical differences between them. △ Less

Submitted 26 November, 2025; originally announced November 2025.

Comments: 16 pages, 7 figures, and 3 tables. Accepted for publication in ApJ

A Method for Gamma-Ray Energy Spectrum Inversion and Correction

Authors: Zhi-Qiang Ding, Xin-Qiao Li, Da-Li Zhang, Zheng-Hua An, Zhen-Xia Zhang, Roberto Battiston, Roberto Iuppa, Zhuo Li, Yan-Qiu Zhang, Yan Huang, Chao Zheng, Yan-Bing Xu, Xiao-Yun Zhao, Lu Wang, Ping Wang, Hong Lu

Abstract: Accurate spectral analysis of high-energy astrophysical sources often relies on comparing observed data to incident spectral models convolved with the instrument response. However, for Gamma-Ray Bursts and other high-energy transient events observed at high count rates, significant distortions (e.g., pile-up, dead time, and large signal trailing) are introduced, complicating this analysis. We pres… ▽ More Accurate spectral analysis of high-energy astrophysical sources often relies on comparing observed data to incident spectral models convolved with the instrument response. However, for Gamma-Ray Bursts and other high-energy transient events observed at high count rates, significant distortions (e.g., pile-up, dead time, and large signal trailing) are introduced, complicating this analysis. We present a method framework to address the model dependence problem, especially to solve the problem of energy spectrum distortion caused by instrument signal pile-up due to high counting rate and high-rate effects, applicable to X-ray, gamma-ray, and particle detectors. Our approach combines physics-based Monte Carlo (MC) simulations with a model-independent spectral inversion technique. The MC simulations quantify instrumental effects and enable correction of the distorted spectrum. Subsequently, the inversion step reconstructs the incident spectrum using an inverse response matrix approach, conceptually equivalent to deconvolving the detector response. The inversion employs a Convolutional Neural Network, selected for its numerical stability and effective handling of complex detector responses. Validation using simulations across diverse input spectra demonstrates high fidelity. Specifically, for 27 different parameter sets of the brightest gamma-ray bursts, goodness-of-fit tests confirm the reconstructed spectra are in excellent statistical agreement with the input spectra, and residuals are typically within $\pm 2σ$. This method enables precise analysis of intense transients and other high-flux events, overcoming limitations imposed by instrumental effects in traditional analyses. △ Less

Submitted 19 November, 2025; originally announced November 2025.

Comments: The Astrophysical Journal has accepted

Measuring the magnetic fields in the chromospheres of low-mass stars

Authors: Tianqi Cang, Pascal Petit, Jean-François Donati, Hui Tian, Jianning Fu, Hao Li, Stefano Bellotti, Xueying Hu, Xiaoyu Ma, Arturo Lopez Ariste, Keyu Xing, Julien Morin, Hongpeng Lu, Weikai Zong

Abstract: Magnetic fields in the upper atmospheres of solar-like stars are believed to provide an enormous amount of energy to power the hot coronae and drive large-scale eruptions that could impact the habitability of planetary systems around these stars. However, these magnetic fields have never been routinely measured on stars beyond the solar system. Through decade-long spectropolarimetric observations,… ▽ More Magnetic fields in the upper atmospheres of solar-like stars are believed to provide an enormous amount of energy to power the hot coronae and drive large-scale eruptions that could impact the habitability of planetary systems around these stars. However, these magnetic fields have never been routinely measured on stars beyond the solar system. Through decade-long spectropolarimetric observations, we have now achieved the measurements of magnetic fields in the lower and middle chromospheres of three M-dwarfs. Our results indicate that the line-of-sight component of the chromospheric magnetic fields can reach up to hundreds of Gauss, whose sign frequently opposes that of the photospheric field. The measurements highlight the magnetic field complexity and the variation with height close to the surface of these M-dwarfs. They provide critical constraints on the energy budget responsible for heating and eruptions of stellar upper atmospheres, and enable assessments of how stellar magnetic activity may affect exoplanet environments. △ Less

Submitted 11 November, 2025; originally announced November 2025.

Comments: 12 figs, 3 tabs, accepted by Nature Communications

Constraint on the Physical Origin of Gamma-Ray Burst Prompt Emission via Its Nondetected Diffuse Neutrino Emission

Authors: Yang-Dong-Jun Ou, Hou-Jun Lü, Jia-Ming Chen, Ben-Yang Zhu, En-Wei Liang

Abstract: The physical origin of prompt emission in gamma-ray bursts (GRBs) remains an open question since it has been studied more than half a century. Three alternative models (i.e. dissipative photosphere, internal shock, and Internal-Collision-induced MAgnetic Reconnection and Turbulence, ICMART) have been proposed to interpret the observations of GRB prompt emission, but none of them can fully interpre… ▽ More The physical origin of prompt emission in gamma-ray bursts (GRBs) remains an open question since it has been studied more than half a century. Three alternative models (i.e. dissipative photosphere, internal shock, and Internal-Collision-induced MAgnetic Reconnection and Turbulence, ICMART) have been proposed to interpret the observations of GRB prompt emission, but none of them can fully interpret all of the observational data collected so far. The question is what is the fraction of these three theoretical models in the prompt emission of GRBs. In this paper, we propose to utilize an innovative method and constrain the fraction of GRB prompt emission models via its nondetected diffuse neutrinos. By adopting two methods (e.g., summing up the individual GRB contributions and assumed luminosity functions of GRB) to calculate diffuse neutrino flux of GRBs for given the benchmark parameters of $Γ=300$ and $\varepsilon_{p} \text{/} \varepsilon_{e}=10$, both approaches indicate that most GRBs should be originated from the ICMART model. Moreover, we find that the fractions of the dissipative photosphere model, the internal shock model, and the ICMART model are constrained to be [0, $0.5\%$], [0, $1.1\%$], and [$98.9\%$, 1], respectively, for the method of summing up the individual GRB contributions. For the method of luminosity functions, the fractions of above three models are constrained to be [0, $6.1\%$], [0, $8.2\%$], and [$91.8\%$, 1], respectively. However, such fractions of different models are also dependent on the parameters of $Γ$ and $\varepsilon_{p} \text{/} \varepsilon_{e}$. △ Less

Submitted 26 November, 2025; v1 submitted 26 October, 2025; originally announced October 2025.

Comments: 9 pages, 4 figures, and 2 tables, matched with the published verison

Journal ref: ApJ,994:91 (2025)

Starspots as the origin of ultrafast drifting radio bursts from an active M dwarf

Authors: Jiale Zhang, Hui Tian, Stefano Bellotti, Tianqi Cang, Joseph R. Callingham, Harish K. Vedantham, Bin Chen, Sijie Yu, Philippe Zarka, Corentin K. Louis, Peng Jiang, Hongpeng Lu, Yang Gao, Jinghai Sun, Hengqian Gan, Hui Li, Chun Sun, Zheng Lei, Menglin Huang

Abstract: Detecting coherent radio bursts from nearby M dwarfs provides opportunities for exploring their magnetic activity and interaction with orbiting exoplanets. However, it remains uncertain if the emission is related to flare-like activity similar to the Sun or magnetospheric process akin to magnetized planets. Using observations (1.0 - 1.5 GHz) taken by the Five-hundred-meter Aperture Spherical radio… ▽ More Detecting coherent radio bursts from nearby M dwarfs provides opportunities for exploring their magnetic activity and interaction with orbiting exoplanets. However, it remains uncertain if the emission is related to flare-like activity similar to the Sun or magnetospheric process akin to magnetized planets. Using observations (1.0 - 1.5 GHz) taken by the Five-hundred-meter Aperture Spherical radio Telescope, we found a type of millisecond-scale radio bursts with exceptionally high frequency drift rates ($\sim 8\;\rm{GHz\;s^{-1}}$) from an active M dwarf, AD Leo. The ultrafast drift rates point to a source region with a notably low magnetic scale height ($<0.15\; r_\star$, $r_\star$ as the stellar radius), a feature not expected in a commonly assumed dipole-like global field but highly possible in localized strong-field structures, i.e. starspots. Our findings suggest that a concentrated magnetic field above starspots could be responsible for some of the most intense radio bursts from M dwarfs, supporting a solar-like electron acceleration mechanism. △ Less

Submitted 20 October, 2025; originally announced October 2025.

Comments: Published on Science Advances, the authors' version

Journal ref: Science Advances (2025) 11, 42

The color evolution of magnetar-powered kilonova emission in merging neutron star-neutron star systems

Authors: Suo-Ning Wang, Hou-Jun Lü, Xiao-Xuan Liu, Jared Rice, Jia Ren, En-Wei Liang

Abstract: The first direct detection of the gravitational wave (GW) event GW170817 and its electromagnetic (EM) counterpart open a new window for studying of multi-messenger astronomy. However, how to identify the remnant of binary neutron star (NS) merger via EM radiation remain an open question. In this paper, we propose a method of color evolution of kilonova emission to identify its progenitors. We assu… ▽ More The first direct detection of the gravitational wave (GW) event GW170817 and its electromagnetic (EM) counterpart open a new window for studying of multi-messenger astronomy. However, how to identify the remnant of binary neutron star (NS) merger via EM radiation remain an open question. In this paper, we propose a method of color evolution of kilonova emission to identify its progenitors. We assume that the energy of the kilonova is contributed from radioactive decay, magnetar spin-down, and pulsar wind nebula (PWN). The color evolution of kilonova emission associated with short GRB is significant when the spectrum is thermal emission, while it tends towards a constant when the spectrum is non-thermal radiation. On the other hand, if the central engine is a black hole (BH) which is promptly generated by the NS-NS merger or NS-BH merger, then the kilonova is powered only by the radioactive decay. There is no color evolution at the beginning before the peak of kilonova emission, but is significantly and rapidly increasing after the peak. On the contrary, if the central engine is a magnetar or stable NS, the kilonova emission is contributed from radioactive decay, magnetar, and PWN. The color evolution after the peak of kilonova emission is complex behavior which depends on the rotational energy and spin-down time-scale of magnetar, and finally tend to a constant in the late state. △ Less

Submitted 21 August, 2025; originally announced August 2025.

Comments: 11 pages and 4 figure. Accepted for publication in MNRAS

Signature of a magnetar central engine with precession motion in the X-ray emission of GRB 220711B

Authors: Ying-Ze Shan, Xing Yang, Hou-Jun Lü, Jared Rice, Hao-Yu Yuan, Xue-Zhao Chang, Zhao Joseph Zhang, Le Zou, En-Wei Liang

Abstract: The $γ$-ray light curve of long-duration GRB 220711B, is characterized by a multi-peaked structure with a duration lasting $\sim$105 seconds. More interestingly, the X-ray afterglow light curve is composed of a plateau emission smoothly connected with a $\sim t^{-2}$ segment overlapping some flares followed by an extremely steep decay. By analysing the light curves of both prompt emission and X-ra… ▽ More The $γ$-ray light curve of long-duration GRB 220711B, is characterized by a multi-peaked structure with a duration lasting $\sim$105 seconds. More interestingly, the X-ray afterglow light curve is composed of a plateau emission smoothly connected with a $\sim t^{-2}$ segment overlapping some flares followed by an extremely steep decay. By analysing the light curves of both prompt emission and X-ray afterglow, no high-confidence-level quasi-periodic oscillation (QPO) signals are found in the light curves of the prompt emission (e.g., BAT and GBM), but it is found that a QPO signal at $\sim$ 50 s above 6$σ$ confidence level indeed exist in the X-ray afterglow. Here, we propose that a supra-massive magnetar as the central engine of GRB 220711B with precession motion is a good interpretation of the features of the X-ray emission. The initial plateau emission and followed decay segment, as well as the extremely steep-decay segment, are consistent with the physical process of supra-massive magnetar spin-down and then collapse into black hole. Moreover, the QPO signal in the X-ray emission can be explained as an effect of the precession motion of the magnetar. If this is the case, one can derive various magnetar parameters such as the initial period ($P_{\rm{0}}$) and surface magnetic field strength ($B_{\rm{p}}$) within a pseudo-redshift range of [1.08, 4.27]. By considering beaming corrections with jet opening angle $5^{\circ}$, we find that $P_{\rm{0}}$ and $B_{\rm{p}}$ lie within the range of [1.87, 6.25] ms and [$1.47\times 10^{16}$, $3.09\times 10^{16}$] G, respectively. The parameter of $B_{\rm{p}}$ is slightly larger than that of other typical long-duration GRBs, but $P_{\rm{0}}$ fall in a reasonable range. △ Less

Submitted 31 July, 2025; originally announced August 2025.

Comments: 14 pages, 7 figures, and 2 Tables, Accepted for publication in MNRAS

Neutron Star-White Dwarf Merger as One Possible Optional Source of Kilonova-like Emission: Implications for GRB 211211A

Authors: Xiao-Xuan Liu, Hou-Jun Lü, Qiu-Hong Chen, Zhao-Wei Du, En-Wei Liang

Abstract: Long-duration GRB 211211A, which lacked an associated supernova at such a low redshift $z=0.076$, but was associated with a possible kilonova emission, has attracted great attention. The neutron star-white dwarf (NS-WD) merger is proposed as a possible progenitor of GRB 211211A, and it could naturally explain the long duration of the prompt emission. However, the NS-WD merger is not an ideal site… ▽ More Long-duration GRB 211211A, which lacked an associated supernova at such a low redshift $z=0.076$, but was associated with a possible kilonova emission, has attracted great attention. The neutron star-white dwarf (NS-WD) merger is proposed as a possible progenitor of GRB 211211A, and it could naturally explain the long duration of the prompt emission. However, the NS-WD merger is not an ideal site for producing heavy elements via r-process nucleosynthesis. In this Letter, we investigate the heavy elements produced in NS-WD mergers based on numerical simulations of nucleosynthesis via SkyNet, and then calculate the resulting kilonova-like emission to compare with the solidly observed case of possible kilonova emission associated with GRB 211211A. By adopting three models (i.e., Model-A, Model-B, and Model-C) from \cite{2023ApJ...956...71K} at different temperatures ($T=4$ GK, 5 GK, and 6 GK), which are treated as free parameters, we find that the mass number of the heaviest element produced in our simulations is less than 90 ($A< 90$). Moreover, by comparing the calculated kilonova-like emission with the afterglow-subtracted observations of the possible kilonova associated with GRB 211211A, it is found that the merger of NS and WD cannot be ruled out as the origin of GRB 211211A to produce the possible kilonova emission if the remnant of the WD-NS merger is a supramassive or stable magnetar. Otherwise, it is difficult to explain the early possible kilonova emission following GRB 211211A by invoking the merger of a WD and an NS. △ Less

Submitted 24 July, 2025; v1 submitted 6 July, 2025; originally announced July 2025.

Comments: 11 pages, 3 figures,and 1 Table, accepted for publication in ApJL,and matched with the published verison

Journal ref: ApJL, 988:L46, 2025

How to identify the object with mass range of $(2.2-3)M_\odot$ in the merger of compact star systems

Authors: ZhaoWei Du, HouJun Lü, Xiaoxuan Liu, XiLong Fan, EnWei Liang

Abstract: High-frequency gravitational-wave (GW) radiation has been detected by LIGO-Virgo-KAGRA in the merger of compact stars. However, two GW events, GW190814 and GW200210, the mass of one companion object falls into the mass region of $(2.2-3)\rm~M_\odot$, and how to identify such object (e.g., as a low-mass black hole (BH) or a massive neutron star (NS)) remains an open question. In this paper, we prop… ▽ More High-frequency gravitational-wave (GW) radiation has been detected by LIGO-Virgo-KAGRA in the merger of compact stars. However, two GW events, GW190814 and GW200210, the mass of one companion object falls into the mass region of $(2.2-3)\rm~M_\odot$, and how to identify such object (e.g., as a low-mass black hole (BH) or a massive neutron star (NS)) remains an open question. In this paper, we propose a method to identify the mystery compact object (MCO) with the mass region of $(2.2-3)\rm~M_\odot$ in a binary system via the possible electromagnetic (EM) radiations before and after the mergers. A multi-band EM emission can be produced with $L\propto(-t)^{7/4}$ (or $L\propto(-t)^{-5/4}$) during the inspiral phase due to the BH battery (or interaction magnetospheres) mechanism, and a bright (or dark) kilonova emission is powered by radioactive decay with ejecta mass ratio $q>1.7$ (or $q<1.7$) during the post-merge state when MCO is as a low-mass BH (or massive NS) to merger with NS. Moreover, by considering the merger system between MCO and a BH when MCO is a massive NS, we find that it requires the BH with high spin (e.g., $a\sim0.8-0.99$) to make sure the tidal disruption event (TDE) occurred, and a multi-band precursor emission and bright kilonova emission can also be produced during the inspiral phase and post-merge state, respectively. In any case, no matter which mechanism we adopt, such precursor emissions are too weak to be detected by most current telescopes unless the distance is close enough. △ Less

Submitted 10 April, 2025; originally announced April 2025.

Comments: 16 pages, 5 figures,and 1 table, accepted for publication in ApJ

Signature of strange star as the central engine of GRB 240529A

Authors: Xiao Tian, HouJun Lü, WenJun Tan, ShaoLin Xiong, HaoYu Yuan, WenYuan Yu, ShuQing Zhong, WenLong Zhang, EnWei Liang

Abstract: GRB 240529A is a long-duration gamma-ray burst (GRB) whose light curve of prompt emission is composed of a triple-episode structure, separated by quiescent gaps of tens to hundreds of seconds. More interestingly, its X-ray light curve of afterglow exhibits two-plateau emissions, namely, an internal plateau emission that is smoothly connected with a $\sim t^{-0.1}$ segment and followed by a… ▽ More GRB 240529A is a long-duration gamma-ray burst (GRB) whose light curve of prompt emission is composed of a triple-episode structure, separated by quiescent gaps of tens to hundreds of seconds. More interestingly, its X-ray light curve of afterglow exhibits two-plateau emissions, namely, an internal plateau emission that is smoothly connected with a $\sim t^{-0.1}$ segment and followed by a $\sim t^{-2}$ power-law decay. The three episodes in the prompt emission, together with two plateau emissions in X-ray, are unique in the Swift era. They are very difficult to explain with the standard internal/external shock model by invoking a black hole central engine. However, it could be consistent with the prediction of a supramassive magnetar as the central engine, the physical process of phase transition from magnetar to strange star, as well as the cooling and spin-down of the strange star. In this paper, we propose that the first- and second-episode emissions in the prompt $γ-$ray of GRB 240529A are from the jet emission of a massive star collapsing into a supramassive magnetar and the re-activity of central engine, respectively. Then, the third-episode emission of prompt is attributed to the phase transition from a magnetar to a strange star. Finally, the first- and second-plateau emissions of the X-ray afterglow are powered by the cooling and spin-down of the strange star, respectively. The observational data of each component of GRB 240529A are roughly coincident with the estimations of the above physical picture. △ Less

Submitted 17 February, 2025; originally announced February 2025.

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

Discovery of high-frequency quasi-periodic oscillation in short-duration gamma-ray bursts

Authors: Xing Yang, Hou-Jun Lü, Jared Rice, En-Wei Liang

Abstract: Rapidly rotating newborn magnetars, which originate from binary neutron star (NS) mergers and serve as the central engines of short gamma-ray bursts (GRBs), may leave some imprints on their prompt gamma-ray light curves even though they are far from their radiating fireballs. A high-frequency quasi-periodic oscillation (QPO) would be a unique feature for the magnetar central engine, especially a h… ▽ More Rapidly rotating newborn magnetars, which originate from binary neutron star (NS) mergers and serve as the central engines of short gamma-ray bursts (GRBs), may leave some imprints on their prompt gamma-ray light curves even though they are far from their radiating fireballs. A high-frequency quasi-periodic oscillation (QPO) would be a unique feature for the magnetar central engine, especially a hypermassive magnetar. By conducting a systematic analysis of the prompt gamma-ray light curves from 605 short GRBs observed by {\em Fermi}/Gamma-ray Burst Monitor, we have identified such QPO signals in three GRBs (e.g. GRB 120323A, GRB 181222B, and GRB 190606A). The QPOs that peaked at $1258^{+6}_{-6}$ Hz for GRB 120323A, $623^{+4}_{-4}$ Hz for GRB 181222B, and $1410^{+4}_{-5}$ Hz for GRB 190606A are all with a confidence level above 5.2 $σ$. The high-frequency QPO signals of those three short GRBs may be caused by a hypermassive magnetar acting as the central engine in a binary NS merger of a binary NS. △ Less

Submitted 19 February, 2025; v1 submitted 23 January, 2025; originally announced January 2025.

Comments: 11 pages, 1 Table, 7 Figures, accepted for publication in MNRAS, and matched with the published verison

Two repeated quasi-periodic oscillations in the FSRQ S5 1044+71 observed by TESS

Authors: Jun-Jie Wang, Ting-Feng Yi, Yangwei Zhang, He Lu, Yuncai Shen, Lisheng Mao, Liang Dong

Abstract: In this work, we report for the first time two repeated quasi-periodic oscillations (QPOs) in the light curve of the Flat Spectrum Radio Quasar (FSRQ) S5 1044+71. This source was observed by the Transiting Exoplanet Survey Satellite (TESS) in multiple sectors. We used the generalized Lomb-Scargle periodogram method and weighted wavelet Z-transform method to search for significant periodic signals.… ▽ More In this work, we report for the first time two repeated quasi-periodic oscillations (QPOs) in the light curve of the Flat Spectrum Radio Quasar (FSRQ) S5 1044+71. This source was observed by the Transiting Exoplanet Survey Satellite (TESS) in multiple sectors. We used the generalized Lomb-Scargle periodogram method and weighted wavelet Z-transform method to search for significant periodic signals. The main results are as follows: We found QPOs of $\sim$ 7.0 days (persisted for 4 cycles, with a significance of $\sim3.5σ$) and $\sim$ 7.3 days (persisted for 5 cycles, with a significance of $\sim3.8σ$) in the light curves of Sector 47 and EP1, respectively. Considering range of error, we consider them to be the same. We discussed two likely models of these rapid quasi-periodic variations: One comes from the jet and the other from the accretion disk. For the first one, we consider kink instability of the jet as a plausible explanation. Second, the QPO is probable to come from the main hot spots in the accretion disk, which is located approximately within the innermost stable circular orbit allowed by general relativity. Based on this model, we estimate the mass of the black hole in S5 1044+71 to be $3.49 \times 10^9 M_{\odot}$. △ Less

Submitted 21 January, 2025; originally announced January 2025.

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

Constraining the Nuclear Equation of State of neutron star via high-frequency quasi-periodic oscillation in short gamma-ray bursts

Authors: Jun-Xiang Huang, Hou-Jun Lü, En-Wei Liang

Abstract: The determination of the equation of state (EOS) of a neutron star (NS) and its maximum mass is very important for understanding the formation and properties of NSs under extreme conditions, but they remain open questions. Short-duration gamma-ray bursts (GRBs) are believed to originate from the merger of binary NSs or giant flares (GFs) of soft gamma repeaters (SGRs). Recently, the high-frequency… ▽ More The determination of the equation of state (EOS) of a neutron star (NS) and its maximum mass is very important for understanding the formation and properties of NSs under extreme conditions, but they remain open questions. Short-duration gamma-ray bursts (GRBs) are believed to originate from the merger of binary NSs or giant flares (GFs) of soft gamma repeaters (SGRs). Recently, the high-frequency quasi-periodic oscillations (QPOs) have been claimed to be identified from two short GRBs (GRB 931101B and GRB 910711). In this paper, we propose that the observed high-frequency QPOs in these two short GRBs result from torsional oscillations in the GFs of SGRs associated with cold NSs, or from radial oscillations of hypermassive NSs as the hot remnants of binary NS mergers, and then to constrain the EOS of NSs. For torsional oscillations, the six selected EOSs (TM1, NL3, APR, SLy4, DDME2, and GM1) of NSs suitable for the zero-temperature condition exhibit significant overlap in mass ranges, suggesting that we cannot constrain the EOS of NSs. For radial oscillations, the six selected EOSs (IUF, TM1, TMA, FSG, BHBLp, and NL3) of NSs suitable for the high-temperature condition cannot be ruled out when redshift is considered. However, it is found that the EOS can only be constrained if the redshift and temperature of the remnant can be measured. △ Less

Submitted 24 February, 2025; v1 submitted 17 January, 2025; originally announced January 2025.

Comments: 16 pages, 2 Tables, 10 Figures, accepted for publication in ApJ, and matched with the published verison

Waveform Reconstruction of Core-Collapse Supernova Gravitational Waves with Improved Multisynchrosqueezing Transform

Authors: Yong Yuan, Ao-Ran Wang, Zhuo-Tao Li, Gang Yu, Hou-Jun Lü, Peng Xu, Xi-Long Fan

Abstract: Gravitational waves (GWs) from core-collapse supernovae (CCSNe) have been proposed as a means to probe the internal physical properties of supernovae. However, due to their complex time-frequency structure, effectively searching for and extracting GW signals from CCSNe remains an unsolved challenge. In this paper, we apply the improved multisynchrosqueezing transform (IMSST) method to reconstruct… ▽ More Gravitational waves (GWs) from core-collapse supernovae (CCSNe) have been proposed as a means to probe the internal physical properties of supernovae. However, due to their complex time-frequency structure, effectively searching for and extracting GW signals from CCSNe remains an unsolved challenge. In this paper, we apply the improved multisynchrosqueezing transform (IMSST) method to reconstruct simulated GW data based on the advanced LIGO (aLIGO) and Einstein Telescope (ET) detectors. These data are generated by the magnetorotational and neutrino-driven mechanisms, and we use the match score as the criterion for evaluating the quality of the reconstruction. To assess whether the reconstructed waveforms correspond to true GW signals, we calculate the false alarm probability of reconstruction (FAPR). For GW sources located at 10 kpc and datasets where the waveform amplitudes are normalized to $5 \times 10^{-21}$ observed by aLIGO, FAPR are $2.1 \times 10^{-2}$ and $6.2 \times 10^{-3}$, respectively. For GW sources at 100 kpc and with waveform amplitudes normalized to $5 \times 10^{-21}$ observed by ET, FAPR are $1.3 \times 10^{-1}$ and $1.5 \times 10^{-2}$, respectively. When the gravitational wave strain reaches $7 \times 10^{-21}$ and the match score threshold is set to 0.75, the IMSST method achieves maximum reconstruction distances of approximately 37 kpc and 317 kpc for aLIGO and ET, respectively. Finally, we compared the performance of IMSST and STFT in waveform reconstruction based on the ET. The results show that the maximum reconstructable distance using STFT is 186 kpc. △ Less

Submitted 8 December, 2024; originally announced December 2024.

Ground electron calibration of the Gamma-ray Transient Monitor onboard DRO-A Satellite

Authors: Pei-Yi Feng, Zheng-Hua An, Yu-Hui Li, Qi Le, Da-Li Zhang, Xin-Qiao Li, Shao-Lin Xiong, Cong-Zhan Liu, Wei-Bin Liu, Jian-Li Wang, Bing-Lin Deng, He Xu, Hong Lu

Abstract: The Gamma-Ray Transient Monitor (GTM) is an all-sky monitor onboard the Distant Retrograde Orbit-A (DRO-A) satellite, with the scientific objective of detecting gamma-ray bursts in the energy range of 20 keV to 1 MeV. The GTM is equipped with five Gamma-Ray Transient Probes (GTPs), utilizing silicon photomultiplier (SiPM) arrays coupled with NaI(Tl) scintillators for signal readout. To test the pe… ▽ More The Gamma-Ray Transient Monitor (GTM) is an all-sky monitor onboard the Distant Retrograde Orbit-A (DRO-A) satellite, with the scientific objective of detecting gamma-ray bursts in the energy range of 20 keV to 1 MeV. The GTM is equipped with five Gamma-Ray Transient Probes (GTPs), utilizing silicon photomultiplier (SiPM) arrays coupled with NaI(Tl) scintillators for signal readout. To test the performance of the GTP in detecting electrons, we independently developed a continuous-energy-tunable, low-current, quasi-single-electron accelerator, and used this facility for ground-based electron calibration of the GTP. This paper provides a detailed description of the operational principles of the unique electron accelerator and comprehensively presents the process and results of electron calibration for the GTP. The calibration results indicate that the dead time for normal signals is less than 4 $μ$s, while for overflow signals, it is approximately 70 $μ$s, consistent with the design specifications. The GTP's time-recording capability is working correctly, accurately recording overflow events. The GTP responds normally to electrons in the 0.4-1.4 MeV energy range. The ground-based electron calibration validates the design of the GTP and enhances the probe's mass model, laying the foundation for payload development, in-orbit observation strategies, and scientific data analysis. △ Less

Submitted 28 November, 2024; originally announced November 2024.

Comments: 14 pages, 16 figures

Estimating Stellar Atmospheric Parameters and [α/Fe] for LAMOST O-M type Stars Using a Spectral Emulator

Authors: Jun-chao Liang, A-Li Luo, Yin-Bi Li, Xiao-Xiao Ma, Shuo Li, Shu-Guo Ma, Hai-Ling Lu, Yun-Jin Zhang, Bing Du, Xiao Kong

Abstract: In this paper, we developed a spectral emulator based on the Mapping Nearby Galaxies at Apache Point Observatory Stellar Library (MaStar) and a grouping optimization strategy to estimate effective temperature (T_eff), surface gravity (log g), metallicity ([Fe/H]) and the abundance of alpha elements with respect to iron ([alpha/Fe]) for O-M-type stars within the Large Sky Area Multi-Object Fiber Sp… ▽ More In this paper, we developed a spectral emulator based on the Mapping Nearby Galaxies at Apache Point Observatory Stellar Library (MaStar) and a grouping optimization strategy to estimate effective temperature (T_eff), surface gravity (log g), metallicity ([Fe/H]) and the abundance of alpha elements with respect to iron ([alpha/Fe]) for O-M-type stars within the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) low-resolution spectra. The primary aim is to use a rapid spectral-fitting method, specifically the spectral emulator with the grouping optimization strategy, to create a comprehensive catalog for stars of all types within LAMOST, addressing the shortcomings in parameter estimations for both cold and hot stars present in the official LAMOST AFGKM-type catalog. This effort is part of our series of studies dedicated to establishing an empirical spectral library for LAMOST. Experimental results demonstrate that our method is effectively applicable to parameter prediction for LAMOST, with the single-machine processing time within $70$ hr. We observed that the internal error dispersions for T_eff, log g, [Fe/H], and [alpha/Fe] across different spectral types lie within the ranges of $15-594$ K, $0.03-0.27$ dex, $0.02-0.10$ dex, and $0.01-0.04$ dex, respectively, indicating a good consistency. A comparative analysis with external data highlighted deficiencies in the official LAMOST catalog and issues with MaStar parameters, as well as potential limitations of our method in processing spectra with strong emission lines and bad pixels. The derived atmospheric parameters as a part of this work are available at https://nadc.china-vo.org/res/r101402/ . △ Less

Submitted 13 November, 2024; originally announced November 2024.

Comments: 34 pages, 28 Figures, published in ApJS

Observation of nonaxisymmetric standard magnetorotational instability induced by a free-shear layer

Authors: Yin Wang, Fatima Ebrahimi, Hongke Lu, Jeremy Goodman, Erik P. Gilson, Hantao Ji

Abstract: The standard magnetorotational instability (SMRI) is widely believed to be responsible for the observed accretion rates in astronomical disks. It is a linear instability triggered in the differentially rotating ionized disk flow by a magnetic field component parallel to the rotation axis. Most studies focus on axisymmetric SMRI in conventional base flows with a Keplerian profile for accretion disk… ▽ More The standard magnetorotational instability (SMRI) is widely believed to be responsible for the observed accretion rates in astronomical disks. It is a linear instability triggered in the differentially rotating ionized disk flow by a magnetic field component parallel to the rotation axis. Most studies focus on axisymmetric SMRI in conventional base flows with a Keplerian profile for accretion disks or an ideal Couette profile for Taylor-Couette flows, since excitation of nonaxisymmetric SMRI in such flows requires a magnetic Reynolds number Rm more than an order of magnitude larger. Here, we report that in a magnetized Taylor-Couette flow, nonaxisymmetric SMRI can be destabilized in a free-shear layer in the base flow at Rm $\gtrsim$ 1, the same threshold as for axisymmetric SMRI. Global linear analysis reveals that the free-shear layer reduces the required Rm, possibly by introducing an extremum in the vorticity of the base flow. Nonlinear simulations validate the results from linear analysis and confirm that a novel instability recently discovered experimentally (Nat. Commun. 13, 4679 (2022)) is the nonaxisymmetric SMRI. Our finding has astronomical implications since free-shear layers are ubiquitous in celestial systems. △ Less

Submitted 4 November, 2024; originally announced November 2024.

The X-ray re-brightening of GRB afterglow revisited: a possible signature from activity of the central engine

Authors: Zhe Yang, Hou-Jun Lü, Xing Yang, Jun Shen, Shuang-Xi Yi

Abstract: Long-duration gamma-ray bursts (GRBs) are thought to be from core collapse of massive stars, and a rapidly spinning magnetar or black hole may be formed as the central engine. The extended emission in the prompt emission, flares and plateaus in X-ray afterglow, are proposed to be as the signature of central engine re-activity. However, the directly evidence from observations of identifying the cen… ▽ More Long-duration gamma-ray bursts (GRBs) are thought to be from core collapse of massive stars, and a rapidly spinning magnetar or black hole may be formed as the central engine. The extended emission in the prompt emission, flares and plateaus in X-ray afterglow, are proposed to be as the signature of central engine re-activity. However, the directly evidence from observations of identifying the central engines remain an open question. In this paper, we systemically search for long-duration GRBs that consist of bumps in X-ray afterglow detected by Swift/XRT, and find that the peak time of the X-ray bumps exhibit bimodal distribution (defined as early and late bumps) with division line at $t=7190$ s. Although we cannot rule out that such a bimodality arises from selection effects. We proposed that the long-duration GRBs with an early (or late) bumps may be originated from the fall-back accretion onto a new-born magnetar (or black hole). By adopting MCMC method to fit the early (or late) bumps of X-ray afterglow with the fall-back accretion of magnetar (or black hole), it is found that the initial surface magnetic filed and period of magnetars for most early bumps are clustered around $5.88\times10^{13}$ G and $1.04$ ms, respectively. Meanwhile, the derived accretion mass of black hole for late bumps is range of $[4\times10^{-4}, 1.8\times10^{-2}]~M_{\odot}$, and the typical fall-back radius is distributed range of $[1.04, 4.23]\times 10^{11}$ cm which is consistent with the typical radius of a Wolf-Rayet star. However, we also find that the fall-back accretion magnetar model is disfavored by the late bumps, but the fall-back accretion of black hole model can not be ruled out to interpret the early bumps of X-ray afterglow. △ Less

Submitted 3 November, 2024; originally announced November 2024.

Comments: 14 pages, 5 figures, accepted for publication in MNRAS

Constraints on primordial black holes from $N_{\text{eff}}$ : scalar induced gravitational waves as an extra radiation component

Authors: Jing-Zhi Zhou, Yu-Ting Kuang, Zhe Chang, H. Lü

Abstract: In June 2023, multiple pulsar timing array collaborations provided evidence for the existence of a stochastic gravitational wave background. Scalar induced gravitational waves (SIGWs), as one of the most likely sources of stochastic gravitational waves, have received widespread attention. When primordial curvature perturbations on small scales are sufficiently large, \acp{PBH} inevitably form, con… ▽ More In June 2023, multiple pulsar timing array collaborations provided evidence for the existence of a stochastic gravitational wave background. Scalar induced gravitational waves (SIGWs), as one of the most likely sources of stochastic gravitational waves, have received widespread attention. When primordial curvature perturbations on small scales are sufficiently large, \acp{PBH} inevitably form, concurrently producing SIGWs with significant observable effects. These SIGWs can serve as an additional radiation component, influencing the relativistic degrees of freedom $N_{\text{eff}}$. Taking into account primordial non-Gaussianity, we study the energy density spectrum of SIGWs up to the third order and use the current observational data of $N_{\text{eff}}$ to constrain small-scale primordial curvature perturbations and the abundance of \acp{PBH}. △ Less

Submitted 5 January, 2025; v1 submitted 13 October, 2024; originally announced October 2024.

Comments: Accepted for publication in ApJ on 2025 January 5

Neutrino constraints and detection prospects from gamma-ray bursts with different jet compositions

Authors: Yang-Dong-Jun Ou, Hou-Jun Lü, Xue-Zhao Chang, Xiao-Xuan Liu, En-Wei Liang

Abstract: The prompt emission mechanism of gamma-ray bursts (GRBs) is a long-standing open question, and GRBs have been considered as potential sources of high-energy neutrinos. Despite many years of search for the neutrino events associated with GRBs from IceCube, there were no results. However, the absence of search results for neutrinos provides a unique opportunity to constrain the parameter space of GR… ▽ More The prompt emission mechanism of gamma-ray bursts (GRBs) is a long-standing open question, and GRBs have been considered as potential sources of high-energy neutrinos. Despite many years of search for the neutrino events associated with GRBs from IceCube, there were no results. However, the absence of search results for neutrinos provides a unique opportunity to constrain the parameter space of GRB jet models. In this paper, we chose four peculiar GRBs with two different types of jet composition to investigate neutrino emission. It is found that only GRB 211211A could be well constrained within the dissipative photosphere model. By adopting the specific parameters of the photosphere, one can obtain $\varepsilon _{\rm p}/\varepsilon_{\rm e}<8$ for $f_{\rm p}>0.2$ from GRB 211211A. For the Internal-collision-induced Magnetic Reconnection and Turbulence (ICMART) model, we can effectively constrain neither GRB 230307A nor GRB 080916C. Moreover, we also investigate the detection prospects of high-energy neutrinos from GRBs and find that it is difficult to detect at least one high-energy neutrino associated with GRBs from the ICMART model even during the IceCube-Gen2 operation. For the GRB 211211A-like events, it is possible to detect at least one neutrino coincident with the gravitational wave during the IceCube-Gen2 operation, if such an event originated from mergers of compact stars within the photosphere dissipation. △ Less

Submitted 23 November, 2024; v1 submitted 12 October, 2024; originally announced October 2024.

Comments: 15 pages, 10 figures, matched with the published verison

Journal ref: The Astrophysical Journal, 976: 174 (10pp), 2024 December 1

The thermal emission in short GRBs with extended emission observed by Fermi/GBM

Authors: Xue-Zhao Chang, Hou-Jun Lü, Xing Yang, Jia-Ming Chen, En-Wei Liang

Abstract: Short gamma-ray bursts (SGRBs) with extended emission (EE) are composed of initial main emission (ME) with a short-hard spike, followed by a long-lasting EE. Whether the ME and EE originated from the same origin or not, as well as the jet composition, remains an open question. In this paper, we present a systematic analysis of 36 gamma ray bursts (GRBs) in our sample, which are identified as the c… ▽ More Short gamma-ray bursts (SGRBs) with extended emission (EE) are composed of initial main emission (ME) with a short-hard spike, followed by a long-lasting EE. Whether the ME and EE originated from the same origin or not, as well as the jet composition, remains an open question. In this paper, we present a systematic analysis of 36 gamma ray bursts (GRBs) in our sample, which are identified as the category of SGRBs with EE as observed by Fermi/Gamma-ray Burst Monitor. By extracting time-integrated spectra of ME and EE with cutoff power-law or Band models for our sample, we find that 20 out of 36 SGRBs have $α$ values that exceed the death line (e.g., -2/3) of synchrotron emission within either ME or EE phases, and we suggest that the quasi-thermal component should exist in the prompt emission. Then, we extract the time-resolved spectra of our samples, but only four GRBs are bright enough to extract the time-resolved spectra. We find that both thermal and nonthermal emissions do exist in the prompt emission of those four bright GRBs, which suggests that a hybrid jet (e.g., matter and Poynting-flux outflow) in GRBs should exist. Moreover, strong positive correlations (e.g., $F_{\rm tot}-Γ$ and $F_{\rm tot}-kT$) in the time-resolved spectra of ME and EE for those four GRBs have been discovered. This indicates that the spectral evolution of both ME and EE seem to share similar behavior, possibly from the same physical origin. △ Less

Submitted 5 November, 2024; v1 submitted 12 October, 2024; originally announced October 2024.

Comments: 30 pages, 9 figures and 11 tables, matched with the published verison

Journal ref: The Astrophysical Journal Supplement Series, 275: 9 (21pp), 2024 November 5

Scalar induced gravitational waves in f(R) gravity

Authors: Jing-Zhi Zhou, Yu-Ting Kuang, Di Wu, Fei-Yu Chen, H. Lü, Zhe Chang

Abstract: We investigate the first and second order cosmological perturbation equations in f(R) modified gravity theory and provide the equation of motion of second order scalar induced gravitational waves. We find that the effects of modified gravity not only change the form of the equation of motion of second order scalar induced gravitational waves but also contribute an additional anisotropic stress ten… ▽ More We investigate the first and second order cosmological perturbation equations in f(R) modified gravity theory and provide the equation of motion of second order scalar induced gravitational waves. We find that the effects of modified gravity not only change the form of the equation of motion of second order scalar induced gravitational waves but also contribute an additional anisotropic stress tensor, composed of first order scalar perturbations, to the source term of the gravitational waves. We calculate the energy density spectrum of second order scalar induced gravitational waves in the HS model. Utilizing current pulsar timing array observational data, we perform a rigorous Bayesian analysis of the parameter space of the HS model. △ Less

Submitted 13 November, 2024; v1 submitted 11 September, 2024; originally announced September 2024.

Comments: Accepted for publication in JCAP on 2024 November 11

Induced gravitational waves for arbitrary higher orders: vertex rules and loop diagrams in cosmological perturbation theory

Authors: Jing-Zhi Zhou, Yu-Ting Kuang, Di Wu, H. Lü, Zhe Chang

Abstract: Gravitational waves induced by primordial perturbations serve as crucial probes for studying the early universe, providing a significant window into potential new physics during cosmic evolution. Due to the potentially large amplitudes of primordial perturbations on small scales, the contributions of high-order cosmological perturbations are highly significant. We propose a vertex approach applica… ▽ More Gravitational waves induced by primordial perturbations serve as crucial probes for studying the early universe, providing a significant window into potential new physics during cosmic evolution. Due to the potentially large amplitudes of primordial perturbations on small scales, the contributions of high-order cosmological perturbations are highly significant. We propose a vertex approach applicable to the study of induced gravitational waves for arbitrary higher orders. Using the vertex approach and tree diagrams, we can directly derive the explicit expressions of higher-order induced gravitational waves without involving the complex and lengthy calculations of higher-order cosmological perturbations. Correlations between different tree diagrams correspond to the loop diagrams of two-point correlation functions of induced gravitational waves. Our investigation reveals that one-particle reducible diagrams impact tensor-scalar induced gravitational waves while leaving scalar induced gravitational waves unaffected. △ Less

Submitted 22 March, 2025; v1 submitted 26 August, 2024; originally announced August 2024.

Comments: Accepted for publication in PRD on 2025 March 22

The Progenitor and Central Engine of short-duration GRB 201006A associated with a coherent radio flash

Authors: Xiao Tian, HouJun Lü, Yong Yuan, Xing Yang, HaoYu Yuan, ShuangXi Yi, WenLong Zhang, EnWei Liang

Abstract: Recently, the detection of a coherent radio flash associated with short-duration GRB 201006A, occurring 76.6 minutes after the burst, has attracted great attention. However, the physical origin of the coherent radio flash remains under debate. By reanalyzing its data observed by Fermi and Swift, we find that an early radio afterglow as the physical origin of the radio flash can be ruled out, but t… ▽ More Recently, the detection of a coherent radio flash associated with short-duration GRB 201006A, occurring 76.6 minutes after the burst, has attracted great attention. However, the physical origin of the coherent radio flash remains under debate. By reanalyzing its data observed by Fermi and Swift, we find that an early radio afterglow as the physical origin of the radio flash can be ruled out, but the coherent radio emission seems to be consistent with the hypothesis of a supramassive magnetar as the central engine collapsing into a black hole. Within this scenario, the derived magnetar surface magnetic field ($B_{\rm p}$) and the initial spin period ($P_{\rm 0}$) fall into a reasonable range but require a preferable low value of $η_{\rm R} = 10^{-7}$ or $10^{-6}$. Moreover, the calculated low-$\varepsilon$ value and $E_{\rm γ,iso}-E_{\rm p}$ correlation of GRB 201006A also supports the progenitor which is from the merger of compact stars. We also discuss the non-detected kilonova emission associated with GRB 201006A, and then compare with its upper limits of optical observations. △ Less

Submitted 11 October, 2024; v1 submitted 13 August, 2024; originally announced August 2024.

Comments: 14 pages, 7 figures, matched with the published verison

Journal ref: The Astronomical Journal, 974:133 (9pp), 2024 October 10

Characterizing the current systems in the Martian ionosphere

Authors: Jiawei Gao, Shibang Li, Anna Mittelholz, Zhaojin Rong, Moa Persson, Zhen Shi, Haoyu Lu, Chi Zhang, Xiaodong Wang, Chuanfei Dong, Lucy Klinger, Jun Cui, Yong Wei, Yongxin Pan

Abstract: When the solar wind interacts with the ionosphere of an unmagnetized planet, it induces currents that form an induced magnetosphere. These currents and their associated magnetic fields play a pivotal role in controlling the movement of charged particles, which is essential for understanding the escape of planetary ions. Unlike the well-documented magnetospheric current systems, the ionospheric cur… ▽ More When the solar wind interacts with the ionosphere of an unmagnetized planet, it induces currents that form an induced magnetosphere. These currents and their associated magnetic fields play a pivotal role in controlling the movement of charged particles, which is essential for understanding the escape of planetary ions. Unlike the well-documented magnetospheric current systems, the ionospheric current systems on unmagnetized planets remain less understood, which constrains the quantification of electrodynamic energy transfer from stars to these planets. Here, utilizing eight years of data from the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, we investigate the global distribution of ionospheric currents on Mars. We have identified two distinct current systems in the ionosphere: one aligns with the solar wind electric field yet exhibits hemispheric asymmetry perpendicular to the electric field direction; the other corresponds to the flow pattern of annually-averaged neutral winds. We propose that these two current systems are driven by the solar wind and atmospheric neutral winds, respectively. Our findings reveal that Martian ionospheric dynamics are influenced by the neutral winds from below and the solar wind from above, highlighting the complex and intriguing nature of current systems on unmagnetized planets. △ Less

Submitted 6 August, 2024; originally announced August 2024.

Comments: 20 pages, 6 figures

GRB 221009A/SN 2022xiw: A Supernova Obscured by a Gamma-Ray Burst Afterglow?

Authors: De-Feng Kong, Xiang-Gao Wang, WeiKang Zheng, Hou-Jun Lü, L. P. Xin, Da-Bin Lin, Jia-Xin Cao, Ming-Xuan Lu, B. Ren, Edgar P. Vidal, J. Y. Wei, En-Wei Liang, Alexei V. Filippenko

Abstract: We present optical photometry for the afterglow of GRB 221009A, in some respects the most extraordinary gamma-ray burst (GRB) ever observed. Good quality in the R-band light curve is obtained, covering 0.32-19.57 days since the Fermi-GBM trigger. We find that a weak bump emerges fromthe declining afterglow at $t \approx 11$ days; a supernova (SN) may be responsible. We use a smooth broken power-la… ▽ More We present optical photometry for the afterglow of GRB 221009A, in some respects the most extraordinary gamma-ray burst (GRB) ever observed. Good quality in the R-band light curve is obtained, covering 0.32-19.57 days since the Fermi-GBM trigger. We find that a weak bump emerges fromthe declining afterglow at $t \approx 11$ days; a supernova (SN) may be responsible. We use a smooth broken power-law and $^{56}\mathrm{Ni}$ model to fit the light curve. The best-fitting results reveal that the SN ejected a total mass of $M_\mathrm{ej} = 3.70 M_\odot$, a $^{56}\mathrm{Ni}$ mass of $M_\mathrm{Ni} = 0.23 M_\odot$, and a kinetic energy of $E_\mathrm{SN,K} = 2.35 \times 10^{52} \mathrm{erg}$. We also compare GRB 221009A with other GRB-SN events based on a GRB-associated SN sample, and find that only SN 2003lw and SN 2011kl can be obviously revealed in the afterglow of GRB 221009A by setting these objects at its distance. This suggests that a supernova (SN 2022xiw) is possibly obscured by the brighter afterglow emission from GRB 221009A. △ Less

Submitted 30 June, 2024; originally announced July 2024.

Generalized Entropy Implies Varying-G: Horizon Area Dependent Field Equations and Black Hole-Cosmology Coupling

Authors: Hengxin Lu, Sofia Di Gennaro, Yen Chin Ong

Abstract: When the Bekenstein-Hawking entropy is modified, ambiguity often arises concerning whether the Hawking temperature or the thermodynamic mass should be modified. The common practice, however, is to keep the black hole solution the same as that in general relativity. On the other hand, if Jacobson's method of deriving Einstein equations from thermodynamic is valid in the general settings, then given… ▽ More When the Bekenstein-Hawking entropy is modified, ambiguity often arises concerning whether the Hawking temperature or the thermodynamic mass should be modified. The common practice, however, is to keep the black hole solution the same as that in general relativity. On the other hand, if Jacobson's method of deriving Einstein equations from thermodynamic is valid in the general settings, then given a generalized entropy one should first derive the corresponding modified gravity, and then look for the compatible black hole solution before investigating its thermodynamics. We comment on some properties and subtleties in this approach. In particular, we point out that generically generalized entropy would lead to a varying effective gravitational "constant" theory, where $G_\text{eff}$ depends on the horizon area. We discuss in what ways such theories are discernible from general relativity despite its seemingly jarring differences, and how to make sense of area-dependent field equations. As a consequence we show that in the Jacobson's approach, the standard quantum gravitational logarithmic correction to Bekenstein-Hawking entropy is equivalent to a running gravitational "constant". A horizon area dependent $G_\text{eff}$ could also lead to a coupling between black hole masses and cosmological expansion, a scenario that has been studied recently in the literature, but so far lacks strong theoretical motivation. In the Tsallis case, we show that the thermodynamic mass for a Schwarzschild black hole is just a constant multiple of its ADM mass, which is considerably simpler than the approach not utilizing the Jacobson's method. △ Less

Submitted 17 January, 2025; v1 submitted 29 June, 2024; originally announced July 2024.

Comments: Minor changes; also added/fixed some references; published version

Journal ref: Annals of Phys. 474 (2025) 169914

Spin, inclination, and magnetic field evolution of magnetar population in vacuum and plasma-filled magnetospheres

Authors: Jun-Xiang Huang, Hou-Jun Lü, Jared Rice, En-Wei Liang

Abstract: Magnetars are potential energy sources or central engines for numerous transient phenomena in the Universe. How newborn magnetars evolve in different environments remains an open question. Based on both observed and candidate magnetars, it is found that the periods of all magnetars or candidates appear as a bimodal distribution, and are defined as the ``long-P'' and ``short-P'' magnetar subclasses… ▽ More Magnetars are potential energy sources or central engines for numerous transient phenomena in the Universe. How newborn magnetars evolve in different environments remains an open question. Based on both observed and candidate magnetars, it is found that the periods of all magnetars or candidates appear as a bimodal distribution, and are defined as the ``long-P'' and ``short-P'' magnetar subclasses, respectively. We find that for the ``short-P'' subclass of magnetars, the $\dot{P}$ values also appear as a bimodal distribution, and therefore can be classified as ``high-$\dot{P}$ short-P'' and ``low-$\dot{P}$ short-P'' magnetar subclasses. In this paper, we use Monte Carlo simulations to generate synthetic magnetar populations and investigate the evolution of the ``high-$\dot{P}$ short-P'' and ``low-$\dot{P}$ short-P'' magnetar subclasses by considering both the magnetar spin and inclination, as well as the decay of their magnetic field within their evolution in both vacuum and plasma-filled magnetospheres. We find that the magnetar evolution is dependent on both spin and magnetic field, but seems to be insensitive to inclination evolution and magnetospheric environment for the ``high-$\dot{P}$ short-P'' sub-class. In comparison for the case of ``high-$\dot{P}$ short-P'', the magnetar evolution is dependent on spin, magnetic field, and inclination evolution, as well as the magnetospheric environment. The best evolution model should be the case of inclination evolution in vacuum with a small value of $\overline{\mathrm{FOM}}$. The differences in the best-fit parameters also suggest that the ``high-$\dot{P}$ short-P'' and ``low-$\dot{P}$ short-P'' magnetar subclasses may be tracking with different evolution channels. △ Less

Submitted 18 June, 2024; v1 submitted 24 May, 2024; originally announced May 2024.

Comments: 18 pages, 2 Tables, and 9 Figures. PRD in press, and matched with the published verison

Cosmology with Higher-Derivative Gravities

Authors: H. Khodabakhshi, M. Farhang, F. Shojai, H. Lü

Abstract: We introduce an ingenious approach to explore cosmological implications of higher-derivative gravity theories. The key novelty lies in the characterization of the additional massive spin-0 modes constructed from Hubble derivatives as an effective density, with the corresponding pressure uniquely determined by energy conservation, while terms with no Hubble derivatives directly alter Friedmann equa… ▽ More We introduce an ingenious approach to explore cosmological implications of higher-derivative gravity theories. The key novelty lies in the characterization of the additional massive spin-0 modes constructed from Hubble derivatives as an effective density, with the corresponding pressure uniquely determined by energy conservation, while terms with no Hubble derivatives directly alter Friedmann equations. This classification of the various high-derivative contributions to Friedmann equations develops insight about their cosmological impacts and is essential for understanding the universe's evolution across energy scales. Various examples of higher-derivative gravity theories illustrate the power of this method in efficiently solving Friedmann equations and exploring new phenomena. Using CMB and BAO data, we apply this method to assess the observational feasibility of wall-bouncing universes, as predicted by scenarios with, e.g., certain third order modifications to general relativity. These models also provide an inflationary phase without the need to introduce extra scalar fields. △ Less

Submitted 14 May, 2024; v1 submitted 5 May, 2024; originally announced May 2024.

Comments: 9 pages, 1 table

Journal ref: Phys.Rev.D 110 (2024) 6, L061504

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

A complex node of the cosmic web associated with the massive galaxy cluster MACS J0600.1-2008

Authors: Lukas J. Furtak, Adi Zitrin, Johan P. Richard, Dominique Eckert, Jack Sayers, Harald Ebeling, Seiji Fujimoto, Nicolas Laporte, David Lagattuta, Marceau Limousin, Guillaume Mahler, Ashish K. Meena, Felipe Andrade-Santos, Brenda L. Frye, Mathilde Jauzac, Anton M. Koekemoer, Kotaro Kohno, Daniel Espada, Harry Lu, Richard Massey, Anna Niemiec

Abstract: MACS J0600.1-2008 (MACS0600) is an X-ray luminous, massive galaxy cluster at $z_{\mathrm{d}}=0.43$, studied previously by the REionization LensIng Cluster Survey (RELICS) and ALMA Lensing Cluster Survey (ALCS) projects which revealed a complex, bimodal mass distribution and an intriguing high-redshift object behind it. Here, we report on the results of a combined analysis of the extended strong le… ▽ More MACS J0600.1-2008 (MACS0600) is an X-ray luminous, massive galaxy cluster at $z_{\mathrm{d}}=0.43$, studied previously by the REionization LensIng Cluster Survey (RELICS) and ALMA Lensing Cluster Survey (ALCS) projects which revealed a complex, bimodal mass distribution and an intriguing high-redshift object behind it. Here, we report on the results of a combined analysis of the extended strong lensing (SL), X-ray, Sunyaev-Zeldovich (SZ), and galaxy luminosity-density properties of this system. Using new JWST and ground-based Gemini-N and Keck data, we obtain 13 new spectroscopic redshifts of multiply imaged galaxies and identify 12 new photometric multiple-image systems and candidates, including two multiply imaged $z\sim7$ objects. Taking advantage of the larger areal coverage, our analysis reveals an additional bimodal, massive SL structure which we measure spectroscopically to lie adjacent to the cluster and whose existence was implied by previous SL-modeling analyses. While based in part on photometric systems identified in ground-based imaging requiring further verification, our extended SL model suggests that the cluster may have the second-largest critical area and effective Einstein radius observed to date, $A_{\mathrm{crit}}\simeq2.16 \mathrm{arcmin}^2$ and $θ_{\mathrm{E}}=49.7''\pm5.0''$ for a source at $z_{\mathrm{s}}=2$, enclosing a total mass of $M(<θ_{\mathrm{E}})=(4.7\pm0.7)\times10^{14} \mathrm{M}_{\odot}$. These results are also supported by the galaxy luminosity distribution, the SZ and X-ray data. Yet another, probably related massive cluster structure, discovered in X-rays $5'$ (1.7 Mpc) further north, suggests that MACS0600 is part of an even larger filamentary structure. This discovery adds to several recent detections of massive structures around SL galaxy clusters and establishes MACS0600 as a prime target for future high-redshift surveys with JWST. △ Less

Submitted 10 August, 2024; v1 submitted 4 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS. Version updated to the accepted paper

The progenitor and central engine of a peculiar GRB 230307A

Authors: ZhaoWei Du, HouJun Lü, Yong Yuan, Xing Yang, EnWei Liang

Abstract: Recently, a lack of supernova-associated long-duration gamma-ray burst (GRB 230307A) at such a low redshift $z=0.065$, but associated with a possible kilonova emission, has attracted great attention. Its heavy element nucleosynthesis and the characteristic of soft X-ray emission suggests that the central engine of GRB 230307A is magnetar which is originated from a binary compact star merger. The c… ▽ More Recently, a lack of supernova-associated long-duration gamma-ray burst (GRB 230307A) at such a low redshift $z=0.065$, but associated with a possible kilonova emission, has attracted great attention. Its heavy element nucleosynthesis and the characteristic of soft X-ray emission suggests that the central engine of GRB 230307A is magnetar which is originated from a binary compact star merger. The calculated lower value of $\varepsilon \sim 0.05$ suggests that the GRB 230307A seems to be with ambiguous progenitor. The lower value of $f_{\rm eff}=1.23$ implies that the GRB 230307A is not likely to be from the effect of "tip of iceberg". We adopt the magnetar central engine model to fit the observed soft X-ray emission with a varying efficiency and find that the parameters constraints of magnetar fall into a reasonable range, i.e., $B<9.4\times10^{15}$ G and $P<2.5$ ms for $Γ_{\rm sat} = 10^3$, and $B<3.6\times10^{15}$ G and $P<1.05$ ms for $Γ_{\rm sat} = 10^4$. Whether the progenitor of GBR 230307A is from the mergers of neutron star - white dwarf (NS - WD) or neutron star - neutron star (NS - NS) remains unknown. △ Less

Submitted 16 February, 2024; originally announced February 2024.

Comments: 7 pages, 1 table, 4 figures. Accepted for publication in ApJ Letters, 962:L27, 2024

Ejecta-circumstellar medium interaction in high-density environment contribution to kilonova emission: Application to GRB 191019A

Authors: Suo-Ning Wang, Hou-Jun Lü, Yong Yuan, Hao-Yu Yuan, Jared Rice, Meng-Hua Chen, En-Wei Liang

Abstract: The nearby long-duration GRB 191019A recently detected by Swift lacks an associated supernova and belongs to a host galaxy with little star formation activity, suggesting that the origin of this burst is the result of a merger of two compact objects with dynamical interactions in a high-density medium of an active galactic nucleus. Given the potential motivation of this event, and given that it oc… ▽ More The nearby long-duration GRB 191019A recently detected by Swift lacks an associated supernova and belongs to a host galaxy with little star formation activity, suggesting that the origin of this burst is the result of a merger of two compact objects with dynamical interactions in a high-density medium of an active galactic nucleus. Given the potential motivation of this event, and given that it occurs in such a high-density environment, the ejecta-circumstellar medium (CSM) interaction cannot be ignored as possibly contributing to the kilonova emission. Here, we theoretically calculate the kilonova emission by considering the contribution of the ejecta-CSM interaction in a high-density environment. We find that the contribution to the kilonova emission from the ejecta-CSM interaction will dominate at a later time, and a smaller ejecta mass will have a stronger kilonova emission from the ejecta-CSM interaction. Moreover, we try to apply it to GRB 191019A, but we find that it is difficult to identify the possible kilonova emission from the observations, due to the contribution of the bright host galaxy. On the other hand, less injected mass (less than $M_{\rm ej}=2\times10^{-5}M_{\odot}$) will be required if one can detect the kilonova emission associated with a GRB 191019A-like event in the future. The {\em r}-process-powered and spin energy contributions from the magnetar are also discussed. △ Less

Submitted 20 February, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: 10 pages, 6 figures, updated the rederence list. ApJ in press, and matched with the published verison

Observational Feasibility of 4D Einstein-Gauss-Bonnet Cosmology: Bouncing and Non-Bouncing Universes

Authors: H. Khodabakhshi, M. Farhang, H. Lü

Abstract: This paper analyzes the possibility of bouncing and non-bouncing universes in the framework of four-dimensional Einstein-Gauss-Bonnet (4D-EGB) gravity, corresponding respectively to negative and positive coupling constants $λ$ of the Gauss-Bonnet term. We also use the Horndeski-type scalar-tensor theory to assess the role of a scalar charge $C$ as a geometrical contribution to the radiation in the… ▽ More This paper analyzes the possibility of bouncing and non-bouncing universes in the framework of four-dimensional Einstein-Gauss-Bonnet (4D-EGB) gravity, corresponding respectively to negative and positive coupling constants $λ$ of the Gauss-Bonnet term. We also use the Horndeski-type scalar-tensor theory to assess the role of a scalar charge $C$ as a geometrical contribution to the radiation in the Universe. We modify the expansion history of the universe to allow for modifications induced by the 4D-EGB gravity. Using Planck measurements of the cosmic microwave background anisotropies as well as various datasets of baryonic acoustic oscillations, we set the upper bounds $λ\le 10^{-16} \text{(km/s/Mpc)}^{-2} $ and $λ\le 10^{-30} \text{(km/s/Mpc)}^{-2} $ for the non-bouncing and bouncing scenarios. The upper limit in the latter case is mainly driven by the requirement to conservatively respect the thermal history at energy scales of the standard model of particle physics. We also find that the contribution of the geometrical radiation-like term of the model cannot exceed 10\% of the current radiation in the Universe. The possibility of an early inflationary phase produced by a single scalar field is also studied and found to be feasible in both bouncing and non-bouncing scenarios. This study shows the feasibility of a bouncing universe, even with normal matter sector, in the 4D-EGB gravity. More theoretical investigation is required to further explore possible observational predictions of the model that can distinguish between general relativity and 4D-EGB gravity. △ Less

Submitted 1 May, 2024; v1 submitted 22 January, 2024; originally announced January 2024.

Comments: 17 pages, 4 figues, 1 table

Journal ref: JCAP 05 (2024) 024

Detector performance of the Gamma-ray Transient Monitor onboard DRO-A Satellite

Authors: Pei-Yi Feng, Zheng-Hua An, Da-Li Zhang, Chen-Wei Wang, Chao Zheng, Sheng Yang, Shao-Lin Xiong, Jia-Cong Liu, Xin-Qiao Li, Ke Gong, Xiao-Jing Liu, Min Gao, Xiang-Yang Wen, Ya-Qing liu, Xiao-Yun Zhao, Fan Zhang, Xi-Lei Sun, Hong Lu

Abstract: Gamma-ray Transient Monitor (GTM) is an all-sky monitor onboard the Distant Retrograde Orbit-A (DRO-A) satellite with the scientific objective of detecting gamma-ray transients ranging from 20 keV to 1 MeV. GTM is equipped with 5 Gamma-ray Transient Probe (GTP) detector modules, utilizing the NaI(Tl) scintillator coupled with a SiPM array. To reduce the SiPM noise, GTP makes use of a dedicated dua… ▽ More Gamma-ray Transient Monitor (GTM) is an all-sky monitor onboard the Distant Retrograde Orbit-A (DRO-A) satellite with the scientific objective of detecting gamma-ray transients ranging from 20 keV to 1 MeV. GTM is equipped with 5 Gamma-ray Transient Probe (GTP) detector modules, utilizing the NaI(Tl) scintillator coupled with a SiPM array. To reduce the SiPM noise, GTP makes use of a dedicated dual-channel coincident readout design. In this work, we firstly studied the impact of different coincidence times on detection efficiency and ultimately selected the 500 ns time coincidence window for offline data processing. To test the performance of GTPs and validate the Monte Carlo simulated energy response, we conducted comprehensive ground calibration tests using Hard X-ray Calibration Facility (HXCF) and radioactive sources, including energy response, detection efficiency, spatial response, bias-voltage response, and temperature dependence. We extensively presented the ground calibration results, and validated the design and mass model of GTP detector. These work paved the road for the in-flight observation and science data analysis. △ Less

Submitted 10 September, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Comments: 15 pages, 25 figures

Journal ref: Sci. China-Phys. Mech. Astron. 67, 111013 (2024)

The Jet Composition of GRB 230307A: Poynting-Flux-Dominated Outflow?

Authors: Zhao-Wei Du, HouJun Lü, Xiaoxuan Liu, EnWei Liang

Abstract: The jet composition of GRB plays an important role in understanding the energy dissipation and radiation mechanisms in GRB physics, but it is poorly constrained from the observational data. Recently, an interesting long-duration GRB 230307A with redshift $z=$0.065 has attracted great attention. The lack of detected thermal emission and mini-structure of prompt emission lightcurve of this burst sug… ▽ More The jet composition of GRB plays an important role in understanding the energy dissipation and radiation mechanisms in GRB physics, but it is poorly constrained from the observational data. Recently, an interesting long-duration GRB 230307A with redshift $z=$0.065 has attracted great attention. The lack of detected thermal emission and mini-structure of prompt emission lightcurve of this burst suggest that the outflow is Poynting-flux-dominated and point towards the ICMART model. In this paper, we invoke two independent methods to investigate the jet composition of GRB 230307A. The high magnetization parameter ($σ>7$ or ever large) for$R_0=10^{10}$ cm that is used to suppress thermal component, strongly suggests that a significant fraction of the outflow energy is likely in a Poynting flux entrained with the baryonic matter. Moreover, it is found that the radiation efficiency of this burst for typical values $ε_e=0.1$ and $ε_B=0.01$ can reach as high as $~50\%$ which disfavors the internal shock model, but is consistent with ICMART model. Finally, a possible unified picture to produce GRB 230307A originated from a compact star merger is also discussed. △ Less

Submitted 12 January, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: 6 pages, 2 figures, 1 table, updated references, and matched with the published verison

Can fallback accretion on magnetar model power the X-ray flares simultaneously observed with gamma-rays of Gamma-ray bursts?

Authors: Wen-Yuan Yu, Hou-Jun Lü, Xing Yang, Lin Lan, Zhe Yang

Abstract: The prompt emission, X-ray plateau, and X-ray flares of Gamma-ray bursts (GRB) are thought to be from internal dissipation, and the magnetar as the central engine with propeller fallback accretion is proposed to interpret the observed phenomena of GRBs. In this paper, by systematically searching for X-ray emission observed by Swift/Xry Telescope, we find that seven robust GRBs include both X-ray f… ▽ More The prompt emission, X-ray plateau, and X-ray flares of Gamma-ray bursts (GRB) are thought to be from internal dissipation, and the magnetar as the central engine with propeller fallback accretion is proposed to interpret the observed phenomena of GRBs. In this paper, by systematically searching for X-ray emission observed by Swift/Xry Telescope, we find that seven robust GRBs include both X-ray flares and plateau emissions with measured redshift. More interestingly, the X-ray flares/bumps for those seven GRBs are simultaneously observed in the gamma-ray band. By adopting the propeller fallback accretion model to fit the observed data, it is found that the free parameters of two GRBs (140512A and 180329B) can be constrained very well, while in the other five cases, more or less, they are not all sufficiently constrained. On the other hand, this requires that the conversion efficiency of the propeller is to be two or three times higher than that of the spindown dipole radiation of the magnetar. If this is the case, it is contradictory to the expectation from the propeller model: namely, a dirtier ejecta should be less efficient in producing gamma-ray emissions. Our results hint that at least the magnetar central engine with propeller fallback accretion model cannot interpret very well both the GRB X-ray flares simultaneously observed in the gamma-ray band and the X-ray flares of GRBs with a high Lorentz factor. △ Less

Submitted 16 January, 2024; v1 submitted 10 January, 2024; originally announced January 2024.

Comments: 13 pages, 5 figures, 1 table, updated references, match with published version

Multiwavelength observation of an active M-dwarf star EV Lac and its stellar flare accompanied by a delayed prominence eruption

Authors: Shun Inoue, Teruaki Enoto, Kosuke Namekata, Yuta Notsu, Satoshi Honda, Hiroyuki Maehara, Jiale Zhang, Hong-Peng Lu, Hiroyuki Uchida, Takeshi Go Tsuru, Daisaku Nogami, Kazunari Shibata

Abstract: We conducted 4-night multiwavelength observations of an active M-dwarf star EV Lac on 2022 October 24$-$27 with simultaneous coverage of soft X-rays (NICER; 0.2$-$12 $\mathrm{keV}$, Swift XRT; 0.2$-$10 $\mathrm{keV}$), near-ultraviolet (Swift UVOT/UVW2; 1600$-$3500 Å), optical photometry (TESS; 6000$-$10000 Å), and optical spectroscopy (Nayuta/MALLS; 6350$-$6800 Å). During the campaign, we detecte… ▽ More We conducted 4-night multiwavelength observations of an active M-dwarf star EV Lac on 2022 October 24$-$27 with simultaneous coverage of soft X-rays (NICER; 0.2$-$12 $\mathrm{keV}$, Swift XRT; 0.2$-$10 $\mathrm{keV}$), near-ultraviolet (Swift UVOT/UVW2; 1600$-$3500 Å), optical photometry (TESS; 6000$-$10000 Å), and optical spectroscopy (Nayuta/MALLS; 6350$-$6800 Å). During the campaign, we detected a flare starting at 12:28 UTC on October 25 with its white-light bolometric energy of $3.4 \times 10^{32}$ erg. At about 1 hour after this flare peak, our $\mathrm{Hα}$ spectrum showed a blue-shifted excess component at its corresponding velocity of $\sim 100 \: \mathrm{km \: s^{-1}}$. This may indicate that the prominence erupted with a 1-hour delay of the flare peak. Furthermore, the simultaneous 20-second cadence near-ultraviolet and white-light curves show gradual and rapid brightening behaviors during the rising phase at this flare. The ratio of flux in NUV to white light at the gradual brightening was $\sim 0.49$, which may suggest that the temperature of the blackbody is low ($< 9000 \: \mathrm{K}$) or the maximum energy flux of a nonthermal electron beam is less than $5\times10^{11} \: \mathrm{erg \: cm^{-2} \: s^{-1}}$. Our simultaneous observations of NUV and white-light flare raise the issue of a simple estimation of UV flux from optical continuum data by using a blackbody model. △ Less

Submitted 30 December, 2023; originally announced January 2024.

Comments: 19 pages, 10 figures. Accepted for publication in Publications of the Astronomical Society of Japan

The Intrinsic Energy Resolution of LaBr$_3$(Ce) Crystal for GECAM

Authors: Pei-Yi Feng, Xi-Lei Sun, Zheng-Hua An, Cheng-Er Wang, Da-Li Zhang, Xin-Qiao Li, Chao Zheng, Shao-Lin Xiong, Hong Lu

Abstract: This study aims to provide an accurate estimation of the intrinsic resolution of LaBr$_3$(Ce) crystal through a combination of experimental and simulation methods. We re-analyzed the data from previous Wide-Angle Compton Coincidence (WACC) and Hard X-ray Calibration Facility (HXCF) experiments, conducted PMT Single-Photoelectron Calibration (SPEC) and radial non-uniformity (also called Spot Scanni… ▽ More This study aims to provide an accurate estimation of the intrinsic resolution of LaBr$_3$(Ce) crystal through a combination of experimental and simulation methods. We re-analyzed the data from previous Wide-Angle Compton Coincidence (WACC) and Hard X-ray Calibration Facility (HXCF) experiments, conducted PMT Single-Photoelectron Calibration (SPEC) and radial non-uniformity (also called Spot Scanning, SS) experiments to acquire new data, and combined these results with Geant4 simulations to isolate the contribution of each physical process to the total energy resolution, thereby allowing for a precise estimation of the scintillator's intrinsic resolution. For 100 keV X-rays, the total energy resolution of LaBr$_3$(Ce) crystal is 3.99% $\pm$ 0.04% (expressed as 1-$σ$), with statistical fluctuations and intrinsic resolution as the main components, contributing 2.47% $\pm$ 0.00% and 3.06% $\pm$ 0.06%, respectively. We identify two main sources of intrinsic resolution: one primarily due to non-proportional scintillation, contributing 2.28% $\pm$ 0.00%, and the other due to fluctuations in the energy transfer process, contributing 2.04% $\pm$ 0.08%. We quantified six components of the total energy resolution and reconstructed the photon response using Geant4. The consistency between the reconstructed relative light yield and the experimental measurements validated the mass model of the LaBr$_3$(Ce) detector used in the simulations. △ Less

Submitted 3 January, 2025; v1 submitted 30 December, 2023; originally announced January 2024.

Comments: 20 pages, 14 figures

The Energy Response of LaBr3(Ce), LaBr3(Ce,Sr) and NaI(Tl) Crystals for GECAM

Authors: Pei-Yi Feng, Xi-Lei Sun, Zheng-Hua An, Yong Deng, Cheng-Er Wang, Huang Jiang, Jun-Jie Li, Da-Li Zhang, Xin-Qiao Li, Shao-Lin Xiong, Chao Zheng, Ke Gong, Sheng Yang, Xiao-Jing Liu, Min Gao, Xiang-Yang Wen, Ya-Qing Liu, Yan-Bing Xu, Xiao-Yun Zhao, Jia-Cong Liu, Fan Zhang, Hong Lu

Abstract: The GECAM series of satellites utilize LaBr3(Ce), LaBr3(Ce,Sr), and NaI(Tl) crystals as sensitive materials for gamma-ray detectors (GRDs). To investigate the non-linearity in the detection of low-energy gamma rays and address errors in the E-C relationship calibration, comprehensive tests and comparative studies of the non-linearity of these three crystals were conducted using Compton electrons,… ▽ More The GECAM series of satellites utilize LaBr3(Ce), LaBr3(Ce,Sr), and NaI(Tl) crystals as sensitive materials for gamma-ray detectors (GRDs). To investigate the non-linearity in the detection of low-energy gamma rays and address errors in the E-C relationship calibration, comprehensive tests and comparative studies of the non-linearity of these three crystals were conducted using Compton electrons, radioactive sources, and mono-energetic X-rays. The non-linearity test results for Compton electrons and X-rays displayed substantial differences, with all three crystals showing higher non-linearity for X-rays and gamma-rays than for Compton electrons. Despite LaBr3(Ce) and LaBr3(Ce,Sr) crystals having higher absolute light yields, they exhibited a noticeable non-linear decrease in light yield, especially at energies below 400 keV. The NaI(Tl) crystal demonstrated excess light output in the 6~200 keV range, reaching a maximum excess of 9.2% at 30 keV in X-ray testing and up to 15.5% at 14 keV during Compton electron testing, indicating a significant advantage in the detection of low-energy gamma rays. Furthermore, this paper explores the underlying causes of the observed non-linearity in these crystals. This study not only elucidates the detector responses of GECAM, but also marks the inaugural comprehensive investigation into the non-linearity of domestically produced lanthanum bromide and sodium iodide crystals. △ Less

Submitted 27 December, 2023; originally announced December 2023.

Comments: 12pages, 16 figures

Constraining the Charge of a Black Hole with Electromagnetic Radiation from a Black Hole-Neutron Star System

Authors: Hao-Yu Yuan, Hou-Jun Lü, Jared Rice, En-Wei Liang

Abstract: Black hole-neutron star (BH-NS) mergers are expected to emit gravitational-wave (GW) and electromagnetic (EM) counterparts when the NS is tidally disrupted or plunges into the BH. Recently, GW 200105 and GW200115 were claimed as originating in BH-NS mergers, even GW 200105 remains in debate. Several optical source candidates are reported to possible associate with the two GW events, but not confir… ▽ More Black hole-neutron star (BH-NS) mergers are expected to emit gravitational-wave (GW) and electromagnetic (EM) counterparts when the NS is tidally disrupted or plunges into the BH. Recently, GW 200105 and GW200115 were claimed as originating in BH-NS mergers, even GW 200105 remains in debate. Several optical source candidates are reported to possible associate with the two GW events, but not confirmed yet. In this work, we assume that the BH is charged (the NS is naturally charged) and try to constrain the charge of the BH by using the possible associated EM emission from the charged BH and NS system working in the inspiral regime. We adopt electric and magnetic dipole radiations for the binaries which power a Poynting-flux-dominated outflow to accelerate electrons. Then, it produces the observed EM radiation via synchrotron radiation. We find that the conversion efficiency in the X-ray band is much higher than that of the ultraviolet (UV), near-infrared, and radio bands. The estimated maximum charge-to-mass ratio (the charge for unit mass) of the BH is $1.12\times 10^{-6}$ and $1.53\times 10^{-6}$ esu for the binary systems of GW200105 and GW200115, respectively, if magnetic field strength $B_{p}\lesssim ~10^{16}$ G and period $P>~1$ ms for the NS spin. △ Less

Submitted 24 September, 2023; originally announced September 2023.

Comments: 13 pages, 4 Tables, and 2 Figures. Accepted for publication in PRD

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

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

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

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

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

Waveform Reconstruction of Core-Collapse Supernovae Gravitational-Waves with Ensemble Empirical Mode Decomposition

Authors: Yong Yuan, Xi-Long Fan, Hou-Jun Lü, Yang-Yi Sun, Kai Lin

Abstract: The gravitational waves (GW) from core-collapse supernovae (CCSN) have been proposed as a probe to investigate physical properties inside of the supernova. However, how to search and extract the GW signals from core-collapse supernovae remains an open question due to its complicated time-frequency structure. In this paper, we apply the Ensemble Empirical Mode Decomposition (EEMD) method to decompo… ▽ More The gravitational waves (GW) from core-collapse supernovae (CCSN) have been proposed as a probe to investigate physical properties inside of the supernova. However, how to search and extract the GW signals from core-collapse supernovae remains an open question due to its complicated time-frequency structure. In this paper, we apply the Ensemble Empirical Mode Decomposition (EEMD) method to decompose and reconstruct simulated GW data generated by magnetorotational mechanism and neutrino-driven mechanism within the advanced LIGO, using the match score as the criterion for assessing the quality of the reconstruction. The results indicate that by decomposing the data, the sum of the first six intrinsic mode functions (IMFs) can be used as the reconstructed waveform. To determine the probability that our reconstructed waveform corresponds to a real GW waveform, we calculate the false alarm probability of reconstruction (FAPR). By setting the threshold of the match score to be 0.75, we obtain FAPR of GW sources at a distance of 5 kpc and 10 kpc to be $6\times10^{-3}$ and $1\times10^{-2}$ respectively. If we normalize the maximum amplitude of the GW signal to $5\times10^{-21}$, the FAPR at this threshold is $4\times10^{-3}$. Furthermore, in our study, the reconstruction distance is not equivalent to the detection distance. When the strain of GW reaches $7 \times 10^{-21}$, and the match score threshold is set at 0.75, we can reconstruct GW waveform up to approximately 36 kpc. △ Less

Submitted 22 February, 2024; v1 submitted 12 September, 2023; originally announced September 2023.

Comments: 9 pages, 6 figures. Accepted by MNRAS

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

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

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

Submitted 26 August, 2023; originally announced August 2023.

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

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

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

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

Submitted 26 August, 2023; originally announced August 2023.

Comments: 11 pages, 4 figures, 1 table

Comprehensive study of the blazars from Fermi-LAT LCR: The log-normal flux distribution and linear RMS-Flux relation

Authors: Na Wang, Ting-Feng Yi, Liang Wang, Li-Sheng Mao, Zhi-Yuan Pu, Gong-Ming Ning, Wei-Tian Huang, He Lu, Shun Zhang, Yu-Tong Chen, Liang Dong

Abstract: Fermi-LAT LCR provide continuous and regularly-sampled gamma-ray light curves, spanning about 14 years, for a large sample of blazars. The log-normal flux distribution and linear RMS-Flux relation of the light curves for a few of Fermi blazar have been examined in previous studies. However, the probability that blazars exhibit log-normal flux distribution and linear RMS-Flux relation in their gamm… ▽ More Fermi-LAT LCR provide continuous and regularly-sampled gamma-ray light curves, spanning about 14 years, for a large sample of blazars. The log-normal flux distribution and linear RMS-Flux relation of the light curves for a few of Fermi blazar have been examined in previous studies. However, the probability that blazars exhibit log-normal flux distribution and linear RMS-Flux relation in their gamma-ray light curves has not been systematically explored. In this study, we comprehensively research on the distribution of gamma-ray flux and the statistical characteristics on a large sample of 1414 variable blazars from the Fermi-LAT LCR catalog, including 572 FSRQs, 477 BL Lacs, and 365 BCUs, and statistically compare their flux distributions with normal and log-normal distributions. The results indicate that the probability of not reject log-normal is 42.05% for the large sample, and there is still 2.05% probability of not reject normality, based on the joint of Kolmogorov-Smirnov, Shapiro-Wilk and Normality tests. We further find that the probability that BL Lacs conforms to the log-normal distribution is higher than that of FSRQs. Besides, after removing sources with less than 200 data points from this large sample, a sample of 549 blazars, which is still a large sample comparing to the previous studies, was obtained. Basing on dividing the light curves into segments every 20 points (or 40 points, or one year), we fitted the linear RMS-Flux relation of this three different sets, and found that the Pearson correlation coefficients are all close to 1 of the most blazars. This result indicates a strong linear correlation between the RMS and the flux of this 549 blazars. The log-normal distribution and linear RMS-Flux relation indicate that the variability of gamma-ray flux for most blazars is non-linear and multiplicative process. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 13pages, 5figures, Accepted for publication in RAA