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.     

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.     

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.     

High-mass X-ray binaries and the spiral structure of the host galaxy

We investigate the manifestation of the spiral structure in the distribution of high-mass X-ray binaries (HMXBs) over the host galaxy. We construct the simplest kinematic model. It shows that the HMXBs should be displaced relative to the spiral structure observed in such traditional star formation rate indicators as the Hα and far-infrared emissions because of their finite lifetimes. Using Chandra observations of M51, we have studied the distribution of X-ray sources relative to the spiral arms of this galaxy observed in Hα. Based on K-band data and background source number counts, we have separated the contributions from high-mass and low-mass X-ray binaries and active galactic nuclei. In agreement with model predictions, the distribution of HMXBs is wider than that of bright H II regions concentrated in the region of ongoing star formation. However, the statistical significance of this result is low, as is the significance of the concentration of the total population of X-ray sources to the spiral arms. We also predict the distribution of HMXBs in our Galaxy in Galactic longitude. The distribution depends on the mean HMXB age and can differ significantly from the distributions of such young objects as ultracompact H II regions.     

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.     

The secondary stars in cataclysmic variables and low-mass X-ray binaries

We critically re-examine the available data on the spectral types, masses and radii of the secondary stars in cataclysmic variables (CVs) and low-mass X-ray binaries (LMXBs), using the new catalogue of Ritter &38; Kolb as a starting point. We find there are 55 reliable spectral type determinations and only 14 reliable mass determinations of CV secondary stars (10 and 5, respectively, in the case of LMXBs). We derive new spectral type–period, mass–radius, mass–period and radius–period relations, and compare them with theoretical predictions. We find that CV secondary stars with orbital periods shorter than 7–8 h are, as a group, indistinguishable from main-sequence stars in detached binaries. We find that it is not valid, however, to estimate the mass from the spectral type of the secondary star in CVs or LMXBs. We find that LMXB secondary stars show some evidence for evolution, with secondary stars which are slightly too large for their mass. We show how the masses and radii of the secondary stars in CVs can be used to test the validity of the disrupted magnetic braking model of CV evolution, but we find that the currently available data are not sufficiently accurate or numerous to allow such an analysis. As well as considering secondary star masses, we also discuss the masses of the white dwarfs in CVs, and find mean values of ⁻ M  = 0.69 ± 0.13 M_⊙ below the period gap, and ⁻ M  = 0.80 ± 0.22 M_⊙ above the period gap.     

On the metallicity dependence of high-mass X-ray binaries

It is commonly assumed that high-mass X-ray binary (HMXB) populations are little affected by metallicity. However, the massive stars making up their progenitor systems depend on metallicity in a number of ways, not least through their winds. We present simulations, well-matched to the observed sample of Galactic HMXBs, which demonstrate that both the number and the mean period of HMXB progenitors can vary with metallicity, with the number increasing by about a factor of 3 between solar and Small Magellanic Cloud (SMC) metallicity. However, the SMC population itself cannot be explained simply by metallicity effects; it requires both that the HMXBs observed therein primarily sample the older end of the HMXB population and that the star formation rate at the time of their formation was very large.     

中子星X射线双星中kHz QPO现象的理论解释

罗西X射线时变探测器(RXTE)在中子星小质量X射线双星中发现了千赫兹准周期振荡现象(kHzQPO)。kHzQPO的频率一般在几百到上千赫兹,其动力学时标与吸积盘最内部区域物质的运动时标一致,因此普遍认为kHz QPO产生于中子星表面附近区域,携带了来自中心中子星及周围强引力场信息,如质量、自转周期、角动量、半径、磁场等。kHz QpO现象的理解为研究强引力场和致密物质状态开启了一扇新的窗口。着重介绍基于kHz QPO的基本现象和相应的理论模型。     

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.     

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.     

The magnetohydrodynamics model of twin kilohertz quasi-periodic oscillations in low-mass X-ray binaries

We suggest an explanation for the twin kilohertz quasi-periodic oscillations (kHz QPOs) in low-mass X-ray binaries (LMXBs) based on magnetohydrodynamics (MHD) oscillation modes in neutron star magnetospheres. Including the effect of the neutron star spin, we derive several MHD wave modes by solving the dispersion equations, and propose that the coupling of the two resonant MHD modes may lead to the twin kHz QPOs. This model naturally relates the upper, lower kHz QPO frequencies with the spin frequencies of the neutron stars, and can well account for the measured data of six LMXBs.     

The population of faint transients in the Galactic Centre

BeppoSAX has detected a population of faint transient X-ray sources in the Galactic Centre. I show that a simple irradiated disc picture gives a consistent fit to the properties of this population, and that it probably consists of low-mass X-ray binaries (LMXBs) that have evolved beyond their minimum orbital periods of ∼80 min. Since all post-minimum systems are transient, and neutron star LMXBs are more common than black hole LMXBs in the Galaxy, the majority of these systems should contain neutron stars, as observed. This picture predicts that the Galactic Centre transients should have orbital periods in the range ∼80–120 min, and that most of them should repeat in the next few years. In this case, the total number of post-minimum transients in the Galaxy would be considerably smaller than the usual estimates of its total LMXB population. I discuss possible reasons for this.     

Late stages of the evolution of close binaries

Close binaries can evolve through various ways of interaction into compact objects (white dwarfs, neutron stars, black holes). Massive binary systems (mass of the primaryM ₁ larger than 14 to 15M ₀) are expected to leave, after the first stage of mass transfer a compact component orbiting a massive star. These systems evolve during subsequent stages into massive X-ray binaries. Systems with initial large periode evolve into Be X-ray binaries.Low mass X-ray sources are probably descendants of lower mass stars, and various channels for their production are indicated. The evolution of massive close binaries is examined in detail and different X-ray stages are discussed. It is argued that a first X-ray stage is followed by a reverse extensive mass transfer, leading to systems like SS 433, Cir X1. During further evolution these systems would become Wolf-Rayet runaways. Due to spiral in these system would then further evolve into ultra short X-ray binaries like Cyg X-3.Finally the explosion of the secondary will in most cases disrupt the system. In an exceptional case the system remains bound, leading to binary pulsars like PSR 1913+16. In such systems the orbit will shrink due to gravitational radiation and finally the two neutron stars will coalesce. It is argued that the millisecond pulsar PSR 1937+214 could be formed in this way.A complete scheme starting from two massive ZAMS stars, ending with a millisecond pulsar is presented.Paper presented at the Lembang-Bamberg IAU Colloquium No. 80 on Double Stars: Physical Properties and Generic Relations, held at Bandung, Indonesia 3–7 June, 1983.     

Neutron star retention and millisecond pulsar production in globular clusters

We investigate the conditions by which neutron star retention in globular clusters is favoured. We find that neutron stars formed in massive binaries are far more likely to be retained. Such binaries are likely to then evolve into contact before encountering other stars, possibly producing a single neutron star after a common envelope phase. A large fraction of the single neutron stars in globular clusters are then likely to exchange into binaries containing moderate-mass main-sequence stars, replacing the lower-mass components of the original systems. These binaries will become intermediate-mass X-ray binaries (IMXBs), once the moderate-mass star evolves off the main sequence, as mass is transferred on to the neutron star, possibly spinning it up in the process. Such systems may be responsible for the population of millisecond pulsars (MSPs) that has been observed in globular clusters. Additionally, the period of mass-transfer (and thus X-ray visibility) in the vast majority of such systems will have occurred 5–10 Gyr ago, thus explaining the observed relative paucity of X-ray binaries today, given the MSP population.     

Are the magnetic fields of millisecond pulsars ∼108 G?

It is generally assumed that the magnetic fields of millisecond pulsars (MSPs) are ～10⁸ G. We argue that this may not be true and the fields may be appreciably greater. We present six evidences for this: (1) The ～10⁸G field estimate is based on magnetic dipole emission losses which is shown to be questionable; (2) The MSPs in low mass X-ray binaries (LMXBs) are claimed to have <10¹¹ G on the basis of a Rayleygh-Taylor instability accretion argument. We show that the accretion argument is questionable and the upper limit 10¹¹ G may be much higher; (3) Low magnetic field neutron stars have difficulty being produced in LMXBs; (4) MSPs may still be accreting indicating a much higher magnetic field; (5) The data that predict ～10⁸ G for MSPs also predict ages on the order of, and greater than, ten billion years, which is much greater than normal pulsars. If the predicted ages are wrong, most likely the predicted ～10⁸ G fields of MSPs are wrong; (6) When magnetic fields are measured directly with cyclotron lines in X-ray binaries, fields ?10⁸ G are indicated. Other scenarios should be investigated. One such scenario is the following. Over 85% of MSPs are confirmed members of a binary. It is possible that all MSPs are in large separation binaries having magnetic fields >10⁸ G with their magnetic dipole emission being balanced by low level accretion from their companions.     

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.     

High energy gamma-ray emission from binaries

The current Cherenkov telescopes together with GLAST are opening up a new window into the physics at work close to black holes and rapidly rotating neutron stars with great breakthrough potential. Very high energy gamma-ray emission up to 10 TeV is now established in several binaries. The radiative output of gamma-ray binaries is in fact dominated by emission above 1–10 MeV. Most are likely powered by the rotational spindown of a young neutron star that generates a highly relativistic wind. The interaction of this pulsar wind with the companion’s stellar wind is responsible for the high energy gamma-ray emission. There are hints that microquasars, accretion-powered binaries emitting relativistic jets, also emit gamma-ray flares that may be linked to the accretion–ejection process. Studying high energy gamma-ray emission from binaries offers good prospects for the study of pulsar winds physics and may bring new insights into the link between accretion and ejection close to black holes.     

Radius of the neutron star magnetosphere during disk accretion

The dependence of the spin frequency derivative \(\dot \nu \) of accreting neutron stars with a strongmagnetic field (X-ray pulsars) on the mass accretion rate (bolometric luminosity, L_bol) has been investigated for eight transient pulsars in binary systems with Be stars. Using data from the Fermi/GBM and Swift/BAT telescopes, we have shown that for seven of the eight systems the dependence \(\dot \nu \) (L_bol) can be fitted by the model of angular momentum transfer through an accretion disk, which predicts the relation \(\dot \nu \) ～ L^6/7_bol. Hysteresis in the dependence \(\dot \nu \) (L_bol) has been confirmed in the system V 0332+53 and has been detected for the first time in the systems KS 1947+300, GRO J1008-57, and 1A 0535+26. Estimates for the radius of the neutron star magnetosphere in all of the investigated systems have been obtained. We show that this quantity varies from pulsar to pulsar and depends strongly on the analytical model and the estimates for the neutron star and binary system parameters.