Skip to content
The following article is Open access

Insight-HXMT, NICER, and NuSTAR Views to the Newly Discovered Black Hole X-Ray Binary Swift J151857.0–572147

, , , , , , , , ,

Published 2024 September 12 © 2024. The Author(s). Published by the American Astronomical Society.
, , Citation Jing-Qiang Peng et al 2024 ApJL 973 L7DOI 10.3847/2041-8213/ad74ec

Download Article PDF
DownloadArticle ePub

You need an eReader or compatible software to experience the benefits of the ePub3 file format.

2041-8205/973/1/L7

Abstract

The systematic properties are largely unknown for the black hole X-ray binary Swift J151857.0–572147 newly discovered in the 2024 outburst. The nature of a black hole can be completely defined by specifying the mass and dimensionless spin parameter. Therefore, accurate measurement of the two fundamental parameters is important for understanding the nature of black holes. The joint spectral fitting of a reflection component with simultaneous observations from Insight-HXMT, NICER, and NuSTAR reveals for the first time a black hole dimensionless spin of ${0.84}_{-0.26}^{+0.17}$ and an inclination angle of ${21.1}_{-3.6}^{+4.5}$ degrees for this system. Monitoring of the soft state by NICER results in disk flux and temperature following ${F}_{\mathrm{disk}}\propto {T}_{\mathrm{in}}^{3.83\pm 0.17}$. For the standard thin disk, ${L}_{\mathrm{disk}}\approx 4\pi {R}_{\mathrm{in}}^{2}\sigma {T}_{\mathrm{in}}^{4}$, so the relationship between the flux and temperature of the disk we measured indicates that the inner radius of the disk is stable and the disk is in the innermost stable circular orbit. With an empirical relation built previously between the black hole outburst profile and the intrinsic power output, the source distance is estimated as 5.8 ± 2.5 kpc according to the outburst profile and peak flux observed by Insight-HXMT and NICER. Finally, a black hole mass of 3.67 ± 1.79–8.07 ± 4.20M can be inferred from a joint diagnostic of the aforementioned parameters measured for this system. This system is also consistent with most black hole X-ray binaries with high spin and a mass in the range of 5–20 M.

Export citation and abstractBibTeXRIS

Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence. Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.

1. Introduction

Black hole (BH) low-mass X-ray binaries consist of a companion star with a mass less than 1M and a BH. Most of the BH low-mass X-ray binaries can only be observed by X-ray telescopes when they enter into outbursts, during which their different spectral states can be observed and traced in the hardness–intensity diagram (HID; J. Homan et al. 2001; R. P. Fender et al. 2004; S. Motta et al. 2012). For the low hard state (LHS), the emission is dominated by nonthermal emission from the corona/jet, accompanied by a relatively low percentage of thermal emission from the disk (T. Belloni et al. 2005). The accretion disk is commonly thought to be truncated, but as of now, there is no definitive observational evidence that the accretion disk is truncated in the LHS (J. M. Miller et al. 2006; M. T. Reynolds & J. M. Miller 2013; P. A. Draghis et al. 2024). As the accretion rate increases, the source passes through the intermediate states (hard and soft intermediate states) to the high soft state (HSS). In the HSS, the inner edge of the disk is expected to lie at the innermost stable circular orbit (ISCO) and dominate the emission in the form of thermal. The nonthermal emission in HSS is relatively weak, usually accounting for less than 25% of the total (A. A. Esin et al. 1997; M. Gierlinski et al. 2008).

The nature of a black hole can be completely defined by specifying two parameters: mass M and dimensionless spin a. Therefore, accurate measurement of the two fundamental parameters is important for understanding the nature of black holes. The spin of a black hole fundamentally changes the geometry of spacetime. The gravitational well of an extreme Kerr black hole of the same mass is significantly deeper compared to a spinless black hole, resulting in a much harder X-ray spectrum and significantly enhancing its efficiency in converting the accreted rest mass into radiant energy. And the spin of black holes plays an important role in the mechanism of jets, core collapse of gamma-ray bursts, gravitational-wave astronomy in predicting the waveforms of merging black holes, and so on (S. E. Woosley 1993; R. D. Blandford & R. L. Znajek 1977; S. N. Zhang et al. 1997; J. E. McClintock et al. 2006; M. Campanelli et al. 2006; J. F. Steiner et al. 2011). The mass of a black hole provides a physical scale. The mass distribution of such black holes can provide important clues to the end stages of the evolution of massive stars.

The BH spin can be inferred from fitting the continuum and reflection spectra. The continuum spectrum approach requires an inner disk to stay around ISCO, independent mass, inclination, and distance measurements, and an assumption about the color correction factor, and many models can be used to estimate BH spin, e.g., kerrbb and kerrbb2 (S. N. Zhang et al. 1997; L.-X. Li et al. 2005), while the reflection-spectrum approach deals with the broadened iron line and the Compton hump around 20 keV (L. W. Brenneman & C. S. Reynolds 2006; J. M. Miller et al. 2009), commonly used models are the relxill family, etc. (T. Dauser et al. 2016). The BH mass measurement is usually a tough task. The most reliable BH mass estimation comes from the mass function obtained from observing the orbital modulation of the line emission from the companion (V. Kalogera & G. Baym 1996; B. E. Tetarenko et al. 2016). BH mass can also be constrained from the merger of the compact objects, and BHs discovered by the Laser Interferometer Gravitational-Wave Observatory and its Virgo interferometer (LIGO/Virgo) cover a rather wide mass range of 3–100M (J. Bodensteiner et al. 2022). For other approaches, J.-Q. Peng et al. (2023, 2024) reported the BH mass of SLX 1746–331 by assuming its HSS reaches 30% LEdd, which results in a BH mass similar to that derived from an empirical relation of BH mass with disk temperature.

Swift J151857.0–572147 is a new BH X-ray binary first observed by SWIFT/BAT. A spectral fit to the combined Swift/XRT data shows a power-law slope with a photon index of 1.78 ± 0.02 (J. A. Kennea et al. 2024). F. J. Cowie et al. (2024) found that the radio spectrum of the source is inverted with a preliminary spectral index of α ∼ +0.5, where flux density as a function of frequency (fν ) is proportional to να . This inverted index and the reported X-ray photon index of +1.78 are consistent with an X-ray binary in the hard state, suggesting that Swift J151857.0–572147 is either a black hole or (radio-bright) neutron star X-ray binary. An extremely bright radio flare was detected by ATCA on March 9, which is usually suggestive of having discrete ejecta from a BH X-ray binary during the transition between hard and soft X-ray states (F. Carotenuto & T. D. Russell 2024). Subsequently, Swift/XRT observation of Swift J151857.0–572147 confirmed a likely transition to the soft state on March 10 (M. Del Santo et al. 2024). The H i spectroscopic distance constraints for Swift J151857.0–572147 give a rather large range of 4.48–15.64 kpc (B. J. Burridge et al. 2024). By taking a distance of 10 kpc, S. Mondal et al. (2024) fitted the joint IXPE and NuSTAR spectra and reported rather broad parameter coverages: 0.6 ± 0.2–0.7 ± 0.2 for BH spin, 38° ± 9°−47° ± 15° for disk inclination, and 9.2 ± 1.6–10.5 ± 1.8M for BH mass. K. Chatterjee et al. (2024) detected type C quasiperiodic oscillations (QPOs) in Swift J151857.0–572147 with Insight-HXMT data.

With the observations from Insight-HXMT, NICER, and NuSTAR, we perform the first detailed investigation of the systematic properties of Swift J151857.0–572147. In Section 2, we describe the observations and data reduction. The detailed results are presented in Section 3. The results are then discussed, and the conclusions are presented in Section 4.

2. Observations and Data Reduction

2.1. NICER

NICER was launched by the Space X Falcon 9 rocket on 2017 June 3 (K. C. Gendreau et al. 2016). NICER has a large effective area and high temporal resolution in the soft X-ray band (0.2–12 keV), which may allow us to study blackbody components at low temperatures better. As shown in Table 1, we selected NICER observations during the soft state of Swift J151857.0–572147, covering a time zone from MJD 60387–MJD 60433.

We use the official software HEASOFT V6.33/NICERDAS v012 with the latest CALDB xti20240216 for NICER data analysis. NICER data are reduced using the standard pipeline tool nicerl2. 8 We extract light curves using nicerl3-lc 9 in 1–4, 4–10, and 1–10 keV. To extract the spectrum, we utilize nicerl3-spect, 10 employing the “nibackgen3C50 11 model to estimate the background for spectral analysis. We select an energy range of 2–10 keV for spectrum fitting. Additionally, nicerl3-spect automatically applies the systematic error using niphasyserr. In the energy range of 0.3–10 keV, the systematic error is about 1.5%.

Table 1. NICER Observations of Swift J151857.0–572147 during the 2024 Outburst

NICERObserved DateExposure Time
ObsID(MJD)(s)
720422011160387.554193
766101010160404.183907
720422011260405.22569
720422011360408.90228
720422011460409.16404
766101010360412.062682
720422011560412.77262
720422011660413.67799
720422011760414.061411
720422011860415.411065
720422011960416.251265
720422012060417.29565
766101010460417.413218
766101010560425.092474
766101010760433.022021

Download table as:  ASCIITypeset image

2.2. NuSTAR

NuSTAR is the first mission to employ two focusing techniques for detecting hard X-rays at 3–79 keV. It was launched on 2012 June 13, at 9 am PDT. (F. A. Harrison et al. 2013). As shown in Table 2, NuSTAR conducted two observations of Swift J151857.0–572147. For NuSTAR data analysis, we employ the software HEASOFT V6.33/NuSTARDAS v2.1.2 with the calibration database (CALDB 20230613) to reduce the data. We extract NuSTAR-filtered data using the standard pipeline program nupipeline. The spectrum is extracted from a 120″ circle region centered on the source and the background is generated from a 60″ circle region away from the source. We utilized FPMA and FPMB 3–78 keV data for spectral analysis.

Table 2. The Quasi-simultaneous NICER, NuSTAR, and Insight-HXMT Observations of Swift J151857.0–572147 during the 2024 Outburst

NICERObserved DateExposure Time
ObsID(MJD)(s)
720422011160387.54193
Insight-HXMTObserved DateExposure TimeExposure Time
  MEHE
ObsID(MJD)(s)(s)
P06143740080260386.737041925
NuSTARObserved DateExposure TimeExposure Time
  FPMAFPMB
ObsID(MJD)(s)(s)
9100131100260386.713,88014,300
9100131100460387.691959449

Download table as:  ASCIITypeset image

2.3. Insight-HXMT

Insight-HXMT, the first Chinese X-ray astronomy satellite, was successfully launched on 2017 June 15 (S. Zhang et al. 2014, 2018; S.-N. Zhang et al. 2020). The satellite is equipped with three scientific payloads: the low-energy X-ray telescope (LE) featuring an SCD detector operating in the 1–15 keV range with an effective area of 384 cm2 (Y. Chen et al. 2020), the medium-energy X-ray telescope (ME) equipped with a Si-PIN detector covering 5–35 keV with an effective area of 952 cm2 (X. Cao et al. 2020), and the high-energy X-ray telescope (HE) utilizing a phoswich NaI(CsI) detector for the 20–250 keV energy range with an effective area of 5100 cm2 (C. Liu et al. 2020).

The Insight-HXMT began observing Swift J151857.0–572147 on 2024 March 4 and ended on March 17. However, due to contamination by Cir X–1, which is located in the field of view at a position close to Swift J151857.0–572147, only the Insight-HXMT HE and ME observations simultaneous with NICER are adopted for joint spectral analysis (Table 2). This is because during the joint observations, Cir X–1 was in dip phase with a rather soft energy spectrum, and hence the expected count rates at ME and HE bands are rather small.

We extract the data from LE, ME, and HE using the Insight-HXMT Data Analysis software HXMTDAS v2.06. The data are filtered with the criteria recommended by the Insight-HXMT Data Reduction Guide v2.06. 12 Due to the relatively high impact of Cir X–1 on LE, the energy bands considered for spectral analysis are ME 8–28 keV and HE 28–100 keV. A 1% systematic error is added to data (J.-Y. Liao et al. 2020), and errors are estimated via Markov Chain Monte Carlo chains with a length of 20,000.

3. Results

3.1. Light Curve and Hardness–Intensity Diagram

Figure 1 shows the light curves of Swift J151857.0–572147 observed by Insight-HXMT and NICER during the 2024 outburst. The red dots represent Insight-HXMT observations and the black dots represent NICER observations. We normalized the light curves to the Crab Nebula by using the count rate ratios between Insight-HXMT and NICER. Since Swift J151857.0–572147 is 0fdg3 away from Cir X–1, both sources are enclosed in the field of view of Insight-HXMT. Fortunately, Cir X–1 has a relatively stable profile for flux evolution over its orbital period of 16.6 days, the overall flux evolution trend at the LE band can be constructed with Insight-HXMT observations outside the outburst of Swift J151857.0–572147. So we extract the Cir X–1 LE light curve at LE 2–12 keV, then normalize it with the counts of Crab in the same energy band and fit it to get the flux evolution trend (see the gray shading in the bottom panel of Figure 1). Then, we subtract off the Cir X–1 and get the LE light curve of the Swift J151857.0–572147 at 2–12 keV.

Figure 1. Refer to the following caption and surrounding text.

Figure 1. The light curves of Swift J151857.0–572147 observed by Insight-HXMT and NICER during the 2024 outburst. The red dots represent Insight-HXMT observations and the black dots represent NICER observations. Top panel: the light curve of Insight-HXMT LE and NICER in 2–12 keV. Middle panel: the light curves are normalized to the Crab. Bottom panel: the light curve of Insight- HXMT after subtracting off the Cir X–1 (a gray area) and the light curve of NICER. The red dotted lines correspond to about 90%, 50%, and 10% of the peak flux, respectively.

Standard image High-resolution image

The HID is built with NICER light curves in 1–4, 4–10, and 1–10 keV. The hardness is defined as the count rate ratio of 4–10 to 1–4 keV, while the intensity takes the count rate of 1–10 keV. As shown in Figure 2, the source was in a soft state during almost the whole NICER observations.

Figure 2. Refer to the following caption and surrounding text.

Figure 2. The NICER HID of Swift J151857.0–572147, where the hardness is defined as the ratio of 4–10 to 1–4 keV count rate. The subplot shows the distribution of the outburst trajectory of Swift J151857.0–572147 after scaling the count rate relative to the trajectory of the 2021 outburst of GX 339–4.

Standard image High-resolution image

3.2. The Spectral Analysis

As shown in Table 2, there are two sets of simultaneous joint observations from Insight-HXMT, NICER, and NuSTAR. One set is from Insight-HXMT and NuSTAR at the time around MJD 60386 and another set is from NICER and NuSTAR at the time around MJD 60387. Accordingly, joint spectral fittings are carried out for each of the observational sets.

We account for interstellar absorption with the tbabs model (J. Wilms et al. 2000), where photoelectric cross sections were given in D. A. Verner et al. (1996). We fit the multitemperature blackbody component of the accretion disk with diskbb (K. Mitsuda et al. 1984). Cutoffpl is used to fit nonthermal emission. Constant is introduced to balance the calibration discrepancies between different telescopes. In order to take into account the differences between the different spectral indices of NICER and NuSTAR, we allow the temperature of the disk and spectral index to vary. At this step, we find a poor-fitting result with χ2/(degrees of freedom (dof)) = 1.34, and a significant reflection component in the residuals (Figure 3). Therefore, we replace cutoffpl with the relxill V2.3 model to fit the reflectance component of the spectrum (T. Dauser et al. 2016), which results in a much improved fit, χ2/(dof) = 1.06 (Figure 4). Therefore, our fitting model M1 is: constanttbabs(diskbb+relxill). We fix Rin at -1 and the disk at ISCO. The break radius and the outer radius of the disk are fixed at 15 and 400, respectively, and both emission indices are fixed at 3. The spectral parameters are shown in Table 3. The BH spin and disk inclination averaged over the two observational sets are ${0.84}_{-0.26}^{+0.17}$ and ${21.1}_{-3.6}^{+4.5}$ degrees, respectively. The reflection fraction is relatively small and is defined in the frame of the primary source as the ratio of intensity emitted toward the disk compared to escaping to infinity. This indicates that the corona is less illuminated by the disk. The normalization of relxill is also smaller, indicating a lower percentage of reflected fluxes. The ionization of the accretion disk and the iron abundance of the material in the accretion disk are also higher, which is related to the lower electron number density (ne = 1015 cm−3) of the relxill. This relatively high metal abundance has also appeared in other BH X-ray binaries such as MAXI J1820 + 070, GX 339–4, and Cyg X–1 (F. Fürst et al. 2015; M. L. Parker et al. 2015; B. You et al. 2021). We replace relxill with the reflection model relxillcp for higher-density disks, so our model M1* is: constant*tbabs(diskbb+relxillcp). The parameters are shown in Table 4. The BH spin and disk inclination averaged over the two observational sets are ${0.87}_{-0.20}^{+0.06}$ and ${21.4}_{-1.9}^{+2.2}$ degrees, respectively, which is consistent with the results obtained for M1. Both the ionization parameters and the iron abundance have also decreased, but the iron abundance is still high because the density of the disk is ne = 1015 cm−3, which is lower than the typical values in the standard thin disk model (N. I. Shakura & R. A. Sunyaev 1973).

Figure 3. Refer to the following caption and surrounding text.

Figure 3. The simultaneous broadband spectrum with the model (constanttbabs(diskbb+cutoffpl)) of Swift J151857.0–572147 is observed from NuSTAR/FPMA (black), NuSTAR/FPMB (red), Insight-HXMT ME (green), and Insight-HXMT HE (blue).

Standard image High-resolution image
Figure 4. Refer to the following caption and surrounding text.

Figure 4. The simultaneous broadband spectrum with the model M1 of Swift J151857.0–572147 is observed from NuSTAR/FPMA (black), NuSTAR/FPMB (red), Insight-HXMT ME (green), and Insight-HXMT HE (blue).

Standard image High-resolution image

Table 3. The Results of Spectral Fitting the Insight-HXMT+NuSTAR and NICER+NuSTAR Data for Model M1

ModelParameter91001311002+HXMTNICER+91001311004
tbabs NH[1022 cm−2] ${3.32}_{-0.01}^{+0.02}$ ${3.37}_{-0.10}^{+0.13}$
diskbb Tin ${0.90}_{-0.01}^{+0.01}$ ${0.90}_{-0.01}^{+0.01}$
 norm ${1320.5}_{-36.2}^{+17.1}$ ${1615.5}_{-113.3}^{+64.4}$
relxill a ${0.86}_{-0.14}^{+0.12}$ ${0.82}_{-0.22}^{+0.13}$
  i [°] ${20.2}_{-0.8}^{+0.5}$ ${21.9}_{-3.5}^{+4.5}$
 Γ ${2.56}_{-0.01}^{+0.01}$ ${2.63}_{-0.01}^{+0.01}$
  Rin −1 a −1 a
  Rout 400 a 400 a
  Rbr 15 a 15 a
 index13 a 3 a
 index23 a 3 a
 logxi ${4.38}_{-0.06}^{+0.10}$ ${4.02}_{-0.21}^{+0.11}$
  Afe ${5.02}_{-0.23}^{+0.13}$ ${3.87}_{-1.53}^{+0.87}$
  Ecut ${111.1}_{-3.3}^{+7.3}$ ${122.2}_{-12.8}^{+21.7}$
 fefl_frac ${0.16}_{-0.01}^{+0.02}$ ${0.21}_{-0.03}^{+0.04}$
 norm[10−2] ${7.67}_{-0.11}^{+0.12}$ ${6.39}_{-0.14}^{+0.17}$
constantcon[NICER]...1 a
 con[NFPMA]1 a ${1.61}_{-0.06}^{+0.06}$
 con[NFPMB] ${0.97}_{-0.01}^{+0.01}$ ${1.57}_{-0.06}^{+0.06}$
 con[ME] ${0.94}_{-0.01}^{+0.01}$
 con[HE] ${0.82}_{-0.04}^{+0.11}$
  χ2/(dof)1.061.01

Note.

a These parameters are fixed during the fitting.

Download table as:  ASCIITypeset image

Table 4. The Results of Spectral Fitting the Insight-HXMT+NuSTAR and NICER+NuSTAR Data for Model M1*

ModelParameter91001311002+HXMTNICER+91001311004
tbabs NH[1022 cm−2] ${2.98}_{-0.09}^{+0.12}$ ${4.44}_{-0.01}^{+0.03}$
diskbb Tin ${0.95}_{-0.01}^{+0.02}$ ${0.84}_{-0.01}^{+0.01}$
 norm ${1000.2}_{-82.3}^{+53.5}$ ${2242.6}_{-100.9}^{+39.7}$
relxill a ${0.87}_{-0.33}^{+0.11}$ ${0.86}_{-0.23}^{+0.06}$
  i [°] ${20.4}_{-2.7}^{+3.1}$ ${22.3}_{-2.7}^{+3.0}$
 Γ ${2.54}_{-0.02}^{+0.02}$ ${2.59}_{-0.01}^{+0.02}$
  Rin −1 a −1 a
  Rout 400 a 400 a
  Rbr 15 a 15 a
 index13 a 3 a
 index23 a 3 a
 logxi ${2.30}_{-0.27}^{+0.03}$ ${3.30}_{-0.02}^{+0.16}$
 logN[cm−3] ${19.00}_{-0.01}^{+0.03}$ ${19.00}_{-0.14}^{+0.03}$
  Afe ${4.02}_{-0.23}^{+0.13}$ ${3.94}_{-0.153}^{+0.24}$
  kTe[keV] ${61.5}_{-13.2}^{+15.1}$ ${67.2}_{-12.8}^{+21.7}$
 fefl_frac ${0.19}_{-0.05}^{+0.06}$ ${0.10}_{-0.02}^{+0.01}$
 norm[10−2] ${6.02}_{-0.28}^{+0.29}$ ${3.49}_{-0.03}^{+0.04}$
constantcon[NICER]...1 a
 con[NFPMA]1 a ${1.25}_{-0.05}^{+0.01}$
 con[NFPMB] ${0.97}_{-0.01}^{+0.01}$ ${1.22}_{-0.05}^{+0.01}$
 con[ME] ${0.94}_{-0.01}^{+0.01}$
 con[HE] ${0.82}_{-0.04}^{+0.11}$
  χ2/(dof)1.061.01

Note.

a These parameters are fixed during the fitting.

Download table as:  ASCIITypeset image

Since the NICER observations are in the soft state of the outburst, we fit the NICER spectrum by removing the relxill in M1. We find that the fit is poor, with χ2/(dof) = 1.41, due to the existence of a significant Compton component in the residuals. By introducing an additional thcomp model to account for the Compton component, the fit is largely improved with χ2/(dof) = 0.85 (Figure 5). So our model M2 is tbabsthcompdiskbb and the fluxes of the disk and nonthermal in the 0.001–100 keV are estimated with cflux. The evolution of the spectral parameters is shown in Figure 6, where the percentage of nonthermal emission is as low as less than 10%, consistent with having a soft state of BHXRBs. Additionally, the normalization of diskbb remains almost unchanged during the outburst.

Figure 5. Refer to the following caption and surrounding text.

Figure 5. The spectral fittings for the data born out of NICER observations of Swift J151857.0–572147.

Standard image High-resolution image
Figure 6. Refer to the following caption and surrounding text.

Figure 6. Evolution of the spectral parameters of Swift J151857.0–572147 in 2024 from model M2: NH is the column density, Γ the low-energy power-law photon index, Tin the temperature of the inner disk, Cov_frac the coverage factor, and rin the inner radius of the disk. Fdisk and Fnon−thermal are disk flux and nonthermal flux, respectively.

Standard image High-resolution image

3.3. Properties of the System and the Compact Object

As shown in Figure 7, we fit the temperature and flux of the disk obtained with M2 using the formula ${F}_{\mathrm{disk}}={\rm{norm}}* {T}_{\mathrm{in}}^{\alpha }$, and we find that ${F}_{\mathrm{disk}}\propto {T}_{\mathrm{in}}^{3.83\pm 0.17}$. For the standard thin disk, the luminosity of the disk is approximately given by ${L}_{\mathrm{disk}}\approx 4\pi {R}_{\mathrm{in}}^{2}\sigma {T}_{\mathrm{in}}^{4}$, where Rin, Tin, and σ represent the inner radius, the temperature of the disk, and the Stefan–Boltzmann constant, respectively. So the formula between the flux and temperature of the disk we measured indicates that the inner radius of the disk is stable and the disk is in the ISCO during the HSS of the outburst. From diskbb, we have $\mathrm{norm}={\left(\tfrac{{R}_{\mathrm{in}}}{{D}_{10}}\right)}^{2}\cos \theta $, where Rin is the inner disk radius, D10 is the distance of the source in units of 10 kpc, and θ is the inclination angle. With an inclination angle of ${21.1}_{-3.6}^{+4.5}$ degrees as aforementioned in the spectral analysis, the inner disk radius can be deduced as ${R}_{\mathrm{in}}={R}_{\mathrm{ISCO}}={53.2}_{-1.08}^{+1.89}{D}_{10}$ km. However, the previous rough estimation of the distance ${4.48}_{-0.47}^{+0.67}-{15.64}_{-0.60}^{+0.77}$ kpc prevents further constraining the Rin.

Figure 7. Refer to the following caption and surrounding text.

Figure 7. The disk unabsorbed flux (0.01–100 keV) versus the disk inner temperature (Tin). Flux is found to vary with temperature in form ${T}_{\mathrm{in}}^{3.83\pm 0.17}$.

Standard image High-resolution image

For the distance measurements, Z. Yan & W. Yu (2015) have statistically investigated the properties of 110 bright X-ray outbursts out of 36 low-mass X-ray binary transients observed by the All-Sky Monitor on RXTE from 1996 to 2011. They derived a number of outburst characteristics, including peak X-ray luminosity, change rate of luminosity on a daily timescale, e-folding rise and decay timescales, outburst duration, and total radiated energy. They built the relationships between the total radiated energy and the peak luminosity as $\mathrm{log}E\,=(-13.26\pm 6.53)+(1.52\pm 0.17)\times \mathrm{log}{L}_{\mathrm{peak}}$, and the relationship between the total radiated energy and the e-folding rise and decay timescales as ${\rm{log}}E=(43.46\pm 0.10)\,+(0.68\pm 0.11)\times {\rm{log}}{\tau }_{{\rm{rise}},10{\rm{ \% }}-90{\rm{ \% }}}$, ${\rm{log}}E=(43.72\pm 0.10)+(0.49\pm 0.13)\times {\rm{log}}{\tau }_{{\rm{rise}},10{\rm{ \% }}-50{\rm{ \% }}}$, ${\rm{log}}E=(43.48\pm 0.12)+(0.60\pm 0.12)\times {\rm{log}}{\tau }_{{\rm{rise}},50{\rm{ \% }}-90{\rm{ \% }}}$, ${\rm{log}}E=(42.65\pm 0.12)+(1.14\pm 0.10)\times {\rm{log}}{\tau }_{{\rm{decay}},10{\rm{ \% }}-90{\rm{ \% }}}$, ${\rm{log}}E=(42.65\pm 0.12)+(1.00\pm 0.11)\times {\rm{log}}{\tau }_{{\rm{decay}},10{\rm{ \% }}-50{\rm{ \% }}}$, and ${\rm{log}}E=(43.09\pm 0.13)+(0.75\pm 0.10)\times {\rm{log}}{\tau }_{{\rm{decay}},50{\rm{ \% }}-90{\rm{ \% }}}$. For Swift J151857.0–572147, we choose F90%, F50%, and F10% from the light curve where the fluxes are about 90%, 50%, and 10% of the peak flux (red dashed lines in Figure 1). The three e-folding decay timescales are then calculated using the formula $\tau ={\rm{\Delta }}t/\mathrm{ln}A$ (Z. Yan & W. Yu 2015) and the corresponding three total radiated energies are obtained from the relationship between the total radiated energy and the e-folding decay timescales (Table 5). The peak flux of Swift J151857.0–572147 is about 1.41 Crab, corresponding to a flux of about 2.6 × 10−8erg cm−2 s−1. As shown in Table 5, we can estimate the three distances based on the three total radiated energies obtained. Accordingly, the average distance can be deduced as 5.8 ± 2.5 kpc. With this distance, we have ${R}_{\mathrm{in}}={30.8}_{-13.3}^{+13.3}$ km. Since the inner disk shall hold at ${R}_{\mathrm{in}}={2.6}_{-0.04}^{+0.04}{R}_{{\rm{g}}}$ for a BH spin of ${0.84}_{-0.26}^{+0.17}$, the mass of the BH harbored in Swift J151857.0–572147 can be finally inferred as ${8.07}_{-3.69}^{+3.70}{M}_{\odot }$.

Table 5. The Parameters Essential for Distance Estimation Obtained by Following the Approach in Z. Yan & W. Yu (2015) for Swift J151857.0–572147 during the 2024 Outburst

  tpeak Fpeak τdecay,10−90% τdecay,10%−50% τdecay,50%−90% E D
 (MJD)(Crab)(days)(days)(days)(1044 erg)kpc
Swift J151857.0-572147603771.14 ± 0.0120.97 ± 0.2326.23 ± 0.4310.77 ± 0.331.43 ± 0.075.94 ± 4.35
 ...............2.18 ± 0.186.74 ± 4.94
 ...............0.73 ± 0.034.72 ± 3.36

Download table as:  ASCIITypeset image

4. Discussion and Conclusion

We have carried out a spectral analysis of the Swift J151857.0–572147 outburst using data from NICER, Insight-HXMT, and NuSTAR. Almost the entire fundamental parameters (e.g., the BH spin, BH mass, disk inclination, source distance, etc.) of this newly discovered X-ray binary (XRB) system are largely uncovered for the first time.

As a newly discovered BH low-mass X-ray binary system, the fundamental parameters of Swift J151857.0–572147 are largely unknown despite a few reports claiming some rough estimations. To simultaneously constrain the spin, inclination, and mass solely with kerrbb is usually challenging and heavily influenced by whether the distance parameter is properly set. BH spin and disk inclination can be measured in the absence of the knowledge of the BH mass and the source distance, in an approach by fitting the disk reflection feature and the broadened iron line. Fortunately, the joint Insight-HXMT, NICER, and NuSTAR observations of Swift J151857.0–572147 provide data with high quality and broad energy coverage, essential for disentangling in the spectrum the reflection component. As a result, the BH spin and the disk inclination Swift J151857.0–572147 are precisely refined as ${0.84}_{-0.26}^{+0.17}$ and ${21.1}_{-3.6}^{+4.5}$ degrees, respectively. As shown in Table 3, an adoption of the flavor of relxill can well account for the reflection component. We notice that an alternative flavor of relxillcp can also fit the spectrum and have little influence on the measurements of BH spin and inclination angle. An unusually high abundance is required from the measured Afe, which is also reported in outbursts of some other BH XRB systems of, e.g., MAXI J1820 + 070, GX 339–4, and Cyg X–1 (F. Fürst et al. 2015; M. L. Parker et al. 2015; B. You et al. 2021). We try a version of relxill with a high disk density (relxllcp) and find the strange abundance value is an artifact of assuming a low density. Another uncertainty comes from the manner we eliminate the calibration difference by allowing the parameters to vary across instruments of NICER and NuSTAR. We try as well instead with constant and find that, while the fitting gets worse, the BH spin and disk inclination can be inferred with fewer changes. Hence, we think the BH spin and inclination angle are measured more or less robust.

Swift J151857.0–572147 was observed to have an outburst on 2024 March 3 (MJD 60372), and enter a soft state on 2024 March 10 (MJD 60379) (J. A. Kennea et al. 2024; M. Del Santo et al. 2024). The peak flux was caught by Insight-HXMT and the soft state was monitored by NICER. Our spectral analysis of the NICER data reveals that in the soft state, the disk flux is proportional to ${T}_{\mathrm{in}}^{3.83\pm 0.17}$, indicating an inner disk staying around ISCO. Since with a BH spin of ${0.84}_{-0.26}^{+0.17}$ and a disk normalization of ${2628.9}_{-90.2}^{+132.9}$ for Swift J151857.0–572147 we have ${R}_{\mathrm{in}}={R}_{\mathrm{ISCO}}={2.6}_{-0.4}^{+0.4}{R}_{{\rm{g}}}={53.2}_{-1.08}^{+1.59}{D}_{10}$ km, the BH mass can in principle be directly inferred once the source distance is known. A distance of ${4.48}_{-0.47}^{+0.67}-{15.64}_{-0.60}^{+0.77}$ kpc claimed previously in B. J. Burridge et al. (2024) highlights the huge uncertainty in distance measurement. As a result, BH mass can fall in a range of ${6.24}_{-1.17}^{+1.35}-{21.77}_{-3.48}^{+3.58}{M}_{\odot }$, with a notable range of uncertainty. Obviously, refined source distance would be essential for further probing the mass of the BH harbored in Swift J151857.0–572147.

Z. Yan & W. Yu (2015) conducted a statistical analysis of 110 bright X-ray outbursts from 36 low-mass X-ray binary transients observed by the All-Sky Monitor on RXTE between 1996 and 2011. They measured several outburst characteristics, including peak X-ray luminosity, daily rate of change of luminosity, e-folding rise and decay timescales, outburst duration, and total radiated energy. Since the peak X-ray luminosity is a function of disk mass, which is expected from the disk instability model (A. R. King & H. Ritter 1998; J.-P. Lasota 2001), disk mass plays a major role in determining outburst properties (W. Yu & J. Dolence 2007). If one can assume as an approximation that the radiation efficiency is constant, the total radiant energy E should correspond to the disk mass accreted during the outburst. Therefore, they found a positive correlation between the total radiative energy and the peak X-ray luminosity, as well as a positive correlation between the total radiative energy and the timescale of e-folding rise or decay in the outburst. With these relations, distances were estimated for SLX 1746–331, SWIFT J1539.2–6227, SWIFT J1842.5–1124, XTE J1755–324, and XTE J2012 + 381. By taking the approach in Z. Yan & W. Yu (2015), and the outburst profile monitored jointly by Insight-HXMT and NICER at 2–12 keV, a distance of 5.8 ± 2.5 kpc can be assigned to Swift J151857.0–572147. Finally, the joint diagnostic upon the systematic properties of Swift J151857.0–572147 ends with a BH mass measured as${8.07}_{-3.69}^{+3.70}{M}_{\odot }$.

We would like to note that such an estimation of BH mass may suffer from a few aspects in measuring the inner disk radius with a simplistic version of diskbb. (1) The diskbb model does not include a zero-torque inner boundary condition. Using a model that includes this important physical effect tends to reduce radii by a factor of about 2.2 (see, e.g., E. R. Zimmerman et al. 2005), which may reduce the BH mass from ∼${8.07}_{-3.69}^{+3.70}{M}_{\odot }$ to ${3.67}_{-1.68}^{+1.68}{M}_{\odot }$. (2) There is a color correction factor that aims to account for scattering in the disk atmosphere. The color correction factor fcol is 1.7 ± 0.2 and is minimally affected by mass and accretion rate, indicating a weak dependence of fcol on black hole spin a (T. Shimura & F. Takahara 1995). (3) Additionally, there are also relativistic corrections to the diskbb results. S. N. Zhang et al. (1997) found that the uncertainty in fcol is 10%, which will cause an error of no more than 20% in the inner radius of the disk. Therefore, the final conservative estimation of BH mass is 3.67 ± 1.79–8.07 ± 4.20M. We find that, with such a BH mass and aforementioned distance, the outburst gets a peak flux of 0.09 ± 0.08–0.19 ± 0.16LEdd, and the disk is expected to be described as a thin disk that is not truncated, and not too thick at its smallest radius.

As addressed in P. A. Draghis et al. (2023, 2024), BH spin has a series of impacts on our understanding of the outburst behavior for emissions in different energy bands and as well the possible originations of the BH. For the former, the inner disk can reach around the event horizon with high radiation efficiency for an extremely spinning BH. For the latter, it seems that so far we have two populations of BH: the group from observing the BH mergers with the gravitational-wave experiments has slow spin, while those observed from XRB outbursts are always highly spinning (P. A. Draghis et al. 2024). The causality of the different BH populations is not yet known and the discovery of a BH spin in Swift J151857.0–572147 enlarges the BH spin sample, which so far consists of only about 40 BHs held with XRBs, and thus helps to disentangle the models at work in understanding the possible formations of these two BH populations.

Acknowledgments

This work is supported by the National Key R&D Program of China (2021YFA0718500), the National Natural Science Foundation of China under grant Nos. 12333007, U1838202, U2038101, U1938103, 12273030, U1938107, 12027803, and 12173103. This work made use of data and software from the Insight-HXMT mission, a project funded by the China National Space Administration (CNSA) and the Chinese Academy of Sciences (CAS). This work was partially supported by the International Partnership Program of the Chinese Academy of Sciences (grant No. 113111KYSB20190020). This research has made use of software provided by data obtained from the High Energy Astrophysics Science Archive Research Center (HEASARC), provided by NASA's Goddard Space Flight Center. L.D.K. is grateful for the financial support provided by the Sino-German (CSC-DAAD) Postdoc Scholarship Program (91839752).

Footnotes

Please wait… references are loading.
10.3847/2041-8213/ad74ec