Rossi X-ray Timing Explorer Observations Research Articles

Measuring the Magnetic Dipole Moment and Magnetospheric Fluctuations of SXP 18.3 with a Kalman Filter

The magnetic dipole moment μ of an accretion-powered pulsar in magnetocentrifugal equilibrium cannot be inferred uniquely from time-averaged pulse period and aperiodic X-ray flux data, because the radiative efficiency η 0 of the accretion is unknown, as are the mass, radius, and distance of the star. The degeneracy associated with the radiative efficiency is circumvented if fluctuations of the pulse period and aperiodic X-ray flux are tracked with a Kalman filter, whereupon μ can be measured uniquely up to the uncertainties in the mass, radius, and distance. Here, the Kalman filter analysis is demonstrated successfully in practice for the first time on Rossi X-ray Timing Explorer observations of the X-ray transient SXP 18.3 in the Small Magellanic Cloud (SMC), which is monitored regularly. The analysis yields μ=8.0−1.2+1.3×1030Gcm3 and η0=0.04−0.01+0.02 , compared to μ=5.0−1.0+1.0×1030Gcm3 as inferred traditionally from time-averaged data assuming η 0 = 1. The analysis also yields time-resolved estimates of two hidden state variables, the mass accretion rate and the Maxwell stress at the disk–magnetosphere boundary. The success of the demonstration confirms that the Kalman filter analysis can be applied in the future to study the magnetic moments and disk–magnetosphere physics of accretion-powered pulsar populations in the SMC and elsewhere.

Investigation into the reflection properties of the neutron star low-mass X-ray binary 4U 1636–53

We present the spectroscopy of the neutron star low-mass X-ray binary 4U 1636–53 using six simultaneous XMM-Newton and Rossi X-ray Timing Explorer observations. We applied different self-consistent reflection models to explore the features when the disk is illuminated by either the corona or the neutron star surface. We found that the spectra could be well fitted by these two types of models, with the derived emissivity index below a typical value of 3. The relative low emissivity can be explained if the neutron star and the corona, working together as an extended illuminator, simultaneously illuminate and ionize the disk. Additionally, the derived ionization parameter in the lamppost geometry is larger than the theoretical prediction. This inconsistency likely suggests that the corona does not emit isotropically in a realistic context. Furthermore, we also found that there is a possible trend between the height of the corona and the normalization of the disk emission. This could be understood either as a variation in the reflected radiation pressure or in the context of a jet base. Finally, we found that the disk is less ionized if it is illuminated by the neutron star, indicating that the illuminating source has significant influence on the physical properties of the disk.

Similarities and differences in accretion flow properties between GRS 1915+105 and IGR J17091-3624: A case study

We perform a comparative spectro-temporal analysis on the variability classes of GRS 1915+105 and IGR J17091-3624 to draw inferences regarding the underlying accretion flow mechanism. The ν, as well as C2 class Rossi X-ray Timing Explorer observation, have been considered for analysis. We investigate the intensity variation of the source in different energy domains that correspond to different components of the accretion flow and infer the relative dominance of these flow components during the dip/flare events. We correlate the dependence of the dynamic photon index (Θ) with intensities in different energy bands and comment on the transition of the source to hard/soft phases during soft dips/flares. We also report the presence of sharp QPOs at ∼7.1 Hz corresponding to both softer and harder domain in the case of ν variability class of GRS 1915+105 and discuss the possible accretion flow configuration it suggests. Sharp QPO around ∼20 mHz is observed in ν and C2 classes of IGR J17091-3624 in low and mid energy band (2.0–6.0 keV and 6.0–15.0 keV), but remains undetected in high energy (15.0–60.0 keV). The 2.5–25.0 keV background-subtracted spectra have also been fitted with TCAF along with a Compton reflection component. A plausible accretion flow mechanism in order to explain the observed variability has been proposed.

Pulsar identification method based on adaptive grey wolf optimization algorithm in X-ray pulsar-based navigations

Fast and accurate identification of pulsar signals is important for X-ray pulsar-based navigation (XNAV). Traditional pulsar signal detection technology based on FFT search and epoch folding (EF) requires a very long time to obtain an appropriate signal-to-noise ratio (SNR) gain, especially for weak pulsar signals with low photon fluxes. This paper proposes an adaptive stochastic resonance (SR) method based on the grey wolf optimizer (GWO) algorithm for fast pulsar identification. The GWO algorithm can optimize the SR parameters to obtain the optimal SR output that matches the input signal. Meanwhile, the initial population of GWO is generated by Chebyshev chaotic maps, which guarantees the diversity of the initial population and enhances the global search capability of GWO. To further improve the performance of SR, we add an extra denoising process, including a high-pass filter and singular value decomposition (SVD) that determines the effective singular components through dispersion entropy (DE), to suppress the noise to a great extent. We use simulation data and Rossi X-ray Timing Explorer (RXTE) observation data to verify the proposed method. The output SNR is greatly improved, and the method can identify pulsars in a short time. The results indicated that the proposed method has great practical value in engineering. Further experiments reveal the importance of SVD, DE and the Chebyshev chaotic map for the algorithm.

The Harmonic Component of the Millihertz Quasi-periodic Oscillations in 4U 1636–53

Abstract We studied the harmonics of the millihertz quasi-periodic oscillations (mHz QPOs) in the neutron star low-mass X-ray binary 4U 1636–53 using the Rossi X-ray Timing Explorer observations. We detected the harmonics of the mHz QPOs in 73 data intervals, with most of them in the transitional spectra state. We found that the ratio between the rms amplitude of the harmonic and that of the fundamental remains constant in a wide range of the fundamental frequency. More importantly, we studied, for the first time, the rms amplitude of the harmonics versus energy in 4U 1636–53 in the 2–5 keV range. We found that the rms amplitude of both the harmonic and the fundamental shows a decreasing trend as the energy increases, which is different from the behaviors reported in QPOs in certain black hole systems. Furthermore, our results suggest that not all observations with mHz QPOs have the harmonic component, although the reason behind this is still unclear.

Light-curve fingerprints: an automated approach to the extraction of X-ray variability patterns with feature aggregation – an example application to GRS 1915+105

ABSTRACT Time series data mining is an important field of research in the era of ‘Big Data’. Next generation astronomical surveys will generate data at unprecedented rates, creating the need for automated methods of data analysis. We propose a method of light-curve characterization that employs a pipeline consisting of a neural network with a long-short term memory variational autoencoder architecture and a Gaussian mixture model. The pipeline performs extraction and aggregation of features from light-curve segments into feature vectors of fixed length that we refer to as light-curve ‘fingerprints’. This representation can be readily used as input of down-stream machine learning algorithms. We demonstrate the proposed method on a data set of Rossi X-ray Timing Explorer observations of the Galactic black hole X-ray binary GRS 1915+105, which was chosen because of its observed complex X-ray variability. We find that the proposed method can generate a representation that characterizes the observations and reflects the presence of distinct classes of GRS 1915+105 X-ray flux variability. We find that this representation can be used to perform efficient classification of light curves. We also present how the representation can be used to quantify the similarity of different light curves, highlighting the problem of the popular classification system of GRS 1915+105 observations, which does not account for intermediate class behaviour.

Power spectral properties of the soft spectral states in four black hole transients

The X-ray variability in the soft X-ray spectral state of black hole binaries is primarily characterized by a power-law noise (PLN), which is thought to originate from the propagation of the modulation in the mass accretion rate of a standard accretion disk flow. Such a PLN has also been revealed in the disk spectral component in the hard and the intermediate states in several black hole binaries. Here we present an investigation of the Rossi X-ray Timing Explorer (RXTE) observations of four black hole transients in which soft spectral states were observed twenty times or more. We show that in the soft spectral state, the PLN index varied in a large range between –1.64 and –0.62, and the fractional rms variability calculated in the 0.01 – 20 Hz frequency range reached as large as 7.67% and as low as 0.83%. Remarkably, we have found evidence of an inclination dependence of the maximal fractional rms variability, the averaged fractional rms variability and the fractional rms variability of the median in the sample based on current knowledge of inclination of black hole binaries. An inclination dependence has only been predicted in early magnetohydrodynamic simulations of isothermal disks limited to a high-frequency regime. In theory, the noise index is related to the physics of inward propagation of disk fluctuations, while the fractional rms amplitude reflects the intrinsic properties of the magnetohydrodynamic nature of the accretion flow. Our results therefore suggest that X-ray variability in the soft state can be used to put constraints on the properties of the accretion flow as well as the inclination of the accretion disk.

The Black Hole Spin in GRS 1915+105, Revisited

Abstract We estimate the black hole spin parameter in GRS 1915+105 using the continuum-fitting method with revised mass and inclination constraints based on the very long baseline interferometric parallax measurement of the distance to this source. We fit Rossi X-ray Timing Explorer observations selected to be accretion-disk-dominated spectral states as described in McClintock et al. (2006) and Middleton et al. (2006), which previously gave discrepant spin estimates with this method. We find that, using the new system parameters, the spin in both data sets increased, providing a best-fit spin of a * = 0.86 for the Middleton et al. data and a poor fit for the McClintock et al. data set, which becomes pegged at the BHSPEC model limit of a * = 0.99. We explore the impact of the uncertainties in the system parameters, showing that the best-fit spin ranges from a * = 0.4 to 0.99 for the Middleton et al. data set and allows reasonable fits to the McClintock et al. data set with near-maximal spin for system distances greater than ∼10 kpc. We discuss the uncertainties and implications of these estimates.

Study on the Energy Limits of kHz QPOs in Sco X-1 with RXTE and Insight-HXMT Observations

Abstract We present a detailed spectral-timing analysis of the Kilohertz quasiperiodic oscillations (kHz QPOs) in Sco X-1 using the data of the Rossi X-ray Timing Explorer (RXTE) and the Hard X-ray Modulation Telescope (Insight-HXMT). The energy band with detectable kHz QPOs is studied for the first time: on the horizontal branch, it is ∼6.89–24.01 and ∼8.68–21.78 keV for the upper and lower kHz QPOs, respectively, detected by the RXTE, and ∼9–27.5 keV for the upper kHz QPOs by the Insight-HXMT; on the lower normal branch, the energy band is narrower. The fractional root mean square (rms) of the kHz QPOs increases with energy at a lower energy, reaches a plateau at about 16 and 20 keV for the lower and upper peaks, and then levels off though with a large uncertainty. The simulation of the deadtime effect of RXTE/PCA shows that the deadtime does not affect much the search of the kHz QPOs but makes the rms amplitude underestimated. No significant QPO is detected below ∼6 keV as shown by the RXTE data, implying that the kHz QPOs do not originate from the blackbody emission of the accretion disk and neutron star surface. In addition, with the combined analysis of the energy spectra and the absolute rms spectra of kHz QPOs, we suggest that the kHz QPOs in Sco X-1 originate from the Comptonization of the inner part of the transition layer, where the rotation sets the frequency and the inward bulk motion makes the spectrum harder.

Long-term Coherent Timing of the Accreting Millisecond Pulsar IGR J17062–6143

Abstract We report on a coherent timing analysis of the 163 Hz accreting millisecond X-ray pulsar IGR J17062–6143. Using data collected with the Neutron Star Interior Composition Explorer and XMM-Newton, we investigated the pulsar evolution over a time span of four years. We obtained a unique phase-coherent timing solution for the stellar spin, finding the source to be spinning up at a rate of (3.77 ± 0.09) × 10−15 Hz s−1. We further find that the 0.4–6 keV pulse fraction varies gradually between 0.5% and 2.5% following a sinusoidal oscillation with a 1210 ± 40 day period. Finally, we supplemented this analysis with an archival Rossi X-ray Timing Explorer observation and obtained a phase-coherent model for the binary orbit spanning 12 yr, yielding an orbital period-derivative measurement of (8.4 ± 2.0) × 10−12 s s−1. This large orbital period derivative is inconsistent with a binary evolution that is dominated by gravitational wave emission and is suggestive of highly nonconservative mass transfer in the binary system.

Very low-frequency oscillations from the 11 Hz pulsar in Terzan 5: frame dragging back on the table.

ABSTRACT We present a re-analysis of 47 Rossi X-ray Timing Explorer observations of the 11 Hz accreting pulsar IGR J17480-2446 in Terzan 5 during its 2010 outburst. We studied the fast-time variability properties of the source and searched for quasi-periodic oscillations (QPOs) in a large-frequency range. General relativity predicts that frame dragging occurs in the vicinity of a spinning compact object and induces the precession of matter orbiting said object. The relativistic precession model predicts that this frame dragging can be observed as QPOs with a characteristic frequency in the light curves of accreting compact objects. Such QPOs have historically been classified as horizontal branch oscillations in neutron star (NS) systems, and for a NS spinning at 11 Hz, these oscillations are expected at frequencies below 1 Hz. However, previous studies of IGR J17480-2446 have classified QPOs at 35–50 Hz as horizontal branch oscillations, thus casting doubts on the frame-dragging nature of such QPOs. Here, we report the detection of seven very low-frequency QPOs, previously undetected, with centroid frequencies below 0.3 Hz, and which can be ascribed to frame dragging. We also discuss the possible nature of the QPOs detected at 35–50 Hz in this alternative scenario.

A minimal timescale for the continuum in 4U 1608-52 and Aql X-1

ABSTRACT Similar to black hole X-ray binary transients, hysteresis-like state transitions are also seen in some neutron-star X-ray binaries. Using a method based on wavelets and light curves constructed from archival Rossi X-ray Timing Explorer observations, we extract a minimal timescale over the complete range of transitions for 4U 1608-52 during the 2002 and 2007 outbursts and the 1999 and 2000 outbursts for Aql X-1. We present evidence for a strong positive correlation between this minimal timescale and a similar timescale extracted from the corresponding power spectra of these sources.

A QPO in NGC 4945 from Archival RXTE Data

Abstract We report the discovery of an ∼6 week quasi-periodic oscillation (QPO) in archival NGC 4945 data observed by the Rossi X-ray Timing Explorer (RXTE) satellite. QPOs are important observables in accretion disks and have been studied extensively in both neutron star (NS) and black hole (BH) X-ray binaries (XRB). QPOs should be present in active galactic nuclei (AGN) if galactic BHs and supermassive black holes are governed by a common set of physical processes. The search for QPOs in AGN has proven difficult because the time scales would be much longer, due to their higher mass. RXTE AGN light curves spanning 1996 to 2011 provide an excellent and perhaps unique opportunity to search for low-frequency QPOs. We investigated the 533 RXTE observations made of the Seyfert-2 AGN, NGC 4945. During a large cluster of observations in 2006–2007 (194 observations, spanning 396 days), the Lomb–Scargle periodogram shows a candidate QPO at 0.274 μHz (period ≈42.2 days). We estimate the uncertainties using the false alarm probability. We discuss the possible identification of this feature with the Lense–Thirring precession period.

Study of long-term flux and photon index distributions of blazars using RXTE observations

ABSTRACT We present a detailed study of flux and index distributions of three blazars [one flat-spectrum radio quasar (FSRQ) and two BL Lacertae objects (BL Lacs)] by using 16 yr of Rossi X-ray Timing Explorer (RXTE) archival data. The three blazars were chosen such that their flux and index distributions have sufficient number of data points (≥90) with relatively less uncertainty $\left(\overline{\sigma _{\rm err}^{2}}/\sigma ^{2} < 0.2\right)$ in light curves. Anderson–Darling (AD) test and histogram fitting show that flux distribution of FSRQ 3C 273 is lognormal, while its photon index distribution is Gaussian. This result is consistent with linear Gaussian perturbation in the particle acceleration time-scale, which produces lognormal distribution in flux. However, for two BL Lacs, viz. Mrk 501 and Mrk 421, AD test shows that their flux distributions are neither Gaussian nor lognormal, and their index distributions are non-normal. The histogram fitting of Mrk 501 and Mrk 421 suggests that their flux distributions are more likely to be a bimodal, and their index distributions are double Gaussian. Since, Sinha et al. had shown that Gaussian distribution of index produces a lognormal distribution in flux, double Gaussian distribution of index in Mrk 501 and Mrk 421 indicates that their flux distributions are probably double lognormal. Observation of double lognormal flux distribution with double Gaussian distribution in index reaffirms two flux states hypothesis. Further, the difference observed in the flux distribution of FSRQ (3C 273) and BL Lacs (Mrk 501 and Mrk 421) at X-rays suggests that the low-energy emitting electrons have a single lognormal flux distribution, while the high-energy ones have a double lognormal flux distribution.

Discovery of an Accretion-rate Independent Absolute RMS Amplitude of Millihertz Quasi-periodic Oscillations in 4U 1636-53

Abstract We investigate the frequency and amplitude of the millihertz quasi-periodic oscillations (mHz QPOs) in the neutron-star low-mass X-ray binary 4U 1636–53 using Rossi X-ray Timing Explorer observations. We find that no mHz QPOs appear when the source is in the hard spectral state. We also find that there is no significant correlation between the frequency and the fractional RMS amplitude of the mHz QPOs. Notwithstanding, for the first time, we find that the absolute rms amplitude of the mHz QPOs is insensitive to the parameter S a , which measures the position of the source in the color–color diagram and is usually assumed to be an increasing function of mass accretion rate. This finding indicates that the transition from marginally stable burning to stable burning or unstable burning could happen very rapidly since, before the transition, the mHz QPOs do not gradually decay as the rate further changes.

XMMU J181227.8–181234: a new ultracompact X-ray binary candidate

We report the discovery of Type I (thermonuclear) X-ray bursts from the transient source XMMU J181227.8–181234 = XTE J1812–182. We found seven X-ray bursts in Rossi X-ray Timing Explorer observations during the 2008 outburst, confirming the source as a neutron star low-mass X-ray binary. Based on the measured burst fluence and the average recurrence time of 1.4|$^{+0.9}_{-0.5}$| h, we deduce that the source is accreting almost pure helium (X ≤ 0.1) fuel. Two bursts occurred just 18 min apart; the first short waiting time bursts observed in a source accreting hydrogen-poor fuel. Taking into consideration the effects on the burst and persistent flux due to the inferred system inclination of 30 ± 10°, we estimate the distance to be 14 ± 2 kpc, where we report the statistical uncertainty but note that there could be up to |$20{{\ \rm per\ cent}}$| variation in the distance due to systematic effects discussed in the paper. The corresponding maximum accretion rate is 0.30 ± 0.05 times the Eddington limit. Based on the low hydrogen content of the accreted fuel and the short average recurrence time, we classify the source as a transient ultracompact low-mass X-ray binary.

A systematic study of the phase difference between QPO harmonics in black hole X-ray binaries

We perform a systematic study of the evolution of the waveform of black hole X-ray binary low-frequency QPOs, by measuring the phase difference between their fundamental and harmonic features. This phase difference has been studied previously for small number of QPO frequencies in individual sources. Here, we present a sample study spanning fourteen sources and a wide range of QPO frequencies. With an automated pipeline, we systematically fit power spectra and calculate phase differences from archival Rossi X-ray Timing Explorer (RXTE) observations. We measure well-defined phase differences over a large range of QPO frequencies for most sources, demonstrating that a QPO for a given source and frequency has a persistent underlying waveform. This confirms the validity of recently developed spectral-timing methods performing phase resolved spectroscopy of the QPO. Furthermore, we evaluate the phase difference as a function of QPO frequency. For Type-B QPOs, we find that the phase difference stays constant with frequency for most sources. We propose a simple jet precession model to explain these constant Type-B QPO phase differences. The phase difference of the Type-C QPO is not constant but systematically evolves with QPO frequency, with the resulting relation being similar for a number of high inclination sources, but more variable for low-inclination sources. We discuss how the evolving phase difference can naturally arise in the framework of precession models for the Type-C QPO, by considering the contributions of a direct and reflected component to the QPO waveform.

Ionisation state of the accretion disc in the neutron-star low-mass X-ray binary 4U 1728−34

We analysed an XMM-Newton plus a simultaneous Rossi X-ray Timing Explorer observation and a separate Suzaku observation of the neutron-star low-mass X-ray binary 4U 1728-34. We fitted the X-ray spectra with the self-consistent reflection model relxill. We found that the inclination angle of 4U 1728-34 is 49 degrees, consistent with the upper limit of 60 degrees deduced from the absence of eclipses or dips in this source. The inclination angle in the fit of the XMM-Newton/RXTE observation is larger than 85 degrees, which may be due to the possible calibration issues of the PN instrument in timing mode. We also found that the thermal emission from the accretion disc is not significant. This could be explained either by the relatively high column density of the interstellar medium along the line of sight to the source, which decreases the number of soft disc photons, or if most of the soft thermal photons from the disc are reprocessed in the corona. The ionisation parameter derived from the fits is larger than the value predicted in the framework of the standard reflection model, wherein the disc is irradiated by an X-ray source above the compact object. This inconsistency suggests that irradiation from the neutron star and its boundary layer may play an important role in the ionisation of the accretion disc, and hence in the reflection component in this source.

A model-independent comparison of the variability of accreting neutron stars and black holes

We use Rossi X-ray Timing Explorer observations to conduct a population study of the timing properties of accretion-powered neutron star Low Mass X-ray Binaries (LMXBs), following a model-independent approach that was originally applied to black hole systems. The ratios of integrated power in four equally spaced Fourier frequency bands allow power spectral shapes to be parameterised with two `power colour' values, providing a simple way of tracking the evolution in timing properties across observations. We show that neutron star LMXBs follow a remarkably similar power spectral evolution to black hole LMXBs, confirming that the broadband noise variability seen in both types of system has a common origin in the accretion flow. Both neutron stars and black holes follow a similar clear track in the power colour-colour diagram as they transition from the hard through soft states. Quantifying the position on this oval track using a single parameter, the power-spectral `hue', we find that the transition in X-ray spectral shape occurs at the same hue for both neutron star and black hole systems. The power colours of Z sources map on to those of soft state black holes, while those of atoll sources correspond to all spectral states. There is no clear dependence of power colour on neutron star spin, or on whether the neutron star is clearly magnetised (determined by ms X-ray pulsations).

On the 2018 Outburst of the Accreting Millisecond X-Ray Pulsar Swift J1756.9–2508 As Seen with NICER

Abstract We report on the coherent timing analysis of the 182 Hz accreting millisecond X-ray pulsar Swift J1756.9−2508 during its 2018 outburst as observed with the Neutron Star Interior Composition Explorer (NICER). Combining our NICER observations with Rossi X-ray Timing Explorer observations of the 2007 and 2009 outbursts, we also studied the long-term spin and orbital evolution of this source. We find that the binary system is well described by a constant orbital period model, with an upper limit on the orbital period derivative of ∣ P ˙ b ∣ < 7.4 × 10 − 13 s s−1. Additionally, we improve upon the source coordinates through astrometric analysis of the pulse arrival times, finding R.A. = 17h56m57.ˢ18 ± 0.ˢ08 and decl. = −25°06′27.″8 ± 3.″5, while simultaneously measuring the long-term spin frequency derivative as ν ˙ = − 7.3 × 10 − 16 Hz s−1. We briefly discuss the implications of these measurements in the context of the wider population of accreting millisecond pulsars.