<figure><img src="image_1.png" /></figure>

# Observations of the Ultra-compact X-Ray Binary 4U 1543-624 in Outburst with NICER, INTEGRAL, Swift, and ATCA

R. M. Ludlam¹, L. Shishkovsky², P. M. Bult³, J. M. Miller¹, A. Zoghbi¹, T. E. Strohmayer³, M. Reynolds¹, L. Natalucci⁴, J. C. A. Miller-Jones⁵, G. K. Jaisawal⁶, S. Guillot⁷,⁸, K. C. Gendreau³, J. A. García⁹,¹⁰, M. Fiocchi⁴, A. C. Fabian¹¹, D. Chakrabarty¹², E. M. Cackett¹³, A. Bahramian⁵, Z. Arzoumanian³, and D. Altamirano¹⁴

¹Department of Astronomy, University of Michigan, 1085 South University Ave, Ann Arbor, MI 48109-1107, USA

²Center for Data Intensive and Time Domain Astronomy, Department of Physics and Astronomy, Michigan State University, East Lansing, MI, USA

³X-ray Astrophysics Laboratory, Astrophysics Science Division, NASA/GSFC, Greenbelt, MD 20771, USA

⁴Istituto Nazionale di Astrofisica, INAF-IAPS, via del Fosso del Cavaliere, 100. I-00133 Roma, Italy

⁵International Centre for Radio Astronomy Research—Curtin University, GPO Box U1987, Perth, WA 6845, Australia

⁶National Space Institute, Technical University of Denmark, Elektrovej 327-328, DK-2800 Lyngby, Denmark

⁷CNRS, IRAP, 9 avenue du Colonel Roche, BP 44346, F-31028 Toulouse Cedex 4, France

⁸Université de Toulouse, CNES, UPS-OMP, F-31028 Toulouse, France

⁹Cahill Center for Astronomy and Astrophysics, California Institute of Technology, Pasadena, CA 91125, USA

¹⁰Remeis Observatory & ECAP, Universität Erlangen-Nürnberg, Sternwartstr. 7, D-96049 Bamberg, Germany

¹¹Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK

¹²MIT Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA

¹³Department of Physics & Astronomy, Wayne State University, 666 West Hancock Street, Detroit, MI 48201, USA

¹⁴Department of Physics & Astronomy, University of Southampton, Highfield, Southampton SO17 1BJ, UK

Received 2019 June 18; revised 2019 July 26; accepted 2019 July 31; published 2019 September 18

## Abstract

We report on X-ray and radio observations of the ultra-compact X-ray binary 4U 1543-624 taken in August 2017 during an enhanced accretion episode. We obtained Neutron Star Interior Composition Explorer (NICER) monitoring of the source over a ~10 day period during which target-of-opportunity observations were also conducted with Swift, INTEGRAL, and the Australia Telescope Compact Array. Emission lines were measured in the NICER X-ray spectrum at ~0.64 keV and ~6.4 keV that correspond to O and Fe, respectively. By modeling these line components, we are able to track changes in the accretion disk throughout this period. The innermost accretion flow appears to move inwards from hundreds of gravitational radii ($R_g = GM/c^2$) at the beginning of the outburst to <8.7 $R_g$ at peak intensity. We do not detect the source in radio, but are able to place a $3\sigma$ upper limit on the flux density at 27 $\mu$Jy beam$^{-1}$. Comparing the radio and X-ray luminosities, we find that the source lies significantly away from the range typical of black holes in the $L_r-L_x$ plane, suggesting a neutron star primary. This adds to the evidence that neutron stars (NSs) do not follow a single track in the $L_r-L_x$ plane, limiting its use in distinguishing between different classes of NSs based on radio and X-ray observations alone.

**Key words:** accretion, accretion disks – stars: neutron

## 1. Introduction

Ultra-compact X-ray binaries (UCXBs) are low-mass X-ray binaries (LMXBs) with a short orbital period of <80 minutes, and typically contain a degenerate stellar companion. 4U 1543-624 is an UCXB with an orbital period of $18.2 \pm 0.1$ minutes determined from optical photometry (Wang & Chakrabarty 2004; Wang et al. 2015). The source was first detected by the *Uhuru* mission (Giacconi et al. 1972). The absence of hydrogen and helium lines, coupled with emission from carbon and oxygen in the optical spectrum of 4U 1543-624, indicates a degenerate donor star such as a carbon-oxygen white dwarf (Nelemans et al. 2003).

The nature of the compact accretor in this system is unknown, but is often assumed to be a neutron star (NS). The Monitor of All-sky X-ray Image (MAXI) recently detected a Type-I X-ray burst (Serino et al. 2018) consistent with the position of 4U 1543-624,<fnref n="15" /> although the instrument's spatial resolution of 1.⁵° (Matsuoka et al. 2009) prevents a definitive conclusion regarding the event's origin because there are other X-ray sources nearby. If the burst indeed arose from

4U 1543-624, then this would irrefutably confirm that the system contains an NS. Additionally, this would place 4U 1543-624 within a subset of unique UCXBs that show helium-powered X-ray bursts even though they have degenerate C/O donors and their optical spectra lack evidence of He (4U 0614+091: Kuulkers et al. 2010, IGR J17062-6143: Strohmayer et al. 2018b).

Further support that 4U 1543-624 contains a NS can be obtained from observing the source in both radio and X-rays. For LMXBs, the general view is that the X-ray emission tracks the accretion inflow closest to the compact object, while radio emission represents an outflow in the form of a relativistic jet (Fender et al. 2003). The exact mechanism creating these jets is not confirmed and may vary with different systemic properties, but is largely accepted to involve the acceleration of material by the magnetic field of the compact object or inner accretion disk, creating the synchrotron radio emission that has been observed in many LMXBs (e.g., Blandford & Znajek 1977; Blandford & Payne 1982; Tudor et al. 2017). Observing the radio and X-ray luminosity simultaneously can be advantageous in discerning between an NS or black hole (BH) accretor, as NSs are typically ~30 times fainter in the radio than BHs at the same X-ray luminosity (Migliari & Fender 2006). However,

<footnote n="15">¹⁵ http://maxi.riken.jp/alert/novae/8151027492/8151027492.htm</footnote>
some caution should be maintained in utilizing this diagnostic tool given that certain classes of NSs, namely transitional and accreting millisecond X-ray pulsars (AMXPs), can be nearly as radio-luminous as the radio-faint BHs (Deller et al. 2015).

In many LMXBs, hard X-ray photons incident on the surrounding material in the accretion disk will result in emission features from the photons being reprocessed. These emission features are broadened by strong Doppler, special relativistic, and general relativistic effects from proximity to the compact accretor (Fabian et al. 2000). The strength of these effects can therefore be used to determine how close the accretion disk is to the NS or BH.

4U 1543−624 was first reported to exhibit a broad O VIII Lyα emission feature at ~0.7 keV by Madej & Jonker (2011) when observed with the high-resolution spectrographs onboard *XMM-Newton* and *Chandra*. Previous evidence for an emission feature near 0.7 keV existed in ASCA and *BeppoSAX* observations (Juett et al. 2001; Schultz 2003), but these early spectra could also be described without invoking O emission if there was an overabundance in neon that increased absorption along the line of sight (Juett et al. 2001). This generated uncertainty surrounding the identification of this feature. A similar O feature was reported with *Chandra* and *XMM-Newton* by Juett & Chakrabarty (2003), but early calibration uncertainties hindered a definitive claim.

In addition to O VIII Lyα, there was also evidence of a possible Fe Kα feature at ~6.6 keV concurrent with the detection of O VIII Lyα (Madej & Jonker 2011). This was not the first time that the Fe line had been detected in this system, because it was present in prior *RXTE* and *EXOSAT* data (Schultz 2003), but Madej & Jonker (2011) were the first to claim that both features were present in the X-ray spectrum. The O VIII feature was found to be more intense relative to the continuum than Fe K. For a typical accretion disk of solar abundance, O VIII is the second strongest feature after Fe K. In order for O to be the most prominent feature, the disk would need a significant overabundance of oxygen (Ballantyne et al. 2002), further supporting that the donor star in the system is a CO or ONe white dwarf (Madej & Jonker 2011).

Madej et al. (2014) presented an X-ray spectral analysis using a preliminary version of a new reflection model tailored to accommodate the typical elemental abundances in UCXBs known as XILLVER<sub>CO</sub>. This model has negligible H and He abundances, variable abundance of C and O, and 10 times solar abundance for all other elemental species, but had only a limited number of grid points (i.e., large steps between parameter values) at the time. Spectral modeling using this initial XILLVER<sub>CO</sub> grid on 4U 1543−624 indicated an inclination of *i* ~ 65° and an inner disk radius <7.4 *R*<sub>g</sub> (where *R*<sub>g</sub> = *GM*/*c*<sup>2</sup>).

In early 2017 August, the *Swift/Burst Alert Telescope* (BAT) detected increased activity from 4U 1543−624 in the 15–50 keV energy band (Ludlam et al. 2017). We requested *Neutron Star Interior Composition Explorer (NICER)* monitoring of the source in order to track any changes in the system as the X-ray intensity increased. Additionally, we secured target-of-opportunity observations with *Swift, International Gamma-Ray Astrophysics Laboratory (INTEGRAL)*, and the Australia Telescope Compact Array (ATCA). We present the results of our X-ray observations throughout this period of enhanced accretion and reanalyze the radio data that was initially reported in Ludlam et al. (2017) below.

## 2. Observations and Data Reduction

### 2.1. NICER

NICER (Gendreau et al. 2012) onboard the International Space Station is comprised of 56 “concentrator” optics and silicon drift detector pairs that each collect X-rays in the 0.2–12 keV range. The 52 operational detectors provide a collecting area of ~1900 cm<sup>2</sup> at 1.5 keV. NICER observed 4U 1543−624 over a ~10 day period in 2017 August (ObsIDs 1050060101-1050060113) for a cumulative exposure time of ~86.2 ks. Data were reduced using NICERDAS 2018-10-07_V005. We created good time intervals (GTIs) using NIMAKETIME to select events that occurred when the space weather index KP < 5 to avoid periods of time with high particle background due to geomagnetic storms and magnetic cutoff rigidity of COR_SAX > 4 to ensure that high particle radiation intervals coincident with the Earth’s auroral zones are removed. These GTIs were applied to the data using NIEX-TRACT-EVENTS selecting pulse-invariant (PI) energy channels between 25 and 1200, inclusive, and EVENT_FLAGS=bxxx1x000.

The event files for each observation were loaded into XSELECT and combined to create light curves and spectra. Light curves were extracted in five different energy bands (super soft: 0.5–1.1 keV, soft: 1.1–2.0 keV, intermediate: 2.0–3.8 keV, hard: 3.8–6.8 keV, and full: 0.5–6.8 keV) and binned to 128 s for color analysis as per Bult et al. (2018). Figure 1 shows the NICER light curve, soft and hard color evolution, color intensity, and color–color diagrams. We extract five spectra based on changes in the light curve and color diagram, which are labeled A–E. These will be referred to as intervals A–E in later sections. The spectra are normalized to instrumental residuals created from the Crab Nebula (see Ludlam et al. 2018 for more details). Backgrounds were generated based upon the same filtering criteria using *RXTE* “blank sky” field 5 (Jahoda et al. 2006). The cleaned source spectra have exposure times of 11.4, 8.45, 10.3, 10.2, and 7.96 ks for intervals A-E, respectively, after filtering. We use the publicly available ARF and RMF instrument response files from 2017 June 1 for these observations.

### 2.2. Swift

4U 1543-624 has been observed on numerous occasions by the *Neil Gehrels Swift Observatory (Swift)*. There was a single X-Ray Telescope (XRT) observation (ObsID: 00010238010) in Windowed Timing mode on 2017 August 18 that overlapped with the NICER data during interval A. The raw *Swift* observation files were downloaded from HEASARC and reprocessed using *XRTPIPELINE* and *CALDB* version 20180710. Although the observation is short (~200 s), the source is detected with a net count rate of 18.05 ± 0.32 ct s⁻¹, which is well below the threshold for photon pile-up in Windowed Timing mode.

### 2.3. INTEGRAL

The INTEGRAL observations occurred during interval E of the NICER monitoring for 23.9 ks. The INTEGRAL/Imager on-Board the INTEGRAL Satellite (IBIS; Ubertini et al. 2003) data for this observation were processed using the standard off-line Scientific Analysis (OSA v10.2) software released by the INTEGRAL Scientific Data Centre (Courvoisier et al. 2003) in
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Panel (a): the 0.5–6.8 keV NICER light curve of 4U 1543−624 starting on 2017 August 16. The vertical dashed purple lines indicate the time intervals in which spectra were extracted. Data were binned to 128 s. Panel (b): the soft color (1.1–2.0 keV/0.5–1.1 keV) evolution of the source. Panel (c): the hard color (3.8–6.8 keV/2.0–3.8 keV) evolution throughout the outburst. Panels (d) and (e): soft color and hard color intensity diagrams. Panel (f): soft color vs. hard color diagram. The single outlier point in interval A is an artifact of the 128 s time binning.</figcaption></figure>

order to obtain a spectrum. These runs were performed with the AVES cluster, designed to optimize performance and disk storage needed for the INTEGRAL data analysis by Federici et al. (2010).

## 2.4. ATCA

4U 1543-624 was observed by the ATCA on 2017 August 23, from 11:24:59.9 to 16:09:19.9 UTC, with ≈3.25 hr on source (project code CX392). This falls within interval E of the NICER monitoring. The 4 cm CABB receiver was used, which was set up with two frequency subbands of 2048 MHz bandwidth each, centered at 5.5 and 9.0 GHz. The ATCA array was in the 1.5A configuration, which has a maximum baseline length of 4.47 km. The flux and bandpass calibrator used for this observation was 1934-638, while 1554-64 was used as a phase calibrator.

Data for each frequency band were separately reduced according to standard procedures with Miriad (Sault et al. 1995). The target visibilities were then imported into the Common Astronomy Software Application (CASA; McMullin et al. 2007) to be imaged. Each 5.5 and 9.0 GHz baseband was imaged separately, with an additional image created combining both data sets in the uv-plane having a central frequency of 7.25 GHz. We used the CASA task `tclean`, selecting Briggs weighting with a robust parameter of 1, as well as nterms =2 to account for non-zero spectral indices in other field sources.

The synthesized beams at 5.5 GHz, 7.25 GHz, and 9.0 GHz were $6.3^\prime \times 2.1^\prime$, $5.4^\prime \times 1.8^\prime$, and $4.1^\prime \times 1.4^\prime$, respectively. The local sensitivities achieved were: 13 µJy beam⁻¹ at 5.5 GHz; 9 µJy beam⁻¹ at 7.25 GHz; 10 µJy beam⁻¹ at 9.0 GHz. These are slightly higher than the expected theoretical sensitivities, due to several bright field sources in the primary beam side-lobes. In any case, 4U 1543−624 is clearly not detected in any of our images. We derive a 3σ upper limit

on the radio luminosity of 4U 1543−624 using the corresponding 3σ sensitivity of our deepest image (7.25 GHz): 27 µJy beam⁻¹.

# 3. X-Ray Analysis and Results

## 3.1. Spectral

We use XSPEC (Arnaud 1996) version 12.10.1 for our spectral analysis. Parameter uncertainties are reported at the 90% confidence level. These are determined from Monte Carlo Markov Chains of length $10^6$ with an equal burn-in length using the “chain” command. NICER spectra are modeled in the 0.4–9.5 keV band because the effective collecting area drops sharply outside of this region (Ludlam et al. 2018). The Swift/XRT data are considered in the 0.3–8.5 keV band because the spectrum becomes quickly background-dominated above this energy. The INTEGRAL/IBIS observation is fit in the 26–100 keV band. In order to account for differences in calibration between missions, a multiplicative constant is allowed to float between the Swift and NICER data in interval A and INTEGRAL and NICER during interval E.

We use TBABS to account for the neutral absorption column along the line of sight, which is tied between intervals. We concurrently fit the continuum for each spectrum using a variety of models that have been used previously for 4U 1543−624. The parameter values and normalization of the continuum components are allowed to vary in order to account for changes in the spectral shape as the source intensity increases. Using an absorbed double thermal component model, TBABS\*(DISKBB+BBODY), similar to Ng et al. (2010), provides a reduced $\chi^2/\nu = 3.51$ and fails to model the INTEGRAL component in interval E. If instead we use an absorbed power-law and blackbody model as per Juett & Chakrabarty (2003), this improves the fit to $\chi^2/\nu = 2.11$. If
instead we use a Comptonization component, NTHCOMP, in place of the power law similar to Schultz (2003), we obtain a comparable fit of $\chi_{\nu}^{2} = 2.11$ from TBABS*(NTHCOMP+BBODY), assuming that the seed photons originate from the accretion disk. Conversely, if we assume the seed photons originate from the blackbody component, TBABS*(NTHCOMP+DISKBB), then the fit becomes worse ($\chi_{\nu}^{2} = 3.66$). An absorbed cutoff power-law and blackbody model from Madej & Jonker (2011) further improves the overall fit to $\chi_{\nu}^{2} = 2.07$.

The model that is able to best describe the NICER and INTEGRAL data is the hybrid model for NSs developed by Lin et al. (2007). This is composed of a multi-temperature blackbody (DISKBB) for the accretion disk, a single-temperature blackbody (BBODY) for a boundary layer or emission from the surface of the NS, and a power-law component (POWERLAW) to account for emission from a coronal region. We apply the same continuum description to each interval given that the hard color in Figure 1 shows very little change throughout the outburst. This continuum model provides a reduced $\chi_{\nu}^{2}$ of 1.78. Table 1 provides the parameter values for continuum fitting (Model 1) for intervals A–E. Each of these components are statistically needed at >27σ level of confidence via an F-test.

Switching out the single-temperature thermal component for Comptonization, NTHCOMP, provided a significantly worse fit ($\chi_{\nu}^{2} = 2.35$). Moreover, the Comptonized component tended toward the shape of a blackbody with a high optical depth of $\tau \sim 9$. Alternatively, if we switch the power law out for the Comptonization component, still assuming the seed photons originate from the single-temperature blackbody, we obtain a marginally better fit over Model 1 in Table 1 ($\chi_{\nu}^{2} = 1.75$). However, the electron temperature tends to the hard limit of $kT_{e} = 1000$ keV in all cases in order to properly describe the INTEGRAL data. Fixing $kT_{e}$ at lower values provides a worse $\chi^{2}$ and fails to fit the high-energy data from the INTEGRAL observation. For this reason, we proceed with our analysis using the continuum description of Lin et al. (2007).

There is an excess in emission in the NICER spectra between 0.6–0.7 keV and 6.0–7.0 keV that can be attributed to the O and Fe K emission lines that have previously been detected in 4U 1543–624. We do not detect any emission features in the Swift/XRT spectrum due to the low signal-to-noise from the short exposure and modest count rate. Figure 2 shows the ratio of the NICER data to the respective continuum model for intervals A–E. A clear evolution in the oxygen line profile with time can be seen, whereas any change in the Fe line profile between intervals is not significantly detected. This may be due to there being three times less NICER collecting area in the Fe K band than in the O band. In Figure 3, we show the change in the O line profile throughout the outburst by dividing the spectrum in each interval by the time-averaged spectrum. There is a deficit of emission (with respect to the time-averaged profile) between 0.7 and 0.8 keV that corresponds to the apparent O VII edge in interval A that disappears as the source intensity increases. This is in agreement with the behavior shown in Figure 2. The disappearance of the O VII edge may be due to a change in the ionization state of the material with time.

To assess the broadening of the features between 0.6–0.7 keV and 6.0–7.0 keV, we employ simple Gaussian components. The line widths are between $1.6 \times 10^{-2}$ keV < $\sigma_{O} < 10.2 \times 10^{-2}$ keV and $0.74$ keV < $\sigma_{Fe} < 1.89$ keV, respectively. The O line width agrees with the values

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Ratio of NICER data to continuum model in each interval for the O band (top) and Fe K band (bottom). The continuum is modeled by an absorbed multi-temperature accretion disk, single-temperature blackbody, and power law. The apparent absorption feature between 0.7 and 0.8 keV in interval A is likely due to an O VII edge that disappears as the ionization state of the material changes (see text for more detail).</figcaption></figure>

previously reported in Madej & Jonker (2011), while the Fe component values are typical of the broadening seen in other accreting NS LMXBs (e.g., Cackett et al. 2012). In order to obtain physical constraints from these features, we add two DISKLINE (Fabian et al. 1989) components to model the emission line broadening due to accretion around a compact object with dimensionless spin of $a = 0$ (where $a = cJ/GM^{2}$). The inclination (i) and emissivity index (q) are tied between all line components for all intervals. The outer disk radius is fixed at 1000 $R_{g}$. The line energy, inner disk radius, and normalization for each line are allowed to vary between intervals in order to track any qualitative changes in the accretion disk. This provides a significant improvement in the overall fit ($\Delta\chi^{2} = 2402.95$ for 32 dof). Values for each parameter are presented under Model 2 in Table 1. Figure 4 shows the broadband spectral model and residuals for interval E. This is representative of the overall model applied to each interval, but we only show one epoch for clarity.

The line energy for O is slightly lower than the laboratory value for O VIII Lyα (0.654 keV), but this could be due to a blending with O VII emission as NICER has an energy resolution of ~80 eV at 1 keV and/or partially due to the
<table><caption>**Table 1**

Spectral Modeling of X-Ray Observations</caption><thead><tr><th>Model 1</th><th>Parameter</th><th>Interval A</th><th>Interval B</th><th>Interval C</th><th>Interval D</th><th>Interval E</th></tr></thead><tbody><tr><td>TBABS</td><td>N<sub>H</sub> (10<sup>21</sup> cm<sup>-2</sup>)</td><td></td><td></td><td>2.92<sup>+0.02</sup><sub>-0.03</sub></td><td></td><td></td></tr><tr><td>CONSTANT</td><td>C</td><td>0.93 ± 0.03</td><td>...</td><td>...</td><td>...</td><td>0.33 ± 0.03</td></tr><tr><td>DISKBB</td><td>kT (10<sup>-2</sup> keV)</td><td>4.9<sup>+0.3</sup><sub>-0.4</sub></td><td>8.7<sup>+2.8</sup><sub>-0.8</sub></td><td>9.3<sup>+0.4</sup><sub>-0.6</sub></td><td>12.2<sup>+0.4</sup><sub>-0.3</sub></td><td>16.1 ± 0.3</td></tr><tr><td></td><td>norm<sub>disk</sub> (10<sup>5</sup>)</td><td>150<sup>+222</sup><sub>-70</sub></td><td>1.9 ± 1.5</td><td>1.2<sup>+0.6</sup><sub>-0.3</sub></td><td>0.36 ± 0.06</td><td>0.20 ± 0.02</td></tr><tr><td>BBODY</td><td>kT (keV)</td><td>0.66 ± 0.01</td><td>0.67<sup>+0.02</sup><sub>-0.04</sub></td><td>0.73 ± 0.01</td><td>0.74<sup>+0.01</sup><sub>-0.02</sub></td><td>0.65<sup>+0.02</sup><sub>-0.01</sub></td></tr><tr><td></td><td>norm<sub>bb</sub> (10<sup>-3</sup>)</td><td>1.04<sup>+0.04</sup><sub>-0.03</sub></td><td>0.89<sup>+0.05</sup><sub>-0.03</sub></td><td>1.06 ± 0.04</td><td>1.17 ± 0.04</td><td>1.22 ± 0.06</td></tr><tr><td>POWERLAW</td><td>Γ</td><td>1.74 ± 0.01</td><td>1.73<sup>+0.01</sup><sub>-0.02</sub></td><td>1.77 ± 0.01</td><td>1.79 ± 0.01</td><td>1.71 ± 0.01</td></tr><tr><td></td><td>norm<sub>pl</sub> (10<sup>-1</sup>)</td><td>1.46 ± 0.01</td><td>1.58<sup>+0.01</sup><sub>-0.03</sub></td><td>1.61 ± 0.01</td><td>1.66<sup>+0.01</sup><sub>-0.02</sub></td><td>1.59 ± 0.03</td></tr><tr><td></td><td>χ<sup>2</sup> (dof)</td><td></td><td></td><td>8306.01 (4668)</td><td></td><td></td></tr><tr><td>Model 2</td><td>Parameter</td><td>Interval A</td><td>Interval B</td><td>Interval C</td><td>Interval D</td><td>Interval E</td></tr><tr><td>TBABS</td><td>N<sub>H</sub> (10<sup>21</sup> cm<sup>-2</sup>)</td><td></td><td></td><td>3.06<sup>+0.03</sup><sub>-0.02</sub></td><td></td><td></td></tr><tr><td>CONSTANT</td><td>C</td><td>0.93 ± 0.03</td><td>...</td><td>...</td><td>...</td><td>0.44<sup>+0.05</sup><sub>-0.04</sub></td></tr><tr><td>DISKBB</td><td>kT (10<sup>-2</sup> keV)</td><td>4.6 ± 0.2</td><td>5.7<sup>+0.3</sup><sub>-0.2</sub></td><td>5.6 ± 0.2</td><td>5.7 ± 0.2</td><td>10.6<sup>+0.2</sup><sub>-0.6</sub></td></tr><tr><td></td><td>norm<sub>disk</sub> (10<sup>5</sup>)</td><td>456<sup>+278</sup><sub>-150</sub></td><td>68<sup>+25</sup><sub>-24</sub></td><td>88<sup>+36</sup><sub>-19</sub></td><td>81<sup>+24</sup><sub>-17</sub></td><td>1.4<sup>+0.5</sup><sub>-0.1</sub></td></tr><tr><td>BBODY</td><td>kT (keV)</td><td>0.69 ± 0.01</td><td>0.74 ± 0.02</td><td>0.80<sup>+0.02</sup><sub>-0.01</sub></td><td>0.83<sup>+0.02</sup><sub>-0.01</sub></td><td>0.81<sup>+0.04</sup><sub>-0.02</sub></td></tr><tr><td></td><td>norm<sub>bb</sub> (10<sup>-3</sup>)</td><td>1.07<sup>+0.03</sup><sub>-0.05</sub></td><td>0.94<sup>+0.06</sup><sub>-0.09</sub></td><td>1.24<sup>+0.06</sup><sub>-0.08</sub></td><td>1.36<sup>+0.09</sup><sub>-0.07</sub></td><td>0.98<sup>+0.08</sup><sub>-0.09</sub></td></tr><tr><td>POWERLAW</td><td>Γ</td><td>1.77 ± 0.01</td><td>1.78 ± 0.01</td><td>1.83 ± 0.01</td><td>1.86 ± 0.01</td><td>1.82 ± 0.01</td></tr><tr><td></td><td>norm<sub>pl</sub> (10<sup>-1</sup>)</td><td>1.51<sup>+0.02</sup><sub>-0.01</sub></td><td>1.78<sup>+0.02</sup><sub>-0.01</sub></td><td>1.66 ± 0.01</td><td>1.74 ± 0.01</td><td>1.86<sup>+0.02</sup><sub>-0.01</sub></td></tr><tr><td>DISKLINE<sub>1</sub></td><td>E<sub>O</sub> (keV)</td><td>0.637<sup>+0.004</sup><sub>-0.009</sub></td><td>0.637<sup>+0.010</sup><sub>-0.008</sub></td><td>0.646<sup>+0.008</sup><sub>-0.010</sub></td><td>0.632<sup>+0.006</sup><sub>-0.008</sub></td><td>0.659<sup>+0.005</sup><sub>-0.008</sub></td></tr><tr><td></td><td>q</td><td></td><td></td><td>-2.36<sup>+0.05</sup><sub>-0.07</sub></td><td></td><td></td></tr><tr><td></td><td>R<sub>in</sub> (R<sub>g</sub>)</td><td>942<sup>+30</sup><sub>-754</sub></td><td>46<sup>+65</sup><sub>-19</sub></td><td>44<sup>+39</sup><sub>-19</sub></td><td>7.9<sup>+2.5</sup><sub>-1.3</sub></td><td>7.3<sup>+1.4</sup><sub>-1.0</sub></td></tr><tr><td></td><td>i (°)</td><td></td><td></td><td>74<sup>+6</sup><sub>-4</sub></td><td></td><td></td></tr><tr><td></td><td>norm<sub>line</sub> (10<sup>-2</sup>)</td><td>0.19<sup>+0.04</sup><sub>-0.02</sub></td><td>0.40<sup>+0.08</sup><sub>-0.07</sub></td><td>0.45<sup>+0.10</sup><sub>-0.06</sub></td><td>1.2 ± 0.1</td><td>1.5<sup>+0.2</sup><sub>-0.1</sub></td></tr><tr><td></td><td>EW<sub>O</sub> (eV)</td><td>5.0<sup>+1.1</sup><sub>-0.5</sub></td><td>10<sup>+2.0</sup><sub>-1.75</sub></td><td>11<sup>+2.4</sup><sub>-1.5</sub></td><td>27 ± 2.3</td><td>33<sup>+4.4</sup><sub>-2.2</sub></td></tr><tr><td>DISKLINE<sub>2</sub></td><td>E<sub>Fe</sub> (keV)</td><td>6.42<sup>+0.28</sup><sub>-0.01</sub></td><td>6.7 ± 0.2</td><td>6.41<sup>+0.13</sup><sub>-0.01</sub></td><td>6.41<sup>+0.06</sup><sub>-0.01</sub></td><td>6.48<sup>+0.27</sup><sub>-0.08</sub></td></tr><tr><td></td><td>R<sub>in</sub> (R<sub>g</sub>)</td><td>28<sup>+26</sup><sub>-2</sub></td><td>17 ± 11</td><td>9<sup>+9</sup><sub>-2</sub></td><td>12<sup>+7</sup><sub>-5</sub></td><td>10.2<sup>+12.5</sup><sub>-3.4</sub></td></tr><tr><td></td><td>norm<sub>line</sub> (10<sup>-3</sup>)</td><td>0.6<sup>+0.1</sup><sub>-0.2</sub></td><td>0.9<sup>+0.3</sup><sub>-0.5</sub></td><td>1.2 ± 0.3</td><td>1.3 ± 0.3</td><td>0.8 ± 0.3</td></tr><tr><td></td><td>EW<sub>Fe</sub> (eV)</td><td>105<sup>+18</sup><sub>-35</sub></td><td>160<sup>+53</sup><sub>-89</sub></td><td>210 ± 53</td><td>228 ± 53</td><td>123 ± 46</td></tr><tr><td></td><td>χ<sup>2</sup> (dof)</td><td></td><td></td><td>5903.06 (4636)</td><td></td><td></td></tr><tr><td>Model 3</td><td>Parameter</td><td>Interval A</td><td>Interval B</td><td>Interval C</td><td>Interval D</td><td>Interval E</td></tr><tr><td>TBABS</td><td>N<sub>H</sub> (10<sup>21</sup> cm<sup>-2</sup>)</td><td></td><td></td><td>3.05<sup>+0.04</sup><sub>-0.02</sub></td><td></td><td></td></tr><tr><td>CONSTANT</td><td>C</td><td>0.93 ± 0.03</td><td>...</td><td>...</td><td>...</td><td>0.45<sup>+0.05</sup><sub>-0.04</sub></td></tr><tr><td>DISKBB</td><td>kT (10<sup>-2</sup> keV)</td><td>4.7<sup>+0.1</sup><sub>-0.2</sub></td><td>5.7<sup>+0.6</sup><sub>-0.4</sub></td><td>5.5<sup>+0.1</sup><sub>-0.3</sub></td><td>5.6 ± 0.2</td><td>10.1<sup>+0.2</sup><sub>-0.4</sub></td></tr><tr><td></td><td>norm<sub>disk</sub> (10<sup>5</sup>)</td><td>400<sup>+200</sup><sub>-90</sub></td><td>74<sup>+74</sup><sub>-43</sub></td><td>102<sup>+77</sup><sub>-22</sub></td><td>97<sup>+47</sup><sub>-26</sub></td><td>1.74<sup>+0.5</sup><sub>-0.2</sub></td></tr><tr><td>BBODY</td><td>kT (keV)</td><td>0.68<sup>+0.02</sup><sub>-0.01</sub></td><td>0.72 ± 0.02</td><td>0.79<sup>+0.02</sup><sub>-0.01</sub></td><td>0.83<sup>+2</sup><sub>-0.01</sub></td><td>0.81<sup>+0.03</sup><sub>-0.01</sub></td></tr><tr><td></td><td>norm<sub>bb</sub> (10<sup>-3</sup>)</td><td>1.00<sup>+0.04</sup><sub>-0.03</sub></td><td>0.88<sup>+0.05</sup><sub>-0.06</sub></td><td>1.17<sup>+0.07</sup><sub>-0.04</sub></td><td>1.40<sup>+0.8</sup><sub>-0.06</sub></td><td>1.00<sup>+0.08</sup><sub>-0.07</sub></td></tr><tr><td>POWERLAW</td><td>Γ</td><td>1.76 ± 0.01</td><td>1.77 ± 0.01</td><td>1.82 ± 0.01</td><td>1.87 ± 0.01</td><td>1.83 ± 0.01</td></tr><tr><td></td><td>norm<sub>pl</sub> (10<sup>-1</sup>)</td><td>1.51<sup>+0.02</sup><sub>-0.01</sub></td><td>1.64<sup>+0.02</sup><sub>-0.01</sub></td><td>1.66<sup>+0.02</sup><sub>-0.01</sub></td><td>1.74<sup>+0.02</sup><sub>-0.01</sub></td><td>1.87<sup>+0.02</sup><sub>-0.01</sub></td></tr><tr><td>DISKLINE<sub>1</sub></td><td>E<sub>O</sub> (keV)</td><td>0.637<sup>+0.003</sup><sub>-0.017</sub></td><td>0.63<sup>+0.01</sup><sub>-0.02</sub></td><td>0.63 ± 0.01</td><td>0.63 ± 0.01</td><td>0.66 ± 0.01</td></tr><tr><td></td><td>q</td><td></td><td></td><td>-2.34<sup>+0.05</sup><sub>-0.07</sub></td><td></td><td></td></tr><tr><td></td><td>R<sub>in</sub> (R<sub>g</sub>)</td><td>515 ± 455</td><td>32<sup>+22</sup><sub>-16</sub></td><td>18<sup>+11</sup><sub>-1</sub></td><td>6.94<sup>+3.0</sup><sub>-9.4</sub></td><td>6.1<sup>+2.1</sup><sub>-0.1</sub></td></tr><tr><td></td><td>i (°)</td><td></td><td></td><td>70<sup>+10</sup><sub>-1</sub></td><td></td><td></td></tr><tr><td></td><td>norm<sub>line</sub> (10<sup>-2</sup>)</td><td>0.20<sup>+0.07</sup><sub>-0.02</sub></td><td>0.4 ± 0.1</td><td>0.55<sup>+0.09</sup><sub>-0.06</sub></td><td>1.2 ± 0.1</td><td>1.6 ± 0.1</td></tr><tr><td></td><td>EW<sub>O</sub> (eV)</td><td>5.5<sup>+1.9</sup><sub>-0.5</sub></td><td>11 ± 3</td><td>13.3<sup>+2.2</sup><sub>-1.4</sub></td><td>27 ± 2</td><td>35 ± 2</td></tr><tr><td>DISKLINE<sub>2</sub></td><td>E<sub>Fe</sub> (keV)</td><td>6.41<sup>+0.28</sup><sub>-0.01</sub></td><td>6.7 ± 0.2</td><td>6.41<sup>+0.20</sup><sub>-0.01</sub></td><td>6.41<sup>+0.05</sup><sub>-0.01</sub></td><td>6.41<sup>+0.21</sup><sub>-0.01</sub></td></tr><tr><td></td><td>norm<sub>line</sub> (10<sup>-3</sup>)</td><td>0.2<sup>+0.3</sup><sub>-0.1</sub></td><td>0.4<sup>+0.2</sup><sub>-0.1</sub></td><td>0.8 ± 0.2</td><td>1.5 ± 0.2</td><td>0.9 ± 0.3</td></tr><tr><td></td><td>EW<sub>Fe</sub> (eV)</td><td>26<sup>+39</sup><sub>-13</sub></td><td>75<sup>+37</sup><sub>-19</sub></td><td>142 ± 36</td><td>267 ± 36</td><td>143 ± 48</td></tr><tr><td></td><td>χ<sup>2</sup> (dof)</td><td></td><td></td><td>5936.04 (4641)</td><td></td><td></td></tr></tbody></table>

Note. Errors are reported at the 90% confidence level and calculated from Markov Chain Monte Carlo of chain length 10⁶. NICER is fit in the 0.4–9.5 keV energy band, INTEGRAL is fit in the 26–100 keV band, and Swift is considered in the 0.3–8.5 keV band. The χ² (dof) reported refers to the simultaneous fitting of all spectra. The multiplicative constant is used on the Swift data in Interval A and INTEGRAL data in Interval E with all other parameters tied to the NICER data in that respective interval. The outer disk radius is fixed at 1000 R_g. Inclination (i) and emissivity index (q) are tied between the two disk line components. R_in is tied between each DISKLINE component in Model 3.

current uncertainty in the gain (<10 eV). The Fe line component is consistent with a blend of Fe XXV with lower ionization states in intervals A, B, and potentially E. Intervals C

and D are inconsistent with He-like Fe XXV. This change in the ionization state of Fe line may be due to the spectral softening of the power-law component with time. This can lead to less
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Ratio of the NICER spectrum in each interval normalized to the time average spectrum. This demonstrates a clear change in the O line profile that is independent of the calibration correction. The teal point indicates the errors in the *x* and *y* directions. Values larger than unity indicate that emission is stronger in that region, whereas values less than unity indicate that the emission is weaker than the mean O profile.</figcaption></figure>

Fe-ionizing flux in the harder energy band. The O line component requires less energy to become ionized and therefore remains consistent with the H-like charge state throughout the outburst with some contribution from He-like O early on.

The inner disk radius appears to move inward as the intensity increases from interval A → E. Although the radii inferred from the O feature and Fe feature differ, they agree at the $3\sigma$ level of confidence, again indicating an evolution in $R_{\text{in}}$. Figure 5 shows the O and Fe line profile in velocity space during interval C. These profiles exhibit comparable broadening, suggesting that they do indeed arise from a similar location in the accretion disk. We only show a single interval for clarity, but the other intervals show analogous behavior. The inclination is between $i = 70^\circ–80^\circ$, which is consistent with the range of inclination values obtained in Madej & Jonker (2011) and Madej et al. (2014). We also tie $R_{\text{in}}$ between each DISKLINE component (Model 3 in Table 1). The fit is $3.9\sigma$ worse than when we allow the O and Fe components to vary in radius (Model 2), but the inclination and overall trend of the inner accretion disk moving toward the compact object remains. Rather than asserting that these lines arise from the same area of the disk, we conservatively proceed with the individual measurements, though this does not exclude the possibility that the lines originate from the same region.

We find a multiplicative constant between NICER and Swift of $C = 0.93 \pm 0.03$ for each model in Table 1. The constant between NICER and INTEGRAL is smaller ($C \sim 0.30–0.50$), but this is likely due to detector-based flux offsets. The low value could also partially be due to spectral evolution between the two fitting bands. However, alternative models, such as a broken or cutoff power-law, do not offer statistical improvements and yield some fairly extreme parameters. There are currently no direct measurements of the cross-calibration constant between INTEGRAL and NICER. Future investigations using simultaneous observations of sources with NICER, NuSTAR, and INTEGRAL will likely be able to provide a better determination of the cross-calibration constant as there would be joint energy coverage between all three missions.

<figure><img src="image_5.png" /><figcaption>**Figure 4.** Spectral modeling of NICER (orange) and INTEGRAL (teal) data in interval E with the residuals divided by the errors shown in the lower panel. The dashed line indicates the multi-temperature blackbody component. The dotted-dashed line corresponds to the single-temperature blackbody. The dot-dot-dashed line indicates the power-law component. The solid lines indicate the DISKLINE components for O and Fe emission. This is representative of the broadband modeling applied to each interval, but only one interval is shown for simplicity. Data were rebinned for clarity.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Figure 5.** The Fe and O line profiles in NICER from interval C plotted in velocity space, referred to the line centroid energies reported in Table 1. The lines show similar broadening profiles indicating a common origin within the accretion disk.</figcaption></figure>

Since the oxygen line profile shows a clear evolution with time, we examine the O line flux versus continuum component flux in intervals A–E in order to discern what the line is responding to within the system, similar to what was done in Lin et al. (2010) for 4U 1705−44. Figure 6 shows the bolometric oxygen line flux versus the unabsorbed 0.57–100 keV flux for the multi-temperature disk blackbody, single-temperature blackbody, and power-law component. This energy band is chosen to encompass all ionizing radiation above the ionization threshold energy for O VII that could possible contribute to the flux of the line. The purple dashed line indicates the ratio of 10% of the line flux to the respective continuum component. The disk blackbody component shows a positive trend with the oxygen line flux as the source increases in intensity.
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Oxygen or iron line flux vs. each continuum component for the five different time intervals shown in Figure 1. The dashed purple line represents the trend expected if the O or Fe line flux is 10% of the continuum component. The bottom right panel shows the O line flux vs. the Fe line flux throughout the outburst; the dotted-dashed line indicates a slope of one as reference. Both the O and Fe line fluxes increase with time, with the exception of Fe in interval E.</figcaption></figure>

A simple Spearman rank test, a measure of the statistical dependence between the line flux and continuum component, returns a correlation coefficient of 0.9 (with an error of 0.3 and *p*-value of 0.037). While the *p*-value is not very significant, as it is limited by the available points, the other components have a Spearman rank coefficient consistent with zero and high *p*-values. This implies that the change in the O line shape is correlated with changes in the accretion disk component,

perhaps due to the combination of (1) a change in the ionization state of the material as the accretion disk temperature increases (see Table 1), (2) more line emission as the disk moves inwards and a larger surface area is able to be irradiated by the external ionizing flux, and/or (3) a change in accretion geometry.

We also examine which continuum components the Fe line is responding to throughout the outburst. We use the unabsorbed flux in the 6.4–100 keV band for each continuum
component. The thermal disk component does not contribute to the flux in the Fe band, therefore we only examine the correlation with the single-temperature blackbody and power law in Figure 6. The bolometric flux estimates for the broadened Fe K component are more uncertain than for the O line because fewer counts are available. There appears to be a smaller change in the Fe flux over the different intervals, but the Fe line is more correlated with the blackbody component than the changes in the power law. Moreover, Figure 6 shows a comparison of the change in the O line flux and Fe line flux throughout the outburst. There is an overall increase in bolometric line flux for each as the source intensity increases, with the exception of the Fe component in interval E. In any case, the O component serves as a better diagnostic tool for NICER in UCXBs since there is a smaller signal-to-noise ratio in the Fe K band as the effective area declines rapidly, as well as the attenuation of Fe K emission by the overabundance of O in the accretion disk material (Koliopanos et al. 2013).

The physical conditions of the line-emitting plasma can be characterized by the ionization parameter $\xi = L/n r^2$, where L is the ionizing luminosity from 1 to 1000 Ryd, n is the number density of the plasma, and r is the distance of the plasma from the ionizing radiation source taken to be $R_{in}$. We can estimate a plausible change in $\xi$ at the beginning and end of the 10 day NICER monitoring period studying the evolution of the O emission features. Figure 13 of Kallman & Bautista (2001) shows the ionization balance of an optically thin photoionized plasma at high density, as a function of $\xi$. Although we expect the accretion disk to be optically thick so we cannot use the exact values for $\xi$, this still allows us to qualitatively estimate how much $\xi$ should change in order to produce the evolution in spectral features that we observe. Given (1) the presence of an O VII edge while O VIII is the dominant ion in interval A and (2) that O VII is no longer present by interval E and the Fe K line is formally consistent with a blending of ionized Fe XXV with neutral species (Kallman & Bautista 2001; Kallman et al. 2004), we determine that $\Delta(\log \xi) \sim 0.3-1.2$ (cgs).

From the limit on $\Delta(\log \xi)$, we can look into the change in the density of the material in the line-emitting region using the following equation:

$$ \Delta(\log n) \approx \log \left( \frac{F_E \xi_A r_A^2}{F_A \xi_E r_E^2} \right). \qquad (1) $$

The ratio of the unabsorbed ionizing flux in interval E to interval A is $F_E/F_A = 1.22$. Using the lower limit on inner disk radius in interval A (188 $R_g$) and upper limit in interval E (8.7 $R_g$), as well as the estimate for the change in ionization, we determine that the density increased by 1.6–2.5 orders of magnitude over the course of the outburst. This change in density is large but could still fall within the range of density assumed within reflection models ($n \sim 10^{15}-10^{19}$: García et al. 2016).

Moreover, we can use the equation for density of a standard thin disk in Frank et al. (2002) as a cross check for the range of values obtained from Equation (1): $\rho = 3.1 \times 10^{-8} \alpha^{-7/10} M_{16}^{11/20} m_1^{5/8} R_{10}^{-15/8} f^{11/5}$ g cm$^{-3}$, where $\alpha$ is the viscosity parameter, $M_{16}$ is the mass accretion rate in units of $10^{16}$ g s$^{-1}$, $m_1$ is the mass of the compact object in solar masses, $R_{10}$ is the radius in units of $10^{10}$ cm, and $f = [1 - R_{NS}/R_{in}]^{1/4}$. $\Delta(\log \rho)$ and $\Delta(\log n)$ should change by the same amount because they are related by a constant. Assuming

that the viscosity parameter ($\alpha$) is constant, the change in density becomes:

$$ \Delta(\log \rho) \approx \log \left[ \left( \frac{F_E}{F_A} \right)^{11/20} \left( \frac{f_E}{f_A} \right)^{11/5} \left( \frac{r_A^2}{r_E^2} \right)^{-15/8} \right]. \qquad (2) $$

For an NS with a 10 km radius and the same values for flux and inner disk radius used in Equation (1), the density in the line-emitting region increases by ~2.5 orders of magnitude as the disk moves closer to the NS and the system becomes more luminous. This is independent of the change in ionization and in agreement with our previous estimate of 1.6–2.5. It is important to note that Equation (1) depends on the estimated change in $\xi$ and strongly on the change in inner disk radius ($\propto r^2$). Therefore, these estimates should be regarded with a degree of caution. For example, if we instead use the values reported in Table 1 for the Fe line component of 28 $R_g$ in interval A and 10.2 $R_g$ in interval E, we obtain a smaller change of $\Delta(\log n) \sim -0.23-0.66$.

Regardless, without fully self-consistent reflection modeling, we are unable to place further constraints on the change in ionization or disk density. We tried fitting the data with an updated version of XILLVERCO that has more grid points, but have identified limitations within the model, such as a steep dependence on small changes in $kT$, that resulted in unsatisfactory fits ($\chi^2_{\nu} > 80$). Therefore, we chose to not include them here. We are currently revising and updating the existing XILLVERCO models and their application will be performed in a future investigation utilizing a sample of UCXBs.

## 3.2. Timing

From a 0.3–9.5 keV light curve with 1/8192 s time resolution, we construct averaged power spectra for each interval according to usual methods (see, e.g., van der Klis 1989; Bult et al. 2018), normalizing the power to units of fractional rms with respect to the total source count rate. The resulting power spectra are fit with a sum of Lorentzian profiles (Belloni et al. 2002) and labeled according to the atoll naming conventions (van der Klis 2006). Note that the frequencies are also consistent with quasi-periodic oscillations (QPOs) seen in stellar mass BHs and a different naming scheme could have been applied (Klein-Wolt & van der Klis 2008). However, since 4U 1543–624 has exhibited a tentative Type-I X-ray burst, we opt for interpreting these features in the context of an NS.

We find that the power spectra of 4U 1543–624 are well described by the sum of four or five Lorentzian profiles. The best-fit parameters are shown in Table 2 and the individual power spectra and their best-fit model are shown in Figure 7. The narrow QPO seen at high frequencies in intervals A and B (>100 Hz, labeled “hHz” for hectohertz) is marginally significant at a detection level of 3.1 and 3.0, respectively. A visual inspection of interval D shows a similar narrow feature at a frequency of 250 Hz, however, with a signal-to-noise ratio of ~1.3 this feature is not significantly detected. We find that as the source evolves to higher count rate (interval A through E), the power spectrum components gradually move to higher frequencies, while decreasing in fractional rms amplitude. This trend is consistent with the expected behavior of a hard/ intermediate state accreting X-ray binary.
<table><caption>**Table 2**
Power Spectrum Fit Parameters</caption><thead><tr><th></th><th>Frequency<br>(Hz)</th><th>Quality<br>Factor</th><th>Fractional<br>Amplitude<br>(% rms)</th><th>χ<sup>2</sup>/dof</th></tr></thead><tbody><tr><td colspan="5"><strong>Interval A</strong></td></tr><tr><td>break2</td><td>0.16(0.03)</td><td>0.28(0.09)</td><td>6.2(0.8)</td><td rowspan="5">253/204</td></tr><tr><td>break</td><td>0.658(0.019)</td><td>0.47(0.08)</td><td>12.4(0.6)</td></tr><tr><td>hump</td><td>3.89(0.07)</td><td>0.59(0.05)</td><td>16.3(0.5)</td></tr><tr><td>low</td><td>32.6(6.5)</td><td>0.0 (fixed)</td><td>12.1(0.6)</td></tr><tr><td>hHz</td><td>143.8(3.9)</td><td>6.8(3.9)</td><td>4.9(1.0)</td></tr><tr><td colspan="5"><strong>Interval B</strong></td></tr><tr><td>break</td><td>0.87(0.05)</td><td>0.08(0.03)</td><td>15.5(0.4)</td><td rowspan="5">224/208</td></tr><tr><td>hump</td><td>5.18(1.4)</td><td>0.58(0.08)</td><td>14.6(0.7)</td></tr><tr><td>low</td><td>42(12)</td><td>0.0 (fixed)</td><td>11.2(0.7)</td></tr><tr><td>hHz</td><td>245(12)</td><td>3.2(1.3)</td><td>7.1(0.9)</td></tr><tr><td colspan="5"><strong>Interval C</strong></td></tr><tr><td>break</td><td>0.94(0.04)</td><td>0.10(0.02)</td><td>15.5(0.3)</td><td rowspan="5">264/185</td></tr><tr><td>hump</td><td>5.96(1.5)</td><td>0.56(0.07)</td><td>14.5(0.6)</td></tr><tr><td>low</td><td>58.0 (4.0)</td><td>0.0 (fixed)</td><td>10.8(2.7)</td></tr><tr><td>hHz</td><td>157.0 (23)</td><td>1.0(1.0)</td><td>8.0(4.0)</td></tr><tr><td colspan="5"><strong>Interval D</strong></td></tr><tr><td>break2</td><td>0.32(0.06)</td><td>0.33(0.10)</td><td>5.9(1.1)</td><td rowspan="5">233/207</td></tr><tr><td>break</td><td>1.19(0.04)</td><td>0.44(0.10)</td><td>12.4(1.0)</td></tr><tr><td>hump</td><td>7.3(0.3)</td><td>0.29(0.06)</td><td>16.2(0.5)</td></tr><tr><td>hHz</td><td>367(12)</td><td>0.0 (fixed)</td><td>12.8(1.3)</td></tr><tr><td colspan="5"><strong>Interval E</strong></td></tr><tr><td>break2</td><td>0.47(0.08)</td><td>0.27(0.08)</td><td>6.5(0.9)</td><td rowspan="5">227/206</td></tr><tr><td>break</td><td>1.74(0.06)</td><td>0.54(0.11)</td><td>11.3(0.9)</td></tr><tr><td>hump</td><td>9.7(0.4)</td><td>0.33(0.08)</td><td>14.7(0.6)</td></tr><tr><td>hHz</td><td>223.0 (42)</td><td>0.4(0.3)</td><td>11.2(1.2)</td></tr></tbody></table>

**Note.** Best fit parameter values for the multi-Lorentzian models describing the power spectra of 4U 1543–624. The values in parentheses indicate $1\sigma$ uncertainties.

Lastly, we searched for coherent pulsations in the data during the 10 day period of NICER monitoring. This would provide further evidence that the compact object is an NS and classify the source as an X-ray pulsar. Both the individual and the total power spectra were searched, but no apparent coherent oscillations were found. It is possible that the Doppler modulation associated with the binary orbital motion sup- presses the pulse signal to below our detection sensitivity (see, e.g., Strohmayer et al. 2018b), however, a full acceleration search is beyond the scope of this work.

## 4. Discussion

We present evolution in the X-ray spectrum of 4U 1543–624 during a period of enhanced accretion through NICER monitoring over a ~10 day period. We divided the data into five intervals that were labeled A–E. The soft color decreased with time while the hard color remained constant, suggesting an evolving low-energy thermal component within the system. Additionally, we obtained observations with Swift that were contemporaneous with interval A, as well as observations with INTEGRAL and ATCA in interval E. The INTEGRAL observation played an important role in determining the

appropriate continuum model for 4U 1543–624, since we were able to rule out models that were unable to describe the NICER and INTEGRAL passbands simultaneously.

The NICER spectra exhibit emission lines from O at ~0.64 keV and Fe K at ~6.4 keV throughout the outburst. The O line showed a clear evolution in the shape of its profile as the source intensity increased. The Fe line profile does not appear to change as dramatically, which is likely due to a combination of (1) there being three times less NICER collecting area in the Fe band in comparison to the O region and (2) the line emission being dampened from the over- abundance in C/O from the companion (Koliopanos et al. 2013). Both show a similar line profile when plotted in velocity space, indicating that they arise from a similar location within the disk. With the advantage of the large collecting area of NICER below 1.5 keV, we are able to track changes in the accretion disk using the O line component. The innermost accretion disk moves from hundreds of $R_g$ initially to <8.7 $R_g$ at peak intensity.

The evolution in the power spectrum is consistent with an X-ray binary in the hard/intermediate state. As expected, the QPO frequencies move as the source intensity changes (Hansinger & van der Klis 1989), aside from the hHz feature, which remains consistent within a factor of 2 (van Straaten et al. 2003). The upper kHz QPO is often associated with the inner accretion flow, and although we are unable to detect any kHz QPOs, it is known that the frequency of kHz QPOs moves in conjunction with the lower QPO features (van der Klis 2004). Since the QPO features are moving toward higher frequency with time and the inner disk radius is decreasing, it stands to reason that the process controlling the QPO frequency is moving with the inner radius. In this regard, the mechanisms producing the QPOs are very likely dependent upon the inner disk radius. However, without measuring the upper kHz QPO frequency we are unable to test if the kHz QPO and reflection features are produced in the same region of the accretion disk. Figure 8 shows the inner disk radius inferred from DISKLINE fitting versus the change in the O line flux, break frequency, and hump frequency with time.

The Chandra observation of 4U 1543–624 analyzed in Madej & Jonker (2011) and Madej et al. (2014) occurred at a similar flux as the NICER observations at peak intensity (i.e., $F_{0.5-10.0 keV} \sim 1 \times 10^{-9} erg s^{-1} cm^{-2}$). The inner disk radius inferred from the O line in the final two intervals (D: $R_{in} = 7.9_{-1.3}^{+2.5} R_g$, E: $R_{in} = 7.3_{-1.0}^{+1.4} R_g$) agree well with the measurements from Madej et al. (2014) using Chandra/LETGS and the initial grid of XILLVERCO ($R_{in} < 7.4 R_g$). The equivalent width of the O lines in interval D (~27 eV) and E (~33 eV) also agree with the range reported in Madej & Jonker (2011) of 27–35 eV. Moreover, the inclination of $i \sim 70^{\circ}-80^{\circ}$ agrees with the upper limits reported in both studies. Our results are subject to change as our understanding of the NICER instrument calibration improves, but we expect that this will only have a minor impact and the agreement with previous studies is reassuring.

Although often assumed to be an NS the identity of the central object of 4U 1543–624 is ultimately unknown. Constraining the mass function of the binary could conclus- ively confirm the presence of a BH or massive NS, but would require repeated spectroscopy of the faint counterpart— currently impossible in a very short period system. The Type-I X-ray burst of Serino et al. (2018) would confirm that
<figure><img src="image_8.png" /><figcaption>**Figure 7.** Power density spectra of 4U 1543−624 in intervals A-E with their best-fit model components, as well as a comparison of each interval. The power spectrum tends toward higher frequency and decreases in fractional rms amplitude as the source increases in intensity. This is in agreement with the behavior expected for an X-ray binary in the hard/intermediate state.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 8.** The change in inner disk radius measured from the O line vs. the O line flux, break frequency, and hump frequency throughout the outburst. Horizontal dashed purple line indicates 6 $R_g$ or the innermost stable circular orbit for $a = 0$. The inner accretion disk moves toward the compact object and the QPOs move to higher frequency as the source intensity increases (see Figure 1).</figcaption></figure>

the compact accretor is an NS, but the resolution of the MAXI instrument (1.°5) precludes a definite association with 4U 1543−624 as there is another known X-ray source within the MAXI beam. X-ray pulsations would also establish that 4U 1543−624 contained an NS, but none were found in our inspection of the NICER data. Finally, inspection of our simultaneous ATCA radio and NICER observations in interval E could constrain the nature and behavior of 4U 1543−624 in the context of the known X-ray and radio properties of other LMXB systems.

The empirically derived $L_r-L_x$ fundamental plane displays the relationship between radio and X-ray luminosity for BHs and different classes of NS LMXBs, and aids in the understanding of the inflow/outflow efficiency of these systems (Gallo et al. 2018). In order to place 4U 1543−624 on the $L_r-L_x$ plane, we need to convert the radio flux density and X-ray flux into luminosities, which depends on the distance to the source. Wang & Chakrabarty (2004) originally estimated a distance of ≈7 kpc by calculating a systemic mass accretion rate using the orbital period, and then scaling the corresponding X-ray luminosity to that previously observed (Juett & Chakrabarty 2003). The Type-I X-ray burst reported in Serino et al. (2018) places the source at a distance of 9.2 ± 2.3 kpc, assuming the empirical Eddington limit of $3.8 \times 10^{38}$ erg s⁻¹ (Kuulkers et al. 2003). Bailer-Jones et al. (2018) infer a limit on the distance to the optical companion of 4U 1543−624 of $d = 3.3_{-1.9}^{+3.3}$ kpc using Gaia DR2 parallax measurements with the assumption of an exponentially decreasing space density prior (Bailer-Jones 2015; Astraatmadja & Bailer-Jones 2016). Given the uncertainty surrounding the source location, we choose to encompass the entire range of distance estimates, 1.4 kpc < *d* < 11.5 kpc, in our luminosity calculations.

Assuming a flat radio spectrum, we convert our flux density upper limit to its corresponding 5.0 GHz luminosity with $L_r = 4\pi d^2\nu S_\nu$, where $S_\nu = 27 \ \mu\text{Jy}$, $\nu = 5.0$ GHz, and *d* is the distance to 4U 1543−624. We then convert the NICER unabsorbed 1.0–10 keV flux of ~9.34 × 10⁻¹⁰ erg s⁻¹ cm⁻² into an X-ray luminosity. By comparing the X-ray-radio luminosity to other known accreting LMXBs, we can determine if the properties of 4U 1543−624 closely resemble those of another class of LMXBs. Figure 9 shows the region (in red) that 4U 1543−624 traces out on the $L_r-L_x$ plane based upon the current estimates of the distance to the source. It is immediately clear that 4U 1543−624 has a radio luminosity

<figure><img src="image_10.png" /><figcaption>**Figure 9.** Radio—X-ray luminosity plane of accreting compact objects adapted from Bahramian et al. (2018). The dashed gray line indicates the best-fit trends for BHs ($L_r \propto L_x^{0.6}$; Gallo et al. 2014). The purple dotted-dashed lines indicate the proposed correlations for NSs (Migliari & Fender 2006) for atoll and Z sources ($L_r \propto L_x^{0.6}$) and for NSs in the hard state ($L_r \propto L_x^{1.4}$). The black points are quiescent/hard state BHs from the available literature (Gallo et al. 2003, 2014; Corbel et al. 2004, 2013; Brocksopp et al. 2005, 2013; Calvelo et al. 2010; Coriat et al. 2011b; Rodriguez et al. 2011; Paizis et al. 2011; Jonker et al. 2012; Russell et al. 2015; Plotkin et al. 2017; Ribó et al. 2017; Dinçer et al. 2018). The light and dark blue squares are NSs in hard states while orange pentagons are NSs in soft states (Rutledge et al. 1998; Migliari & Fender 2006; Tudose et al. 2009; Coriat et al. 2011a; Migliari et al. 2011; Tetarenko et al. 2016; Gusinskaia et al. 2017; Tudor et al. 2017); the dark blue is used to differentiate the UCXB 4U 0614+091 from the other hard state NSs. The magenta stars are accreting millisecond X-ray pulsars (AMXPs: Migliari et al. 2005, 2011; Tudor et al. 2016, 2017; Tetarenko et al. 2017; Strohmayer et al. 2018a; van den Eijnden et al. 2018a). The green triangles are transitional millisecond pulsars (tMSPs) when they have switched from rotation powered to accretion powered mode (Hill et al. 2011; Papitto et al. 2013; Deller et al. 2015; Bogdanov et al. 2018). The purple upside-down triangles are Z sources (Migliari & Fender 2006). The shaded red region indicates where 4U 1543−624 resides based upon the X-ray flux in interval E, the 27 µJy upper limit on the radio flux density, and the uncertainty on distance. The source lies two orders of magnitude below the trend for black holes.</figcaption></figure>

that is (a minimum of) two orders of magnitude below what is typically observed for BHs. Although the $L_r-L_x$ plane should be used with caution in making definitive claims about an LMXB’s nature (e.g., IGR J17591−2342: Ferrigno et al. 2018;
Russell et al. 2018), we consider its position to be highly suggestive that the central object is some class of NS.

The utility of the $L_r-L_x$ plane for separating NSs by class is somewhat more limited. We consider a Z type origin unlikely as 4U 1543–624 is at least a magnitude fainter in both radio and X-ray luminosity than the faintest Z source, and displays none of the fast X-ray variability typical in such systems (Migliari & Fender 2006). The X-ray properties of 4U 1543−624 alone point away from a soft state system, as its spectrum has been observed to be persistently hard. The behavior of accreting tMSPs is more difficult to characterize as there are few confirmed tMSP systems. Presently, tMSPs are among the most radio-bright NS systems and, if they obey their current observed $L_r-L_x$ trend (Deller et al. 2015), are expected to be at least an order of magnitude brighter in the radio than our upper limit for 4U 1543−624. Furthermore, they exhibit many other unique phenomena, e.g., optical pulsations, radio pulsations, rapid X-ray and radio flaring, gamma-ray emission—none of which have been observed in 4U 1543−624 (Hill et al. 2011; Papitto et al. 2013; Bogdanov & Halpern 2015; Bogdanov et al. 2015; Zampieri et al. 2019).

To further aid our interpretation, we consider 4U 0614+091, another UCXB that exhibits similar X-ray spectral properties to 4U 1543−624. Both show the presence of a prominent O VIII feature and a broad Fe K line (Madej et al. 2014; Ludlam et al. 2019). Unlike 4U 1543−624, 4U 0614+091 has been detected in the radio band. The source is located on the $L_r-L_x$ plane (dark blue squares in Figure 9) near $L_{r,5}$ GHz ~ $10^{28}$ erg s⁻¹ and $L_{x,1-10}$ keV ~ $10^{36}$ erg s⁻¹ (Migliari et al. 2005). For the same X-ray luminosity, 4U 1543−624 would be more than five times less radio-luminous than 4U 0614+091. The upper limit on the magnetic field in 4U 0614+091 is $B \le 14.5 \times 10^8$ G at the poles (Ludlam et al. 2019). Using the same procedure, we can estimate a conservative upper limit on the B-field in 4U 1543−624 for comparison. The unabsorbed flux in the 0.5–50 keV band is $2.01 \times 10^{-9}$ erg s⁻¹ cm⁻² in interval E and the upper limit on $R_{in}$ is 8.7 $R_g$. This provides an upper limit of $B \le 3.2 \times 10^8$ G at $d = 11.5$ kpc or $B \le 0.4 \times 10^8$ G at $d = 1.42$ kpc.

The difference in the magnetic field strength of 4U 1543−624 and 4U 0614+091 could contribute to the discrepancy in radio luminosity between these two UCXBs, but how magnetic field strength affects radio emission in NS binaries is poorly understood. Recent studies of both AMXPs and non-pulsating NS LMXBs do not find obvious trends in their radio emission with respect to magnetic field strength or other systemic properties (Tetarenko et al. 2016; Tudor et al. 2017), although AMXPs overall have slightly higher radio luminosities than typical NS systems. It is suggested that an interplay between NS spin and magnetic field strength could have an important role in the efficiency and method of jet production (van den Eijnden et al. 2018b), but further investigation on this point is required.

Although the absence of X-ray pulsations in 4U 1543−624 favors a typical hard state NS scenario, we cannot rule out that the source contains an AMXP whose pulsations are shielded by the accretion geometry (Lamb et al. 2009) or suppressed below our detection threshold by the short-period orbital motion (Strohmayer et al. 2018b). Regardless, the radio and X-ray relationship of 4U 1543−624 strongly favors an NS over a BH and adds to the evidence that NSs do not follow a single track in the $L_r-L_x$ plane, limiting its use in distinguishing between

different classes of NSs based on radio and X-ray observations only.

This work was supported by NASA through the NICER mission and the Astrophysics Explorers Program, and made use of data and software provided by the High Energy Astrophysics Science Archive Research Center (HEASARC). The Australia Telescope Compact Array is part of the Australia Telescope National Facility which is funded by the Australian Government for operation as a National Facility managed by CSIRO. R.M.L. acknowledges partial funding from a NASA Earth and Space Science Fellowship. E.M.C. gratefully acknowledges support through CAREER award number AST-1351222. M.F. acknowledges ASI financial/programmatic support via ASI-INAF agreement n. 2013-025.R1 and ASI-INAF N. 2017-14-H.O. J.A.G. acknowledges support from Chandra Theory grant TM8-19003X and from the Alexander von Humboldt Foundation. J.C.A. M.-J. is the recipient of an Australian Research Council Future Fellowship (FT 140101082), funded by the Australian Government.

## ORCID iDs

R. M. Ludlam https://orcid.org/0000-0002-8961-939X

A. Zoghbi https://orcid.org/0000-0002-0572-9613

T. E. Strohmayer https://orcid.org/0000-0001-7681-5845

G. K. Jaisawal https://orcid.org/0000-0002-6789-2723

S. Guillot https://orcid.org/0000-0002-6449-106X

M. Fiocchi https://orcid.org/0000-0001-5697-6019

A. C. Fabian https://orcid.org/0000-0002-9378-4072

E. M. Cackett https://orcid.org/0000-0002-8294-9281

A. Bahramian https://orcid.org/0000-0003-2506-6041

## References

Arnaud, K. A. 1996, in ASP Conf. Ser. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby & J. Barnes (San Francisco, CA: ASP), 17

Astraatmadja, T. L., & Bailer-Jones, C. A. L. 2016, ApJ, 832, 137

Bahramian, A., Miller-Jones, J. C. A., Strayder, J., et al. 2018, Radio/X-ray Correlation Database for X-ray binaries, Zenodo, doi:10.5281/zenodo.1252036

Bailer-Jones, C. A. L. 2015, PASP, 127, 994

Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., et al. 2018, AJ, 156, 58

Ballantyne, D. R., Ross, R. R., & Fabian, A. C. 2002, MNRAS, 336, 867

Belloni, T., Psaltis, D., & van der Klis, M. 2002, ApJ, 572, 392

Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883

Blandford, R. D., & Znajek, R. L. 1977, MNRAS, 179, 433

Bogdanov, S., Archibald, A. M., Bassa, C., et al. 2015, ApJ, 806, 148

Bogdanov, S., Deller, A. T., Miller-Jones, J. C. A., et al. 2018, ApJ, 856, 54

Bogdanov, S., & Halpern, J. P. 2015, ApJL, 803, L27

Brocksopp, C., Corbel, S., Fender, R. P., et al. 2005, MNRAS, 356, 125

Brocksopp, C., Corbel, S., Tzioumis, A., et al. 2013, MNRAS, 432, 931

Bult, P., Altamirano, D., Arzoumanian, Z., et al. 2018, ApJL, 860, L9

Cackett, E. M., Miller, J. M., Reis, R. C., Fabian, A. C., & Barret, D. 2012, ApJ, 755, 27

Calvelo, D., Fender, R. P., Russell, D. M., et al. 2010, MNRAS, 409, 839

Corbel, S., Coriat, M., Brocksopp, C., et al. 2013, MNRAS, 428, 2500

Corbel, S., Fender, R. P., Tomsick, J. A., Tzioumis, A. K., & Tingay, S. 2004, ApJ, 617, 1272

Coriat, M., Corbel, S., Prat, L., et al. 2011a, MNRAS, 414, 677

Coriat, M., Tzioumis, T., Corbel, S., et al. 2011b, A Tel, 3119, 1

Courvoisier, T. J.-L., Walter, R., Beckmann, V., et al. 2003, A&A, 411, L53

Deller, A. T., Moldon, J., Miller-Jones, J. C. A., et al. 2015, ApJ, 809, 13

Dinçer, T., Bailyn, C. D., Miller-Jones, J. C. A., Buxton, M., & MacDonald, R. K. D. 2018, ApJ, 852, 4

Fabian, A. C., Iwasawa, K., Reynolds, C. S., & Young, A. J. 2000, PASP, 112, 1145

Fabian, A. C., Rees, M. J., Stella, L., & White, N. E. 1989, MNRAS, 238, 729
Federici, M., Martino, B. L., & Ubertini, P. 2010, ExA, 34, 105

Fender, R. P., Gallo, E., & Jonker, P. G. 2003, MNRAS, 343, L99

Ferrigno, C., Bozzo, E., Sanna, A., et al. 2018, ATel, 11957, 1

Frank, J., King, A., & Raine, D. 2002, Accretion Power in Astrophysics (3rd ed.; Cambridge: Cambridge Univ. Press)

Gallo, E., Degenaar, N., & van den Eijnden, J. 2018, MNRAS, 478, L132

Gallo, E., Fender, R. P., & Pooley, G. G. 2003, MNRAS, 344, 60

Gallo, E., Miller-Jones, J. C. A., Russell, D. M., et al. 2014, MNRAS, 445, 290

García, J. A., Fabian, A. C., Kallman, T. R., et al. 2016, MNRAS, 462, 751

Gendreau, K. C., Arzoumanian, Z., & Okajima, T. 2012, Proc. SPIE, 8443, 13

Giacconi, R., Murray, S., Gursky, H., et al. 1972, ApJ, 178, 281

Gusinskaia, N. V., Deller, A. T., Hessels, J. W. T., et al. 2017, MNRAS, 470, 1871

Hansinger, G., & van der Klis, M. 1989, A&A, 225, 79

Hill, A. B., Szostek, A., Corbel, S., et al. 2011, MNRAS, 415, 235

Jahoda, K., Markwardt, C. B., Radeva, Y., et al. 2006, ApJS, 163, 401

Jonker, P. G., Miller-Jones, J. C. A., Homan, J., et al. 2012, MNRAS, 423, 3308

Juett, A. M., & Chakrabarty, D. 2003, ApJ, 599, 498

Juett, A. M., Psaltis, D., & Chakrabarty, D. 2001, ApJL, 560, L59

Kallman, T., & Bautista, M. 2001, ApJS, 133, 221

Kallman, T. R., Palmeri, P., Bautista, M. A., Mendoza, C., & Krolik, J. H. 2004, ApJS, 155, 165

Klein-Wolt, M., & van der Klis, M. 2008, ApJ, 675, 1407

Koliopanos, F., Golfanov, M., & Bildsten, L. 2013, MNRAS, 432, 1264

Kuulkers, E., den Hartog, P. R., in't Zand, J. J. M., et al. 2003, A&A, 399, 663

Kuulkers, E., in't Zand, J. J. M., Atteia, J.-L., et al. 2010, A&A, 514, A65

Lamb, F. K., Boutloukos, S., Van Wassenhove, S., et al. 2009, ApJ, 706, 427

Lin, D., Remillard, R. A., & Homan, J. 2007, ApJ, 667, 1073

Lin, D., Remillard, R. A., & Homan, J. 2010, ApJ, 719, 1350

Ludlam, R. M., Miller, J. M., Arzoumanian, Z., et al. 2018, ApJL, 858, L5

Ludlam, R. M., Miller, J. M., Barret, D., et al. 2019, ApJ, 873, 99

Ludlam, R. M., Miller, J. M., Miller-Jones, J., & Reynolds, M. 2017, ATel, 10690, 1

Madej, O. K., García, J., Jonker, P. G., et al. 2014, MNRAS, 442, 1157

Madej, O. K., & Jonker, P. G. 2011, MNRAS, 412, L11

Matsuoka, M., Kawasaki, K., Ueno, S., et al. 2009, PASJ, 61, 999

McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in ASP Conf. Ser. 376, Astronomical Data Analysis Software and Systems XVI, ed. R. A. Shaw, F. Hill, & D. J. Bell (San Francisco, CA: ASP), 127

Migliari, S., & Fender, R. 2006, MNRAS, 366, 79

Migliari, S., Fender, R. P., & van der Klis, M. 2005, MNRAS, 363, 112

Migliari, S., Miller-Jones, J. C. A., & Russell, D. M. 2011, MNRAS, 415, 2407

Nelemans, G., Jonker, P. G., Marsh, T. R., & van der Klis, M. 2003, MNRAS, 348, L7

Ng, C., Díaz Trigo, M., Cadolle Bel, M., & Migliari, S. 2010, A&A, 522, A96

Paizis, A., Nowak, M. A., Wilms, J., et al. 2011, ApJ, 738, 183

Papitto, A., Ferrigno, C., Bozzo, E., et al. 2013, Natur, 501, 517

Plotkin, R. M., Bright, J., Miller-Jones, J. C. A., et al. 2017, ApJ, 848, 92

Ribó, M., Munar-Adrover, P., Paredes, J. M., et al. 2017, ApJ, 835, 33

Rodriguez, J, Corbel, S., Caballero, I., et al. 2011, A&A, 533, L4

Russell, T., Degenaar, N., Wijnands, R., & van den Eijnden, J. 2018, ATel, 11954, 1

Russell, T. D., Miller-Jones, J. C. A., Curran, P. A., et al. 2015, MNRAS, 450, 1745

Rutledge, R., Moore, C., Fox, D., Lewin, W., & van Paradijs, J. 1998, ATel, 8, 1

Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in ASP Conf. Ser. 77, Astronomical Data Analysis Software and Systems IV, ed. R. A. Shaw, H. E. Payne, & J. J. E. Hayes (San Francisco, CA: ASP), 433

Schultz, J. 2003, A&A, 397, 249

Serino, M., Hiaka, K., Negoro, H., et al. 2018, ATel, 11302, 1

Strohmayer, T. E., Arzoumanian, Z., Bogdanov, S., et al. 2018a, ApJL, 858, L13

Strohmayer, T. E., Ray, P. S., Gendreau, K. C., et al. 2018b, ATel, 11520, 1

Tetarenko, A. J., Bahramian, A., Sivakoff, G. R., et al. 2016, MNRAS, 460, 345

Tetarenko, A. J., Bahramian, A., Sivakoff, G. R., et al. 2017, ATel, 10894, 1

Tetarenko, A. J., Bahramian, A., Wijnands, R., et al. 2016, ApJ, 854, 125

Tudor, V., Bahramian, A., Sivakoff, G., et al. 2016, ATel, 8914, 1

Tudor, V., Miller-Jones, J. C. A., Patruno, A., et al. 2017, MNRAS, 470, 324

Tudose, V., Fender, R. P., Linares, M., Maitra, D., & van der Klis, M. 2009, MNRAS, 400, 2111

Ubertini, P., Lebrun, F., Di Cocco, G., et al. 2003, A&A, 411, L131

van den Eijnden, J., Degenaar, N., Russell, T. D., et al. 2018a, ATel, 11487, 1

van den Eijnden, J., Degenaar, N., Russell, T. D., et al. 2018b, Natur, 562, 233

van der Klis, M. 1989, in Timing Neutron Stars, Vol. 262, ed. H. Ögelman & E. P. J. van den Heuvel (Dordrecht: Kluwer)

van der Klis, M. 2004, arXiv:astro-ph/0410551

van der Klis, M. 2006, Rapid X-ray Variability, Compact Stellar X-ray Sources (Cambridge Univ. Press), 39

van Straaten, S., van der Klis, M., & Méndez, M. 2003, ApJ, 596, 1155

Wang, Z., & Chakrabarty, D. 2004, ApJL, 616, L139

Wang, Z., Tziamtzis, A., Kaplan, D. L., & Chakrabarty, D. 2015, PASA, 32, 35

Zampieri, L., Burtovoi, A., Fiori, M., et al. 2019, MNRAS, 485, L109