Researcher profile

Li-Ming Song

Li-Ming Song contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.HE astro-ph.IM

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author

15works

0followers

2topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Brightest GRB flare observed in GRB 221009A: bridge the last gap between flare and prompt emission in GRB

Flares are usually observed during the afterglow phase of Gamma-Ray Bursts (GRBs) in soft X-ray, optical and radio bands, but rarely in gamma-ray band. Despite the extraordinary brightness, GECAM-C has accurately measured both the bright prompt emission and flare emission of GRB 221009A without instrumental effects, offering a good opportunity to study the relation between them. In this work, we present a comprehensive analysis of flare emission of GRB 221009A, which is composed of a series of flares. Among them, we identify an exceptionally bright flare with a record-breaking isotropic energy $E_{\rm iso} = 1.82 \times 10^{53}$ erg of GRB flares. It exhibits the highest peak energy ever detected in GRB flares, $E_{\rm peak} \sim 300$ keV, making it a genuine gamma-ray flare. It also shows rapid rise and decay timescales, significantly shorter than those of typical X-ray flares observed in soft X-ray or optical band, but comparable to those observed in prompt emissions. Despite these exceptional properties, the flare shares several common properties with typical GRB flares. We note that this is the first observation of a GRB flare in the keV-MeV band with sufficiently high temporal resolution and high statistics, which bridges the last gap between prompt emission and flare.

preprint2026arXiv

On the Ultra-Long Gamma-Ray Transient GRB 250702B/EP250702

GRB 250702B/EP250702a is an interesting long-duration gamma-ray transient whose nature is in debate. To obtain a full picture in gamma-ray band, we implement a comprehensive targeted search of burst emission in a wide window of 30 days jointly with Insight-HXMT, GECAM and Fermi/GBM data within the ETJASMIN framework. In gamma-ray band, we find there is a 50-second precursor about 25 hours before the 4-hour main burst, which generally consists of 4 emission episodes. Remarkably, we find that the soft X-ray emission (after the main burst) decays as a power-law with start time aligning with the last episode of main emission and index of -5/3 perfectly consistent with the canonical prediction of fallback accretion. We conclude that the properties of precursor, main burst and the following soft X-ray emission strongly support the atypical collapsar Ultra-Long Gamma-Ray Burst (ULGRB) scenario rather than the Tidal Disruption Event (TDE), and all these gamma-ray and soft X-ray emission probably originate from relativistic jet whose luminosity is dominated by the fallback accretion rate during the death collapse of a supergiant star.

preprint2025arXiv

LUNCH: A Lightweight Unified Deep-Learning Framework for General Transients Classification in High-Energy Time-Domain Astronomy

The increasing data volume of high-energy space monitors necessitates real-time, automated transient classification for multi-messenger follow-up. Conventional methods rely on empirical features like hardness ratios and reliable localization, which are not always precisely available during early detection. We developed the Lightweight Unified Neural Classifier for High-energy Transients (LUNCH) - an end-to-end deep-learning framework that performs general transient classification directly from raw multi-band light curves, eliminating the need for background subtraction or source localization. Its dual-scale architecture fuses long- and short-scale temporal evolution adaptively. Evaluated on 15 years of Fermi/GBM triggers, the optimal model achieves 97.23% accuracy when trained on complete energy spectra. A lightweight version using only three broad energy bands retains 95.07% accuracy, demonstrating that coarse spectral information fused with temporal context enables robust discrimination. The system significantly outperforms the GBM in-flight classifier on three months of independent test data. Feature visualization reveals well-separated class clusters, confirming physical interpretability. LUNCH combines high accuracy, low computational cost, and instrument-agnostic inputs, offering a practical solution for real-time in-flight processing that enables timely triggers for immediate multi-wavelength and multi-messenger follow-up observations in future time-domain missions.

preprint2022arXiv

Determination of QPO properties in the presence of strong broad-band noise: a case study on the data of MAXI J1820+070

Accurate calculation of the phase lags of quasi-periodic oscillations (QPOs) will provide insight into their origin. In this paper we investigate the phase lag correction method which has been applied to calculate the intrinsic phase lags of the QPOs in MAXI J1820+070. We find that the traditional additive model between BBN and QPOs in the time domain is rejected, but the convolution model is accepted. By introducing a convolution mechanism in the time domain, the Fourier cross-spectrum analysis shows that the phase lags between QPOs components in different energy bands will have a simple linear relationship with the phase lags between the total signals, so that the intrinsic phase lags of the QPOs can be obtained by linear correction. The power density spectrum (PDS) thus requires a multiplicative model to interpret the data. We briefly discuss a physical scenario for interpreting the convolution. In this scenario, the corona acts as a low-pass filter, the Green's function containing the noise is convolved with the QPOs to form the low-frequency part of the PDS, while the high-frequency part requires an additive component. We use a multiplicative PDS model to fit the data observed by Insight-HXMT. The overall fitting results are similar compared to the traditional additive PDS model. Neither the width nor the centroid frequency of the QPOs obtained from each of the two PDS models were significantly different, except for the r.m.s. of the QPOs. Our work thus provides a new perspective on the coupling of noise and QPOs.

preprint2022arXiv

Gamma-Ray Polarimetry of the Crab Pulsar Observed by POLAR

The X/$γ$ ray polarimetry of the Crab pulsar/nebula is believed to hold crucial information on their emission models. In the past, several missions have shown evidence of polarized emission from the Crab. The significance of these measurements remains however limited. New measurements are therefore required. POLAR is a wide Field of View Compton-scattering polarimeter (sensitive in 50-500 keV) onboard the Chinese spacelab Tiangong-2 which took data from September 2016 to April 2017. Although not designed to perform polarization measurements of pulsars, we present here a novel method which can be applied to POLAR as well as that of other wide Field of View polarimeters. The novel polarimetric joint-fitting method for the Crab pulsar observations with POLAR, allows us to obtain constraining measurements of the pulsar component. The best fitted values and corresponding 1$σ$ deviations for the averaged phase interval: (PD=$14\substack{+15 \\ -10}$\%, PA=$108\substack{+33 \\ -54} ^{\circ}$), for Peak 1: (PD=$17\substack{+18 \\ -12}$\%, PA=$174\substack{+39 \\ -36} ^{\circ}$) and for Peak 2: (PD=$16\substack{+16 \\ -11}$\%, PA=$78\substack{+39 \\ -30} ^{\circ}$). Further more, the 3$σ$ upper limits on the polarization degree are for the averaged phase interval (55\%), Peak 1 (66\%) and Peak 2 (57\%). Finally, to illustrate the capabilities of this method in the future, we simulated two years observation to the Crab pulsar with POLAR-2. The results show that POLAR-2 is able to confirm the emission to be polarized with $5σ$ and $4σ$ confidence level if the Crab pulsar is polarized at $20\%$ and $10\%$ respectively.

preprint2022arXiv

GRB 210121A: Observation of photospheric emissions from different regimes and the evolution of the outflow

GRB 210121A was observed by Insight-HXMT, Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM), Fermi Gamma-ray Burst Monitor (Fermi/GBM) on Jan 21st 2021. In this work, photospheric emission from a structured jet is preferred to interpret the prompt emission phase of GRB 210121A and emissions from different regimes are observed on-axis. Particularly, the emission from the intermediate photosphere is first observed in the first 1.3 s of the prompt emission, while those from the other part are dominant by the saturated regime, and offers an alternative explanation compared with the previous work. Moreover, the emissions with considering the intermediate photosphere can well interpret the changes on low-energy photon index $α$ during the pulses. Besides, the evolution of the outflow is extracted from time-resolved analysis, and a correlation of $Γ_0 \propto L^{0.25\pm0.05}_0$ is obtained, which implies that the jet may be mainly launched by neutrino annihilation in a hyper-accretion disk.

preprint2022arXiv

iPTF14hls in the circumstellar medium interaction model: A promising candidate for a pulsational pair-instability supernova

iPTF14hls is a luminous Type II supernova (SN) with a bumpy light curve that remains debated for its origin. It maintains roughly a constant effective temperature and luminosity since discovery for about 600 days, followed by a slow decay. On $\sim 1000$\ days post discovery the light curve transitions to a very steep decline. A spectrum taken during this steep decline phase shows clear signatures of shock interaction with dense circumstellar medium (CSM). Here we explore the possibility of iPTF14hls as an interaction-powered SN. The light curve of iPTF14hls can be fitted with wind-like CSMs. Analytic modeling indicates that iPTF14hls may have undertaken six episodes of mass loss during the last $\sim 200\mathrm{yr}$. Assuming that the 1954 eruption triggered the last mass-loss episode, the stellar-wind velocity is determined to be $40-70\mathrm{km}\mathrm{s}^{-1}$, depending on different models. Mass loss rates are in the range $% 0.4-3.3M_{\odot }\mathrm{yr}^{-1}$. The inferred total mass of ejecta and CSMs ($M_{\mathrm{ej}}+M_{\mathrm{CSMs}}\simeq 245M_{\odot }$) supports the idea that iPTF14hls may be a candidate for a (pulsational) pair-instability SN. Discovery and observations of more similar stellar explosions will help understand these peculiar SNe.

preprint2022arXiv

The First Insight-HXMT Gamma-Ray Burst Catalog: The First Four Years

The Hard X-ray Modulation Telescope (Insight-HXMT), is China's first X-ray astronomy satellite launched on June 15, 2017. The anti-coincidence CsI detectors of the High Energy X-ray telescope (HE) onboard Insight-HXMT could serve as an all-sky gamma-ray monitor in about 0.2-3 MeV. In its first four years of operation, Insight-HXMT has detected 322 Gamma-Ray Bursts (GRBs) by offline search pipeline including blind search and targeted search. For the GOLDEN sample of Insight-HXMT GRBs, joint analyses were performed with other GRB missions, including Fermi Gamma-ray Burst Monitor (Fermi/GBM), Swift Burst Alert Telescope (Swift/BAT) and Gravitational wave high-energy Electromagnetic Counterpart All-sky Monitor (GECAM). It shows that Insight-HXMT can provide better constraint on GRB spectrum at higher energy band. The properties of Insight-HXMT GRBs are reported in detail, including their trigger time, duration, spectral parameters, peak fluxes of different time scales and fluence. This catalog is an official product of the Insight-HXMT GRB team.

preprint2020arXiv

A peculiar cyclotron line near 16 keV detected in the 2015 outburst of 4U 0115+63?

In 2015 October, the Be/X-ray binary 4U 0115+63 underwent a type II outburst, reaching an X-ray luminosity of $\sim$ {\DP{38}\es}. During the outburst, Nuclear Spectroscopic Telescope Array ({\it NuSTAR}) performed two Target of Opportunity observations. Using the broad-band spectra from {\it NuSTAR} (3$-$79 keV), we have detected multiple cyclotron lines of the source, i.e., $\sim$ 12, 16, 22 and 33/35 keV. Obviously, the 16 keV line is not a harmonic component of the 12 keV line. As described by the phase-dependent equivalent widths of these cyclotron lines, the 16 keV and 12 keV lines are two different fundamental lines. In our work, we apply the two-poles cyclotron line model to the observation, i.e., the two line sets are formed at the same altitude ($\sim$ 0.2 km over the NS surface) of different magnetic poles, with $\sim$ {1.1\TDP{12}} and {1.4\TDP{12}} G in two poles, respectively.

preprint2020arXiv

Background Model for the High-Energy Telescope of Insight-HXMT

Accurate background estimation is essential for spectral and temporal analysis in astrophysics. In this work, we construct the in-orbit background model for the High-Energy Telescope (HE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Based on the two-year blank sky observations of Insight-HXMT/HE, we first investigate the basic properties of the background and find that both the background spectral shape and intensity have long-term evolution at different geographical sites. The entire earth globe is then divided into small grids, each with a typical area of 5x5 square degrees in geographical coordinate system. For each grid, an empirical function is used to describe the long-term evolution of each channel of the background spectrum; the intensity of the background can be variable and a modification factor is introduced to account for this variability by measuring the contemporary flux of the blind detector. For a given pointing observation, the background model is accomplished by integrating over the grids that are passed by the track of the satellite in each orbit. Such a background model is tested with both the blank sky observations and campaigns for observations of a series of celestial sources. The results show an average systematic error of 1.5% for the background energy spectrum (26-100 keV) under a typical exposure of 8 ks, and <3% for background light curve estimation (30-150 keV). Therefore, the background model introduced in this paper is included in the Insight-HXMT software as a standard part specialized for both spectral and temporal analyses.

preprint2020arXiv

Background Model for the Low-Energy Telescope of Insight-HXMT

With more than 150 blank sky observations at high Galactic latitude, we make a systematic study to the background of the Low Energy Telescope (LE) of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). Both the on-ground simulation and the in-orbit observation indicate that the background spectrum mainly has two components. One is the particle background that dominates above 7 keV and its spectral shape is consistent in every geographical locations. Another is the diffuse X-ray background that dominates below 7 keV and has a stable spectrum less dependent of the sky region. The particle background spectral shape can be obtained from the blind detector data of all the blank sky observations, and the particle background intensity can be measured by the blind detector at 10-12.5 keV. The diffuse X-ray background in the high Galactic latitude can also be obtained from the blank sky spectra after subtracting the particle background. Based on these characteristics, we develop the background model for both the spectrum and the light curve. The systematic error for the background spectrum is investigated with different exposures (T_exp). For the spectrum with T_exp=1 ks, the average systematic errors in 1-7 keV and 1-10 keV are 4.2% and 3.7%, respectively. We also perform the systematic error analyses of the background light curves with different energy bands and time bins. The results show that the systematic errors for the light curves with different time bins are <8% in 1-10 keV.

preprint2020arXiv

Discovery of oscillations above 200 keV in a black hole X-ray binary with Insight-HXMT

Low-frequency quasi-periodic oscillations (LFQPOs) are commonly found in black hole X-ray binaries, and their origin is still under debate. The properties of LFQPOs at high energies (above 30 keV) are closely related to the nature of the accretion flow in the innermost regions, and thus play a crucial role in critically testing various theoretical models. The Hard X-ray Modulation Telescope (Insight-HXMT) is capable of detecting emissions above 30 keV, and is therefore an ideal instrument to do so. Here we report the discovery of LFQPOs above 200 keV in the new black hole MAXI J1820+070 in the X-ray hard state, which allows us to understand the behaviours of LFQPOs at hundreds of kiloelectronvolts. The phase lag of the LFQPO is constant around zero below 30 keV, and becomes a soft lag (that is, the high-energy photons arrive first) above 30 keV. The soft lag gradually increases with energy and reaches ~0.9s in the 150-200 keV band. The detection at energies above 200 keV, the large soft lag and the energy-related behaviors of the LFQPO pose a great challenge for most currently existing models, but suggest that the LFQPO probably originates from the precession of a small-scale jet.

preprint2020arXiv

The Background Model of the Medium Energy X-ray telescope of Insight-HXMT

The Medium Energy X-ray Telescope (ME) is one of the main payloads of the Hard X-ray Modulation Telescope (dubbed as Insight-HXMT). The background of Insight-HXMT/ME is mainly caused by the environmental charged particles and the background intensity is modulated remarkably by the geomagnetic field, as well as the geographical location. At the same geographical location, the background spectral shape is stable but the intensity varies with the level of the environmental charged particles. In this paper, we develop a model to estimate the ME background based on the ME database that is established with the two-year blank sky observations of the high Galactic latitude. In this model, the entire geographical area covered by Insight-HXMT is divided into grids of $5^{\circ}\times5^{\circ}$ in geographical coordinate system. For each grid, the background spectral shape can be obtained from the background database and the intensity can be corrected by the contemporary count rate of the blind FOV detectors. Thus the background spectrum can be obtained by accumulating the background of all the grids passed by Insight-HXMT during the effective observational time. The model test with the blank sky observations shows that the systematic error of the background estimation in $8.9-44.0$ keV is $\sim1.3\%$ for a pointing observation with an average exposure $\sim5.5$ ks. We also find that the systematic error is anti-correlated with the exposure, which indicates the systematic error is partly contributed by the statistical error of count rate measured by the blind FOV detectors.

preprint2020arXiv

The influence of the Insight-HXMT/LE time response on timing analysis

LE is the low energy telescope of Insight-HXMT. It uses swept charge devices (SCDs) to detect soft X-ray photons. The time response of LE is caused by the structure of SCDs. With theoretical analysis and Monte Carlo simulations we discuss the influence of LE time response (LTR) on the timing analysis from three aspects: the power spectral density, the pulse profile and the time lag. After the LTR, the value of power spectral density monotonously decreases with the increasing frequency. The power spectral density of a sinusoidal signal reduces by a half at frequency 536 Hz. The corresponding frequency for QPO signals is 458 Hz. The Root mean square (RMS) of QPOs holds the similar behaviour. After the LTR, the centroid frequency and full width at half maxima (FWHM) of QPOs signals do not change. The LTR reduces the RMS of pulse profiles and shifts the pulse phase. In the time domain, the LTR only reduces the peak value of the crosscorrelation function while it does not change the peak position. Thus it will not affect the result of the time lag. When considering the time lag obtained from two instruments and one among them is LE, a 1.18 ms lag is expected caused by the LTR. The time lag calculated in the frequency domain is the same as that in the time domain.

preprint2020arXiv

The POLAR Gamma-Ray Burst Polarization Catalog

Despite over 50 years of research, many open questions remain about the origin and nature of GRBs. Polarization measurements of the prompt emission of these extreme phenomena have long been thought to be the key to answering a range of these questions. The POLAR detector was designed to produce the first set of detailed and reliable polarization measurements in an energy range of approximately 50-500 keV. During late 2016 and early 2017, POLAR detected a total of 55 GRBs. Analysis results of 5 of these GRBs have been reported in the past. The results were found to be consistent with a low or unpolarized flux. However, previous reports by other collaborations found high levels of polarization. We study the polarization for all the 14 GRBs observed by POLAR for which statistically robust inferences are possible. Additionally, time-resolved polarization studies are performed on GRBs with sufficient apparent flux. A publicly available polarization analysis tool, developed within the 3ML framework, was used to produce statistically robust results. The method allows to combine spectral and polarimetric data from POLAR with spectral data from the Fermi GBM and Swift-BAT to jointly model the spectral and polarimetric parameters. The time integrated analysis finds all results to be compatible with a low or zero polarization with the caveat that, when time-resolved analysis is possible within individual pulses, we observe moderate polarization with a rapidly changing polarization angle. Thus, time-integrated polarization results, while pointing to lower polarization are potentially an artifact of summing over the changing polarization signal and thus, washing out the true moderate polarization. Therefore, we caution against over interpretation of any time-integrated results and encourage one to wait for more detailed polarization measurements from forthcoming missions such as POLAR-2 and LEAP.

Institution

Affiliation not imported yet

This author record came from a source that does not expose affiliation metadata. Once the author claims the profile or we enrich the record from another provider, this section will link to the concrete institution.

Topic footprint