Accretion-powered Millisecond Pulsar Research Articles

From Feast to Famine: A Systematic Study of Accretion onto Oblique Pulsars with 3D GRMHD Simulations

Disk-fed accretion onto neutron stars can power a wide range of astrophysical sources ranging from X-ray binaries, to accretion-powered millisecond pulsars, ultraluminous X-ray sources, and gamma-ray bursts. A crucial parameter controlling the gas–magnetosphere interaction is the strength of the stellar dipole. In addition, coherent X-ray pulsations in many neutron star systems indicate that the star's dipole moment is oblique relative to its rotation axis. Therefore, it is critical to systematically explore the 2D parameter space of the star's magnetic field strength and obliquity, which is what this work does, for the first time, in the framework of 3D general-relativistic magnetohydrodynamics. If the accretion disk carries its own vertical magnetic field, this introduces an additional factor: the relative polarity of the disk and stellar magnetic fields. We find that depending on the strength of the stellar dipole and the star–disk relative polarity, the neutron star's jet power can either increase or decrease with increasing obliquity. For weak dipole strength (equivalently, high accretion rate), the parallel polarity results in a positive correlation between jet power and obliquity, whereas the antiparallel orientation displays the opposite trend. For stronger dipoles, the relative-polarity effect disappears, and jet power always decreases with increasing obliquity. The influence of the relative polarity gradually disappears as obliquity increases. Highly oblique pulsars tend to have an increased magnetospheric radius, a lower mass accretion rate, and enter the propeller regime at lower magnetic moments than aligned stars.

Polarized radiation from an accretion shock in accreting millisecond pulsars using exact Compton scattering formalism

Pulse profiles of accreting millisecond pulsars can be used to determine neutron star (NS) parameters, such as their masses and radii, and therefore provide constraints on the equation of state of cold dense matter. Information obtained by the Imaging X-ray Polarimetry Explorer (IXPE) can be used to decipher pulsar inclination and magnetic obliquity, providing ever tighter constraints on other parameters. In this paper, we develop a new emission model for accretion-powered millisecond pulsars based on thermal Comptonization in an accretion shock above the NS surface. The shock structure was approximated by an isothermal plane-parallel slab and the Stokes parameters of the emergent radiation were computed as a function of the zenith angle and energy for different values of the electron temperature, the Thomson optical depth of the slab, and the temperature of the seed blackbody photons. We show that our Compton scattering model leads to a significantly lower polarization degree of the emitted radiation compared to the previously used Thomson scattering model. We computed a large grid of shock models, which can be combined with pulse profile modeling techniques both with and without polarization included. In this work, we used the relativistic rotating vector model for the oblate NS in order to produce the observed Stokes parameters as a function of the pulsar phase. Furthermore, we simulated the data to be produced by IXPE and obtained constraints on model parameters using nested sampling. The developed methods can also be used in the analysis of the data from future satellites, such as the enhanced X-ray Timing and Polarimetry mission.

Neutron star parameter constraints for accretion-powered millisecond pulsars from the simulated IXPE data

We have simulated the X-ray polarization data that can be obtained with the Imaging X-ray Polarimetry Explorer, when observing accretion-powered millisecond pulsars. We estimated the necessary exposure times for SAX J1808.4−3658 in order to obtain different accuracy in the measured time-dependent Stokes profiles integrated over all energy channels. We found that the measured relative errors strongly depend on the relative configuration of the observer and the emitting hotspot. The improvement in the minimum relative error in Stokes Q and U parameters as a function of observing time t scales as 1/√t, and it spans the range from 30–90% with a 200 ks exposure time to 20–60% with a 500 ks exposure time (in the case of data binned in 19 phase bins). The simulated data were also used to predict how accurate measurements of the geometrical parameters of the neutron star can be made when modelling only Q and U parameters, but not the flux. We found that the observer inclination and the hotspot co-latitude could be determined with better than 10° accuracy for most of the cases we considered. In addition, we show that the position of a secondary hotspot can also be constrained when the spot is not obscured by an accretion disc. These measurements can be used to further constrain the neutron star mass and radius when combined with modelling of the X-ray pulse profile.

NS parameters from IXPE data of AMPs

AbstractIXPE is a NASA/ASI Small Explorer Mission. It will probe the X-ray polarization properties of celestial sources. In particular, for Accretion-powered Millisecond Pulsars (AMPs), IXPE can provide us with unique information on their geometry. These information together with pulse shape modelling will strongly boost the achievable sensitivity on measuring the AMPs mass and the radius. As a case-study, we simulated an observation of SAXJ1808.4-3658 and studied the accuracy that can be achieved in the measured time-dependent Stokes profiles. From these data we estimated how well IXPE will be able to constraint NS geometrical parameters, such as the inclination and hot-spot co-latitude angles.

Optical pulsations from a transitional millisecond pulsar

Weakly magnetic, millisecond spinning neutron stars attain their very fast rotation through a 1E8-1E9 yr long phase during which they undergo disk-accretion of matter from a low mass companion star. They can be detected as accretion-powered millisecond X-ray pulsars if towards the end of this phase their magnetic field is still strong enough to channel the accreting matter towards the magnetic poles. When mass transfer is much reduced or ceases altogether, pulsed emission generated by particle acceleration in the magnetosphere and powered by the rotation of the neutron star is observed, preferentially in the radio and gamma-ray bands. A few transitional millisecond pulsars that swing between an accretion-powered X-ray pulsar regime and a rotationally-powered radio pulsar regime in response to variations of the mass in-flow rate have been recently identified. Here we report the detection of optical pulsations from a transitional pulsar, the first ever from a millisecond spinning neutron star. The pulsations were observed when the pulsar was surrounded by an accretion disk and originated inside the magnetosphere or within a few hundreds of kilometres from it. Energy arguments rule out reprocessing of accretion-powered X-ray emission and argue against a process related to accretion onto the pulsar polar caps; synchrotron emission of electrons in a rotation-powered pulsar magnetosphere seems more likely.

Searching for X-ray Pulsations from Neutron Stars Using NICER

AbstractThe Neutron Star Interior Composition Explorer (NICER) presents an exciting new capability for exploring the modulation properties of X-ray emitting neutron stars, including large area, low background, extremely precise absolute event time stamps, superb low-energy response and flexible scheduling. The Pulsation Searches and Multiwavelength Coordination working group has designed a 2.5 Ms observing program to search for emission and characterize the modulation properties of about 30 known or suspected neutron star sources across a number of source categories. A key early goal will be to search for pulsations from millisecond pulsars that might exhibit thermal pulsations from the surface suitable for pulse profile modeling to constrain the neutron star equation of state. In addition, we will search for pulsations from transitional millisecond pulsars, isolated neutron stars, low-mass X-ray binaries (LMXBs), accretion-powered millisecond pulsars, central compact objects and other sources. We present our science plan and initial results from the first months of the NICER mission, including the discovery of pulsations from the millisecond pulsar J1231–1411.

Observational constraints on neutron star masses and radii

Precise and reliable measurements of the masses and radii of neutron stars with a variety of masses would provide valuable guidance for improving models of the properties of cold matter with densities above the saturation density of nuclear matter. Several different approaches for measuring the masses and radii of neutron stars have been tried or proposed, including analyzing the X-ray fluxes and spectra of the emission from neutron stars in quiescent low-mass X-ray binary systems and thermonuclear burst sources; fitting the energy-dependent X-ray waveforms of rotation-powered millisecond pulsars, burst oscillations with millisecond periods, and accretion-powered millisecond pulsars; and modeling the gravitational radiation waveforms of coalescing double neutron star and neutron star -- black hole binary systems. We describe the strengths and weaknesses of these approaches, most of which currently have substantial systematic errors, and discuss the prospects for decreasing the systematic errors in each method.

MULTIBAND STUDIES OF THE OPTICAL PERIODIC MODULATION IN THE X-RAY BINARY SAX J1808.4–3658 DURING ITS QUIESCENCE AND 2008 OUTBURST

We report on time-resolved optical imaging of the X-ray binary SAX J1808.4-3658 during its quiescent state and 2008 outburst. The binary, containing an accretion-powered millisecond pulsar, has a large sinusoidal-like modulation in its quiescent optical emission. We employ a Markov chain Monte Carlo technique to fit our multi-band light curve data in quiescence with an irradiated star model, and derive a tight constraint of $50^{+6}_{-5}$ deg on the inclination angle $i$ of the binary system. The pulsar and its companion are constrained to have masses of $0.97^{+0.31}_{-0.22} M_{\sun}$ and $0.04^{+0.02}_{-0.01} M_{\sun}$ (both 1$\sigma$ ranges), respectively. The dependence of these results on the measurements of the companion's projected radial velocity is discussed. We also find that the accretion disk had nearly constant optical fluxes over a $\sim$500 day period in the quiescent state our data covered, but started brightening 1.5 months before the 2008 outburst. Variations in modulation during the outburst were detected in our four observations made 7-12 days after the start of the outburst, and a sinusoidal-like modulation with 0.2 mag amplitude changed to have a smaller amplitude of 0.1 mag. The modulation variations are discussed. We estimate the albedo of the companion during its quiescence and the outburst, which was approximately 0 and 0.8 (for isotropic emission), respectively. This large difference probably provides additional evidence that the neutron star in the binary turns on as a radio pulsar in quiescence.

Chandra X-ray and Gemini near-infrared observations of the eclipsing millisecond pulsar SWIFT J1749.4−2807 in quiescence

We report on Chandra X-ray and Gemini-North near-infrared K-band observations of the eclipsing accretion-powered millisecond X-ray pulsar SWIFT J1749.4-2807 in quiescence. Using the Chandra observation we derive a source position of Right Ascencion: 17:49:31.73 and Declination:-28:08:05.08. The position is accurate to 0.6" (90 per cent confidence). We find one source at a magnitude K=18.44+-0.03 with a position fully consistent with the accurate Chandra X-ray localization and a second source at K=19.2+-0.1 that falls close to the edge of the error circle in the deep K-band images. The presence of a few weaker sources as suggested by previous H-band observations presented in the literature cannot the ruled out. There is marginal evidence that the brighter of the these two sources is variable. Follow-up spectroscopy of this potential counterpart will show if this source is the true counterpart to SWIFT J1749.4-2807. If so, baring the presence of complicating effects such as heating of the mass-donor star, it would allow for the mass of the neutron star to be measured through the measurement of periodic radial velocity variations.

The 2009 outburst of accreting millisecond pulsar IGR J17511-3057 as observed by Swift and RXTE

The twelfth accretion-powered millisecond pulsar, IGR J17511-3057, was discovered in September 2009. In this work we study its spectral and timing properties during the 2009 outburst based on Swift and RXTE data. Our spectral analysis of the source indicates only slight spectral shape evolution during the entire outburst. The equivalent width of the iron line and the apparent area of the blackbody emission associated with the hotspot at the stellar surface both decrease significantly during the outburst. This is consistent with a gradual receding of the accretion disc as the accretion rate drops. The pulse profile analysis shows absence of dramatic shape evolution with a moderate decrease in pulse amplitude. This behaviour might result from a movement of the accretion column footprint towards the magnetic pole as the disc retreats. The time lag between the soft and the hard energy pulses increase by a factor of two during the outburst. A physical displacement of the centroid of the accretion shock relative to the blackbody spot or changes in the emissivity pattern of the Comptonization component related to the variations of the accretion column structure could cause this evolution. We have found that IGR J17511-3057 demonstrates outburst stages similar to those seen in SAX J1808.4-3658. A transition from the "slow decay" into the "rapid drop" stage, associated with the dramatic flux decrease, is also accompanied by a pulse phase shift which could result from an appearance of the secondary spot due to the increasing inner disc radius.

Timing of the accreting millisecond pulsar IGR J17511-3057

Context. Timing analysis of accretion-powered millisecond pulsars (AMPs) is a powerful tool for probing the physics of compact objects. The recently discovered IGR J17511-3057 was the twelfth discovered of the 13 AMPs known. The Rossi XTE satellite provided an extensive coverage of the 25 days-long observation of the source outburst.

DISCOVERY OF AN ACCRETING MILLISECOND PULSAR IN THE ECLIPSING BINARY SYSTEM SWIFT J1749.4–2807

We report on the discovery and the timing analysis of the first eclipsing accretion-powered millisecond X-ray pulsar (AMXP): SWIFT J1749.4–2807. The neutron star rotates at a frequency of ∼517.9 Hz and is in a binary system with an orbital period of 8.8 hr and a projected semimajor axis of ∼1.90 lt-s. Assuming a neutron star between 0.8 and 2.2 M☉ and using the mass function of the system and the eclipse half-angle, we constrain the mass of the companion and the inclination of the system to be in the ∼0.46–0.81 M☉ and ∼ 744–773 range, respectively. To date, this is the tightest constraint on the orbital inclination of any AMXP. As in other AMXPs, the pulse profile shows harmonic content up to the third overtone. However, this is the first AMXP to show a first overtone with rms amplitudes between ∼6% and ∼23%, which is the strongest ever seen and which can be more than two times stronger than the fundamental. The fact that SWIFT J1749.4–2807 is an eclipsing system that shows uncommonly strong harmonic content suggests that it might be the best source to date to set constraints on neutron star properties including compactness and geometry.

ON THE NATURE OF PULSE PROFILE VARIATIONS AND TIMING NOISE IN ACCRETING MILLISECOND PULSARS

Timing noise in the data on accretion-powered millisecond pulsars (AMP) appears as irregular pulse phase jumps on timescales from hours to weeks. A large systematic phase drift is also observed in the first discovered AMP SAX J1808.4-3658. To study the origin of these timing features, we use here the data of the well studied 2002 outburst of SAX J1808.4-3658. We develop first a model for pulse profile formation accounting for the screening of the antipodal emitting spot by the accretion disk. We demonstrate that the variations of the visibility of the antipodal spot associated with the receding accretion disk cause a systematic shift in Fourier phases, observed together with the changes in the pulse form. We show that a strong secondary maximum can be observed only in a narrow intervals of inner disk radii, which explains the very short appearance of the double-peaked profiles in SAX J1808.4-3658. By directly fitting the pulse profile shapes with our model, we find that the main parameters of the emitting spot such as its mean latitude and longitude as well as the emissivity pattern change irregularly causing small shifts in pulse phase, and the strong profile variations are caused by the increasing inner disk radius. We finally notice that significant variations in the pulse profiles in the 2002 and 2008 outbursts of SAX J1808.4-3658 happen at fluxes differing by a factor of 2, which can be explained if the inner disk radius is not a simple function of the accretion rate, but depends on the previous history.

A MODEL FOR THE WAVEFORM BEHAVIOR OF ACCRETING MILLISECOND X-RAY PULSARS: NEARLY ALIGNED MAGNETIC FIELDS AND MOVING EMISSION REGIONS

We investigate further a model of the accreting millisecond X-ray pulsars we proposed earlier. In this model, the X-ray-emitting regions of these pulsars are near their spin axes but move. This is to be expected if the magnetic poles of these stars are close to their spin axes, so that accreting gas is channeled there. As the accretion rate and the structure of the inner disk vary, gas is channeled along different field lines to different locations on the stellar surface, causing the X-ray-emitting areas to move. We show that this 'nearly aligned moving spot model' can explain many properties of the accreting millisecond X-ray pulsars, including their generally low oscillation amplitudes and nearly sinusoidal waveforms; the variability of their pulse amplitudes, shapes, and phases; the correlations in this variability; and the similarity of the accretion- and nuclear-powered pulse shapes and phases in some. It may also explain why accretion-powered millisecond pulsars are difficult to detect, why some are intermittent, and why all detected so far are transients. This model can be tested by comparing with observations the waveform changes it predicts, including the changes with accretion rate.

ORIGIN OF INTERMITTENT ACCRETION-POWERED X-RAY OSCILLATIONS IN NEUTRON STARS WITH MILLISECOND SPIN PERIODS

We have shown previously that many of the properties of persistent accretion-powered millisecond pulsars can be understood if their X-ray emitting areas are near their spin axes and move as the accretion rate and structure of the inner disk vary. Here we show that this "nearly aligned moving spot model" may also explain the intermittent accretion-powered pulsations that have been detected in three weakly magnetic accreting neutron stars. We show that movement of the emitting area from very close to the spin axis to about 10 degrees away can increase the fractional rms amplitude from less than about 0.5 percent, which is usually undetectable with current instruments, to a few percent, which is easily detectable. The second harmonic of the spin frequency usually would not be detected, in agreement with observations. The model produces intermittently detectable oscillations for a range of emitting area sizes and beaming patterns, stellar masses and radii, and viewing directions. Intermittent oscillations are more likely in stars that are more compact. In addition to explaining the sudden appearance of accretion-powered millisecond oscillations in some neutron stars with millisecond spin periods, the model explains why accretion-powered millisecond oscillations are relatively rare and predicts that the persistent accretion-powered millisecond oscillations of other stars may become undetectable for brief intervals. It suggests why millisecond oscillations are frequently detected during the X-ray bursts of some neutron stars but not others and suggests mechanisms that could explain the occasional temporal association of intermittent accretion-powered oscillations with thermonuclear X-ray bursts.

ERRATUM: “DISCOVERY OF THE ACCRETION-POWERED MILLISECOND PULSAR SWIFT J1756.9-2508 WITH A LOW-MASS COMPANION” (2007, ApJ, 668, L147)

We report on the discovery by the Swift Gamma-Ray Burst Explorer of the eighth known transient accretion-powered millisecond pulsar, SWIFT J1756.9-2508, as part of routine observations with the Swift Burst Alert Telescope hard X-ray transient monitor. The pulsar was subsequently observed by both the X-Ray Telescope on Swift and the Rossi X-Ray Timing Explorer Proportional Counter Array. It has a spin frequency of 182 Hz (5.5 ms) and an orbital period of 54.7 minutes. The minimum companion mass is between 0.0067 and 0.0086 solar masses, depending on the mass of the neutron star, and the upper limit on the mass is 0.030 solar masses (95% confidence level). Such a low mass is inconsistent with brown dwarf models, and comparison with white dwarf models suggests that the companion is a He-dominated donor whose thermal cooling has been at least modestly slowed by irradiation from the accretion flux. No X-ray bursts, dips, eclipses or quasi-periodic oscillations were detected. The current outburst lasted approximately 13 days and no earlier outbursts were found in archival data.

DISCOVERY OF BURST OSCILLATIONS IN THE INTERMITTENT ACCRETION-POWERED MILLISECOND PULSAR HETE J1900.1-2455

We report the discovery of burst oscillations from the intermittent accretion-powered millisecond pulsar (AMP) HETE J1900.1-2455, with a frequency approximately 1 Hz below the known spin frequency. The burst oscillation properties are far more similar to those of the non-AMPs and Aql X-1 (an intermittent AMP with a far lower duty cycle), than those of the AMPs SAX J1808.4-3658 and XTE J1814-338. We discuss the implications for models of the burst oscillation and intermittency mechanisms.

THE LUMINOSITY AND ENERGY DEPENDENCE OF PULSE PHASE LAGS IN THE ACCRETION-POWERED MILLISECOND PULSAR SAX J1808.4–3658

Soft phase lags, in which X-ray pulses in lower energy bands arrive later than pulses in higher energy bands, have been observed in nearly all accretion-powered millisecond pulsars, but their origin remains an open question. In a study of the 2.5 ms accretion-powered pulsar SAX J1808.4-3658, we report that the magnitude of these lags is strongly dependent on the accretion rate. During the brightest stage of the outbursts from this source, the lags increase in magnitude as the accretion rate drops; when the outbursts enter their dimmer flaring-tail stage, the relationship reverses. We evaluate this complex dependence in the context of two theoretical models for the lags, one relying on the scattering of photons by the accretion disk and the other invoking a two-component model for the photon emission. In both cases, the turnover suggests that we are observing the source transitioning into the "propeller" accretion regime.

AN ACCURATE DETERMINATION OF THE OPTICAL PERIODIC MODULATION IN THE X-RAY BINARY SAX J1808.4–3658

We report on optical imaging of the X-ray binary SAX J1808.4-3658 with the 8-m Gemini South Telescope. The binary, containing an accretion-powered millisecond pulsar, appears to have a large periodic modulation in its quiescent optical emission. In order to clarify the origin of this modulation, we obtained three time-resolved $r'$-band light curves (LCs) of the source in five days. The LCs can be described by a sinusoid, and the long time-span between them allows us to determine optical period P=7251.9 s and phase 0.671 at MJD 54599.0 (TDB; phase 0.0 corresponds to the ascending node of the pulsar orbit), with uncertainties of 2.8 s and 0.008 (90 % confidence), respectively. This periodicity is highly consistent with the X-ray orbital ephemeris. By considering this consistency and the sinusoidal shape of the LCs, we rule out the possibility of the modulation arising from the accretion disk. Our study supports the previous suggestion that the X-ray pulsar becomes rotationally powered in quiescence, with its energy output irradiating the companion star, causing the optical modulation. While it has also been suggested that the accretion disk would be evaporated by the pulsar, we argue that the disk exists and gives rise to the persistent optical emission. The existence of the disk can be verified by long-term, multi-wavelength optical monitoring of the source in quiescence, as an increasing flux and spectral changes from the source would be expected based on the standard disk instability model.

XMM-Newton Observations of IGR J00291+5934: Signs of a Thermal Spectral Component during Quiescence

We present X-ray observations of the transient accretion-powered millisecond pulsar IGR J00291+5934 during quiescence. IGR J00291+5934 is the first source among accretion-powered millisecond pulsars to show signs of a thermal component in its quiescent spectrum. Fitting this component with a neutron star atmosphere or a blackbody model we obtain soft temperatures (~64 and ~110 eV, respectively). As in other sources of this class a hard spectral component is also present, comprising more than 60% of the unabsorbed 0.5-10 keV flux. Interpreting the soft component as cooling emission from the neutron star, we can conclude that the compact object can be spun up to millisecond periods by accreting only 0.2 M☉.