Frederic A. Rasio

Dynamical Outcomes of Planet-Planet Scattering

Observations in the past decade have revealed extrasolar planets with a wide range of orbital semimajor axes and eccentricities. Based on the present understanding of planet formation via core accretion and oligarchic growth, we expect that giant planets often form in closely packed configurations. While the protoplanets are embedded in a protoplanetary gas disk, dissipation can prevent eccentricity growth and suppress instabilities from becoming manifest. However, once the disk dissipates, eccentricities can grow rapidly, leading to close encounters between planets. Strong planet-planet gravitational scattering could produce both high eccentricities and, after tidal circularization, very short period planets, as observed in the exoplanet population. We present new results for this scenario based on extensive dynamical integrations of systems containing three giant planets, both with and without residual gas disks. We assign the initial planetary masses and orbits in a realistic manner following the core accretion model of planet formation. We show that, with realistic initial conditions, planet-planet scattering can reproduce quite well the observed eccentricity distribution. Our results also make testable predictions for the orbital inclinations of short-period giant planets formed via strong planet scattering followed by tidal circularization.

SECULAR EVOLUTION OF HIERARCHICAL TRIPLE STAR SYSTEMS

We derive octupole-level secular perturbation equations for hierarchical triple systems, using classical Hamiltonian perturbation techniques. Our equations describe the secular evolution of the orbital eccentricities and inclinations over timescales that are long compared to the orbital periods. By extending previous work done to leading (quadrupole) order to octupole level (i.e., including terms of order α3, where α ≡ a1/a2 < 1 is the ratio of semimajor axes), we obtain expressions that are applicable to a much wider range of parameters. In particular, our results can be applied to high-inclination as well as coplanar systems, and our expressions are valid for almost all mass ratios for which the system is in a stable hierarchical configuration. In contrast, the standard quadrupole-level theory of Kozai gives a vanishing result in the limit of zero relative inclination. The classical planetary perturbation theory, while valid to all orders in α, applies only to orbits of low-mass objects orbiting a common central mass, with low eccentricities and low relative inclinations. For triple systems containing a close inner binary, we also discuss the possible interaction between the classical Newtonian perturbations and the general relativistic precession of the inner orbit. In some cases we show that this interaction can lead to resonances and a significant increase in the maximum amplitude of eccentricity perturbations. We establish the validity of our analytic expressions by providing detailed comparisons with the results of direct numerical integrations of the three-body problem obtained for a large number of representative cases. In addition, we show that our expressions reduce correctly to previously published analytic results obtained in various limiting regimes. We also discuss applications of the theory in the context of several observed triple systems of current interest, including the millisecond pulsar PSR B1620-26 in M4, the giant planet in 16 Cygni, and the protostellar binary TMR-1.

Hot Jupiters from secular planet-planet interactions

About 25 per cent of ‘hot Jupiters’ (extrasolar Jovian-mass planets with close-in orbits) are actually orbiting counter to the spin direction of the star. Perturbations from a distant binary star companion can produce high inclinations, but cannot explain orbits that are retrograde with respect to the total angular momentum of the system. Such orbits in a stellar context can be produced through secular (that is, long term) perturbations in hierarchical triple-star systems. Here we report a similar analysis of planetary bodies, including both octupole-order effects and tidal friction, and find that we can produce hot Jupiters in orbits that are retrograde with respect to the total angular momentum. With distant stellar mass perturbers, such an outcome is not possible. With planetary perturbers, the inner orbits angular momentum component parallel to the total angular momentum need not be constant. In fact, as we show here, it can even change sign, leading to a retrograde orbit. A brief excursion to very high eccentricity during the chaotic evolution of the inner orbit allows planet–star tidal interactions to rapidly circularize that orbit, decoupling the planets and forming a retrograde hot Jupiter.

Tidal decay of close planetary orbits

The 4.2-day orbit of the newly discovered planet around 51~Pegasi is formally unstable to tidal dissipation. However, the orbital decay time in this system is longer than the main-sequence lifetime of the central star. Given our best current understanding of tidal interactions, a planet of Jupiters mass around a solar-like star could have dynamically survived in an orbit with a period as short as

Origins of Eccentric Extrasolar Planets: Testing the Planet-Planet Scattering Model

\sim10\,

Compact Object Modeling with the StarTrack Population Synthesis Code

hr. Since radial velocities increase with decreasing period, we would expect to find those planets close to the tidal limit first and, unless this is a very unusual system, we would expect to find many more. We also consider the tidal stability of planets around more evolved stars and we re-examine in particular the question of whether the Earth can dynamically survive the red-giant phase in the evolution of the Sun.

Formation of Massive Black Holes in Dense Star Clusters. I. Mass Segregation and Core Collapse

In planetary systems with two or more giant planets, dynamical instabilities can lead to collisions or ejections through strong planet-planet scattering. Previous studies for initial conditions with two equal-mass planets revealed two discrepancies between the results of simulations and the observed orbits of exoplanets: potentially frequent collisions between giant planets and a narrow distribution of final eccentricities. We show that simulations with two unequal-mass planets starting on nearly circular orbits predict fewer collisions and a broader range of final eccentricities. Thus, the two-planet scattering model can reproduce the observed eccentricities with a plausible distribution of planet mass ratios. The model also predicts a maximum eccentricity of 0.8, independent of the distribution of planet mass ratios, provided that both planets are initially placed on nearly circular orbits. This compares favorably with observations and will be tested by future planet discoveries. Moreover, the combination of planet-planet scattering and tidal circularization may explain the existence of some giant planets with very short period orbits. Orbital migration due to planet scattering could play an important role in explaining the distribution of orbital periods found by radial velocity surveys. We also reexamine and discuss various possible correlations between eccentricities and other properties of exoplanets. We find that radial velocity observations are consistent with planet eccentricities being correlated with the ratio of the escape velocity from the planets surface relative to the escape velocity from the host star at the planets location. We demonstrate that the observed distribution of planet masses, periods, and eccentricities can provide constraints for models of planet formation and evolution.

Structure and Evolution of Nearby Stars with Planets. II. Physical Properties of ~1000 Cool Stars from the SPOCS Catalog

We present a comprehensive description of the population synthesis code StarTrack. The original code has been significantly modified and updated. Special emphasis is placed here on processes leading to the formation and further evolution of compact objects (white dwarfs, neutron stars, and black holes). Both single and binary star populations are considered. The code now incorporates detailed calculations of all mass transfer phases, a full implementation of orbital evolution due to tides, as well as the most recent estimates of magnetic braking. This updated version of StarTrack can be used for a wide variety of problems, with relevance to observations with many current and planned observatories, e.g., studies of X-ray binaries (Chandra, XMM-Newton), gravitational radiation sources (LIGO, LISA), and gamma-ray burst progenitors (HETE-II, Swift). The code has already been used in studies of Galactic and extragalactic X-ray binary populations, black holes in young star clusters, Type Ia supernova progenitors, and double compact object populations. Here we describe in detail the input physics, we present the code calibration and tests, and we outline our current studies in the context of X-ray binary populations.

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

We review possible dynamical formation processes for central massive black holes in dense star clusters. We focus on the early dynamical evolution of young clusters containing a few thousand to a few million stars. One natural formation path for a central seed black hole in these systems involves the development of the Spitzer instability, through which the most massive stars can drive the cluster to core collapse in a very short time. The sudden increase in the core density then leads to a runaway collision process and the formation of a very massive merger remnant, which must then collapse to a black hole. Alternatively, if the most massive stars end their lives before core collapse, a central cluster of stellar-mass black holes is formed. This cluster will likely evaporate before reaching the highly relativistic state necessary to drive a runaway merger process through gravitational radiation, thereby avoiding the formation of a central massive black hole. We summarize the conditions under which these different paths will be followed, and present the results of recent numerical simulations demonstrating the process of rapid core collapse and runaway collisions between massive stars.

Stellar collisions during binary-binary and binary-single star interactions

We derive detailed theoretical models for 1074 nearby stars from the SPOCS (Spectroscopic Properties of Cool Stars) Catalog. The California and Carnegie Planet Search has obtained high-quality (R 70,000-90,000, S/N 300-500) echelle spectra of over 1000 nearby stars taken with the Hamilton spectrograph at Lick Observatory, the HIRES spectrograph at Keck, and UCLES at the Anglo Australian Observatory. A uniform analysis of the high-resolution spectra has yielded precise stellar parameters (Teff, log g, v sin i, [M/H], and individual elemental abundances for Fe, Ni, Si, Na, and Ti), enabling systematic error analyses and accurate theoretical stellar modeling. We have created a large database of theoretical stellar evolution tracks using the Yale Stellar Evolution Code (YREC) to match the observed parameters of the SPOCS stars. Our very dense grids of evolutionary tracks eliminate the need for interpolation between stellar evolutionary tracks and allow precise determinations of physical stellar parameters (mass, age, radius, size and mass of the convective zone, surface gravity, etc.). Combining our stellar models with the observed stellar atmospheric parameters and uncertainties, we compute the likelihood for each set of stellar model parameters separated by uniform time steps along the stellar evolutionary tracks. The computed likelihoods are used for a Bayesian analysis to derive posterior probability distribution functions for the physical stellar parameters of interest. We provide a catalog of physical parameters for 1074 stars that are based on a uniform set of high-quality spectral observations, a uniform spectral reduction procedure, and a uniform set of stellar evolutionary models. We explore this catalog for various possible correlations between stellar and planetary properties, which may help constrain the formation and dynamical histories of other planetary systems.