Melvyn B. Davies

ULTRALUMINOUS X-RAY SOURCES IN EXTERNAL GALAXIES

We investigate models for the class of ultraluminous nonnuclear X-ray sources (i.e., ultraluminous compact X-ray sources [ULXs]) seen in a number of galaxies and probably associated with star-forming regions. Models in which the X-ray emission is assumed to be isotropic run into several difficulties. In particular, the formation of sufficient numbers of the required ultramassive black hole X-ray binaries is problematic, and the likely transient behavior of the resulting systems is not in good accord with observation. The assumption of mild X-ray beaming suggests instead that ULXs may represent a short-lived but extremely common stage in the evolution of a wide class of X-ray binaries. The best candidate for this is the phase of thermal-timescale mass transfer that is inevitable in many intermediate- and high-mass X-ray binaries. This in turn suggests a link with the Galactic microquasars. The short lifetimes of high-mass X-ray binaries would explain the association of ULXs with episodes of star formation. These considerations still allow the possibility that individual ULXs may contain extremely massive black holes.

High-resolution calculations of merging neutron stars - III. Gamma-ray bursts

Recent three-dimensional, high-resolution simulations of neutron star coalescences are analysed to assess whether short gamma-ray bursts (GRBs) could originate from such encounters. The two most popular modes of energy extraction - namely the annihilation of and magnetohydrodynamic processes - are explored in order to investigate their viability in launching GRBs. We find that annihilation can provide the necessary stresses to drive a highly relativistic expansion. However, unless the outflow is beamed into less than 1 per cent of the solid angle, this mechanism may fail to explain the apparent isotropized energies implied for short GRBs at cosmological distances. We argue that the energetic, neutrino-driven wind that accompanies the merger event will have enough pressure to provide adequate collimation to the -annihilation-driven jet, thereby comfortably satisfying constraints on event rate and apparent luminosity. We also assess magnetic mechanisms to transform the available energy into a GRB. If the central object does not collapse immediately into a black hole, it will be convective and it is expected to act as an effective large scale dynamo, amplifying the seed magnetic fields to a few times 1017 G within a small fraction of a second. The associated spindown time-scale is 0.2 s, coinciding with the typical duration of a short GRB. The efficiencies of the various assessed magnetic processes are high enough to produce isotropized luminosities in excess of 1052 erg s-1 even without beaming. (Less)

Merging white dwarf / black hole binaries and gamma-ray bursts

The merger of compact binaries, especially black holes and neutron stars, is frequently invoked to explain gamma-ray bursts (GRBs). In this paper, we present three-dimensional hydrodynamical simulations of the relatively neglected mergers of white dwarfs and black holes. During the merger, the white dwarf is tidally disrupted and sheared into an accretion disk. Nuclear reactions are followed, and the energy release is negligible. Peak accretion rates are ~0.05 M☉ s-1 (less for lower mass white dwarfs) and last for approximately a minute. Many of the disk parameters can be explained by a simple analytic model that we derive and compare to our simulations. This model can be used to predict accretion rates for white dwarf and black hole (or neutron star) masses that are not simulated here. Although the mergers studied here create disks with larger radii and longer accretion times than those from the merger of double neutron stars, a larger fraction of the white dwarfs mass becomes part of the disk. Thus the merger of a white dwarf and a black hole could produce a long-duration GRB. The event rate of these mergers may be as high as 10-6 yr-1 per galaxy.

Blue straggler production in globular clusters

Recent Hubble Space Telescope observations of a large sample of globular clusters reveal that every cluster contains between 40 and 400 blue stragglers.The population does not correlate with either stellar collision rate (as would be expected if all blue stragglers were formed via collisions) or total mass (as would be expected if all blue stragglers were formed via the unhindered evolution of a subset of the stellar population). In this paper, we support the idea that blue stragglers are made through both channels. The number produced via collisions tends to increase with cluster mass. In this paper we show how the current population produced from primordial binaries decreases with increasing cluster mass;exchange encounters with third, single stars in the most massivec lusters tend to reduce the fraction of binaries containing a primary close to the current turn-off mass. Rather, their primaries tend to be somewhat more massive (~1-3 Msolar) and have evolved off the main sequence, filling their Roche lobes in the past, often converting their secondaries into blue stragglers (but more than 1 Gyr or so ago and thus they are no longer visible as blue stragglers). We show that this decline in the primordial blue straggler population is likely to be offset by the increase in the number of blue stragglers produced via collisions. The predicted total blue straggler population is therefore relatively independent of cluster mass, thus matching the observed population. This result does not depend on any particular assumed blue straggler lifetime. (Less)

Relative Frequencies of Blue Stragglers in Galactic Globular Clusters: Constraints for the Formation Mechanisms

We discuss the main properties of the Galactic globular cluster (GC)blue straggler stars (BSSs), as inferred from our new catalog containingnearly 3000 BSSs. The catalog has been extracted from thephotometrically homogeneous V versus (B-V) color-magnitude diagrams(CMDs) of 56 GCs, based on Wide Field Planetary Camera 2 images of theircentral cores. In our analysis, we used consistent relative distancesbased on the same photometry and calibration. The number of BSSs hasbeen normalized to obtain relative frequencies (FBSS) andspecific densities (NS) using different stellar populationsextracted from the CMD. The cluster FBSS is significantlysmaller than the relative frequency of field BSSs. We find a significantanticorrelation between the BSS relative frequency in a cluster and itstotal absolute luminosity (mass). There is no statistically significanttrend between the BSS frequency and the expected collision rate. Thevalue of FBSS does not depend on other cluster parameters,apart from a mild dependence on the central density. Post-core-collapseclusters act like normal clusters as far as the BSS frequency isconcerned. We also show that the BSS luminosity function for the mostluminous clusters is significantly different, with a brighter peak andextending to brighter luminosities than in the less luminous clusters.These results imply that the efficiency of BSS production mechanisms andtheir relative importance vary with the cluster mass.Based on observations with the NASA/ESA Hubble Space Telescope, obtainedat the Space Telescope Science Institute, which is operated by theAssociation of Universities for Research in Astronomy, Inc., under NASAcontract NAS 5-26555. (Less)

Merging neutron stars. 1. Initial results for coalescence of noncorotating systems

We present three-dimensional Newtonian simulations of the coalescence of two neutron stars, using a smoothed particle hydrodynamics (SPH) code. We begin the simulations with the two stars in a hard, circular binary, and have them spiral together as angular momentum is lost through gravitational radiation at the rate predicted by modeling the system as two point masses. We model the neutron stars as hard polytropes (gamma = 2.4) of equal mass, and investigate the effect of the initial spin of the two stars on the coalescence. The process of coalescence, from initial contact to the formation of an axially symmetric object, takes only a few orbital periods. Some of the material from the two neutron stars is shed, forming a thick disk around the central, coalesced object. The mass of this disk is dependent on the initial neutron star spins; higher spin rates result in greater mass loss and thus more massive disks. For spin rates that are most likely to be applicable to real systems, the central coalesced object has a mass of 2.4 solar mass, which is tantalizingly close to the maximum mass allowed by any neutron star equation of state for an object that is supported in part by rotation. Using a realistic nuclear equation of state, we estimate the temperature of the material after the coalescence. We find that the central object is at a temperature of approximately 10 MeV, while the disk is heated by shocks to a temperature of 2 to 4 MeV.

Gamma-Ray Bursts, Supernova Kicks, and Gravitational Radiation

We suggest that the collapsing core of a massive rotating star may fragment to produce two or more compact objects. Their coalescence under gravitational radiation gives the resulting black hole or neutron star a significant kick velocity, which may explain those observed in pulsars. A gamma-ray burst can result only when this kick is small. Thus, only a small fraction of core-collapse supernovae produce gamma-ray bursts. The burst may be delayed significantly (hours to days) after the supernova, as suggested by recent observations. If our picture is correct, core-collapse supernovae should be significant sources of gravitational radiation with a chirp signal similar to a coalescing neutron star binary.

Post-Newtonian Smoothed Particle Hydrodynamics

We introduce an adaptation of the well-known tree+SPH numerical scheme to post-Newtonian (PN) hydrodynamics and gravity. Our code solves the (0 + 1 + 2.5)PN equations. These equations include Newtonian hydrodynamics and gravity (0PN), the first-order relativistic corrections to those (1PN), and the lowest order gravitational radiation terms (2.5PN). We test various aspects of our code using analytically solvable test problems. We then proceed to study the 1PN effects on binary neutron star coalescence by comparing calculations with and without the 1PN terms. We find that the effect of the 1PN terms is rather small. The largest effect arises with a stiff equation of state for which the maximum rest mass density increases. This could induce black hole formation. The gravitational wave luminosity is also affected.We introduce an adaptation of the well known Tree+SPH numerical scheme to Post Newtonian (PN) hydrodynamics and gravity. Our code solves the (0+1+2.5)PN equations. These equations include Newtonian hydrodynamics and gravity (0PN), the first order relativistic corrections to those (1PN) and the lowest order gravitational radiation terms (2.5PN). We test various aspects of our code using analytically solvable test problems. We then proceed to study the 1PN effects on binary neutron star coalescence by comparing calculations with and without the 1PN terms. We find that the effect of the 1PN terms is rather small. The largest effect arises with a stiff equation of state for which the maximum rest mass density increases. This could induce black hole formation. The gravitational wave luminosity is also affected. Subject headings: gravitation,hydrodynamic,relativity,stars: neutron

The minimum orbital period in thermal time-scale mass transfer

We show that the usual picture of supersoft X-ray binary evolution as driven by conservative thermal time-scale mass transfer cannot explain the short orbital periods of RX J0537.7–7034 (3.5 h) and 1E 0035.4–7230 (4.1 h). Non-conservative evolution may produce such periods, but requires very significant mass loss, and is highly constrained.

Formation of the binary pulsars J1141-6545 and B2303+46

The binaries PSR J1141-6545 and PSR B2303+46 each appear to contain a white dwarf that formed before the neutron star. We describe an evolutionary pathway to produce these two systems. In this scenario, the primary transfers its envelope on to the secondary, which is then the more massive of the two stars, and indeed sufficiently massive later to produce a neutron star via a supernova. The core of the primary produces a massive white dwarf, which enters into a common envelope with the core of the secondary when the latter evolves off the main sequence. During the common-envelope phase, the white dwarf and the core of the secondary spiral together as the envelope is ejected. The evolutionary histories of PSR J1141-6545 and PSR B2303+46 differ after this phase. In the case of PSR J1141-6545, the secondary (now a helium star) evolves into contact transferring its envelope on to the white dwarf. We propose that the vast majority of this material is, in fact, ejected from the system. The remains of the secondary then explode as a supernova, producing a neutron star. Generally the white dwarf and neutron star will remain bound in tight, often eccentric, systems resembling PSR J1141-6545. These systems will spiral in and merge on a relatively short time-scale and may make a significant contribution to the population of gamma-ray burst progenitors. In PSR B2303+46, the helium-star secondary and white dwarf never come into contact. Rather the helium star loses its envelope via a wind, which increases the binary separation slightly. Only a small fraction of such systems will remain bound when the neutron star is formed (as the systems are wider). Those systems that are broken up will produce a population of high-velocity white dwarfs and neutron stars.