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.     

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.     

Supernova kicks and misaligned Be star binaries

Be stars are rapidly spinning B stars surrounded by an outflowing disc of gas in Keplerian rotation. Be star/X-ray binary systems contain a Be star and a neutron star. They are found to have non-zero eccentricities and there is evidence that some systems have a misalignment between the spin axis of the star and the spin axis of the binary orbit. The eccentricities in these systems are caused by a kick to the neutron star during the supernova that formed it. Such kicks would also give rise to misalignments. In this paper, we investigate the extent to which the same kick distribution can give rise to both the observed eccentricity distribution and the observed misalignments. We find that a Maxwellian distribution of velocity kicks with a low velocity dispersion,  σ_k≈ 15 km s⁻¹  , is consistent with the observed eccentricity distribution but is hard to reconcile with the observed misalignments, typically   i ≥ 25°  . Alternatively, a higher velocity kick distribution,  σ_k= 265 km s⁻¹  , is consistent with the observed misalignments but not with the observed eccentricities, unless post-supernova circularization of the binary orbits has taken place. We discuss briefly how this might be achieved.     

The enigma of the oldest 'nova': the central star and nebula of CK Vul

CK Vul is classified as, amongst others, the slowest known nova, a hibernating nova or a very late thermal pulse object. Following its eruption in ad 1670, the star remained visible for 2 yr. A 15-arcsec nebula was discovered in the 1980s, but the star itself has not been detected since the eruption. We here present radio images which reveal a 0.1-arcsec radio source with a flux of 1.5 mJy at 5 GHz. Deep Hα images show a bipolar nebula with a longest extension of 70 arcsec, with the previously known compact nebula at its waist. The emission-line ratios show that the gas is shock-ionized, at velocities  >100 km s⁻¹  . Dust emission yields an envelope mass of  ∼5 × 10⁻² M_⊙  . Echelle spectra indicate outflow velocities up to 360 km s⁻¹. From a comparison of images obtained in 1991 and 2004 we find evidence for expansion of the nebula, consistent with an origin in the 1670 explosion; the measured expansion is centred on the radio source. No optical or infrared counterpart is found at the position of the radio source. The radio emission is interpreted as thermal free–free emission from gas with   T _e∼ 10⁴ K  . The radio source may be due to a remnant circumbinary disc, similar to those seen in some binary post-AGB stars. We discuss possible classifications of this unique outburst, including that of a sub-Chandrasekhar mass supernova, a nova eruption on a cool, low-mass white dwarf or a thermal pulse induced by accretion from a circumbinary disc.     

The violent past of Cygnus X-2

Cygnus X-2 appears to be the descendant of an intermediate-mass X-ray binary (IMXB). Using Mazzitelli's stellar code we compute detailed evolutionary sequences for the system and find that its prehistory is sensitive to stellar input parameters, in particular the amount of core overshooting during the main-sequence phase. With standard assumptions for convective overshooting a case B mass transfer starting with a 3.5-M_⊙ donor star is the most likely evolutionary solution for Cygnus X-2. This makes the currently observed state rather short-lived, of order 3 Myr, and requires a formation rate > 10⁻⁷–10⁻⁶ yr⁻¹ of such systems in the Galaxy. Our calculations show that neutron star IMXBs with initially more massive donors (≳4 M_⊙) encounter a delayed dynamical instability; they are unlikely to survive this rapid mass transfer phase. We determine limits for the age and initial parameters of Cygnus X-2 and calculate possible dynamical orbits of the system in a realistic Galactic potential, given its observed radial velocity. We find trajectories which are consistent with a progenitor binary on a circular orbit in the Galactic plane inside the solar circle that received a kick velocity ≤200 km s⁻¹ at the birth of the neutron star. The simulations suggest that about 7 per cent of IMXBs receiving an arbitrary kick velocity from a standard kick velocity spectrum would end up in an orbit similar to Cygnus X-2, while about 10 per cent of them reach yet larger Galactocentric distances.     

Hyperfast pulsars as the remnants of massive stars ejected from young star clusters

Recent proper motion and parallax measurements for the pulsar PSR B1508+55 indicate a transverse velocity of  ∼1100 km s⁻¹  , which exceeds earlier measurements for any neutron star. The spin-down characteristics of PSR B1508+55 are typical for a non-recycled pulsar, which implies that the velocity of the pulsar cannot have originated from the second supernova disruption of a massive binary system. The high velocity of PSR B1508+55 can be accounted for by assuming that it received a kick at birth or that the neutron star was accelerated after its formation in the supernova explosion. We propose an explanation for the origin of hyperfast neutron stars based on the hypothesis that they could be the remnants of a symmetric supernova explosion of a high-velocity massive star which attained its peculiar velocity (similar to that of the pulsar) in the course of a strong dynamical three- or four-body encounter in the core of dense young star cluster. To check this hypothesis, we investigated three dynamical processes involving close encounters between: (i) two hard massive binaries, (ii) a hard binary and an intermediate-mass black hole (IMBH) and (iii) a single stars and a hard binary IMBH. We find that main-sequence O-type stars cannot be ejected from young massive star clusters with peculiar velocities high enough to explain the origin of hyperfast neutron stars, but lower mass main-sequence stars or the stripped helium cores of massive stars could be accelerated to hypervelocities. Our explanation for the origin of hyperfast pulsars requires a very dense stellar environment of the order of  10⁶– 10⁷ stars pc⁻³  . Although such high densities may exist during the core collapse of young massive star clusters, we caution that they have never been observed.     

Time-series Paschen-β spectroscopy of SU Aurigae

We present time-series echelle spectra of the Paβ line of the T Tauri star SU Aur, observed over three consecutive nights. The line shows strong variability (∼10 per cent) over the velocity range 100–420 km s⁻¹ in the red broad absorption component, and weaker variability (∼2 per cent) over the velocity range  −200–0 km s⁻¹  in the blue wing. The variability in the velocity range  −200–0 km s⁻¹  is correlated with that in  200–400 km s⁻¹  , and the variability in these velocity ranges anticorrelates with that in  0–100 km s⁻¹  . The mean spectrum from the second night shows the suggestion of a blueshifted absorption component at about  −150 km s⁻¹  , similar to that found in the Hα and Hβ lines. We find the position of the subpeak in the red absorption component changes steadily with time, and its motion modulates at half the rotational period. We also find that the modulation of the line equivalent width is possibly associated with a half and a third of the rotational period, which is consistent with the surface Doppler images of SU Aur. Radiative transfer models of a rotationally modulated Paβ line, produced in the shock-heated magnetospheric accretion flow, are also presented. Models with a magnetic dipole offset reproduce the overall characteristics of the observed line variability, including the line equivalent width and the motion of the subpeak in the red absorption trough.     

Inflow and outflow from the accretion disc of the microquasar SS 433: UKIRT spectroscopy

A succession of near-infrared (near-IR) spectroscopic observations, taken nightly throughout an entire cycle of SS 433's orbit, reveal (i) the persistent signature of SS 433's accretion disc, having a rotation speed of  ∼500 km s⁻¹  , (ii) the presence of circumbinary disc recently discovered at optical wavelengths by Blundell, Bowler & Schmidtobreick (2008) and (iii) a much faster outflow than has previously been measured for the disc wind, with a terminal velocity of  ∼1500 km s⁻¹  . The increased wind terminal velocity results in a mass-loss rate of  ∼10⁻⁴ M_⊙ yr⁻¹  . These, together with the newly (upwardly) determined masses for the components of the SS 433 system, result in an accurate diagnosis of the extent to which SS 433 has super-Eddington flows. Our observations imply that the size of the companion star is comparable with the semiminor axis of the orbit which is given by     , where e is the eccentricity. Our relatively spectral resolution at these near-IR wavelengths has enabled us to deconstruct the different components that comprise the Brackett-γ (Brγ) line in this binary system, and their physical origins. With this line being dominated throughout our series of observations by the disc wind, and the accretion disc itself being only a minority (∼15 per cent) contribution, we caution against use of the unresolved Brγ line intensity as an 'accretion signature' in X-ray binaries or microquasars in any quantitative way.     

The circumstellar medium of the peculiar supernova SN1997ab

We report the detection of the slow-moving wind into which the compact supernova remnant SN 1997ab is expanding. Echelle spectroscopy provides clear evidence for a well-resolved narrow (full width at zero intensity, FWZI ∼180 km s⁻¹) P Cygni profile, both in Hα and Hβ, superimposed on the broad emission lines of this compact supernova remnant. From theoretical arguments we know that the broad and strong emission lines imply a circumstellar density ( n  ≥ 10⁷ cm⁻³). This, together with our detection, implies a massive and slow stellar wind experienced by the progenitor star shortly prior to the explosion.     

Radial velocities of pulsating subdwarf B stars: KPD 2109+4401 and PB 8783

High-speed spectroscopy of two pulsating subdwarf B stars, KPD 2109+4401 and PB 8783, is presented. Radial motions are detected with the same frequencies as reported from photometric observations and with amplitudes of ∼2 km s⁻¹ in two or more independent modes. These represent the first direct observations of surface motion arising from multimode non-radial oscillations in subdwarf B stars. In the case of the sdB+F binary PB 8783, the velocities of both components are resolved; high-frequency oscillations are found only in the sdB star and not the F star. There also appears to be evidence for mutual motion of the binary components. If confirmed, it implies that the F-type companion is ≳1.2 times more massive than the sdB star, while the amplitude of the F-star acceleration over 4 h would constrain the orbital period to lie between 0.5 and 3.2 d.     

The mass of the neutron star in Centaurus X-3

We report new radial velocity observations of V779 Cen, the optical companion to the X-ray pulsar Cen X-3. Two sets of results at two epochs yield very different radial velocity amplitudes. We demonstrate there are problems with the first set, not least that they are incompatible with the observed duration of the X-ray eclipse for all inclination angles. The anomalously high radial velocities are probably a result of changes in the outflow behaviour of the companion star. Although there is no reason to doubt the results from the second epoch when viewed in isolation, given the anomalous radial velocities of the first epoch, they must be treated with caution. Using these data, the semi-amplitude of the resulting radial velocity curve is found to be 24.4±4.1 km s⁻¹. Given the accurately measured semi-amplitude of the orbit of the pulsar, 414.3±0.9 km s⁻¹, the mass ratio of the system is 0.059±0.010. The inclination of the system is found to be 702±27, assuming that the optical component fills its Roche lobe, and that the system is in synchronous rotation. Hence the mass of the neutron star is 1.21±0.21 M_⊙, and the mass of the optical companion is 20.5±0.7 M_⊙. This is a smaller uncertainty than previously reported values, and is consistent with the canonical neutron star mass of 1.4 M_⊙.
In addition, we use our spectra to determine the spectral class of V779 Cen to be O6-7II-III.     

Magnetized relativistic jets and long-duration GRBs from magnetar spin-down during core-collapse supernovae

We use ideal axisymmetric relativistic magnetohydrodynamic simulations to calculate the spin-down of a newly formed millisecond,   B ∼ 10¹⁵ G  , magnetar and its interaction with the surrounding stellar envelope during a core-collapse supernova (SN) explosion. The mass, angular momentum and rotational energy lost by the neutron star are determined self-consistently given the thermal properties of the cooling neutron star's atmosphere and the wind's interaction with the surrounding star. The magnetar drives a relativistic magnetized wind into a cavity created by the outgoing SN shock. For high spin-down powers  (∼10⁵¹–10⁵² erg s⁻¹)  , the magnetar wind is superfast at almost all latitudes, while for lower spin-down powers  (∼10⁵⁰ erg s⁻¹)  , the wind is subfast but still super-Alfvénic. In all cases, the rates at which the neutron star loses mass, angular momentum and energy are very similar to the corresponding free wind values (≲30 per cent differences), in spite of the causal contact between the neutron star and the stellar envelope. In addition, in all cases that we consider, the magnetar drives a collimated  (∼5–10°)  relativistic jet out along the rotation axis of the star. Nearly all of the spin-down power of the neutron star escapes via this polar jet, rather than being transferred to the more spherical SN explosion. The properties of this relativistic jet and its expected late-time evolution in the magnetar model are broadly consistent with observations of long duration gamma-ray bursts (GRBs) and their associated broad-lined Type Ic SN.     

A statistical study of 233 pulsar proper motions

We present and analyse a catalogue of 233 pulsars with proper motion measurements. The sample contains a wide variety of pulsars including recycled objects and those associated with globular clusters or supernova remnants. After taking the most precise proper motions for those pulsars for which multiple measurements are available, the majority of the proper motions (58 per cent) are derived from pulsar timing methods, 41 per cent using interferometers and the remaining 1 per cent using optical telescopes. Many of the one-dimensional (1D) and two-dimensional (2D) speeds (referring to speeds measured in one coordinate only and the magnitudes of the transverse velocities, respectively) derived from these measurements are somewhat lower than earlier estimates because of the use of the most recent electron density model in determining pulsar distances. The mean 1D speeds for the normal and recycled pulsars are 152(10) and 54(6) km s⁻¹, respectively. The corresponding mean 2D speeds are 246(22) and 87(13) km s⁻¹. PSRs B2011+38 and B2224+64 have the highest inferred 2D speeds of  ∼1600 km s⁻¹  . We study the mean speeds for different subsamples and find that, in general, they agree with previous results. Applying a novel deconvolution technique to the sample of 73 pulsars with characteristic ages less than 3 Myr, we find the mean three-dimensional (3D) pulsar birth velocity to be 400(40) km s⁻¹. The distribution of velocities is well described by a Maxwellian distribution with  1D rms σ= 265 km s⁻¹  . There is no evidence for a bimodal velocity distribution. The proper motions for PSRs B1830−08 and B2334+61 are consistent with their proposed associations with the supernova remnants W41 and G114.3+0.3, respectively.     

Spectroscopic solution of the star QX Cas

High-resolution spectroscopic observations around the Hα line of the binary star QX Cas covering the whole orbital period are presented. Our radial velocity solution, the first ever determined, requires an eccentric orbit with the following orbital parameters: eccentricity,   e = 0.22 ± 0.01  ; longitude of periastron,  ω= 45°± 5°  ; semi-amplitudes of the radial velocity curves of the primary and secondary stars,   K ₁ sin  i = 125.8 ± 0.9 km s⁻¹  and   K ₂ sin  i = 144.8 ± 1.1 km s⁻¹  ; gamma velocity,   V ₀= 65.1 ± 0.5 km s⁻¹  ; and mass ratio,   q = 0.869 ± 0.013  . The corresponding lower limits of the masses of the components and their separation are         , and   a sin  i = 31.34 ± 0.48 R_⊙  .     

Radial velocity study of the pulsating subdwarf B star and possible Type Ia supernova progenitor KPD 1930+2752

We report the results of a high-time-resolution radial velocity study of the subdwarf B star and possible Type Ia supernova progenitor KPD 1930+2752. There were no significant peaks in the power spectrum of the velocity curve above our detection limit, about 4 km s⁻¹, at the frequencies where peaks arising from pulsation were present in the photometric data of previous researchers. We report an orbital velocity amplitude,  348.5±1 km s^-1  , in agreement with that reported by previous investigators. We find an orbital period of   P =0.095 093 08±0.000 000 15 d  based on our data and the ephemeris of Maxted et al.     

The high-velocity molecular stream in W51B: a ring structure with compact H ii regions

We present ¹³ CO J  = 1 − 0 line observations of the H  ii region complex W51B located in the high-velocity (HV) stream. These observations reveal a filamentary and clumpy structure in the molecular gas. The mean local standard of rest (LSR) velocity ∼ + 65 km s⁻¹ of the molecular gas in this region is greater than the maximum velocities allowed by kinematic Galactic rotation curves. The size and mass of the molecular clouds are ∼ 48 × 17 pc² and ∼ 2.4 × 10⁵ M⊙ respectively. In a position–velocity diagram, molecular gas in the southern part comprises a redshifted ring structure with v _LSR=+ 60 to +73 km s⁻¹. The velocity gradient of this ring is ∼ 0.5 km s⁻¹ pc⁻¹, and the mass is ∼ 6.2 × 10⁴ M⊙. If we assume that the ring is expanding with a uniform velocity, the expansion velocity, radius and kinetic energy are ∼ 7 km s⁻¹, ∼ 13 pc and ∼ 3.0 × 10 ⁴⁹ erg respectively. The kinetic energy and mass spectrum of the ring could be explained by an expanding cylindrical cloud with a centrally condensed mass distribution. The locations of two compact H  ii regions, G49.0−0.3 and G48.9−0.3, coincide with the two molecular clumps in this ring. We discuss star formation, and the mechanism that produced the ring structure.     

A source of high-velocity white dwarfs

We investigate whether the recently observed population of high-velocity white dwarfs can be derived from a population of binaries residing initially within the thin disc of the Galaxy. In particular, we consider binaries where the primary is sufficiently massive to explode as a Type II supernova. A large fraction of such binaries are broken up when the primary then explodes as a supernova, owing to the combined effects of the mass loss from the primary and the kick received by the neutron star on its formation. For binaries where the primary evolves to fill its Roche lobe, mass transfer from the primary leads to the onset of a common envelope phase during which the secondary and the core of the primary spiral together as the envelope is ejected. Such binaries are the progenitors of X-ray binaries if they are not broken up when the primary explodes. For those systems that are broken up, a large number of the secondaries receive kick velocities ∼100–200 km s⁻¹ and subsequently evolve into white dwarfs. We compute trajectories within the Galactic potential for this population of stars and relate the birth rate of these stars over the entire Galaxy to those seen locally with high velocities relative to the local standard of rest (LSR) . We show that for a reasonable set of assumptions concerning the Galactic supernova rate and the binary population, our model produces a local number density of high-velocity white dwarfs compatible with that inferred from observations. We therefore propose that a population of white dwarfs originating in the thin disc may make a significant contribution to the observed population of high-velocity white dwarfs.     

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  .     

Kinematics of massive stars in the Small Magellanic Cloud

We present radial velocities for 2045 stars in the Small Magellanic Cloud (SMC), obtained from the 2dF survey by Evans et al. The great majority of these stars are of OBA type, tracing the dynamics of the young stellar population. Dividing the sample into ad hoc 'bar' and 'wing' samples (north and south, respectively, of the line:  δ=−77°50'+[4α]'  , where α is the right ascension in minutes of time) we find that the velocities in the SMC bar show a gradient of 26.3 ± 1.6 km s⁻¹ deg⁻¹ at a position angle of 126°± 4°. The derived gradient in the bar is robust to the adopted line of demarcation between the two samples. The largest redshifts are found in the SMC wing, in which the velocity distribution appears distinct from that in the bar, most probably a consequence of the interaction between the Magellanic Clouds that is predicted to have occurred 0.2 Gyr ago. The mean velocity for all stars in the sample is +172.0 ± 0.2 km s⁻¹ (redshifted by ∼20 km s⁻¹ when compared to published results for older populations), with a velocity dispersion of 30 km s⁻¹.     

The mass of the white dwarf in the old nova BT Mon

We present spectrophotometry of the eclipsing old nova BT Mon (Nova Mon 1939). By detecting weak absorption features from the secondary star, we find its radial velocity semi-amplitude to be K _R = 205 ± 5 km s⁻¹ and its rotational velocity to be v  sin  i  = 138 ± 5 km s⁻¹. We also measure the radial velocity semi-amplitude of the primary star to be K _R = 170 ± 10 km s⁻¹. From these parameters we obtain a mass of 1.04 ± 0.06 M⊙ for the white dwarf primary star and a mass of 0.87 ⊙ 0.06 M⊙ for the G8 V secondary star. The inclination of the system is found to be 82^°_.2 ± 32^°_.2 and we estimate that the system lies at a distance of 1700 ± 300 pc. The high mass of the white dwarf and our finding that BT Mon was probably a fast nova together constitute a new piece of evidence in favour of the thermonuclear runaway model of classical nova outbursts. The emission lines are single-peaked throughout the orbital cycle, showing absorption around phase 0.5, high-velocity S-wave components and large phase offsets in their radial velocity curves. In each of these respects, BT Mon is similar to the SW Sex stars. We also find quasi-periodic flaring in the trailed spectra, which makes BT Mon a candidate intermediate polar.