Light Curves Of Gamma-ray Bursts Research Articles

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.

One Fits All: A Unified Synchrotron Model Explains GRBs with FRED-shape Pulses

The analysis of gamma-ray burst (GRB) spectra often relies on empirical models lacking a distinct physical explanation. Previous attempts to couple physical models with observed data focus on individual burst studies, fitting models to segmented spectra with independent physical parameters. However, these approaches typically neglect to explain the time evolution of observed spectra. In this study, we propose a novel approach by incorporating the synchrotron radiation model to provide a self-consistent explanation for a selection of single-pulse GRBs. Our study comprehensively tests the synchrotron model under a unified physical condition, such as a single injection event of electrons. By tracing the evolution of cooling electrons in a decaying magnetic field, our model predicts time-dependent observed spectra that align well with the data. Using a single set of physical parameters, our model successfully fits all time-resolved spectra within each burst. Our model suggests that the rising phase of the GRB light curve results from the increasing number of radiating electrons, while the declining phase is attributed to the curvature effect, electron cooling, and the decaying magnetic field. Our model provides a straightforward interpretation of the peak energy’s evolution, linked to the decline of the magnetic field and electron cooling due to the expansion of the GRB emission region. Our findings strongly support the notion that spectral and temporal evolution in GRB pulses originates from the expansion of the GRB emission region, with an initial radius of approximately 1015 cm, and synchrotron radiation as the underlying emission mechanism.

What absorbs the early TeV photons of GRB 221009A?

ABSTRACT The tera-electronvolt (TeV) light curve of gamma-ray burst 221009A exhibits an unprecedentedly rapid rise at its start. This phenomenon could be due to the strong absorption of photons and electrons within the emitting region. As the external shock expands outward and the radius increases, the volume of matter also grows, resulting in a gradual reduction of the optical depth for TeV photons. We investigate several potential explanations for the early TeV light curves. We calculate the optical depth for TeV photons, considering both annihilation with lower energy photons in the external shock and their scattering by electrons generated through the cascading of the TeV emission. Even under the favourable assumptions, we have determined that the optical depths for these processes are orders of magnitude too small to explain the observed light curve. Additional sources of absorbers, including electrons in the ejecta or external shock, also do not result in sufficient optical depths. Hence, the cause behind the early peculiar TeV light curve remains unclear.

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.

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.

The signature of refreshed shocks in the afterglow of GRB 030329

ABSTRACT GRB 030329 displays one clear and, possibly, multiple less intense fast-rising (Δt/t ∼ 0.3) jumps in its optical afterglow light curve. The decay rate of the optical light curve remains the same before and after the flux jumps. This may be the signature of energy injection into the shocked material at the front of the jet. In this study, we model the Gamma-ray Burst (GRB) ejecta as a series of shells. We follow the dynamical evolution of the ejecta as it interacts with itself (i.e. internal shocks) and with the circumburst medium (i.e. external forward and reverse shocks), and calculate the emission from each shock event assuming synchrotron emission. We confirm the viability of the proposed model in which the jumps in the optical afterglow light curve of GRB 030329 are produced via refreshed shocks. The refreshed shocks may be the signatures of collisions between earlier ejected material with an average Lorentz factor $\bar{\Gamma }\gtrsim 100$ and later ejected material with $\bar{\Gamma } \sim 10$ once the early material has decelerated due to interaction with the circumburst medium. We show that even if the late material is ejected with a spread of Lorentz factors, internal shocks naturally produce a narrow distribution of Lorentz factors (ΔΓ/Γ ≲ 0.1), which is a necessary condition to produce the observed quick rise times of the jumps. These results imply a phase of internal shocks at some point in the dynamical evolution of the ejecta, which requires a low magnetization in the outflow.

A Stochastic Approach to Reconstruct Gamma-Ray-burst Light Curves

Gamma-ray bursts (GRBs), as they are observed at high redshift (z = 9.4), are vital to cosmological studies and investigating Population III stars. To tackle these studies, we need correlations among relevant GRB variables with the requirement of small uncertainties on their variables. Thus, we must have good coverage of GRB light curves (LCs). However, gaps in the LC hinder the precise determination of GRB properties and are often unavoidable. Therefore, extensive categorization of GRB LCs remains a hurdle. We address LC gaps using a stochastic reconstruction, wherein we fit two preexisting models (the Willingale model; W07; and a broken power law; BPL) to the observed LC, then use the distribution of flux residuals from the original data to generate data to fill in the temporal gaps. We also demonstrate a model-independent LC reconstruction via Gaussian processes. At 10% noise, the uncertainty of the end time of the plateau, its correspondent flux, and the temporal decay index after the plateau decreases by 33.3%, 35.03%, and 43.32% on average for the W07, and by 33.3%, 30.78%, 43.9% for the BPL, respectively. The uncertainty of the slope of the plateau decreases by 14.76% in the BPL. After using the Gaussian process technique, we see similar trends of a decrease in uncertainty for all model parameters for both the W07 and BPL models. These improvements are essential for the application of GRBs as standard candles in cosmology, for the investigation of theoretical models, and for inferring the redshift of GRBs with future machine-learning analyses.

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.

Afterglow Light Curves from Off-Axis GRB Jets in Stratified Circumburst Medium

We study the gamma-ray burst (GRB) afterglow light curves produced by an off-axis jet in a stratified circumburst medium and summarize the temporal indices of the coasting phase, the deceleration phase, the Newtonian phase, and the deep Newtonian phase for various viewing angles and power-law indices of medium density. Generally, the afterglow light curves of off-axis GRBs in the homogeneous interstellar medium have a steep rise arising due to jet deceleration. In the stratified medium, the flux rises is more shallow but peaks earlier for the same viewing angle due to faster deceleration of the jet running into the denser stratified medium, compared with the case of the interstellar medium (ISM). Observations of off-axis bursts will possibly increase over the coming years due to the arrival of the multi-messenger era and the forthcoming surveys in multiple bands. The temporal indices of off-axis afterglows derived in the paper will provide a reference for comparison with the observations and can diagnose the circumburst environment. The numerical code calculating the afterglow light curve also can be used to fit the multi-wavelength light curves.

The Synchrotron Low-energy Spectrum Arising from the Cooling of Electrons in Gamma-Ray Bursts

This work is a continuation of a previous effort (Panaitescu) to study the cooling of relativistic electrons through radiation (synchrotron and self-Compton) emission and adiabatic losses, with application to the spectra and light curves of the synchrotron gamma-ray burst (GRB) produced by such cooling electrons. Here, we derive the low-energy slope β LE of a GRB pulse-integrated spectrum and quantify the implications of the measured distribution of β LE. Radiative processes that produce soft integrated spectra can accommodate the harder slopes measured by CGRO/BATSE and Fermi/GBM only if the magnetic field lifetime t B is shorter than the time during which the typical GRB electrons cool to radiate below 1–10 keV, which is less than (at most) 10 radiative cooling timescales t rad of the typical GRB electron. In this case, there is a one-to-one correspondence between t B and β LE. To account for low-energy slopes β LE > −3/4, the adiabatic electron-cooling requires a similar restriction on t B . In this case, the diversity of slopes arises mostly from how the electron-injection rate varies with time (temporal power-law injection rates yield power-law low-energy GRB spectra) and not from the magnetic field timescale.

Comprehension of Jet Physics from the Analysis of the Gamma‐Ray Burst Afterglow Parameter Distributions

AbstractThe discovery of a short gamma‐ray burst (GRB), GRB170817A, in association with a gravitational wave (GW) and a bright kilonova started a new era in the high‐energy astrophysics. The observation of GRB170817A and more recently, GRB200826A and GRB211211A, a short and a long burst, respectively, with a possible kilonova, reinforce the concern about new ways of classification. For this reason, a new machine learning technique is applied to Swift‐BAT data, searching for morphological similarities in the light curves. The resulting map is characterized by two distinct groups, although still correlated with standard T90 duration. Since a jet viewed off‐axis could explain the emission from GRB170817A, the modeling of this kind of sources is of great importance. A public code called JetFit, based on the “boosted fireball” model, is applied to fit Swift‐XRT afterglow light curves of short and long GRBs, with known red‐shift, from 2005 to 2021. JetFit does not model the flaring activity. For this purpose, a new technique to remove the time flaring phases, is developed. This analysis provides a comprehensive study of the prompt and of the afterglow phase, trough the study of the best‐fit parameters.

On the Jet Structures of GRB 050820A and GRB 070125

We present the broadband numerical modeling of afterglows for two remarkably bright long gamma-ray bursts (GRBs), GRB 050820A and GRB 070125, with a wide range of observations from the radio band to the X-ray band. In our work, we fit light curves and constrain physical parameters using a standard forward shock model from the afterglowpy Python package, considering different jet structures and the jet lateral expansion. For GRB 050820A, the constrained jet is close to a top-hat jet with an extremely small half opening angle of about 0.015 rad, and the circumburst matter density is as small as 10−7 cm−3, which suggests that this peculiar long GRB might originate from metal-poor stars with low mass-loss rates. To explain the late time optical light curves of GRB 070125, the effects of the lateral expansion and the participation factor of electrons that are accelerated by the shock have to be taken into account. The constrained results for GRB 070125 show that the jet is also close to a top-hat jet with a half opening angle of about 0.1 rad, the viewing angle is about 0.05 rad, the circumburst density is about 10 cm−3, and the participation factor is about 0.1. The jet energy of the two bursts is required to be ∼1051–1052 erg, which can be produced by a millisecond magnetar or a hyper-accreting black hole.

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.

Gamma-Ray Burst Prompt Emission Spectrum and E p Evolution Patterns in the ICMART Model

Abstract In this paper, we simulate the prompt emission light curve, spectrum, and E p evolution patterns of gamma-ray bursts (GRBs) within the framework of the Internal-collision-induced Magnetic Reconnection and Turbulence (ICMART) model. We show that this model can produce a Band-shaped spectrum, whose parameters (E p, α, β) could follow the typical distribution of GRB observations, as long as the magnetic field and the electron acceleration process in the emission region are under appropriate conditions. On the other hand, we show that for one ICMART event, E p evolution is always a hard-to-soft pattern. However, a GRB light curve is usually composed of multiple ICMART events that are fundamentally driven by the erratic activity of the GRB central engine. In this case, we find that if one individual broad pulse in the GRB light curve is composed of multiple ICMART events, the overall E p evolution could be disguised as an intensity-tracking pattern. Therefore, mixed E p evolution patterns can coexist in the same burst, with a variety of combined patterns. Our results support the ICMART model as a competitive model to explain the main properties of GRB prompt emission. The possible challenges faced by the ICMART model are also discussed in detail.

Instrumental Tip-of-the-iceberg Effects on the Prompt Emission of Swift/BAT Gamma-ray Bursts

Abstract The observed durations of prompt gamma-ray emission from gamma-ray bursts (GRBs) are often used to infer the progenitors and energetics of the sources. Inaccurate duration measurements will have a significant impact on constraining the processes powering the bursts. The “tip-of-the-iceberg” effect describes how the observed signal is lost into background noise; lower instrument sensitivity leads to higher measurement bias. In this study, we investigate how observing conditions, such as the number of enabled detectors, background level, and incident angle of the source relative to the detector plane, affect the measured duration of GRB prompt emission observed with the Burst Alert Telescope on board the Neil Gehrels Swift Observatory (Swift/BAT). We generate “simple-pulse” light curves from an analytical fast rise exponential decay function and from a sample of eight real GRB light curves. We fold these through the Swift/BAT instrument response function to simulate light curves Swift/BAT would have observed for specific observing conditions. We find duration measurements are highly sensitive to observing conditions and the incident angle of the source has the highest impact on measurement bias. In most cases duration measurements of synthetic light curves are significantly shorter than the true burst duration. For the majority of our sample, the percentage of duration measurements consistent with the true duration is as low as ∼25%–45%. In this article, we provide quantification of the tip-of-the-iceberg effect on GRB light curves due to Swift/BAT instrumental effects for several unique light curves.

GRB 191016A: The onset of the forward shock and evidence of late energy injection

ABSTRACTWe present optical and near-infrared (NIR) photometric observations of GRB 191016 with the COATLI,DDOTI, and RATIR ground-based telescopes over the first three nights. We present the temporal evolution of the optical afterglow and describe five different stages that were not completely characterized in previous works, mainly due to scarcity of data points to accurately fit the different components of the optical emission. After the end of the prompt gamma-ray emission, we observed the afterglow rise slowly in the optical and NIR wavelengths and peak at around T + 1450 s in all filters. This was followed by an early decay, a clear plateau from T + 5000 s to T + 11 000 s, and then a regular late decay. We also present evidence of the jet break at later times, with a temporal index in good agreement with the temporal slope obtained from X-ray observations. Although many of the features observed in the optical light curves of gamma-ray bursts are usually well explained by a reverse shock (RS) or forward shock (FS), the shallowness of the optical rise and enhanced peak emission in the GRB 191016A afterglow is not well fitted by only a FS or a RS. We propose a theoretical model which considers both of these components and combines an evolving FS with a later embedded RS and a subsequent late energy injection from the central engine activity. We use this model to successfully explain the temporal evolution of the light curves and discuss its implications on the fireball properties.

On the Emergence of a New Instability during Core Collapse of Very Massive Stars

Abstract The process of uniform supernovae (SNe) explosions is well investigated for all their types. However, observational data suggests that the SNe may be not spherically symmetric. Modern multidimensional simulations of SNe demonstrate development of hydrodynamical instabilities during the explosion phase. But the configuration of a star and inhomogeneities prior to explosion could strongly affect how the SN develops. A number of papers on numerical modeling of pair-instability SNe explosion considered the case when thermonuclear energy in the central region of a massive star is introduced by a series of several hot spots. It leads to the appearance of many fragments of hot matter behind the divergence shock wave. An observable manifestation of this may be the presence of peaks on light curves of gamma-ray bursts associated with explosions of massive stars. The physical nature of such inhomogeneities is not evident and the number and size of spots is conjecture. In this work, we study the possibility of formation of these inhomogeneities at the stage of core collapse (CC) in a massive star. To check this assumption, we chose an analytic self-similar model of CC and investigated the stability of solutions obtained from it with respect to small multidimensional perturbations. It shows there are no conditions where the collapse of a very massive star may remain stable, although, for a less massive star, it is possible. Using the relations obtained, we found characteristic features of developing instability, thereby making it possible to estimate the amount and characteristic size of the inhomogeneities.