Spectral and Timing Analysis of the accretion-powered pulsar 4U 1626-67 observed with Suzaku and NuSTAR

We present an analysis of the spectral shape and pulse profile of the accretion-powered pulsar 4U 1626-67 observed with Suzaku and NuSTAR during a spin-up state. The pulsar, which experienced a torque reversal to spin-up in 2008, has a spin period of 7.7 s. Comparing the phase-averaged spectra obtained with Suzaku in 2010 and with NuSTAR in 2015, we find that the spectral shape changed between the two observations: the 3-10 keV flux increased by 5% while the 30-60 keV flux decreased significantly by 35%. Phase-averaged and phase-resolved spectral analysis shows that the continuum spectrum observed by NuSTAR is well described by an empirical NPEX continuum with an added broad Gaussian emission component around the spectral peak at 20 keV. Taken together with the observed Pdot value obtained from Fermi/GBM, we conclude that the spectral change between the Suzaku and NuSTAR observations was likely caused by an increase of the accretion rate. We also report the possible detection of asymmetry in the profile of the fundamental cyclotron line. Furthermore, we present a study of the energy-resolved pulse profiles using a new relativistic ray tracing code, where we perform a simultaneous fit to the pulse profiles assuming a two-column geometry with a mixed pencil- and fan-beam emission pattern. The resulting pulse profile decompositions enable us to obtain geometrical parameters of accretion columns (inclination, azimuthal and polar angles) and a fiducial set of beam patterns. This information is important to validate the theoretical predictions from radiation transfer in a strong magnetic field.Comment: 19 pages, 14 figures, Accepted for publication in ApJ on May 5, 201

The giant outburst of 4U 0115+634 in 2011 with Suzaku and RXTE - Minimizing cyclotron line biases

We present an analysis of X-ray spectra of the high-mass X-ray binary 4U 0115+634 as observed with Suzaku and RXTE in 2011 July, during the fading phase of a giant X-ray outburst. We used a continuum model consisting of an absorbed cutoff power law and an ad hoc Gaussian emission feature centered around 8.5 keV, which we attribute to cyclotron emission. Our results are consistent with a fundamental cyclotron absorption line centered at ∼10.2 keV for all observed flux ranges. At the same time we rule out significant influence of the 8.5 kev Gaussian on the parameters of the cyclotron resonant scattering feature, which are not consistent with the cyclotron line energies or the depths of previously reported flux-dependent descriptions. We also show that some continuum models can lead to artificial line-like residuals in the analyzed spectra, which are then misinterpreted as unphysically strong cyclotron lines. Specifically, our results do not support the existence of a previously claimed additional cyclotron feature at ∼15 keV. Apart from these features, we find for the first time evidence for a He-like Fe XXV emission line at ∼6.7 keV and weak H-like Fe XXVI emission close to ∼7.0 keV.We acknowledge funding by the European Space Agency under contract number C4000115860/15/NL/IB, by the Bundesministerium für Wirtschaft und Technologie under Deutsches Zentrum für Luft- und Raumfahrt grants 50OR0808, 50OR0905, 50OR1113, and 50OR1207, and by the Deutscher Akademischer Austauschdienst. MTW is supported by the NASA Astrophysical Data Analysis Program and the Chief of Naval Research. VG is supported through the Margarethe von Wrangell fellowship by the ESF and the Ministry of Science, Research and the Arts Baden-Württemberg. SMN and JMT acknowledge Spanish Ministerio de Ciencia, Tecnología e Innovación (MICINN) through the grant ESP2016-76683-C3-1-R and ESP2017-85691-P, respectively

The Transient Accereting X-Ray Pulsar XTE J1946+274: Stability of the X-Ray Properties at Low Flux and Updated Orbital Solution

We present a timing and spectral analysis of the X-ray pulsar XTE J1946+274 observed with Suzaku during an outburst decline in 2010 October and compare with previous results. XTE J1946+274 is a transient X-ray binary consisting of a Be-type star and a neutron star with a 15.75 s pulse period in a 172 days orbit with 2–3 outbursts per orbit during phases of activity. We improve the orbital solution using data from multiple instruments. The X-ray spectrum can be described by an absorbed Fermi–Dirac cut-off power-law model along with a narrow Fe Kα line at 6.4 keV and a weak Cyclotron Resonance Scattering Feature (CRSF) at ~35 keV. The Suzaku data are consistent with the previously observed continuum flux versus iron line flux correlation expected from fluorescence emission along the line of sight. However, the observed iron line flux is slightly higher, indicating the possibility of a higher iron abundance or the presence of non-uniform material. We argue that the source most likely has only been observed in the subcritical (non-radiation dominated) state since its pulse profile is stable over all observed luminosities and the energy of the CRSF is approximately the same at the highest (~5 × 10^(37) erg s^(−1)) and lowest (~5 × 10^(36) erg s^(−1)) observed 3–60 keV luminosities

Spectral and Timing Analysis of the Accretion-powered Pulsar 4U 1626–67 Observed with Suzaku and NuSTAR

We present an analysis of the spectral shape and pulse profile of the accretion-powered pulsar 4U 1626−67 observed with Suzaku and Nuclear Spectroscopic Telescope Array (NuSTAR) during a spin-up state. The pulsar, which experienced a torque reversal to spin-up in 2008, has a spin period of ~7.7 s. Comparing the phase-averaged spectra obtained with Suzaku in 2010 and with NuSTAR in 2015, we find that the spectral shape changed between the two observations: the 3–10 keV flux increased by ~5%, while the 30–60 keV flux decreased significantly by ~35%. Phase-averaged and phase-resolved spectral analysis shows that the continuum spectrum observed by NuSTAR is well described by an empirical negative and positive power law times exponential continuum with an added broad Gaussian emission component around the spectral peak at ~20 keV. Taken together with the observed Ṗ value obtained from the Fermi/gamma-ray burst monitor data, we conclude that the spectral change between the Suzaku and NuSTAR observations was likely caused by an increase in the accretion rate. We also report the possible detection of asymmetry in the profile of the fundamental cyclotron line. Furthermore, we present a study of the energy-resolved pulse profiles using a new relativistic ray tracing code, where we perform a simultaneous fit to the pulse profiles assuming a two-column geometry with a mixed pencil- and fan-beam emission pattern. The resulting pulse profile decompositions enable us to obtain geometrical parameters of accretion columns (inclination, azimuthal and polar angles) and a fiducial set of beam patterns. This information is important to validate the theoretical predictions from radiation transfer in a strong magnetic field

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, for example, a cutoff power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. Previous works have presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this solution that is both fast and accurate enough to be fitted to observed spectra is available as a model in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during the χ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3. The derived accretion-flow structure features a small column radius of ∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies