Spectral Changes in the Hyperluminous Pulsar in NGC 5907 as a Function of Super-Orbital Phase

We present broad-band, multi-epoch X-ray spectroscopy of the pulsating ultra-luminous X-ray source (ULX) in NGC 5907. Simultaneous XMM-Newton and NuSTAR data from 2014 are best described by a multi-color black-body model with a temperature gradient as a function of accretion disk radius significantly flatter than expected for a standard thin accretion disk (T(r) ~ r^{-p}, with p=0.608^{+0.014}_{-0.012}). Additionally, we detect a hard power-law tail at energies above 10 keV, which we interpret as being due to Comptonization. We compare this observation to archival XMM-Newton, Chandra, and NuSTAR data from 2003, 2012, and 2013, and investigate possible spectral changes as a function of phase over the 78d super-orbital period of this source. We find that observations taken around phases 0.3-0.4 show very similar temperature profiles, even though the observed flux varies significantly, while one observation taken around phase 0 has a significantly steeper profile. We discuss these findings in light of the recent discovery that the compact object is a neutron star and show that precession of the accretion disk or the neutron star can self-consistently explain most observed phenomena.Comment: 7 pages, 5 figures, submitted to ApJ; comments welcom

Multiple cyclotron line-forming regions in GX 301-2

We present two observations of the high-mass X-ray binary GX 301-2 with NuSTAR, taken at different orbital phases and different luminosities. We find that the continuum is well described by typical phenomenological models, like a very strongly absorbed NPEX model. However, for a statistically acceptable description of the hard X-ray spectrum we require two cyclotron resonant scattering features (CRSF), one at ~35 keV and the other at ~50 keV. Even though both features strongly overlap, the good resolution and sensitivity of NuSTAR allows us to disentangle them at >=99.9% significance. This is the first time that two CRSFs are seen in GX 301-2. We find that the CRSFs are very likely independently formed, as their energies are not harmonically related and, if it were a single line, the deviation from a Gaussian shape would be very large. We compare our results to archival Suzaku data and find that our model also provides a good fit to those data. We study the behavior of the continuum as well as the CRSF parameters as function of pulse phase in seven phase bins. We find that the energy of the 35 keV CRSF varies smoothly as function of phase, between 30-38 keV. To explain this variation, we apply a simple model of the accretion column, taking the altitude of the line-forming region, the velocity of the in-falling material, and the resulting relativistic effects into account. We find that in this model the observed energy variation can be explained simply due to a variation of the projected velocity and beaming factor of the line forming region towards us.Comment: 18 pages, 10 figures, accepted for publication in A&

A Study of the 20 Day Superorbital Modulation in the High-Mass X-ray Binary IGR J16493-4348

We report on Nuclear Spectroscopic Telescope Array (NuSTAR), Neil Gehrels Swift Observatory (Swift) X-ray Telescope (XRT) and Swift Burst Alert Telescope (BAT) observations of IGR J16493-4348, a wind-fed Supergiant X-ray Binary (SGXB) showing significant superorbital variability. From a discrete Fourier transform of the BAT light curve, we refine its superorbital period to be 20.058 $\pm$ 0.007 days. The BAT dynamic power spectrum and a fractional root mean square analysis both show strong variations in the amplitude of the superorbital modulation, but no observed changes in the period were found. The superorbital modulation is significantly weaker between MJD 55,700 and MJD 56,300. The joint NuSTAR and XRT observations, which were performed near the minimum and maximum of one cycle of the 20 day superorbital modulation, show that the flux increases by more than a factor of two between superorbital minimum and maximum. We find no significant changes in the 3-50 keV pulse profiles between superorbital minimum and maximum, which suggests a similar accretion regime. Modeling the pulse-phase averaged spectra we find a possible Fe K$\alpha$ emission line at 6.4 keV at superorbital maximum. The feature is not significant at superorbital minimum. While we do not observe any significant differences between the pulse-phase averaged spectral continua apart from the overall flux change, we find that the hardness ratio near the broad main peak of the pulse profile increases from superorbital minimum to maximum. This suggests the spectral shape hardens with increasing luminosity. We discuss different mechanisms that might drive the observed superorbital modulation.Comment: 17 pages, 14 figures, 3 tables, accepted for publication in The Astrophysical Journal on 2019 May 1

Spectral and Timing Nature of the Symbiotic X-ray Binary 4U 1954+319: The Slowest Rotating Neutron Star in an X-ray Binary System

The symbiotic X-ray binary 4U 1954+319 is a rare system hosting a peculiar neutron star (NS) and an M-type optical companion. Its ~5.4h NS spin period is the longest among all known accretion-powered pulsars and exhibited large (~7%) fluctuations over 8 years. A spin trend transition was detected with Swift/BAT around an X-ray brightening in 2012. The source was in quiescent and bright states before and after this outburst based on 60 ks Suzaku observations in 2011 and 2012. The observed continuum is well described by a Comptonized model with the addition of a narrow 6.4 keV Fe Kalpha line during the outburst. Spectral similarities to slowly rotating pulsars in high-mass X-ray binaries, its high pulsed fraction (~60-80%), and the location in the Corbet diagram favor high B-field (>~1e+12 G) over a weak field as in low-mass X-ray binaries. The observed low X-ray luminosity (1e+33-1e+35 erg/s), probable wide orbit, and a slow stellar wind of this SyXB make quasi-spherical accretion in the subsonic settling regime a plausible model. Assuming a ~1e+13 G NS, this scheme can explain the ~5.4 h equilibrium rotation without employing the magnetar-like field (~1e+16 G) required in the disk accretion case. The time-scales of multiple irregular flares (~50 s) can also be attributed to the free-fall time from the Alfven shell for a ~1e+13 G field. A physical interpretation of SyXBs beyond the canonical binary classifications is discussed.Comment: 20 pages, 18 figures, 3 tables, accepted for publication in the Astrophysical Journa

Disentangling the Complex Broadband X-ray Spectrum of IRAS 13197-1627 with NuSTAR, XMM-Newton and Suzaku

We present results from a coordinated $XMM$-$Newton$+$NuSTAR$ observation of the type 1.8 Seyfert galaxy IRAS 13197-1627. This is a highly complex source, with strong contributions from relativistic reflection from the inner accretion disk, neutral absorption and further reprocessing by more distant material, and ionised absorption from an outflow. We undertake a detailed spectral analysis combining the broadband coverage provided by $XMM$-$Newton$+$NuSTAR$ with a multi-epoch approach incorporating archival observations performed by $XMM$-$Newton$ and $Suzaku$. Our focus is on characterising the reflection from the inner accretion disk, which previous works have suggested may dominate the AGN emission, and constraining the black hole spin. Using lamppost disk reflection models, we find that the results for the inner disk are largely insensitive to assumptions regarding the geometry of the distant reprocessor and the precise form of the illuminating X-ray continuum. However, these results do depend on the treatment of the iron abundance of the distant absorber/reprocessor. The multi-epoch data favour a scenario in which the AGN is chemically homogeneous, and we find that a rapidly rotating black hole is preferred, with $a^* \geq 0.7$, but a slowly-rotating black hole is not strongly excluded. In addition to the results for the inner disk, we also find that both the neutral and ionised absorbers vary from epoch to epoch, implying that both have some degree of inhomogeneity in their structure.Comment: 15 pages, 7 figures, accepted for publication in MNRA