Spin evolution of neutron stars in binary systems

I review our understanding of the evolution of the spin periods of neutron stars in binary stellar systems, from their birth as fast, spin-powered pulsars, through their middle life as accretion-powered pulsars, upto their recycling or “rebirth” as spin-powered pulsars with relatively low magnetic fields and fast rotation. I discuss how the new-born neutron star is spun down by electromagnetic and “propeller” torques, until accretion of matter from the companion star begins, and the neutron star becomes an accretion-powered X-ray pulsar. Detailed observations of massive radio pulsar binaries like PSR 1259-63 will yield valuable information about this phase of initial spindown. I indicate how the spin of the neutron star then evolves under accretion torques during the subsequent phase as an accretion-powered pulsar. Finally, I describe how the neutron star is spun up to short periods again during the subsequent phase of recycling, with the accompanying reduction in the stellar magnetic field, the origins of which are still not completely understood.     

The neutron stars of Soft X–ray Transients

Summary. Soft X–ray Transients (SXRTs) have long been suspected to contain old, weakly magnetic neutron stars that have been spun up by accretion torques. After reviewing their observational properties, we analyse the different regimes that likely characterise the neutron stars in these systems across the very large range of mass inflow rates, from the peak of the outbursts to the quiescent emission. While it is clear that close to the outburst maxima accretion onto the neutron star surface takes place, as the mass inflow rate decreases, accretion might stop at the magnetospheric boundary because of the centrifugal barrier provided by the neutron star. For low enough mass inflow rates (and sufficiently short rotation periods), the radio pulsar mechanism might turn on and sweep the inflowing matter away. The origin of the quiescent emission, observed in a number of SXRTs at a level of , plays a crucial role in constraining the neutron star magnetic field and spin period. Accretion onto the neutron star surface is an unlikely mechanism for the quiescent emission of SXRTs, as it requires very low magnetic fields and/or long spin periods. Thermal radiation from a cooling neutron star surface in between the outbursts can be ruled out as the only cause of the quiescent emission. We find that accretion onto the neutron star magnetosphere and shock emission powered by an enshrouded radio pulsar provide far more plausible models. In the latter case the range of allowed neutron star spin periods and magnetic fields is consistent with the values recently inferred from the properties of kHz quasi-periodic oscillation in low mass X–ray binaries. If quiescent SXRTs contain enshrouded radio pulsars, they provide a missing link between X–ray binaries and millisecond pulsars. Received 4 November 1997; Accepted 15 April 1998     

On the torque exerted by a warped,magnetically threaded accretion disk

Most astrophysical accretion disks are likely to be warped.In X-ray binaries,the spin evolution of an accreting neutron star is critically dependent on the interaction between the neutron star magnetic field and the accretion disk.There have been extensive investigations on the accretion torque exerted by a coplanar disk that is magnetically threaded by the magnetic field lines from the neutron stars,but relevant works on warped/tilted accretion disks are still lacking.In this paper we develop a simplified twocomponent model,in which the disk is comprised of an inner coplanar part and an outer,tilted part.Based on standard assumption on the formation and evolution of the toroidal magnetic field component,we derive the dimensionless torque and show that a warped/titled disk is more likely to spin up the neutron star compared with a coplanar disk.We also discuss the possible influence of various initial parameters on the torque.     

Accretion spin-up of low-magnetic neutron stars

The present paper is concerned with the spin-up of low-magnetic neutron stars by the accretion of matter onto the star. Calculations have been made for the evolution of the rotation of a neutron star and applied to different stellar models. It is shown that the existence of a millisecond pulsar imposes no restriction on any of the equations of state considered. However, constraints would arise with the possible discovery of third-octave pulsars (with frequencies in excess of 1000 Hz). Predictions are made as to the distribution of bursters over the orbital periods of neutron stars (about half of these having similar orbital periods). It is demonstrated that in the case of continued accretion onto a star, after it has acquired the critical angular frequency allowing no diviation from axial symmetry, specific accretion disks can be formed with a smooth transition into a star. The specific angular momentum is computed for a neutron star for the instant of the attainment of the Oppenheimer-Volkoff limit.     

Magnetic and spin evolution of neutron stars in close binaries

The evolution of neutron stars in close binary systems with a low-mass companion is considered, assuming the magnetic field to be confined within the solid crust. We adopt the standard scenario for the evolution in a close binary system, in which the neutron star passes through four evolutionary phases ('isolated pulsar'–'propeller'– accretion from the wind of a companion – accretion resulting from Roche-lobe overflow). Calculations have been performed for a great variety of parameters characterizing the properties of both the neutron star and the low-mass companion. We find that neutron stars with more or less standard magnetic field and spin period that are processed in low-mass binaries can evolve to low-field rapidly rotating pulsars. Even if the main-sequence life of a companion is as long as 10¹⁰ yr, the neutron star can maintain a relatively strong magnetic field to the end of the accretion phase. The model that is considered can account well for the origin of millisecond pulsars.     

Gamma-Ray Emission from Be/X-ray Binaries

Be/X-ray binaries are systems formed by a massive Be star and a magnetized neutron star, usually in an eccentric orbit. The Be star has strong equatorial winds occasionally forming a circumstellar disk. When the neutron star intersects the disk the accretion rate dramatically increases and a transient accretion disk can be formed around the compact object. This disk can last longer than a single orbit in the case of major outbursts. If the disk rotates faster than the neutron star, the Cheng-Ruderman mechanism can produce a current of relativistic protons that would impact onto the disk surface, producing gamma-rays from neutral pion decays and initiating electromagnetic cascades inside the disk. In this paper we present calculations of the evolution of the disk parameters during both major and minor X-ray events, and we discuss the generation of gamma-ray emission at different energies within a variety of models that include both screened and unscreened disks.     

On the progenitors of core-collapse supernovae

Theory holds that a star born with an initial mass between about 8 and 140 times the mass of the Sun will end its life through the catastrophic gravitational collapse of its iron core to a neutron star or black hole. This core collapse process is thought to usually be accompanied by the ejection of the star’s envelope as a supernova. This established theory is now being tested observationally, with over three dozen core-collapse supernovae having had the properties of their progenitor stars directly measured through the examination of high-resolution images taken prior to the explosion. Here I review what has been learned from these studies and briefly examine the potential impact on stellar evolution theory, the existence of “failed supernovae”, and our understanding of the core-collapse explosion mechanism.     

Disk luminosity and angular momentum for accreting,weak field neutron stars in the ‘Slow’ rotation approximation

For accretion on to neutron stars possessing weak surface magnetic fields and substantial rotation rates (corresponding to the secular instability limit), we calculate the disk and surface layer luminosities general relativistically using the Hartle & Thorne formalism, and illustrate these quantities for a set of representative neutron star equations of state. We also discuss the related problem of the angular momentum evolution of such neutron stars and give a quantitative estimate for this accretion driven change in angular momentum. Rotation always increases the disk luminosity and reduces the rate of angular momentum evolution. These effects have relevance for observations of low-mass X-ray binaries.     

Analysing the atolls: X-ray spectral transitions of accreting neutron stars

We systematically analyse all the available X-ray spectra of disc accreting neutron stars (atolls and millisecond pulsars) from the RXTE data base. We show that while all these have similar spectral evolution as a function of mass accretion rate, there are also subtle differences. There are two different types of hard/soft transition, those where the spectrum softens at all energies, leading to a diagonal track on a colour–colour diagram, and those where only the higher energy spectrum softens, giving a vertical track. The luminosity at which the transition occurs is correlated with this spectral behaviour, with the vertical transition at   L / L _Edd∼ 0.02  while the diagonal one is at ∼0.1. Superimposed on this is the well-known hysteresis effect, but we show that classic, large-scale hysteresis occurs only in the outbursting sources, indicating that its origin is in the dramatic rate of change of mass accretion rate during the disc instability. We show that the long-term mass accretion rate correlates with the transition behaviour, and speculate that this is due to the magnetic field being able to emerge from the neutron star surface for low average mass accretion rates. While this is not strong enough to collimate the flow except in the millisecond pulsars, its presence may affect the inner accretion flow by changing the properties of the jet.     

Interpretation of the light curve for the young binary (?) system UY Aur

Using visual, photographic, and photoelectric measurements, we have constructed a historical light curve for the young binary system UY Aur on an interval longer than 100 yr. About a quarter of all magnitude estimates have been obtained for the first time from photographic plates of the Sternberg Astronomical Institute and Harvard College Observatory Astronomical Plate Stacks. Analysis of the light curve and the magnitude dependences of the polarization and color has led us to the following conclusions. Cyclic variations in the seasonally mean brightness of the binary’s primary component UY Aur A with a period of ≃16.3 yr occurred from the mid-1920s to the mid-1940s and after 1986. The variations are caused by the change in the rate of disk accretion onto the star attributable to the motion of the hypothetical companion UY Aur C around the primary star in an orbit with a semimajor axis of ≃ 6 AU. From the early 1950s to the mid-1980s, the periodicity of the seasonally mean variations was not noticeable due to nonperiodic eclipses of UY Aur A by gas-dust clouds. Between 1945 and 1974, another gas-dust cloud obscured and still obscures the component UY Aur B, causing its mean optical brightness to drop by several magnitudes. The role of the clouds that caused an almost simultaneous eclipse of the stars, whose separation in projection onto the celestial sphere exceeds 100 AU, is played by the denser and puffed-up regions of their accretion disks. These regions are the result of a dynamical interaction between the binary’s stars and the outer circumbinary accretion disk. The extinction variations with time are attributable to orbital motion of the binary’s stars and azimuthal inhomogeneity of the clump regions in the disks. A number of observational tests are suggested to verify our conclusions.     

A possible detection mechanism for pre-low-mass X-ray binaries

This work presents a possible detection mechanism for close, detached, neutron star–red dwarf binaries, which are expected to be the evolutionary precursors of low-mass X-ray binaries (LMXBs). Although this pre-low-mass X-ray binary (pre-LMXB) phase of evolution is predicted theoretically, as yet no such systems have been identified observationally. The calculations presented here suggest that the X-ray luminosity of neutron star wind accretion in a pre-LMXB system can be expected to exceed the intrinsic X-ray luminosity of the red dwarf secondary star. Furthermore, the temperature of the radiation emitted from the neutron star wind accretion process is expected, within the confines of a reasonable set of conditions, to lie within the detection range of X-ray satellites. Sources with X-ray luminosities greater than that expected for a red dwarf star, but the positions of which coincide with that of a red dwarf star, are then candidate pre-LMXB systems. These candidate systems should be surveyed for the radial velocity shifts that would occur as a result of the orbital motion of a red dwarf star within a close binary system containing a high-mass compact object.     

Flux-Vortex Pinning and Neutron Star Evolution

G. Srinivasan et al. (1990) proposed a simple and elegant explanation for the reduction of the neutron star magnetic dipole moment during binary evolution leading to low mass X-ray binaries and eventually to millisecond pulsars: Quantized vortex lines in the neutron star core superfluid will pin against the quantized flux lines of the proton superconductor. As the neutron star spins down in the wind accretion phase of binary evolution, outward motion of vortex lines will reduce the dipole magnetic moment in proportion to the rotation rate. The presence of a toroidal array of flux lines makes this mechanism inevitable and independent of the angle between the rotation and magnetic axes. The incompressibility of the flux-line array (Abrikosov lattice) determines the epoch when the mechanism will be effective throughout the neutron star. Flux vortex pinning will not be effective during the initial young radio pulsar phase. It will, however, be effective and reduce the dipole moment in proportion with the rotation rate during the epoch of spindown by wind accretion as proposed by Srinivasan et al. The mechanism operates also in the presence of vortex creep.     

Ohm’s law for plasma in general relativity and Cowling’s theorem

The general-relativistic Ohm’s law for a two-component plasma which includes the gravitomagnetic force terms even in the case of quasi-neutrality has been derived. The equations that describe the electromagnetic processes in a plasma surrounding a neutron star are obtained by using the general relativistic form of Maxwell equations in a geometry of slow rotating gravitational object. In addition to the general-relativistic effect first discussed by Khanna and Camenzind (Astron. Astrophys. 307:665, 1996) we predict a mechanism of the generation of azimuthal current under the general relativistic effect of dragging of inertial frames on radial current in a plasma around neutron star. The azimuthal current being proportional to the angular velocity ω of the dragging of inertial frames can give valuable contribution on the evolution of the stellar magnetic field if ω exceeds 2.7×10¹⁷(n/σ) s⁻¹ (n is the number density of the charged particles, σ is the conductivity of plasma). Thus in general relativity a rotating neutron star, embedded in plasma, can in principle generate axial-symmetric magnetic fields even in axisymmetry. However, classical Cowling’s antidynamo theorem, according to which a stationary axial-symmetric magnetic field can not be sustained against ohmic diffusion, has to be hold in the general-relativistic case for the typical plasma being responsible for the rotating neutron star.     

A Two—Temperature Supernova Fallback Disk Model for Anomalous X—ray Pulsars

We present a case study of the relevance of the radially pulsational instability of a two-temperature accretion disk around a neutron star to anomalous X-ray pulsars (AXPs). Our estimates are based on the approximation that such a neutron star disk with mass in the range of 10^-6-10^-5M⊙ is formed by supernova fallback. We derive several peculiar properties of the accretion disk instability: a narrow interval of X-ray pulse periods; lower X-ray luminosities; a period derivative and an evolution time scale. All these results are in good agreement with the observations of the AXPs.     

Evolution of massive binary stars in the LMC and its implications for radio pulsar population

The Hertzsprung-Russell diagram of the Large Magellanic Cloud compiled recently by Fitzpatrick & Garmany (1990) shows that there are a number of supergiant stars immediately redward of the main sequence although theoretical models of massive stars with normal hydrogen abundance predict that the region 4.5 ≤ logT _eff ≤ 4.3 should be un-populated (“gap”). Supergiants having surface enrichment of helium acquired for example from a previous phase of accretion from a binary companion, however, evolve in a way so that the evolved models and observed data are consistent — an observation first made by Tuchman & Wheeler (1990). We compare the available optical data on OB supergiants with computed evolutionary tracks of massive stars of metallicity relevant to the LMC with and without helium-enriched envelopes and conclude that a large fraction (≈ 60 per cent) of supergiant stars may occur in binaries. As these less evolved binaries will later evolve into massive X-ray binaries, the observed number and orbital period distribution of the latter can constrain the evolutionary scenarios of the supergiant binaries. The distributions of post main sequence binaries and closely related systems like WR + O stars are bimodal-consisting of close and wide binaries in which the latter type is numerically dominating. When the primary star explodes as a supernova leaving behind a neutron star, the system receives a kick and in some cases can lead to runaway O-stars. We calculate the expected space velocity distribution for these systems. After the second supernova explosion, the binaries in most cases, will be disrupted leading to two runaway neutron stars. In between the two explosions, the first born neutron star’s spin evolution will be affected by accretion of mass from the companion star. We determine the steady-state spin and radio luminosity distributions of single pulsars born from the massive stars under some simple assumptions. Due to their great distance, only the brightest radio pulsars may be detected in a flux-limited survey of the LMC. A small but significant number of observable single radio pulsars arising out of the disrupted massive binaries may appear in the short spin period range. Most pulsars will have a low velocity of ejection and therefore may cluster around the OB associations in the LMC.     

A new look at the origin of the 6.67 hr period X-ray pulsar 1E 161348-5055

The point X-ray source 1E 161348-5055 is observed to display pulsations with the period 6.67?hr and $|\dot{P}| \leq1.6 \times10^{-9}\,{\rm s\,s^{-1}}$ . It is associated with the supernova remnant RCW?103 and is widely believed to be a ～2000?yr old neutron star. Observations give no evidence for the star to be a member of a binary system. Nevertheless, it resembles an accretion-powered pulsar with the magnetospheric radius ～3000?km and the mass-accretion rate $\sim 10^{14}\,{\rm g\,s^{-1}}$ . This situation could be described in terms of accretion from a (residual) fossil disk established from the material falling back towards the star after its birth. However, current fall-back accretion scenarios encounter major difficulties explaining an extremely long spin period of the young neutron star. We show that the problems can be avoided if the accreting material is magnetized. The star in this case is surrounded by a fossil magnetic slab in which the material is confined by the magnetic field of the accretion flow itself. We find that the surface magnetic field of the neutron star within this scenario is ～10¹²?G and that a presence of $\gtrsim10^{-7}\,{\rm M_{\odot}}$ magnetic slab would be sufficient to explain the origin and current state of the pulsar.     

The evolution of the magnetic fields of neutron stars

Observational evidence, and theoretical models of the magnetic field evolution of neutron stars is discussed. Observational data indicates that the magnetic field of a neutron star decays significantly only if it has been a member of a close interacting binary. Theoretically, the magnetic field evolution has been related to the processing of a neutron star in a binary system through the spin evolution of the neutron star, and also through the accretion of matter on the neutron star surface. I describe two specific models, one in which magnetic flux is expelled from the superconducting core during spin-down, via a copuling between Abrikosov fluxoids and Onsager-Feynman vortices; and another in which the compression and heating of the stellar crust by the accreted mass drastically reduces the ohmic decay time scale of a magnetic field configuration confined entirely to the crust. General remarks about the behaviour of the crustal field under ohmic diffusion are also made.     

A model of two-stream non-radial accretion for binary X-ray pulsars

The general case of non-radial accretion is assumed to occur in real binary systems containing X-ray pulsars. The structure and the stability of the magnetosphere, the interaction between the magnetosphere and accreted matter, as well as evolution of neutron star in close binary system are examined within the framework of the two-stream model of nonradial accretion onto a magnetized neutron star. Observable parameters of X-ray pulsars are explained in terms of the model considered.     

Evolution of neutron stars in high-mass X-ray binaries

We consider the evolution of neutron stars during the X-ray phase of high-mass binaries. Calculations are performed assuming a crustal origin of the magnetic field. A strong wind from the companion can significantly influence the magnetic and spin behaviour of a neutron star even during the main-sequence life of the companion. In the course of evolution, the neutron star passes through four evolutionary phases ('isolated pulsar', propeller, wind accretion, and Roche lobe overflow). The model considered can naturally account for the observed magnetic fields and spin periods of neutron stars, as well as the existence of pulsating and non-pulsating X-ray sources in high-mass binaries. Calculations also predict the existence of a particular sort of high-mass binary with a secondary that fills its Roche lobe and a neutron star that does not accrete the overflowing matter because of fast spin.     

Forced oscillations in a hydrodynamical accretion disk and QPOs

This is the second of a series of papers aimed to look for an explanation on the generation of high frequency quasi-periodic oscillations (QPOs) in accretion disks around neutron star, black hole, and white dwarf binaries. The model is inspired by the general idea of a resonance mechanism in the accretion disk oscillations as was already pointed out by Abramowicz and Klu’zniak (2001). In a first paper (P'etri, 2005a, paper I), we showed that a rotating misaligned magnetic field of a neutron star gives rise to some resonances close to the inner edge of the accretion disk. In this second paper, we suggest that this process does also exist for an asymmetry in the gravitational potential of the compact object. We prove that the same physics applies, at least in the linear stage of the response to the disturbance in the system. This kind of asymmetry is well suited for neutron stars or white dwarfs possessing an inhomogeneous interior allowing for a deviation from a perfectly spherically symmetric gravitational field. After a discussion on the magnitude of this deformation applied to neutron stars, we show by a linear analysis that the disk initially in a cylindrically symmetric stationary state is subject to {three kinds of resonances: a corotation resonance, a Lindblad resonance due to a driven force and a parametric resonance}. In a second part, we focus on the linear response of a thin accretion disk in the 2D limit. {Waves are launched at the aforementioned resonance positions and propagate in some permitted regions inside the disk, according to the dispersion relation obtained by a WKB analysis}. In a last part, these results are confirmed and extended via non linear hydrodynamical numerical simulations performed with a pseudo-spectral code solving Euler's equations in a 2D cylindrical coordinate frame. {We found that for a weak potential perturbation, the Lindblad resonance is the only effective mechanism producing a significant density fluctuation}. In a last step, we replaced the Newtonian potential by the so called logarithmically modified pseudo-Newtonian potential in order to take into account some general-relativistic effects like the innermost stable circular orbit (ISCO). The latter potential is better suited to describe the close vicinity of a neutron star or a black hole. However, from a qualitative point of view, the resonance conditions remain the same. The highest kHz QPOs are then interpreted as the orbital frequency of the disk at locations where the response to the resonances are maximal. It is also found that strong gravity is not required to excite the resonances.