X-ray binary unified analysis: pulse/RV application to Vela X1: GP Velorum

A procedure for simultaneous analysis of multiple kinds of observations is developed for binaries that contain X-ray pulsars. A wide variety of observation types might be included, although we consider only velocity: pulse and light: velocity: pulse cases at present. The model operates with equipotentials and can accommodate non-synchronous rotation and eccentric orbits. The duration condition imposed by observed X-ray eclipses is incorporated as an embedded constraint, so that only solutions consistent with the eclipse duration are found. Relations needed for the method of differential corrections in least-squares solutions are specified, and we apply the Marquardt scheme for improvement of solution convergence. Parameters of Vela X1: GP Velorum are obtained from seven combinations of pulse and radial velocity data. The estimated masses are thereby put into some perspective, especially for the neutron star. The mass of the supergiant ( m  sin³  i ) now ranges only ±5 per cent from the mean of seven results, although the uncertainty in i makes the actual mass range larger. We discuss the determinacy of certain parameters and their historical consistency. For example, the systemic radial velocity differs among spectroscopic data sets by more than 20 km s⁻¹. We show that, contrary to published assertions, the orbital inclinations of high-mass X-ray binaries are essentially indeterminable by the usual methods, although lower limits to inclination may be meaningful for some examples.     

High-velocity stars from decay of small stellar systems

In this paper we present numerical results on the decay of small stellar systems under different initial conditions (multiplicity 3 ≤  N  ≤ 10, and various mass spectra, initial velocities and initial configurations). The numerical treatment uses the CHAIN1 code (Mikkola &38; Aarseth). Particular attention is paid to the distribution of high-velocity escapers: we define these as stars with velocity above 30 km s⁻¹. These numerical experiments confirm that small N -body systems are dynamically unstable and produce cascades of escapers in the process of their decay. It is shown that the fraction of stars that escape from small dense stellar systems with an escape velocity greater than 30 km s⁻¹ is ∼1 per cent for all systems treated here. This relatively small fraction must be considered in relation to the rate of star formation in the Galaxy in small groups: this could explain some moderately high-velocity stars observed in the Galactic disc and possibly some young stars with relatively high metallicity in the thick disc.     

Stream overflow and dwarf nova eruptions

We consider the effects of accretion stream overflow on the viscous dynamics of accretion discs in dwarf novae. If the stream from the secondary star is geometrically thick enough, some fraction of its material can flow over and under the disc. The mass and specific angular momentum of the stream are then deposited not only at the point of collision with the outer disc, but also at those radii in the inner disc with geometric heights that are large enough to intercept the residual stream, or near the radius where the disc has the same specific angular momentum as the stream. The overflowing stream can alter the behaviour of heating fronts and cooling fronts in the disc. If the mass fraction of the overflowing stream is of order tens of per cent, the deposition of mass in the inner parts of the disc is sufficient to change the character of the eruption light curves significantly.     

The destructive effects of binary encounters on red giants in the Galactic Centre

We consider the destructive effects of encounters between binaries and red giant stars in the Galactic Centre. Such encounters may explain the observed depletion of luminous red giants within the central 0.2 pc of the galaxy. We consider encounters involving 2- and 8-M⊙ red giants, and thus span the range of stellar masses contributing to the most luminous red giants observed in the Galactic Centre. To explore the phase space of encounters thoroughly, we simulate 18 × 10³ encounters using a modified four-body code in which the red giant core and components of the binary are treated as point masses, and where the envelope configuration is assumed to remain static throughout the encounter. We then rerun a small number of encounters with a smoothed particle hydrodynamics (SPH) code to confirm the reliability of conclusions drawn from the four-body runs. We see two possible pathways to red giant destruction. A large fraction of encounters lead to the formation of common-envelope systems, where two compact objects (drawn from the red giant core and the components of the original binary) form a binary within a common gaseous envelope, whilst the third body is ejected. The destruction of the red giant will then follow when the envelope is ejected as the binary hardens. In a smaller number of encounters, the intruding binary passes through the star and ejects the red giant core from the envelope. The red giant envelope will then disperse on short time-scales. We compute the time-scales for both of these processes to occur in the Galactic Centre for a variety of binary populations.     

Spectroscopic observations of the eclipsing polar MN Hya (RX J0929–24)

We present low–medium resolution optical spectroscopy of the eclipsing AM Her system MN Hya (RX J0929–24). We determine the magnetic field strength at the primary accretion region of the white dwarf to be 42 MG from the spacing of cyclotron features visible during π ∼ 0.4–0.7. From spectra taken during the eclipse we find that the secondary has an M3–4 spectral type. Combined with the eclipse photometry of Sekiguchi, Nakada &38; Bassett and an estimate of the interstellar extinction we find a distance of ∼300–700 pc. We find unusual line variations at π ∼ 0.9: Hα is seen in absorption and emission. This is at the same point in the orbital phase at which a prominent absorption dip is seen in soft X-rays.     

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.     

Numerical simulations of protostellar encounters — II. Coplanar disc–disc encounters

It is expected that an average protostar will undergo at least one impulsive interaction with a neighbouring protostar whilst a large fraction of its mass is still in a massive, extended disc. Such interactions must have a significant impact upon the evolution of the protostars and their discs.   We have carried out a series of simulations of coplanar encounters between two stars, each possessing a massive circumstellar disc, using an SPH code that models gravitational, hydrodynamic and viscous forces. We find that during a coplanar encounter, disc material is swept up into a shock layer between the two interacting stars, and the layer then fragments to produce new protostellar condensations. The truncated remains of the discs may subsequently fragment; and the outer regions of the discs may be thrown off to form circumbinary disc-like structures around the stars. Thus coplanar disc–disc encounters lead efficiently to the formation of multiple star systems and small- N clusters, including substellar objects.     

A timing formula for main-sequence star binary pulsars

In binary radio pulsars with a main-sequence star companion, the spin-induced quadrupole moment of the companion gives rise to a precession of the binary orbit. As a first approximation one can model the secular evolution caused by this classical spin-orbit coupling by linear-in-time changes of the longitude of periastron and the projected semi-major axis of the pulsar orbit. This simple representation of the precession of the orbit neglects two important aspects of the orbital dynamics of a binary pulsar with an oblate companion. First, the quasiperiodic effects along the orbit, owing to the anisotropic 1/ r ³ nature of the quadrupole potential. Secondly, the long-term secular evolution of the binary orbit, which leads to an evolution of the longitude of periastron and the projected semi-major axis, which is non-linear in time.   In this paper a simple timing formula for binary radio pulsars with a main-sequence star companion is presented which models the short-term secular and most of the short-term periodic effects caused by the classical spin-orbit coupling. I also give extensions of the timing formula that account for long-term secular changes in the binary pulsar motion. It is shown that the short-term periodic effects are important for the timing observations of the binary pulsar PSR B1259–63. The long-term secular effects are likely to become important in the next few years of timing observations of the binary pulsar PSR J0045–7319. They could help to restrict or even determine the moments of inertia of the companion star and thus probe its internal structure.   Finally, I reinvestigate the spin-orbit precession of the binary pulsar PSR J0045–7319 since the analysis given in the literature is based on an incorrect expression for the precession of the longitude of periastron. A lower limit of 20° for the inclination of the B star with respect to the orbital plane is derived.     

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.     

Gamma-ray bursts from stellar remnants: probing the Universe at high redshift

A gamma-ray burst (GRB) releases an amount of energy similar to that of a supernova explosion, which combined with its rapid variability suggests an origin related to neutron stars or black holes. Since these compact stellar remnants form from the most massive stars not long after their birth, GRBs should trace the star formation rate in the Universe; we show that the GRB flux distribution is consistent with this. Because of the strong evolution of the star formation rate with redshift, it follows that the dimmest known bursts have z  ∼ 6, much above the value usually quoted and beyond the most distant quasars. This explains the absence of bright galaxies in well-studied GRB error boxes. The increased distances imply a peak luminosity of 8.3 × 10⁵¹ erg s⁻¹ and a rate density of 0.025 per million years per galaxy. These values are 20 times higher and 150 times lower, respectively, than are implied by fits with non-evolving GRB rates. This means either that GRBs are caused by a much rarer phenomenon than mergers of binary neutron stars, or that their gamma-ray emission is often invisible to us due to beaming. Precise burst locations from optical transients will discriminate between the various models for GRBs from stellar deaths, because the distance between progenitor birth place and burst varies greatly among them. The dimmest GRBs are then the most distant known objects, and may probe the Universe at an age when the first stars were forming.