We study the evolution of newborn neutron stars in high-mass X-ray binaries interacting with a wind-fed super-Eddington disk. The inner disk is regularized to a radiation-dominated quasi-spherical configuration for which we calculate the inner radius of the disk, the total luminosity of the system and the torque acting on the neutron star accordingly, following the evolution of the system through the ejector and early propeller stages. We find that the systems with $B \gtrsim 10^{13}$ G pass through a short ($\sim 20\,{\rm yr}$) ejector stage appearing as supernova impostors followed by a propeller stage lasting $\sim 10^3\,{\rm yr}$. In the super-critical propeller stage the system is still bright ($L\sim 10^{40}\,{\rm erg\, s^{-1}}$) due to the spindown power and therefore appears as an ultra-luminous X-ray source (ULX). The system evolves into pulsating ULX (PULX) when the neutron star spins down to a period ($P\sim 1$ s) allowing for accretion onto its surface to commence. Systems with lower magnetic fields, $B \sim 10^{11}$ G, pass through a long ($10^5\,{\rm yr}$) super-critical propeller stage with luminosities similar to those of the ultra-luminous super-soft sources (ULS), $L \lesssim 10^{40}\,{\rm erg\, s^{-1}}$. The equilibrium periods of these systems in the accretion stage is about $10\,{\rm ms}$, which is much smaller than the typical period range of PULX observed to date. Such systems could have a larger population, but their pulsations would be elusive due to the smaller size of the magnetosphere. Our results suggest that the ULS and some nonpulsating ULX are rapidly spinning and highly magnetized young neutron stars at the super-critical propeller stage.

Ultra-luminous X-ray sources (ULXs) are off-nuclear point sources in nearby galaxies with luminosities well exceeding the Eddington limit for stellar-mass objects. It has been recognized after the discovery of pulsating ULXs (PULXs) that a fraction of these sources could be accreting neutron stars in high-mass X-ray binaries (HMXBs) though the majority of ULXs are lacking in coherent pulsations. The earliest stage of some HMXBs may harbor rapidly rotating neutron stars propelling out the matter transferred by the massive companion. The spin-down power transferred by the neutron-star magnetosphere to the accretion disk at this stage can well exceed the Eddington luminosities and the system appears as a non-pulsating ULX. In this picture, PULXs appear as super-critical mass-accreting descendants of non-pulsating ULXs. We present this evolutionary scenario within a self-consistent model of magnetosphere-disk interaction and discuss the implications of our results on the spin and magnetic field of the neutron star.

The neutron stars in low-mass X-ray binaries (LMXBs) are usually thought to be weakly magnetized objects accreting matter from their low-mass companions in the form of a disk. Albeit weak as compared to those in young neutron-star systems, the neutron-star magnetospheres in LMXBs can play an important role in determining the correlations between spectral and temporal properties. Parallel tracks appearing in the plane of kilohertz (kHz) quasi-periodic oscillation (QPO) frequency versus X-ray flux can be used as a tool to study the magnetosphere-disk interaction in neutron-star LMXBs. For dynamically important weak fields, the formation of a non-Keplerian magnetic boundary layer at the innermost disk truncated near the surface of the neutron star is highly likely. Such a boundary region may harbor oscillatory modes of frequencies in the kHz range. We generate parallel tracks using the boundary region model of kHz QPOs. We also present the direct application of our model to the reproduction of the observed parallel tracks of individual sources such as 4U 1608--52, 4U 1636--53, and Aql X-1. We reveal how the radial width of the boundary layer must vary in the long-term flux evolution of each source to regenerate the parallel tracks. The run of the radial width looks similar for different sources and can be fitted by a generic model function describing the average steady behavior of the boundary region in the long term. The parallel tracks then correspond to the possible quasi-steady states the source can occupy around the average trend.

Using the spin and flux evolution of the accreting pulsar 4U 1626$-$67 across the 2008 torque reversal, we determine the fastness parameter dependence of the dimensionless torque acting on the pulsar. We find that the dimensionless torque is qualitatively different from the existing models: it is concave-up across the torque equilibrium whereas the existing torque models predict a concave-down (convex) relation with the fastness parameter. We show that the dimensionless torque has a cubic dependence on the fastness parameter near the torque equilibrium. We also find that the torque can not attain large values away from the equilibrium, either in the positive or the negative side, but saturates at limited values. The spin-down torque can attain a 2.5 times larger magnitude at the saturation limit than the spin-up torque. From the evolution of the frequency of quasi-periodic oscillations of 4U 1626$-$67 across the torque reversal of 1990, we determine the critical fastness parameter corresponding to torque equilibrium to be $\omega_{\rm c} \simeq 0.75$ within the framework of the beat frequency model and boundary region model for reasonable values of the model parameters. We find that the disc magnetosphere interaction becomes unstable when the inner radius approaches the corotation radius as predicted by some models, though with a longer timescale. We also find that there is an unstable regime that is triggered when the fastness parameter is 0.8 times the critical fastness parameter ($\omega=0.6$ for $\omega_{\rm c} \simeq 0.75$) possibly associated with an instability observed in numerical simulations.

We study the dependence of kHz quasi-periodic oscillation (QPO) frequency on accretion-related parameters in the ensemble of neutron star low-mass X-ray binaries. Based on the mass accretion rate, $\dot{M}$, and the magnetic field strength, $B$, on the surface of the neutron star, we find a correlation between the lower kHz QPO frequency and $\dot{M}/B^{2}$. The correlation holds in the current ensemble of Z and atoll sources and therefore can explain the lack of correlation between the kHz QPO frequency and X-ray luminosity in the same ensemble. The average run of lower kHz QPO frequencies throughout the correlation can be described by a power-law fit to source data. The simple power-law, however, cannot describe the frequency distribution in an individual source. The model function fit to frequency data, on the other hand, can account for the observed distribution of lower kHz QPO frequencies in the case of individual sources as well as the ensemble of sources. The model function depends on the basic length scales such as the magnetospheric radius and the radial width of the boundary region, both of which are expected to vary with $\dot{M}$ to determine the QPO frequencies. In addition to modifying the length scales and hence the QPO frequencies, the variation in $\dot{M}$, being sufficiently large, may also lead to distinct accretion regimes, which would be characterized by Z and atoll phases.

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

The kHz QPOs observed in a neutron star low mass X-ray binary are likely to be produced in the innermost regions of accretion disk around the neutron star. The rotational dynamics of the inner disk can be characterized by the presence of either sub-Keplerian or super-Keplerian accretion flow depending on the relative fastness of the neutron spin as compared to the Keplerian frequency at the inner disk radius. Within the magnetosphere-disk interaction model, the frequency difference between the two kHz QPOs observed in a given source can be estimated to be slightly higher than or nearly around the neutron star spin frequency if the neutron star is a slow rotator and less than the stellar spin frequency if the neutron star is a fast rotator.

We apply the global mode analysis, which has been recently developed for the modeling of kHz quasi-periodic oscillations (QPOs) from neutron stars, to the inner region of an accretion disk around a rotating black hole. Within a pseudo-Newtonian approach that keeps the ratio of the radial epicyclic frequency $\kappa$ to the orbital frequency $\Omega$ the same as the corresponding ratio for a Kerr black hole we determine the innermost disk region where the hydrodynamic modes grow in amplitude. We find that the radiation flux emerging from the inner disk has the highest values within the same region. Using the flux weighted averages of the frequency bands over this region we identify the growing modes with highest frequency branches $\Omega +\kappa $ and $\Omega $ to be the plausible candidates for the high-frequency QPO pairs observed in black hole systems. The observed frequency ratio around 1.5 can therefore be understood naturally in terms of the global free oscillations in the innermost region of a viscous accretion disk around a black hole without invoking a particular resonance to produce black hole QPOs. Although the frequency ratio $<\Omega +\kappa>/<\Omega>$ is found to be not sensitive to the black hole's spin which is good for explaining the high-frequency QPOs it may work as a limited diagnostic of the spin parameter to distinguish black holes with very large spin from the slowly rotating ones. Within our model we estimate the frequency ratio of a high-frequency QPO pair to be greater than 1.5 if the black hole is a slow rotator. For fast rotating black holes, we expect the same ratio to be less than 1.5.

We show that the period clustering of anomalous X-ray pulsars (AXPs) and soft gamma-ray repeaters (SGRs), their X-ray luminosities, ages and statistics can be explained with fallback disks with large initial specific angular momentum. The disk evolution models are developed by comparison to self-similar analytical models. The initial disk mass and angular momentum set the viscous timescale. An efficient torque, with (1 - w^2) dependence on the fastness parameter w leads to period clustering in the observed AXP-SGR period range under a wide range of initial conditions. The timescale t_0 for the early evolution of the fallback disk, and the final stages of fallback disk evolution, when the disk becomes passive, are the crucial determinants of the evolution. The disk becomes passive at temperatures around 100 K, which provides a natural cutoff for the X-ray luminosity and defines the end of evolution in the observable AXP and SGR phase. This low value for the minimum temperature for active disk turbulence indicates that the fallback disks are active up to a large radius greater than ~10^{12} cm. We find that transient AXPs and SGRs are likely to be older than their persistent cousins. A fallback disk with mass transfer rates corresponding to the low quiescent X-ray luminosities of the transient sources in early evolutionary phases would have a relatively lower initial mass, such that the mass-flow rate in the disk is not sufficient for the inner disk to penetrate into the light cylinder of the young neutron star, making mass accretion onto the neutron star impossible. The transient AXP phase therefore must start later. The model results imply that the transient AXP/SGRs, although older, are likely to be similar in number to persistent sources (abridged).

The observational characteristics of quasi-periodic oscillations (QPOs) from accreting neutron stars strongly indicate the oscillatory modes in the innermost regions of accretion disks as a likely source of the QPOs. The inner regions of accretion disks around neutron stars can harbor very high frequency modes related to the radial epicyclic frequency $\kappa $. The degeneracy of $\kappa $ with the orbital frequency $\Omega $ is removed in a non-Keplerian boundary or transition zone near the magnetopause between the disk and the compact object. We show, by analyzing the global hydrodynamic modes of long wavelength in the boundary layers of viscous accretion disks, that the fastest growing mode frequencies are associated with frequency bands around $\kappa $ and $\kappa \pm \Omega $. The maximum growth rates are achieved near the radius where the orbital frequency $\Omega $ is maximum. The global hydrodynamic parameters such as the surface density profile and the radial drift velocity determine which modes of free oscillations will grow at a given particular radius in the boundary layer. In accordance with the peak separation between kHz QPOs observed in neutron-star sources, the difference frequency between two consecutive bands of the fastest growing modes is always related to the spin frequency of the neutron star. This is a natural outcome of the boundary condition imposed by the rotating magnetosphere on the boundary region of the inner disk.

We show that: (i) the long-term X-ray outburst light curve of the transient AXP XTE J1810-197 can be accounted for by a fallback disk that is evolving towards quiescence through a disk instability after having been heated by a soft gamma-ray burst, (ii) the spin-frequency evolution of this source in the same period can also be explained by the disk torque acting on the magnetosphere of the neutron star, (iii) most significantly, recently observed pulsed-radio emission from this source coincides with the epoch of minimum X-ray luminosity. This is natural in terms of a fallback disk model, as the accretion power becomes so low that it is not sufficient to suppress the beamed radio emission from XTE J1810-197.

The recent detection of the anomalous X-ray pulsar (AXP) 4U 0142+61 in the mid-infrared with the {\it Spitzer} observatory by Z.Wang and coworkers constitutes the first instance of a disk around an AXP. We show, by analyzing earlier optical and near-IR data together with the recent data, that the overall broadband data set can be reproduced by a single model of an irradiated and viscously heated disk.

The anomalous X-ray pulsar 4U 0142+61 was recently detected in the mid infrared bands with the SPITZER Observatory (Wang, Chakrabarty & Kaplan 2006). This observation is the first instance for a disk around an AXP. From a reanalysis of optical and infrared data, we show that the observations indicate that the disk is likely to be an active disk rather than a passive dust disk beyond the light cylinder, as proposed in the discovering paper. Furthermore, we show that the irradiated accretion disk model can also account for all the optical and infrared observations of the anomalous X-ray pulsars in the persistent state.

We investigate the rotational dynamics of disk accretion around a strongly magnetized neutron star with an aligned dipole field. The magnetospheric field is assumed to thread the disk plasma both inside and outside the corotation radius. As a result of disk-star interaction, the magnetic torque on the disk affects the structure of accretion flow to yield the observed spin-up or spin-down rates for a source of given fastness, magnetic field strength, and mass accretion rate. Within the model we obtain a prescription for the dynamical viscosity of such magnetically modified solutions for a Keplerian disk. We then use this prescription to find a model solution for the rotation rate profile throughout the entire disk including the non-Keplerian inner disk. We find that the non-Keplerian angular velocity transition region is not necessarily narrow for a source of given spin state. The boundary layer approximation, as in the standard magnetically threaded disk model, holds only in the case of dynamical viscosity decreasing all the way to the innermost edge of the disk. These results are applied to several observed disk-fed X-ray pulsars that have exhibited quasi-periodic oscillations (QPOs). The QPO frequencies provide a constraint on the fastness parameter and enable one to determine uniquely the width of the angular velocity transition zone for each source within model assumptions. We discuss the implications of these results on the value of the critical fastness parameter for a magnetized star in spin equilibrium. Applications of our model are also made with relevant parameters from recent numerical simulations of quasi-stationary disk-magnetized star interactions.