The origin of subdwarf B stars - II

Abstract

We have carried out a detailed binary population synthesis (BPS) study of the formation of subdwarf B (sdB) stars and related objects (sdO, sdOB stars) using the latest version of the BPS code developed by Han and co-workers. We systematically investigate the importance of the five main evolutionary channels in which the sdB stars form after one or two common-envelope (CE) phases, one or two phases of stable Roche lobe overflow (RLOF) or as the result of the merger of two helium white dwarfs (WDs). Our best BPS model can satisfactorily explain the main observational characteristics of sdB stars, in particular their distributions in the orbital period-minimum companion mass (log P-Mcomp) diagram and in the effective temperature-surface gravity (Teff-log g) diagram, their distributions of orbital period, log (gθ4) (θ = 5040 K /Teff) and mass function, their binary fraction and the fraction of sdB binaries with WD companions, their birth rates and their space density. We obtain a Galactic formation rate for sdB stars of 0.014-0.063 yr-1 with a best estimate of ∼0.05 yr-1 and a total number in the Galaxy of 2.4-9.5 × 106 with a best estimate of ∼6 × 10 6; half of these may be missing in observational surveys owing to selection effects. The intrinsic binary fraction is 76-89 per cent, although the observed frequency may be substantially lower owing to the selection effects. The first CE ejection channel, the first stable RLOF channel and the merger channel are intrinsically the most important channels, although observational selection effects tend to increase the relative importance of the second CE ejection and merger channels. We also predict a distribution of masses for sdB stars that is wider than is commonly assumed and that some sdB stars have companions of spectral type as early as B. The percentage of A-type stars with sdB companions can in principle be used to constrain some of the important parameters in the binary evolution model. We conclude that (i) the first RLOF phase needs to be more stable than is commonly assumed, either because the critical mass ratio qcrit for dynamical mass transfer is higher or because of tiddly enhanced stellar wind mass loss; (ii) mass transfer in the first stable RLOF phase is non-conservative, and the mass lost from the system takes away a specific angular momentum similar to that of the system; and (iii) common-envelope ejection is very efficient.