## OBSERVATIONS OF THE ULTRA-COMPACT X-RAY BINARY 4U 1543 -624 IN OUTBURST WITH NICER INTEGRAL , , Swift , AND ATCA

R. M. LUDLAM 1 , L. SHISHKOVSKY 2 , P. M. BULT 3 , J. M. MILLER 1 , A. ZOGHBI 1 , T. E. STROHMAYER 3 , M. REYNOLDS 1 , L. NATALUCCI 4 , J. C. A. MILLER-JONES 5 , G. K. JAISAWAL 6 , S. GUILLOT 7,8 , K. C. GENDREAU 3 , J. A. GARC´ IA 9,10 , M. FIOCCHI 4 , A. C. FABIAN 11 , D. CHAKRABARTY 12 , E. M. CACKETT 13 , A. BAHRAMIAN 5 , Z. ARZOUMANIAN 3 , D. ALTAMIRANO 14

Draft version August 2, 2019

## 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 NICER monitoring of the source over a ∼ 10 day period during which target-of-opportunity observations were also conducted with Swift , INTEGRAL , and ATCA. 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 &lt; . 8 7 R g at peak intensity. We do not detect the source in radio, but are able to place a 3 σ upper limit on the flux density at 27 µ 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 (NS) primary. This 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 alone.

Keywords: accretion, accretion disks - stars: neutron - stars: individual (4U 1543 -624) - X-rays: binaries

## 1. INTRODUCTION

Ultra-compact X-ray binaries (UCXBs) are low-mass Xray binaries (LMXBs) with a short orbital period of &lt; 80 minutes, and typically contain a degenerate stellar companion. 4U 1543 -624 is an UCXB with an orbital period of 18 2 . ± 0 1 . minutes determined from optical photometry (Wang &amp; 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).

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

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 1 , although the instrument's spatial resolution of 1.5 ◦ (Matsuoka et al. 2009) prevents a definitive conclusion regarding the event's origin since there are other X-ray sources nearby. If the burst indeed arose from 4U 1543 -624, then this would irrefutably confirm the system contains a 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).

- 2 Center for Data Intensive and Time Domain Astronomy, Department of Physics and Astronomy, Michigan State University, East Lansing, MI, USA
- 4 Istituto Nazionale di Astrofisica, INAF-IAPS, via del Fosso del Cavaliere, 100. 00133 Roma, Italy
- 3 X-ray Astrophysics Laboratory, Astrophysics Science Division, NASA/GSFC, Greenbelt, MD 20771, USA
- 5 International Centre for Radio Astronomy Research - Curtin University, GPO Box U1987, Perth, WA 6845, Australia
- 6 National Space Institute, Technical University of Denmark, Elektrovej 327-328, DK-2800 Lyngby, Denmark
- 7 CNRS, IRAP, 9 avenue du Colonel Roche, BP 44346, F-31028 Toulouse Cedex 4, France
- Universit´ e de Toulouse, CNES, UPS-OMP, F-31028 Toulouse, France 9 Cahill Center for Astronomy and Astrophysics, California Institute of Technology, Pasadena, CA 91125, USA
- 8
- 10 Remeis Observatory &amp; ECAP, Universit¨ at Erlangen-N¨ urnberg, Sternwartstr. 7, D-96049, Bamberg, Germany
- 11 Institute of Astronomy, Madingley Road, Cambridge CB3 0HA, UK
- 12 MIT Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- 13 Department of Physics &amp; Astronomy, Wayne State University, 666 West Hancock Street, Detroit, MI 48201, USA
- 14 Department of Physics &amp; Astronomy, University of Southampton, Highfield, Southampton SO17 1BJ, UK

Further support that 4U 1543 -624 contains a NS can be obtained from observing the source in both radio and Xrays. 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 &amp; Znajek 1977; Blandford &amp; Payne 1982; Tudor et al. 2017). Observing the radio and X-ray luminosity simultaneously can be advantageous in discerning between a NS or BH accretor, since NSs are typically ∼ 30 times fainter in the radio than BHs at the same X-ray luminosity (Migliari &amp; Fender 2006). However, some caution should be maintained in utilizing this diagnostic tool given that certain classes of NSs, namely transitional and accreting millisecond X-ray pul-

350

Intensity (cts/s)

Figure 1. Panel (a): The 0 5 . -6 8 . keV NICER light-curve of 4U 1543 -624 starting on 2017 Aug 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 versus hard color diagram. The single outlier point in interval A is an artifact of the 128 s time binning.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00539v1-figure-c4e49d7f-1eed-4a5f-a438-b73641ada9db.jpg)

sars (tMXPs and AMXPS), can be nearly as radio-luminous as the radio-faint black holes (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.

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 &amp; Jonker 2011).

4U 1543 -624 was first reported to exhibit a broad O VIII Ly α emission feature at ∼ 0 7 . keV by Madej &amp; 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 XMMNewton by Juett &amp; 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 &amp; Jonker 2011). This was not the first time that the Fe line had been detected in this system, since it was present in prior RXTE and EXOSAT data (Schultz 2003), but Madej &amp; 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

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 CO . 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 CO grid on 4U 1543 -624 indicated an inclination of i ∼ 65 ◦ and an inner disk radius &lt; 7 4 . R g (where R g = GM/c 2 ).

In early August 2017, the Swift /BAT detected increased activity from 4U 1543 -624 in the 15 -50 keV energy band (Ludlam et al. 2017). We requested NICER monitoring of the source in order to track any changes in the system as the X-ray intensity increased. Additionally, we secured target-ofopportunity observations with Swift , INTEGRAL , and 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

The Neutron Star Interior Composition Explorer ( 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 2 at 1.5 keV. NICER observed 4U 1543 -624 over a ∼ 10 day period in 2017 August (ObsIDs 10500601011050060113) for a cumulative exposure time of ∼ 86 2 . ks.

Table 1 Spectral Modeling of X-ray Observations

| Model 1               | Parameter                                                             | Interval A                                                                                                                                     | Interval B                                                                                         | Interval C                                                                                                                                     | Interval D                                                                   | Interval E                                                                                                                      |
|-----------------------|-----------------------------------------------------------------------|------------------------------------------------------------------------------------------------------------------------------------------------|----------------------------------------------------------------------------------------------------|------------------------------------------------------------------------------------------------------------------------------------------------|------------------------------------------------------------------------------|---------------------------------------------------------------------------------------------------------------------------------|
| TBABS                 | N H ( 10 21 cm - 2 )                                                  |                                                                                                                                                |                                                                                                    | 2 . 92 +0 . 02 - 0 . 03                                                                                                                        |                                                                              |                                                                                                                                 |
| CONSTANT              | C                                                                     | 0 . 93 ± 0 . 03                                                                                                                                | ...                                                                                                | ...                                                                                                                                            | ...                                                                          | 0 . 33 ± 0 . 03                                                                                                                 |
| DISKBB                | kT ( 10 - 2 keV)                                                      | 4 . 9 +0 . 3 - 0 . 4                                                                                                                           | 8 . 7 +2 . 8 - 0 . 8                                                                               | 9 . 3 +0 . 4 - 0 . 6                                                                                                                           | 12 . 2 +0 . 4 - 0 . 3                                                        | 16 . 1 ± 0 . 3                                                                                                                  |
|                       | norm disk ( 10 5 )                                                    | 150 +222 - 70                                                                                                                                  | 1 . 9 ± 1 . 5                                                                                      | 1 . 2 +0 . 6 - 0 . 3                                                                                                                           | 0 . 36 ± 0 . 06                                                              | 0 . 20 ± 0 . 02                                                                                                                 |
| BBODY                 | kT (keV)                                                              | 0 . 66 ± 0 . 01                                                                                                                                | 0 . 67 +0 . 02 - 0 . 04                                                                            | 0 . 73 ± 0 . 01                                                                                                                                | 0 . 74 +0 . 01 - 0 . 02                                                      | 0 . 65 +0 . 02 - 0 . 01                                                                                                         |
|                       | norm bb ( 10 - 3 )                                                    | 1 . 04 +0 . 04 - 0 . 03                                                                                                                        | 0 . 89 +0 . 05 - 0 . 03                                                                            | 1 . 06 ± 0 . 04                                                                                                                                | 1 . 17 ± 0 . 04                                                              | 1 . 22 ± 0 . 06                                                                                                                 |
| POWERLAW              | Γ                                                                     | 1 . 74 ± 0 . 01                                                                                                                                | 1 . 73 +0 . 01 - 0 . 02                                                                            | 1 . 77 ± 0 . 01                                                                                                                                | 1 . 79 ± 0 . 01                                                              | 1 . 71 ± 0 . 01                                                                                                                 |
|                       | norm pl ( 10 - 1 ) χ 2 (dof)                                          | 1 . 46 ± 0 . 01                                                                                                                                | 1 . 58 +0 . 01 - 0 . 03                                                                            | 1 . 61 ± 0 . 01 8306.01 (4668)                                                                                                                 | 1 . 66 +0 . 01 - 0 . 02                                                      | 1 . 59 ± 0 . 03                                                                                                                 |
| Model 2               | Parameter                                                             | Interval A                                                                                                                                     | Interval B                                                                                         | Interval C                                                                                                                                     | Interval D                                                                   | Interval E                                                                                                                      |
| TBABS CONSTANT DISKBB | N H ( 10 21 cm - 2 ) C kT ( 10 - 2 keV) norm disk ( 10 5 )            | 0 . 93 ± 0 . 03 4 . 6 ± 0 . 2 456 +278 - 150 0 . 69 ± 0 . 01 1 . 07 +0 . 03 - 0 . 05 1 . 77 ± 0 . 1 . 51 +0 . - 0 . 0 . 637 +0 . - 0 . 942 +30 | ... 5 . 7 +0 . 3 - 0 . 2 68 +25 - 24 - 0 . 01 . 637 +0 . - 0 . 46 +65 - 19 0 . 40 +0 . 08 - 0 . 07 | 3 . 06 +0 . 03 - 0 . 02 ... 5 . 6 ± 0 . 2 88 +36 - 19 - 0 . 010 - 2 . 36 +0 . 05 - 0 . 07 44 +39 - 19 74 +6 - 4 0 . 45 +0 . 10 - 0 . 06 +2 . 4 | ... 5 . 7 ± 0 . 2 81 +24 - 17 7 . 9 +2 . 5 - 1 . 3                           | 0 . 44 +0 . 05 - 0 . 04 10 . 6 +0 . 2 - 0 . 6 1 . 4 +0 . 5 - 0 . 1 7 . 3 - 1 . 0 1 . 5 +0 . 2 - 0 . 1 33 +4 . 4 - 2 . 2 +0 . 27 |
| DISKLINE              | ◦ O E Fe (keV) R in ( R g ) norm line ( 10 - 3 ) EW Fe (eV) χ 2 (dof) | 0 . 19 +0 . - 0 . 5 . 0 +1 . 1 - 0 . 5 6 . 42 +0 . -                                                                                           | 10 +2 . 0 - 1 . 75 6 . 7 ± 0 . 2 17 ± 11 0 . 9 +0 . 3 - 0 . 5                                      | 11 - 1 . 5 6 . 41 +0 . 13 - 0 . 01 9 +6 - 2 1 . 2 ± 0 . 3                                                                                      | 1 . 2 ± 0 . 1 27 ± 2 . 3 6 . 41 +0 . 06 - 0 . 01 +7                          |                                                                                                                                 |
| 1                     | Parameter                                                             | 28 0 . 01 28 +26 - 2 +0 . 1                                                                                                                    |                                                                                                    | 210 ± 53                                                                                                                                       | 12 - 5                                                                       | 6 . 48 - 0 . 08 10 . 2 +12 . 5 - 3 . 4                                                                                          |
| BBODY POWERLAW        | kT (keV) norm bb ( 10 - 3 ) Γ norm pl ( 10 - 1 ) E O (keV)            | 01 02 01                                                                                                                                       | 0 . 74 ± 0 . 02 0 . 94 +0 . 06 - 0 . 09 1 . 78 ± 0 . 01                                            | 0 . 80 +0 . 02 - 0 . 01 1 . 24 +0 . 06 - 0 . 08 1 . 83 ± 0 . 01 1 . 66 ± 0 . 01 +0 . 008                                                       | 0 . 83 +0 . 02 - 0 . 01 1 . 36 +0 . 09 - 0 . 07 1 . 86 ± 0 . 01 1 . 74 ± 0 . | 0 . 81 +0 . 04 - 0 . 02 0 . 98 +0 . 08 - 0 . 09 1 . 82 ± 0 . 1 . 86 +0 . 02 - 0 . 01                                            |
|                       |                                                                       |                                                                                                                                                | +0 . 02                                                                                            | 0 . 646                                                                                                                                        | 01 +0 .                                                                      | 01 0 . 659 +0 . 005 - 0 . 008                                                                                                   |
|                       | q R in ( R g )                                                        | 004 009                                                                                                                                        | 1 . 78 0 010 008                                                                                   |                                                                                                                                                | 0 . 632 006 - 0 . 008                                                        | +1 . 4                                                                                                                          |
| DISKLINE 2            |                                                                       | 0 . 6 - 0 . 2 105 +18 - 35                                                                                                                     | 160 +53 - 89                                                                                       |                                                                                                                                                | 1 . 3 ± 0 . 3 228 ± 53                                                       | 0 . 8 ± 0 . 3 123 ± 46                                                                                                          |
|                       | i ( ) norm line ( 10 - 2 ) EW (eV)                                    | - 754 04 02                                                                                                                                    |                                                                                                    |                                                                                                                                                |                                                                              |                                                                                                                                 |
| Model 3               | 21 -                                                                  |                                                                                                                                                | Interval B                                                                                         | 5903.06 (4636) Interval C +0 . 04                                                                                                              | Interval D                                                                   | Interval E                                                                                                                      |
| TBABS                 | N H ( 10 cm 2 ) C kT ( 10 - 2 keV)                                    | Interval A 0 . 93 ± 0 . 03 4 . 7 +0 . 1                                                                                                        | ...                                                                                                | 3 . 05 - 0 . 02                                                                                                                                | ...                                                                          | 0 . 45 +0 . 05                                                                                                                  |
| BBODY                 | norm disk ( 10 5 )                                                    | - 0 . 2 400 +200 - 90 0 . 68 +0 . 02 - 0 . 01                                                                                                  | 74 +74 - 43 0 . 72 ± 0 . 0 . 88 +0 . 05 - 0 . 06                                                   | 5 . 5 - 0 . 3 102 +77 - 22 0 . 79 +0 . 02                                                                                                      | 5 . 6 ± 0 . 97 +47 - 26                                                      | 10 . 1 - 0 . 4 1 . 74 +0 . 5 - 0 . 2 0 . 81 +0 . 03                                                                             |
| DISKBB                | kT (keV) norm bb ( 10 - 3 )                                           | 1 . 00 +0 . 04 - 0 . 03                                                                                                                        | 5 . 7 +0 . 6 - 0 . 4 02                                                                            | ...                                                                                                                                            | 0 . 83 +2 - 0 . 01                                                           | - 0 . 01 1 . 00 +0 . 08                                                                                                         |
| CONSTANT              |                                                                       |                                                                                                                                                | 1 . 77 ± 0 . 01 1 . 64 +0 . 02 - 0 . 01 0 . 63 +0 . 01 - 0 . 02                                    | +0 . 1 - 0 . 01                                                                                                                                | 2                                                                            | - 0 . 04 +0 . 2                                                                                                                 |
| POWERLAW              | Γ norm pl ( 10 - 1 ) E O (keV)                                        | 1 . 76 ± 0 . 01 0 . 637 +0 . 003 - 0 . 017 515 ± 455                                                                                           | 32 +22 - 16                                                                                        | 1 . 17 +0 . 07 - 0 . 04 1 . 82 ± 0 . 01 1 . 66 +0 . 02 - 0 . 01                                                                                | 1 . 40 +0 . 8 - 0 . 06 1 . 87 ± 0 . 01                                       |                                                                                                                                 |
| DISKLINE 1            | q                                                                     | 1 . 51 +0 . 02 - 0 . 01                                                                                                                        |                                                                                                    | 0 . 63 ± 0 . 01 +0 . 05                                                                                                                        | 1 . 74 +0 . 02 - 0 . 01 0 . 63 ± 0 . 01                                      | - 0 . 07 1 . 83 ± 0 . 01 1 . 87 +0 . 02 - 0 . 01                                                                                |
|                       | R in ( R g ) ◦                                                        | 0 . 20 +0 . 07                                                                                                                                 |                                                                                                    | - 2 . 34 - 0 . 07 18 +11 - 5                                                                                                                   | 6 . 94 +3 . 0 - 9 . 4                                                        | 0 . 66 ± 0 . 01 6 . 1 +2 . 1 - 0 . 1                                                                                            |
|                       | i ( )                                                                 | - 0 . 02                                                                                                                                       |                                                                                                    | 70 +10                                                                                                                                         |                                                                              |                                                                                                                                 |
|                       | - 2                                                                   |                                                                                                                                                |                                                                                                    | - 1 0 . 55 +0 . 09 - 0 . 06                                                                                                                    | 1 . 2 ± 0 . 1                                                                |                                                                                                                                 |
|                       | norm line ( 10 ) EW O (eV)                                            | 5 . 5 +1 . 9                                                                                                                                   | 0 . 4 ± 0 . 1                                                                                      | +2 . 2                                                                                                                                         |                                                                              | 1 . 6 ± 0 .                                                                                                                     |
| DISKLINE              |                                                                       | - 0 . 5 +0 . 28                                                                                                                                | 11 ± 3                                                                                             | 13 . 3 - 1 . 4                                                                                                                                 | 27 ± 2                                                                       | 1 35 ± 2                                                                                                                        |
| 2                     |                                                                       |                                                                                                                                                |                                                                                                    |                                                                                                                                                |                                                                              | +0 .                                                                                                                            |
|                       |                                                                       |                                                                                                                                                |                                                                                                    |                                                                                                                                                | +0 . 05                                                                      |                                                                                                                                 |
|                       | E Fe (keV)                                                            | 6 . 41                                                                                                                                         |                                                                                                    | +0 . 20                                                                                                                                        |                                                                              |                                                                                                                                 |
|                       | norm line ( 10 - 3 )                                                  | - 0 . 01 0 . 2 +0 . 3                                                                                                                          | 6 . 7 ± 0 . 2 0 . 4 +0 . 2                                                                         | 6 . 41 - 0 . 01 0 . 8 ± 0 . 2                                                                                                                  | 6 . 41 - 0 . 01 1 . 5 ± 0 . 2                                                | 6 . 41 21 - 0 . 01 0 . 9 ± 0 . 3                                                                                                |
|                       | EW Fe (eV) χ 2                                                        | - 0 . 1 26 +39 - 13                                                                                                                            | - 0 . 1 75 +37 - 19                                                                                | 142 ± 36                                                                                                                                       | 267 ± 36                                                                     | 143 ± 48                                                                                                                        |
|                       | (dof)                                                                 |                                                                                                                                                |                                                                                                    | 5936.04 (4641)                                                                                                                                 |                                                                              |                                                                                                                                 |

Note.- Errors are reported at the 90% confidence level and calculated from Markov Chain Monte Carlo (MCMC) of chain length 10 6 . 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 χ 2 (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 ( ) and emissivity index ( i q ) are tied between the two disk line components. R in is tied between each DISKLINE component in Model 3.

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 &lt; 5 to avoid periods of time with high particle background due to geomagnetic storms and magnetic cut-off rigidity of COR SAX &gt; 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 NIEXTRACT-EVENTS se-

lecting 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. Lightcurves were extracted in 5 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, colorintensity, 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 ks, 8.45 ks, 10.3 ks, 10.2 ks, 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 XRT observation (ObsID: 00010238010) in Windowed Timing mode on 2017 Aug 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 /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 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 Australia Telescope Compact Array (ATCA) on 2017 Aug 23, from 11:24:59.9 to 16:09:19.9 UTC, with ≈ 3.25 hrs on source (project code CX392). This falls within interval E of the NICER monitoring. The 4cm 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 GHz and 9.0 GHz baseband was imaged separately, with an additional image created combining both datasets in the uv-plane having a central frequency of

1.15

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 . -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).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00539v1-figure-f90a7bf2-158b-410b-8aff-20fcace99ada.jpg)

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 ′′ × 2.1 ′′ , 5.4 ′′ × 1.8 ′′ , and 4.1 ′′ × 1.4 ′′ , respectively. The local sensitivities achieved were: 13 µ Jy beam -1 at 5.5 GHz; 9 µ Jy beam -1 at 7.25 GHz; 10 µ Jy beam -1 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 -1 .

## 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

1.10

Ratio

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.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00539v1-figure-f91e269c-e434-4fb8-b72a-21c99f0009e5.jpg)

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 dot-dashed line corresponds to the single-temperature blackbody. The dotdot-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.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00539v1-figure-fb336f6e-1dee-4657-a5dc-7af917a828ac.jpg)

using the 'chain' command. NICER spectra are modeled in the 0 4 . -9 5 . keV band since 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 since 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

Ratio

0

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.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00539v1-figure-3dc469ab-fc60-464f-976a-ca6606712f26.jpg)

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 χ 2 ν = 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 &amp; Chakrabarty (2003), this improves the fit to χ 2 ν = 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 χ 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 ( χ 2 ν = 3 66 . ). An absorbed cut-off power-law and blackbody model from Madej &amp; Jonker (2011) further improves the overall fit to χ 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 χ 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 &gt; 27 σ level of confidence via an F-test.

Switching out the single-temperature thermal component for Comptonization, NTHCOMP, provided a significantly worse fit ( χ 2 ν = 2 35 . ). Moreover, the Comptonized component tended towards the shape of a blackbody with a high optical depth of τ ∼ 9 . Alternatively, if we switch the powerlaw 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 ( χ 2 ν = 1 75 . ). However, the electron temperature tends to the

Figure 6. Oxygen or iron line flux versus 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 versus the Fe line flux throughout the outburst; the dot-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.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00539v1-figure-1bf8b4f1-d997-4129-a345-ef09c8ccafb0.jpg)

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 χ 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-tonoise 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 . -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 . × 10 -2 keV &lt; σ O &lt; 10 2 . × 10 -2 keV and 0 74 keV . &lt; σ Fe &lt; 1 89 . keV, respectively. The O line width agrees with the values previously reported in Madej &amp; 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 ( ∆ χ 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 since NICER has an energy resolution of ∼ 80 eV at 1 keV and/or partially due to the current uncertainty in the gain ( &lt; 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 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 inten- sity increases from interval A → E. Although the radii inferred from the O feature and Fe feature differ, they agree at the 3 σ level of confidence, again indicating an evolution in R 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 ◦ -80 ◦ , which is consistent with the range of inclination values obtained in Madej &amp; Jonker (2011) and Madej et al. (2014). We also tie R in between each DISKLINE component (Model 3 in Table 1). The fit is 3 9 . σ 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 towards 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 . ± 0 03 . for each model in Table 1. The constant between NICER and INTEGRAL is smaller ( C ∼ 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 since there would be joint energy coverage between all three missions.

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.

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 unab-

sorbed 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 powerlaw 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 ξ = L/nr 2 , where L is the ionizing luminosity from 1 -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 ξ at the beginning and end of the 10 day NICER monitoring period studying the evolution of the O emission features. Figure 13 of Kallman &amp; Bautista (2001) shows the ionization balance of an optically-thin photoionized plasma at high density, as a function of ξ . Although we expect the accretion disk to be optically thick so we cannot use the exact values for ξ , this still allows us to qualitatively estimate how much ξ 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 &amp; Bautista 2001; Kallman et al. 2004), we determine that ∆(log ξ ) ∼ 0 3 . -1 2 . (cgs).

From the limit on ∆(log ξ ) , we can look into the change in the density of the material in the line emitting region using the following equation:

$$\Delta ( \log n ) \simeq \log \left ( \frac { F _ { E } } { F _ { A } } \frac { \xi _ { A } } { \xi _ { E } } \frac { r _ { A } ^ { 2 } } { r _ { E } ^ { 2 } } \right ) \quad \ \ ( 1 ) \quad \ \ & \text{dif} \\ \text{hil}$$

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 ∼ 10 15 -10 19 : Garc´ ıa et al. 2016).

Moreover, we can use the equation for density of a standard thin disk in Frank, King, &amp; Raine (2002) as a cross check for the range of values obtained from Equation 1: ρ = 3 1 . × 10 -8 α -7 / 10 ˙ M 11 / 20 16 m 5 / 8 1 R -15 / 8 10 f 11 / 5 g cm -3 , where α 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 . ∆(log ρ ) and ∆(log n ) should change by the same amount since they are related by a constant. Assuming that the viscosity parameter ( α ) is constant, the change in density becomes:

$$\intertext { \colon a } \text{$\mathfrak { u } $} \Delta ( \log \rho ) \simeq \log \left [ \left ( \frac { F _ { E } } { F _ { A } } \right ) ^ { 1 1 / 2 0 } \left ( \frac { f _ { E } } { f _ { A } } \right ) ^ { 1 1 / 5 } \left ( \frac { r _ { A } ^ { 2 } } { r _ { E } ^ { 2 } } \right ) ^ { - 1 5 / 8 } \right ] _ { ( 2 ) } \\ \intertext { \mathfrak { n } t } \text{ for a NS with a 10 km radius and the same values for flux}$$

For a 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 ξ and strongly on the change in inner disk radius ( ∝ 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 ∆(log n ) ∼ -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 XILLVER CO 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 ( χ 2 ν &gt; 80 ). Therefore, we chose to not include them here. We are currently revising and updating the existing XILLVER CO 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 QPOs seen in stellar mass BHs and a different naming scheme could have been applied (Klein-Wolt &amp; 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 a 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 ( &gt; 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.

Lastly, we searched for coherent pulsations in the data during the 10 day period of NICER monitoring. This would pro-

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.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00539v1-figure-1c5bab50-aa43-455e-a19b-c5ba5a2e3565.jpg)

Figure 8. The change in inner disk radius measured from the O line versus 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 towards the compact object and the QPOs move to higher frequency as the source intensity increases (see Figure 1).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00539v1-figure-25817215-fbaf-49fe-9f9a-be5a2cbe9a32.jpg)

Table 2 Power spectrum fit parameters

|                           | Frequency (Hz)                                          | Quality factor                                        | Fractional amplitude (% rms)                    | χ 2 / dof   |
|---------------------------|---------------------------------------------------------|-------------------------------------------------------|-------------------------------------------------|-------------|
| Interval A                | Interval A                                              |                                                       |                                                 |             |
| break2 break hump low hHz | 0.16(0.03) 0.658(0.019) 3.89(0.07) 32.6(6.5) 143.8(3.9) | 0.28(0.09) 0.47(0.08) 0.59(0.05) 0.0 (fixed) 6.8(3.9) | 6.2(0.8) 12.4(0.6) 16.3(0.5) 12.1(0.6) 4.9(1.0) | 253 / 204   |
| Interval B                | Interval B                                              |                                                       |                                                 |             |
| break hump low hHz        | 0.87(0.05) 5.18(1.4) 42(12) 245(12)                     | 0.08(0.03) 0.58(0.08) 0.0 (fixed) 3.2(1.3)            | 15.5(0.4) 14.6(0.7) 11.2(0.7) 7.1(0.9)          | 224 / 208   |
| Interval C                | Interval C                                              |                                                       |                                                 |             |
| break hump low hHz        | 0.94(0.04) 5.96(1.5) 58.0 (4.0) 157.0 (23)              | 0.10(0.02) 0.56(0.07) 0.0 (fixed) 1.0(1.0)            | 15.5(0.3) 14.5(0.6) 10.8(2.7) 8.0(4.0)          | 264 / 185   |
| break2 break hump         | 0.32(0.06) 1.19(0.04) 7.3(0.3)                          | 0.33(0.10) 0.44(0.10) 0.29(0.06) 0.0 (fixed)          | 5.9(1.1) 12.4(1.0)                              | 233 / 207   |
| Interval D                | Interval D                                              |                                                       |                                                 |             |
| hHz                       | 367(12)                                                 |                                                       | 16.2(0.5) 12.8(1.3)                             |             |
| Interval E                | Interval E                                              |                                                       |                                                 |             |
| break2 break hump hHz     | 0.47(0.08) 1.74(0.06) 9.7(0.4) 223.0 (42)               | 0.27(0.08) 0.54(0.11) 0.33(0.08) 0.4(0.3)             | 6.5(0.9) 11.3(0.9) 14.7(0.6) 11.2(1.2)          | 227 / 206   |

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

vide further evidence that the compact object is a 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 suppresses 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 overabundance 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 &lt; . 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 &amp; 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 towards 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 de-

Figure 9. Radio - X-ray luminosity plane of accreting compact objects adapted from Bahramian et al. (2018). The dashed grey line indicates the best-fit trends for BHs ( L r ∝ L 0 6 . x ; Gallo et al. 2014). The purple dot-dashed lines indicate the proposed correlations for NSs (Migliari &amp; Fender 2006) for atoll and Z sources ( L r ∝ L 0 6 . x ) and for NSs in the hard state ( L r ∝ L 1 4 . x ). The black points are quiescent/hard state BHs from the available literature (Calvelo et al. 2010; Coriat et al. 2011b; Rodriguez et al. 2011; Paizis et al. 2011; Jonker et al. 2012; Corbel et al. 2004, 2013; Brocksopp et al. 2005, 2013; Gallo et al. 2003, 2014; Russell et al. 2015; Plotkin et al. 2017; Rib´ et al. 2017; Dincer et al. 2018). o ¸ The light and dark blue squares are NSs in hard states while orange pentagons are NSs in soft states (Rutledge et al. 1998; Migliari &amp; Fender 2006; Tudose et al. 2009; Migliari et al. 2011; Coriat et al. 2011a; Tetarenko et al. 2016; Tudor et al. 2017; Gusinskaia 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; van den Eijnden et al. 2018a; Strohmayer 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 &amp; 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.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00539v1-figure-023dd2a8-7c9f-4d77-84d4-a13bb6527dd9.jpg)

pendent 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 &amp; 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 . ∼ 1 × 10 -9 ergs s -1 cm -2 ). The inner disk radius inferred from the O line in the final two intervals (D: R in = 7 9 . +2 5 . -1 3 . R g , E: R in = 7 3 . +1 4 . -1 0 . R g ) agree well with the measurements from Madej et al. (2014) using Chandra /LETGS and the initial grid of XILLVER CO ( R in &lt; 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 &amp; Jonker (2011) of 27 -35 eV. Moreover, the inclination of i ∼ 70 ◦ -80 ◦ 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.

central object of 4U 1543 -624 is ultimately unknown. Constraining the mass function of the binary could conclusively 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 the compact accretor is a 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 a 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.

Although often assumed to be a NS, the identity of the

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 Xray flux into luminosities, which depends on the distance to

the source. Wang &amp; 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 &amp; 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 . × 10 38 erg s -1 (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 . +3 3 . -1 9 . kpc using Gaia DR2 parallax measurements with the assumption of an exponentially decreasing space density prior (Bailer-Jones et al. 2015; Astraatmadja &amp; Bailer-Jones 2016). Given the uncertainty surrounding the source location, we choose to encompass the entire range of distance estimates, 1.4 kpc &lt; d &lt; 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 πd νS 2 ν , where S ν = 27 µ Jy, ν = 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 -10 ergs s -1 cm -2 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 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: Russell et al. 2018, Ferrigno 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 &amp; 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 &amp; 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 -1 and L x, 1 -10 keV ∼ 10 36 erg s -1 (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 ≤ 14 5 . × 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 . × 10 -9 ergs s -1 cm -2 in interval E and the upper limit on R in is 8.7 R g . This provides an upper limit of B ≤ 3 2 . × 10 8 G at d = 11 5 . kpc or B ≤ 0 4 . × 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 (Tudor et al. 2017; Tetarenko et al. 2018), 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 Xray relationship of 4U 1543 -624 strongly favors a 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.0. J.A.G. acknowledges support from Chandra Theory grant TM819003X 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.

## REFERENCES

Astraatmadja, T. L., &amp; Bailer-Jones, C. A. L. 2016, ApJ, 832, 137 Arnaud, K. A. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby &amp; J. Barnes, 17 Bahramian, A., Miller-Jones, J. C. A., Strayder, J., et al. 2018, https://doi.org/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., &amp; Fabian A. C., 2002, MNRAS, 336, 867 Belloni, T., Psaltis, D., &amp; van der Klis, M. 2002, ApJ, 572, 392 Blandford, R. D., &amp; Znajek, R. L. 1977, MNRAS, 179, 433. Blandford, R. D., &amp; Payne, D. G. 1982, MNRAS, 199, 883. Bogdanov, S., &amp; Halpern, J. P. 2015, ApJ, 803, L27. 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 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., &amp; Barret, D. 2012, ApJ, 755, 27

Calvelo, D., Fender, R. P., Russell, D. M., et al. 2010, MNRAS, 409, 839 Coriat, M., Tzioumis, T., Corbel, S., et al. 2011a, ATel, 3119, 1 Coriat, M., Corbel, S., Prat, L., et al. 2011b, MNRAS, 414, 677 Corbel, S., Fender, R. P., Tomsick, J. A., Tzioumis, A. K., Tingay, S. 2004,

- ApJ, 617, 1272

Corbel, S., Coriat, M., Brocksopp, C., et al. 2013, MNRAS, 428, 2500 Courvoisier T. J.-L., Walter R., Beckmann V. et al. 2003, A&amp;A, 411, L53 Deller A. T., Moldon, J., Miller-Jones, J. C. A., et al. 2015, ApJ, 809, 13 Dinc ¸er, T., Bailyn, C. D., Miller-Jones, J. C. A., Buxton, M., MacDonald, R. K. D. 2018, ApJ, 852, 4

- Fabian, A. C., Rees, M. J., Stella, L., &amp; White, N. E. 1989, MNRAS, 238, 729
- Fabian A. C., Iwasawa K., Reynolds C. S., &amp; Young A. J., 2000, PASP, 112, 1145
- Federici M., Martino B. L., &amp; Ubertini P. 2010, Experimental Astronomy, Volume 34, Issue 1, pp.105-121

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

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

- Frank, J., King, A., &amp; Raine, D., 2002, Accretion Power in Astrophysics (3rd ed.; Cambridge: Cambridge Univ. Press)
- Gallo, E., Fender, R. P., &amp; Pooley, G. G. 2003, MNRAS, 344, 60

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

290

- Gallo, E., Degenaar, N., &amp; van den Eijnden, J. 2018, MNRAS, 478, L132 Garc´ ıa, J. A., Fabian, A. C., Kallman, T. R., et al. 2016, MNRAS, 462, 751 Gendreau, K. C., Arzoumanian, Z., &amp; Okajima, T. 2012, Proc. SPIE, 8443, 13
- Giacconi, R., Murray, S., Gursky, H., Hellogg, E., Schreier, E., &amp; Tananbaum, H. 1972, ApJ, 178, 281
- Gusinskaia, N. V., Deller, A. T., Hessels, J. W. T., et al. 2017, MNRAS, 470, 1871

Hansinger, G., &amp; van der Klis, M., 1989, A&amp;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., Psaltis D., &amp; Chakrabarty D., 2001, ApJ, 560, L59

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

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

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

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

- Koliopanos, F., Golfanov, M., &amp; Bildsten, L. 2013, MNRAS, 432, 1264 Kuulkers, E., den Hartog, P. R., in't Zand, J. J. M., et al. 2003, A&amp;A, 399, 663

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

Lamb, F. K., Boutloukos, S., Van Wassenhove, S., et al. 2009, ApJ, 706, 427 Lin, D., Remillard, R. A., &amp; Homan, J. 2007, ApJ, 667, 1073 Lin, D., Remillard, R. A., &amp; Homan, J. 2010, ApJ, 719,1350 Ludlam, R. M., Miller, J. M., Miller-Jones, J., &amp; Reynolds, M. 2017, ATel,

- 10690, 1
- 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 Madej, O. K., &amp; Jonker, P. G. 2011, MNRAS, 412, L11
- Madej, O. K., Garc´ ıa, J., Jonker, P. G., et al. 2014, MNRAS, 442, 1157 Matsuoka, M., Kawasaki, K., Ueno, S., et al. 2009, PASJ, 61, 999. McMullin, J. P., Waters, B., Schiebel, D., Young, W., &amp; Golap, K. 2007, ASPC, 376, 127
- Migliari S., Fender R. P., &amp; van der Klis M., 2005, MNRAS, 363, 112 Migliari, S., &amp; Fender, R. 2006, MNRAS, 366, 79
- Migliari, S., Miller-Jones,J. C. A., &amp; Russell, D. M. 2011, MNRAS, 415, 2407
- Miller, J. M., Ludlam, R. M., Reynolds, M. T., Kuulkers, E., Ferringo, C., &amp; Bozzo, E. 2017, ATel, 10719, 1
- Nelemans, G., Jonker, P. G., Marsh, T. R., &amp; van der Klis, M. 2003, MNRAS, 348, L7

Nelemans, G., &amp; Jonker, P. G. 2010, New Astronomy Reviews, 54, 87 Ng C., D´ ıaz Trigo M., Cadolle Bel M., Migliari S., 2010, A&amp;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, Nature, 501, 517 Plotkin, R. M., Bright, J., Miller-Jones, J. C. A., et al. 2017, ApJ, 848, 92 Rib´ o, M., Munar-Adrover, P., Paredes, J. M., et al. 2017, ApJ, 835, 33 Rodriguez, J, Corbel, S., Caballero, I., et al. 2011, A&amp;A, 533, L4 Russell, T. D., Miller-Jones, J. C. A., Curran, P. A., et al. 2015, MNRAS, 450, 1745

- Russell, T., Degenaar, N., Wijnands, R., &amp; van den Eijnden, J. 2018, ATel, 11954, 1
- Rutledge, R., Moore, C., Fox, D., Lewin, W., &amp; van Paradijs, J. 1998, ATel, 8, 1
- Sault, R. J., Teuben, P. J., &amp; Wright, M. C. H. 1995, ASPC, 77, 433 Serino, M., Hiaka, K., Negoro, H., et al. 2018, ATel, 11302, 1 Schultz J., 2003, A&amp;A, 397, 249
- Strohmayer, T.E., Ray, P. S., Gendreau, K. C., et al. 2018a, ATel, 11520, 1 Strohmayer, T.E., Arzoumanian, Z., Bogdanov, S., et al. 2018b, ApJ858, L13 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 Tetarenko, A. J., Bahramian, A., Wijnands, R., et al. 2016, ApJ, 854, 125 Tudose V., Fender R. P., Linares M., Maitra D., &amp; van der Klis M., 2009, MNRAS, 400, 2111

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 Ubertini P., Lebrun F., Di Cocco G. et al. 2003, A&amp;A, 411, L131 Wang, Z., &amp; Chakrabarty, D. 2004, ApJL, 616, L139

- Wang, Z., Tziamtzis, A., Kaplan, D. L., &amp; Chakrabarty, D., 2015, PASA, 32, 35
- 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, Nature, 562, 233
- van der Klis, M. 1989, in NATO ASI Series, Vol. 262, Timing Neutron Stars, ed. H. Ogelman &amp; 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 University Press, p. 39

van Straaten, S., van der Klis, M., &amp; M´ endez, M., 2003, ApJ, 596, 1155

Zampieri, L., Burtovoi, A., Fiori, M., et al. 2019, MNRAS, 485, L109.
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 D. Altamirano
Observations of the Ultra-compact X-ray Binary 4U 1543-624 in Outburst with NICER, INTEGRAL, Swift, and ATCA
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 NICER monitoring of the source over a $\sim10$ day period during which
target-of-opportunity observations were also conducted with Swift, INTEGRAL,
and ATCA. Emission lines were measured in the NICER X-ray spectrum at
$\sim0.64$ keV and $\sim6.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 (NS) primary. This 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 alone.