Free Access
Issue
A&A
Volume 652, August 2021
Article Number A89
Number of page(s) 12
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/202038531
Published online 13 August 2021

© ESO 2021

1. Introduction

X-ray pulsars (XRPs) are rotating, strongly magnetized neutron stars (NSs) accreting matter from a non-degenerate companion star. The accreting material is funneled by the magnetic field from either an accretion disc or a strong stellar wind to the NS surface, producing pulsed X-ray emission. The X-ray spectra of some of the known XRPs show absorption-like features, the so-called cyclotron resonant scattering features (CRSFs; see e.g. Staubert et al. 2019), which are caused by resonant scattering of photons by the electrons in the accretion stream. Such features, if detected, provide the only direct method to estimate the strength of the NS magnetic field.

4U 1901+03 is a high-mass X-ray binary (HXMB) which was first detected by Uhuru and Vela 5B during a giant outburst in 1970–1971 (Forman et al. 1976; Priedhorsky & Terrell 1984). The source then remained in quiescence until February 2003 when it was detected again by the Rossi X-ray Timing Explorer/All Sky Monitor (RXTE/ASM) during another giant outburst. This detection triggered an extensive observational campaign (Galloway et al. 2003, 2005), which allowed the source to be studied in a broad range of luminosities up to ∼1038 erg s−1 (for the assumed distance of 10 kpc; Galloway et al. 2005), and pulsations were detected with a period of around 2.763 s, clearly establishing the source as an XRP. The pulsar timing also allowed us to determine the orbital parameters of the binary system: the orbital period Porb = 22.58 d, the eccentricity e = 0.036, and the projected semi-major axis ax sin i = 107 lt-s (Galloway et al. 2005).

Although no optical counterpart has been identified, the source was classified as a transient Be/X-ray pulsar based on its orbital parameters and observed outbursts (Galloway et al. 2005). Using the optical and infrared (IR) observations, Reig & Milonaki (2016) were unable to identify obvious candidates for the optical counterpart, largely because of the poor X-ray localisation. Later in 2019, Mereminskiy et al. (2019) and Hemphill et al. (2019) improved the source localisation significantly using the Swift/XRT and Chandra/HRC data. This allowed the authors to propose the IR/optical candidate star J190339.39+031215.6 in the UKIDS-DR6 catalogue as the optical counterpart to 4U 1901+03. The optical spectrum of the star was subsequently obtained and an early spectral type was confirmed, consistent with the Be classification (Strader et al. 2019; McCollum & Laine 2019).

According to the Gaia measurements, the distance to the source is 3.0 1.1 + 2.0 $ ^{+2.0}_{-1.1} $ kpc (Bailer-Jones et al. 2018). However, Strader et al. (2019) argued that the parallax measurement of the optical counterpart in Gaia DR2 is not significant and thus this estimate is not constraining. Furthermore, the PS1 reddening maps (Green et al. 2018), high X-ray absorption, and optical reddening in the direction of the source also point to a larger distance of > 12 kpc (Strader et al. 2019). A similarly large distance of d = 12.4 ± 0.2 kpc was recently derived by Tuo et al. (2020) based on the torque models and the pulse profile evolution during the 2019 outburst.

The pulse profile of 4U 1901+03 is relatively complex, and demonstrates variable structure depending both on the X-ray luminosity and photon energy (Chen et al. 2008; Lei et al. 2009). Based on the NICER and Insight-HXMT data obtained during the 2019 outburst, Ji et al. (2020) discovered short flares in the light curve of the source with significant variability of the pulse profile in the flares and persistent emission. The observed behaviour was interpreted as changes in the beam pattern due to transitions between sub- and supercritical accretion regimes and the onset of an accretion column.

Assuming the spin equilibrium, Galloway et al. (2005) estimated the magnetic field strength of the NS to be B ≈ 0.3 × 1012(d/10 kpc)−6/7 G. Later, James et al. (2011) interpreted the discovered quasi-periodic oscillations (QPOs) at around 0.135 Hz in the context of the beat-frequency model (Alpar & Shaham 1985; Miller et al. 1998). Assuming that the QPO appears at the beat frequency between the pulsar spin frequency and the Keplerian frequency at the magnetospheric radius, the authors obtained the same value of B ≈ 0.3 × 1012 G (assuming distance of 10 kpc). On the other hand, the spin evolution of 4U 1901+03 during the 2019 outburst led Tuo et al. (2020) to argue in favour of ten-times stronger field B ≈ 4.3 × 1012 G, leaving questions as to the NS magnetic field open.

The X-ray spectrum of 4U 1901+03 has a shape similar to those of other accreting pulsars. Galloway et al. (2005) modelled the RXTE data with a thermal Comptonization model. We note that although some residuals around 10 keV were still visible when applying this model (see Fig. 2 in Galloway et al. 2005), the authors did not claim a detection of a cyclotron resonant scattering feature (CRSF) in the spectrum at this energy. Later, Reig & Milonaki (2016) reanalysed the RXTE observations concluding that it is not possible to fit the X-ray spectrum without a Gaussian absorption component that accounts for the absorption-like feature around 10 keV. The authors also found a positive correlation between the energy of the absorption feature and X-ray flux. This makes 4U 1901+03 one of a handful of objects displaying such correlation (Staubert et al. 2019), and allowed Reig & Milonaki (2016) to tentatively suggest that the 10 keV feature is indeed associated with a cyclotron line. This also implies that 4U 1901+03 remained in the subcritical accretion regime during the 2003 outburst. In addition, Coley et al. (2019) recently reported negative residuals from their best-fit model around 30 keV using NuSTAR observations during the declining phase of the 2019 outburst, suggesting that it could be a cyclotron absorption line.

Here we report a detailed broad-band spectral analysis of 4U 1901+03 in a wide range of X-ray luminosities using independent high-quality data obtained by the NuSTAR and Insight-HXMT observatories during the 2019 outburst. Specifically, we aim to clarify the origin of both the 10- and 30 keV absorption features.

2. Observations and data reduction

The 2019 outburst of the source has been extensively monitored by several instruments. Here we focus on the analysis of data obtained by the NuSTAR and Insight-HXMT observatories both of which cover a broad range of energies essential for spectral analysis. A summary of the observations is provided in Table 1 and Fig. 1. The latter shows the Swift/BAT long-term light curve with the dates of the Insight-HXMT and NuSTAR observations indicated. Details of the data reduction for both instruments are presented below.

thumbnail Fig. 1.

15–50 keV X-ray light curve obtained by Swift/BAT during the source outburst. The orange and blue vertical solid lines show the dates of HXMT and NuSTAR observations, respectively.

Table 1.

Observation log of 4U 1901+03 during its 2019 outburst with Insight-HXMT and NuSTAR.

2.1. Insight-HXMT observations

Hard X-ray Modulation Telescope (HXMT; also known as Insight-HXMT), is China’s first X-ray telescope and was launched on 2017 June 15 (Li 2007; Zhang et al. 2020). The telescope has a large effective area in hard X-rays, and covers a broad energy band 1–250 keV. The three slat-collimated instruments onboard are the low-energy X-ray telescope (LE), the medium-energy X-ray telescope (ME), and the high-energy X-ray telescope (HE) operating in the energy ranges 1–15, 5–30, and 20–250 keV, respectively (Zhang et al. 2020; Liu et al. 2020; Cao et al. 2020; Chen et al. 2020).

During the declining phase of the outburst, 4U 1901+03 was observed 20 times by Insight-HXMT over the period from MJD 58563 to MJD 58616. This long duration monitoring provided us with a good opportunity to study the spectral evolution of the source at different luminosities. We note that the background starts to dominate at energies above ∼50 keV even close to the peak of the outburst, and so the HE data unfortunately provide no useful information. Therefore, for the present study we decided to use only LE and ME data covering the energy range between 1 and 30 keV. In order to reduce the HXMT data and extract the source spectra, we performed the standard data-reduction procedure explained in the Insight-HXMT official user guide1 using the software HXMTSOFT v2.01 with a CALDB v2.01 (Li et al. 2020; Liao et al. 2020a,b; Guo et al. 2020).

2.2. NuSTAR observations

The Nuclear Spectroscopic Telescope Array (NuSTAR) is the first focusing hard X-ray observatory and is equipped with two identical co-aligned X-ray telescopes, Focal Plane Module A and B (FPMA and FPMB; Harrison et al. 2013). The instruments cover a wide energy range 3–79 keV and provide an imaging resolution of 18″ (full width at half maximum, FWHM) and a spectral energy resolution of 400 eV (FWHM) at 10 keV. 4U 1901+03 was observed by NuSTAR four times during its recent outburst in a four-month period from February to May 2019 (see Table 1 for more details). The first observation (ObsID 90501305001) was performed at the rising phase of the outburst while the others were done at the declining phase (see Fig. 1). The uniform distribution of the observations over the outburst provides us with a good opportunity to study the evolution of the source properties as a function of luminosity and compare these results to the ones obtained by Insight-HXMT. From now on, we refer to the four NuSTAR observations as NuObs1, NuObs2, NuObs3, and NuObs4, named in the chronological order.

In order to reduce the raw data and extract the source light curves and spectra, we performed the standard data reduction procedure following the NuSTAR official user guides2. We used the NuSTAR Data Analysis Software NUSTARDAS v2.0.0 with a CALDB version 20201130. A source-centred circular region with a radius of 120″ was chosen to extract the source photons from both FPMA and FPMB. The background was extracted similarly from a source-free circular region of the same size. The spectra collected from NuSTAR/FPMA and FMPB were simultaneously fitted in XSPEC in the energy range 4–79 keV. In addition, some of the Insight-HXMT data sets are performed quasi-simultaneously with two NuSTAR observations: Insight-HXMT P0211006008 with NuObs3 and Insight-HXMT P0211006022 with NuObs4. The joint temporal and spectral studies on Her X-1 using data from Insight-HXMT and NuSTAR showed that the two observatories revealed highly consistent results (Xiao et al. 2019), indicating no calibration problems. Therefore, to find the spectral model that best describes the broadband continuum spectrum of the source, we define four data sets as NuObs1 (#1), NuObs2 (#2), NuObs3+P0211006008 (#3), and NuObs4+P0211006022 (#4) to be used for the detailed spectroscopy.

In order to use χ2 statistics, the Insight-HXMT and NuSTAR spectra were each grouped to have at least 30 and 25 counts in each energy bin, respectively, for phased-average and phased-resolved spectral analyses. For all Insight-HXMT and NuSTAR observations, we used HEASOFT 6.283 and XSPEC 12.11.0m4 to perform the spectral analysis.

The pulse phase-resolved spectral analysis also requires determination of the timing properties of the source. To obtain the pulse profile of each observation, we extracted the light curves from the same regions and then applied the barycentric correction to them using the standard BARYCORR task (see Sect. 3.2).

3. Results

3.1. Pulse phase-averaged spectroscopy

Similarly to other accreting XRPs, the spectrum of 4U 1901+03 has a smooth power-law-like shape with a cutoff at high energies (see e.g., Filippova et al. 2005). We attempted to use different physical and phenomenological continuum models such as: a power law with a high-energy cutoff model CUTOFFPL described by the equation N(E) = KE−αexp(−E/β); a Comptonization model COMPTT (Titarchuk 1994); a power law with a high-energy cutoff model PO × HIGHECUT, with the latter defined by the equation M(E) = exp[(Ec − E)/Ef] for (E ≥ Ec) and M(E) = 1 for (E ≤ Ec); a NPEX model, representing the sum of power laws of negative and positive slopes given by equation N(E) = (A1Eα1 + A2E+α2) exp(−E/kT) (Mihara et al. 1998); and a power-law model with the Fermi-Dirac cutoff FDCUT described by equation N(E) = APLE−Γ/[exp((E − Ecut)/Efold)+1] (Tanaka et al. 1986). We also included a Gaussian emission component GAU corresponding to the fluorescent iron line at 6.4 keV and modified the continuum by the photoelectric absorption to account for interstellar and intrinsic absorption TBABS (Wilms et al. 2000) with default abundances from Anders & Grevesse (1989) and atomic cross-sections adopted from Verner et al. (1996). A cross-normalisation constant was also included to account for minor differences in the absolute calibration of the effective areas of the individual instruments. We note that none of the mentioned models used in this analysis give a physical description of the X-ray emission from XRPs as the spectral formation is much more complex in these systems. In recent years, several physically motivated models have been proposed by different authors (see e.g., Becker & Wolff 2007; Farinelli et al. 2012; West et al. 2017; Gornostaev 2021); however, applying such models to fit the spectrum of 4U 1901+03is beyond the scope of this study.

Application of the models listed above lead to fits of unacceptable quality. For example, for the data set #3 we obtain χ2 values of 6781 (for 2957 d.o.f.), 3927 (2956), 5808 (2955), 3720 (2954), and 4425 (2956), respectively, for CUTOFFPL, COMPTT, PO × HIGHECUT, NPEX, and PO × FDCUT. The fits produced residuals around 10 and 30 keV (see Fig. 2-left). Therefore, the source spectrum is more complex than predicted by any of the above-mentioned models. A detailed description of the fitting results for individual models is provided below.

thumbnail Fig. 2.

Left panels: X-ray spectra of 4U 1901+03 obtained by Insight-HXMT (LE and ME detectors in purple and pink for third data set and orange and yellow for fourth data set, respectively) and NuSTAR (from the first to the last chronologically in black, red, green, and blue, respectively. Only FPMA data are shown here for clarity). During the outburst together with the composite model TBABS × (COMPTT + GAUSSIAN) (solid grey curves) shown in panel a. The dashed curves with the same colours show the iron line component. Residuals for the four data sets using the corresponding composite model are shown in units of standard deviations (panels b–e). Right panels: same X-ray spectra as in panel a fitted with the same composite model but with a GABS component added to account for the absorption feature at 10 keV (panel a′). Residuals of the same spectra using the corresponding composite model are shown in panels b′–e′.

3.1.1. Continuum model CUTOFFPL

To be able to make a comparison between our results and previous findings from the literature, in the first step we adopted the model of a power law with a high-energy exponential cutoff (CUTOFFPL in XSPEC) with the additional Gaussian absorption component GABS to fit the absorption structure around 10 keV. Such a model was used by Reig & Milonaki (2016) to fit the RXTE spectra of the source. The composite model TBABS × (GAU+CUTOFFPL × GABS) revealed acceptable fits with χ2(d.o.f.) values of 2102 (1668), 1720 (1568), 3081 (2954), and 3522 (3090), respectively, for our four data sets. The fit parameters and the corresponding uncertainties are given in Table 2. All the uncertainties in the spectral analysis are given at 68.3% (1σ) confidence level everywhere in this paper. We also applied this model to all 20 Insight-HXMT/LE+ME spectra to study the evolution of spectral parameters as a function of X-ray flux in the energy range 1–30 keV. Of particular interest is the dependence of the putative 10 keV absorption feature centroid energy on flux. However, when the results from two telescopes –where all parameters can be fitted independently– are compared, we find that there are inconsistencies between the values of the parameters Eabs1, Γ, NH, and Ecut most probably caused by the different composition of the data sets, that is, a continuous broadband spectrum for NuSTAR in contrast to two joint bands in Insight-HXMT and different energy bands covered by the instruments. As a result, the cross-normalisation constant used in the Insight-HXMT data adds one additional free parameter.

Table 2.

Best-fit parameters for the model TBABS × (GAU+CUTOFFPL × GABS) for the four data sets of pulse-averaged spectra during the 2019 outburst.

In order to obtain more consistent results, we used two combined spectra (NuSTAR + Insight-HXMT), corresponding to data set #3 and data set #4, and found the values of parameters NH and Ecut, which cannot be constrained by the fit using the Insight-HXMT data alone. Therefore, for the NuObs1 and NuObs2 and the remaining 18 Insight-HXMT observations, we fixed the NH and Ecut to values of 2.4 × 1022 cm−2 and 6.66 keV, respectively, which were obtained by averaging the corresponding best-fit values for data sets #3 and #4 to avoid this inconsistency. As shown in Fig. 3, the energy of the 10 keV feature appears to be correlated with X-ray flux, which is in full agreement with the RXTE results (Reig & Milonaki 2016). The photon index also shows some evolution, decreasing with increasing flux, similarly to the RXTE data.

thumbnail Fig. 3.

Evolution of the 10 keV absorption feature energy Eabs1(panel a), its width σabs1(panel b), the photon index Γ(panel c), and the reduced χ2 of the fits (panel d) as a function of the 1–30 keV flux obtained from the Insight-HXMT data (black points). Red and blue dots show the same parameters derived from two single NuSTAR observations and two combined observations (NuSTAR+Insight-HXMT), respectively. We used the model TBABS × (GAU+CUTOFFPL × GABS) here with the NH and Ecut fixed at 2.4 × 1022 cm−2 and 6.66 keV, respectively (see the text).

As we see from Fig. 3, the 10 keV absorption feature centroid energy correlates with the flux from the source. At the same time, the width of this feature demonstrates the opposite behaviour. In the lowest available data set (NuObs4+P0211006022), the line energy and width reach 8.3 and 5.4 keV, respectively. Similar results were obtained from the RXTE observations of the source in which the 10 keV absorption feature centroid energy (width) dropped (rised) remarkably from E ∼ 11.8 (σ∼ 2.5) keV at flux ∼8 × 10−9 erg s−1 cm−2 to E ∼ 6.3 (σ ∼ 5) keV at flux ∼6 × 10−10 erg s−1 cm−2 (see Fig. 7 in Reig & Milonaki 2016). Obviously, in such a situation, the absorption component cannot be considered to be a narrow feature but rather acts as a continuum spectral component. To illustrate the problem, in Fig. 4 we show the 4U 1901+03 spectrum from NuObs4 with the 10 keV feature depth (Dabs1) set to zero. This casts very strong doubt on the interpretation of this feature as a CRSF and/or the use of the continuum model CUTOFFPL proposed by Reig & Milonaki (2016). Therefore, for further analysis we decided to use the broadly accepted Comptonization model COMPTT instead of CUTOFFPL to fit the source spectra.

thumbnail Fig. 4.

Panel a: source spectrum obtained during NuObs4 together with the model TBABS × (GAU+CUTOFFPL × GABS) where the absorption line depth was set to zero (solid line). Panel b: corresponding residuals.

3.1.2. Continuum model COMPTT

The composite model TBABS×(GAU+COMPTT×GABS) also revealed acceptable fits with χ2(d.o.f.) values of 2178 (1668), 1756 (1568), 3167 (2952), and 3579 (3089) for four of our data sets. The corresponding fitted spectra are shown in Fig. 2a′–e′). However, we note that the obtained description of the spectra is still not ideal, with some negative residuals around 30 keV clearly seen in all data sets (especially in data set #3; see Fig. 2d′). This points towards an additional absorption feature, as was already suggested by Coley et al. (2019). Use of the continuum models such as CUTOFFPL, NPEX, and PO × HIGHECUT did not help to get rid of this absorption-like feature (see Fig. 5). The fits revealed χ2/d.o.f. of 1731/1527, 1689/1524, and 1933/1526, respectively, for NuObs3, that is, they have not resulted in a significantly improved description of the spectrum. Therefore, it can be inferred that the residuals are not caused by the use of an inappropriate continuum model. Our best-fit model was then modified by the inclusion of a second Gaussian absorption component GABS.

thumbnail Fig. 5.

Panel a: X-ray spectrum of 4U 1901+03 taken from NuObs3 together with the composite model TBABS × (GAU + CONTINUUM × GABS) (solid lines) where CONTINUUM stands for CUTOFFPL, NPEX and PO × HIGHECUT and shown by blue, red and green curves, respectively. The corresponding residuals are shown in panels b, c and d. Only FPMA data are shown for the sake of clarity.

Accordingly, the following spectral analysis was done using model TBABS×(GAU+COMPTT×GABS×GABS) which gave better χ2 (d.o.f.) of 2085 (1665), 1615 (1564), 2963 (2949), and 3417 (3087), respectively, for the four data sets in comparison to the one with no absorption model at 30 keV. The corresponding best-fit spectra and residuals are shown in Fig. 6 and the best-fit parameters are given in Table 3. Unfortunately, the presence of the 30 keV absorption feature cannot be either confirmed or ruled out based on the Insight-HXMT data, because it appears exactly at the edge of the working energy band of the ME (see, e.g. Fig. 2). We note that for our four data sets the neutral column density NH varies in the range (0.3–1.4) × 1022 cm−2 which is in good agreement with the Galactic mean value of 1.2 × 1022 cm−2 in this direction (Willingale et al. 2013). However, it is lower than the value NH ∼ 3.5 × 1022 cm−2 previously reported by Reig & Milonaki (2016), who used the CUTOFFPL model (without the second GABS) to fit the spectra obtained during the 2003 outburst.

thumbnail Fig. 6.

Same as Fig. 2, but only for the NuSTAR observations together with the best-fit model (TBABS × (GAU+COMPTT × GABS × GABS)). Lighter colours represent the data obtained from the NuSTAR/FPMB module.

Table 3.

Best-fit parameters for the model TBABS × (GAU+COMPTT × GABS × GABS) for the four data sets of pulse-averaged spectra during the 2019 outburst.

Finally, in order to estimate the significance of the additional Gaussian absorption component at 30 keV in the NuSTAR observations, we used the SIMFTEST5 script in XSPEC. For this, we simulated 104 spectra and estimated the probability of a false detection. As a result, for all NuSTAR observations fitted by the best-fit model, we obtained the upper limit probability of 10−4 with the null hypothesis that the absorption component is caused by statistical fluctuations, suggesting the 30 keV feature has a significance of at least ∼4σ. We also tried to roughly estimate the F-test probability using FTEST6 task in XSPEC although it may not be a correct method to calculate the significance level of a line-like feature (Protassov et al. 2002). The F-test probability for all four of our data sets is 4.4 × 10−21, 9.6 × 10−27, 4.1 × 10−42, and 5.7 × 10−32, respectively.

In order to investigate the evolution of the COMPTT model on the flux in more detail, we applied it to all 20 Insight-HXMT/LE+ME spectra, excluding the 30 keV absorption feature from the total model. Interestingly, no absorption feature at 10 keV was required to fit the data in contrast to the CUTOFFPL model. However, we let this component remain in the model in order to be able to compare the spectral evolution using both Insight-HXMT and NuSTAR observations. To do this, we fixed the parameters of the 10 keV Gaussian absorption component to the best-fit values obtained from the data set #3 (NuObs3+P0211006008; see Table 3) representing the average luminosity state of the source.

Figure 7 indicates the evolution of photoelectric absorption (NH), plasma temperature (kT), input photon temperature (T0), and plasma optical depth (τp) as a function of flux. As can be seen, the NH value for both Insight-HXMT and NuSTAR observations shows a negative correlation with the X-ray flux varying in the range (0.3–1.5) × 1022 cm−2. The kT and τp values do not show clear correlation with the flux. At the same time, the T0 shows a very prominent positive correlation with flux.

thumbnail Fig. 7.

Evolution of the interstellar absorption NH (panel a), the plasma temperature kT (panel b), the input photon temperature T0 (panel c), the plasma optical depth τp (panel d), and the reduced χ2 of the fits (panel e) as a function of the 1–30 keV flux obtained from the Insight-HXMT spectra using the TBABS × (GAU+CUTOFFPL × GABS) model. Red and blue dots show the same parameters derived from two single NuSTAR observations and two combined observations (NuSTAR+Insight-HXMT), respectively.

Although we used the same composite model to fit the Insight-HXMT and NuSTAR observations, there is some discrepancy in the results obtained with these two observatories, as seen in the figure. Most probably, the observed differences are caused by the different composition of the data sets from different instruments, as discussed above.

3.1.3. Two-component continuum model

Another possible explanation for the observed absorption-like feature around 10 keV is the two-component spectrum. Originally, this idea was proposed for XRP X Persei to explain a similar feature around 20–30 keV (Di Salvo et al. 1998; Doroshenko et al. 2012), which alternatively was also interpreted as a broad CRSF (Coburn et al. 2001; Lutovinov et al. 2012). More recently, using observations of two XRPs, A 0535+262 and GX 304–1, in their low-luminosity states, Tsygankov et al. (2019a,b) showed that the appearance of the double-hump spectrum is not caused by the cyclotron absorption lines, but may be attributed to the emission of cyclotron photons in the NS atmosphere caused by collisional excitation of electrons to upper Landau levels and consequent Comptonization by hot electron gas (Mushtukov et al. 2021). We emphasise that the two-component spectrum appeared at comparatively low luminosities (∼1035 erg s−1) for the sources mentioned above, and so the validity of this scenario is questionable for 4U 1901+03. However, spectra of 1A 0535+262 at higher luminosities (∼1036 erg s−1) were found to be well described with the same model (Tsygankov et al. 2019a) even if the presence of both components was not obvious. We note also that a similar shaped spectrum can also be produced by other physical mechanisms, for instance by scattering of the emission from the accretion column of the NS surface (Ferrigno et al. 2011). However, it is important to note that such a complex continuum shape may easily be mistakenly interpreted as a cyclotron absorption feature in the standard XRP spectrum.

In order to verify the applicability of this scenario to 4U 1901+03, we tried to describe the NuSTAR spectra using a model consisting of the two separate components. First, we adopted the same two-component model consisting of two Comptonization models, COMPTT, as previously used to describe the double-hump spectra of XRPs A 0535+262 and GX 304–1 in their low-luminosity states (Tsygankov et al. 2019a,b). We then modified it by adding photoelectric absorption (TBABS) and an iron line emission component (GAU). For data set #3, this model resulted in χ2 = 3306 for 2951 d.o.f., but left residuals around 30 keV. Therefore, we added a Gaussian absorption component (GABS) to account for the 30 keV absorption feature which improved the fit significantly to a χ2 = 2970 for 2948 d.o.f.

The unfolded spectrum obtained from data set #3 fitted with a model consisting of two COMPTT components is shown in Fig. 8 and the corresponding fit parameters are listed in Table 4. As shown in Fig. 8b, the two-component composite model fits the source spectrum without leaving any residuals associated with the absorption feature at 10 keV. Application of this model to the remaining three data sets gives a good fit for NuObs4 (χ2/d.o.f. = 3901/3087), but leaves some insignificant residuals around 10 keV for the brighter observations NuObs1 (χ2 (d.o.f.) = 2089 (1663)) and NuObs2 (χ2/d.o.f. = 1674/1563). The non-ideal quality of the fit can be explained by the fact that the simple COMPTT model is not physically motivated for either component in this kind of spectra.

thumbnail Fig. 8.

Panel a: X-ray spectrum of 4U 1901+03 collected during NuObs3 together with the two-component models TBABS × (GAU + COMPTT + COMPTT × GABS) and TBABS × (GAU + PO × HIGHECUT + PO × HIGHECUT × GABS) where the corresponding model components are shown by the dashed curves in red and blue, respectively. Panel b: residuals for the former composite model (COMPTT). Panel c: residuals for the latter composite model (PO×HIGHECUT). The color code is same as for Fig. 2.

Table 4.

Best-fit parameters from the fit to the NuObs3 spectrum with the two-component models.

To verify this assumption, we adopted the same model that Di Salvo et al. (1998) used to fit the X Persei spectrum consisting of two (PO × HIGHECUT) components. The modelled spectrum and the corresponding fit parameters for data set #3 are shown in Fig. 8 and Table 4. We emphasise that for this composite model, the 30-keV absorption feature also needs to be taken into consideration. We note that the continuum model HIGHECUT introduces artificial absorption-line like structures because of the break in the derivative at the cutoff energy. Therefore, in order to get rid of such artificial residuals, we introduced a Gaussian absorption component (GABS) to the model (see Coburn et al. 2002) in which the centroid energy is fixed to the position of the cutoff energy determined by HIGHECUT. The width of the component is also set to 10% of the centroid energy. This model allowed us to fit all four data sets without leaving any residuals associated with the absorption feature at 10 keV for χ2(d.o.f.) of 2059 (1664), 1635 (1564), 2949 (2949), and 3492 (3086), respectively.

3.2. Phase-resolved spectroscopy

Phase-resolved spectroscopy may serve as a unique source of information about the spatial properties of the emitting region of the NS. Using the barycentrically corrected X-ray light curves, we obtained the spin periods for four NuSTAR observations following the standard epoch folding technique (Leahy et al. 1983) using the EFSEARCH routines from the FTOOL package (Table 5). The epoch Tmax of the maximal flux in the pulse profile used as a zero phase for the phase-resolved spectral analysis is presented in the same table for each observation. The spin periods and the corresponding uncertainties were obtained from a simulation of 103 light curves using the count rates and corresponding errors from the original light curves (see e.g., Boldin et al. 2013, and references therein). We then obtained the pulse profiles of all four NuSTAR observations in the full energy band (Fig. 9). We applied the same procedure to the Insight-HXMT P021106008 observation in order to be able to perform a simultaneous phase-resolved analysis using NuObs3+P021106008.

thumbnail Fig. 9.

Pulse profiles of 4U 1901+03 from four NuSTAR observations in black (NuObs1), red (NuObs2), green (NuObs3), and blue (NuObs4) in the 3–79 keV energy band. The purple pulse profile was obtained from Insight-HXMT P021106008 in the energy band 1–30 keV. Black dashed lines show the phase segments we used to extract the phase-resolved spectra.

Table 5.

Timing properties of 4U 1901+03.

A detailed study of the pulse profile evolution of 4U 1901+03 during the 2019 outburst was performed by Ji et al. (2020) and Tuo et al. (2020). We then performed a phase-resolved spectroscopy for 4U 1901+03 using all four data sets. We find that the dependence of the model parameters on the spin phase is very similar for all the observations and does not depend on the source luminosity. We therefore discuss the results of the phase-resolved analysis mostly for data set #3 (NuObs3+P021106008) where the 30 keV feature is most prominent. Taking into the account the available counting statistics, we considered ten equally spaced phase bins for this analysis as shown in Fig. 9. Each spectrum was fitted with our best-fit model. However, we had to fix the width of the 30 keV absorption feature σabs2 to its phase-averaged values 5.8 keV as it was not constrained in the fits due to lower count statistics at high energies.

The NH remains constant around 1.2 × 1022 cm−2, indicating no dependence on the pulse phase (see Fig. 10). The soft photon temperature T0 shows insignificant variations between 1 and 1.2 keV. The plasma temperature kT and the optical depth τp are clearly dependent on the pulse phases, showing negative and positive correlations with the X-ray flux, respectively. The position Eabs1 and width σabs1 of the first absorption feature remain almost constant around 10 keV and 2–3 keV, respectively. In contrast to that, parameters of the 30-keV absorption feature strongly depend on the pulse phase. In order to illustrate how these variations depend on the source luminosity, we show in Fig. 11 the results of our spectral analysis for all NuSTAR data. As seen from the figure, the position of the 30-keV feature Eabs2 shows a positive correlation with the X-ray flux. At the same time, the optical depth of the line demonstrates different behaviours at different luminosities. In the first two bright observations, the optical depth is almost correlated with X-ray flux. However, in the latter two observations, it shows a negative correlation with flux in the main maximum and a positive correlation in the main minimum.

thumbnail Fig. 10.

Variations of the best-fit parameters of model TBABS × (GAU+COMPTT × GABS × GABS) as a function of pulse phase obtained from the NuObs3+P0211006008. Panels a–i: observed X-ray flux in the energy band 1–79 keV in units of 10−9 erg s−1 cm−2, first absorption feature energy in keV, width in keV, and its optical depth, neutral hydrogen column density NH in units of 1022 cm−2, input soft photon temperate in keV, plasma temperature in keV, plasma optical depth, and the reduced χ2, respectively. The values of the red points are fixed at the best-fit values and the blue dot represents the 3σ upper limit of the first absorption feature optical depth. The second absorption feature parameters are shown in Fig. 11.

thumbnail Fig. 11.

Second absorption feature centroid energy (in keV) and its optical depth variations as a function of pulse phase obtained from NuObs1, NuObs2, NuObs3+P0211006008, and NuObs4. The pulse profile of each observation is shown in grey. For each observation, the width of the second absorption feature σabs2 was fixed to its phase-averaged value.

4. Discussion and conclusions

Using the broadband NuSTAR data, we find that both the 10 and 30 keV features appear significant when applying conventional continuum models. However, in contrast to previous studies (Reig & Milonaki 2016; Beri et al. 2021), we favour the adoption of the Comptonization model COMPTT to describe the broadband X-ray spectrum of the source instead of the CUTOFFPL continuum model as we find this latter to be affected by systematic bias as explained below.

Particularly, by using CUTOFFPL we were able to confirm that the energy Eabs1 of the 10 keV feature demonstrates a positive correlation with the X-ray flux as was previously reported for the RXTE data (see the Fig. 7 in Reig & Milonaki 2016). This implies that, assuming the 10 keV feature is a CRSF, 4U 1901+03 remains in a subcritical accretion regime during the outburst with the luminosity below the critical one (Basko & Sunyaev 1976; Becker et al. 2012; Mushtukov et al. 2015) which was also discussed by Reig & Milonaki (2016). Assuming critical luminosity to be around 1037 erg s−1, the distance to the source cannot be more than ∼3–4 kpc, which is consistent with the Gaia measurements (Bailer-Jones et al. 2018). Taking this distance into account, the X-ray luminosities observed by NuSTAR for the source lie in the range (2.8–8.3) × 1036 erg s−1. Therefore, interpretation of the 10-keV feature as a CRSF in the subcritical pulsar requires a small distance to the source of below ∼4 kpc.

However, a more favourable large distance above 12 kpc (Strader et al. 2019; Reig & Milonaki 2016; Galloway et al. 2005; Tuo et al. 2020) results in X-ray luminosities in the range (4.5–13.3) × 1037 erg s−1 for the four NuSTAR observations, which is above the critical luminosity for a XRP with a typical magnetic field. This is at odds with the idea that the 10-keV feature has a magnetic origin because, in the case of bright XRPs, a negative CRSF energy–luminosity correlation is expected (Mihara et al. 1998; Tsygankov et al. 2006; Becker et al. 2012). Moreover, Mushtukov et al. (2015) demonstrated that for a XRP with magnetic field around 1012 G, the critical luminosity may be as low as a few ×1036 erg s−1. In this case, according to the NuSTAR observations, 4U 1901+03 should be in the supercritical accretion regime during the outburst even at a distance of 3 kpc. Therefore, the observed properties of the 10-keV feature (i.e., positive correlation of Eabs1 with X-ray luminosity) is hard to explain in terms of the cyclotron line for any available estimates of the distance. It is worth noting that the NuSTAR spectra may show an instrumental feature around 10 keV as seen in some cases. For instance, in Be-transient XRP KS 1947+300, the 12.5 keV absorption feature was long debated until recently when Doroshenko et al. (2020) drew attention to the instrumental effects causing the appearance of the feature in the NuSTAR spectra. The instrumental systematic effect in NuSTAR observations results from the tungsten L-edge of the telescope optics coupled with small uncertainties in the energy scale offset (Madsen et al. 2015). However, the existence of the absorption like feature around 10 keV is confirmed by the RXTE/PCA (Reig & Milonaki 2016), ruling out an instrumental origin.

Moreover, we noticed that the combination of the continuum model CUTOFFPL with a Gaussian absorption component GABS tends to modify the continuum itself rather than producing a narrow absorption feature at 10 keV. As a result, the width of the feature at low flux becomes comparable to its centroid energy (see Figs. 3 and 4). This behaviour was not discussed previously in the literature and casts serious doubt on the CRSF nature of the feature.

Taking into account all the problems mentioned above, we applied COMPTT to the continuum, which resulted in a significantly different description of the spectrum. Particularly, we did not detect an absorption feature at 10 keV in any of the 20 Insight-HXMT observations. This is in line with the results of Reig & Milonaki (2016) who also admitted that the Comptonization models cause the 10-keV absorption feature to appear insignificant in the RXTE/PCA data. However, the feature was required in all NuSTAR observations. Unfortunately, the Comptonization model also showed the same problem as CUTOFFPL by revealing a large width σabs1 = 4.4 keV for the centroid energy Eabs1 = 9 keV in the NuObs4 which was obtained when the source showed low luminosity. Interestingly, the centroid energy Eabs1 of the 10 keV feature remains almost constant at around 10.6 keV during the outburst as seen in the first three NuSTAR observations. This clearly demonstrates that some properties of the 10 keV feature depend on the continuum model.

We demonstrate that the position of the 10 keV feature is stable along the pulse phase. On the other hand, phase-resolved spectroscopy reveals a strong dependence of the 30 keV feature on the pulse phases in all luminosity states of the source, which is typical behaviour for CRSFs. It is worth mentioning here that the centroid energies of the two absorption features do not appear to be harmonically related. Although not a typical behaviour for CRSFs, this can be explained by assuming that both features are associated with fundamental cyclotron lines originating at different altitudes above the NS. Particularly, considering the strong dependence of the magnetic field strength on the height of the emission region (B ∝ h−3), the 10 keV absorption line may originate at ∼4.5 km (for a NS radius of R* = 10 km) above the NS surface. This scenario may also be combined with the so-called offset magnetic dipole resulting in the presence of two independent fundamental cyclotron lines (see e.g., Iyer et al. 2015).

Finally, we suggest another possible explanation for the 10 keV absorption feature. We demonstrate that the two-component continuum provides an acceptable description of the broadband spectra of 4U 1901+03 with no need to introduce an absorption 10 keV feature into the model, despite the two-component COMPTT model leaving some residuals around 10 keV in three out of four data sets. The purpose of using COMPTT in this context as a phenomenological model was to be able to compare the results with the previous observational studies. However, physically motivated theoretical models for the emission from highly magnetised NS atmospheres at low mass-accretion rates have appeared recently and demonstrate very good consistency with the observational data on other low-luminosity XRPs (Mushtukov et al. 2021; Sokolova-Lapa et al. 2021).

It should be noted that transformation of a typical cutoff power-law-shaped spectrum observed in XRPs at high luminosities to a double-component spectrum happens at low luminosities of ∼1034−36 erg s−1 (Doroshenko et al. 2012; Tsygankov et al. 2019a,b; Lutovinov et al. 2021). Therefore, this explanation may also require a small distance to the source. At the same time, evidence for a two-component continuum starts to appear in some bright XRPs as well (Doroshenko et al. 2020). We emphasise that even in the case of the two-component continuum model, the 30 keV feature still needs to be taken into account.

It is worth noting that a magnetic origin of the 30 keV feature in the source spectrum was supported by the investigation of the pulse profile shape as a function of energy performed by Beri et al. (2021). The authors discovered an abrupt change in the pulse profiles around this energy that is a typical behaviour around CRSF energy (see e.g., Tsygankov et al. 2006). At the same time, no significant deviation from the general trend around 10 keV was found. This fact also supports speculation that the 10 keV feature may have a non-magnetic origin.


Acknowledgments

This work was supported by the grant 14.W03.31.0021 of the Ministry of Science and Higher Education of the Russian Federation. We also acknowledge the support from the Finnish Cultural Foundation project 00200764 (AN), the Academy of Finland travel grants 317552, 322779, 324550, 331951, and 333112 (AN, SST, JP), the German Academic Exchange Service (DAAD) travel grants 57405000 and 57525212 (LJ, VD), the National Key R&D Program of China (2016YFA0400800), and the National Natural Science Foundation of China grants U1838201, U1838202, U1938103 and 11733009 (LJ, YL, SZ, SZ, LJ). This work made use of data from the Insight-HXMT mission, a project funded by China National Space Administration (CNSA) and the Chinese Academy of Sciences (CAS).

References

  1. Alpar, M. A., & Shaham, J. 1985, Nature, 316, 239 [NASA ADS] [CrossRef] [Google Scholar]
  2. Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta., 53, 197 [NASA ADS] [CrossRef] [Google Scholar]
  3. Bailer-Jones, C. A. L., Rybizki, J., Fouesneau, M., Mantelet, G., & Andrae, R. 2018, AJ, 156, 58 [NASA ADS] [CrossRef] [Google Scholar]
  4. Basko, M. M., & Sunyaev, R. A. 1976, MNRAS, 175, 395 [NASA ADS] [CrossRef] [Google Scholar]
  5. Becker, P. A., & Wolff, M. T. 2007, ApJ, 654, 435 [NASA ADS] [CrossRef] [Google Scholar]
  6. Becker, P. A., Klochkov, D., Schönherr, G., et al. 2012, A&A, 544, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Beri, A., Girdhar, T., Iyer, N. K., & Maitra, C. 2021, MNRAS, 500, 1350 [CrossRef] [Google Scholar]
  8. Boldin, P. A., Tsygankov, S. S., & Lutovinov, A. A. 2013, Astron. Lett., 39, 375 [NASA ADS] [CrossRef] [Google Scholar]
  9. Cao, X., Jiang, W., Meng, B., et al. 2020, Sci. China Phys. Mech. Astron., 63, 249504 [NASA ADS] [CrossRef] [Google Scholar]
  10. Chen, W., Qu, J.-L., Zhang, S., Zhang, F., & Zhang, G.-B. 2008, Chin. Astron. Astrophys., 32, 241 [NASA ADS] [CrossRef] [Google Scholar]
  11. Chen, Y., Cui, W., Li, W., et al. 2020, Sci. China Phys. Mech. Astron., 63, 249505 [NASA ADS] [CrossRef] [Google Scholar]
  12. Coburn, W., Heindl, W. A., Gruber, D. E., et al. 2001, ApJ, 552, 738 [NASA ADS] [CrossRef] [Google Scholar]
  13. Coburn, W., Heindl, W. A., Rothschild, R. E., et al. 2002, ApJ, 580, 394 [NASA ADS] [CrossRef] [Google Scholar]
  14. Coley, J. B., Fuerst, F., Hemphill, P., et al. 2019, ATel., 12684, 1 [Google Scholar]
  15. Di Salvo, T., Burderi, L., Robba, N. R., & Guainazzi, M. 1998, ApJ, 509, 897 [NASA ADS] [CrossRef] [Google Scholar]
  16. Doroshenko, V., Santangelo, A., Kreykenbohm, I., & Doroshenko, R. 2012, A&A, 540, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Doroshenko, R., Piraino, S., Doroshenko, V., & Santangelo, A. 2020, MNRAS, 493, 3442 [CrossRef] [Google Scholar]
  18. Farinelli, R., Ceccobello, C., Romano, P., & Titarchuk, L. 2012, A&A, 538, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Ferrigno, C., Falanga, M., Bozzo, E., et al. 2011, A&A, 532, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Filippova, E. V., Tsygankov, S. S., Lutovinov, A. A., & Sunyaev, R. A. 2005, Astron. Lett., 31, 729 [NASA ADS] [CrossRef] [Google Scholar]
  21. Forman, W., Jones, C., & Tananbaum, H. 1976, ApJ, 206, L29 [NASA ADS] [CrossRef] [Google Scholar]
  22. Galloway, D., Remillard, R., Morgan, E., & Swank, J. 2003, IAU circ., 8070, 2 [Google Scholar]
  23. Galloway, D. K., Wang, Z., & Morgan, E. H. 2005, ApJ, 635, 1217 [NASA ADS] [CrossRef] [Google Scholar]
  24. Gornostaev, M. I. 2021, MNRAS, 501, 564 [CrossRef] [Google Scholar]
  25. Green, G. M., Schlafly, E. F., Finkbeiner, D., et al. 2018, MNRAS, 478, 651 [NASA ADS] [CrossRef] [Google Scholar]
  26. Guo, C.-C., Liao, J.-Y., Zhang, S., et al. 2020, J. High Energy Astrophys., 27, 44 [CrossRef] [Google Scholar]
  27. Harrison, F. A., Craig, W. W., Christensen, F. E., et al. 2013, ApJ, 770, 103 [Google Scholar]
  28. Hemphill, P., Coley, J., Fuerst, F., et al. 2019, ATel., 12556, 1 [Google Scholar]
  29. Iyer, N., Mukherjee, D., Dewangan, G. C., Bhattacharya, D., & Seetha, S. 2015, MNRAS, 454, 741 [NASA ADS] [CrossRef] [Google Scholar]
  30. James, M., Paul, B., Devasia, J., & Indulekha, K. 2011, MNRAS, 410, 1489 [NASA ADS] [CrossRef] [Google Scholar]
  31. Ji, L., Ducci, L., Santangelo, A., et al. 2020, MNRAS, 493, 5680 [CrossRef] [Google Scholar]
  32. Leahy, D. A., Darbro, W., Elsner, R. F., et al. 1983, ApJ, 266, 160 [NASA ADS] [CrossRef] [Google Scholar]
  33. Lei, Y.-J., Chen, W., Qu, J.-L., et al. 2009, ApJ, 707, 1016 [NASA ADS] [CrossRef] [Google Scholar]
  34. Li, T.-P. 2007, Nucl. Phys. B Proc. Suppl., 166, 131 [NASA ADS] [CrossRef] [Google Scholar]
  35. Li, X., Li, X., Tan, Y., et al. 2020, J. High Energy Astrophys., 27, 64 [CrossRef] [Google Scholar]
  36. Liao, J.-Y., Zhang, S., Lu, X.-F., et al. 2020a, J. High Energy Astrophys., 27, 14 [CrossRef] [Google Scholar]
  37. Liao, J.-Y., Zhang, S., Chen, Y., et al. 2020b, J. High Energy Astrophys., 27, 24 [CrossRef] [Google Scholar]
  38. Liu, C., Zhang, Y., Li, X., et al. 2020, Sci. China Phys. Mech. Astron., 63, 249503 [NASA ADS] [CrossRef] [Google Scholar]
  39. Lutovinov, A., Tsygankov, S., & Chernyakova, M. 2012, MNRAS, 423, 1978 [NASA ADS] [CrossRef] [Google Scholar]
  40. Lutovinov, A., Tsygankov, S., Molkov, S., et al. 2021, ApJ, 912, 17 [Google Scholar]
  41. Madsen, K. K., Harrison, F. A., Markwardt, C. B., et al. 2015, ApJS, 220, 8 [NASA ADS] [CrossRef] [Google Scholar]
  42. McCollum, B., & Laine, S. 2019, ATel., 12560, 1 [Google Scholar]
  43. Mereminskiy, I. A., Lutovinov, A. A., Tsygankov, S. S., Semena, A. N., & Shtykovskiy, A. E. 2019, ATel., 12514, 1 [Google Scholar]
  44. Mihara, T., Makishima, K., & Nagase, F. 1998, Adv. Space Res., 22, 987 [NASA ADS] [CrossRef] [Google Scholar]
  45. Miller, M. C., Lamb, F. K., & Psaltis, D. 1998, ApJ, 508, 791 [NASA ADS] [CrossRef] [Google Scholar]
  46. Mushtukov, A. A., Suleimanov, V. F., Tsygankov, S. S., & Poutanen, J. 2015, MNRAS, 447, 1847 [NASA ADS] [CrossRef] [Google Scholar]
  47. Mushtukov, A. A., Suleimanov, V. F., Tsygankov, S. S., & Portegies Zwart, S. 2021, MNRAS, 503, 5193 [Google Scholar]
  48. Priedhorsky, W. C., & Terrell, J. 1984, ApJ, 280, 661 [NASA ADS] [CrossRef] [Google Scholar]
  49. Protassov, R., van Dyk, D. A., Connors, A., Kashyap, V. L., & Siemiginowska, A. 2002, ApJ, 571, 545 [NASA ADS] [CrossRef] [Google Scholar]
  50. Reig, P., & Milonaki, F. 2016, A&A, 594, A45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Sokolova-Lapa, E., Gornostaev, M., Wilms, J., et al. 2021, A&A, 651, A12 [CrossRef] [EDP Sciences] [Google Scholar]
  52. Staubert, R., Trümper, J., Kendziorra, E., et al. 2019, A&A, 622, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Strader, J., Chomiuk, L., Swihart, S., & Aydi, E. 2019, ATel., 12554, 1 [Google Scholar]
  54. Tanaka, Y. 1986, in Radiation Hydrodynamics in Stars and Compact Objects, eds. D. Mihalas, & K. H. A. Winkler (Berlin: Springer-Verlag), Lect. Notes Phys., 255, 198 [NASA ADS] [CrossRef] [Google Scholar]
  55. Titarchuk, L. 1994, ApJ, 434, 570 [NASA ADS] [CrossRef] [Google Scholar]
  56. Tsygankov, S. S., Lutovinov, A. A., Churazov, E. M., & Sunyaev, R. A. 2006, MNRAS, 371, 19 [NASA ADS] [CrossRef] [Google Scholar]
  57. Tsygankov, S. S., Doroshenko, V., Mushtukov, A. A., et al. 2019a, MNRAS, 487, L30 [CrossRef] [Google Scholar]
  58. Tsygankov, S. S., Rouco Escorial, A., Suleimanov, V. F., et al. 2019b, MNRAS, 483, L144 [NASA ADS] [CrossRef] [Google Scholar]
  59. Tuo, Y. L., Ji, L., Tsygankov, S. S., et al. 2020, J. High Energy Astrophys., 27, 38 [CrossRef] [Google Scholar]
  60. Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 465, 487 [Google Scholar]
  61. West, B. F., Wolfram, K. D., & Becker, P. A. 2017, ApJ, 835, 130 [NASA ADS] [CrossRef] [Google Scholar]
  62. Willingale, R., Starling, R. L. C., Beardmore, A. P., Tanvir, N. R., & O’Brien, P. T. 2013, MNRAS, 431, 394 [NASA ADS] [CrossRef] [Google Scholar]
  63. Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 [Google Scholar]
  64. Xiao, G. C., Ji, L., Staubert, R., et al. 2019, J. High Energy Astrophys., 23, 29 [CrossRef] [Google Scholar]
  65. Zhang, S.-N., Li, T., Lu, F., et al. 2020, Sci. China Phys. Mech. Astron., 63, 249502 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Observation log of 4U 1901+03 during its 2019 outburst with Insight-HXMT and NuSTAR.

Table 2.

Best-fit parameters for the model TBABS × (GAU+CUTOFFPL × GABS) for the four data sets of pulse-averaged spectra during the 2019 outburst.

Table 3.

Best-fit parameters for the model TBABS × (GAU+COMPTT × GABS × GABS) for the four data sets of pulse-averaged spectra during the 2019 outburst.

Table 4.

Best-fit parameters from the fit to the NuObs3 spectrum with the two-component models.

Table 5.

Timing properties of 4U 1901+03.

All Figures

thumbnail Fig. 1.

15–50 keV X-ray light curve obtained by Swift/BAT during the source outburst. The orange and blue vertical solid lines show the dates of HXMT and NuSTAR observations, respectively.

In the text
thumbnail Fig. 2.

Left panels: X-ray spectra of 4U 1901+03 obtained by Insight-HXMT (LE and ME detectors in purple and pink for third data set and orange and yellow for fourth data set, respectively) and NuSTAR (from the first to the last chronologically in black, red, green, and blue, respectively. Only FPMA data are shown here for clarity). During the outburst together with the composite model TBABS × (COMPTT + GAUSSIAN) (solid grey curves) shown in panel a. The dashed curves with the same colours show the iron line component. Residuals for the four data sets using the corresponding composite model are shown in units of standard deviations (panels b–e). Right panels: same X-ray spectra as in panel a fitted with the same composite model but with a GABS component added to account for the absorption feature at 10 keV (panel a′). Residuals of the same spectra using the corresponding composite model are shown in panels b′–e′.

In the text
thumbnail Fig. 3.

Evolution of the 10 keV absorption feature energy Eabs1(panel a), its width σabs1(panel b), the photon index Γ(panel c), and the reduced χ2 of the fits (panel d) as a function of the 1–30 keV flux obtained from the Insight-HXMT data (black points). Red and blue dots show the same parameters derived from two single NuSTAR observations and two combined observations (NuSTAR+Insight-HXMT), respectively. We used the model TBABS × (GAU+CUTOFFPL × GABS) here with the NH and Ecut fixed at 2.4 × 1022 cm−2 and 6.66 keV, respectively (see the text).

In the text
thumbnail Fig. 4.

Panel a: source spectrum obtained during NuObs4 together with the model TBABS × (GAU+CUTOFFPL × GABS) where the absorption line depth was set to zero (solid line). Panel b: corresponding residuals.

In the text
thumbnail Fig. 5.

Panel a: X-ray spectrum of 4U 1901+03 taken from NuObs3 together with the composite model TBABS × (GAU + CONTINUUM × GABS) (solid lines) where CONTINUUM stands for CUTOFFPL, NPEX and PO × HIGHECUT and shown by blue, red and green curves, respectively. The corresponding residuals are shown in panels b, c and d. Only FPMA data are shown for the sake of clarity.

In the text
thumbnail Fig. 6.

Same as Fig. 2, but only for the NuSTAR observations together with the best-fit model (TBABS × (GAU+COMPTT × GABS × GABS)). Lighter colours represent the data obtained from the NuSTAR/FPMB module.

In the text
thumbnail Fig. 7.

Evolution of the interstellar absorption NH (panel a), the plasma temperature kT (panel b), the input photon temperature T0 (panel c), the plasma optical depth τp (panel d), and the reduced χ2 of the fits (panel e) as a function of the 1–30 keV flux obtained from the Insight-HXMT spectra using the TBABS × (GAU+CUTOFFPL × GABS) model. Red and blue dots show the same parameters derived from two single NuSTAR observations and two combined observations (NuSTAR+Insight-HXMT), respectively.

In the text
thumbnail Fig. 8.

Panel a: X-ray spectrum of 4U 1901+03 collected during NuObs3 together with the two-component models TBABS × (GAU + COMPTT + COMPTT × GABS) and TBABS × (GAU + PO × HIGHECUT + PO × HIGHECUT × GABS) where the corresponding model components are shown by the dashed curves in red and blue, respectively. Panel b: residuals for the former composite model (COMPTT). Panel c: residuals for the latter composite model (PO×HIGHECUT). The color code is same as for Fig. 2.

In the text
thumbnail Fig. 9.

Pulse profiles of 4U 1901+03 from four NuSTAR observations in black (NuObs1), red (NuObs2), green (NuObs3), and blue (NuObs4) in the 3–79 keV energy band. The purple pulse profile was obtained from Insight-HXMT P021106008 in the energy band 1–30 keV. Black dashed lines show the phase segments we used to extract the phase-resolved spectra.

In the text
thumbnail Fig. 10.

Variations of the best-fit parameters of model TBABS × (GAU+COMPTT × GABS × GABS) as a function of pulse phase obtained from the NuObs3+P0211006008. Panels a–i: observed X-ray flux in the energy band 1–79 keV in units of 10−9 erg s−1 cm−2, first absorption feature energy in keV, width in keV, and its optical depth, neutral hydrogen column density NH in units of 1022 cm−2, input soft photon temperate in keV, plasma temperature in keV, plasma optical depth, and the reduced χ2, respectively. The values of the red points are fixed at the best-fit values and the blue dot represents the 3σ upper limit of the first absorption feature optical depth. The second absorption feature parameters are shown in Fig. 11.

In the text
thumbnail Fig. 11.

Second absorption feature centroid energy (in keV) and its optical depth variations as a function of pulse phase obtained from NuObs1, NuObs2, NuObs3+P0211006008, and NuObs4. The pulse profile of each observation is shown in grey. For each observation, the width of the second absorption feature σabs2 was fixed to its phase-averaged value.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.