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.     

Massive X-ray binaries: Their physics and evolution

We review various aspects of the evolutionary history of massive X-ray binaries. It is expected that moderately massive close binaries evolve to Be X-ray binaries, while very massive systems evolve to standard X-ray binaries.The compact objects are formed through supernova explosions. The fairly low galactic latitudes of those systems indicate that the explosion should, in general, not have accelerated the system to a velocity larger than 50kms^–1. This implies that the mass of the exploding stars is in general less than 5 to 6M .After the explosion, tidal forces will circularize the orbit of short period systems. Even if the tidal evolution has been completed, the expansion of the optical star during the course of its evolution will continously disturb the stability of the orbit. Short period systems with large mass ratio may eventually become tidally unstable. Cen X-3 may be an example of such a system. The predicted rate of the orbital period decrease of Cen X-3 is in agreement with the observed rate.A way to represent the rotational and magnetic evolution of neutron stars in close binary systems is presented. The observed distribution of the pulsation periods of X-ray pulsars with Be companions is consistent with initial magnetic fields of 10¹²–10¹³ G of the neutron stars. We suggest that the fast X-ray pulsars 4U 0115+63 and A 0538-66 are young neutron stars, while Cen X-3 and SMC X-1 are recycled pulsars.The evolutionary relationship between massive X-ray binaries, binary pulsars, and millisecond pulsars is also discussed.Invited paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.     

Massive stars

We describe the present state of massive star research seen from the viewpoint of stellar evolution, with special emphasis on close binaries. Statistics of massive close binaries are reasonably complete for the Solar neighbourhood. We defend the thesis that within our knowledge, many scientific results where the effects of binaries are not included, have an academic value, but may be far from reality. In chapter I, we summarize general observations of massive stars where we focus on the HR diagram, stellar wind mass loss rates, the stellar surface chemistry, rotation, circumstellar environments, supernovae. Close binaries can not be studied separately from single stars and vice versa. First, the evolution of single stars is discussed (chapter I). We refer to new calculations with updated stellar wind mass loss rate formalisms and conclusions are proposed resulting from a comparison with representative observations. Massive binaries are considered in chapter II. Basic processes are briefly described, i.e. the Roche lobe overflow and mass transfer, the common envelope process, the spiral-in process in binaries with extreme mass ratio, the effects of mass accretion and the merging process, the implications of the (asymmetric) supernova explosion of one of the components on the orbital parameters of the binary. Evolutionary computations of interacting close binaries are discussed and general conclusions are drawn. The enormous amount of observational data of massive binaries is summarized. We separately consider the non-evolved and evolved systems. The latter class includes the semi-detached and contact binaries, the WR binaries, the X-ray binaries, the runaways, the single and binary pulsars. A general comparison between theoretical evolution and observations is combined with a discussion of specially interesting binaries: the evolved binaries HD 163181, HD 12323, HD 14633, HD 193516, HD 25638, HD 209481, Per and Sgr; the WR+OB binary V444 Cyg; the high mass X-ray binaries Vela X-1, Wray 977, Cyg X-1; the low mass X-ray binaries Her X-1 and those with a black hole candidate; the runaway Pup, the WR+compact companion candidates Cyg X-3, HD 50896 and HD 197406. We finally propose an overall evolutionary model of massive close binaries as a function of primary mass, mass ratio and orbital period. Chapter III deals with massive star population synthesis with a realistic population of binaries. We discuss the massive close binary frequency, mass ratio and period distribution, the observations that allow to constrain possible asymmetries during the supernova explosion of a massive star. We focuss on the comparison between observed star numbers (as a function of metallicity) and theoretically predicted numbers of stellar populations in regions of continuous star formation and in starburst regions. Special attention is given to the O-type star/WR star/red supergiant star population, the pulsar and binary pulsar population, the supernova rates. Received 17 July 1998     

Frequency of Be and Ae stars in spectroscopic binaries

The frequency of Be and Ae stars in spectroscopic binaries and Algol binaries is examined based on the available catalogue data, and compared with that of stars in theBright Star Catalogue. It is shown that in binary systems Be and Ae stars reveal different statistical behaviours. The frequency of Be stars in spectroscopic binaries shows a notable peak in its orbital-period distribution in the range of 100–300 days, suggesting a statistical group separated from Algol systems. In contrast, most of Ae stars in spectroscopic binaries belong to the Algol systems.     

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.     

Nearby young single black holes

We consider nearby young black holes formed after supernova explosions in close binaries whose secondary components are currently observed as the so-called runaway stars. Using data on runaway stars and making reasonable assumptions about the mechanisms of supernova explosion and binary breakup, we estimate the present positions of nearbyyoung black holes. For two objects, we obtained relatively small error regions (～50–100 deg²). The possibility of detecting these nearby young black holes is discussed.     

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.     

Influence of Neutron Star Parameters on Evolutions of Different Types of Pulsar

Influences of the mass, moment of inertia, rotation, absence of stability in the atmosphere and some other parameters of neutron stars on the evolution of pulsars are examined. It is shown that the locations and evolutions of soft gamma repeaters, anomalous X-ray pulsars and other types of pulsar on the period versus period derivative diagram can be explained adopting values of B < 10¹⁴ G for these objects if they have smaller mass (e.g. about 0.5 Solar mass) compared to the conventionally adopted values of mass. This approach gives the possibility to explain many properties of different types of pulsar.     

Pulsars: Space kinematic characteristics and their connection with supernova remnants

The galactic distribution of pulsars and shell remnants of supernovae (SN) as investigated on the basis of newly-estimated parameters. Special attention was paid to taking into account all possible selection effects and an attempt was made to reveal a statistically-pure ensemble of objects. On the basis of this ensemble we studied the radial andz-distribution of pulsars and supernova remnants (SNR).It is shown that the radial distribution of both objects is considered to have an annular structure with the maximum surface density at a distance of 4.5–5.3 kpc (if the distance of the Sun from the galactic centre is assumed to be equal to 8.5 kpc). The scle-height of the progenitors of SNRs is not more than 110 pc and only 15% of the SNRs, whose progenitors may also be massive runaway stars, are situated at 300 pc. The mean application of the pulsars is not more than 300 pc which agrees with the hypothesis about the genetic connection with type-II SN outbursts at the kinematic aget _k5×10⁶ yr and thez-component of spatial velocity beingV _z=120 km s^–1.The possible precursors of type-I SNRs by the shape of their radial distribution in late spirals and the various model calculations given in the literature, are also discussed.     

On the origin of high-velocity runaway stars

We explore the hypothesis that some high-velocity runaway stars attain their peculiar velocities in the course of exchange encounters between hard massive binaries and a very massive star (either an ordinary  50–100 M_⊙  star or a more massive one, formed through runaway mergers of ordinary stars in the core of a young massive star cluster). In this process, one of the binary components becomes gravitationally bound to the very massive star, while the second one is ejected, sometimes with a high speed. We performed three-body scattering experiments and found that early B-type stars (the progenitors of the majority of neutron stars) can be ejected with velocities of  ≳200–400 km s⁻¹  (typical of pulsars), while  3–4 M_⊙  stars can attain velocities of  ≳300–400 km s⁻¹  (typical of the bound population of halo late B-type stars). We also found that the ejected stars can occasionally attain velocities exceeding the Milky Ways's escape velocity.     

On evolution of contact eclipsing binaries

It is shown that during contact eclipsing binaries evolution under the influence of stellar wind, magnetic stellar wind and with matter transfer by gas flow, in binary stellar systems there may take place a process of star merger (low mass stars) within 10⁵–10⁷ yr and a fast increase of distance between stars of massive binaries. W UMa-type stars are a finite evolutionary stage of very close and low mass binary pairs. As for contact systems of early spectral types (CE-systems), they are more varied in evolution.     

A near-infrared survey of IRAS sources in the South Galactic Cap

The paper describes the JHK colours of late-type stars which were investigated as part of a survey of South Galactic Cap (b < -30°) IRAS sources selected on the basis of their 12/25µm flux ratios as high mass-loss candidates. Near-infrared two-colour diagrams provide an effective technique for distinguishing between various groups of late-type stars. Such diagrams are also useful in indicating which stars are likely to be peculiar and worthy of more detailed study. The late-type stars isolated by this survey comprise: 61 Mira variables (3 of which are carbon stars with very thick shells), 3 young stellar objects, 4 interacting binaries, 2 semi-regular carbon variables and 154 oxygen-rich giants.     

The HYPER-MUCHFUSS project—target selection and analysis

The HYPER-MUCHFUSS project targets a population of high velocity subluminous B stars to discover either close binaries with massive unseen companions or hyper-velocity stars. Our starting point is the enormous database of SDSS. We preselected sdO/B candidates by colour and classified them by visual inspection of their spectra. We measured the radial velocity from the coadded SDSS spectra, which serves as first epoch measurement. Stars with high Galactic rest-frame velocities were selected and second epoch observations were obtained starting in 2007 at several sites. For the brighter targets we also included the SDSS individual spectra as additional information. In the course of our survey we observed 88 out of 265 stars from our target list. We discovered 39 HVS candidates as well as 49 close binaries. In addition we analysed all single spectra of sdBs from SDSS and found 120 close binaries. For the targets with constant RVs we performed a proper motion analysis with the highest possible accuracy from the available digitised photographic plates. Together with the analysed spectra and the calculation of the spectroscopic distance, we calculated complete trajectories and deduced the origins of these stars. Targets with high RV variability on short timescales were selected for follow-up. Numerical simulations based on the period and companion mass distribution of the known sdB binary sample were carried out to optimise the target selection and single out candidate binaries with massive companions. The follow-up campaign using WHT/ISIS and CAHA-3.5m/TWIN started in 2009.     

Evolution of massive close binaries and formation of neutron stars and black holes

Main results of computations of evolution for massive close binaries (10M +9.4M , 16M +15M , 32M +30M , 64M +60M ) up to oxygen exhaustion in the core are described. Mass exchange starting in core hydrogen, shell hydrogen and core helium burning stages was studied. Computations were performed assuming both the Ledoux and Schwarzschild stability criteria for semiconvection. The influence of UFI-neutrino emission on evolution of close binaries was investigated. The results obtained allow to outline the following evolutionary chain: two detached Main-Sequence stars — mass exchange — Wolf-Rayet star or blue supergiant plus main sequence star — explosion of the initially more massive star appearing as a supernova event — collapsed or neutron star plus Main-Sequence star, that may be observed as a runaway star — mass exchange leading to X-rays emission — collapsed or neutron star plus WR-star or blue supergiant — second explosion of supernova that preferentially disrupts the system and gives birth to two single high spatial velocity pulsars.Numerical estimates concerning the number and properties of WR-stars, pulsars and X-ray sources are presented. The results are in favour of the existence of UFI-neutrino and of the Ledoux criterion for describing semiconvection. Properties of several well-known X-ray sources and the binary pulsar are discussed on base of evolutionary chain of close binaries.     

The HYPER-MUCHFUSS project—the constant high-velocity population

The HYPER-MUCHFUSS (HYPER-velocity stars or Massive Unseen Companions of Hot Faint Underluminous Stars from SDSS) project targets a population of high-velocity subluminous B stars to discover either close binaries with massive unseen companions or hyper-velocity stars. We re-observed high-velocity subdwarf selected candidates from the SDSS spectroscopic Data Release 6. Starting in 2007 we used several instruments and have now reached a completion level of 33% (from 265 targets), whereas we found at least 16 close binaries. Here we present results for two of our 39 hyper-velocity star candidates. From the available Digitized Sky Surveys photographic plates we measured a significant proper motion for 14 stars. Combining this information with a detailed spectroscopic analysis allows for the first time a complete determination of the 3D-trajectories for a high-velocity sample. We present our preliminary results for the two subdwarfs J1644+4523 and J1211+1437. Assuming the Standard Allen and Santillan (Rev. Mex. Astron. Astrofis. 22:255, 1991) potential the first one is bound and originates in the central region of the Galaxy. The subdwarf B star J1211+1437 is possibly unbound and seems to originate in the Galactic rim. We also performed numerical kinematical experiments with increased dark matter halo mass. and found that the origin of J1644+4523 in the central region is not changed but the time-of-flight is drastically shortened. J1211+1437 would be bound and probably belongs to population II.     

Massive early-type emission-line stars: an attempt to distinguish a new group of stellar objects

A new sample of possibly massive early-type emission-line stars (METELS) based on the previous lists of peculiar Be stars is presented. It consists of 36 objects divided amongst supergiants, possible binaries, and candidates to the list. The central stars are probably more massive than 10M . Two new relations allowing idientification of possible binaries among the objects are proposed.     

High-luminosity single carbon stars in stellar and galactic evolution

Summary In the solar neighborhood, approximately half of all intermediate mass main sequence stars with initially between 1 and about 5 Mbecome carbon stars with luminosities near 10⁴ L for typically less than 10⁶ years. These high luminosity carbon stars lose mass at rates nearly always in excess of 10^–7 M yr^–1 and sometimes in excess of 10^–5 M yr^–1. Locally, close to half of the mass returned into the interstellar medium by intermediate mass stars before they become white dwarfs is during the carbon star phase. A much greater fraction of lower metallicity stars become carbon-rich before they evolve into planetary nebulae than do higher metallicity stars; therefore, carbon stars are much more importan t in the outer than in the inner Galaxy.     

A direct imaging search for close stellar and sub‐stellar companions to young nearby stars

A total of 28 young nearby stars (ages ≤60 Myr) have been observed in the K_s‐band with the adaptive optics imager Naos‐Conica of the Very Large Telescope at the Paranal Observatory in Chile. Among the targets are ten visual binaries and one triple system at distances between 10 and 130 pc, all previously known. During a first observing epoch a total of 20 faint stellar or sub‐stellar companion‐candidates were detected around seven of the targets. These fields, as well as most of the stellar binaries, were re‐observed with the same instrument during a second epoch, about one year later. We present the astrometric observations of all binaries. Their analysis revealed that all stellar binaries are co‐moving. In two cases (HD 119022 AB and FG Aqr B/C) indications for significant orbital motions were found. However, all sub‐stellar companion candidates turned out to be non‐moving background objects except PZ Tel which is part of this project but whose results were published elsewhere. Detection limits were determined for all targets, and limiting masses were derived adopting three different age values; they turn out to be less than 10 Jupiter masses in most cases, well below the brown dwarf mass range. The fraction of stellar multiplicity and of the sub‐stellar companion occurrence in the star forming regions in Chamaeleon are compared to the statistics of our search, and possible reasons for the observed differences are discussed. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Relative astrometry and near‐infrared differential photometry of nearby southern orbital binaries with adaptive optics

We present relative astrometry and differential photometry measurements for a sample of nearby southern orbital binaries making use of the technique of Adaptive Optics. The observations were made in December 2000, with the ADONIS camera mounted at the 3.6‐m ESO telescope from La Silla Observatory, equipped with the broad‐band near‐infrared filters (J ‐, H ‐, K ‐passbands). Our sample contains stars which do not fit very well the empirical mean mass‐luminosity relation (according to our previous study), but for which accurate parallaxes (determined by the Hipparcos satellite) and high‐quality orbits were available thanks to many previous efforts. We derived accurate positions and J, H, K magnitudes of the individual components of those binaries. The individual stellar components have near‐infrared colour indices well grouped in those plots and are comparable to standard single stars. The data reduction procedure used for deriving those results is described in detail. It is based on a least‐squares fit of Moffat‐Lorentz profiles in direct imaging for well‐resolved systems and on Fourier analysis for very close pairs. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)