Eccentricities of Double Neutron Star Binaries

Recent pulsar surveys have increased the number of observed double neutron stars (DNS) in our galaxy enough so that observable trends in their properties are starting to emerge. In particular, it has been noted that the majority of DNS have eccentricities less than 0.3, which are surprisingly low for binaries that survive a supernova explosion that we believe imparts a significant kick to the neutron star. To investigate this trend, we generate many different theoretical distributions of DNS eccentricities using Monte Carlo population synthesis methods. We determine which eccentricity distributions are most consistent with the observed sample of DNS binaries. In agreement with Chaurasia & Bailes (2005), assuming all double neutron stars are equally as probable to be discovered as binary pulsars, we find that highly eccentric, coalescing DNS are less likely to be observed because of their accelerated orbital evolution due to gravitational wave emission and possible early mergers. Based on our results for coalescing DNS, we also find that models with vanishingly or moderately small kicks (sigma < about 50 km/s) are inconsistent with the current observed sample of such DNS. We discuss the implications of our conclusions for DNS merger rate estimates of interest to ground-based gravitational-wave interferometers. We find that, although orbital evolution due to gravitational radiation affects the eccentricity distribution of the observed sample, the associated upwards correction factor to merger rate estimates is rather small (typically 10-40%).Comment: 9 pages, 8 figures, accepted by ApJ. Figures reduced and some content changed, references adde

The Galactic Population of Low- and Intermediate-Mass X-ray Binaries

(abridged) We present the first study that combines binary population synthesis in the Galactic disk and detailed evolutionary calculations of low- and intermediate-mass X-ray binaries (L/IMXBs). We show that the formation probability of IMXBs with initial donor masses of 1.5--4 Msun is typically >~5 times higher than that of standard LMXBs, and suggest that the majority of the observed systems may have descended from IMXBs. Distributions at the current epoch of the orbital periods, donor masses, and mass accretion rates have been computed, as have orbital-period distributions of BMPs. Several significant discrepancies between the theoretical and observed distributions are discussed. The orbital-period distribution of observed BMPs strongly favors cases where the envelope of the neutron-star progenitor is more easily ejected during the common-envelope phase. However, this leads to a >~100-fold overproduction of the theoretical number of luminous X-ray sources relative to the total observed number of LMXBs. X-ray irradiation of the donor star may result in a dramatic reduction in the X-ray active lifetime of L/IMXBs, thus possibly resolving the overproduction problem, as well as the long-standing BMP/LMXB birthrate problem.Comment: 12 pages, emulateapj, submitted to Ap

A New Way to Detect Massive Black Holes in Galaxies: The Stellar Remnants of Tidal Disruption

We point out that the tidal disruption of a giant may leave a luminous (10^35-10^39 ergs/s), hot (10-100 eV) stellar core. The ``supersoft'' source detected by Chandra at the center of M31 may be such a core; whether or not it is, the observations have shown that such a core is detectable, even in the center of a galaxy. We therefore explore the range of expected observational signatures and how they may be used to (1) test the hypothesis that the M31 source is a remnant of tidal stripping and (2) discover evidence of black holes and disruption events in other galaxies.Comment: Four pages with 1 figure. Appeared in ApJL (2001, 551, L37

Modelling the water budget and the riverflows of the Maritsa basin in Bulgaria

International audienceA soil-vegetation-atmosphere transfer model coupled with a macroscale distributed hydrological model was used in order to simulate the water cycle for a large region in Bulgaria. To do so, an atmospheric forcing was built for two hydrological years (1 October 1995 to 30 September 1997), at an eight km resolution. It was based on the data available at the National Institute of Meteorology and Hydrology (NIMH) of Bulgaria. Atmospheric parameters were carefully checked and interpolated with a high level of detail in space and time (3-h step). Comparing computed Penman evapotranspiration versus observed pan evaporation validated the quality of the implemented forcing. The impact of the human activities on the rivers (especially hydropower or irrigation) was taken into account. Some improvements of the hydrometeorological model were made: for better simulation of summer riverflow, two additional reservoirs were added to simulate the slow component of the runoff. Those reservoirs were calibrated using the observed data of the 1st year, while the 2nd year was used for validation. 56 hydrologic stations and 12 dams were used for the model calibration while 41 rivergages were used for the validation of the model. The results compare well with the daily-observed discharges, with good results obtained over more than 25% of the rivergages. The simulated snow depth was compared to daily measurements at 174 stations and the evolution of the snow water equivalent was validated at 5 sites. The process of melting and refreezing of snow was found to be important on this region. The comparison of the normalized values of simulated versus measured soil moisture showed good correlation. The surface water budget shows large spatial variations due to the elevation influence on the precipitations, soil properties and vegetation variability. An inter annual difference was observed in the water cycle as the first year was more influenced by Mediterranean climate, while the second year was characterised by continental influence. Energy budget shows a dominating sensible heat component in summer, due to the fact that the water stress limits the evaporation. This study is a first step for the implementation of an operational hydrometeorological model that could be used for real time monitoring and forecast the water budget and the riverflow of Bulgaria

Equipotential Surfaces and Lagrangian points in Non-synchronous, Eccentric Binary and Planetary Systems

We investigate the existence and properties of equipotential surfaces and Lagrangian points in non-synchronous, eccentric binary star and planetary systems under the assumption of quasi-static equilibrium. We adopt a binary potential that accounts for non-synchronous rotation and eccentric orbits, and calculate the positions of the Lagrangian points as functions of the mass ratio, the degree of asynchronism, the orbital eccentricity, and the position of the stars or planets in their relative orbit. We find that the geometry of the equipotential surfaces may facilitate non-conservative mass transfer in non-synchronous, eccentric binary star and planetary systems, especially if the component stars or planets are rotating super-synchronously at the periastron of their relative orbit. We also calculate the volume-equivalent radius of the Roche lobe as a function of the four parameters mentioned above. Contrary to common practice, we find that replacing the radius of a circular orbit in the fitting formula of Eggleton (1983) with the instantaneous distance between the components of eccentric binary or planetary systems does not always lead to a good approximation to the volume-equivalent radius of the Roche-lobe. We therefore provide generalized analytic fitting formulae for the volume-equivalent Roche lobe radius appropriate for non-synchronous, eccentric binary star and planetary systems. These formulae are accurate to better than 1% throughout the relevant 2-dimensional parameter space that covers a dynamic range of 16 and 6 orders of magnitude in the two dimensions.Comment: 12 pages, 10 figures, 2 Tables, Accepted by the Astrophysical Journa

The Search for Stellar Companions to Exoplanet Host Stars Using the CHARA Array

Most exoplanets have been discovered via radial velocity studies, which are inherently insensitive to orbital inclination. Interferometric observations will show evidence of a stellar companion if it sufficiently bright, regardless of the inclination. Using the CHARA Array, we observed 22 exoplanet host stars to search for stellar companions in low-inclination orbits that may be masquerading as planetary systems. While no definitive stellar companions were discovered, it was possible to rule out certain secondary spectral types for each exoplanet system observed by studying the errors in the diameter fit to calibrated visibilities and by searching for separated fringe packets.Comment: 26 pages, 5 tables, 8 figure

The Role of Helium Stars in the Formation of Double Neutron Stars

We have calculated the evolution of 60 model binary systems consisting of helium stars in the mass range of M_He= 2.5-6Msun with a 1.4Msun neutron star companion to investigate the formation of double neutron star systems.Orbital periods ranging from 0.09 to 2 days are considered, corresponding to Roche lobe overflow starting from the helium main sequence to after the ignition of carbon burning in the core. We have also examined the evolution into a common envelope phase via secular instability, delayed dynamical instability, and the consequence of matter filling the neutron star's Roche lobe. The survival of some close He-star neutron-star binaries through the last mass transfer episode (either dynamically stable or unstable mass transfer phase) leads to the formation of extremely short-period double neutron star systems (with P<~0.1days). In addition, we find that systems throughout the entire calculated mass range can evolve into a common envelope phase, depending on the orbital period at the onset of mass transfer. The critical orbital period below which common envelope evolution occurs generally increases with M_He. In addition, a common envelope phase may occur during a short time for systems characterized by orbital periods of 0.1-0.5 days at low He-star masses (~ 2.6-3.3Msun). The existence of a short-period population of double neutron stars increases the predicted detection rate of inspiral events by ground-based gravitational-wave detectors and impacts their merger location in host galaxies and their possible role as gamma-ray burst progenitors. We use a set of population synthesis calculations and investigate the implications of the mass-transfer results for the orbital properties of DNS populations.Comment: 30 pages, Latex (AASTeX), 1 table, 8 figures. To appear in ApJ, v592 n1 July 20, 200

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, and (iii) a single star and a hard binary intermediate-mass black hole. 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^6 -10^7 stars pc^{-3}. Although such high densities may exist during the core collapse of young massive star clusters, we caution that they have never been observed.Comment: 11 pages, 6 figures, 1 table, accepted to MNRA

Evolution of Neutron-Star, Carbon-Oxygen White-Dwarf Binaries

At least one, but more likely two or more, eccentric neutron-star, carbon-oxygen white-dwarf binaries with an unrecycled pulsar have been observed. According to the standard scenario for evolving neutron stars which are recycled in common envelope evolution we expect to observe \gsim 50 such circular neutron star-carbon oxygen white dwarf binaries, since their formation rate is roughly equal to that of the eccentric binaries and the time over which they can be observed is two orders of magnitude longer, as we shall outline. We observe at most one or two such circular binaries and from that we conclude that the standard scenario must be revised. Introducing hypercritical accretion into common envelope evolution removes the discrepancy by converting the neutron star into a black hole which does not emit radio waves, and therefore would not be observed.Comment: 25 pages, 1 figure, accepted in Ap

The late stages of evolution of helium star-neutron star binaries and the formation of double neutron star systems

With a view to understanding the formation of double neutron-stars (DNS), we investigate the late stages of evolution of helium stars with masses of 2.8 - 6.4 Msun in binary systems with a 1.4 Msun neutron-star companion. We found that mass transfer from 2.8 - 3.3 Msun helium stars and from 3.3 - 3.8 Msun in very close orbits (P_orb > 0.25d) will end up in a common-envelope (CE) and spiral-in phase due to the development of a convective helium envelope. If the neutron star has sufficient time to complete the spiraling-in process before the core collapses, the system will produce very tight DNSs (P_orb ~ 0.01d) with a merger timescale of the order of 1 Myr or less. These systems would have important consequences for the detection rate of GWR and for the understanding of GRB progenitors. On the other hand, if the time left until the explosion is shorter than the orbital-decay timescale, the system will undergo a SN explosion during the CE phase. Helium stars with masses 3.3 - 3.8 Msun in wider orbits (P_orb > 0.25d) and those more massive than 3.8 Msun do not go through CE evolution. The remnants of these massive helium stars are DNSs with periods in the range of 0.1 - 1 d. This suggests that this range of mass includes the progenitors of the galactic DNSs with close orbits (B1913+16 and B1534+12). A minimum kick velocity of 70 km/s and 0 km/s (for B1913+16 and B1534+12, respectively) must have been imparted at the birth of the pulsar's companion. The DNSs with wider orbits (J1518+4904 and probably J1811-1736) are produced from helium star-neutron star binaries which avoid RLOF, with the helium star more massive than 2.5 Msun. For these systems the minimum kick velocities are 50 km/s and 10 km/s (for J1518+4904 and J1811-1736, respectively).Comment: 16 pages, latex, 12 figures, accepted for publication in MNRA