Engineering topics

Rasio, Frederic A.

Publications and source records attributed to Rasio, Frederic A..

At least 19 records

NASA NTRS1995

On blue straggler information by direct collisions of main sequence stars

We report the results of new smoothed particle hydrodynamics calculations of parabolic collisions between main-sequence (MS) stars. The stars are assumed to be close the MS turnoff point in a globular cluster and are therefore modeled as n = 3, Gamma = 5/3 polytropes. We find that the high degree of central mass concentration in these stars has a profound effect on the hydrodynamics. In particular, very little hydrodynamic mixing occurs between the dense, helium-rich inner cores and the outer envelopes. As a result, and in contrast to what has been assumed in previous studies, blue stragglers formed by direct stellar collisions are not necessarily expected to have anomalously high helium abundances in their envelopes or to have their cores replenished with fresh hydrogen fuel.

Lombardi, James, C. jr.↗

NASA NTRS1995

The orbital eccentricities of binary millisecond pulsars in globular clusters

Low-mass binary millisecond pulsars (LMBPs) are born with very small orbital eccentricities, typically of order e(sub i) approximately 10(exp -6) to 10(exp -3). In globular clusters, however, higher eccentricities e(sub f) much greater than e(sub i) can be induced by dynamical interactions with passing stars. Here we show that the cross section for this process is much larger than previously estimated. This is becuse, even for initially circular binaries, the induced eccentricity e(sub f) for an encounter with pericenter separation r(sub p) beyond a few times the binary semimajor axis a declines only as a power law (e(sub f) varies as (r(sub p)/a)(exp -5/2), and not as an exponential. We find that all currently known LMBPs in clusters were probably affected by interactions, with their current eccentricities typically greater than at birth by an order of magnitude or more.

Rasio, Frederic A.↗

NASA NTRS1995

The minimum mass ratio of W Ursae Majoris binaries

The minimum mass ratio for tidal stability of a contact binary containing two unevolved main-sequence stars is calculated to be q(sub min) approximately =0.09 in the case of a mostly radiative primary, and it is higher if an appreciable fraction of the mass lies in a convective envelope. At least one observed system, AW UMa, has a mass ratio just below this value (q = 0.075), implying that, if the system is stable, the primary must be slightly evolved and must have a very shallow convective envelope. Contact binaries with mass ratios significantly below that of AW UMa should not be observed, since they are tidally unstable and quickly merge into a single, rapidly rotating object, on a timescale approximately 10(exp 3)-10(exp 4) yr.

Rasio, Frederic A.↗

NASA NTRS1995

Hydrodynamics of binary coalescence. 2: Polytropes with gamma = 5/3

We present a new numerical study of the equilibrium and stability properties of close binary systems. We use the smoothed-particle hydrodynamics (SPH) technique both to construct accurate equilibrium configurations in three dimensions and to follow their hydrodynamic evolution. We adopt a simple polytropic equation of state p = K(sub rho)(exp gamma) with gamma = 5/3 and K = constant within each star, applicable to low-mass degenerate dwarfs as well as low-mass main-sequence stars. For degenerate configurations, we set the two polytropic constants equal, K = K prime, independent of the mass ratio. For main-sequence stars, we adjust K and K prime so as to obtain a simple mass-radius relation of the form R/R prime = M/M prime, where R prime and M prime are the radius and mass of the secondary. Along a sequence of binary equilibrium configurations for two identical stars, we demonstrate the existence of both secular and dynamical instabilities, confirming directly the results of recent analytic work. We use the SPH method to calculate the nonlinear development of the dynamical instability and to determine the final fate of the system. We find that the two stars merge together into a single, rapidly rotating object in just a few orbital periods. Equilibrium sequences are also constructed for systems containing two nonidentical stars. These sequences terminate at a Roche limit, which we can determine very accurately using SPH. For two low-mass main-sequence stars with mass ratio q approximately less than 0.4 we find that the (synchronized) Roche limit configuration is secularly unstable. We discuss the implications of our results for the evolution of double white-dwarf systems and W Ursae Majoris binaries.

Rasio, Frederic A.↗

NASA NTRS1995

Binary-binary interactions and the formation of the PSR B1620-26 triple system in M4

The hierarchical triple system containing the millisecond pulsar PSR B1620-26 in M4 is the first star system ever detected in a globular cluster. Such systems should form in globular clusters as a result of dynamical interactions between binaries. We propose that the triple system containing PSR B1620-26 formed through an exchange interaction between a wide primordial binary and a pre-existing binary millisecond pulsar. This scenario would have the advantage of reconciling the approximately 10(exp 9) yr timing age of the pulsar with the much shorter lifetime of the triple system in the core of M4.

Rasio, Frederic A.↗

NASA NTRS1994

Hydrodynamics of rotating stars and close binary interactions: Compressible ellipsoid models

We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids, and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest, such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.

Lai, Dong↗

NASA NTRS1994

A possible optical counterpart for the triple pulsar system PSR B1620-26

We have searched deep optical images of the globular cluster M4 for a possible optical counterpart for the third body in the PSR B1620-26 system. We identify as a possible optical counterpart a star located within 0.3 sec of the nominal pulsar position with V approximately equal to 20. This magnitude is consistent with a main-sequence cluster member with M approximately equal to 0.45 solar mass. However, this star may be blended with fainter stars, and the probability of a chance superposition is nonnegligible.

Bailyn, Charles D.↗

NASA NTRS1994

Collapse of primordial gas clouds and the formation of quasar black holes

The formation of quasar black holes during the hydrodynamic collapse of protogalactic gas clouds is discussed. The dissipational collapse and long-term dynamical evolution of these systems is analyzed using three-dimensional numerical simulations. The calculations focus on the final collapse stages of the inner baryonic component and therefore ignore the presence of dark matter. Two types of initial conditions are considered: uniformly rotating spherical clouds, and iirotational ellipsoidal clouds. In both cases the clouds are initially cold, homogeneous, and not far from rotational support (T/(absolute value of W) approximately equals 0.1). Although the details of the dynamical evolution depend sensitively on the initial conditions, the qualitative features of the final configurations do not. Most of the gas is found to fragment into small dense clumps, that eventually make up a spheroidal component resembling a galactic bulge. About 5% of the initial mass remains in the form of a smooth disk of gas supported by rotation in the gravitational potential potential well of the outer spheroid. If a central seed black hole of mass approximately greater than 10(exp 6) solar mass forms, it can grow by steady accretion from the disk and reach a typical quasar black hole mass approximately 10(exp 8) solar mass in less than 5 x 10(exp 8) yr. In the absence of a sufficiently massive seed, dynamical instabilities in a strongly self-gravitating inner region of the disk will inhibit steady accretion of gas and may prevent the immediate formation of quasar.

Loeb, Abraham↗

NASA NTRS1994

Hydrodynamics of binary coalescence. 1: Polytropes with stiff equations of state

We performed a series of three-dimensional hydrodynamic calculations of binary coalescence using the smoothed particle hydrodynamics (SPH) method. The initial conditions are exact polytropic equilibrium configurations on the verge of dynamical instability. We consider synchronized equilbria only and concentrate on stiff equations of state, with adiabatic Gamma greater than 5/3. We assume that the polytropic constants (K identically equal to P/(rho(exp Gamma)) are the same for both components. These conditions apply well to models of neutron star binaries. Accordingly, we discuss our results in the context of the Laser Interferometer Gravitational-Wave Observatory (LIGO) project, and we calculate the emission of gravitational radiation in the quadruple approximation. The fully nonlinear development of the instability is followed using SPH until a new equilibrium configuration is reached by the system. We find that the properties of this final configuration depend sensitively on both the compressibility and mass ratio. An axisymmetric merged configuration is always produced when the adiabatic exponent Gamma approximately less than 2.3. As a consequence, the emission of gravitational radiation shuts off abruptly right after the onset of dynamical instability. In contrast, triaxial merged configurations are obtained when Gamma approximately greater than 2.3, and the system continues to emit gravitational waves after the final coalescence. Systems with mass ratios q not equal to 1 typically become dynamically unstable before the onset of mass transfer. Stable mass transfer from one neutron star to another in a close binary is therefore probably ruled out. For a mass ratio q approximately less than 0.5, however, dynamical mass transfer can temporarily retard the coalescence by causing a rapid reexpansion of the binary into a new, slightly eccentric but dynamically stable orbit. The maximum amplitude h(sub max) and peak luminosity L(sub max) of the gravitational waves emitted during the final coalescence are nearly independent of Gamma, but depend sensitively on the mass ratio q. The approximate scalings we find are h(sub max) varies as q(exp 2) and L(sub max) varies as q(exp 6) for q close to unity.

Rasio, Frederic A.↗

NASA NTRS1994

Is there a planet int he PSR 1620-26 triple system?

The unusually large eccentricity (e(sub 1) = 0.025) of the low-mass binary millisecond pular PSR B1620-26 can be explained naturally as arising from the secular perturbation of a second, more distant companion. Such a triple configuration has been proposed recently as the most likely cause of the anomalous second period derivative of the pulsar. The current timing data are consistent with a second-companion mass m(sub 2) as low as approximately 10(exp -3) solar mass, i.e., comparable to that of Jupiter. However, if the eccentricity is indeed produced by secular perturbations, then the second companion must be another star, most likely of mass m(sub 2) less than or approximately equal to 1 solar mass and in a very eccentric (e(sub 2) greater than or approximately equal to 0.5) orbit of period P(sub 2) approximately equal 10(exp 2) - 10(exp 3) yr. A second companion of planetary mass cannot induce the observed eccentricity. Independent of the mass of the second companion, small changes in the binary pulsar's orbit should become detectable with just a few additional years of timing data. This detection would provide direct confirmation of the triple nature of the system, and an accurate measurement of the effects would place important new constraints on the orbital parameters.

Rasio, Frederic A.↗

NASA NTRS1994

Equilibrium, stability, and orbital evolution of close binary systems

We present a new analytic study of the equilibrium and stability properties of close binary systems containing polytropic components. Our method is based on the use of ellipsoidal trial functions in an energy variational principle. We consider both synchronized and nonsynchronized systems, constructing the compressible generalizations of the classical Darwin and Darwin-Riemann configurations. Our method can be applied to a wide variety of binary models where the stellar masses, radii, spins, entropies, and polytropic indices are all allowed to vary over wide ranges and independently for each component. We find that both secular and dynamical instabilities can develop before a Roche limit or contact is reached along a sequence of models with decreasing binary separation. High incompressibility always makes a given binary system more susceptible to these instabilities, but the dependence on the mass ratio is more complicated. As simple applications, we construct models of double degenerate systems and of low-mass main-sequence star binaries. We also discuss the orbital evoltuion of close binary systems under the combined influence of fluid viscosity and secular angular momentum losses from processes like gravitational radiation. We show that the existence of global fluid instabilities can have a profound effect on the terminal evolution of coalescing binaries. The validity of our analytic solutions is examined by means of detailed comparisons with the results of recent numerical fluid calculations in three dimensions.

Lai, Dong↗

NASA NTRS1994

Hydrodynamic instability and coalescence of binary neutron stars

We study the importance of hydrodynamic effects on the evolution of coalescing binary neutron stars. Using an approximate energy functional constructed from equilibrium solutions for polytropic binary configuration, we incorporate hydrodynamic effects into the calculation of the orbital decay driven by gravitational wave emission. In particular, we follow the transition between the quasi-static, secular decay of the orbit at large separation and the rapid dynamical evolution of confirgurations approaching contact. We show that a purely Newtonian hydrodynamic instability can significantly accelerate the coalescence at small separation. Such an instability occurs in all close binary configurations containing sufficiently incompressible stars. Calculations are performed for various neutron star masses, radii, and spins. The influence of the stiffness of the equation of state is also explored by varying the effective polytopic index. Typically, we find that the radial infall velocity just prior to contact is about 10% of the tangential orbital velocity. Once the stability limit is reached, the final evolution only takes another orbit. Post-Newtonian effects can move the stability limit to a larger binary separation, and may induce an even larger radial velocity. We also consider the possibiltiy of mass transfer form one neutron star to the other. We show that stable mass transfer is unlikely except when the mass of one of the components is very small (M is less than or approximately 0.4 solar mass) and the viscosity is high enough to maintain corotation. Otherwise, either the two stars come into contact or the dynamical instability sets in before a Roche limit can be reached.

Lai, Dong↗

NASA NTRS1993

Ellipsoidal figures of equilibrium - Compressible models

The results of Chandrasekhar (1969) are generalized to polytropes, using a formalism based on ellipsoidal energy variational principle to construct approximate stellar equilibrium solutions and study their stability. After reviewing the energy variational method and describing the approach, several equivalent stability conditions are established and secular vs. dynamical instabilities are discussed. Then, the equilibrium structure equations are derived for isolated, rotating polytropes, and axisymmetric configurations (compressible Maclaurin spheroids) are considered. Particular attention is given to triaxial configurations, either in a state of uniform rotation (generalizing the classical Jacobi ellipsoids) or with internal fluid motions of uniform vorticity (the compressible analogues of Riemann-S ellipsoids) and to the stability of these single star configurations. The compressible generalizations of the Roche and Roche-Riemann problems for a polytrope in orbit about a point-mass companion are solved, and the generalized Darwin problem for two identical polytropes in a binary is considered.

Lai, Dong↗

NASA NTRS1993

Hydrodynamic stellar interactions in dense star clusters

Highly detailed HST observations of globular-cluster cores and galactic nuclei motivate new theoretical studies of the violent dynamical processes which govern the evolution of these very dense stellar systems. These processes include close stellar encounters and direct physical collisions between stars. Such hydrodynamic stellar interactions are thought to explain the large populations of blue stragglers, millisecond pulsars, X-ray binaries, and other peculiar sources observed in globular clusters. Three-dimensional hydrodynamics techniques now make it possible to perform realistic numerical simulations of these interactions. The results, when combined with those of N-body simulations of stellar dynamics, should provide for the first time a realistic description of dense star clusters. Here I review briefly current theoretical work on hydrodynamic stellar interactions, emphasizing its relevance to recent observations.

Rasio, Frederic A.↗

NASA NTRS1993

Collisions and close encounters between massive main-sequence stars

Smooth-particle hydrodynamics is used to investigate collisions and close encounters between two massive main-sequence stars. Approximate scaling relations and fitting formulas are given for the amount of mass loss and for the critical impact parameters for capture or merging. The more massive stars, which have smaller ratios of specific binding energy to the square of escape velocity, are more easily disrupted in collisions. As the relative velocity increases, the orbital energy of the colliding stars can only be dissipated by shock heating, and the critical capture radius decreases much faster than predicted by linear theory. The average fractional mass loss per collision in a cluster is found not to depend sensitively on the stellar velocity dispersion.

Lai, Dong↗

NASA NTRS1993

Hydrodynamic instability and coalescence of close binary systems

The meaning of hydrodynamic instability of close binary systems is clarified, and the resulting large (but finite) value of the radial infall velocity near the end of the coalescence is estimated. Emphasis is placed on the case where the orbital decay is driven by gravitational radiation. A simple calculation of the orbital decay is presented for a simple generic model of a close binary.

Lai, Dong↗

NASA NTRS1993

Orbital evolution of the PSR1257+12 planetary system

The detection of orbital perturbation effects in the PSR1257+12 timing data would provide irrefutable confirmation that planets are indeed orbiting the pulsar. Here we give an overview of how perturbation effects are expected to affect the orbital elements of the two planets over the next few years. In particular, we give a simple calculation of resonant perturbations, including the nonlinear effects which could be important if sin i less than about 0.1. We also present a new analysis of the effects of close encounters between the two planets and we discuss their detectability.

Rasio, Frederic A.↗

NASA NTRS1992

Hydrodynamical evolution of coalescing binary neutron stars

The hydrodynamics of the final merging of two neutron stars and the corresponding gravitational wave emission is studied in detail. Various test calculations are presented, including the compressible Roche and Darwin problems and the head-on collision of two polytropes. A complete coalescence calculation is presented for the simplest case of two identical neutron stars, represented by Gamma = 2 polytropes, in a circular orbit, with their spins aligned and synchronized with the orbital rotation.

Rasio, Frederic A.↗