The Evolution of Blue Stragglers Formed Via Stellar Collisions

We have used the results of recent smoothed particle hydrodynamic simulations of colliding stars to create models appropriate for input into a stellar evolution code. In evolving these models, we find that little or no surface convection occurs, precluding angular momentum loss via a magnetically-driven stellar wind as a viable mechanism for slowing rapidly rotating blue stragglers which have been formed by collisions. Angular momentum transfer to either a circumstellar disk (possibly collisional ejecta) or a nearby companion are plausible mechanisms for explaining the observed low rotation velocities of blue stragglers. Under the assumption that the blue stragglers seen in NGC 6397 and 47 Tuc have been created solely by collisions, we find that the majority of these blue stragglers cannot have been highly mixed by convection or meridional circulation currents at anytime during their evolution. Also, on the basis of the agreement between the predictions of our non-rotating models and the observed blue straggler distribution, the evolution of blue stragglers is apparently not dominated by the effects of rotation.Comment: 36 pages, including 1 table and 7 postscript figures (LaTeX2e). Also avaliable at http://astrowww.phys.uvic.ca/~ouellet/ . Accepted for publication in A

Merger of binary neutron stars with realistic equations of state in full general relativity

We present numerical results of three-dimensional simulations for the merger of binary neutron stars (BNSs) in full general relativity. Hybrid equations of state (EOSs) are adopted to mimic realistic nuclear EOSs. In this approach, we divide the EOSs into two parts, i.e., the thermal part and the cold part. For the cold part, we assign a fitting formula for realistic EOSs of cold nuclear matter slightly modifying the formula developed by Haensel and Potekhin. We adopt the SLy and FPS EOSs for which the maximum allowed ADM mass of cold and spherical neutron stars (NSs) is ~ 2.04Mo and 1.80Mo, respectively. Simulations are performed for BNSs of the total ADM mass in the range between 2.4Mo and 2.8Mo with the rest-mass ratio Q_M to be in the range 0.9 < Q_M < 1. It is found that if the total ADM mass of the system is larger than a threshold M_{thr}, a black hole (BH) is promptly formed in the merger irrespective of the mass ratios. In the other case, the outcome is a hypermassive NS of a large ellipticity, which results from the large adiabatic index of the realistic EOSs adopted. The value of M_{thr} depends on the EOS: ~ 2.7Mo and ~ 2.5Mo for the SLy and FPS EOSs, respectively. Gravitational waves are computed in terms of a gauge-invariant wave extraction technique. In the formation of the hypermassive NS, quasiperiodic gravitational waves of a large amplitude and of frequency between 3 and 4 kHz are emitted. The estimated emission time scale is < 100 ms, after which the hypermassive NS collapses to a BH. Because of the long emission time, the effective amplitude may be large enough to be detected by advanced laser interferometric gravitational wave detectors if the distance to the source is smaller than ~ 100 Mpc.Comment: Typos corrected, 2 references and comments on them added, 26 pages, 54 Postscript figures, Phys.Rev.D in pres

Gravitational Radiation from Nonaxisymmetric Instability in a Rotating Star

We present the first calculations of the gravitational radiation produced by nonaxisymmetric dynamical instability in a rapidly rotating compact star. The star deforms into a bar shape, shedding $\sim 4\%$ of its mass and $\sim 17\%$ of its angular momentum. The gravitational radiation is calculated in the quadrupole approximation. For a mass $M \sim 1.4$ M$_{\odot}$ and radius $R \sim 10$ km, the gravitational waves have frequency $\sim 4$ kHz and amplitude $h \sim 2 \times 10^{-22}$ at the distance of the Virgo Cluster. They carry off energy $\Delta E/M \sim 0.1\%$ and radiate angular momentum $\Delta J/J \sim 0.7\%$.Comment: 16 pages, LaTeX with REVTEX macros, reprints available - send mailing address to [email protected]. Published: PRL 72, 1314 (1994

The two hybrid B-type pulsators: Nu Eridani and 12 Lacertae

The rich oscillation spectra determined for the two stars, Nu Eridani and 12 Lacertae, present an interesting challenge to stellar modelling. The stars are hybrid objects showing a number of modes at frequencies typical for Beta Cep stars but also one mode at frequency typical for SPB stars. We construct seismic models of these stars considering uncertainties in opacity and element distribution. We also present estimate of the interior rotation rate and address the matter of mode excitation. We use both the OP and OPAL opacity data and find significant difference in the results. Uncertainty in these data remains a major obstacle in precise modelling of the objects and, in particular, in estimating the overshooting distance. We find evidence for significant rotation rate increase between envelope and core in the two stars. Instability of low-frequency g-modes was found in seismic models of Nu Eri built with the OP data, but at frequencies higher than those measured in the star. No such instability was found in models of 12 Lac. We do not have yet a satisfactory explanation for low frequency modes. Some enhancement of opacity in the driving zone is required but we argue that it cannot be achieved by the iron accumulation, as it has been proposed.Comment: 9 pages, 8 figures, accepted for publication in MNRA

Dynamical Tide in Solar-Type Binaries

Circularization of late-type main-sequence binaries is usually attributed to turbulent convection, while that of early-type binaries is explained by resonant excitation of g modes. We show that the latter mechanism operates in solar-type stars also and is at least as effective as convection, despite inefficient damping of g modes in the radiative core. The maximum period at which this mechanism can circularize a binary composed of solar-type stars in 10 Gyr is as low as 3 days, if the modes are damped by radiative diffusion only and g-mode resonances are fixed; or as high as 6 days, if one allows for evolution of the resonances and for nonlinear damping near inner turning points. Even the larger theoretical period falls short of the observed transition period by a factor two.Comment: 17 pages, 2 postscript figures, uses aaspp4.sty. Submitted to Ap

Gravitational Waves from the Dynamical Bar Instability in a Rapidly Rotating Star

A rapidly rotating, axisymmetric star can be dynamically unstable to an m=2 "bar" mode that transforms the star from a disk shape to an elongated bar. The fate of such a bar-shaped star is uncertain. Some previous numerical studies indicate that the bar is short lived, lasting for only a few bar-rotation periods, while other studies suggest that the bar is relatively long lived. This paper contains the results of a numerical simulation of a rapidly rotating gamma=5/3 fluid star. The simulation shows that the bar shape is long lived: once the bar is established, the star retains this shape for more than 10 bar-rotation periods, through the end of the simulation. The results are consistent with the conjecture that a star will retain its bar shape indefinitely on a dynamical time scale, as long as its rotation rate exceeds the threshold for secular bar instability. The results are described in terms of a low density neutron star, but can be scaled to represent, for example, a burned-out stellar core that is prevented from complete collapse by centrifugal forces. Estimates for the gravitational-wave signal indicate that a dynamically unstable neutron star in our galaxy can be detected easily by the first generation of ground based gravitational-wave detectors. The signal for an unstable neutron star in the Virgo cluster might be seen by the planned advanced detectors. The Newtonian/quadrupole approximation is used throughout this work.Comment: Expanded version to be published in Phys. Rev. D: 13 pages, REVTeX, 13 figures, 9 TeX input file

Faraday resonance in dynamical bar instability of differentially rotating stars

We investigate the nonlinear behaviour of the dynamically unstable rotating star for the bar mode by three-dimensional hydrodynamics in Newtonian gravity. We find that an oscillation along the rotation axis is induced throughout the growth of the unstable bar mode, and that its characteristic frequency is twice as that of the bar mode, which oscillates mainly along the equatorial plane. A possibility to observe Faraday resonance in gravitational waves is demonstrated and discussed.Comment: 13 pages with 13 figures, revtex4.cls. Accepted for publication in the Physical Review

(In)finite extent of stationary perfect fluids in Newtonian theory

For stationary, barotropic fluids in Newtonian gravity we give simple criteria on the equation of state and the "law of motion" which guarantee finite or infinite extent of the fluid region (providing a priori estimates for the corresponding stationary Newton-Euler system). Under more restrictive conditions, we can also exclude the presence of "hollow" configurations. Our main result, which does not assume axial symmetry, uses the virial theorem as the key ingredient and generalises a known result in the static case. In the axially symmetric case stronger results are obtained and examples are discussed.Comment: Corrections according to the version accepted by Ann. Henri Poincar

A Robust Measure of Tidal Circularization in Coeval Binary Populations: The solar-type spectroscopic Binary Population in The Open Cluster M35

We present a new homogeneous sample of 32 spectroscopic binary orbits in the young (~ 150 Myr) main-sequence open cluster M35. The distribution of orbital eccentricity vs. orbital period (e-log(P)) displays a distinct transition from eccentric to circular orbits at an orbital period of ~ 10 days. The transition is due to tidal circularization of the closest binaries. The population of binary orbits in M35 provide a significantly improved constraint on the rate of tidal circularization at an age of 150 Myr. We propose a new and more robust diagnostic of the degree of tidal circularization in a binary population based on a functional fit to the e-log(P) distribution. We call this new measure the tidal circularization period. The tidal circularization period of a binary population represents the orbital period at which a binary orbit with the most frequent initial orbital eccentricity circularizes (defined as e = 0.01) at the age of the population. We determine the tidal circularizationperiod for M35 as well as for 7 additional binary populations spanning ages from the pre main-sequence (~ 3 Myr) to late main-sequence (~ 10 Gyr), and use Monte Carlo error analysis to determine the uncertainties on the derived circularization periods. We conclude that current theories of tidal circularization cannot account for the distribution of tidal circularization periods with population age.Comment: 37 pages, 9 figures, to be published in The Astrophysical Journal, February 200

Accurate evolutions of unequal-mass neutron-star binaries: properties of the torus and short GRB engines

We present new results from accurate and fully general-relativistic simulations of the coalescence of unmagnetized binary neutron stars with various mass ratios. The evolution of the stars is followed through the inspiral phase, the merger and prompt collapse to a black hole, up until the appearance of a thick accretion disk, which is studied as it enters and remains in a regime of quasi-steady accretion. Although a simple ideal-fluid equation of state with \Gamma=2 is used, this work presents a systematic study within a fully general relativistic framework of the properties of the resulting black-hole--torus system produced by the merger of unequal-mass binaries. More specifically, we show that: (1) The mass of the torus increases considerably with the mass asymmetry and equal-mass binaries do not produce significant tori if they have a total baryonic mass M_tot >~ 3.7 M_sun; (2) Tori with masses M_tor ~ 0.2 M_sun are measured for binaries with M_tot ~ 3.4 M_sun and mass ratios q ~ 0.75-0.85; (3) The mass of the torus can be estimated by the simple expression M_tor(q, M_tot) = [c_1 (1-q) + c_2](M_max-M_tot), involving the maximum mass for the binaries and coefficients constrained from the simulations, and suggesting that the tori can have masses as large as M_tor ~ 0.35 M_sun for M_tot ~ 2.8 M_sun and q ~ 0.75-0.85; (4) Using a novel technique to analyze the evolution of the tori we find no evidence for the onset of non-axisymmetric instabilities and that very little, if any, of their mass is unbound; (5) Finally, for all the binaries considered we compute the complete gravitational waveforms and the recoils imparted to the black holes, discussing the prospects of detection of these sources for a number of present and future detectors.Comment: 35 pages; small changes to match the published versio