Burst Light Curves Research Articles

Experimentally well-constrained masses of 27P and 27S: Implications for studies of explosive binary systems

The mass of 27P is expected to impact the X-ray burst (XRB) model predictions of burst light curves and the composition of the burst ashes, but large uncertainties and inconsistencies still exist in the reported 27P masses. We have used the β-decay spectroscopy of 27S to determine the most precise mass excess of 27P to date to be −659(9) keV, which is 63 keV (2.3σ) higher and a factor of 3 more precise than the value recommended in the 2016 Atomic Mass Evaluation. Based on the new 27P mass, the Si26(p,γ)27P reaction rate and its uncertainty were recalculated using Monte Carlo techniques. We also estimated the previously unknown mass excess of 27S to be 17678(77) keV, based on the measured β-delayed two-proton energy and the Coulomb displacement energy relations. The impact of these well-constrained masses and reaction rates on the modeling of the explosive astrophysical scenarios has been investigated by post-processing XRB and hydrodynamic nova models. Compared to the model calculations based on the masses and rates from databases, the abundance of A=26 in the burst ashes is increased by a factor of 2.4, while no substantial change was found in the XRB energy generation rate or the light curve. Our calculation also suggests that 27S is not a significant waiting point in the rapid proton capture process, and the change of the Si26(p,γ)27P reaction rate is not sufficiently large to affect the conclusion previously drawn on the nova contribution to the synthesis of galactic 26Al.

Spallation-altered Accreted Compositions for X-Ray Bursts: Impact on Ignition Conditions and Burst Ashes

Abstract Model predictions of X-ray burst ashes and light curves depend on the composition of the material accreted from the companion star, in particular the abundance of CNO elements. It has previously been pointed out that spallation in the atmosphere of the accreting neutron star can destroy heavy elements efficiently. In this work we study this spallation using a realistic reaction network that follows the complete spallation cascade and takes into account not only destruction, but also production of elements by the spallation of heavier species. We find an increased survival probability of heavier elements compared to previous studies, resulting in significantly higher CNO abundances. We provide resulting compositions as a function of accretion rate, and explore their impact on 1D multi-zone X-ray burst models. We find significant changes in the composition of the burst ashes, which will affect the thermal and compositional structure of accreted neutron star crusts.

Transient detection capabilities of small satellite gamma‐ray detectors

AbstractThe new, small satellite‐based gamma‐ray detectors, like Cubesats Applied for MEasuring and Localizing Transients, will provide a new way to detect gamma transients in the multimessenger era. The efficiency and the detection capabilities of such a system will be compared with current missions, for example, Fermi Gamma‐ray Burst Monitor (GBM). We used the Fermi GBM's observed short gamma‐ray burst light curves aggregated from observed discrete detector event for the simulation input. The corresponding direction‐dependent detector response matrices were used to generate photon counts and light curves around a simulated event, enabling to determine the statistics. This method can be used in the future for trigger algorithm and detector system development, and also to estimate the efficiency of the data analysis pipeline regarding the observable gamma‐ray bursts' parameters as well as other electromagnetic transients.

Consistent Modeling of GS 1826-24 X-Ray Bursts for Multiple Accretion Rates Demonstrates the Possibility of Constraining rp-process Reaction Rates

Abstract Type-I X-ray burst light curves encode unique information about the structure of accreting neutron stars and the nuclear reaction rates of the rp-process that powers bursts. Using the first model calculations of hydrogen/helium-burning bursts for a large range of astrophysical conditions performed with the code MESA, this work shows that simultaneous model–observation comparisons for bursts from several accretion rates are required to remove degeneracies in astrophysical conditions that otherwise reproduce bursts for a single and that such consistent multi-epoch modeling could possibly limit the 15O(α, γ)19Ne reaction rate. Comparisons to the 1998, 2000, and 2007 bursting epochs of the neutron star GS 1826-24 show that must be larger than previously inferred and that the shallow heating in this source must be below 0.5 MeV/u, providing a new method to constrain the shallow heating mechanism in the outer layers of accreting neutron stars. Features of the light curve rise are used to demonstrate that a lower limit could likely be placed on the 15O(α, γ) reaction rate, demonstrating the possibility of constraining nuclear reaction rates with X-ray burst light curves.

Short GRBs: Opening Angles, Local Neutron Star Merger Rate, and Off-axis Events for GRB/GW Association

Abstract The jet breaks in the afterglow light curves of short gamma-ray bursts (SGRBs), rarely detected so far, are crucial for estimating the half-opening angles of the ejecta (θ j) and hence the neutron star merger rate. In this work, we report the detection of jet decline behaviors in GRB 150424A and GRB 160821B, and find θ j ∼ 0.1 rad. Together with five events reported before 2015 and three others “identified” recently (GRB 050709, GRB 060614, and GRB 140903A), we have a sample consisting of nine SGRBs and one long-short GRB with reasonably estimated θ j. In particular, three Swift bursts in the sample have redshifts z ≤ 0.2, with which we estimate the local neutron star merger rate density to be or if the narrowly beamed GRB 061201 is excluded. Inspired by the typical θ j ∼ 0.1 rad found currently, we further investigate whether the off-beam GRBs (in the uniform jet model) or the off-axis events (in the structured jet model) can significantly enhance the GRB/GW association. For the former, the enhancement is at most moderate, while for the latter the enhancement can be much greater and a high GRB/GW association probability of ∼10% is possible. We also show that the data of GRB 160821B may contain a macronova/kilonova emission component with a temperature of ∼3100 K at ∼3.6 days after the burst and more data are needed to ultimately clarify.

The Signature of the Blandford-Znajek Mechanism in GRB Light Curves

In 1977, Blandford and Znajek showed how the spin energy of a rotating black hole may be extracted electromagnetically through a magnetic field that threads the black hole horizon. A characteristic feature of this mechanism is that, under certain fairly general conditions, the energy loss rate decays exponentially. We looked precisely for such behavior in the X-ray light curves of Long and Ultra Long duration Gamma-Ray Bursts (GRBs) observed with the XRT instrument on board the Swift satellite, and found that almost 30% of XRT light curves show an exponential decay before they reach the afterglow plateau. A similar behavior (Fast Rise Exponential Decay-FRED) was observed in γ -rays with the BATSE instrument aboard the CGRO satellite. We consider both of these findings as the signature of the Blandford-Znajek mechanism in action in the central engine of GRBs.

Millihertz quasi-periodic oscillations in 4U 1636−53 associated with bursts with positive convexity only

We investigated the convexity of all type I X-ray bursts with millihertz quasi-periodic oscillations (mHz QPOs) in 4U 1636-53 using archival observations with the Rossi X-ray Timing Explorer. We found that, at a 3.5 $\sigma$ confidence level, in all 39 cases in which the mHz QPOs disappeared at the time of an X-ray burst, the convexity of the burst is positive. The convexity measures the shape of the rising part of the burst light curve and, according to recent models, it is related to the ignition site of bursts on the neutron star surface. This finding suggests that in 4U 1636-53 these 39 bursts and the marginally-stable nuclear burning process responsible for the mHz QPOs take place at the neutron-star equator. This scenario could explain the inconsistency between the high accretion rate required for triggering mHz QPOs in theoretical models and the relatively low accretion rate derived from observations.

THREE-DIMENSIONAL SIMULATIONS OF LONG DURATION GAMMA-RAY BURST JETS: TIMESCALES FROM VARIABLE ENGINES

ABSTRACT Gamma-ray burst (GRB) light curves are characterized by marked variability, each showing unique properties. The origin of this variability, at least for a fraction of long GRBs, may be the result of an unsteady central engine. It is thus important to study the effects that an episodic central engine has on the jet propagation and, eventually, on the prompt emission within the collapsar scenario. Thus, in this study we follow the interaction of pulsed outflows with their progenitor stars with hydrodynamic numerical simulations in both two and three dimensions. We show that the propagation of unsteady jets is affected by the interaction with the progenitor material well after the break-out time, especially for jets with long quiescent times comparable to or larger than a second. We also show that this interaction can lead to an asymmetric behavior in which pulse durations and quiescent periods are systematically different. After the pulsed jets drill through the progenitor and the interstellar medium, we find that, on average, the quiescent epochs last longer than the pulses (even in simulations with symmetrical active and quiescent engine times). This could explain the asymmetry detected in the light curves of long quiescent time GRBs.

Primordial Black Holes: Observational characteristics of the final evaporation

Many early universe theories predict the creation of Primordial Black Holes (PBHs). PBHs could have masses ranging from the Planck mass to 105 solar masses or higher depending on the size of the universe at formation. A Black Hole (BH) has a Hawking temperature which is inversely proportional to its mass. Hence a sufficiently small BH will quasi-thermally radiate particles at an ever-increasing rate as emission lowers its mass and raises its temperature. The final moments of this evaporation phase should be explosive and its description is dependent on the particle physics model. In this work we investigate the final few seconds of BH evaporation, using the Standard Model and incorporating the most recent Large Hadron Collider (LHC) results, and provide a new parameterization for the instantaneous emission spectrum. We calculate for the first time energy-dependent PBH burst light curves in the GeV/TeV energy range. Moreover, we explore PBH burst search methods and potential observational PBH burst signatures. We have found a unique signature in the PBH burst light curves that may be detectable by GeV/TeV gamma-ray observatories such as the High Altitude Water Cerenkov (HAWC) observatory. The implications of beyond the Standard Model theories on the PBH burst observational characteristics are also discussed, including potential sensitivity of the instantaneous photon detection rate to a squark threshold in the 5–10 TeV range.

CSI 2264: CHARACTERIZING YOUNG STARS IN NGC 2264 WITH STOCHASTICALLY VARYING LIGHT CURVES*

ABSTRACT We provide CoRoT and Spitzer light curves and other supporting data for 17 classical T Tauri stars in NGC 2264 whose CoRoT light curves exemplify the “stochastic” light curve class as defined in 2014 by Cody et al. The most probable physical mechanism to explain the optical variability within this light curve class is time-dependent mass accretion onto the stellar photosphere, producing transient hot spots. Where we have appropriate spectral data, we show that the veiling variability in these stars is consistent in both amplitude and timescale with the optical light curve morphology. The veiling variability is also well-correlated with the strength of the He i 6678 Å emission line, predicted by models to arise in accretion shocks on or near the stellar photosphere. Stars with accretion burst light curve morphology also have variable mass accretion. The stochastic and accretion burst light curves can both be explained by a simple model of randomly occurring flux bursts, with the stochastic light curve class having a higher frequency of lower amplitude events. Members of the stochastic light curve class have only moderate mass accretion rates. Their Hα profiles usually have blueshifted absorption features, probably originating in a disk wind. The lack of periodic signatures in the light curves suggests that little of the variability is due to long-lived hot spots rotating into or out of our line of sight; instead, the primary driver of the observed photometric variability is likely to be instabilities in the inner disk that lead to variable mass accretion.

REACTION RATES OF 64Ge(p,γ)65As AND 65As(p,γ)66Se AND THE EXTENT OF NUCLEOSYNTHESIS IN TYPE I X-RAY BURSTS

ABSTRACT The extent of nucleosynthesis in models of type I X-ray bursts (XRBs) and the associated impact on the energy released in these explosive events are sensitive to nuclear masses and reaction rates around the 64Ge waiting point. Using the well known mass of 64Ge, the recently measured 65As mass, and large-scale shell model calculations, we have determined new thermonuclear rates of the 64Ge(p,γ)65As and 65As(p,γ)66Se reactions with reliable uncertainties. The new reaction rates differ significantly from previously published rates. Using the new data, we analyze the impact of the new rates and the remaining nuclear physics uncertainties on the 64Ge waiting point in a number of representative one-zone XRB models. We find that in contrast to previous work, when all relevant uncertainties are considered, a strong 64Ge rp-process waiting point cannot be ruled out. The nuclear physics uncertainties strongly affect XRB model predictions of the synthesis of 64Zn, the synthesis of nuclei beyond A = 64, the energy generation, and the burst light curve. We also identify key nuclear uncertainties that need to be addressed to determine the role of the 64Ge waiting point in XRBs. These include the remaining uncertainty in the 65As mass, the uncertainty of the 66Se mass, and the remaining uncertainty in the 65As(p,γ)66Se reaction rate, which mainly originates from uncertain resonance energies.

DETECTING THE SUPERNOVA BREAKOUT BURST IN TERRESTRIAL NEUTRINO DETECTORS

ABSTRACT We calculate the distance-dependent performance of a few representative terrestrial neutrino detectors in detecting and measuring the properties of the ν e breakout burst light curve in a Galactic core-collapse supernova. The breakout burst is a signature phenomenon of core collapse and offers a probe into the stellar core through collapse and bounce. We examine cases of no neutrino oscillations and oscillations due to normal and inverted neutrino-mass hierarchies. For the normal hierarchy, other neutrino flavors emitted by the supernova overwhelm the ν e signal, making a detection of the breakout burst difficult. For the inverted hierarchy (IH), some detectors at some distances should be able to see the ν e breakout burst peak and measure its properties. For the IH, the maximum luminosity of the breakout burst can be measured at 10 kpc to accuracies of ∼30% for Hyper-Kamiokande (Hyper-K) and ∼60% for the Deep Underground Neutrino Experiment (DUNE). Super-Kamiokande (Super-K) and Jiangmen Underground Neutrino Observatory (JUNO) lack the mass needed to make an accurate measurement. For the IH, the time of the maximum luminosity of the breakout burst can be measured in Hyper-K to an accuracy of ∼3 ms at 7 kpc, in DUNE to ∼2 ms at 4 kpc, and JUNO and Super-K can measure the time of maximum luminosity to an accuracy of ∼2 ms at 1 kpc. Detector backgrounds in IceCube render a measurement of the ν e breakout burst unlikely. For the IH, a measurement of the maximum luminosity of the breakout burst could be used to differentiate between nuclear equations of state.

Update on the GRB universal scaling EX,iso–Eγ,iso–Epk with 10 years of Swift data

From a comprehensive statistical analysis of $Swift$ X-ray light-curves of gamma-ray bursts (GRBs) collected from December 2004 to the end of 2010, we found a three-parameter correlation between the isotropic energy emitted in the rest frame 1-10$^4$ keV energy band during the prompt emission (E$_{\rm{\gamma,iso}}$), the rest frame peak of the prompt emission energy spectrum (E$_{\rm{pk}}$), and the X-ray energy emitted in the rest frame 0.3-30 keV observed energy band (E$_{\rm{X,iso}}$), computed excluding the contribution of the flares. In this paper, we update this correlation with the data collected until June 2014, expanding the sample size with $\sim$35% more objects, where the number of short GRBs doubled. With this larger sample we confirm the existence of a universal correlation that connects the prompt and afterglow properties of long and short GRBs. We show that this correlation does not depend on the X-ray light-curve morphology and that further analysis is necessary to firmly exclude possible biases derived by redshift measurements. In addition we discuss about the behavior of the peculiar objects as ultra-long GRBs and we propose the existence of an intermediate group between long and short GRBs. Interestingly, two GRBs with uncertain classification fall into this category. Finally, we discuss the physics underlying this correlation, in the contest of the efficiency of conversion of the prompt $\gamma$-ray emission energy into the kinetic energy of the afterglow, the photosferic model, and the cannonball model.

Effect of rapid evolution of magnetic tilt angle on a newborn magnetar's dipole radiation**Supported by the National Natural Science Foundation of China.

We study the electromagnetic radiation from a newborn magnetar whose magnetic tilt angle decreases rapidly. We calculate the evolution of the angular spin frequency, the perpendicular component of the surface magnetic field strength, and the energy loss rate through magnetic dipole radiation. We show that the spin-down of the magnetar experiences two stages characterized by two different timescales. The apparent magnetic field decreases with the decrease of the tilt angle. We further show that the energy loss rate of the magnetar is very different from that in the case of a fixed tilt angle. The evolution of the energy loss rate is consistent with the overall light curves of gamma-ray bursts which show a plateau structure in their afterglow stage. Our model supports the idea that some gamma-ray bursts with a plateau phase in their afterglow stage may originate from newborn millisecond magnetars.

Broad-band modelling of short gamma-ray bursts with energy injection from magnetar spin-down and its implications for radio detectability

The magnetar model has been proposed to explain the apparent energy injection in the X-ray light curves of short gamma-ray bursts (SGRBs), but its implications across the full broadband spectrum are not well explored. We investigate the broadband modelling of four SGRBs with evidence for energy injection in their X-ray light curves, applying a physically motivated model in which a newly formed magnetar injects energy into a forward shock as it loses angular momentum along open field lines. By performing an order of magnitude search for the underlying physical parameters in the blast wave, we constrain the characteristic break frequencies of the synchrotron spectrum against their manifestations in the available multi-wavelength observations for each burst. The application of the magnetar energy injection profile restricts the successful matches to a limited family of models that are self-consistent within the magnetic dipole spin-down framework.We produce synthetic light curves that describe how the radio signatures of these SGRBs ought to have looked given the restrictions imposed by the available data, and discuss the detectability of these signatures with present-day and near-future radio telescopes. Our results show that both the Atacama Large Millimetre Array and the upgraded Very Large Array are now sensitive enough to detect the radio signature within two weeks of trigger in most SGRBs, assuming our sample is representative of the population as a whole. We also find that the upcoming Square Kilometre Array will be sensitive to depths greater than those of our lower limit predictions.

MEPSA: A flexible peak search algorithm designed for uniformly spaced time series

Abstract We present a novel algorithm aimed at identifying peaks within a uniformly sampled time series affected by uncorrelated Gaussian noise. The algorithm, called “MEPSA” (multiple excess peak search algorithm), essentially scans the time series at different timescales by comparing a given peak candidate with a variable number of adjacent bins. While this has originally been conceived for the analysis of gamma-ray burst light (GRB) curves, its usage can be readily extended to other astrophysical transient phenomena, whose activity is recorded through different surveys. We tested and validated it through simulated featureless profiles as well as simulated GRB time profiles. We showcase the algorithm’s potential by comparing with the popular algorithm by Li and Fenimore, that is frequently adopted in the literature. Thanks to its high flexibility, the mask of excess patterns used by MEPSA can be tailored and optimised to the kind of data to be analysed without modifying the code. The C code is made publicly available.

β-delayed proton emission of 69Kr and the 68Se rp-process waiting point

The slow β-decay of the 68Se waiting point in the astrophysical rp-process can in principle be bypassed by a sequential two proton capture. The rate of this reaction sequence depends exponentially on the 69Br proton separation energy. We studied β-delayed proton emission of 69Kr and extracted a proton separation energy of −641(42) keV. In addition, we determined a 69Kr β-decay half-life of T1/2=28(1) ms and an excitation energy of 3153(55) keV of the 69Kr isobaric analog state in 69Br. X-ray burst model calculations show that regardless of the values of other uncertain masses in the region, the new Sp(Br69) allows for a reaction flow via Se68(2p,γ) of at most 20%. Uncertainties are sufficiently reduced to conclude that Se68(2p,γ) has at best a very small effect on burst light curve and composition, and that 68Se is a strong rp-process waiting point. Our results also exclude the possibility of a suggested longer lived, so far unobserved, 69Br ground state.

Investigating signatures of cosmological time dilation in duration measures of prompt gamma-ray burst light curves

We study the evolution with redshift of three measures of gamma-ray burst (GRB) duration ($T_{\rm 90}$, $T_{\rm 50}$ and $T_{\rm R45}$) in a fixed rest frame energy band for a sample of 232 Swift/BAT detected GRBs. Binning the data in redshift we demonstrate a trend of increasing duration with increasing redshift that can be modelled with a power-law for all three measures. Comparing redshift defined subsets of rest-frame duration reveals that the observed distributions of these durations are broadly consistent with cosmological time dilation. To ascertain if this is an instrumental effect, a similar analysis of Fermi/GBM data for the 57 bursts detected by both instruments is conducted, but inconclusive due to small number statistics. We then investigate under-populated regions of the duration redshift parameter space. We propose that the lack of low-redshift, long duration GRBs is a physical effect due to the sample being volume limited at such redshifts. However, we also find that the high-redshift, short duration region of parameter space suffers from censorship as any Swift GRB sample is fundamentally defined by trigger criteria determined in the observer frame energy band of Swift/BAT. As a result, we find that the significance of any evidence for cosmological time dilation in our sample of duration measures typically reduces to $<2\sigma$.

GAMMA-RAY BURST FLARES: X-RAY FLARING. II.

We present a catalog of 497 flaring periods found in gamma-ray burst (GRB) light curves taken from the online XRT GRB Catalogue. We analyzed 680 individual light curves using a flare detection method developed and used on our UV/optical GRB Flare Catalog. The method makes use of the Bayesian Information Criterion to analyze the residuals of fitted GRB light curves and statistically determines the optimal fit to the light curve residuals in attempt to identify any additional features. These features, which we classify as flares, are identified by iteratively adding additional `breaks' to the light curve. We find evidence of flaring in 310 of the analyzed light curves. For those light curves with flares, we find an average number of ~1.5 flares per GRB. As with the UV/optical, flaring in our sample is generally confined to the first 1000 s of the afterglow, but can be detected to beyond 10^5 s. Only ~50% of the detected flares follow the `classical' definition of \Delta t/t << 1, with many of the largest flares exhibiting \Delta t/t ~ 1 and some exceeding this value.

UNCOVERING THE INTRINSIC VARIABILITY OF GAMMA-RAY BURSTS

We develop a robust technique to determine the minimum variability timescale for Gamma-ray Burst light curves, utilizing Haar wavelets. Our approach averages over the data for a given GRB, providing an aggregate measure of signal variation while also retaining sensitivity to narrow pulses within complicated time-series. In contrast to previous studies using wavelets, which simply define the minimum timescale in reference to the measurement noise floor, our approach identifies the signature of temporally-smooth features in the wavelet scaleogram and then additionally identifies a break in the scaleogram on longer timescales as signature of a true, temporally-unsmooth light curve feature or features. We apply our technique to the large sample of Swift GRB Gamma-ray light curves and for the first time -- due to the presence of a large number of GRBs with measured redshift -- determine the distribution of minimum variability timescales in the source frame. We find a median minimum timescale for long-duration GRBs in the source frame of $\Delta t_{\rm min}=0.5$ s, with the shortest timescale found being on the order of 10 ms. This short timescale suggests a compact central engine ($3 \times 10^3$ km). We discuss further implications for the GRB fireball model and present a tantalizing correlation between minimum timescale and redshift, which may in part be due to cosmological time-dilation.