Mass transfer in eccentric binary stars

The concept of Roche lobe overflow is fundamental to the theory of interacting binaries. Based on potential theory, it is dependent on all the relevant material corotating in a single frame of reference. Therefore if the mass losing star is asynchronous with the orbital motion or the orbit is eccentric, the simple theory no longer applies and no exact analytical treatment has been found. We use an analytic approximation whose predictions are largely justified by smoothed particle hydrodynamic simulations (SPH). We present SPH simulations of binary systems with the same semi-major axis   a = 5.55 R_⊙  , masses   M ₁= 1 M_⊙, M ₂= 2 M_⊙  and radius   R ₁= 0.89 R_⊙  for the primary star but with different eccentricities   e = 0.4, 0.5, 0.6  and 0.7. In each case the secondary star is treated as a point mass. When   e = 0.4  no mass is lost from the primary while at   e = 0.7  catastrophic mass transfer, partly through the L₂ point, takes place near periastron. This would probably lead to common-envelope evolution if star 1 were a giant or to coalescence for a main-sequence star. In between, at   e ≥ 0.5  , some mass is lost through the L₁ point from the primary close to periastron. However, rather than being all accreted by the secondary, some of the stream appears to leave the system. Our results indicate that the radius of the Roche lobe is similar to circular binaries when calculated for the separation and angular velocity at periastron. Part of the mass loss occurs through the L₂ point.     

KPD 0422 + 5421: a new short-period subdwarf B/white dwarf binary

The subdwarf B (sdB) star KPD 0422 + 5421 was discovered to be a single-lined spectroscopic binary with a period of P  = 0.090 1795 ± (3 × 10⁻⁷) d (2 h 10 min). The U B light curves display an ellipsoidal modulation with amplitudes of ≈ 0.02 mag. The sdB star contributes nearly all of the observed flux. This and the absence of any reflection effect suggest that the unseen companion star is small (i.e. R _comp ≈ 0.01 R_⊙) and therefore degenerate. We modelled the U B light curves and derived i  = 78.05° ± 0.50° and a mass ratio of q  =  M _comp/ M _sdB = 0.87 ± 0.15. The sdB star fills 69 per cent of its Roche lobe. These quantities may be combined with the mass function of the companion [ f ( M ) = 0.126 ± 0.028 M_⊙] to derive M _sdB = 0.72 ± 0.26 M_⊙ and M _comp = 0.62 ±  0.18 M_⊙. We used model spectra to derive the effective temperature, surface gravity and helium abundance of the sdB star. We found T _eff = 25 000 ± 1500 K, log g  = 5.4 ± 0.1 and [He/H] = −1.0. With a period of 2 h 10 min, KPD 0422 + 5421 has one of the shortest known orbital periods of a detached binary. This system is also one of only a few known binaries that contain a subdwarf B star and a white dwarf. Thus KPD 0422 + 5421 represents a relatively unobserved, and short-lived, stage of binary star evolution.     

Physical parameters of the O6.5V+B1V eclipsing binary system LS 1135

The 'All Sky Automated Survey' (ASAS) photometric observations of LS 1135, an O-type single-lined binary (SB1) system with an orbital period of 2.7 d, show that the system is also eclipsing performing a numerical model of this binary based on the Wilson–Devinney method. We obtained an orbital inclination     . With this value of the inclination, we deduced masses   M ₁∼ 30 ± 1 M_⊙  and   M ₂∼ 9 ± 1 M_⊙  , and radii   R ₁∼ 12 ± 1 R_⊙  and   R ₂∼ 5 ± 1 R_⊙  for primary and secondary components, respectively. Both the components are well inside their respective Roche lobes. Fixing the T _eff of the primary to the value corresponding to its spectral type (O6.5V), the T _eff obtained for the secondary component corresponds approximately to a spectral type of B1V. The mass ratio   M ₂/ M ₁∼ 0.3  is among the lowest known values for spectroscopic binaries with O-type components.     

The radial velocity of the companion star in the low-mass X-ray binary 2S 0921–630: limits on the mass of the compact object

In this paper we report on optical spectroscopic observations of the low-mass X-ray binary 2S 0921–630 obtained with the Very Large Telescope. We found sinusoidal radial velocity variations of the companion star with a semi-amplitude of  99.1 ± 3.1 km s⁻¹  modulated on a period of 9.006 ± 0.007 d, consistent with the orbital period found previously for this source, and a systemic velocity of  44.4 ± 2.4 km s⁻¹  . Owing to X-ray irradiation, the centre of light measured by the absorption lines from the companion star is probably shifted with respect to the centre of mass. We try to correct for this using the so-called K -correction. Conservatively applying the maximum correction possible and using the previously measured rotational velocity of the companion star, we find a lower limit to the mass of the compact object in 2S 0921–630 of   M_X sin³ i > 1.90 ± 0.25 M_⊙  (1σ errors). The inclination in this system is well constrained since partial eclipses have been observed in X-ray and optical bands. For inclinations in the range  60° < i < 90°  we find  1.90 ± 0.25 < M_X < 2.9 ± 0.4 M_⊙  . However, using this maximum K -correction we find that the ratio between the mass of the companion star and that of the compact object, q , is 1.32 ± 0.37, implying super-Eddington mass-transfer rates; however, evidence for that has not been found in 2S 0921–630. We conclude that the compact object in 2S 0921–630 is either a (massive) neutron star or a low-mass black hole.     

Roche tomography of cataclysmic variables – III. Star-spots on AE Aqr

We present a Roche tomography reconstruction of the secondary star in the cataclysmic variable AE Aqr. The tomogram reveals several surface inhomogeneities that are due to the presence of large, cool star-spots. In addition to a number of lower latitude spots, the maps also show the presence of a large, high-latitude spot similar to that seen in Doppler images of rapidly rotating isolated stars, and a relative paucity of spots at a latitude of 40°. In total, we estimate that some 18 per cent of the Northern hemisphere of AE Aqr is spotted.
We have also applied the entropy landscape technique to determine accurate parameters for the binary system. We obtain optimal masses   M ₁= 0.74 M_⊙, M ₂= 0.50 M_⊙  , a systemic velocity  γ=−63 km s⁻¹  and an orbital inclination   i = 66°  .
Given that this is the first study to successfully image star-spots on the secondary star in a cataclysmic variable, we discuss the role that further studies of this kind may play in our understanding of these binaries.     

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 late spectral types of cataclysmic variable secondaries

We investigate why the spectral type of most cataclysmic variable (CV) secondaries is significantly later than that of a zero-age main-sequence (ZAMS) star with the same mean density. Using improved stellar input physics, tested against observations of low-mass stars at the bottom of the main sequence, we calculate the secular evolution of CVs with low-mass donors. We consider sequences with different mass transfer rates and with a different degree of nuclear evolution of the donor prior to mass transfer.
Systems near the upper edge of the gap ( P ∼3–6 h) can be reproduced by models with a wide range of mass transfer rates from 1.5×10⁻⁹ M_⊙ yr⁻¹ to 10⁻⁸ M_⊙ yr⁻¹. Evolutionary sequences with a small transfer rate and donors that are substantially evolved off the ZAMS (central hydrogen content 0.05–0.5) reproduce CVs with late spectral types above P ≳6 h. Systems with the most discrepant (late) spectral type should have the smallest donor mass at any given P .
Consistency with the period gap suggests that the mass transfer rate increases with decreasing donor mass for evolved sequences above the period gap. In this case, a single-parameter family of sequences with varying X _c and increasing mass transfer rate reproduces the full range of observed spectral types. This would imply that CVs with such evolved secondaries dominate the CV population.     

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.     

Orbital solution and fundamental parameters of σ Scorpii

The first orbital solution for the spectroscopic pair in the multiple star system σ Scorpii, determined from measurements with the Sydney University Stellar Interferometer, is presented. The primary component is of β Cephei variable type and has been one of the most intensively studied examples of its class. The orbital solution, when combined with radial velocity results found in the literature, yields a distance of  174⁺²³₋₁₈ pc  , which is consistent with, but more accurate than the Hipparcos value. For the primary component we determine  18.4 ± 5.4 M_⊙, −4.12 ± 0.34 mag  and  12.7 ± 1.8 R_⊙  for the mass, absolute visual magnitude and radius, respectively. A B1 dwarf spectral type and luminosity class for the secondary is proposed from the mass determination of  11.9 ± 3.1 M_⊙  and the estimated system age of 10 Myr.     

Modelling the binary progenitor of Supernova 1993J

We have developed a detailed stellar evolution code capable of following the simultaneous evolution of both stars in a binary system, together with their orbital properties. To demonstrate the capabilities of the code, we investigate potential progenitors for the Type IIb Supernova 1993J, which is believed to have been an interacting binary system prior to its primary exploding. We use our detailed binary stellar evolution code to model this system to determine the possible range of primary and secondary masses that could have produced the observed characteristics of this system, with particular reference to the secondary. Using the luminosities and temperatures for both stars (as determined by Maund et al.) and the remaining mass of the hydrogen envelope of the primary at the time of explosion, we find that if mass transfer is 100 per cent efficient, the observations can be reproduced by a system consisting of a  15 M_⊙  primary and a  14 M_⊙  secondary in an orbit with an initial period of 2100 days. With a mass transfer efficiency of 50 per cent, a more massive system consisting of a  17 M_⊙  primary and a  16 M_⊙  secondary in an initial orbit of 2360 days is needed. We also investigate some of the uncertainties in the evolution, including the effects of tidal interaction, convective overshooting and thermohaline mixing.     

Mass transfer from a giant star to a main-sequence companion and its contribution to long-orbital-period blue stragglers

Binary population synthesis shows that mass transfer from a giant star to a main-sequence (MS) companion may account for some observed long-orbital-period blue stragglers. However, little attention is paid to this blue straggler formation scenario as dynamical instability often happens when the mass donor is a giant star. In this paper, we have studied the critical mass ratio, q _c, for dynamically stable mass transfer from a giant star to a MS companion using detailed evolution calculations. The results show that a more evolved star is generally less stable for Roche lobe overflow. Meanwhile,   q _c  almost linearly increases with the amount of the mass and angular momentum lost during mass transfer, but has little dependance on stellar wind. To conveniently use the result, we give a fit of q _c as a function of the stellar radius at the onset of Roche lobe overflow and of the mass-transfer efficiency during the Roche lobe overflow.
To examine the formation of blue stragglers from mass transfer between giants and MS stars, we have performed Monte Carlo simulations with various q _c. The simulations show that some binaries with the mass donor on the first giant branch may contribute to blue stragglers with q _c obtained in this paper but will not from previous q _c. Meanwhile, from our q _c, blue stragglers from the mass transfer between an asymptotic giant branch star and a MS companion may be more numerous and have a wider range of orbital periods than those from the other q _c.     

The case and fate of HD 75767 – neutron star or supernova?

We report the discovery of the nearby  ( d = 24 pc)  HD 75767 as an eight billion year old quadruple system consisting of a distant M dwarf pair, HD 75767 C–D, in orbit around the known short-period   P = 10.25 d  single-lined binary HD 75767 A–B, the primary of which is a solar-like G star. On the reasonable assumption of synchronous orbital rotation as well as rotational and orbital coplanarity for the inner pair, we get   M _B= 0.96 M_⊙  for the unseen HD 75767 B, that is, the case of a massive white dwarf. Upon future evolution, mass transfer towards HD 75767 B will render the   M _A= 0.96 M_⊙  G-type primary, now a turnoff star, to become a helium white dwarf of   M _A∼ 0.33 M_⊙  . Depending on the mass accretion rate, accretion efficiency and composition of the massive white dwarf, this in turn may result in a collapse of HD 75767 B with the formation of a millisecond pulsar, i.e. the creation of a low-mass binary pulsar (LMBP), or, instead, a Type Ia supernova explosion and the complete disruption of HD 75767 B. Irrespective of which scenario applies, we point to the importance of the distant M dwarfs as the likely agents for the formation of the inner, short-period HD 75767 A–B pair, and hence a path that particularly avoids preceding phases of common envelope evolution.     

Introducing a hybrid radiative transfer method for smoothed particle hydrodynamics

A new means of incorporating radiative transfer into smoothed particle hydrodynamics (SPH) is introduced, which builds on the success of two previous methods – the polytropic cooling approximation as devised by Stamatellos et al. and flux-limited diffusion. This hybrid method preserves the strengths of its individual components, while removing the need for atmosphere matching or other boundary conditions to marry optically thick and optically thin regions. The code uses a non-trivial equation of state to calculate temperatures and opacities of SPH particles, which captures the effects of H₂ dissociation, H⁰ ionization, He⁰ and He⁺ ionization, ice evaporation, dust sublimation, molecular absorption, bound-free and free–free transitions and electron scattering. The method is tested in several scenarios, including (i) the evolution of a  0.07 M_⊙  protoplanetary disc surrounding a  0.5 M_⊙  star; (ii) the collapse of a  1 M_⊙  protostellar cloud and (iii) the thermal relaxation of temperature fluctuations in a static homogeneous sphere.     

Optical spectroscopy of the quiescent counterpart to EXO 0748−676

We present phase resolved optical spectroscopy and X-ray timing of the neutron star X-ray binary EXO 0748−676 after the source returned to quiescence in the autumn of 2008. The X-ray light curve displays eclipses consistent in orbital period, orbital phase and duration with the predictions and measurements before the return to quiescence. Hα and He  i emission lines are present in the optical spectra and show the signature of the orbit of the binary companion, placing a lower limit on the radial velocity semi-amplitude of   K ₂ > 405 km s⁻¹  . Both the flux in the continuum and the emission lines show orbital modulations, indicating that we observe the hemisphere of the binary companion that is being irradiated by the neutron star. Effects due to this irradiation preclude a direct measurement of the radial velocity semi-amplitude of the binary companion; in fact, no stellar absorption lines are seen in the spectrum. Nevertheless, our observations place a stringent lower limit on the neutron star mass of   M ₁ > 1.27 M_⊙  . For the canonical neutron star mass of   M ₁= 1.4 M_⊙  , the mass ratio is constrained to  0.075 < q < 0.105  .     

Single and binary evolution of Population III stars and their supernova explosions

We present stellar evolution calculations for Population III stars for both single- and binary-star evolutions. Our models include 10- and  16.5-M_⊙  single stars and a  10-M_⊙  model star that undergoes an episode of accretion resulting in a final mass of  16.1 M_⊙  . For comparison, we present the evolution of a solar heavy element abundance model. We use the structure from late-stage evolution models to calculate simulated supernova light curves. Light curve comparisons are made between accretion and non-accretion progenitor models, and models for single-star evolution of comparable masses. Where possible, we make comparisons to previous works. Similar investigations have been carried out, but primarily for solar or near-solar heavy metal abundance stars and not including both the evolution and the supernova explosions in one work.     

The Jeans mass and the origin of the knee in the IMF

We use numerical simulations of the fragmentation of a  1000 M_⊙  molecular cloud and the formation of a stellar cluster to study how the initial conditions for star formation affect the resulting initial mass function (IMF). In particular, we are interested in the relation between the thermal Jeans mass in a cloud and the knee of the IMF, i.e. the mass separating the region with a flat IMF slope from that typified by a steeper, Salpeter-like, slope. In three isothermal simulations with   M _Jeans= 1, 2  and  5 M_⊙  , the number of stars formed, at comparable dynamical times, scales roughly with the number of initial Jeans masses in the cloud. The mean stellar mass also increases (though less than linearly) with the initial Jeans mass in the cloud. It is found that the IMF in each case displays a prominent knee, located roughly at the mass scale of the initial Jeans mass. Thus clouds with higher initial Jeans masses produce IMFs which are shallow to higher masses. This implies that a universal IMF requires a physical mechanism that sets the Jeans mass to be near  1 M_⊙  . Simulations including a barotropic equation of state as suggested by Larson, with cooling at low densities followed by gentle heating at higher densities, are able to produce realistic IMFs with the knee located at  ≈1 M_⊙  , even with an initial   M _Jeans= 5 M_⊙  . We therefore suggest that the observed universality of the IMF in the local Universe does not require any fine tuning of the initial conditions in star forming clouds but is instead imprinted by details of the cooling physics of the collapsing gas.     

Optical spectroscopy and photometry of SAX J1808.4−3658 in outburst

We present phase resolved optical spectroscopy and photometry of V4580 Sagittarii, the optical counterpart to the accretion powered millisecond pulsar SAX J1808.4−3658, obtained during the 2008 September/October outburst. Doppler tomography of the N  iii λ4640.64 Bowen blend emission line reveals a focused spot of emission at a location consistent with the secondary star. The velocity of this emission occurs at  324 ± 15 km s⁻¹  ; applying a ' K -correction', we find the velocity of the secondary star projected on to the line of sight to be  370 ± 40 km s⁻¹  . Based on existing pulse timing measurements, this constrains the mass ratio of the system to be  0.044^+0.005_−0.004  , and the mass function for the pulsar to be  0.44^+0.16_−0.13 M_⊙  . Combining this mass function with various inclination estimates from other authors, we find no evidence to suggest that the neutron star in SAX J1808.4−3658 is more massive than the canonical value of  1.4 M_⊙  . Our optical light curves exhibit a possible superhump modulation, expected for a system with such a low mass ratio. The equivalent width of the Ca  ii H and K interstellar absorption lines suggest that the distance to the source is ∼2.5 kpc. This is consistent with previous distance estimates based on type-I X-ray bursts which assume cosmic abundances of hydrogen, but lower than more recent estimates which assume helium-rich bursts.     

The possible companions of young radio pulsars

We discuss the formation of pulsars with massive companions in eccentric orbits. We demonstrate that the probability for a non-recycled radio pulsar to have a white dwarf as a companion is comparable to that of having an old neutron star as a companion. Special emphasis is given to PSR B1820−11 and PSR B2303+46. Based on population synthesis calculations we argue that PSR B1820−11 and PSR B2303+46 could very well be accompanied by white dwarfs with mass ≳1.1 M_⊙. For PSR B1820−11, however, we cannot exclude the possibility that its companion is a main-sequence star with a mass between ∼0.7 M_⊙ and ∼5 M_⊙.     

Stellar-mass black hole binaries as ultraluminous X-ray sources

Ultraluminous X-ray sources (ULXs) with   L _x > 10³⁹ erg s⁻¹  have been discovered in great numbers in external galaxies with ROSAT , Chandra and XMM-Newton . The central question regarding this important class of sources is whether they represent an extension in the luminosity function of binary X-ray sources containing neutron stars and stellar-mass black holes (BHs), or a new class of objects, e.g. systems containing intermediate-mass BHs  (100–1000 M_⊙)  . We have carried out a theoretical study to test whether a large fraction of the ULXs, especially those in galaxies with recent star formation activity, can be explained with binary systems containing stellar-mass BHs. To this end, we have applied a unique set of binary evolution models for BH X-ray binaries, coupled to a binary population synthesis code, to model the ULXs observed in external galaxies. We find that for donor stars with initial masses  ≳10 M_⊙  the mass transfer driven by the normal nuclear evolution of the donor star is sufficient to potentially power most ULXs. This is the case during core hydrogen burning and, to an even more pronounced degree, while the donor star ascends the giant branch, although the latter phases last only ∼5 per cent of the main-sequence phase. We show that with only a modest violation of the Eddington limit, e.g. a factor of ∼10, both the numbers and properties of the majority of the ULXs can be reproduced. One of our conclusions is that if stellar-mass BH binaries account for a significant fraction of ULXs in star-forming galaxies, then the rate of formation of such systems is  ∼3 × 10⁻⁷ yr⁻¹  normalized to a core-collapse supernova rate of 0.01 yr⁻¹.     

On the formation and evolution of black hole binaries

We present the results of a systematic study of the formation and evolution of binaries containing black holes and normal-star companions with a wide range of masses. We first reexamine the standard formation scenario for close black hole binaries, where the progenitor system, a binary with at least one massive component, experienced a common-envelope phase and where the spiral-in of the companion in the envelope of the massive star caused the ejection of the envelope. We estimate the formation rates for different companion masses and different assumptions about the common-envelope structure and other model parameters. We find that black hole binaries with intermediate- and high-mass secondaries can form for a wide range of assumptions, while black hole binaries with low-mass secondaries can only form with apparently unrealistic assumptions (in agreement with previous studies).
We then present detailed binary evolution sequences for black hole binaries with secondaries of 2 to 17 M_⊙ and demonstrate that in these systems the black hole can accrete appreciably even if accretion is Eddington-limited (up to 7 M_⊙ for an initial black hole mass of 10 M_⊙) and that the black holes can be spun up significantly in the process. We discuss the implications of these calculations for well-studied black hole binaries (in particular GRS 1915+105) and ultraluminous X-ray sources of which GRS 1915+105 appears to represent a typical Galactic counterpart. We also present a detailed evolutionary model for Cygnus X-1, a massive black hole binary, which suggests that at present the system is most likely in a wind mass-transfer phase following an earlier Roche-lobe overflow phase. Finally, we discuss how some of the assumptions in the standard model could be relaxed to allow the formation of low-mass, short-period black hole binaries, which appear to be very abundant in nature.