Burst Light Curves Research Articles

A Comprehensive Analysis of Insight-HXMT Gamma-Ray Burst Data. I. Power Density Spectrum

The power density spectrum (PDS) is a powerful tool to study light curves of gamma-ray bursts (GRBs). We show the average PDS and individual PDS analysis with GRB data from the Hard X-ray Modulation Telescope (also named Insight-HXMT). The values of the power-law index of the average PDS ( αP¯ ) for long GRBs (LGRBs) vary from 1.58 to 1.29 (for 100–245, 245–600, and 600–2000 keV). The Insight-HXMT data allow us to extend the energy of the LGRBs up to 2000 keV, and a relation between αP¯ and energy E is obtained: αP¯∝E−0.09 (8–2000 keV). We first systematically investigate the average PDS and individual PDSs for short GRBs (SGRBs), and obtain αP¯∝E−0.07 (8–1000 keV), where the values of αP¯ vary from 1.86 to 1.34. The distribution of the power-law index of an individual PDS (α) of an SGRB is consistent with that of an LGRB, and the α value for the group with a dominant timescale (the bent power law) is higher than that for the group with no dominant timescale (the single power law). Both LGRBs and SGRBs show similar α and αP¯ , which indicates that they may be the result of similar stochastic processes. Typical values of the dominant timescale τ for LGRBs and SGRBs are 1.58 s and 0.02 s, respectively. It seems that τ varies in proportion to the duration of GRBs T 90, with a relation τ∝T900.86 . The GRB light curve may result from superposing a number of pulses with different timescales. No periodic or quasi-periodic signal above the 3σ significance threshold is found in our sample.

Insights from the Gaussian Process Method for the Fast Radio Burst–associated X-Ray Burst of SGR 1935+2154

The Gaussian process method is employed to analyze the light curves of bursts detected by Insight-HXMT, Neutron Star Interior Composition Explorer (NICER), and Gravitational Wave High-energy Electromagnetic Counterpart All-sky Monitor from SGR 1935+2154 between 2020 and 2022. It is found that a stochastically driven damped simple harmonic oscillator (SHO) is necessary to capture the characteristics of the X-ray bursts (XRBs). A variability timescale of the XRBs, corresponding to the broken frequencies in the SHO power spectral densities (PSDs), is extracted. In particular, a high broken frequency of 35 Hz where the index of the SHO PSD changes from −4 to −2 is constrained by the HXMT-HE burst associated with fast radio burst (FRB) 200428. It is suggested that the corresponding timescale of 0.03 s could be the retarding timescale of the system driven by some energy release, and the production of the HE photon should be quasi-simultaneous with the response. The other special event is a NICER burst with a retarding timescale of 1/(39 Hz) ≈ 0.02 s. In the normal XRBs, no retarding timescale is constrained; a long relax/equilibrium timescale (corresponding to a broken frequency of 1–10 Hz, where the index of the SHO PSD changes from −4/−2 to 0 in the SHO PSD) is obtained. The results indicate that the FRB-associated HXMT-HE XRB could be produced immediately when the system is responding to the energy disturbance, far before the equilibrium state.

The accretion burst of the massive young stellar object G323.46−0.08

Context. Accretion bursts from low-mass young stellar objects (YSOs) have been known for many decades. In recent years, the first accretion bursts of massive YSOs (MYSOs) have been observed. These phases of intense protostellar growth are of particular importance for studying massive star formation. Bursts of MYSOs are accompanied by flares of Class II methanol masers (hereafter masers), which are caused by an increase in exciting mid-infrared (MIR) emission. They can lead to long-lasting thermal afterglows of the dust continuum radiation visible at infrared (IR) and (sub)millimeter (hereafter (sub)mm) wavelengths. Furthermore, they might cause a scattered light echo. The G323.46−0.08 (hereafter G323) event, which shows all these features, extends the small sample of known MYSO bursts. Aims. Maser observations of the MYSO G323 show evidence of a flare, which was presumed to be caused by an accretion burst. This should be verified with IR data. We used time-dependent radiative transfer (TDRT) to characterize the heating and cooling timescales for eruptive MYSOs and to infer the main burst parameters. Methods. Burst light curves, as well as the pre-burst spectral energy distribution (SED) were established from archival IR data. The properties of the MYSO, including its circumstellar disk and envelope, were derived by using static radiative transfer modeling of pre-burst data. For the first time, TDRT was used to predict the temporal evolution of the SED. Observations with SOFIA/HAWC+ were performed to constrain the burst energy from the strength of the thermal afterglow. Image subtraction and ratioing were applied to reveal the light echo. Results. The G323 accretion burst is confirmed. It reached its peak in late 2013/early 2014 with a Ks-band increase of ∼2.5 mag. Both Ks-band and integrated maser flux densities follow an exponential decay. TDRT indicates that the duration of the thermal afterglow in the far-infrared (FIR) can exceed the burst duration by years. The latter was proved by SOFIA observations, which indicate a flux increase of (14.2 ± 4.6)% at 70 μm and (8.5 ± 6.1)% at 160 μm in 2022 (2 yr after the burst ended). A one-sided light echo emerged that was propagating into the interstellar medium. Conclusions. The burst origin of the G323 maser flare has been verified. TDRT simulations revealed the strong influence of the burst energetics and the local dust distribution on the strength and duration of the afterglow. The G323 burst is probably the most energetic MYSO burst that has been observed so far. Within 8.4 yr, an energy of (0.9−0.8+2.5) × 1047 erg was released. The short timescale points to the accretion of a compact body, while the burst energy corresponds to an accumulated mass of at least (7−6+20) MJup and possibly even more if the protostar is bloated. In this case, the accretion event might have triggered protostellar pulsations, which give rise to the observed maser periodicity. The associated IR light echo is the second observed from a MYSO burst.

Long gamma-ray burst light curves as the result of a common stochastic pulse--avalanche process

The complexity and variety exhibited by the light curves of long gamma-ray bursts (GRBs) enclose a wealth of information that has not yet been fully deciphered. Despite the tremendous advance in the knowledge of the energetics, structure, and composition of the relativistic jet that results from the core collapse of the progenitor star, the nature of the inner engine, how it powers the relativistic outflow, and the dissipation mechanisms remain open issues. A promising way to gain insights is describing GRB light curves as the result of a common stochastic process. In the Burst And Transient Source Experiment (BATSE) era, a stochastic pulse avalanche model was proposed and tested through the comparison of ensemble-average properties of simulated and real light curves. Here our aim was to revive and further test this model. We applied it to two independent datasets, BATSE and Swift /BAT, through a machine learning approach: the model parameters are optimised using a genetic algorithm. The average properties were successfully reproduced. Notwithstanding the different populations and passbands of both datasets, the corresponding optimal parameters are interestingly similar. In particular, for both sets the dynamics appear to be close to a critical state, which is key to reproducing the observed variety of time profiles. Our results propel the avalanche character in a critical regime as a key trait of the energy release in GRB engines, which underpins some kind of instability.

On the Detection and Characterization of Quasiperiodic Oscillations in Astronomical Time Series: Gamma-Ray Burst X-Ray Light Curves as a Test Case

The study of temporal properties of variable sources can elucidate their physical processes. In this context, we present a critical study comparing three approaches to periodic or quasiperiodic behavior: Gaussian process, power spectrum, and wavelet analysis, using celerite, Lomb–Scargle periodograms, and weighted wavelet Z-transforms, respectively. We use 15 Swift X-ray Telescope light curves of short gamma-ray bursts (sGRBs) as examples. A comprehensive analysis of two sGRB X-ray light curves is performed. The results reveal the importance of artifacts, largely in the form of false quasiperiodic oscillation signals, possibly introduced by preprocessing (such as detrending) or other aspects of the analysis. The exploration described in this paper can be helpful for future studies of variability in gamma-ray bursts, active galactic nuclei, and other astronomical sources.

Insight-HXMT observations of thermonuclear X-ray bursts in 4U 1636−53

ABSTRACT We conducted an analysis of 45 bursts observed from 4U 1636−53. To investigate the mechanism behind the light-curve profiles and the impact of thermonuclear X-ray bursts on the accretion environment in accreting neutron star low-mass X-ray binaries. This analysis employed both light-curve and time-resolved spectroscopy methodologies, with data collected by the Insight-Hard X-ray Modulation Telescope instrument. We found that 30 bursts exhibited similar light-curve profiles and were predominantly in the hard state, and two photospheric radius expansion (PRE) bursts were in the soft state. The light curves of most bursts did not follow a single exponential decay but displayed a dual-exponential behaviour. The initial exponent had a duration of approximately 6 s. We utilized both the standard method and the ‘fa’ method to fit the burst spectra. The majority of the ‘fa’ values exceeded 1, indicating an enhancement of the persistent emission during the burst. Under the two Comptonization components assumption, we suggest that the scattering of burst photons by the inner corona may mainly contribute to the persistent emission enhancement. We also observed an inverse correlation between the maximum fa and the persistent emission flux in the non-PRE burst. This anticorrelation suggests that when the accretion rate is lower, there is a greater enhancement of persistent emission during the burst peak. The prediction based on Poynting–Robertson drag (P–R drag) aligns with this observed anticorrelation.

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

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

Quasiperiodic Peak Energy Oscillations in X-Ray Bursts from SGR J1935+2154

Magnetars are young neutron stars powered by the strongest magnetic fields in the Universe (1013–15 G). Their transient X-ray emission usually manifests as short (a few hundred milliseconds), bright, energetic (∼1040–41 erg) X-ray bursts. Since its discovery in 2014, SGR J1935+2154 has become one of the most prolific magnetars, exhibiting very active bursting episodes and other fascinating events, such as pulse timing antiglitches and fast radio bursts. Here we present evidence for possible 42 Hz (24 ms) quasiperiodic oscillations in the ν F ν spectrum peak energy (E p ) identified in a unique burst detected with the Fermi Gamma-ray Burst Monitor in 2022 January. While quasiperiodic oscillations have been previously reported in the intensity of magnetar burst light curves, quasiperiodic oscillations in E p have not. We also find an additional event from the same outburst that appears to exhibit a similar character in E p , albeit of lower statistical quality. For these two exceptional transients, such E p oscillations can be explained by magnetospheric density and pressure perturbations. For burst-emitting plasma consisting purely of e + e − pairs, these acoustic modes propagate along a highly magnetized flux tube of length up to around L ∼ 130 neutron star radii, with L being lower if ions are present in the emission zone. Detailed time-resolved analyses of other magnetar bursts are encouraged to evaluate the rarity of these events and their underlying mechanisms.

Evidence for a luminosity–decay correlation in GRB GeV light curves

ABSTRACT Correlations between intrinsic properties of gamma-ray burst (GRB) light curves provide clues to the nature of the central engine, the jet, and a possible means to standardize GRBs for cosmological use. Here, we report on the discovery of a correlation between the intrinsic early-time luminosity, LG,10 s, measured at rest frame 10 s, and the average decay rate measured from rest frame 10 s onward, $\alpha _{\mathrm{G,avg\gt 10\, s}}$, in a sample of 13 Fermi Large Area Telescope long GRB light curves. We note that our selection criteria, in particular the requirement for a redshift to construct luminosity light curves, naturally limits our sample to energetic GRBs. A Spearman’s rank correlation gives a coefficient of –0.74, corresponding to a confidence level of 99.6 per cent, indicating that brighter afterglows decay faster than less luminous ones. Assuming a linear relation with log(LG,10s), we find $\alpha _{\mathrm{G,avg\gt 10\, s}}$$= -0.31_{-0.09}^{+0.12}\log ($LG,10s$) + 14.43_{-5.97}^{+4.55}$. The slope of −0.31 is consistent at 1σ with previously identified correlations in the optical/ultraviolet and X-ray light curves. We speculate that differences in the rate at which energy is released by the central engine or differences in observer viewing angle may be responsible for the correlation.

The Temporal Symmetrical and Translational Structure in Gamma-Ray Burst Light Curves

Tremendous information is hidden in the light curve of a gamma-ray burst (GRB). Based on Compton Gamma Ray Observatory/Burst And Transient Source Experiment (BATSE) data, Hakkila found a majority of GRBs can be characterized by a smooth, single-peaked component superposed with a temporally symmetrical residual structure, i.e., a mirror feature for the fast-varying component. In this study, we conduct a similar analysis on the same data, as well as on Fermi/Gamma-Ray Burst Monitor data. We obtained a similar conclusion, which is that most GRBs have this symmetrical fast-varying component. Furthermore, we chose an alternative model to characterize the smooth component and used a three-parameter model to identify the residual, i.e., the fast component. By choosing 226 BATSE GRBs based on a few criteria, we checked the time-symmetrical feature and time-translational feature for the fast components and found the ratio is roughly 1:1. We propose that both features could come from the structure of the ejected shells. In the future, the Square Kilometre Array might be able to observe the early radio emission from the collision of the shells.

Observing Supernova Neutrino Light Curves with Super-Kamiokande. IV. Development of SPECIAL BLEND: A New Public Analysis Code for Supernova Neutrinos

Supernova neutrinos are invaluable signals that offer information about the interior of supernovae. Because a nearby supernova can occur at any time, preparing for future supernova neutrino observation is an urgent task. For the prompt analysis of supernova neutrinos, we have developed a new analysis code, the “Supernova Parameter Estimation Code based on Insight on Analytic Late-time Burst Light curve at Earth Neutrino Detector” (SPECIAL BLEND). This code estimates the parameters of supernovae based on an analytic model of supernova neutrinos from the proto-neutron star cooling phase. For easy availability to the community, this code is public and easily runs in web environments. SPECIAL BLEND can estimate the parameters better than the analysis pipeline we developed in a previous paper. By using SPECIAL BLEND, we can estimate the supernova parameters within 10% precision up to ∼20 and ∼60 kpc (Large Magellanic Cloud contained) with Super-Kamiokande and Hyper-Kamiokande, respectively.

The Imprint of Convection on Type I X-Ray Bursts: Pauses in Photospheric Radius Expansion Lightcurves

Motivated by the recent observation by NICER of a type I X-ray burst from SAX J1808.4–3658 with a distinct “pause” feature during its rise, we show that bursts which ignite in a helium layer underneath a hydrogen-rich shell naturally give rise to such pauses, as long as enough energy is produced to eject the outer layers of the envelope by super-Eddington winds. The length of the pause is determined by the extent of the convection generated after ignition, while the rate of change of luminosity following the pause is set by the hydrogen gradient left behind by convection. Using the MESA stellar evolution code, we simulate the accumulation, nuclear burning, and convective mixing prior to and throughout the ignition of the burst, followed by the hydrodynamic wind. We show that the results are sensitive to the treatment of convection adopted within the code. In particular, the efficiency of mixing at the H/He interface plays a key role in determining the shape of the lightcurve. The data from SAX J1808.4–3658 favor strong mixing scenarios. Multidimensional simulations will be needed to properly model the interaction between convection and nuclear burning during these bursts, which will then enable a new way to use X-ray burst lightcurves to study neutron star surfaces.

A Concept of Assessment of LIV Tests with THESEUS Using the Gamma-Ray Bursts Detected by Fermi/GBM

According to Einstein’s special relativity theory, the speed of light in a vacuum is constant for all observers. However, quantum gravity effects could introduce its dispersion depending on the energy of photons. The investigation of the spectral lags between the gamma-ray burst (GRB) light curves recorded in distinct energy ranges could shed light on this phenomenon: the lags could reflect the variation of the speed of light if it is linearlydependent on the photon energy and a function of the GRB redshift. We propose a methodology to start investigating the dispersion law of light propagation in a vacuum using GRB light curves. This technique is intended to be fully exploited using the GRB data collected with THESEUS.

The Impacts of Neutron Star Structure and Base Heating on Type I X-Ray Bursts and Code Comparison

Type I X-ray bursts are rapidly brightening phenomena triggered by thermonuclear burning on the accreting layers of a neutron star (NS). The light curves represent the physical properties of NSs and the nuclear reactions on the proton-rich nuclei. The numerical treatments of the accreting NS and physics of the NS interior are not established, which shows uncertainty in modeling for observed X-ray light curves. In this study, we investigate theoretical X-ray burst models compared with burst light curves with GS 1826-24 observations. We focus on the impacts of the NS mass and radius and base heating on the NS surface using the MESA code. We find a monotonic correlation between the NS mass and the parameters of the light curve. The higher the mass, the longer the recurrence time and the greater the peak luminosity. While the larger the radius, the longer the recurrence time, the peak luminosity remains nearly constant. In the case of increasing base heating, both the recurrence time and peak luminosity decrease. We also examine the above results with a different numerical code, HERES, based on general relativity and consider the central NS. We find that the burst rate, energy, and strength are almost the same in two X-ray burst codes by adjusting the base heat parameter in MESA (the relative errors ≲5%), while the duration and rise times are significantly different between (the relative error is possibly ∼50%). The peak luminosity and the e-folding time change irregularly between two codes for different accretion rates.

TexAT detector upgrade for 14O([formula omitted])17F cross section measurement

A direct cross-section measurement of the 14O(α,p)17F reaction is important to understand the light curves of x-ray bursts. The measurement will be performed using the Texas Active Target TPC version 2 (TexAT_v2). The TexAT_v2 aims at measuring lower energy protons from the reaction than the original TexAT. Newly developed silicon and CsI(Tl) detector arrays are added at the left, right and bottom of a modified field cage to increase its detection efficiency. This paper describes the overall specifications and two commissioning experiments performed at Texas A&M University.

First Direct Measurement Constraining the ^{34}Ar(α,p)^{37}K Reaction Cross Section for Mixed Hydrogen and Helium Burning in Accreting Neutron Stars.

The rate of the final step in the astrophysical αp process, the ^{34}Ar(α,p)^{37}K reaction, suffers from large uncertainties due to a lack of experimental data, despite having a considerable impact on the observable light curves of x-ray bursts and the composition of the ashes of hydrogen and helium burning on accreting neutron stars. We present the first direct measurement constraining the ^{34}Ar(α,p)^{37}K reaction cross section, using the Jet Experiments in Nuclear Structure and Astrophysics gas jet target. The combined cross section for the ^{34}Ar,Cl(α,p)^{37}K,Ar reaction is found to agree well with Hauser-Feshbach predictions. The ^{34}Ar(α,2p)^{36}Ar cross section, which can be exclusively attributed to the ^{34}Ar beam component, also agrees to within the typical uncertainties quoted for statistical models. This indicates the applicability of the statistical model for predicting astrophysical (α,p) reaction rates in this part of the αp process, in contrast to earlier findings from indirect reaction studies indicating orders-of-magnitude discrepancies. This removes a significant uncertainty in models of hydrogen and helium burning on accreting neutron stars.

Mass measurements show slowdown of rapid proton capture process at waiting-point nucleus 64Ge

AbstractX-ray bursts are among the brightest stellar objects frequently observed in the sky by space-based telescopes. A type-I X-ray burst is understood as a violent thermonuclear explosion on the surface of a neutron star, accreting matter from a companion star in a binary system. The bursts are powered by a nuclear reaction sequence known as the rapid proton capture process (rp process), which involves hundreds of exotic neutron-deficient nuclides. At so-called waiting-point nuclides, the process stalls until a slower β+ decay enables a bypass. One of the handful of rp process waiting-point nuclides is 64Ge, which plays a decisive role in matter flow and therefore the produced X-ray flux. Here we report precision measurements of the masses of 63Ge, 64,65As and 66,67Se—the relevant nuclear masses around the waiting-point 64Ge—and use them as inputs for X-ray burst model calculations. We obtain the X-ray burst light curve to constrain the neutron-star compactness, and suggest that the distance to the X-ray burster GS 1826–24 needs to be increased by about 6.5% to match astronomical observations. The nucleosynthesis results affect the thermal structure of accreting neutron stars, which will subsequently modify the calculations of associated observables.

The spectral analysis and study of GRB 120709A, a burst with three distinct emission episodes

The reason for quiescent intervals and isolated emission episodes in gamma-ray burst (GRB) light curves is still not well understood. Probing these questions further requires the study of the spectral characteristics and properties of the outflow of such GRBs. GRB 120709A is a burst with three distinct emission episodes, separated by quiescent periods lasting several seconds. The episodes are comparable with each other in terms of peak flux and duration. The GRB was observed by the Fermi LAT and GBM instruments, covering over six orders of magnitude of energy. In this work, we carry out the time-integrated and time-resolved spectral analysis of GRB 120709A. We use these results to study the temporal variation of its parameters. In particular, we find that the peak energy of the spectrum of each episode traces the flux of those episodes. We also study and parameterise the correlation between the isotropic luminosity and intrinsic peak energy. Finally, we determine the parameters of the outflow based on gamma-gamma opacity calculations. We find that if we require the MeV emission to occur above the photosphere, then the ratio of the GeV and MeV emission radii can not be greater than ∼1.2. This implies an internal origin of the GeV emission. This also allows us to estimate the average value of the Lorentz factor to be in the ∼100–200 range. We also evaluate the photospheric and MeV emission radii and find them to be in the (4.5–5.3) ×1012 cm and (2.1–5.2) ×1014 cm ranges, respectively. We note the limitations of this exercise due to the relative faintness of the GRB which prevents us from doing a very fine time-resolved analysis or fitting the spectra with multi-component models. We also discuss GRB 120709A in the context of models for distinct emission episodes and quiescent intervals.

Multi-messenger detection prospects of gamma-ray burst afterglows with optical jumps

Some afterglow light curves of gamma-ray bursts (GRBs) exhibit very complex temporal and spectral features, such as a sudden intensity jump about one hour after the prompt emission in the optical band. We assume that this feature is due to the late collision of two relativistic shells and investigate the corresponding high-energy neutrino emission within a multi-messenger framework, while contrasting our findings with the ones from the classic afterglow model. For a constant density circumburst medium, the total number of emitted neutrinos can increase by about an order of magnitude when an optical jump occurs with respect to the self-similar afterglow scenario. By exploring the detection prospects with the IceCube Neutrino Observatory and future radio arrays such as IceCube-Gen2 radio, RNO-G and GRAND200k, as well as the POEMMA spacecraft, we conclude that the detection of neutrinos with IceCube-Gen2 radio could enable us to constrain the fraction of GRB afterglows with a jump as well as the properties of the circumburst medium. We also investigate the neutrino signal expected for the afterglows of GRB 100621A and a GRB 130427A-like burst with an optical jump. The detection of neutrinos from GRB afterglows could be crucial to explore the yet-to-be unveiled mechanism powering the optical jumps.

Early evolution of a newborn magnetar with strong precession motion in GRB 180620A

ABSTRACT The observed early X-ray plateau in the afterglow lightcurves of some gamma-ray bursts (GRBs) is attributed to the dipole radiations (DRs) of a newborn magnetar. A quasi-periodic oscillation (QPO) signal in the plateau would be strong evidence of the magnetar precession motion. By making a time-frequency domain analysis for the X-ray afterglow lightcurve of GRB 180620A, we find a QPO signal of ∼650 s in its early X-ray plateau. We fit the lightcurve with a magnetar precession model by adopting the Markov chain Monte Carlo algorithm. The observed lightcurve and the QPO signal are well represented with our model. The derived magnetic field strength of the magnetar is $B_{\rm p}= (1.02^{+0.59}_{-0.61})\times 10^{15}$ G. It rapidly spins down with angular velocity evolving as Ωs ∝ (1 + t/τsd)−0.96, where τsd = 9430 s. Its precession velocity evolution is even faster than Ωs, i.e. Ωp ∝ (1 + t/τp)−2.18 ± 0.11, where τp = 2239 ± 206 s. The inferred braking index is n = 2.04. We argue that the extra energy loss via the magnetospheric processes results in its rapid spin-down, low braking index, and strong precession motion of the magnetar.