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.     

The Flared Disc Project: RXTE and ASCA observations of X 1822−371

We present archival Rossi X-ray Timing Explorer ( RXTE ) and simultaneous Advanced Satellite for Cosmology and Astrophysics ( ASCA ) data of the eclipsing low mass X-ray binary (LMXB) X 1822−371. Our spectral analysis shows that a variety of simple models can fit the spectra relatively well. Of these models, we explore two in detail through phase-resolved fits. These two models represent the case of a very optically thick and a very optically thin corona. While systematic residuals remain at high energies, the overall spectral shape is well approximated. The same two basic models are fitted to the X-ray light curve, which shows sinusoidal modulations interpreted as absorption by an opaque disc rim of varying height. The geometry we infer from these fits is consistent with previous studies: the disc rim reaches out to the tidal truncation radius, while the radius of the corona (approximated as spherical) is very close to the circularization radius. Timing analysis of the RXTE data shows a time-lag from hard to soft consistent with the coronal size inferred from the fits. Neither the spectra nor the light curve fits allow us to rule out either model, leaving a key ingredient of the X 1822−371 puzzle unsolved. Furthermore, while previous studies were consistent with the central object being a 1.4 M_⊙ neutron star, which has been adopted as the best guess scenario for this system, our light curve fits show that a white dwarf or black hole primary can work just as well. Based on previously published estimates of the orbital evolution of X 1822−371, however, we suggest that this system contains either a neutron star or a low mass (≲2.5 M_⊙) black hole and is in a transitional state of duration shortward of 10⁷ yr.     

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⁻¹.     

An 'outside-in' outburst of Aql X–1

We present optical spectroscopy and optical and infrared photometry of the neutron star soft X-ray transient Aql X–1 during its X-ray outburst of 1997 August. By modelling the X-ray, optical and IR light curves, we find a 3-d delay between the IR and X-ray rise times, analogous to the UV–optical delay seen in dwarf novae outbursts and black hole X-ray transients. We interpret this delay as the signature of an 'outside-in' outburst, in which a thermal instability in the outer disc propagates inward. This outburst is the first of this type definitively identified in a neutron star X-ray transient.     

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.     

An XMM–Newton observation of the neutron star X-ray transient 2S 1803−245 in quiescence

We observed the neutron star X-ray transient 2S 1803−245 in quiescence with the X-ray satellite XMM–Newton , but did not detect it. An analysis of the X-ray bursts observed during the 1998 outburst of 2S 1803−245 gives an upper limit to the distance of ≤7.3 kpc, leading to an upper limit on the quiescent 0.5–10 keV X-ray luminosity of  ≤2.8 × 10³² erg s⁻¹  (3σ). Since the expected orbital period of 2S 1803−245 is several hours, this limit is not much higher than those observed for the quiescent black hole transients with similar orbital periods.     

The population of faint transients in the Galactic Centre

BeppoSAX has detected a population of faint transient X-ray sources in the Galactic Centre. I show that a simple irradiated disc picture gives a consistent fit to the properties of this population, and that it probably consists of low-mass X-ray binaries (LMXBs) that have evolved beyond their minimum orbital periods of ∼80 min. Since all post-minimum systems are transient, and neutron star LMXBs are more common than black hole LMXBs in the Galaxy, the majority of these systems should contain neutron stars, as observed. This picture predicts that the Galactic Centre transients should have orbital periods in the range ∼80–120 min, and that most of them should repeat in the next few years. In this case, the total number of post-minimum transients in the Galaxy would be considerably smaller than the usual estimates of its total LMXB population. I discuss possible reasons for this.     

Continuous frequency spectrum of the global hydromagnetic oscillations of a magnetically confined mountain on an accreting neutron star

We compute the continuous part of the ideal-magnetohydrodynamic (ideal-MHD) frequency spectrum of a polar mountain produced by magnetic burial on an accreting neutron star. Applying the formalism developed by Hellsten & Spies, extended to include gravity, we solve the singular eigenvalue problem subject to line-tying boundary conditions. This spectrum divides into an Alfvén part and a cusp part. The eigenfunctions are chirped and anharmonic with an exponential envelope, and the eigenfrequencies cover the whole spectrum above a minimum ω_low. For equilibria with accreted mass  1.2 × 10⁻⁶≲ M _a/M_⊙≲ 1.7 × 10⁻⁴  and surface magnetic fields  10¹¹≲ B _*/G ≲ 10¹³, ω_low  is approximately independent of   B _*  , and increases with M _a. The results are consistent with the Alfvén spectrum excited in numerical simulations with the zeus-mp solver. The spectrum is modified substantially by the Coriolis force in neutron stars spinning faster than ∼100 Hz. The implications for gravitational-wave searches for low-mass X-ray binaries are considered briefly.     

Quasi-periodic oscillations in low-mass X-ray binaries and constraints on the equation of state of neutron star matter

Recently discovered quasi-periodic oscillations in the X-ray brightness of low-mass X-ray binaries are used to derive constraints on the mass of the neutron star component and the equation of state of neutron star matter. The observations are compared with models of rapidly rotating neutron stars which are calculated by means of an exact numerical method in full relativity. For the equations of state we select a broad collection of models representing different assumptions about the many-body structure and the complexity of the composition of superdense matter. The mass constraints differ from their values in the approximate treatment by ∼10 per cent. Under the assumption that the maximum frequency of the quasi-periodic oscillations originates from the innermost stable orbit, the mass of the neutron star is in the range M ∼1.92–2.25 M_⊙. The quasi-periodic oscillation in the Atoll-source 4U 1820−30 in particular is only consistent with equations of state that are rather stiff at high densities, which is explainable, so far, only with pure nucleonic/leptonic composition. This interpretation contradicts the hypothesis that the protoneutron star formed in SN 1987A collapsed to a black hole, since this would demand a maximum neutron star mass below 1.6 M_⊙. The recently suggested identification of quasi-periodic oscillations with frequencies of about 10 Hz with the Lense–Thirring precession of the accretion disc is found to be inconsistent with the models studied in this work, unless it is assumed that the first overtone of the precession is observed.     

Optical spectroscopy and Doppler tomography of Cygnus X-2

We present phase resolved optical spectroscopy and Doppler tomography of V1341 Cygni, the optical counterpart to the neutron star low-mass X-ray binary (LMXB) Cygnus X-2 (Cyg X-2). We derive a radial velocity (RV) curve for the secondary star, finding a projected RV semi-amplitude of   K ₂= 79 ± 3 km s⁻¹  , leading to a mass function of  0.51 ± 0.06 M_⊙, ∼30  per cent lower than the previous estimate. We tentatively attribute the lower value of K ₂ (compared to that obtained by other authors) to variations in the X-ray irradiation of the secondary star at different epochs of observations. The limited phase coverage and/or longer timebase of previous observations may also contribute to the difference in K ₂. Our value for the mass function implies a primary mass of  1.5 ± 0.3 M_⊙  , somewhat lower than previous dynamical estimates, but consistent with the value found by analysis of type-I X-ray bursts from this system. Our Doppler tomography of the broad He  ii λ4686 line reveals that most of the emission from this line is produced on the irradiated face of the donor star, with little emission from the accretion disc. In contrast, the Doppler tomogram of the N  iii λ4640.64 Bowen blend line shows bright emission from near the gas stream/accretion disc impact region, with fainter emission from the gas stream and secondary star. This is the first LMXB for which the Bowen blend is dominated by emission from the gas stream/accretion disc impact region, without comparable emission from the secondary star. This has implications for the interpretation of Bowen blend Doppler tomograms of other LMXBs for which the ephemeris may not be accurately known.     

The ultraviolet line spectrum of the soft X-ray transient XTE J1118+480: a CNO-processed core exposed

We compare ultraviolet (UV) spectra of the recent soft X-ray transients XTE J1118+480 and XTE J1859+226. The emission line strengths in XTE J1118+480 strongly suggest that the accreting material has been CNO processed. We show that this system must have come into contact with a secondary star of about 1.5 M_⊙, and an orbital period ∼15 h, very close to the bifurcation value at which the nuclear and angular momentum loss time-scales are similar. Subsequent evolution to the current period of 4.1 h was driven by angular momentum loss. In passing through a period of 7.75 h the secondary star would have shown essentially normal surface abundances. XTE J1118+480 could thus represent a slightly later evolutionary stage of A0620-00. We briefly discuss the broad Ly α absorption wings in XTE J1118+480.     

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.     

The mass of the neutron star in the binary millisecond pulsar PSR J 1012 + 5307

We have measured the radial velocity variation of the white dwarf secondary in the binary system containing the millisecond pulsar PSR J 1012 + 5307. Combined with the orbital parameters of the radio pulsar, we infer a mass ratio q (≡ M ₁/ M ₂) = 10.5 ± 0.5. Our optical spectroscopy has also allowed us to determine the mass of the white dwarf companion by fitting the spectrum to a grid of DA model atmospheres: we estimate M ₂ = 0.16 ± 0.02 M⊙, and hence the mass of the neutron star is 1.64 ± 0.22 M⊙, where the error is dominated by that of M ₂. The orbital inclination is 52 ± 4°. For an initial neutron star mass of ∼ 1.4 M⊙, only a few tenths of a solar mass at most has been successfully accreted over the lifetime of the progenitor low-mass X-ray binary. If the initial mass of the secondary was ∼ 1 M⊙, our result suggests that the mass transfer may have been non-conservative.     

On the effect of overshooting as predicted by the modelling of the pre-main-sequence evolution of a 2 M⊙ star

We discuss the effects of convective overshooting in the pre-main-sequence (PMS) evolution of intermediate-mass stars, by analysing in detail the early evolution towards the main sequence of a  2 M_⊙  stellar model. These effects can be extremely important in the end of the PMS, when the abundances in CNO elements approach the equilibrium in the centre. We provide a possible physical explanation on why a moderate amount of overshooting produces, as the star approaches the zero-age main-sequence, an extra loop in the evolutionary tracks on the Hertzsprung–Russell diagram.
An interesting feature is that there is a very well defined amount of overshooting (for a given stellar mass and chemical composition) beyond which a loop is produced. For smaller amounts of overshooting such a loop does not take place and the evolutionary tracks are similar to those found in the literature. The amount of overshooting needed to produce the loop decreases with stellar mass.
We discuss the underlining physical reasons for the behaviour predicted by the evolution models and argue that it provides a crucial observational test for convective overshooting in the core of intermediate-mass stars.     

Finding Faint Intermediate-Mass Black Holes in the Radio Band

We discuss the prospects for detecting faint intermediate-mass black holes, such as those predicted to exist in the cores of globular clusters and dwarf spheroidal galaxies. We briefly summarize the difficulties of stellar dynamical searches, then show that recently discovered relations between black hole mass, X-ray luminosity and radio luminosity imply that in most cases, these black holes should be more easily detected in the radio than in the X-rays. Finally, we show upper limits from some radio observations of globular clusters, and discuss the possibility that the radio source in the core of the Ursa Minor dwarf spheroidal galaxy might be a ∼10,000–100,000 M⊙ black hole.     

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  .     

The determination of progenitor masses for planetary nebulae

We have evaluated the likely progenitor masses M _PG of nebulae having elliptical, circular and bipolar morphologies, using observed ratios between the populations of these sources, and deduced central star mass functions. We find that most bipolar nebulae (BPNe) are likely to arise from progenitors having mass M _PG>2.3 M_⊙ and spectral types earlier than A3.2, whilst circular sources are associated with progenitors of mass 1.0 M_⊙< M _PG<1.2 M_⊙ and spectral range G1.9–F7.8 . Elliptical sources arise from intermediate-mass progenitors. The procedures employed to determine these values are relatively insensitive to uncertainties in scaleheights and population ratios, and completely insensitive to uncertainties in the distance scale. They are, however, dependent upon the precise forms adopted for the initial–final and central star mass functions, and we discuss the sensitivity of M _PG to uncertainties in these functions.     

Spectral variability in Cygnus X-3

We model the broad-band X-ray spectrum of Cyg X-3 in all states displayed by this source as observed by the Rossi X-ray Timing Explorer . From our models, we derive for the first time unabsorbed spectral shapes and luminosities for the full range of spectral states. We interpret the unabsorbed spectra in terms of Comptonization by a hybrid electron distribution and strong Compton reflection. We study the spectral evolution and compare with other black hole as well as neutron star sources. We show that a neutron star accretor is not consistent with the spectral evolution as a function of L _E and especially not with the transition to a hard state. Our results point to the compact object in Cyg X-3 being a massive,  ∼30 M_⊙  black hole.     

Black hole mergers: can gas discs solve the 'final parsec' problem?

We compute the effect of an orbiting gas disc in promoting the coalescence of a central supermassive black hole binary. Unlike earlier studies, we consider a finite mass of gas with explicit time dependence: we do not assume that the gas necessarily adopts a steady state or a spatially constant accretion rate, i.e. that the merging black hole was somehow inserted into a pre-existing accretion disc. We consider the tidal torque of the binary on the disc, and the binary's gravitational radiation. We study the effects of star formation in the gas disc in a simple energy feedback framework.
The disc spectrum differs in detail from that found before. In particular, tidal torques from the secondary black hole heat the edges of the gap, creating bright rims around the secondary. These rims do not in practice have uniform brightness either in azimuth or time, but can on average account for as much as 50 per cent of the integrated light from the disc. This may lead to detectable high-photon-energy variability on the relatively long orbital time-scale of the secondary black hole, and thus offer a prospective signature of a coalescing black hole binary.
We also find that the disc can drive the binary to merger on a reasonable time-scale only if its mass is at least comparable with that of the secondary black hole, and if the initial binary separation is relatively small, i.e.   a ₀≲ 0.05  pc. Star formation complicates the merger further by removing mass from the disc. In the feedback model we consider, this sets an effective limit to the disc mass. As a result, binary merging is unlikely unless the black hole mass ratio is ≲0.001. Gas discs thus appear not to be an effective solution to the 'last parsec' problem for a significant class of mergers.     

A comparison of white dwarf mass estimators using the intermediate polar XY Ari

We determine the mass of the white dwarf in the eclipsing intermediate polar XY Ari following the method given in Cropper, Ramsay &38; Wu using a multitemperature bremsstrahlung model. By fitting X-ray spectra from Ginga RXTE and ASCA we find that the mean of the best fits to the data taken using different detectors is M _wd = 1.28 ± 0.04 M⊙. This figure is too high to be consistent with the mass of the white dwarf found by Hellier from X-ray eclipse timings. There are also small systematic differences between the masses derived using different X-ray satellites.