X-ray Bursts Research Articles

Monster Radiative Shocks in the Perturbed Magnetospheres of Neutron Stars

Magnetospheres of neutron stars can be perturbed by star quakes, interaction in a binary system, or sudden collapse of the star. The perturbations are typically in the kilohertz band and excite magnetohydrodynamic waves. We show that compressive magnetospheric waves steepen into monster shocks, possibly the strongest shocks in the Universe. The shocks are radiative, i.e., the plasma energy is radiated before it crosses the shock. As the kilohertz wave with the radiative shock expands through the magnetosphere, it produces a bright X-ray burst. Then, it launches an approximately adiabatic blast wave, which will expand far from the neutron star. These results suggest a new mechanism for X-ray bursts from magnetars and support the connection of magnetar X-ray activity with fast radio bursts. Similar shocks may occur in magnetized neutron-star binaries before they merge, generating an X-ray precursor of the merger. Powerful radiative shocks are also predicted in the magnetosphere of a neutron star when it collapses into a black hole, producing a bright X-ray transient.

Quantifying chaos and randomness in magnetar bursts

ABSTRACT In this study, we explore the dynamical stability of magnetar bursts within the context of the chaos-randomness phase space for the first time, aiming to uncover unique behaviours compared to various astrophysical transients, including fast radio bursts (FRBs). We analyse burst energy time series data from active magnetar sources SGR J1550−5418 and SGR J1935+2154, focusing on burst arrival time and energy differences between consecutive events. We find a distinct separation in the time domain, where magnetar bursts exhibit significantly lower randomness compared to FRBs, solar flares, and earthquakes, with a slightly higher degree of chaos. In the energy domain, magnetar bursts exhibit a broad consistency with other phenomena, primarily due to the wide distribution of chaos-randomness observed across different bursts and sources. Intriguingly, contrary to expectations from the FRB-magnetar connection, the arrival time patterns of magnetar bursts in our analysis do not exhibit significant proximity to repeating FRBs in the chaos-randomness plane. This finding may challenge the hypothesis that FRBs are associated with typical magnetar bursts but indirectly supports the evidence that FRBs may primarily be linked to special magnetar bursts like peculiar X-ray bursts from SGR J1935+2154 observed simultaneously with Galactic FRB 200428.

NICER Observations of Thermonuclear Bursts from 4U 1728-34: Detection of Oscillations prior to the Onset of Two Bursts

We present temporal and time-resolved spectral analyses of all the thermonuclear X-ray bursts observed from the neutron star low-mass X-ray binary 4U 1728−34 with NICER from 2017 June to 2019 September. In total, we detected 11 X-ray bursts from the source and performed time-resolved spectroscopy. Unlike some of the earlier results for other bursting sources from NICER, our spectral results indicate that the use of a scaling factor for the persistent emission is not statistically necessary. This is primarily a result of the strong interstellar absorption in the line of sight toward 4U 1728−34, which causes the count rates to be significantly lower at low energies. We also searched for burst oscillations and detected modulations in six different bursts at around the previously known burst oscillation frequency of 363 Hz. Finally, we report the detection of oscillations prior to two bursts at 356 and 359 Hz, respectively. This is the first time in the literature where burst oscillations are detected before the rapid rise in X-ray flux, from any known burster. These oscillations disappear as soon as the burst starts to rise and occur at a somewhat lower frequency than the oscillations we detect during the bursts.

Insight-HXMT observations on thermonuclear X-ray bursts from 4U 1608–52 in 2022: The accretion rate dependent anisotropy of burst emission

Thermonuclear X-ray bursts occur on the surface of an accreting neutron star (NS), and their characteristics and interplay with the surrounding circumstance could be a clue to understand the nature of the NS and accretion process. For this purpose, Insight-HXMT has performed high cadence observations on the bright thermonuclear X-ray burster–4U 1608–52 during its outburst in July and August 2022; nine bursts were detected, including seven bursts with the photospheric radius expansion (PRE). Time-resolved spectroscopy of the bright PRE bursts reveals that an enhancement of accretion rate or the Comptonization of the burst emission by the corona could reduce the residuals when fitting their spectra with the conventional model–blackbody. The inferred energy increment rate of the burst photon gained from the corona is up to ∼40%, even though the bursts have different peak fluxes and locate at different accretion rates. Moreover, the flux shortage of the rising PRE is observed in the bursts at a high mass accretion rate, but not for the burst with a faint persistent emission, which has been predicted theoretically but first observed in this work. If the flux shortage is due to the disk obscuration, i.e., the burst emission is anisotropic, the phenomenon above could indicate that the anisotropy of the burst emission is accretion rate dependent, which could also be evidence of the truncated disk in the low/hard state.

Observing white dwarf tidal stripping with TianQin gravitational wave observatory

ABSTRACT Recently discovered regular X-ray bursts known as quasi-periodic eruptions have a proposed model that suggests a tidal stripping white dwarf inspiralling into the galaxy’s central black hole on an eccentric orbit. According to this model, the interaction of the stripping white dwarf with the central black hole would also emit gravitational wave signals, their detection can help explore the formation mechanism of quasi-periodic eruptions and facilitate multimessenger observations. In this paper, we investigated the horizon distance of TianQin on this type of gravitation wave signal and found it can be set to 200 Mpc. We also find that those stripping white dwarf model sources with central black hole mass within $10^4 \!-\! 10^{5.5}\, \mathrm{M}_\odot$ are more likely to be detected by TianQin. We assessed the parameter estimation precision of TianQin on those stripping white dwarf model sources. Our result shows that, even in the worst case, TianQin can determine the central black hole mass, the white dwarf mass, the central black hole spin, and the orbital initial eccentricity with a precision of 10−2. In the optimistic case, TianQin can determine the central black hole mass and the white dwarf mass with a precision of 10−7, determine the central black hole spin with a precision of 10−5, and determine the orbital initial eccentricity with a precision of 10−8. Moreover, TianQin can determine the luminosity distance with a precision of 10−1 and determine the sky localization with a precision of 10−2–10 $\rm deg^2$.

AstroSat observation of the magnetar SGR J1830−0645 during its first detected X-ray outburst

ABSTRACT We present here timing and spectral analyses of SGR J1830−0645 based on an AstroSat observation carried out on 2020 October 16, about a week after the onset of its first detected X-ray outburst. Using data taken with the Soft X-ray Telescope (SXT) and Large Area X-ray Proportional Counter (LAXPC), we have detected 0.9–10 keV coherent pulsations at a period of ∼10.4 s. The pulse profiles were single-peaked, asymmetric, and consisted of minor peaks attributable to hotspots on the neutron star surface. The pulsed fraction evolved significantly with energy, increasing to energies around 5 keV with a steep drop thereafter. The 0.9–25 keV SXT–LAXPC energy spectrum is best described with two thermal components having temperatures ∼0.46 and ∼1.1 keV (emission radii of ∼2.4 and ∼0.65 km, respectively, assuming a distance of 4 kpc) along with a power-law component having a photon index of ∼0.39. We report the detection of 67 X-ray bursts having an average duration of ∼33 ms. The brightest burst lasted for about 90 ms and had a 3–25 keV fluence of ∼5 × 10−9 erg cm−2.

Statistical Properties of X-Ray Bursts from SGR J1935+2154 Detected by Insight-HXMT

As one class of the most important objects in the universe, magnetars can produce a lot of different frequency bursts including X-ray bursts. In Cai et al., 75 X-ray bursts produced by magnetar SGR J1935+2154 during an active period in 2020 are published, including the duration and net photon counts of each burst, and waiting time based on the trigger time difference. In this paper, we utilize the power-law model, , to fit the cumulative distributions of these parameters. It can be found that all the cumulative distributions can be well fitted, which can be interpreted by a self-organizing criticality theory. Furthermore, we check whether this phenomenon still exists in different energy bands and find that there is no obvious evolution. These findings further confirm that the X-ray bursts from magnetars are likely to be generated by some self-organizing critical process, which can be explained by a possible magnetic reconnection scenario in magnetars.

Mass-loss and composition of wind ejecta in type I X-ray bursts

Context. X-Ray bursts (XRBs) are powerful thermonuclear events on the surface of accreting neutron stars (NSs) where nucleosynthesis of intermediate-mass elements occurs. The high surface gravity of an NS prevents the ejection of material by the XRB thermonuclear explosion. However, the predicted and observed XRB luminosities sometimes exceed Eddington’s value, and some of the material may escape by means of stellar wind. Aims. This work aims to determine the mass-loss and chemical composition of the material ejected through radiation-driven winds and its significance for Galactic abundances. It also reports on the evolution of observational quantities during the wind phase, which can help constrain the mass-radius relation in NSs. Methods. A non-relativistic radiative wind model was implemented, with modern opacity tables and treatment of the critical point, and linked through a new technique to a series of XRB hydrodynamic simulations that include over 300 isotopes. This allowed us to construct a quasi-stationary time evolution of the wind during the XRB. Results. In the models we studied, the total mass ejected by the wind was about 6 × 1019 g; the average ejected mass per burst represented 2.6% of the accreted mass between bursts, with 0.1% of the envelope mass ejected per burst; and approximately 90% of the ejecta was composed by 60Ni, 64Zn, 68Ge, and 58Ni. The ejected material also contained a small fraction (10−4 − 10−5) of some light p-nuclei, but not enough to account for their Galactic abundances. Additionally, the observable magnitudes during the wind phase showed remarkable correlations, partly due to the fact that photospheric luminosity stays close to the Eddington limit. Some of these correlations involve wind parameters, such as energy and mass outflows, that are determined by the conditions at the base of the wind envelope. Conclusions. The simulations resulted in the first realistic quantification of mass-loss for each isotope synthesized in the XRB. The photospheric correlations we found could be used to link observable magnitudes to the physics of the innermost parts of the envelope, close to its interface with the NS crust. This is a promising result regarding the issue of NS radius determination.

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.

Sensitivity of He Flames in X-Ray Bursts to Nuclear Physics

Through the use of axisymmetric 2D hydrodynamic simulations, we further investigate laterally propagating flames in X-ray bursts (XRBs). Our aim is to understand the sensitivity of a propagating helium flame to different nuclear physics. Using the Castro simulation code, we confirm the phenomenon of enhanced energy generation shortly after a flame is established by adding 12C(p, γ)13N(α, p)16O to the network, in agreement with the past literature. This sudden outburst of energy leads to a short accelerating phase, causing a drastic alteration in the overall dynamics of the flame in XRBs. Furthermore, we investigate the influence of different plasma screening routines on the propagation of the XRB flame. We finally examine the performance of simplified spectral deferred correction, a novel approach to hydrodynamics and reaction coupling incorporated in Castro, as an alternative to operator splitting.

Spectral analysis of the LMXB XTE J1810−189 with NICER data

ABSTRACT X-ray Timing Explorer J1810−189 is a low-mass X-ray binary transient system hosting a neutron star, which underwent a three-month-long outburst in 2020. In order to study its spectral evolution during this outburst, we analysed all the available observations performed by Neutron star Interior Composition Explorer, in the 1–10 keV energy band. First, we fitted the spectra with a thermal Comptonization model. Our analysis revealed the lack of a significant direct emission from a blackbody-like component, therefore we calculated the optical depth of the Comptonizing region, deriving an upper limit of 4.5, which suggests the presence of a moderately thick corona. We also attempted to fit the spectrum with an alternative model, i.e. a cold Comptonized emission from a disc and a direct thermal component from the neutron star, finding a similarly good fit. The source did not enter a full high luminosity/soft state throughout the outburst, with a photon index ranging from ∼1.7 to ∼2.2, and an average unabsorbed flux in the 1–10 keV band of ∼3.6 × 10−10 erg cm−2 s−1. We searched for the presence of Fe K-shell emission lines in the range ∼6.4–7 keV, significantly detecting a broad component only in a couple of observations. Finally, we conducted a time-resolved spectral analysis of the detected type-I X-ray burst, observed during the outburst, finding no evidence of a photospheric radius expansion. The type-I burst duration suggests a mix of H/He fuel.

Study of the ^{22}Mg Waiting Point Relevant for X-Ray Burst Nucleosynthesis via the ^{22}Mg(α,p)^{25}Al Reaction.

The ^{22}Mg(α,p)^{25}Al reaction rate has been identified as a major source of uncertainty for understanding the nucleosynthesis flow in Type-I x-ray bursts. We report a direct measurement of the energy- and angle-integrated cross sections of this reaction in a 3.3-6.9MeV center-of-mass energy range using the MUlti-Sampling Ionization Chamber (MUSIC). The new ^{22}Mg(α,p)^{25}Al reaction rate is a factor of ∼4 higher than the previous direct measurement of this reaction within temperatures relevant for x-ray bursts, resulting in the ^{22}Mg waiting point of x-ray burst nucleosynthesis flow to be significantly bypassed via the (α,p) reaction.

Impact of neutron star spin on Poynting–Robertson drag during a type I X-ray burst

ABSTRACT External irradiation of a neutron star (NS) accretion disc induces Poynting–Robertson (PR) drag, removing angular momentum and increasing the mass accretion rate. Recent simulations show PR drag significantly enhancing the mass accretion rate during Type I X-ray bursts, which could explain X-ray spectral features such as an increase in the persistent emission and a soft excess. However, prograde spin of the NS is expected to weaken PR drag, challenging its importance during bursts. Here, we study the effect of spin on PR drag during X-ray bursts. We run four simulations, with two assuming a non-spinning NS and two using a spin parameter of a* = 0.2, corresponding to a rotation frequency of 500 Hz. For each scenario, we simulate the disc evolution subject to an X-ray burst and compare it to the evolution found with no burst. PR drag drains the inner disc region during a burst, moving the inner disc radius outwards by ≈1.6 km in the a* = 0 and by ≈2.2 km in the a* = 0.2 simulation. The burst enhances the mass accretion rate across the innermost stable circular orbit ≈7.9 times when the NS is not spinning and ≈11.2 times when it is spinning. The explanation for this seemingly contradictory result is that the disc is closer to the NS when a* = 0.2, and the resulting stronger irradiating flux offsets the weakening effect of spin on the PR drag. Hence, PR drag remains a viable explanation for the increased persistent emission and soft excess observed during X-ray bursts in spinning NS systems.

Improved thermonuclear rate of42Ti(p,γ)43V and its astrophysical implication in therpprocess

Context.Accurate42Ti(p,γ)43V reaction rates are crucial for understanding the nucleosynthesis path of the rapid capture process (rpprocess) that occurs in X-ray bursts.Aims.We aim to improve the thermonuclear rates of42Ti(p,γ)43V based on more complete resonance information and a more accurate direct component, together with the recently released nuclear masses data. We also explore the impact of the newly obtained rates on therpprocess.Methods.We reevaluated the reaction rate of42Ti(p,γ)43V by the sum of the isolated resonance contribution instead of the Hauser-Feshbach statistical model. We used a Monte Carlo method to derive the associated uncertainties of new rates. The nucleosynthesis simulations were performed via the NuGrid post-processing code ppn.Results.The new rates differ from previous estimations due to the use of a series of updated resonance parameters and a direct S factor. Compared with the previous results from the Hauser-Feshbach statistical model, which assumes compound nucleus43V with a sufficiently high-level density in the energy region of astrophysical interest, large differences exist over the entire temperature region ofrp-process interest, up to two orders of magnitude. We consistently calculated the photodisintegration rate using our new nuclear masses via the detailed balance principle, and found the discrepancies among the different reverse rates are much larger than those for the forward rate, up to ten orders of magnitude at the temperature of 108K. Using a trajectory with a peak temperature of 1.95×109K, we performed therp-process nucleosynthesis simulations to investigate the impact of the new rates. Our calculations show that the adoption of the new forward and reverse rates result in abundance variations for Sc and Ca of 128% and 49%, respectively, compared to the variations for the statistical model rates. On the other hand, the overall abundance pattern is not significantly affected. The results of using new rates also confirm that therp-process path does not bypass the isotope43V.Conclusions.Our study found that the Hauser-Feshbach statistical model is inappropriate to the reaction rate evaluation for42Ti(p,γ)43V. The adoption of the new rates confirms that the reaction path of42Ti(p,γ)43V(p,γ)44Cr(β+)44V is a key branch of therpprocess in X-ray bursts.

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.

The Minimum Variation Timescales of X-Ray Bursts from SGR J1935+2154

The minimum variation timescale (MVT) of soft gamma-ray repeaters can be an important probe to estimate the emission region in pulsar-like models, as well as the Lorentz factor and radius of the possible relativistic jet in gamma-ray burst (GRB)–like models, thus revealing their progenitors and physical mechanisms. In this work, we systematically study the MVTs of hundreds of X-ray bursts (XRBs) from SGR J1935+2154 observed by Insight-HXMT, GECAM, and Fermi/Gamma-ray Burst Monitor from 2014 July to 2022 January through the Bayesian block algorithm. We find that the MVTs peak at ∼2 ms, corresponding to a light-travel time size of about 600 km, which supports the magnetospheric origin in pulsar-like models. The shock radius and the Lorentz factor of the jet are also constrained in GRB-like models. Interestingly, the MVT of the XRB associated with FRB 200428 is ∼70 ms, which is longer than that of most bursts and implies its special radiation mechanism. In addition, the median of MVTs is 7 ms, shorter than the median MVTs of 40 ms and 480 ms for short GRBs or long GRBs, respectively. However, the MVT is independent of duration, similar to GRBs. Finally, we investigate the energy dependence of MVT and suggest that there is a marginal evidence for a power-law relationship like GRBs, but the rate of variation is at least about an order of magnitude smaller. These features may provide an approach to identify bursts with a magnetar origin.

Sub-second infrared variability from the archetypal accreting neutron star 4U 1728−34

ABSTRACT We report on the first simultaneous high-time resolution X-ray and infrared (IR) observations of a neutron star low mass X-ray binary in its hard state. We performed $\approx 2\,$ h of simultaneous observations of 4U 1728−34 using HAWK-I@VLT, XMM–Newton, and NuSTAR. The source displayed significant X-ray and IR variability down to sub-second time-scales. By measuring the cross-correlation function between the IR and X-ray lightcurves, we discovered a significant correlation with an IR lead of $\approx 30 \!-\! 40\,$ ms with respect to the X-rays. We analysed the X-ray energy dependence of the lag, finding a marginal increase towards higher energies. Given the sign of the lag, we interpret this as possible evidence of Comptonization from external seed photons. We discuss the origin of the IR seed photons in terms of cyclo-synchrotron radiation from an extended hot flow. Finally, we also observed the IR counterpart of a type-I X-ray burst, with a delay of $\approx 7.2\,$ s. Although some additional effects may be at play, by assuming that this lag is due to light travel time between the central object and the companion star, we find that 4U 1728−34 must have an orbital period longer than $3\,$ h and an inclination higher than 8○.

EOS-dependent millihertz quasi-periodic oscillation in low-mass X-ray binary

ABSTRACT We studied the frequency and critical mass accretion rate of millihertz quasi-periodic oscillations (mHz QPOs) using a one-zone X-ray burst model. The surface gravity is specified by two kinds of equation of states: neutron star (NS) and strange star (SS). The base flux, Qb, is set in the range of 0–2 MeV nucleon−1. It is found that the frequency of mHz QPO is positively correlated to the surface gravity but negatively to the base heating. The helium mass fraction has a significant influence on the oscillation frequency and luminosity. The observed 7–9 mHz QPOs can be either explained by a heavy NS/light SS with a small base flux or a heavy SS with a large base flux. As base flux increases, the critical mass accretion rate for marginally stable burning is found to be lower. Meanwhile, the impact of metallicity on the properties of mHz QPOs was investigated using one-zone model. It shows that both the frequency and critical mass accretion rate decrease as metallicity increases. An accreted NS/SS with a higher base flux and metallicity, combined with a lower surface gravity and helium mass fraction, could be responsible for the observed critical mass accretion rate ($\dot{m}\simeq 0.3\dot{m}_{\rm Edd}$). The accreted fuel would be in stable burning if base flux is over than ∼2 MeV nucleon−1. This finding suggests that the accreting NSs/SSs in low-mass X-ray binaries showing no type I X-ray bursts possibly have a strong base heating.

Reanalysis of the X-Ray-burst-associated FRB 200428 with Insight-HXMT Observations

A double-peak X-ray burst from the Galactic magnetar SGR J1935+2154 was discovered to be associated with the two radio pulses of FRB 200428 separated by 28.97 ± 0.02 ms. Precise measurements of the timing and spectral properties of the X-ray bursts are helpful for understanding the physical origin of fast radio bursts. In this paper, we have reconstructed some information about the hard X-ray events, which were lost because the High Energy X-ray Telescope (HE) on board the Insight-HXMT mission was saturated by this extremely bright burst, and used the information to improve the temporal and spectral analyses of the X-ray burst. The arrival times of the two X-ray peaks by fitting the new Insight-HXMT/HE lightcurve with multi-Gaussian profiles are 2.77 ± 0.45 ms and 34.30 ± 0.56 ms after the first peak of FRB 200428, respectively, while these two parameters are 2.57 ± 0.52 ms and 32.5 ± 1.4 ms if the fitting profile is a fast-rise and exponential decay function. The spectrum of the two X-ray peaks could be described by a cutoff power law with cutoff energy ∼60 keV and photon index ∼1.4; the latter is softer than the result obtained in Li et al. when the two X-ray peaks appeared.

Sub-nanosecond time-resolved radiation measurement using x-ray focusing crystal spectrometers

In this paper, we describe a technique using a crystal spectrometer, a silicon-diode detector, and a filtered photoconductive detector to monitor photon energies in the L-shell (0.9–1 keV) and K-shell regimes for nickel and copper hybrid X-pinch x-ray sources. The detectors, system cabling, and an 8 GHz digital oscilloscope in combination enable time resolution better than 200 ps for photoconductive detectors and 700 ps for silicon-diode detectors of the K- and L-shell radiation signals, respectively. We substantially improve the relative timing of signals obtained using the oscilloscope by using an x-ray streak camera with a crystal spectrometer to monitor the L-shell line spectra and, separately, the K-shell line spectra relative to the continuum burst to better than 17 ps time resolution. This combination of instruments enabled and validated a new method by which plasma conditions in nickel and copper X-pinches can be assessed immediately before and after the ∼30 ps continuum x-ray burst produced by 370 kA hybrid X-pinches. In general, the method described here can be applied to observe otherwise highly filter-absorbed radiation in the presence of a broad spectrum of higher energy radiation by combining x-ray crystals and detectors.