---
title: Variable Persistent Emission Method
url: https://www.emergentmind.com/topics/variable-persistent-emission-method
type: topic
---

# Variable Persistent Emission Method

Searching arXiv for the core papers on the variable persistent emission method and closely related burst spectroscopy applications.
Variable persistent emission method denotes an X-ray timing and spectral-analysis framework in which the persistent component is not assumed to remain fixed during intervals traditionally treated as burst-only emission. In thermonuclear-burst spectroscopy, the method replaces subtraction of an immutable pre-burst spectrum with a model in which the persistent contribution is multiplied by a free scaling factor \(f_a\) or, in Comptonization-based implementations, by a time-dependent normalization. In magnetar analyses, a related formulation quantifies over-Poisson variability in the persistent light curve through an RMS statistic and interprets that variability with a micro-burst model. Across these uses, the common methodological move is to elevate persistent emission from a fixed background term to an explicitly variable observable [1501.02070][1801.07407].

## 1. Canonical \(f_a\) formulation in thermonuclear-burst spectroscopy

In the burst-spectroscopy formulation, let \(P(E)\) be the best-fit model to the pre-burst persistent spectrum, including Galactic absorption; let \(B(E;T_{\rm bb},K_{\rm bb})\) be a simple blackbody describing the burst emission; and let \(b_{\rm inst}(E)\) denote instrumental background. The variable-persistent-emission model is

$$
S(E)=A(E)\cdot B(E;T_{\rm bb},K_{\rm bb})+f_a\cdot P(E)+b_{\rm inst}(E),
$$

with \(A(E)=\exp[-n_H\,\sigma(E)]\). In shorthand, the method is often written as

$$
\text{total\_spectrum}=f_a\times \text{persistent\_spectrum}+\text{burst\_spectrum}.
$$

Within this formalism, \(f_a\) quantifies whether accretion-powered emission brightens \((f_a>1)\) or dims \((f_a<1)\) during the nuclear flash. The standard approach for time-resolved X-ray spectral analysis of thermonuclear bursts instead subtracts the entire pre-burst emission as background and fits the residual burst spectrum with a single blackbody, thereby assuming both a constant persistent spectral shape and a persistent normalization fixed exactly at the pre-burst level. Introducing \(f_a\) relaxes only the normalization constraint: the shape of \(P(E)\) remains frozen, but its intensity may vary on burst timescales [1501.02070].

This distinction is methodologically important. Under the standard subtraction approach, any real burst-driven change in disk or coronal emission is forced into the burst blackbody fit or into residual structure. Under the \(f_a\) approach, such variability is absorbed into an explicit parameter. A plausible implication is that the method functions simultaneously as a fitting improvement and as a diagnostic of burst–accretion-flow coupling.

## 2. Fitting workflow, parameter control, and statistical validation

The procedure begins with pre-burst modeling. For each burst, a 16 s pre-burst spectrum is accumulated and fitted with a suite of candidate XSPEC models, such as absorbed disk plus Comptonization forms, with the model of minimum \(\chi_\nu^2\) chosen to define \(P(E)\). Time-resolved burst spectroscopy then proceeds by dividing the burst into short intervals, for example \(0.25\) s to a few seconds, while ensuring \(>10^3\) counts. Each interval is fitted with the composite model \(S(E)=A\cdot B(T_{\rm bb},K_{\rm bb})+f_a\cdot P(E)+b_{\rm inst}(E)\), keeping the shape of \(P(E)\) fixed and allowing \(T_{\rm bb}\), \(K_{\rm bb}\), and \(f_a\) to vary freely. A uniform prior is effectively adopted for \(f_a\) over a broad range such as \(-100\) to \(+100\), although physically one typically enforces \(f_a\ge 0\) [1501.02070].

Model selection is performed by comparing goodness-of-fit with and without a free \(f_a\). A typical \(\Delta\chi^2\gg 1\) for one additional degree of freedom indicates significant improvement, and in the RXTE sample the Bayes factor favoring the variable-\(f_a\) approach over the standard method is approximately \(64\), even after penalizing the extra parameter. Procedural validation includes checking that \(f_a\) returns to approximately unity in pre-burst and late-tail intervals. Practical recommendations include extracting a high-quality pre-burst spectrum with at least \(10^3\) counts, fitting multiple plausible persistent models, using sufficiently short time bins to follow rapid \(f_a\) evolution while retaining \(S/N\gtrsim 30\), and adopting thresholds such as \(\Delta\chi^2>4\) for one extra parameter at \(95\%\) confidence when claiming a significant \(f_a\) deviation [1501.02070].

The central statistical point is that the method does not merely add flexibility. It tests a specific null hypothesis—constant persistent normalization—and evaluates whether the data justify replacing that null by a variable normalization while preserving the pre-burst spectral shape.

## 3. Empirical behavior of \(f_a\) and its physical interpretation

A large RXTE reanalysis of \(332\) photospheric radius expansion bursts from \(40\) sources found that, for the majority of spectra, the best-fit value of \(f_a\) is significantly greater than \(1\), indicating that the persistent emission typically increases during a burst. Elevated \(f_a\) values were measured not only during the radius-expansion interval but also in the cooling tail. The modified model yields a lower average value of the \(\chi^2\) fit statistic, although not yet to the level of formal statistical consistency for all spectra. In the same study, an inverse correlation of \(f_a\) with the persistent flux was measured, consistent with theoretical models of disk response [1501.02070].

In the broader burst sample summarized for the method, typical peak values of \(f_a\) in PRE bursts reach approximately \(2\)–\(6\) times the pre-burst level, while in non-PRE bursts they commonly rise by factors of \(2\)–\(4\). The characteristic time profile begins near \(f_a\sim 1\) before the burst, rises during the burst rise, sometimes peaks around photospheric touchdown, and then decays back toward unity in the cooling tail. Non-PRE bursts show a strong positive correlation between instantaneous burst flux \(F_{\rm burst}\) and \(f_a\), with Kendall \(\tau>0\) at \(>3\sigma\) in \(\gtrsim 65\%\) of cases. When normalized by the persistent-to-Eddington ratio \(\gamma\equiv F_{\rm pers}/F_{\rm Edd}\), the product \(f_a\gamma\) is empirically bounded above by \(0.84\) for PRE bursts and \(0.37\) for non-PRE bursts [1501.02070].

The usual physical interpretation is Poynting–Robertson drag. During a bright burst, intense radial photon flux from the neutron-star surface can remove angular momentum from inner-disk material and transiently raise the mass flow onto the star. In this reading, the instantaneous accretion rate is written as

$$
\dot{M}(t)=f_a(t)\,\dot{M}_0,
$$

and simple analytic estimates relate \(f_a-1\sim L_{\rm burst}/L_{\rm Edd}\). At very high burst luminosity approaching \(L_{\rm Edd}\), the disk may be temporarily evacuated or partially disrupted, so the observed spectral signatures can be more complex than a pure normalization change. The method therefore supports, but does not uniquely prove, an accretion-rate interpretation [1501.02070].

## 4. Extensions to superbursts and instrument-specific implementations

During the 2021 superburst of 4U 1820–30, time-resolved spectra from NICER and MAXI were modeled with \(\texttt{tbabs}\times(\texttt{bbodyrad}+\texttt{compTT})\), and in some tail intervals with \(\texttt{tbabs}\times(\texttt{bbodyrad}+\texttt{compTT})\times\texttt{gabs}\). Here \(\texttt{compTT}\) in disk geometry carries free parameters \(kT_0\), \(kT_e\), \(\tau\), and a normalization \(A_C(t)\); the normalization serves as the tracer of variable persistent emission and is identified with the bolometric Comptonization flux \(F_C(t)\). NICER burst-tail spectra were extracted in \(\Delta t=100\) s bins over \(0.5\)–\(10\) keV, with \(1\sigma\) errors derived by the \(\Delta\chi^2=1\) criterion. The recovered persistent flux followed a logistic form,

$$
F_C(t)=\frac{F_{\rm max}}{1+e^{-k(t-t_0)}},
$$

with \(F_{\rm max}=(7.74\pm0.01)\times10^{-9}\ {\rm erg\,s^{-1}\,cm^{-2}}\), \(k=(2.38\pm0.01)\ {\rm hr}^{-1}\), and \(t_0=(4.41\pm0.01)\ {\rm hr}\). The associated \(10\)–\(90\%\) rise time is approximately \(1.8\) hr. The minimum persistent flux at the superburst peak was estimated as \(F_{C,\min}\approx2.2\times10^{-13}\ {\rm erg\,s^{-1}\,cm^{-2}}\), implying \(F_{C,\min}/F_{\rm max}\sim 3\times10^{-5}\), described as nearly complete quenching. Comparison of the superburst total energy, \(\sim3.5\times10^{42}\) erg, with the gravitational binding energy of disk material between \(6\,r_g\) and \(1000\,r_g\), \(\sim6\times10^{39}\) erg, suggests that radiation pressure or Poynting–Robertson drag can evacuate the inner disk; the subsequent recovery timescale is consistent with standard \(\alpha\)-disk viscous times of \(1\)–\(3\) hr for \(\alpha\simeq0.1\)–0.2 and \(R\sim1000\,r_g\) [2412.05785].

The same event also exhibited a transient absorption line that shifted from \(4.15\) to \(3.62\) keV in the \(5\)–\(10\) hr interval. Assigning it to Ar XVIII K\(\alpha\) with rest energy \(E_0=4.15\) keV gives a gravitational redshift

$$
1+z=\frac{E_0}{E_{\rm obs}(t)},
$$

and hence

$$
1+z=\left[1-\frac{2GM}{Rc^2}\right]^{-1/2}.
$$

For \(E_{\rm obs}=3.62\) keV, one obtains \(z\approx0.146\) and \(R\approx17\) km for \(M=1.4\,M_\odot\). The absorption feature was interpreted as likely originating in the inner accretion disk rather than in burst emission from the neutron-star surface, and its evolution suggested inward recovery of the disk [2412.05785].

A NuSTAR study of 4U 1323–62 provides a source-specific \(f_a\) implementation using a pre-burst persistent model \(\texttt{const*TBabs*edge*nthcomp}\) and a burst model \(\texttt{TBabs}\,[f_a\times(\texttt{edge*nthcomp})+\texttt{bbodyrad}]\). In the pre-burst fit, \(N_H\equiv1.2\times10^{22}\ {\rm cm}^{-2}\) was frozen, the absorption edge was found at \(E_{\rm edge}=7.87\pm0.11\) keV with \(\tau_{\rm edge}=0.12\pm0.02\), and the Comptonization parameters were \(\Gamma=1.75\pm0.01\), \(kT_e=32^{+30}_{-8}\) keV, and seed-\(kT_{\rm bb}=1.55\pm0.09\) keV. During burst fits, all persistent-shape parameters were frozen, while only \(f_a\), burst \(kT_{\rm bb}\), and burst normalization were allowed to vary over \(3\)–\(20\) keV. For three bursts divided into five segments \(S1\)–\(S5\) of \(20\) s, except a final \(40\) s segment, all three showed \(f_a\) rising from a pre-burst value near unity to a maximum in \(S2\), then declining toward quiescence. The largest reported enhancement was \(f_a=6.00^{+2.28}_{-2.30}\) in burst B2. In that study, the method improved the fit by \(\Delta\chi^2\approx-136\) with only one extra parameter over a simple blackbody model, recovered average apparent blackbody radii of \(1.5\)–\(3.5\) km, and yielded \(f_a\) values spanning approximately \(1\)–\(8\) [2511.03172].

## 5. RMS-based persistent-emission variability in magnetars

In magnetar work, the expression “variable persistent emission” denotes a different but related methodology. For a background-subtracted light curve with counts \(x_i\), counting errors \(\delta x_i\), and \(N\) bins, the dimensionless RMS intensity variation is defined by

$$
R=\sqrt{\frac{1}{N-1}\left(\sum_i (x_i-\bar{x})^2-\sum_i \delta x_i^2\right)}\Big/\bar{x},
$$

where \(\bar{x}=(1/N)\sum_i x_i\). The subtraction of \(\sum_i\delta x_i^2\) removes the variance expected from counting statistics, so the residual RMS measures intrinsic over-Poisson source variability. Nakagawa et al. proposed that the persistent X-ray emission of magnetars is the superposition of numerous short, \(\lesssim1\) ms micro-bursts of various fluences \(S_c\), with cumulative number–fluence relation

$$
N_c(>S_c)=A_c\,S_c^\alpha,
$$

where \(\alpha\approx -1.1\) from observations of strong bursts, over \(S_1\le S_c\le S_2\) with \(S_1\sim10^{-14}\ {\rm erg\,cm^{-2}}\) and \(S_2\sim10^{-11}\ {\rm erg\,cm^{-2}}\). The associated probability density is

$$
P(S_c)=\frac{-\alpha\,S_c^{\alpha-1}}{S_2^\alpha-S_1^\alpha},
$$

and the expected fractional RMS due solely to micro-burst statistics is

$$
R_M=\sqrt{\sigma_c^2-\sigma_p^2}/S_a.
$$

Inserting \(\alpha=-1.1\), the quoted \(S_1\) and \(S_2\), and a typical persistent X-ray flux of \(10^{-11}\ {\rm erg\,cm^{-2}\,s^{-1}}\), implying \(S_a\approx7.5\times10^{-13}\ {\rm erg\,cm^{-2}}\) in \(0.2\)–\(12\) keV, gives \(R_M\approx14\%\), consistent with observed values of approximately \(5\)–\(20\%\) [1801.07407].

The observational implementation used Suzaku XIS \(0.2\)–\(12\) keV light curves with \(8\) s bins, or \(32\) s for multi-band analysis, and HXD-PIN \(10\)–\(70\) keV light curves with \(128\) s bins. Bright bursts were removed by flagging bins above \(\lambda+5\sigma\), visually verifying them, and recomputing the RMS as \(R'\). Across \(11\) magnetars and \(22\) observations, significant excess RMS intensity variations were found in all \(11\) objects. In four magnetars, corresponding to six observations, the RMS increased clearly toward higher energy bands; in those cases \(R(E)\) rose above the soft–hard crossover of approximately \(4\)–\(8\) keV and tracked the hard power-law component rather than the thermal blackbody. The authors interpreted these results as evidence that persistent emission and burst emission have identical emission mechanisms and that the soft thermal component and hard X-ray component are emitted from different regions far apart from each other. Monte Carlo checks further indicated that spin modulation on \(2\)–\(12\) s periods and day-scale flux drifts do not bias the RMS when binning of at least \(8\) s is used [1801.07407].

## 6. Advantages, limitations, and interpretive boundaries

The principal advantage of the method in burst spectroscopy is improved fit fidelity. Allowing \(f_a\) to vary lowers the average \(\chi^2\) statistic relative to the standard background-subtraction approach, and in the NuSTAR application it removed high-energy residuals while improving the fit by \(\Delta\chi^2\approx-136\) with only one additional parameter. It also yields more reliable blackbody temperatures and radii because the persistent continuum is no longer forced into the burst tail residuals [1501.02070][2511.03172].

Its chief limitation is structural: the method assumes that the *shape* of the persistent spectrum does not change, only its normalization. Rapid coronal cooling, observed above approximately \(30\) keV, or disk ionization changes could violate that assumption. Within the approximate \(2.5\)–\(20\) keV PCA band and on timescales \(\lesssim100\) s, however, no significant shape changes were detected outside bursts. Parameter degeneracy is another persistent concern. \(f_a\) is often strongly covariant with the blackbody normalization, leading to larger uncertainties in both; if the persistent spectrum closely mimics a blackbody, \(f_a\) may be poorly constrained. The method is also sensitive to the quality of the pre-burst fit, so errors in \(N_H\), seed temperature, edge energy, or continuum choice can propagate directly into the inferred burst parameters [1501.02070][2511.03172].

Interpretively, \(f_a\) should not automatically be identified with a pure change in accretion rate. The physical reading \(\dot{M}(t)=f_a(t)\dot{M}_0\) is explicitly conditional on other processes—such as reflection, corona collapse, or more general changes in the Comptonizing medium—being ruled out or modeled. The superburst application strengthens the case that persistent emission can also *decrease* dramatically, since the Comptonization component in 4U 1820–30 was inferred to be nearly completely quenched before recovering on a viscous timescale [2412.05785]. In magnetar work, similarly, excess persistent-emission variability is not attributed to counting noise because the Poisson term is explicitly subtracted and the residual RMS is tested against instrumental and timing-systematics checks [1801.07407].

Taken together, these studies establish variable persistent emission as a methodological category rather than a single code path. In bursting low-mass X-ray binaries it is primarily a spectral-decomposition strategy centered on \(f_a\) or its Comptonization-normalization analogue; in magnetars it is an RMS-based variability formalism tied to a micro-burst hypothesis. What unifies these uses is the rejection of a strictly static view of persistent X-ray emission during high-energy activity.

Source: https://www.emergentmind.com/topics/variable-persistent-emission-method