Photometry of four binary subdwarf B stars and the nature of their unseen companion stars

We present light curves of four binary subdwarf B stars (sdB), Ton 245, Feige 11, PG 1432+159 and PG 1017−086. We also present new spectroscopic data for PG 1017−086 from which we derive its orbital period,   P =0.073 d  , and the mass function,   f _m=0.0010±0.0002 M_⊙.  This is the shortest period for an sdB binary measured to date. The values of P and f _m for the other sdB binaries have been published elsewhere. We are able to exclude the possibility that the unseen companion stars to Ton 245, Feige 11 and PG 1432+159 are main-sequence stars or subgiant stars from the absence of a sinusoidal signal, which would be caused by the irradiation of such a companion star, i.e. they show no reflection effect. The unseen companion stars in these binaries are likely to be white dwarf stars. In contrast, the reflection effect in PG 1017−086 is clearly seen. The lack of eclipses in this binary combined with other data suggests that the companion is a low-mass M-dwarf or, perhaps, a brown dwarf.     

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.     

The formation of the black hole in the X-ray binary system V404 Cyg

Using new and archival radio data, we have measured the proper motion of the black hole X-ray binary V404 Cyg to be  9.2 ± 0.3 mas yr⁻¹  . Combined with the systemic radial velocity from the literature, we derive the full three-dimensional heliocentric space velocity of the system, which we use to calculate a peculiar velocity in the range 47–102 km s⁻¹, with a best-fitting value of 64 km s⁻¹. We consider possible explanations for the observed peculiar velocity and find that the black hole cannot have formed via direct collapse. A natal supernova is required, in which either significant mass  (∼11 M_⊙)  was lost, giving rise to a symmetric Blaauw kick of up to ∼65 km s⁻¹, or, more probably, asymmetries in the supernova led to an additional kick out of the orbital plane of the binary system. In the case of a purely symmetric kick, the black hole must have been formed with a mass  ∼9 M_⊙  , since when it has accreted  0.5–1.5 M_⊙  from its companion.     

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

Cooling curves and initial models for low-mass white dwarfs (<0.25 M⊙) with helium cores

We present a detailed calculation of the evolution of low-mass (<0.25 M_⊙) helium white dwarfs. These white dwarfs (the optical companions to binary millisecond pulsars) are formed via long-term, low-mass binary evolution. After detachment from the Roche lobe, the hot helium cores have a rather thick hydrogen layer with mass between 0.01 and 0.06 M_⊙. As a result of mixing between the core and outer envelope, the surface hydrogen content ( X _surf) is 0.5–0.35 , depending on the initial value of the heavy element Z and the initial secondary mass. We found that the majority of our computed models experience one or two hydrogen shell flashes. We found that the mass of the helium dwarf in which the hydrogen shell flash occurs depends on the chemical composition. The minimum helium white dwarf mass in which a hydrogen flash takes place is 0.213 M_⊙ ( Z =0.003), 0.198 M_⊙ ( Z =0.01), 0.192 M_⊙ ( Z =0.02) or 0.183 M_⊙ ( Z =0.03). The duration of the flashes (independent of chemical composition) is between a few ×10⁶ and a few ×10⁷ yr. In several flashes the white dwarf radius will increase so much that it forces the model to fill its Roche lobe again. Our calculations show that the cooling history of the helium white dwarf depends dramatically on the thickness of the hydrogen layer. We show that the transition from a cooling white dwarf with a temporarily stable hydrogen-burning shell to a cooling white dwarf in which almost all residual hydrogen is lost in a few thermal flashes (via Roche lobe overflow) occurs between 0.183 and 0.213 M_⊙ (depending on the heavy element value).     

The imprint of massive black hole formation models on the LISA data stream

The formation, merging and accretion history of massive black holes (MBHs) along the hierarchical build-up of cosmic structures leaves a unique imprint on the background of gravitational waves (GWs) at mHz frequencies. We study here, by means of dedicated simulations of black hole build-up, the possibility of constraining different models of black hole cosmic evolution using future GW space-borne missions, such as LISA . We consider two main scenarios for black hole formation, namely, one where seeds are light (  ≃10² M_⊙  , remnant of Population III stars) and one where seeds are heavy (  ≳10⁴ M_⊙  , direct collapse). In all the models we have investigated, MBH binary coalescences do not produce a stochastic GW background, but rather, a set of individual resolved events. Detection of several hundreds merging events in a 3-yr LISA mission will be the sign of a heavy seed scenario with efficient formation of black hole seeds in a large fraction of high-redshift haloes. At the other extreme, a low event rate, about a few tens in 3 yr, is peculiar of scenarios where either the seeds are light, and many coalescences do not fall into the LISA band, or seeds are massive, but rare. In this case a decisive diagnostic is provided by the shape of the mass distribution of detected events. Light binaries  ( m < 10⁴ M_⊙)  are predicted in a fairly large number in Population III remnant models, but are totally absent in direct collapse models. Finally, a further, helpful diagnostic of black hole formation models lies in the distribution of the mass ratios in binary coalescences. While heavy seed models predict that most of the detected events involve equal-mass binaries, in the case of light seeds, mass ratios are equally distributed in the range 0.1–1.     

Two-Micron All-Sky Survey J01542930+0053266: a new eclipsing M dwarf binary system

We report on Two-Micron All-Sky Survey (2MASS) J01542930+0053266, a faint eclipsing system composed of two M dwarfs. The variability of this system was originally discovered during a pilot study of the 2MASS Calibration Point Source Working Data base. Additional photometry from the Sloan Digital Sky Survey yields an eight-passband light curve from which we derive an orbital period of  2.639 0157 ± 0.000 0016  d. Spectroscopic followup confirms our photometric classification of the system, which is likely composed of M0 and M1 dwarfs. Radial velocity measurements allow us to derive the masses  (M₁= 0.66 ± 0.03 M_⊙; M₂= 0.62 ± 0.03 M_⊙)  and radii  (R₁= 0.64 ± 0.08 R_⊙; R₂= 0.61 ± 0.09 R_⊙)  of the components, which are consistent with empirical mass–radius relationships for low-mass stars in binary systems. We perform Monte Carlo simulations of the light curves which allow us to uncover complicated degeneracies between the system parameters. Both stars show evidence of Hα emission, something not common in early-type M dwarfs. This suggests that binarity may influence the magnetic activity properties of low-mass stars; activity in the binary may persist long after the dynamos in their isolated counterparts have decayed, yielding a new potential foreground of flaring activity for next generation variability surveys.     

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.     

PSR J1012+5307: younger than it looks?

Lorimer et al. have recently reported that the spin-down age (∼7 × 10⁹ yr) of the low-mass binary pulsar PSR J1012+5307 is much higher than the cooling age (3 × 10⁸ yr) of its white dwarf companion. The proposed solutions for this discrepancy are outlined and discussed. In particular, the revised cooling age estimate proposed by Alberts et al. agrees with data from other low-mass binary pulsar systems if a transition to the 'classical' cooling regime occurs between ∼0.14 and ∼0.28 M_⊙. If this transition is excluded, PSR J1012+5307 seems to have finished its accretion phase far from the spin-up line.     

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.     

Gravitational wave emission from a companion black hole in the presence of an accretion disc around a super-massive Kerr black hole

Gravitational wave signal characteristics from a binary black hole system in which the companion moves through the accretion disc of the primary are studied. We chose the primary to be a super-massive  ( M = 10⁸ M_⊙)  Kerr black hole and the companion to be a massive black hole  ( M = 10⁵ M_⊙)  to clearly demonstrate the effects. We show that the drag exerted on the companion by the disc is sufficient to reduce the coalescence time of the binary. The drag is primarily due to the fact that the accretion disc on a black hole deviates from a Keplerian disc and becomes sub-Keplerian due to inner boundary condition on the black hole horizon. We consider two types of accretion rates on to the companion. The companion is deeply immersed inside the disc and it can accrete at the Bondi rate which depends on the instantaneous density of the disc. However, an accretion disc can also form around the smaller black hole and it can accrete at its Eddington rate. Thus, this case is also studied and the results are compared. We find that the effect of the disc will be significant in reducing the coalescence time and one needs to incorporate this while interpreting gravitational wave signals emitted from such a binary system.     

Low-mass binary-induced outflows from asymptotic giant branch stars

A significant fraction of planetary nebulae (PNe) and protoplanetary nebulae (PPNe) exhibit aspherical, axisymmetric structures, many of which are highly collimated. The origin of these structures is not entirely understood, however, recent evidence suggests that many observed PNe harbour binary systems, which may play a role in their shaping. In an effort to understand how binaries may produce such asymmetries, we study the effect of low-mass  (<0.3 M_⊙)  companions (planets, brown dwarfs and low-mass main-sequence stars) embedded into the envelope of a  3.0-M_⊙  star during three epochs of its evolution [red giant branch, asymptotic giant branch (AGB), interpulse AGB]. We find that common envelope evolution can lead to three qualitatively different consequences: (i) direct ejection of envelope material resulting in a predominately equatorial outflow, (ii) spin-up of the envelope resulting in the possibility of powering an explosive dynamo-driven jet and (iii) tidal shredding of the companion into a disc which facilitates a disc-driven jet. We study how these features depend on the secondary's mass and discuss observational consequences.     

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.     

High-resolution optical and infrared spectroscopic observations of Cir X-1

We present new optical and infrared (IR) observations of Cir X-1 taken near apastron. Both sets of spectra show asymmetric emission lines. Archival optical observations show that an asymmetric H α emission line has been in evidence for the past 20 years, although the shape of the line has changed significantly. We present an eccentric ( e ∼0.7–0.9) low-mass binary model, where the system consists of a neutron star orbiting around a (sub)giant companion star of 3–5 M_⊙. We suggest that the broad components of the emission lines arise in a high-velocity, optically thick flow near the neutron star, while the narrow components of the optical and the IR lines arise near the companion star and a heated ejecta shell surrounding the binary respectively. In this model, the velocity of the narrow component reflects the space velocity of the binary; the implied radial velocity (+430 km s⁻¹ after correcting for Galactic rotation) is the highest velocity known for an X-ray binary.     

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 census of nuclear activity of late-type galaxies in the Virgo cluster

The first spectroscopic census of active galactic nuclei (AGNs) associated with late-type galaxies in the Virgo cluster was carried out by observing 213 out of a complete set of 237 galaxies more massive than   M _dyn > 10^8.5 M_⊙  . Among them, 77 are classified as AGNs [including 21 transition objects, 47 low-ionization nuclear emission regions (LINERs) and nine Seyferts] and comprise 32 per cent of the late-type galaxies in Virgo. Due to spectroscopic incompleteness, at most 21 AGNs are missed in the survey, so that the fraction would increase up to 41 per cent. Using corollary near-infrared observations that enable us to estimate galaxy dynamical masses, it is found that AGNs are hosted exclusively in massive galaxies, i.e.   M _dyn≳ 10¹⁰ M_⊙  . Their frequency increases steeply with the dynamical mass from zero at   M _dyn≈ 10^9.5 M_⊙  to virtually 1 at   M _dyn > 10^11.5 M_⊙  . These frequencies are consistent with those of low-luminosity AGNs found in the general field by the Sloan Digital Sky Survey. Massive galaxies that harbour AGNs commonly show conspicuous r -band star-like nuclear enhancements. Conversely, they often, but not necessarily, contain massive bulges. A few well-known AGNs (e.g. M61, M100, NGC 4535) are found in massive Sc galaxies with little or no bulge. The AGN fraction seems to be only marginally sensitive to galaxy environment. We infer the black hole masses using the known scaling relations of quiescent black holes. No black holes lighter than  ∼10⁶ M_⊙  are found active in our sample.     

Circinus X-1: survivor of a highly asymmetric supernova

We have analysed the kinematical parameters of Cir X-1 to constrain the nature of its companion star, the eccentricity of the binary and the pre-supernova parameter space. We argue that the companion is most likely to be a low-mass (≲2.0 M_⊙) unevolved star and that the eccentricity of the orbit is 0.94±0.04. We have evaluated the dynamical effects of the supernova explosion and we find it must have been asymmetric. On average , we find that a kick of ∼740 km s⁻¹ is needed to account for the recently measured radial velocity of +430 km s⁻¹ (Johnston, Fender & Wu) for this extreme system. The corresponding minimum kick velocity is ∼500 km s⁻¹. This is the largest kick needed to explain the motion of any observed binary system. If Cir X-1 is associated with the supernova remnant G321.9-0.3 then we find a limiting minimum age of this remnant of ∼60 000 yr. Furthermore, we predict that the companion star has lost ∼10 per cent of its mass as a result of stripping and ablation from the impact of the supernova shell shortly after the explosion.     

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  .     

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.     

Long Secondary Periods in variable red giants

We present a study of a sample of Large Magellanic Cloud red giants exhibiting Long Secondary Periods (LSPs). We use radial velocities obtained from VLT spectral observations and MACHO and OGLE light curves to examine properties of the stars and to evaluate models for the cause of LSPs. This sample is much larger than the combined previous studies of Hinkle et al. and Wood, Olivier & Kawaler.
Binary and pulsation models have enjoyed much support in recent years. Assuming stellar pulsation, we calculate from the velocity curves that the typical fractional radius change over an LSP cycle is greater than 30 per cent. This should lead to large changes in T _eff that are not observed. Also, the small light amplitude of these stars seems inconsistent with the radius amplitude. We conclude that pulsation is not a likely explanation for the LSPs. The main alternative, physical movement of the star – binary motion – also has severe problems. If the velocity variations are due to binary motion, the distribution of the angle of periastron in our large sample of stars has a probability of  1.4 × 10⁻³  that it comes from randomly aligned binary orbits. In addition, we calculate a typical companion mass of  0.09 M_⊙  . Less than 1 per cent of low-mass main-sequence stars have companions near this mass  (0.06–0.12 M_⊙)  whereas ∼25–50 per cent of low-mass red giants end up with LSPs. We are unable to find a suitable model for the LSPs and conclude by listing their known properties.