Forming short-period Wolf-Rayet X-ray binaries and double black holes through stable mass transfer

We show that black-hole High-Mass X-ray Binaries (HMXBs) with O- or B-type donor stars and relatively short orbital periods, of order one week to several months may survive spiral in, to then form Wolf-Rayet (WR) X-ray binaries with orbital periods of order a day to a few days; while in systems where the compact star is a neutron star, HMXBs with these orbital periods never survive spiral-in. We therefore predict that WR X-ray binaries can only harbor black holes. The reason why black-hole HMXBs with these orbital periods may survive spiral in is: the combination of a radiative envelope of the donor star, and a high mass of the compact star. In this case, when the donor begins to overflow its Roche lobe, the systems are able to spiral in slowly with stable Roche-lobe overflow, as is shown by the system SS433. In this case the transferred mass is ejected from the vicinity of the compact star (so-called "isotropic re-emission" mass loss mode, or "SS433-like mass loss"), leading to gradual spiral-in. If the mass ratio of donor and black hole is $>3.5$, these systems will go into CE evolution and are less likely to survive. If they survive, they produce WR X-ray binaries with orbital periods of a few hours to one day. Several of the well-known WR+O binaries in our Galaxy and the Magellanic Clouds, with orbital periods in the range between a week and several months, are expected to evolve into close WR-Black-Hole binaries,which may later produce close double black holes. The galactic formation rate of double black holes resulting from such systems is still uncertain, as it depends on several poorly known factors in this evolutionary picture. It might possibly be as high as $\sim 10^{-5}$ per year.Comment: MNRAS in pres

Mutation supply and the repeatability of selection for antibiotic resistance

Whether evolution can be predicted is a key question in evolutionary biology. Here we set out to better understand the repeatability of evolution. We explored experimentally the effect of mutation supply and the strength of selective pressure on the repeatability of selection from standing genetic variation. Different sizes of mutant libraries of an antibiotic resistance gene, TEM-1 $\beta$-lactamase in Escherichia coli, were subjected to different antibiotic concentrations. We determined whether populations went extinct or survived, and sequenced the TEM gene of the surviving populations. The distribution of mutations per allele in our mutant libraries- generated by error-prone PCR- followed a Poisson distribution. Extinction patterns could be explained by a simple stochastic model that assumed the sampling of beneficial mutations was key for survival. In most surviving populations, alleles containing at least one known large-effect beneficial mutation were present. These genotype data also support a model which only invokes sampling effects to describe the occurrence of alleles containing large-effect driver mutations. Hence, evolution is largely predictable given cursory knowledge of mutational fitness effects, the mutation rate and population size. There were no clear trends in the repeatability of selected mutants when we considered all mutations present. However, when only known large-effect mutations were considered, the outcome of selection is less repeatable for large libraries, in contrast to expectations. Furthermore, we show experimentally that alleles carrying multiple mutations selected from large libraries confer higher resistance levels relative to alleles with only a known large-effect mutation, suggesting that the scarcity of high-resistance alleles carrying multiple mutations may contribute to the decrease in repeatability at large library sizes.Comment: 31pages, 9 figure

Expected Coalescence Rate of Double Neutron Stars for Ground Based Interferometers

In this paper we present new estimates of the coalescence rate of neutron star binaries in the local universe and we discuss its consequences for the first generations of ground based interferometers. Our approach based on both evolutionary and statistical methods gives a galactic merging rate of 1.7 10$^{-5}$ yr$^{-1}$, in the range of previous estimates 10$^{-6}$ - 10$^{-4}$ yr$^{-1}$. The local rate which includes the contribution of elliptical galaxies is two times higher, in the order of 3.4 10$^{-5}$ yr$^{-1}$. We predict one detection every 148 and 125 years with initial VIRGO and LIGO, and up to 6 events per year with their advanced configuration. Our recent detection rate estimates from investigations on VIRGO future improvements are quoted.Comment: talk given at the GWDAW9 (Annecy, 2004) to be published in CQ

The Evolution of Globular Clusters in the Galaxy

We investigate the evolution of globular clusters using N-body calculations and anisotropic Fokker-Planck (FP) calculations. The models include a mass spectrum, mass loss due to stellar evolution, and the tidal field of the parent galaxy. Recent N-body calculations have revealed a serious discrepancy between the results of N-body calculations and isotropic FP calculations. The main reason for the discrepancy is an oversimplified treatment of the tidal field employed in the isotropic FP models. In this paper we perform a series of calculations with anisotropic FP models with a better treatment of the tidal boundary and compare these with N-body calculations. The new tidal boundary condition in our FP model includes one free parameter. We find that a single value of this parameter gives satisfactory agreement between the N-body and FP models over a wide range of initial conditions. Using the improved FP model, we carry out an extensive survey of the evolution of globular clusters over a wide range of initial conditions varying the slope of the mass function, the central concentration, and the relaxation time. The evolution of clusters is followed up to the moment of core collapse or the disruption of the clusters in the tidal field of the parent galaxy. In general, our model clusters, calculated with the anisotropic FP model with the improved treatment for the tidal boundary, live longer than isotropic models. The difference in the lifetime between the isotropic and anisotropic models is particularly large when the effect of mass loss via stellar evolution is rather significant. On the other hand the difference is small for relaxation- dominated clusters which initially have steep mass functions and high central concentrations.Comment: 36 pages, 11 figures, LaTeX; added figures and tables; accepted by Ap

A platform independent communication library for distributed computing

We present MPWide, a platform independent communication library for performing message passing between supercomputers. Our library couples several local MPI applications through a long distance network using, for example, optical links. The implementation is deliberately kept light-weight, platform independent and the library can be installed and used without administrative privileges. The only requirements are a C++ compiler and at least one open port to a wide area network on each site. In this paper we present the library, describe the user interface, present performance tests and apply MPWide in a large scale cosmological N-body simulation on a network of two computers, one in Amsterdam and the other in Tokyo

On the Theory of Gamma Ray Bursts and Hypernovae: The Black Hole Soft X-ray Transient Sources

We show that a common evolutionary history can produce the black hole binaries in the Galaxy in which the black holes have masses of ~ 5-10 M_sun. In with low-mass, <~ 2.5 M_sun, ZAMS (zero age main sequence) companions, the latter remain in main sequence during the active stage of soft X-ray transients (SXTs), most of them being of K or M classification. In two intermediate cases, IL Lupi and Nova Scorpii with ZAMS ~ 2.5 M_sun companions the orbits are greatly widened because of large mass loss in the explosion forming the black hole, and whereas these companions are in late main sequence evolution, they are close to evolving. Binaries with companion ZAMS masses >~ 3 M_sun are initially "silent" until the companion begins evolving across the Herzsprung gap. We provide evidence that the narrower, shorter period binaries, with companions now in main sequence, are fossil remnants of gamma ray bursters (GRBs). We also show that the GRB is generally accompanied by a hypernova explosion (a very energetic supernova explosion). We further show that the binaries with evolved companions are good models for some of the ultraluminous X-ray sources (ULXs) recently seen by Chandra in other galaxies. The great regularity in our evolutionary history, especially the fact that most of the companions of ZAMS mass <~ 2.5 M_sun remain in main sequences as K or M stars can be explained by the mass loss in common envelope evolution to be Case C; i.g., to occur only after core He burning has finished. Since our argument for Case C mass transfer is not generally understood in the community, we add an appendix, showing that with certain assumptions which we outline we can reproduce the regularities in the evolution of black hole binaries by Case C mass transfer.Comment: 59 pages, 12 figures, review articl

MRI for the evaluation of gastric physiology

Radiolog