Nuclear Spectroscopic Telescope Array Data Research Articles

Weak secondary cyclotron line in eclipsing high-mass X-ray binary Cen X-3

ABSTRACT We report the time-resolved spectroscopy results from two observations of Centaurus X-3, over one binary orbit with AstroSat and two binary orbits with Nuclear Spectroscopic Telescope Array (NuSTAR). NuSTAR covered two intensity states where the light curve showed transition in count rate from first to second binary orbit by a factor of ∼3. A phenomenological model comprising of partially absorbed power law with smoothed high energy cutoff, cyclotron absorption ∼24 and 6.4 keV iron emission gave good fit for AstroSat observation. NuSTAR spectra required two additional emission components, a broad one ∼5.7 keV and a narrow one ∼6.9 keV. A weak secondary absorption feature at ∼11.6 and ∼ 14.5 keV was seen in the residuals of the spectral fit for AstroSat and NuSTAR data, respectively. The secondary absorption energy showed no correlation with the cutoff energy. Its strength varied within 0.1–0.6 keV with its width ∼1.6 keV. Its energy and optical depth showed linear positive correlation with the fundamental cyclotron line energy and depth, respectively. The cyclotron line energy showed anticorrelation to flux described by a power law with negative index and the secondary absorption also showed similar trend to flux. Depth of secondary absorption was ∼45 per cent and centroid energy was ∼54 per cent of fundamental. Depth and energy ratio of secondary to fundamental lied within 2σ deviation from 0.5. We suggest this secondary absorption to be a redshifted dipolar cyclotron resonance feature exhibiting sub-harmonic behaviour.

Accretion Spin-up and a Strong Magnetic Field in the Slow-spinning Be X-Ray Binary MAXI J0655-013

We present Monitor of All-sky X-ray Image (MAXI) and Nuclear Spectroscopic Telescope Array (NuSTAR) observations of the Be X-ray binary, MAXI J0655−013, in outburst. NuSTAR observed the source once early in the outburst, when spectral analysis yields a bolometric (0.1–100 keV), unabsorbed source luminosity of L bol = 5.6 × 1036 erg s−1, and a second time 54 days later, by which time the luminosity had dropped to L bol = 4 × 1034 erg s−1 after first undergoing a dramatic increase. Timing analysis of the NuSTAR data reveals a neutron star spin period of 1129.09 ± 0.04 s during the first observation, which decreased to 1085 ± 1 s by the time of the second observation, indicating spin-up due to accretion throughout the outburst. Furthermore, during the first NuSTAR observation, we observed quasi-periodic oscillations (QPOs) with centroid frequency ν 0 = 89 ± 1 mHz, which exhibited a second harmonic feature. By combining the MAXI and NuSTAR data with pulse period measurements reported by Fermi/GBM, we are able to show that apparent flaring behavior in the MAXI light curve is an artifact introduced by uneven sampling of the pulse profile, which has a large pulsed fraction. Finally, we estimate the magnetic field strength at the neutron star surface via three independent methods, invoking a tentative cyclotron resonance scattering feature at 44 keV, QPO production at the inner edge of the accretion disk, and spin-up via interaction of the neutron star magnetic field with accreting material. Each of these result in a significantly different value. We discuss the strengths and weaknesses of each method and infer that MAXI J0655−013 is likely to have a high surface magnetic field strength, B s > 1013 G.

Achromatic rapid flares in hard X-rays in the γ-ray binary LS I + 61 303

ABSTRACT We report on the presence of very rapid hard X-ray variability in the γ-ray binary LS I + 61 303. The results were obtained by analysing Nuclear Spectroscopic Telescope Array (NuSTAR) data, which show two achromatic strong flares on ks time-scales before apastron. The Swift/Burst Alert Telescope orbital X-ray light curve is also presented, and the NuSTAR data are put in the context of the system orbit. The spectrum and estimated physical conditions of the emitting region indicate that the radiation is synchrotron emission from relativistic electrons, likely produced in a shocked pulsar wind. The achromaticity suggests that losses are dominated by escape or adiabatic cooling in a relativistic flow, and the overall behaviour in hard X-rays can be explained by abrupt changes in the size of the emitting region and/or its motion relative to the line of sight, with Doppler boosting potentially being a prominent effect. The rapid changes of the emitter could be the result of different situations such as quick changes in the intra-binary shock, variations in the re-accelerated shocked pulsar wind outside the binary, or strong fluctuations in the location and size of the Coriolis shock region. Although future multiwavelength observations are needed to further constrain the physical properties of the high-energy emitter, this work already provides important insight into the complex dynamics and radiation processes in LS I + 61 303.

Broadband X-Ray Spectroscopy of the Pulsar Wind Nebula in HESS J1640-465

We present updated measurements of the X-ray properties of the pulsar wind nebula associated with the TeV γ-ray source HESS J1640-465 derived from Chandra and Nuclear Spectroscopic Telescope Array data. We report a high N H value along the line of sight, consistent with previous work, which led us to incorporate the effects of dust scattering in our spectral analysis. Due to uncertainties in the dust scattering, we report a range of values for the PWN properties (photon index and unabsorbed flux). In addition, we fit the broadband spectrum of this source and found evidence for spectral softening and decreasing unasborbed flux as we go to higher photon energies. We then used a one-zone time-dependent evolutionary model to reproduce the dynamical and multiwavelength spectral properties of our source. Our model suggests a short spin-down timescale, a relatively higher than average magnetized pulsar wind, a strong pulsar wind nebula magnetic field and maximum electron energy up to PeV, suggesting HESS J1640-465 could be a PeVatron candidate.

The First Survey of Quiet Sun Features Observed in Hard X-Rays with NuSTAR

We present the first survey of quiet Sun features observed in hard X-rays (HXRs), using the Nuclear Spectroscopic Telescope ARray (NuSTAR), a HXR focusing optics telescope. The recent solar minimum, combined with NuSTAR’s high sensitivity, has presented a unique opportunity to perform the first HXR imaging spectroscopy on a range of features in the quiet Sun. By studying the HXR emission of these features, we can detect or constrain the presence of high temperature (> 5 MK) or non-thermal sources, to help understand how they relate to larger, more energetic solar phenomena, and determine their contribution to heating the solar atmosphere. We report on several features observed in the 28 September 2018 NuSTAR full-disk quiet Sun mosaics, the first of the NuSTAR quiet Sun observing campaigns, which mostly include steady features of X-ray bright points and an emerging flux region, which later evolved into an active region, as well as a short-lived jet. We find that the features’ HXR spectra are well fitted with isothermal models with temperatures ranging between 2.0 – 3.2 MK. Combining the NuSTAR data with softer X-ray emission from Hinode/XRT and EUV from SDO/AIA, we recover the differential emission measures, confirming little significant emission above 4 MK. The NuSTAR HXR spectra allow us to constrain the possible non-thermal emission that would still be consistent with a null HXR detection. We found that for only one of the features (the jet) was there a potential non-thermal upper limit capable of powering the heating observed. However, even here, the non-thermal electron distribution had to be very steep (effectively mono-energetic) with a low energy cut-off between 3 – 4 keV.

The study of thermonuclear X-ray bursts in accreting millisecond pulsar MAXI J1816−195 with NuSTAR and NICER

ABSTRACT The millisecond pulsar MAXI J1816−195 was recently discovered in an outburst by the Monitor of All-sky X-ray Image (MAXI) in 2022 May. We study different properties of the pulsar using data from the Nuclear Spectroscopic Telescope Array (NuSTAR) and the Neutron Star Interior Composition Explorer (NICER) observations. The unstable burning of accreted material on the surface of neutron stars induces thermonuclear (Type-I) bursts. Several such thermonuclear bursts have been detected by MAXI J1816−195 during its outburst. We investigate the evolution of the burst profiles with flux and energy using NuSTAR and NICER observations. During the NuSTAR observation, a total of four bursts were detected from the source. The duration of each burst is around ∼30 s and the ratio of peak to persistent count rate is ∼26 as seen from the NuSTAR data. The burst profiles are modelled using a sharp linear rise and exponential decay function to determine the burst timing parameters. The burst profiles show a relatively long tail at lower energies. The broad-band time-resolved spectra during the burst periods are successfully modelled with a combination of an absorbed blackbody along with a non-thermal component to account for the persistent emission. From our modelling results, we are able to estimate the maximum apparent emitting area of the blackbody of the neutron star to be ∼12.5 km during the peak of the outburst and the maximum distance to the object to be 8.7 kpc. Our findings for the mass accretion rate and the α factor indicate the stable burning of hydrogen via a hot CNO cycle during the bursts.

Flux and spectral variability of Mrk 421 during its moderate activity state usingNuSTAR: Possible accretion disc contribution?

Context.The X-ray emission in BL Lac objects is believed to be dominated by synchrotron emission from their relativistic jets. However, when the jet emission is not strong, one could expect signatures of X-ray emission from inverse Compton scattering of accretion disc photons by hot and energetic electrons in the corona. Moreover, the observed X-ray variability can also originate in the disc, and get propagated and amplified by the jet.Aims.Here, we present results on the BL Lac object Mrk 421 using the Nuclear Spectroscopic Telescope Array data acquired during 2017 when the source was in a moderate X-ray brightness state. For comparison with high jet activity state, we also considered one epoch of data in April 2013 when the source was in a very high X-ray brightness state. Our aim is to explore the possibility of the signature of accretion disc emission in the overall X-ray emission from Mrk 421 and also examine changes in accretion parameters considering their contribution to spectral variations.Methods.We divided each epoch of data into different segments in order to find small-scale variability. Data for all segments were fitted using a simple power-law model. We also fitted the full epoch data using the two component advective flow (TCAF) model to extract the accretion flow parameters. Furthermore, we estimated the X-ray flux coming from the different components of the flow using the lowest normalisation method and analysed the relations between them. For consistency, we performed the spectral analysis using models available in the literature.Results.The simple power-law function does not fit the spectra well, and a cutoff needs to be added. The spectral fitting of the data using the TCAF model shows that the data can be explained with a model where (a) the size of the dynamic corona at the base of the jet is from ∼28 to 10 rs, (b) the disc mass accretion rate is from 0.021 to 0.051ṀEdd, (c) the halo mass accretion rate is from 0.22 to 0.35ṀEdd, and (d) the viscosity parameter of the Keplerian accretion disc from 0.18 to 0.25. In the assumed model, the total flux, disc and jet flux correlate with the radio flux observed during these epochs.Conclusions.From the spectral analysis, we conclude that the spectra of all the epochs of Mrk 421 in 2017 are well described by the accretion-disc-based TCAF model. The estimated disc and jet flux relations with radio flux show that accretion disc can contribute to the observed X-ray emission, when X-ray data (that cover a small portion of the broad band spectral energy distribution of Mrk 421) are considered in isolation. However, the present disc-based models are disfavoured with respect to the relativistic jet models when considering the X-ray data in conjunction with data at other wavelengths.

Reconstruction of the NuSTAR point spread function using single-laser metrology

We describe a method by which the metrology system of the Nuclear Spectroscopic Telescope Array (NuSTAR) X-ray space observatory, which uses two lasers to characterize the relative motion of the optics and focal plane benches, can be approximated should one laser fail. The two benches are separated by a 10-m-long rigid mast that undergoes small amounts of thermal flexing that needs to be compensated for to produce a nonblurred image. We analyze the trends of mast motion by archival observation parameters to discover whether the mast motion in future observations can be predicted. We find that, using the solar aspect angle, observation date, and orbital phase, we can simulate the motion of one laser by translating the track produced by the other and applying modifications to the resulting mast aspect solution, allowing for the reconstruction of a minimally distorted point spread function in most cases. We will implement the generation of simulated mast files along with the usual NuSTAR data reduction pipeline for contingency purposes. This work has implications for reducing the risk of implementing laser metrology systems on future missions that use deployable masts to achieve the long focal lengths required in high-energy astronomy by mitigating the impact of a metrology laser failure in the extended phase of a mission.

Compton-thick AGN in the NuSTAR Era VI: The Observed Compton-thick Fraction in the Local Universe

Abstract We present the analysis of simultaneous Nuclear Spectroscopic Telescope Array (NuSTAR) and XMM-Newton data of eight Compton-thick active galactic nuclei (CT-AGN) candidates selected in the Swift-BAT 100 month catalog. This work is part of an ongoing effort to find and characterize all CT-AGN in the Local (z ≤ 0.05) Universe. We used two physically motivated models, MYTorus and borus02, to characterize the sources in the sample, finding five of them to be confirmed CT-AGN. These results represent an increase of ∼19% over the previous NuSTAR-confirmed, BAT-selected CT-AGN at z ≤ 0.05, bringing the total number to 32. This corresponds to an observed fraction of ∼8% of all AGN within this volume-limited sample, although it increases to 20% ± 5% when limiting the sample to z ≤ 0.01. Out of a sample of 48 CT-AGN candidates, selected using BAT and soft (0.3−10 keV) X-ray data, only 24 are confirmed as CT-AGN with the addition of the NuSTAR data. This highlights the importance of NuSTAR when classifying local obscured AGN. We also note that most of the sources in our full sample of 48 Seyfert 2 galaxies with NuSTAR data have significantly different lines of sight and average torus column densities, favoring a patchy torus scenario.

NuSTAR observations of four nearby X-ray faint AGNs: low luminosity or heavy obscuration?

ABSTRACT We present NuSTAR (Nuclear Spectroscopic Telescope Array) observations of four active galactic nuclei (AGNs) located within 15 Mpc. These AGNs, namely ESO 121-G6, NGC 660, NGC 3486, and NGC 5195, have observed X-ray luminosities of L2–10 keV,obs ≲ 1039 erg s−1, classifying them as low-luminosity AGN (LLAGN). We perform broad-band X-ray spectral analysis for the AGN by combining our NuSTAR data with Chandra or XMM–Newton observations to directly measure their column densities (NH) and infer their intrinsic power. We complement our X-ray data with archival and new high-angular resolution mid-infrared (mid-IR) data for all objects, except NGC 5195. Based on our X-ray spectral analysis, we found that both ESO 121-G6 and NGC 660 are heavily obscured (NH > 1023 cm−2; L2–10 keV,int ∼ 1041 erg s−1), and NGC 660 may be Compton thick. We also note that the X-ray flux and spectral slope for ESO 121-G6 have significantly changed over the last decade, indicating significant changes in the obscuration and potentially accretion rate. On the other hand, NGC 3486 and NGC 5195 appear to be unobscured and just mildly obscured, respectively, with L2–10 keV,int < 1039 erg s−1, i.e. genuine LLAGN. Both of the heavily obscured AGNs have Lbol > 1041 erg s−1 and λEdd ≳ 10−3, and are detected in high-angular resolution mid-IR imaging, indicating the presence of obscuring dust on nuclear scale. NGC 3486, however, is undetected in high-resolution mid-IR imaging, and the current data do not provide stringent constraints on the presence or absence of obscuring nuclear dust in the AGN.

First Hard X-Ray Observation of a Compact Symmetric Object: A Broadband X-Ray Study of a Radio Galaxy OQ+208 with NuSTAR and Chandra

Abstract Compact symmetric objects (CSOs) have been observed with Chandra and XMM-Newton to gain insights into the initial stages of a radio source evolution and to probe the black hole activity at the time of relativistic outflow formation. However, there have been no CSO observations to date at the hard X-ray energies (>10 keV), impeding our ability to robustly constrain the properties of the intrinsic X-ray emission and of the medium surrounding the young expanding jets. We present the first hard X-ray observation of a CSO performed with the Nuclear Spectroscopic Telescope Array (NuSTAR). Our target, OQ +208, is detected up to 30 keV, and thus we establish CSOs as a new class of the NuSTAR sources. We analyze the NuSTAR data jointly with our new Chandra and archival XMM-Newton data and find that a young (∼250 yr old) radio jet spanning the length of 10 pc coexists with cold obscuring matter, consistent with a dusty torus, with an equivalent hydrogen column density of N H = 1023–1024 cm−2. The primary X-ray emission is characterized by a photon index of Γ ∼ 1.45 and an intrinsic 0.5–30 keV luminosity of L ≃ 1043 erg s−1. The results of our spectral modeling and broad-line optical classification of the source suggest a porous structure of the obscuring torus. Alternatively, the source may belong to the class of optically unobscured/X-ray-obscured active galactic nucleus. The observed X-ray emission is too weak compared to that predicted by the expanding radio lobes model, leaving an accretion disk corona or jets as the possible origins of the X-ray emission from this young radio galaxy.

XRB continuum fitting with sensitive high-energy X-ray detectors

The launch of the Nuclear Spectroscopic Telescope Array (NuSTAR) heralded a new era of sensitive high energy X-ray spectroscopy for X-ray binaries (XRBs). In this paper we show how multiple physical parameters can be measured from the accretion disk spectrum when the high-energy side of the disk spectrum can be measured precisely using NuSTAR. This immediately makes two exciting developments possible. If the mass and distance of the source are known, the continuum fitting method can be used to calculate the spin and inner disk inclination independently of the iron line fitting method. If the mass and distance are unknown, the two methods can be combined to constrain these values to a narrow region of parameter space. In this paper we perform extensive simulations to establish the reliability of these techniques. We find that with high quality spectra, spin and inclination can indeed be simultaneously measured using the disk spectrum. These measurements are much more precise at higher spin values, where the relativistic effects are stronger. The inclusion of a soft X-ray snapshot observation alongside the NuSTAR data significantly improves the reliability, particularly for lower temperature disks, as it gives a greatly improved measurement of the disk peak. High signal to noise data are not necessary for this, as measuring the peak temperature is relatively easy. We discuss the impact of systematic effects on this technique, and the implications of our results such as robust measurements of accretion disk warps and XRB mass surveys.

Compton-thick AGNs in the NuSTAR Era

Abstract We present the 2–100 keV spectral analysis of 30 candidate Compton-thick-(CT-)active galactic nuclei (AGNs) selected in the Swift-Burst Alert Telescope (BAT) 100 month survey. The average redshift of these objects is , and they all lie within ∼500 Mpc. We used the MyTorus model to perform X-ray spectral fittings both without and with the contribution of the Nuclear Spectroscopic Telescope Array (NuSTAR) data in the 3–50 keV energy range. When the NuSTAR data are added to the fit, 13 out of 30 of these objects (43% of the whole sample) have intrinsic absorption N H < 1024 cm−2 at the >3σ confidence level, i.e., they are reclassified from Compton thick to Compton thin. Consequently, we infer an overall observed fraction of the CT-AGN, with respect to the whole AGN population, lower than the one reported in previous works, as low as ∼4%. We find evidence that this overestimation of N H is likely due to the low quality of a subsample of spectra, either in the 2–10 keV band or in the Swift-BAT one.

The Sleeping Monster: NuSTAR Observations of SGR 1806–20, 11 Years After the Giant Flare

Abstract We report the analysis of five Nuclear Spectroscopic Telescope Array (NuSTAR) observations of SGR 1806−20 spread over a year from 2015 April to 2016 April, more than 11 years following its giant flare (GF) of 2004. The source spin frequency during the NuSTAR observations follows a linear trend with a frequency derivative Hz s−1, implying a surface dipole equatorial magnetic field G. Thus, SGR 1806−20 has finally returned to its historical minimum torque level measured between 1993 and 1998. The source showed strong timing noise for at least 12 years starting in 2000, with increasing one order of magnitude between 2005 and 2011, following its 2004 major bursting episode and GF. SGR 1806−20 has not shown strong transient activity since 2009, and we do not find short bursts in the NuSTAR data. The pulse profile is complex with a pulsed fraction of with no indication of energy dependence. The NuSTAR spectra are well fit with an absorbed blackbody, , plus a power law, . We find no evidence for variability among the five observations, indicating that SGR 1806−20 has reached a persistent and potentially its quiescent X-ray flux level after its 2004 major bursting episode. Extrapolating the NuSTAR model to lower energies, we find that the 0.5–10 keV flux decay follows an exponential form with a characteristic timescale days. Interestingly, the NuSTAR flux in this energy range is a factor of ∼2 weaker than the long-term average measured between 1993 and 2003, a behavior also exhibited in SGR 1900+14. We discuss our findings in the context of the magnetar model.

NUSTAR OBSERVATIONS OF THE YOUNG, ENERGETIC RADIO PULSAR PSR B1509–58

ABSTRACT We report on Nuclear Spectroscopic Telescope Array (NuSTAR) hard X-ray observations of the young rotation-powered radio pulsar PSR B1509−59 in the supernova remnant MSH 15−52. We confirm the previously reported curvature in the hard X-ray spectrum, showing that a log parabolic model provides a statistically superior fit to the spectrum compared with the standard power law. The log parabolic model describes the NuSTAR data, as well as previously published γ-ray data obtained with COMPTEL and AGILE, all together spanning 3 keV through 500 MeV. Our spectral modeling allows us to constrain the peak of the broadband high energy spectrum to be at 2.6 ± 0.8 MeV, an improvement of nearly an order of magnitude in precision over previous measurements. In addition, we calculate NuSTAR spectra in 26 pulse phase bins and confirm previously reported variations of photon indices with phase. Finally, we measure the pulsed fraction of PSR B1509−58 in the hard X-ray energy band for the first time. Using the energy resolved pulsed fraction results, we estimate that the pulsar’s off-pulse emission has a photon index value between 1.26 and 1.96. Our results support a model in which the pulsar’s lack of GeV emission is due to viewing geometry, with the X-rays originating from synchrotron emission from secondary pairs in the magnetosphere.

NuSTAROBSERVATIONS OF THE COMPTON-THICK ACTIVE GALACTIC NUCLEUS AND ULTRALUMINOUS X-RAY SOURCE CANDIDATE IN NGC 5643

We present two Nuclear Spectroscopic Telescope Array (NuSTAR) observations of the local Seyfert 2 active galactic nucleus (AGN) and an ultraluminous X-ray source (ULX) candidate in NGC 5643. Together with archival data from Chandra, XMM-Newton, and Swift-BAT, we perform a high-quality broadband spectral analysis of the AGN over two decades in energy (∼0.5–100 keV). Previous X-ray observations suggested that the AGN is obscured by a Compton-thick (CT) column of obscuring gas along our line of sight. However, the lack of high-quality ≳10 keV observations, together with the presence of a nearby X-ray luminous source, NGC 5643 X–1, have left significant uncertainties in the characterization of the nuclear spectrum. NuSTAR now enables the AGN and NGC 5643 X–1 to be separately resolved above 10 keV for the first time and allows a direct measurement of the absorbing column density toward the nucleus. The new data show that the nucleus is indeed obscured by a CT column of NH ≳ 5 × 1024 cm−2. The range of 2–10 keV absorption-corrected luminosity inferred from the best-fitting models is L2–10,int = (0.8–1.7) × 1042 erg s−1, consistent with that predicted from multiwavelength intrinsic luminosity indicators. In addition, we also study the NuSTAR data for NGC 5643 X–1 and show that it exhibits evidence of a spectral cutoff at energy E ∼ 10 keV, similar to that seen in other ULXs observed by NuSTAR. Along with the evidence for significant X-ray luminosity variations in the 3–8 keV band from 2003 to 2014, our results further strengthen the ULX classification of NGC 5643 X–1.

NuSTAROBSERVATIONS OF THE MAGNETAR 1E 2259+586

We report on new broad band spectral and temporal observations of the magnetar 1E 2259+586, which is located in the supernova remnant CTB 109. Our data were obtained simultaneously with the Nuclear Spectroscopic Telescope Array (NuSTAR) and Swift, and cover the energy range from 0.5-79 keV. We present pulse profiles in various energy bands and compare them to previous RXTE results. The NuSTAR data show pulsations above 20 keV for the first time and we report evidence that one of the pulses in the double-peaked pulse profile shifts position with energy. The pulsed fraction of the magnetar is shown to increase strongly with energy. Our spectral analysis reveals that the soft X-ray spectrum is well characterized by an absorbed double-blackbody or blackbody plus power-law model in agreement with previous reports. Our new hard X-ray data, however, suggests that an additional component, such as a power-law, is needed to describe the NuSTAR and Swift spectrum. We also fit the data with the recently developed coronal outflow model by Beloborodov for hard X-ray emission from magnetars. The outflow from a ring on the magnetar surface is statistically preferred over outflow from a polar cap.