Photospheric Radius Expansion and a Double-peaked Type-I X-Ray Burst from GRS 1741.9–2853

Accretion onto compact objects produces an array of phenomena, the study of which can illuminate the nature of these objects and the process of accretion itself. Type-I X-ray bursts represent a particularly dramatic consequence of accretion, wherein the flux of an accreting neutron star (NS) is observed to increase by an order of magnitude or more in the span of only a few seconds or less. In a typical type-I burst, this fast rise in flux is followed by an exponential decay that can take anywhere from tens to thousands of seconds to return to a persistent flux level. Type-I X-ray bursts can provide a wealth of information about the compact stellar remnants from which they originate. The observed fast-rise exponential-decay behavior is thought to be due to the runaway thermonuclear burning of accreted hydrogen and helium. This is possible only if the accretor has a surface onto which this material can accumulate in a thin shell that is unable to effectively cool during burning. Therefore, the detection of a type-I X-ray burst confirms the nature of the accreting source as an NS rather than a black hole. Additionally, in some cases, type-I bursts allow for the measurement of other important parameters such as the distance to the source or the rotation period of the NS. For a recent review of type-I X-ray bursts, refer to Galloway & Keek (2021).

During a type-I X-ray burst, the source spectrum can be described by the addition of a blackbody component on top of whichever model best describes the persistent emission. Analyzing how the spectrum changes throughout a type-I burst is necessary in some cases to distinguish between changes in flux that are energy independent and those that vary across different energy ranges, reflecting changes in the blackbody radius and temperature (Tawara et al. 1984). The latter case is often evidence for photospheric radius expansion (PRE). During PRE, the blackbody that models the burst emission is shown to rapidly increase in radius while decreasing in temperature as the luminosity remains constant at the peak of the burst. This behavior has been attributed to radiation pressure lifting material from the NS surface upon reaching the Eddington limit. Due to the relatively narrow range of acceptable values for NS mass and radius, bursts that demonstrate PRE are often considered "standard candles," which can be used to estimate the distance to a bursting source (van Paradijs 1978; Lewin et al. 1984). Some sources show significant scatter in peak luminosity across many PRE events, however (Kuulkers et al. 2003). As such, these events may serve as approximate standard candles, and source properties such as disk inclination and the composition of burning material must be taken into account when calculating the Eddington luminosity.

Additionally, timing analysis before, during, and after type-I bursts can reveal oscillatory behavior related to burning on the NS surface. Coherent high-frequency oscillations have been observed during some type-I bursts, most often during the decay in brightness (Watts 2012). These oscillations have been attributed to brightness asymmetries induced by burning on the NS surface, implying that the frequency of burst oscillations corresponds to the rotational frequency of an NS. In some systems, quasi-periodic oscillations (QPOs) at millihertz frequencies have also been observed in the persistent emission leading up to type-I bursts. This may indicate that as material accumulates, oscillatory nuclear burning occurs on the surface prior to the onset of the runaway burning that causes type-I bursts.

GRS 1741.9–2853 is a low-mass X-ray binary that was first reported to exhibit type-I X-ray bursts in 1999 (Cocchi et al. 1999). The source resides near the Galactic Center, about 10' from Sgr A* at α = 17^h45^m02^s, (Muno et al. 2003). As it accretes matter from its main-sequence companion, it transitions between a low-flux hard spectral state during which angular momentum transfer in the cold disk is inefficient, leading to low mass accretion rates, and a high-flux soft spectral state when the accretion disk enters a hot, efficiently accreting state (Lasota 2001). GRS 1741.9–2853 is a member of the very faint class of transients, meaning that it reaches a peak of 2–10 keV luminosity in the range of 10³⁴−10³⁶ erg s⁻¹ (Wijnands et al. 2005). The first detection of GRS 1741.9–2853 was made by the GRANAT satellite and reported in 1990, but was attributed to the nearby source 1E 1741.7–2850 (Mandrou 1990). Later analysis resolved the source (Sunyaev 1990; Syunyaev et al. 1991). GRS 1741.9–2853 undergoes periods of outburst every ∼2 yr, with typical outbursts ranging in duration from 5 to 15 weeks, resulting in a duty cycle of ∼10%.

An investigation of the bursting behavior of GRS 1741.9–2853 was presented by Trap et al. (2009). The authors analyzed 15 type-I bursts observed during two periods of outburst in 2005 and 2007 by INTEGRAL JEM-X, Swift BAT and XRT, and XMM-Newton. The bursts had a typical recurrence time between 79 ks and 2 Ms. From this analysis they were able to infer a distance to the source of 7 kpc, and they determined that the bursting behavior was consistent with pure He burning. During another period of outburst in 2013, the source was observed serendipitously during two observations by the Nuclear Spectroscopic Telescope Array (NuSTAR) (Barrière et al. 2015). One type-I X-ray burst lasting 800 s was observed, which was consistent with mixed H/He burning. It exhibited mild PRE, allowing the authors to infer a distance to the source that was consistent with the previous estimate of 7 kpc.

Again in 2020 April the source was observed by the XRT instrument aboard the Neil Gehrels Swift Observatory to be increasing in brightness, rising from ∼10⁻² counts s⁻¹ to ∼9 × 10⁻² counts s⁻¹ over the course of two days, indicating that the source was entering a period of outburst (Degenaar et al. 2020). This detection was followed up by NuSTAR on 2020 May 7.

We report on the detection of two type-I X-ray bursts during the 2020 May NuSTAR observation of GRS 1741.9–2853 in outburst, as well as spectral and timing analyses of the persistent emission. We begin by describing the observation and the methods of data reduction and analysis in Section 2. Next, we present an analysis of the NuSTAR light curve, including modeling of each of the two type-I bursts in Section 3, the first of which is shown to peak twice. In Section 4 we present the persistent and burst spectra and demonstrate evidence for PRE during the second burst. For completeness, in Section 5 we present an analysis of the timing features of both the persistent and burst emission, including a search for QPOs and burst oscillations. We end with a discussion of our results in Section 6.

NuSTAR, launched in 2012 June, is the first high-energy focusing X-ray telescope (Harrison et al. 2013). It is composed of two focal planes, each paired with a set of focusing optics with a focal length of 10 m. The focal planes, FPMA and FPMB, are each composed of four cadmium zinc telluride detectors attached to custom readout electronics.

GRS 1741.9–2853 was observed by NuSTAR on 2020 May 7 for a total exposure time of 28 ks (OBSID 90601317002). We extracted the data products using NuSTARDAS version 2.0.0 and CALDB version 20200826. Photon arrival times were shifted to the barycenter of the solar system to eliminate residuals due to the orbital motion of NuSTAR. Instrumental effects on the photon arrival times due to changes in the temperature of the onboard clock were corrected using clockfile v110 generated on 2020 September 12. For the production of source spectra as well as cospectra, we used circular extraction regions with radius 100'', centered on the source using automatic centroid detection in DS9. To produce background spectra, we used extraction regions of radius 60'', residing on the same detector as the source region but sufficiently removed from it so as to exclude source counts. Additionally, we found that in this case, the automated pipeline overcorrected for changes in the multilayer insulation (Madsen et al. 2020), resulting in inconsistencies between FPMA and FPMB below 5 keV. Therefore, we chose to revert to the previous FPMA ancillary response file. ¹⁰

Timing analyses, such as the production of cospectra and calculation of epoch-folding statistics, were performed using the Python package Stingray (Huppenkothen et al. 2016). Spectral modeling was performed using the X-ray spectral analysis package Xspec (v12.11.1 Arnaud 1996). Spectra were binned using the variable binning algorithm described by Kaastra & Bleeker (2016), which takes into account both the number of photons in a given energy bin as well as the average energy of the photons in that bin. We note that this binning procedure does not allow one to specify an exact energy range, meaning that the binned spectra may not extend all the way down to 3 keV. We restricted our spectral analysis to an energy range of 3–40 keV for persistent emission and a range of 3–20 keV for burst emission. At higher photon energies, background counts begin to contribute significantly to the overall spectra. All spectral fitting was performed using the Cash statistic (Cash 1979), but throughout the paper we present chi-squared fit statistics in order to provide an idea of the quality of various fits using a formalism that is easy to interpret and gives a clear comparison between models. For the purpose of readability, spectra shown in figures have been further rebinned such that each bin has at least 5σ significance. All light curves shown in this paper are the sum of the simultaneous light curves observed by FPMA and FPMB, and have been background subtracted and corrected for variations in the livetimes of the detectors. Errors quoted throughout this paper represent 90% confidence regions unless otherwise stated.

The total (3–78 keV) NuSTAR light curve is shown in Figure 1, where the top panel shows the light curve binned in 0.5 s increments in order to demonstrate the bursting behavior, and the bottom panel shows the light curve binned in 10 s increments and restricted to count rates below 85 count s⁻¹ in order to show the persistent flux level. Throughout the observation, GRS 1741.9–2853 showed little variability aside from two dramatic increases in count rate occurring around 12,700 and 52,600 s. Additionally, the persistent count rate appeared to increase gradually by about 30% from about 50 counts s⁻¹ to about 65 counts s⁻¹.

Zoom In Zoom Out Reset image size

Due to their short duration and dramatic increase in count rate, the features at roughly 12,700 and 52,600 s resemble type-I bursts. Indeed upon closer inspection, they exhibit the fast rise followed by an exponential decay, which are characteristic of type-I bursts, confirming this classification. Figure 2 shows the burst light curves for the entire 3–78 keV NuSTAR band binned into 0.5 s intervals. The bursts not only differ in maximum count rate but in their structure as well. The first burst, hereafter Burst 1, appears to have a longer rise time and exhibits two clear peaks in count rate, while the second burst, hereafter Burst 2, exhibits a much faster rise and a single peak.

Zoom In Zoom Out Reset image size

In order to better understand the structure of the two type-I bursts, we began by fitting the light curves to a simple fast-rise exponential -decay (FRED) model, given by

for t > t₀, where t₀ is the time at burst onset, and τ_R and τ_D are the rise and decay times, respectively, A is a factor that determines the height of the burst, and C is the persistent count rate. From Equation (1), it can be seen that the burst peak occurs at . The height of the burst above the persistent contribution is therefore given by . The fits to the simple FRED model, including residuals, are shown in the top panels of Figure 2, in which the light curves are plotted in logarithmic space in order to better illustrate the structure of each burst.

As expected, Burst 1 shows clear residuals around the peak of the FRED model indicating a double-peaked structure. We therefore proceeded to model Burst 1 as the sum of two "subbursts" in quick succession, each modeled as a FRED burst. The fit to this model is shown in the bottom-left panel of Figure 2. This results in a significantly better fit, eliminating the residuals around the two peaks. Fitting reveals that the first subburst has a significantly slower rise time than the second subburst but a somewhat faster decay time. It can also be seen in Figure 2 that the two subbursts have comparable peak count rates. In Section 6 we discuss the physical mechanisms by which this double-peaked structure could be realized.

The second burst bears more resemblance to a typical type-I burst. The FRED model describes the shape of the burst well for early times, but excess emission can be seen in the tail of the burst. We therefore adopt a phenomenological model similar to the one used by Barrière et al. (2015), in which the rise, peak, and beginning of the decay are modeled by a FRED curve, and the late-time emission is modeled as a power-law decay rather than an exponential decay. This model can be written as

where t_tail is the time at which the emission transitions from an exponential tail to a power-law tail, γ is the power-law index, and such that the exponential tail and the power- law tail have the same count rate at t = t_tail. This model successfully describes the burst emission both at early and late times, as demonstrated by the bottom-right panel of Figure 2.

Table 1 lists the values of the fitted and calculated parameters for Burst 1 and Burst 2. Both subburst 2 and Burst 2 show short rise times and similar decay times. The increase in count rate, which is apparent from the persistent light curve, is reflected by the difference in the fitted values of C between Burst 1 and Burst 2. The magnitude of the first subburst of Burst 1 is not precisely constrained due to degeneracy with other values and due to the limited number of constraining data points as compared to subburst 2, which dominates the light curve beginning about 8 s after the onset of Burst 1. With the onsets of each burst well constrained, we are able to calculate the recurrence time between Burst 1 and Burst 2, Δt = 39,862.2 ± 2.3 s. This recurrence time is less than half the shortest recurrence time previously reported for GRS 1741.9–2853 (Trap et al. 2009). Due to the orbital gaps in the light curve, however, it is possible that the two bursts were not consecutive.

In order to better understand the bursting behavior of GRS 1741.9–2853, it is necessary to characterize both the persistent and burst spectra, as well as to understand how the spectrum changes during each burst. We began by producing good time intervals (GTI) for the persistent emission as well as for each of the two type-I bursts. For the first burst, we defined the GTI to begin at the onset of the first subburst and to end after five decay times (τ_D ) had elapsed following the onset of the second subburst. This results in an interval with length 41 s. For the second burst, we similarly defined the GTI to begin at the onset of the burst and to end after five decay times (τ_D ) had elapsed, resulting in an interval with length 24 s. The GTI for persistent emission was defined such that it excluded the bursts, with a 100 s buffer prior to the onset of each burst, and a 500 s buffer following the end of each burst's GTI, such that in total 1265 s of the observation was excluded. Using these good time intervals, we were able to extract the persistent and burst spectra. When extracting burst spectra, we loosened the event filters by specifying the status expression "STATUS =b0000xx000xxxx000" when running nupipeline. We therefore were able to avoid removing source photons that would otherwise be mistaken for spurious events at such high count rates.

4.1. Persistent Emission

The persistent spectrum is described well by the standard physical picture of an NS surface or boundary layer emitting blackbody radiation, which is Compton upscattered into a power-law-like component by a hot corona. In Xspec, we represented the blackbody emission using bbodyrad and the upscattered emission using nthcomp (Zdziarski et al. 1996; Życki et al. 1999). The seed photon temperature for nthcomp was tied to the temperature of the blackbody. An absorbing column, represented by tbabs, was also applied to the sum of the model components, with molecular abundances described by Wilms et al. (2000) and cross sections described by Verner et al. (1996). Fitting with this model resulted in a reduced chi-squared statistic of (χ² = 408.63; d. o. f. = 249). We observed a clear Fe Kα emission feature in the resulting residuals, so we also included a Gaussian component at ≈6.5 keV. The addition of this component improved the fit significantly, bringing the reduced chi-squared down to (χ² = 293.5; d. o. f. = 246), confirming that the component is necessary to accurately describe the observed spectrum.

Additionally, we found that the Fe line emission could be accounted for by replacing the Gaussian component with a relativistic disk reflection model. We found that the spectrum could be described well by the reflection off an accretion disk of a cutoff power law approximating the nthcomp component using the relxill ¹¹ model (Dauser et al. 2014; García et al. 2014), or by the reflection of the blackbody component using a modified version of relxill, relxillNS (García et al. 2021). Although reflection of the blackbody using relxillNS can account for the line emission, it cannot account for the hard emission on its own, and the nthcomp component is still necessary to accurately model the spectrum up to high photon energies.

For both of our reflection models, the disk emissivity index was assumed to be constant throughout the disk, the disk density was fixed at n = 10¹⁵ cm⁻³, the spin parameter was fixed at a = 0, and the reflection fraction was fixed to a value of −1, such that only the reflected component was modeled by relxill and relxillNS (in other words, we added the reflected component to the nonreflection model described above, rather than replacing the direct blackbody and Comptonized emission with the relxill models). In the case of blackbody reflection, the emissivity index of the disk was fixed to q = 3.5, corresponding to irradiation of the disk by the NS surface or boundary layer (Wilkins 2018), and the temperature of the incident blackbody was tied to that of the bbodyrad component. For the case of power-law reflection, we fixed the emissivity index at a value of q = 3 corresponding to a disk corona geometry. We tied the relxill power-law index to that of nthcomp, and we related the relxill cutoff energy to the nthcomp electron temperature as E_cut = 3kT_e − 2 keV such that the shape of the incident power law roughly approximated the component described by nthcomp.

The best-fit parameters for all three models are listed in Table 2, where the nonreflection model is referred to as Model 1, the reflected blackbody (relxillNS) is referred to as Model 2, and the reflected cutoff power law (relxill) is referred to as Model 3. We find that the three models provide a similarly good fit: the reduced chi-squared statistics for Models 1, 2, and 3 are , 1.17 (d. o. f = 244), and 1.13 (d. o. f = 244), respectively. The observed spectra as well as the residuals and components of the Model 3 are shown in the left panel of Figure 3. The three models are largely consistent with one another, showing little change between their shared components. Additionally, the spectral parameters are consistent with previous measurements (Barrière et al. 2015) with somewhat lower column density and blackbody temperature.

Zoom In Zoom Out Reset image size

Table 2. Spectral Parameters Determined by Fitting the Persistent Spectra in Xspec

Component	Parameter	Model 1	Model 2	Model 3
tbabs	NH (1022 cm−2)
bbodyrad	kT (keV)	1.26 ± 0.01		1.24 ± 0.01
		29 ± 2	20 ± 4
nthcomp	Γ
	kTe (keV)
	Norm (10−3)
Gaussian	E (keV)	6.48 ± 0.07	⋯	⋯
	σ (keV)		⋯	⋯
	K (10−4 photon cm−2 s−1)	4.5 ± 1.0	⋯	⋯
	Equivalent width (eV)		⋯	⋯
relxill/relxillNS	q	⋯	3.5 a	3.0 a
	Inclination ({°})	⋯
	Rin (Rg)	⋯
		⋯
	AFe (Solar)	⋯	3.0 ± 1.9
	Norm (10−4)	⋯
χ2/d. o. f.		293.5/246	284.5/244	274.9/244

Notes. Model 1 corresponds to a Comptonized blackbody without disk reflection, Model 2 corresponds to the relativistic disk reflection of an irradiating blackbody, and Model 3 corresponds to relativistic disk reflection of an irradiating cutoff power law.

^a The values of the disk emissivity index were fixed during fitting.

Download table as:  ASCIITypeset image

From the nonreflection model, we find an apparent blackbody radius of R = 3.8 ± 0.1 km assuming a distance to the source of 7 kpc. Given that the gravitational redshift on the surface of the NS, given by , is significant, the radius measured by fitting to a blackbody is larger than the radius in the emitting frame by a factor of (1 + z)^3/2. Assuming M = 1.4 M_⊙ and R_NS = 10 km, we get 1 + z = 1.31. This yields an actual blackbody radius of 2.5 km, consistent with accretion onto a small band rather than spherical accretion onto the entire surface of the NS. The relativistic reflection models provide similar values for the ionization of the disk, ξ ∼ 10³ erg cm s⁻¹, as well as the iron abundance, A_Fe, which is loosely constrained at about three times the solar value. The iron abundance, however, depends strongly on the disk density, which cannot be varied using relxill and which we could not constrain with relxillNS. We therefore advise against interpreting the iron abundance as the true physical value. Both relxill and relxillNS yield a low inclination of ∼15° and an inner disk radius consistent with ∼50 R_g , corresponding to about 100 km for an NS mass of M = 1.4 M_⊙.

Having modeled the persistent emission, we were able to constrain the persistent flux and corresponding mass accretion rate. From Galloway et al. (2008), the bolometric luminosity, L_bol, is related to the accretion rate per unit area, , by . Solving for , as in Barrière et al. (2015), gives

where F_bol is the bolometric flux, M and R_NS are the NS mass and radius, respectively, d is the distance to the source, and z is the gravitational redshift at the surface of the NS. The values that we measured for the bolometric flux and the resulting calculated values of bolometric luminosity and accretion rate, assuming a distance of d = 7 kpc, are shown in Table 3.

Table 3. Physical Quantities Calculated for the Persistent Emission, Based on the Nonreflection Spectral Model

Quantity	Value	Units
Fbol a	1.20 ± 0.01	10−9 erg cm−2 s−1
Lbol b		1036 erg s−1
c	3.9	103 g cm−2 s−1
d	4.5	percent
e	7.8	10−10 M⊙ yr−1

Notes.

^a Bolometric (0.1–100 keV) unabsorbed flux. ^b Bolometric unabsorbed luminosity assuming a distance of 7 kpc. ^c Mass accretion rate per unit surface area as determined by Equation (3). ^d Ratio of the mass accretion rate to the Eddington limited rate, , assuming M = 1.4 M_⊙ and R_NS = 10 km. ^e Mass accretion rate assuming an NS radius of 10 km.

Download table as:  ASCIITypeset image

4.2. Burst Emission

The burst spectra were modeled by freezing the persistent model (for simplicity, we used the nonreflection model) and adding an additional blackbody component that was allowed to vary in both temperature and normalization. ¹² Rather than pinning the spectral parameters of the persistent model to the values determined for the full observation, we determined the persistent model parameters relevant to each burst using spectra extracted near the time of that burst, excluding the interval of variability between ∼15 and ∼35 ks.

For each burst, we produced GTIs of varying lengths, beginning at the burst onset and ending after several decay times. The GTI lengths were set such that each bin contained a roughly equal number of counts and so that we could achieve good temporal resolution during the first several seconds of the burst while also maintaining enough photons per time bin to enable robust spectral fitting. For each of these intervals, the spectrum was again extracted and fit to the persistent-plus-blackbody model described above, allowing the blackbody temperature and normalization to vary. The resulting parameters, including the temperature, the unabsorbed bolometric luminosity, and the radius of the blackbody (assuming a distance of 7 kpc), are shown in Figure 4. Each data point corresponds to a GTI for which spectra were analyzed.

Zoom In Zoom Out Reset image size

Time-resolved spectroscopy reveals that the double-peaked structure observed in the count rate during Burst 1 is also visible in the luminosity and temperature evolution of the burst. The first subburst reached a blackbody temperature of kT_bb = 2.2 keV and a luminosity of L_bol = 4.4 × 10³⁷ erg s⁻¹, while the second subburst reached a blackbody temperature of kT_bb = 2.3 keV and a luminosity of L_bol = 6.8 × 10³⁷ erg s⁻¹. This analysis confirms that the double-peaked structure is intrinsic to the source. That is, it is not the result of an increase in the fraction of incident photons energies falling below NuSTAR's lower energy limit of 3 keV due to a softening of the spectrum at a constant luminosity. The available data do not suggest that the NS underwent PRE during either of the subbursts that make up Burst 1, but our analysis is limited by our inability to resolve these short subbursts into smaller time bins due to insufficient counts.

We showed in Section 3 that Burst 2 exhibited a very different structure compared Burst 1. This is borne out in our time-resolved spectroscopic analysis as well. The blackbody luminosity shows a single peak, reaching a bolometric flux of (L_peak = 1.7 × 10³⁸ erg s⁻¹, assuming a distance of 7 kpc) within a few seconds of the onset of the burst and a long decay time. On the other hand, the apparent blackbody temperature and radius show variability during the luminosity peak, with the temperature decreasing rapidly from 2.7 to 1.9 keV in the first second of the burst, followed by an increase up to 2.5 keV only a few seconds later. This temperature evolution is accompanied by a rapid expansion and contraction of the apparent blackbody radius, which reaches a maximum of 8.7 km before contracting to about 5 km. We note that the observed values of the blackbody radius and temperature are under- and overestimates, respectively, of their actual physical values. Scattering of the blackbody emission in the NS atmosphere results in a harder emergent spectrum. The overall effect of this scattering is to introduce a color correction factor, f_c, to the effective temperature so that , where T_obs is the temperature determined via modeling of the observed spectrum. Similarly, the effective radius is related to the observed radius by . The color correction factor may vary between 1.4 and 1.7 for luminosities approaching the Eddington limit (Suleimanov et al. 2011). Therefore, at the peak of Burst 2, the blackbody temperature may be a factor of ∼0.6 lower than that observed and the blackbody radius a factor of ∼3 larger.

The behavior described above is the signature of a "touchdown" event wherein the photosphere undergoes expansion due to radiation pressure, then rapidly falls back down to the surface such that luminosity remains constant while the contraction causes the temperature to increase. We therefore conclude that the NS underwent PRE during Burst 2. Because this phenomenon is generally attributed to the balance of radiation pressure and gravitational force at the Eddington limit (L_Edd ∼ 2 × 10³⁸ erg s⁻¹ for an NS with typical mass M =1.4 M_⊙), it is often assumed that PRE allows for the measurement of the source distance. We defer estimation of the distance and further discussion until Section 6.

We searched both the persistent and burst emission for timing features including QPOs and coherent oscillations. We began by producing cross-spectra between FPMA and FPMB for 113 × 256 s intervals included in the persistent GTI described in Section 4. We analyzed the cospectrum, the real part of the cross-spectrum, due to its advantages described by Bachetti et al. (2015), namely the fact that it eliminates contributions to the power spectrum introduced by Poisson noise and dead time, the latter of which can be significant even for moderate NuSTAR count rates. In order to produce the cospectra, we used Stingray to produce light curves using only events within the 100'' source extraction region and with energy between 3 and 78 keV. We binned the light curves with a resolution of 2048⁻¹ s and ensured that the light curves were simultaneous between FPMA and FPMB.

We inspected the resulting dynamical cospectrum as well as the averaged cospectrum to look for coherent signals, particularly around the previously detected signal at 589 Hz (Strohmayer et al. 1997). There is no indication of such a feature, and the cospectrum of the persistent emission does not appear to deviate from a simple power law. We confirmed this visual inspection by performing a chi-squared fit to a power-law model, which yielded a reduced chi-squared of χ²/d. o. f. = 0.98. Nonetheless, we performed an automated search for QPO features by fitting the average cospectrum to a power law with an added Lorentzian feature, such that the total model is given by

where f(ν) is the rms-normalized spectral power density at frequency ν, A is the normalization of the power-law component, α is the power-law index, B is the integrated power underneath the Lorentzian component, ν₀ is the frequency at the peak of the Lorentzian component, and γ gives the width of the Lorentzian and is related to the quality factor by . We scanned ν₀ over the frequency range, bounded on the lower end by the length of the intervals, and bounded on the upper end by the Nyquist frequency, equal to twice the binning frequency. For each value of ν₀, we fit the average cospectrum to the model described by Equation (4), allowing all parameters aside from ν₀ to vary and compared the resulting χ² fit statistic to the statistic attained by fitting the cospectrum to a power law alone. This produced a Δχ² distribution over the range of allowed frequencies. The greatest fit improvement is ∣Δχ²∣ < 12, compared to a baseline statistic of , indicating that the addition of a Lorentzian signal to the underlying power-law spectrum is not warranted. We nonetheless attempted to fit QPO signals at several peaks in the ∣Δχ²∣ distribution, this time allowing ν₀ to vary. We found that none of these signals had significance greater than 1.7σ, where significance is defined as the ratio of the integrated power under the Lorentzian to the error of the power, B/σ_B .

We performed a similar analysis for the burst emission as well as a more focused search for QPOs in two 10 ks intervals preceding each burst. We produced cospectra for 115 × 1 s intervals during Burst 1 and 98 × 1 s intervals during Burst 2 and averaged all 213 cospectra together. We note that for the purposes of timing analyses, we extended each of the burst GTIs by 75 s in order to include the burst decays. For the preburst emission, we produced a total of 34 cospectra, each corresponding to an interval of 256 s, 18 of which precede Burst 1 and 16 of which precede Burst 2. We averaged over each of the preburst epochs separately, and, prompted by the work of Revnivtsev et al. (2001), we filtered the preburst emission to include only photon energies in the range 3–12 keV. Following the QPO search method we described for the persistent emission yields similar results: we do not significantly detect QPOs at any frequency during the bursts nor during the preburst intervals.

We also investigated the burst cospectra for evidence of coherent oscillations corresponding to a spin period. Visual inspection of both the dynamical and averaged (across both bursts) cospectra does not indicate the presence of coherent oscillations. We find a 99% confidence upper limit of 16% for the fractional rms integrated between 588 and 590 Hz, exclusive.

We further searched for burst oscillations by performing a dynamical epoch-folding search on the burst emission. We analyzed each burst in 1 s intervals, stepping through the bursts in 0.5 s increments, resulting in overlapping time bins. For each focal plane module, we produced pulse profiles corresponding to a range of oscillation frequencies near 589 Hz. We folded only events between 3 and 12 keV, again restricting our analysis to events within the source extraction region used described above. The resulting pulse profiles were further corrected for variations in dead time using the method described in Madsen et al. (2015). For each corrected pulse profile, we calculated the statistic (Buccheri et al. 1983) as well as the corresponding probability that each of these statistics was not produced by noise. For each burst, we were thus left with two distributions, one for FPMA and another for FPMB, of the probability of detection of oscillations over a range of frequencies. We did not find any significant signal at any time during the bursts, and those spurious signals that we did observe only appeared in a single focal plane module.

6.1. Explaining the Double-peaked Burst

6.2. Photospheric Radius Expansion

Our analysis of the spectral evolution of Burst 2 revealed that the burning material underwent a process of rapid expansion and contraction. We have interpreted this behavior as evidence for PRE, during which the burning material reached the Eddington luminosity and was therefore momentarily lifted from the NS surface by radiation pressure. The evolution of the blackbody radius, temperature, and luminosity strongly resembles the PRE burst observed by Barrière et al. (2015), with the major difference being the overall shorter timescales observed in Burst 2 compared to the hydrogen-rich burst observed by those authors.

Compared to the burst observed by Barrière et al., the short rise and decay time observed for Burst 2, in addition to the relatively high persistent flux leading up to the burst (see Table 4), may indicate that at the time of observation the source was in the stable H-burning regime. In other words, while the accreted material is likely to contain a mixture of H and He, the accreted H burns stably on the NS surface rather than leading to runaway thermonuclear burning. This in turn leads to pure He bursts (Trap et al. 2009). In addition to the burst structure and persistent flux, the recurrence timescale and burst fluence provide further evidence for pure He burning during Burst 2. The recurrence timescale, τ_rec, and the mass accretion rate, , can be combined to calculate the ignition column: . Next, the energy released per unit mass, _nuc, can be calculated using the relation , where E_b is the total energy radiated during the burst and is related to the burst fluence, f_b, by E_b = 4π d² f_b. We may combine the relations listed above with Equation (3) to obtain an equation for _nuc that does not depend on the distance to the source but rather on the ratio of the observed burst fluence to the integrated persistent flux:

for an NS with M = 1.4M_⊙ and R_NS = 10 km. Due to the orbital gaps in the observation, we were only able to determine an upper limit on the recurrence timescale, τ_rec ≲ 4 × 10⁴ s. Integrating the blackbody flux during Burst 2 gives a burst fluence of f_b = (2.8 ± 0.1) × 10⁻⁷ erg cm⁻². Thus, we arrive at an energy per unit mass of _nuc ≳ 1.4 × 10¹⁸ erg g⁻¹ or an energy per nucleon of Q_nuc ≳ 1.4 MeV nucleon⁻¹. Given Q_nuc =(1.35 + 6.05X)MeV nucleon⁻¹ (Goodwin et al. 2019), this corresponds to a H fraction, X, consistent with zero.

Table 4. Comparison of PRE Bursts Presented in This Work and in Barrière et al. (2015)

Reference	Fpers a	Fpeak b	Xc	LEdd d	de
	(10−11 erg cm−2 s−1)	(10−8 erg cm−2 s−1)	⋯	(1038 erg s−1)	(kpc)
Barrière et al. (2015)			0	3	8.3 ± 0.7
			0.7	1.7	6.3 ± 0.5
This work	130 ± 2		0	3	9.0 ± 0.5
			0.7	1.7	7.0 ± 0.4

Notes.

^a Unabsorbed bolometric flux calculated for the period leading up to the burst. This period corresponds to epoch O4 in Barrière et al. (2015) and to the final four orbits of the observation presented in this work. ^b Unabsorbed bolometric flux calculated for time slice corresponding to the peak of the burst. ^c Hydrogen mass fraction assumed when calculating the Eddington luminosity. ^d Eddington luminosity calculated using Equation (6). ^e Distance calculated assuming that the peak flux corresponds to the Eddington luminosity.

Download table as:  ASCIITypeset image

With this in mind, we may calculate the distance to the source. The Eddington luminosity as measured by a distant observer is given by Lewin et al. (1993):

where κ₀ is the opacity for pure He and is given by κ₀ = 0.2 cm² g⁻¹, M is the mass of the accretor, R is the radius at which the Eddington luminosity is being calculated, T is the temperature of the material at touchdown, and X is the hydrogen mass fraction of the material, which we assume to be zero. We assume a typical NS mass of M = 1.4M_⊙, and because the peak flux is achieved at the time of touchdown when the burning material is assumed to have returned to the NS surface, we similarly assume a typical NS radius of R = 10 km. The temperature, T, that appears in Equation (6) is that of the material in the emission reference frame. In reality, the temperature measured by a distant observer via fitting to a blackbody differs from the temperature in the frame of the emitting material, . Given a temperature at touchdown of kT_∞ = 2.44 keV, a redshift of (1 + z) = 1.31 (see Section 4.1), and a color correction factor of f_c = 1.7, we arrive at an Eddington luminosity of L_Edd,∞ = 3 × 10³⁸ erg s⁻¹. Combined with our measured peak unabsorbed, bolometric flux of , this corresponds to a distance to the source of d = 9.0 ± 0.5 kpc.

This result is consistent with previous estimates that have placed the source at distances ranging from 5 to 9 kpc depending upon the assumed material composition and the method of calculating the Eddington luminosity (Cocchi et al. 1999; Galloway et al. 2008; Trap et al. 2009; Barrière et al. 2015). Our measurement lies on the most distant end of these estimates. This could be the result of a number of effects. As we have shown, the material composition has a significant effect on the Eddington luminosity and therefore on any measurement of the distance. For example, if we assume solar composition, X = 0.7, rather than pure He burning, the resulting distance is d = 7.0 ± 0.4 kpc. In order to remain consistent with our measurements of the burst fluence, this would require a recurrence time that is a factor of ∼4 smaller than the time that elapsed between Burst 1 and Burst 2. This in turn would imply that an additional burst took place during one of the orbital gaps of the observation. It is interesting to note that the peak flux we have observed is nearly consistent with the peak flux presented by Barriere et al., while the structures of the corresponding bursts are very different. This may represent a challenge to the method of inferring the composition of burning material via light-curve analysis alone. Throughout this paper we have assumed isotropic emission when calculating burst and persistent fluxes. However, the inclination of the system can have a significant effect on the actual values of these quantities (Galloway et al. 2003). Spectral modeling of the persistent emission suggests that the inclination of the inner disk is less than 30°. Depending upon the disk geometry, this may result in an overestimate of the burst flux by a factor of 1.5–2, and of the persistent flux by up to a factor of 2.5 (He & Keek 2016). Therefore, the distance may be underestimated, both in this work and in previous work, which also assumed isotropic emission, by a factor of up to .The bursting behavior presented by Trap et al. (2009) indicates a possible correlation between peak flux and burst fluence. We therefore speculate that those bursts could actually have achieved the same peak luminosity (modulo composition effects) that only appeared to vary due to changes in inclination. However, because Trap et al. are unable to claim that the bursts undergo PRE, it may not be appropriate to assume that the bursts reached the Eddington luminosity, meaning that this interpretation may not be accurate. Future investigations of the long-term variability of GRS 1741.9–2853 could reveal whether changing disk inclination has a significant impact on the measured peak flux of PRE bursts.

We have presented light-curve, spectral, and timing analyses of two type-I X-ray bursts as well as persistent emission originating from the source GRS 1741.9–2853 during a period of outburst observed by NuSTAR. This represents the first time the source has been observed by NuSTAR at such a high persistent flux—the source was over 15 times brighter than during the previous NuSTAR observation that took place in 2013—which allowed us to analyze the bursting behavior in an accretion regime previously unexplored by NuSTAR for this source. Additionally, because we observed two bursts, we were able to obtain an upper limit on the recurrence timescale between the bursts, τ_rec ≲ 4 × 10⁴ s. This upper limit is half of the shortest recurrence time previously reported between bursts originating from GRS 1741.9–2853 (Trap et al. 2009), and the measurement of the recurrence time served as an additional piece of evidence that the burning material was H deficient. Aside from the two bursts, the source appeared to show only slight variability. We confirmed this by analyzing the cospectra taken between light curves observed by the two focal plane modules. In all cases we found the cospectra to be relatively featureless, being described well by a single power-law model. Additionally, we performed an epoch-folding analysis that did not produce evidence for coherent burst oscillations.

We found that the spectra of the persistent emission were modeled well by an absorbed blackbody plus upscattered emission originating from a hot corona or accretion flow, in addition to Fe Kα emission, which can be modeled either phenomenologically using a simple Gaussian emission line or self-consistently using relativistic disk reflection models. Having modeled the persistent emission, we were able to perform a time-resolved spectroscopic analysis on each of the type-I bursts. We found that the first burst exhibited a double-peaked structure, and we showed that the structure we observed did not appear to result from the rebrightening of an ongoing burst but rather could be described as two "subbursts" occurring in quick succession. This is the first confirmed case of such a burst originating from GRS 1741.9–2853, placing it within a population of only a handful of multiple-peaked bursters.

The first burst was significantly less bright and energetic than the second burst, which exhibited a canonical FRED structure in addition to an extended tail. Spectroscopic analysis showed that the blackbody describing the emission during the second burst underwent a period of rapid contraction accompanied by a rise in temperature, all occurring at a constant luminosity. We therefore inferred that the source had undergone photospheric radius expansion, or PRE, which typically implies that the source has reached the Eddington luminosity. Based on this assumption, we calculated a distance to the source of d = 9.0 ± 0.5 kpc, which is somewhat larger than, but still consistent with, previous estimates. However, the value of the Eddington luminosity depends strongly on the composition of the burning material. While the observed properties of the burst as well as the persistent flux and recurrence timescale are consistent with pure He burning, it is not possible to determine the composition with absolute certainty. Future observations of type-I X-ray bursts originating from GRS 1741.9–2853 may provide additional opportunities to determine how the composition of burning material varies with persistent flux and how this affects bursting behavior, as well as whether changing disk inclination is responsible for the variation in peak flux between PRE bursts at different times.

J.J. acknowledges support by the Tsinghua Astrophysics Outstanding (TAO) Fellowship and the Tsinghua Shuimu Scholar Programme. R.M.L. acknowledges the support of NASA through Hubble Fellowship Program grant HST-HF2-51440.001. D.J.K.B. acknowledges support from the Royal Society. J.A.G. acknowledges support from NASA ATP grant 80NSSC20K0540 and from the Alexander von Humboldt Foundation. This work was partially supported under NASA contract No. NNG08FD60C and made use of data from the NuSTAR mission, a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory, and funded by the National Aeronautics and Space Administration. We thank the NuSTAR Operations, Software, and Calibration teams for support with the execution and analysis of these observations. This research has made use of the NuSTAR Data Analysis Software (NuSTARDAS), jointly developed by the ASI Science Data Center (ASDC, Italy) and the California Institute of Technology (USA). We would also like to acknowledge the insightful feedback provided by the reviewer, which significantly improved the quality of this work.

Software: Astropy, DS9, Scipy, Stingray, XSpec.