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.     

The bottom magnetic field and magnetosphere evolution of neutron star in low-mass X-ray binary

The accretion-induced neutron star (NS) magnetic field evolution is studied through considering the accretion flow to drag the field lines aside and dilute the polar-field strength, and as a result the equatorial field strength increases, which is buried inside the crust on account of the accretion-induced global compression of star crust. The main conclusions of model are as follows: (i) the polar field decays with increase in the accreted mass; (ii) the bottom magnetic field strength of about 10⁸ G can occur when the NS magnetosphere radius approaches the star radius, and it depends on the accretion rate as     ; and (iii) the NS magnetosphere radius decreases with accretion until it reaches the star radius, and its evolution is little influenced by the initial field and the accretion rate after accreting  ∼0.01 M_⊙  , which implies that the magnetosphere radii of NSs in low-mass X-ray binaries would be homogeneous if they accreted the comparable masses. As an extension, the physical effects of the possible strong magnetic zone in the X-ray NSs and recycled pulsars are discussed. Moreover, the strong magnetic fields in the binary pulsars PSR 1831−00 and PSR 1718−19 after accreting about  0.5 M_⊙  in the binary-accretion phase,  8.7 × 10¹⁰  and  1.28 × 10¹² G  , respectively, can be explained through considering the incomplete frozen flow in the polar zone. As an expectation of the model, the existence of the low magnetic field  (∼3 × 10⁷ G)  NSs or millisecond pulsars is suggested.     

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.     

Neutron Stars in X-ray Binaries and their Environments

Neutron stars in X-ray binary systems are fascinating objects that display a wide range of timing and spectral phenomena in the X-rays. Not only parameters of the neutron stars, like magnetic field strength and spin period evolve in their active binary phase, the neutron stars also affect the binary systems and their immediate surroundings in many ways. Here we discuss some aspects of the interactions of the neutron stars with their environments that are revelaed from their X-ray emission. We discuss some recent developments involving the process of accretion onto high magnetic field neutron stars: accretion stream structure and formation, shape of pulse profile and its changes with accretion torque. Various recent studies of reprocessing of X-rays in the accretion disk surface, vertical structures of the accretion disk and wind of companion star are also discussed here. The X-ray pulsars among the binary neutron stars provide excellent handle to make accurate measurement of the orbital parameters and thus also evolution of the binray orbits that take place over time scale of a fraction of a million years to tens of millions of years. The orbital period evolution of X-ray binaries have shown them to be rather complex systems. Orbital evolution of X-ray binaries can also be carried out from timing of the X-ray eclipses and there have been some surprising results in that direction, including orbital period glitches in two X-ray binaries and possible detection of the most massive circum-binary planet around a Low Mass X-ray Binary.     

Scenarios for the formation of binary and millisecond pulsars – A critical assessment

The evolution of high-and low-mass X-ray binaries (HMXB and LMXB) into different types of binary radio pulsars, the ‘high-mass binary pulsars’(HMBP) and ‘low-mass binary pulsars’ (LMBP) is discussed. The HMXB evolve either into Thorne-Zytkow objects or into short-period binaries consisting of a helium star plus a neutron star (or a black hole), resembling Cygnus X-3. The latter systems evolve (with or without a second common-envelope phase) into close binary pulsars, in which the companion of the pulsar may be a massive white dwarf, a neutron star or a black hole ( some final systems may also consist of two black holes). A considerable fraction of the systems may also be disrupted in the second supernova explosion. We discuss the possible reasons why the observed numbers of double neutron stars and of systems like Cyg X-3 are several orders of magnitude lower than theoretically predicted. It is argued that the observed systems form the tip of an iceberg of much larger populations of unobserved systems, some of which may become observable in the future. As to the LMBP, we consider in some detail the origins of systems with orbital periods in the range 1–20 days. We show that to explain their existence, losses of orbital angular momentum (e.g., by magnetic braking) and in a number of cases: also of mass, have to be taken into account. The masses of the low-mass white dwarf companions in these systems can be predicted accurately. We notice a clear correlation between spin period and orbital period for these systems, as well as a clear correlation between pulsar magnetic field strength and orbital period. These relations strongly suggest that increased amounts of mass accreted by the neutron stars lead to increased decay of their magnetic fields: we suggest a simple way to understand the observed value of the ‘bottom’ field strengths of a few times 10⁸ G. Furthermore, we find that the LMBP-systems in which the pulsar has a strong magnetic field (> 10¹¹ G) have an about two orders of magnitude larger birth rate (i.e., about 4 × 10^-4 yr^-1 in the Galaxy) than the systems with millisecond pulsars (which have B < 10⁹ G). Using the observational fact that neutron stars receive a velocity kick of ∼450 km/s at birth, we find that some 90% of the potential progenitor systems of the strong-field LMBP must have been disrupted in the Supernovae in which their neutron stars were formed. Hence, the formation rate of the progenitors of the strong-field LMBP is of the same order as the galactic supernova rate (4 × 10^-3 yr^-1). This implies that a large fraction of all Supernovae take place in binaries with a close low-mass (< 2.3 M⊙) companion.     

The strong magnetic fields of the X-ray pulsars Hercules X-1 and 4U 1626–67 and their evolutionary scenarios

We consider the magnetic and spin evolution of the X-ray binary pulsars Her X-1 and 4U 1626–67, assuming that their magnetic fields are of crustal origin. We adopt the standard evolutionary model which implies that the neutron star passes through several phases in a binary system ('isolated pulsar' – propeller – wind accretion – Roche lobe overflow). In the framework of the model under consideration, the strong magnetic fields of relatively old pulsars like Her X-1 and 4U 1626–67 can naturally be understood if, at their birth, they had a sufficiently strong magnetic field, ∼3 × 10¹³ G, comparable to the maximal field observed in radio pulsars.     

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     

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 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.     

Modeling the luminosity function of galactic low-mass X-ray binaries

The evolution of the family of binaries with a low-mass star and a compact neutron star companion (low-mass X-ray binaries (LMXBs) with neutron stars) ismodeled by the method of population synthesis. Continuous Roche-lobe filling by the optical star in LMXBs is assumed to be maintained by the removal of orbital angular momentum from the binary by a magnetic stellar wind from the optical star and the radiation of gravitational waves by the binary. The developed model of LMXB evolution has the following significant distinctions: (1) allowance for the effect of the rotational evolution of a magnetized compact remnant on themass transfer scenario in the binary, (2) amore accurate allowance for the response of the donor star to mass loss at the Roche-lobe filling stage. The results of theoretical calculations are shown to be in good agreement with the observed orbital period-X-ray luminosity diagrams for persistent Galactic LMXBs and their X-ray luminosity function. This suggests that the main elements of binary evolution, on the whole, are correctly reflected in the developed code. It is shown that most of the Galactic bulge LMXBs at luminosities L _x > 10³⁷ erg s^?1 should have a post-main-sequence Roche-lobe-filling secondary component (low-mass giants). Almost all of the models considered predict a deficit of LMXBs at X-ray luminosities near ～10^36.5 erg s^?1 due to the transition of the binary from the regime of angular momentum removal by a magnetic stellar wind to the regime of gravitational waves (analogous to the widely known period gap in cataclysmic variables, accreting white dwarfs). At low luminosities, the shape of the model luminosity function for LMXBs is affected significantly by their transient behavior-the accretion rate onto the compact companion is not always equal to the mass transfer rate due to instabilities in the accretion disk around the compact object. The best agreement with observed binaries is achieved in the models suggesting that heavy neutron stars with masses 1.4–1.9M _⊙ can be born.     

The evolution of the core and surface magnetic fields in isolated neutron stars

We apply the model of flux expulsion from the superfluid and superconductive core of a neutron star, developed by Konenkov & Geppert, both to neutron star models based on different equations of state and to different initial magnetic field structures. Initially, when the core and the surface magnetic field are of the same order of magnitude, the rate of flux expulsion from the core is almost independent of the equation of state, and the evolution of the surface field decouples from the core field evolution with increasing stiffness. When the surface field is initially much stronger than the core field, the magnetic and rotational evolution resembles that of a neutron star with a purely crustal field configuration; the only difference is the occurrence of a residual field. In the case of an initially submerged field, significant differences from the standard evolution only occur during the early period of the life of a neutron star, until the field has been re-diffused to the surface. The reminder of the episode of submergence is a correlation of the residual field strength with the submergence depth of the initial field. We discuss the effect of the re-diffusion of the magnetic field on the difference between the real and the active age of young pulsars and on their braking indices. Finally, we estimate the shear stresses built up by the moving fluxoids at the crust–core interface and show that these stresses may cause crust cracking, preferentially in neutron stars with a soft equation of state.     

X-ray binaries

Summary The various types and classes of X-ray binary are reviewed high-lighting recent results. The high mass X-ray binaries (HMXRBs) can be used to probe the nature of the mass loss from the OB star in these systems. Absorption measurements through one orbital cycle of the supergiant system X1700-37 are well modelled by a radiation driven wind and also require a gas stream trailing behind the X-ray source. In Cen X-3 the gas stream is accreted by the X-ray source via an accretion disk. Changes in the gas stream can cause the disk to thicken and the disk to obscure the X-ray source. How close the supergiant is to corotation seems to be as much a critical factor in these systems as how close it is to filling its Roche lobe. In the Be star X-ray binaries a strong correlation between the neutron stars rotation period and its orbital period has been explained as due to the neutron star being immersed in a dense, slow moving equatorial wind from the Be star. For the X-ray pulsars in the transient Be X-ray binaries a centrifugal barrier to accretion is important in determining the X-ray lightcurve and the spin evolution. The X-ray orbital modulations from the low mass X-ray binaries, LMXRBs, include eclipses by the companion and/or periodic dipping behaviour from structure at the edge of the disk. The corresponding optical modulations show a smooth sinusoidal like component and in some cases a sharp eclipse by the companion. The orbital period of the LMXRB XB1916-05 is 1% longer in the optical compared to that given by the X-ray dip period. The optical period has been interpreted as the orbital period, but this seems inconsistent with the well established view of the origin of the X-ray modulations in LMXRB. A new model is presented that assumes the X-ray dip period is the true orbital period. The 5.2 h eclipsing LMXRB XB2129+47 recently entered a low state and optical observations unexpectedly reveal an F star which is too big to fit into the binary. This is probably the first direct evidence that an X-ray binary is part of a hierarchical triple. Finally the class of X-ray binaries containing black hole candidates is reviewed focusing on the value of using X-ray signatures to identify new candidates.     

Orbital period modulation and magnetic cycles in close binaries

We discuss the observed orbital period modulations in close binaries, and focus on the mechanism proposed by Applegate relating the changes of the stellar internal rotation associated with a magnetic activity cycle with the variation of the gravitational quadrupole moment of the active component; the variation of this quadrupole moment in turn forces the orbital motion of the binary stars to follow the activity level of the active star. We generalize this approach by considering the details of this interaction, and develop some illustrative examples in which the problem can be easily solved in analytical form. Starting from such results, we consider the interplay between rotation and magnetic field generation in the framework of different types of dynamo models, which have been proposed to explain solar and stellar activity. We show how the observed orbital period modulation in active binaries may provide new constraints for discriminating between such models. In particular, we study the case of the prototype active binary RS Canum Venaticorum, and suggest that torsional oscillations — driven by a stellar magnetic dynamo — may account for the observed behaviour of this star. Further possible applications of the relationship between magnetic activity and orbital period modulation, related to the recent discovery of binary systems containing a radio pulsar and a convecting upper main-sequence or a late-type low-mass companion, are discussed.     

Evolution of The Spin Periods of Neutron Stars in Low Mass X-ray Binaries

We present a numerical analysis of the spin evolution of neutron stars in low-mass X-ray binaries, trying to explain the discrepancy in the spin period distribution between observations of millisecond pulsars and theoretical results. In our calculations, we take account of possible effects of radiation pressure and irradiation-induced instability on the structure of the disk, and the evolution of the mass transfer rate, respectively. We report the following results: (1) The radiation pressure in the accretion disk leads to a slight increase of spin periods, and the variation of mass transfer rate caused by the neutron star irradiation can shorten the spin-down phase of evolution. (2) The calculated results of the model combining radiation pressure and irradiation show that the accretion is strongly limited by the radiation pressure in the high mass transfer phase. (3) The accreted mass and fastness parameter can affect the number of systems in the equilibrium state.     

Crust–core coupling and r-mode damping in neutron stars: a toy model

r-modes in neutron stars with crusts are damped by viscous friction at the crust–core boundary. The magnitude of this damping, evaluated by Bildsten & Ushomirsky (BU) under the assumption of a perfectly rigid crust, sets the maximum spin frequency for neutron stars spun up by accretion in low-mass X-ray binaries (LMXBs). In this paper we explore the mechanical coupling between the core r-modes and the elastic crust, using a toy model of a constant-density neutron star having a crust with a constant shear modulus. We find that, at spin frequencies in excess of ≈50 Hz, the r-modes strongly penetrate the crust. This reduces the relative motion (slippage) between the crust and the core compared with the rigid-crust limit. We therefore revise down, by as much as a factor of 10²–10³ , the damping rate computed by BU, significantly reducing the maximal possible spin frequency of neutron stars with solid crusts. The dependence of the crust–core slippage on the spin frequency is complicated, and is very sensitive to the physical thickness of the crust. If the crust is sufficiently thick, the curve of the critical spin frequency for the onset of the r-mode instability becomes multivalued for some temperatures; this is related to avoided crossings between the r-mode and higher-order torsional modes in the crust. The critical frequencies are comparable to the observed spins of neutron stars in LMXBs and millisecond pulsars.     

Progenitors with enhanced rotation and the origin of magnetars

Among the dozen known magnetar candidates, there are no binary objects. Given that the fraction of binary neutron stars is estimated to be about 3–10 per cent, it is reasonable to address the question of solitarity of magnetars, to estimate theoretically the fraction of binary objects among them, and to identify the most probable companions. We present population synthesis calculations of massive binary systems. In this study, we adopt the hypothesis that magnetic field of a magnetar is generated at the protoneutron star stage due to a dynamo mechanism, so rapid rotation of the core of a progenitor star is essential. Our goal is to estimate the number of neutron stars originated from progenitors with enhanced rotation. In our calculations, the fraction of neutron stars originating from such progenitors is about 8–9 per cent. This should be considered as an upper limit to the fraction of magnetars, as some of the progenitors can lose momentum. Most of these objects are isolated due to coalescences of components prior to neutron star formation, or due to system disruption after the second supernova explosion. The fraction of such neutron stars in surviving binaries is about 1 per cent or lower. Their most numerous companions are black holes.     

Neutron star spin–kick velocity correlation effect on binary neutron star coalescence rates and spin–orbit misalignment of the components

We study the effect of the neutron star spin–kick velocity alignment observed in young radio pulsars on the coalescence rate of binary neutron stars. Two scenarios are considered for neutron star formation: when the kick is always present, and when it is small or absent if a neutron star is formed in a binary system as a result of electron-capture degenerate core collapse. The effect is shown to be especially strong for large kick amplitudes and tight alignments, reducing the expected galactic rate of binary neutron star coalescence compared to calculations with randomly directed kicks. The spin–kick correlation also leads to a much narrower neutron star spin–orbit misalignment.     

Fallback Disk-Involved Spin-Down of Young Radio Pulsars

Disks originating from supernova fallback have been suggested to surround young neutron stars. Interaction between the disk and the magnetic field of the neutron star may considerably influence the evolution of the star through the so called propeller effect. There are many controversies about the efficiency of the propeller mechanism proposed in the literature. We investigate the fallback disk-involved spin-down of young pulsars. By comparing the simulated and measured results of pulsar evolution, we present some possible constraints on the propeller torques exerted by the disks on neutron stars.     

The spin evolution of nascent neutron stars

The loss of angular momentum owing to unstable r-modes in hot young neutron stars has been proposed as a mechanism for achieving the spin rates inferred for young pulsars. One factor that could have a significant effect on the action of the r-mode instability is fallback of supernova remnant material. The associated accretion torque could potentially counteract any gravitational-wave-induced spin-down, and accretion heating could affect the viscous damping rates and hence the instability. We discuss the effects of various external agents on the r-mode instability scenario within a simple model of supernova fallback on to a hot young magnetized neutron star. We find that the outcome depends strongly on the strength of the magnetic field of the star. Our model is capable of generating spin rates for young neutron stars that accord well with initial spin rates inferred from pulsar observations. The combined action of r-mode instability and fallback appears to cause the spin rates of neutron stars born with very different spin rates to converge, on a time-scale of approximately 1 year. The results suggest that stars with magnetic fields ≤10¹³ G could emit a detectable gravitational wave signal for perhaps several years after the supernova event. Stars with higher fields (magnetars) are unlikely to emit a detectable gravitational wave signal via the r-mode instability. The model also suggests that the r-mode instability could be extremely effective in preventing young neutron stars from going dynamically unstable to the bar-mode.