We describe an analytic model for an evolving protoplanetary disk driven by viscosity and a disk wind. The disk is heated by stellar irradiation and energy generated by viscosity. The evolution is controlled by 3 parameters: (i) the inflow velocity towards the central star at a reference distance and temperature, (ii) the fraction of this inflow caused by the disk wind, and (iii) the mass loss rate via the wind relative to the inward flux in the disk. The model gives the disk midplane temperature and surface density as a function of time and distance from the star. It is intended to provide an efficient way to calculate conditions in a protoplanetary disk for use in simulations of planet formation. In the model, disks dominated by viscosity spread radially while losing mass onto the star. Radial spreading is the main factor reducing the surface density in the inner disk. The disk mass remains substantial at late times. Temperatures in the inner region are high at early times due to strong viscous heating. Disks dominated by a wind undergo much less radial spreading and weaker viscous heating. These disks have a much lower mass at late times than purely viscous disks. When mass loss via a wind is significant, the surface density gradient in the inner disk becomes shallower, and the slope can become positive in extreme cases.

Conventionally, a habitable planet is one that can support liquid water on its surface. Habitability depends on temperature, which is set by insolation and the greenhouse effect, due mainly to CO2 and water vapor. The CO2 level is increased by volcanic outgassing, and decreased by continental and seafloor weathering. Here, I examine the climate evolution of Earth-like planets using a globally averaged climate model that includes both weathering types. Climate is sensitive to the relative contributions of continental and seafloor weathering, even when the total weathering rate is fixed. Climate also depends strongly on the dependence of seafloor weathering on CO2 partial pressure. Both these factors are uncertain. Earth-like planets have two equilibrium climate states: (i) an ice-free state where outgassing is balanced by both weathering types, and (ii) an ice-covered state where outgassing is balanced by seafloor weathering alone. The second of these has not been explored in detail before. For some planets, neither state exists, and the climate cycles between ice-covered and ice-free states. For some other planets, both equilibria exist, and the climate depends on the initial conditions. Insolation increases over time due to stellar evolution, so a planet usually encounters the ice-covered equilibrium first. Such a planet will remain ice-covered, even if the ice-free state appears subsequently, unless the climate receives a large perturbation. The ice-covered equilibrium state covers a large fraction of phase space for Earth-like planets. Many planets conventionally assigned to a star's habitable zone may be rendered uninhabitable as a result.

The solar system's terrestrial planets are thought to have accreted over millions of years out of a sea of smaller embryos and planetesimals. Because it is impossible to know the surface density profile for solids and size frequency distribution in the primordial solar nebula, distinguishing between the various proposed evolutionary schemes has been historically difficult. Nearly all previous simulations of terrestrial planet formation assume that Moon to Mars massed embryos formed throughout the inner solar system during the primordial gas-disk phase. However, validating this assumption through models of embryo accretion is computationally challenging because of the large number of bodies required. Here, we reevaluate this problem with GPU-accelerated, direct N-body simulations of embryo growth starting from r~100 km planetesimals. We find that embryos emerging from the primordial gas phase at a given radial distance already have masses similar to the largest objects at the same semi-major axis in the modern solar system. Thus, Earth and Venus attain ~50% of their modern mass, Mars-massed embryos form in the Mars region, and Ceres-massed objects are prevalent throughout asteroid belt. Consistent with other recent work, our new initial conditions for terrestrial accretion models produce markedly improved solar system analogs when evolved through the giant impact phase of planet formation. However, we still conclude that an additional dynamical mechanism such as giant planet migration is required to prevent Earth-massed Mars analogs from growing.

Previous studies have shown that planets that rotate retrograde (backwards with respect to their orbital motion) generally experience less severe obliquity variations than those that rotate prograde (the same direction as their orbital motion). Here we examine retrograde-rotating planets on eccentric orbits and find a previously unknown secular spin-orbit resonance that can drive significant obliquity variations. This resonance occurs when the frequency of the planet's rotation axis precession becomes commensurate with an orbital eigenfrequency of the planetary system. The planet's eccentricity enables a participating orbital frequency through an interaction in which the apsidal precession of the planet's orbit causes a cyclic nutation of the planet's orbital angular momentum vector. The resulting orbital frequency follows the relationship $f = 2 \dot{\varpi} - \dot{\Omega}$, where $\dot{\varpi}$ and $\dot{\Omega}$ are the rates of the planet's changing longitude of periapsis and ascending node, respectively. We test this mechanism by simulating cases of a simple Earth-Jupiter system, and confirm the predicted resonance. Over the course of 100 Myr, the test Earths with rotation axis precession rates near the predicted resonant frequency experienced pronounced obliquity variations of order $10^\circ$-$30^\circ$. These variations can be significant, and suggest that while retrograde rotation is a stabilizing influence most of the time, retrograde rotators can experience large obliquity variations if they are on eccentric orbits and enter this spin-orbit resonance.

Zechmeister et al. (2009) surveyed 38 nearby M dwarfs from March 2000 to March 2007 with VLT2 and the UVES spectrometer. This data has recently been reanalyzed (Butler et al. 2019), yielding a significant improvement in the Doppler velocity precision. Spurred by this, we have combined the UVES data with velocity sets from HARPS, Magellan/PFS, and Keck/HIRES. Sixteen planet candidates have been uncovered orbiting nine M dwarfs. Five of them are new planets corresponding to radial velocity signals, which are not sensitive to the choice of noise models and are identified in multiple data sets over various timespans. Eight candidate planets require additional observation to be confirmed. We also confirm three previously reported planets. Among the new planets, GJ 180 d and GJ 229A c are super-Earths located in the conservative habitable zones of their host stars. We investigate their dynamical stability using the Monte Carlo approach and find both planetary orbits are robust to the gravitational perturbations of the companion planets. Due to their proximity to the Sun, the angular separation between the host stars and the potentially habitable planets in these two systems is 25 and 59 mas, respectively. They are thus good candidates for future direct imaging by JWST and E-ELT. In addition we find GJ 433 c, a cold super-Neptune belonging to an unexplored population of Neptune-like planets. With a separation of 0.5 as from its host star, GJ 433 c is probably the first realistic candidate for direct imaging of cold Neptunes. A comprehensive survey of these planets is important for the studies of planet formation.

The nature of type Ia supernovae (SNIa) - thermonuclear explosions of white dwarf stars - is an open question in astrophysics. Virtually all existing theoretical models of normal, bright SNIa require the explosion to produce a detonation in order to consume all of stellar material, but the mechanism for the deflagration-to-detonation transition (DDT) remains unclear. We present a unified theory of turbulence-induced DDT that describes the mechanism and conditions for initiating detonation both in unconfined chemical and thermonuclear explosions. The model is validated using experiments with chemical flames and numerical simulations of thermonuclear flames. We use the developed theory to determine criteria for detonation initiation in the single degenerate Chandrasekhar-mass SNIa model, and show that DDT is almost inevitable at densities 10^7 - 10^8 g/cc.

Variations in the axial tilt, or obliquity, of terrestrial planets can affect their climates and therefore their habitability. Kepler-62f is a 1.4 R$_\oplus$ planet orbiting within the habitable zone of its K2 dwarf host star (Borucki et al. 2013). We perform N-body simulations that monitor the evolution of obliquity of Kepler-62f for 10 million year timescales to explore the effects on model assumptions, such as the masses of the Kepler-62 planets and the possibility of outer bodies. Significant obliquity variation occurs when the rotational precession frequency overlaps with one or more of the secular orbital frequencies, but most variations are limited to $\lesssim$10$^\circ$. Moderate variations ($\sim$10$^\circ - 20^\circ$) can occur over a broader range of initial obliquities when the relative nodal longitude ($\Delta\Omega$) overlaps with the frequency and phase of a given secular mode. However, we find that adding outer gas giants on long period orbits ($\gtrsim$ 1000 days) can produce large ($\sim$60$^\circ$) variations in obliquity if Kepler-62f has a very rapid (4 hr) rotation period. The possibility of giant planets on long period orbits impacts the climate and habitability of Kepler-62f through variations in the latitudinal surface flux, where the timescale for large variation can occur on million year timescales.

Of the solar system's four terrestrial planets, the origin of Mercury is perhaps the most mysterious. Modern numerical simulations designed to model the dynamics of terrestrial planet formation systematically fail to replicate Mercury; which possesses just 5% the mass of Earth and the highest orbital eccentricity and inclination among the planets. However, Mercury's large iron-rich core and low volatile inventory stand out among the inner planets, and seem to imply a violent collisional origin. Because most algorithms used for simulating terrestrial accretion do not consider the effects of collisional fragmentation, it has been difficult to test these collisional hypotheses within the larger context of planet formation. Here, we analyze a large suite of terrestrial accretion models that account for the fragmentation of colliding bodies. We find that planets with core mass fractions boosted as a result of repeated hit-and-run collisions are produced in 90% of our simulations. While many of these planets are similar to Mercury in mass, they rarely lie on Mercury-like orbits. Furthermore, we perform an additional batch of simulations designed to specifically test the single giant impact origin scenario. We find less than a 1% probability of simultaneously replicating the Mercury-Venus dynamical spacing and the terrestrial system's degree of orbital excitation after such an event. While dynamical models have made great strides in understanding Mars' low mass, their inability to form accurate Mercury analogs remains a glaring problem.

The solar system's dynamical state can be explained by an orbital instability among the giant planets. A recent model has proposed that the giant planet instability happened during terrestrial planet formation. This scenario has been shown to match the inner solar system by stunting Mars' growth and preventing planet formation in the asteroid belt. Here we present a large sample of new simulations of the "Early Instability" scenario. We use an N-body integration scheme that accounts for collisional fragmentation, and also perform a large set of control simulations that do not include an early giant planet instability. Since the total particle number decreases slower when collisional fragmentation is accounted for, the growing planets' orbits are damped more strongly via dynamical friction and encounters with small bodies that dissipate angular momentum (eg: hit-and-run impacts). Compared with simulations without collisional fragmentation, our fully evolved systems provide better matches to the solar system's terrestrial planets in terms of their compact mass distribution and dynamically cold orbits. Collisional processes also tend to lengthen the dynamical accretion timescales of Earth analogs, and shorten those of Mars analogs. This yields systems with relative growth timescales more consistent with those inferred from isotopic dating. Accounting for fragmentation is thus supremely important for any successful evolutionary model of the inner solar system.

Symplectic integrators separate a problem into parts that can be solved in isolation, alternately advancing these sub-problems to approximate the evolution of the complete system. Problems with a single, dominant mass can use mixed-variable symplectic (MVS) integrators that separate the problem into Keplerian motion of satellites about the primary, and satellite-satellite interactions. Here, we examine T+V algorithms where the problem is separated into kinetic T and potential energy V terms. T+V integrators are typically less efficient than MVS algorithms. This difference is reduced by using different step sizes for primary-satellite and satellite-satellite interactions. The T+V method is improved further using 4th and 6th-order algorithms that include force gradients and symplectic correctors. We describe three 6th-order algorithms, containing 2 or 3 force evaluations per step, that are competitive with MVS in some cases. Round-off errors for T+V integrators can be reduced by several orders of magnitude, at almost no computational cost, using a simple modification that keeps track of accumulated changes in the coordinates and momenta. This makes T+V algorithms desirable for long-term, high-accuracy calculations.

The physics of planet formation is investigated using a population synthesis approach. We develop a simple model for planetary growth including pebble and gas accretion, and orbital migration in an evolving protoplanetary disk. The model is run for a population of 2000 stars with a range of disk masses and radii, and initial protoplanet orbits. The resulting planetary distribution is compared with the observed population of extrasolar planets, and the model parameters are improved iteratively using a particle swarm optimization scheme. The characteristics of the final planetary systems are mainly controlled by the pebble isolation mass, which is the mass of a planet that perturbs nearby gas enough to halt the inward flux of drifting pebbles and stop growth. The pebble isolation mass increases with orbital distance so that giant planet cores can only form in the outer disk. Giants migrate inwards, populating a wide range of final orbital distances. The best model fits have large initial protoplanet masses, short pebble drift timescales, low disk viscosities, and short atmospheric cooling times, all of which promote rapid growth. The model successfully reproduces the observed frequency and distribution of giant planets and brown dwarfs. The fit for super Earths is poorer for single-planet systems, but improves steadily when more protoplanets are included. Although the study was designed to match the extrasolar planet distribution, analogs of the Solar System form in 1-2% of systems that contain at least 4 protoplanets.

In the core accretion model, gas-giant planets first form a solid core, which then accretes gas from a protoplanetary disk when the core exceeds a critical mass. Here, we model the atmosphere of a core that grows by accreting ice-rich pebbles. The ice fraction of pebbles evaporates in warm regions of the atmosphere, saturating it with water vapor. Excess water precipitates to lower altitudes. Beneath an outer radiative region, the atmosphere is convective, following a moist adiabat in saturated regions due to water condensation and precipitation. Atmospheric mass, density and temperature increase with core mass. For nominal model parameters, planets with core masses (ice + rock) between 0.08 and 0.16 Earth masses have surface temperatures between 273 K and 647 K and form an ocean. In more massive planets, water exists as a super-critical convecting fluid mixed with gas from the disk. Typically, the core mass reaches a maximum (the critical mass) as a function of the total mass when the core is 2-5 Earth masses. The critical mass depends in a complicated way on pebble size, mass flux, and dust opacity due to the occasional appearance of multiple core-mass maxima. The core mass for an atmosphere of 50 percent hydrogen and helium may be a more robust indicator of the onset of gas accretion. This mass is typically 1-3 Earth masses for pebbles that are 50 percent ice by mass, increasing with opacity and pebble flux, and decreasing with pebble ice/rock ratio.

Collisional fragmentation is shown to not be a barrier to rocky planet formation at small distances from the host star. Simple analytic arguments demonstrate that rocky planet formation via collisions of homogeneous gravity-dominated bodies is possible down to distances of order the Roche radius ($r_\mathrm{Roche}$). Extensive N-body simulations with initial bodies ${\gtrsim}1700$ km that include plausible models for fragmentation and merging of gravity-dominated bodies confirm this conclusion and demonstrate that rocky planet formation is possible down to ${\sim}$1.1 $r_\mathrm{Roche}$. At smaller distances, tidal effects cause collisions to be too fragmenting to allow mass build-up to a final, dynamically stable planetary system. We argue that even differentiated bodies can accumulate to form planets at distances that are not much larger than $r_\mathrm{Roche}$.

I examine the standard model of planet formation, including pebble accretion, using numerical simulations. Planetary embryos large enough to become giant planets do not form beyond the ice line within a typical disk lifetime unless icy pebbles stick at higher speeds than in experiments using rocky pebbles. Systems like the Solar System (small inner planets, giant outer planets) can form if (i) icy pebbles are stickier than rocky pebbles, and (ii) the planetesimal formation efficiency increases with pebble size, which prevents the formation of massive terrestrial planets. Growth beyond the ice line is dominated by pebble accretion. Most growth occurs early, when the surface density of pebbles is high due to inward drift of pebbles from the outer disk. Growth is much slower after the outer disk is depleted. The outcome is sensitive to the disk radius and turbulence level, which control the lifetime and maximum size of pebbles. The outcome is sensitive to the size of the largest planetesimals since there is a threshold mass for the onset of pebble accretion. The planetesimal formation rate is unimportant provided that some large planetesimals form while pebbles remain abundant. Two outcomes are seen, depending on whether pebble accretion begins while pebbles are still abundant. Either (i) multiple gas giant planets form beyond the ice line, small planets form close to the star, and a Kuiper-belt-like disk of bodies is scattered outwards by the giant planets; or (ii) no giants form and bodies remain Earth-mass or smaller.

The late stages of terrestrial planet formation are dominated by giant impacts that collectively influence the growth, composition and habitability of any planets that form. Hitherto, numerical models designed to explore these late stage collisions have been limited by assuming that all collisions lead to perfect accretion, and many of these studies lack the large number of realizations needed to account for the chaotic nature of N-body systems. We improve on these limitations by performing 280 simulations of planet formation around a Sun-like star, half of which used an N-body algorithm that has recently been modified to include fragmentation and hit-and-run (bouncing) collisions. We find that when fragmentation is included, the final planets formed are comparable in terms of mass and number, however their collision histories differ significantly and the accretion time approximately doubles. We explored impacts onto Earth-like planets which we parameterized in terms of their specific impact energies. Only 15 of our 164 Earth-analogs experienced an impact that was energetic enough to strip an entire atmosphere. To strip about half of an atmosphere requires energies comparable to recent models of the Moon-forming giant impact. Almost all Earth-analogs received at least one impact that met this criteria during the 2 Gyr simulations and the median was three giant impacts. The median time of the final giant impact was 43 Myr after the start of the simulations, leading us to conclude that the time-frame of the Moon-forming impact is typical amongst planetary systems around Sun-like stars.

Venus currently rotates slowly, with its spin controlled by solid-body and atmospheric thermal tides. However, conditions may have been far different 4 billion years ago, when the Sun was fainter and most of the carbon within Venus could have been in solid form, implying a low-mass atmosphere. We investigate how the obliquity would have varied for a hypothetical rapidly rotating Early Venus. The obliquity variation structure of an ensemble of hypothetical Early Venuses is simpler than that Earth would have if it lacked its large Moon (Lissauer et al., 2012), having just one primary chaotic regime at high prograde obliquities. We note an unexpected long-term variability of up to $\pm7^\circ$ for retrograde Venuses. Low-obliquity Venuses show very low total obliquity variability over billion-year timescales -- comparable to that of the real Moon-influenced Earth.

Many features of the outer solar system are replicated in numerical simulations if the giant planets undergo an orbital instability that ejects one or more ice giants. During this instability, Jupiter and Saturn's orbits diverge, crossing their 2:1 mean motion resonance (MMR), and this resonance-crossing can excite the terrestrial planet orbits. Using a large ensemble of simulations of this giant planet instability, we directly model the evolution of the terrestrial planet orbits during this process, paying special attention to systems that reproduce the basic features of the outer planets. In systems that retain four giant planets and finish with Jupiter and Saturn beyond their 2:1 MMR, we find at least an 85% probability that at least one terrestrial planet is lost. Moreover, systems that manage to retain all four terrestrial planets often finish with terrestrial planet eccentricities and inclinations larger than the observed ones. There is less than a ~5% chance that the terrestrial planet orbits will have a level of excitation comparable to the observed orbits. If we factor in the probability that the outer planetary orbits are well-replicated, we find a probability of 1% or less that the orbital architectures of the inner and outer planets are simultaneously reproduced in the same system. These small probabilities raise the prospect that the giant planet instability occurred before the terrestrial planets had formed. This scenario implies that the giant planet instability is not the source of the Late Heavy Bombardment and that terrestrial planet formation finished with the giant planets in their modern configuration.

This chapter concerns the long-term dynamical evolution of planetary systems from both theoretical and observational perspectives. We begin by discussing the planet-planet interactions that take place within our own Solar System. We then describe such interactions in more tightly-packed planetary systems. As planet-planet interactions build up, some systems become dynamically unstable, leading to strong encounters and ultimately either ejections or collisions of planets. After discussing the basic physical processes involved, we consider how these interactions apply to extrasolar planetary systems and explore the constraints provided by observed systems. The presence of a residual planetesimal disc can lead to planetary migration and hence cause instabilities induced by resonance crossing; however, such discs can also stabilise planetary systems. The crowded birth environment of a planetary system can have a significant impact: close encounters and binary companions can act to destabilise systems, or sculpt their properties. In the case of binaries, the Kozai mechanism can place planets on extremely eccentric orbits which may later circularise to produce hot Jupiters.

We conducted a search for very short-period transiting objects in the publicly available Kepler dataset. Our preliminary survey has revealed four planetary candidates, all with orbital periods less than twelve hours. We have analyzed the data for these candidates using photometric models that include transit light curves, ellipsoidal variations, and secondary eclipses to constrain the candidates' radii, masses, and effective temperatures. Even with masses of only a few Earth masses, the candidates' short periods mean they may induce stellar radial velocity signals (a few m/s) detectable by currently operating facilities. The origins of such short-period planets are unclear, but we discuss the possibility that they may be the remnants of disrupted hot Jupiters. Whatever their origins, if confirmed as planets, these candidates would be among the shortest-period planets ever discovered. Such planets would be particularly amenable to discovery by the planned TESS mission.

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $\beta$ Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While $\beta$ Pic b's near-IR colors mimick that of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects -- i.e. $\kappa$ And b and 1RXJ 1609B -- match $\beta$ Pic b's $JHK_{s}L^\prime$ photometry, and its 3.1 $\mu m$ and 5 $\mu m$ photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large ($\sim$ 60 $\mu m$) dust grains fail to reproduce the $\beta$ Pic b spectrum. However, models incorporating thick clouds similar to those found for HR 8799 bcde but also with small (a few microns) modal particle sizes yield fits consistent with the data within uncertainties. Assuming solar abundance models, thick clouds, and small dust particles ($<a>$ = 4 $\mu m$) we derive atmosphere parameters of log(g) = 3.8 $\pm$ 0.2 and $T_{eff}$ = 1575--1650 $K$, an inferred mass of 7$^{+4}_{-3}$ $M_{J}$, and a luminosity of log(L/L$_{\odot}$) $\sim$ -3.80 $\pm$ 0.02. The best-estimated planet radius, $\approx$ 1.65 $\pm$ 0.06 $R_{J}$, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges with constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if $\beta$ Pic b is younger than $\approx$ 7 Myr, consistent with a late formation, well after its host star's birth $\sim$ 12$^{+8}_{-4}$ Myr ago.

We report five new transit epochs of the extrasolar planet OGLE-TR-111b, observed in the v-HIGH and Bessell I bands with the FORS1 and FORS2 at the ESO Very Large Telescope, between April and May 2008. The new transits have been combined with all previously published transit data for this planet to provide a new Transit Timing Variations (TTVs) analysis of its orbit. We discard TTVs with amplitudes larger than 1.5 minutes over a 4-year observation time baseline, in agreement with the recent result by Adams et al.(2010a). Dynamical simulations fully exclude the presence of additional planets in the system with masses greater than 1.3, 0.4 and 0.5 M_earth at the 3:2, 1:2, 2:1 resonances, respectively. We also place an upper limit of about 30 M_earth on the mass of potential second planets in the region between the 3:2 and 1:2 mean-motion resonances.

Five new planets orbiting G and K dwarfs have emerged from the Magellan velocity survey. These companions are jovian-mass planets in eccentric (e \geq 0.24) intermediate and long-period orbits. HD 86226b orbits a solar metallicity G2 dwarf. The MP sin i mass of the planet is 1.5 MJUP, the semi-major axis is 2.6 AU, and the eccentricity 0.73. HD 129445b orbits a metal rich G6 dwarf. The minimum mass of the planet is MP sin i =1.6 MJUP, the semi-major axis is 2.9 AU, and the eccentricity 0.70. HD 164604b orbits a K2 dwarf. The MP sin i mass is 2.7 MJUP, semi-major axis is 1.3 AU, and the eccentricity is 0.24. HD 175167b orbits a metal rich G5 star. The MP sin i mass is 7.8 MJUP, the semi-major axis is 2.4 AU, and the eccentricity 0.54. HD 152079b orbits a G6 dwarf. The MP sin i mass of the planet is 3 MJUP, the semi-major axis is 3.2 AU, and the eccentricity is 0.60.

The Doppler technique measures the reflex radial motion of a star induced by the presence of companions and is the most successful method to detect exoplanets. If several planets are present, their signals will appear combined in the radial motion of the star, leading to potential misinterpretations of the data. Specifically, two planets in 2:1 resonant orbits can mimic the signal of a single planet in an eccentric orbit. We quantify the implications of this statistical degeneracy for a representative sample of the reported single exoplanets with available datasets, finding that 1) around 35 percent of the published eccentric one-planet solutions are statistically indistinguishable from planetary systems in 2:1 orbital resonance, 2) another 40 percent cannot be statistically distinguished from a circular orbital solution and 3) planets with masses comparable to Earth could be hidden in known orbital solutions of eccentric super-Earths and Neptune mass planets.

It has long been assumed that the planet Jupiter acts as a giant shield, significantly lowering the impact rate of small bodies on the Earth. However, until recently, very little work had been carried out examining the role played by Jupiter in determining the frequency of such collisions. In this work, the third of a series of papers, we examine the degree to which the impact rate on Earth resulting from the Oort cloud comets is enhanced or lessened by the presence of a giant planet in a Jupiter-like orbit, in an attempt to more fully understand the impact regime under which life on Earth has developed. Our results show that the presence of a giant planet in a Jupiter-like orbit significantly alters the impact rate of Oort cloud comets on the Earth, decreasing the rate as the mass of the giant increases. The greatest bombardment flus is observed when no giant planet is present.

We report low mass companions orbiting five Solar-type stars that have emerged from the Magellan precision Doppler velocity survey, with minimum (Msini) masses ranging from 1.2 to 25 Mjup. These nearby target stars range from mildly metal-poor to metal-rich, and appear to have low chromospheric activity. The companions to the brightest two of these stars have previously been reported from the CORALIE survey. Four of these companions (HD 48265-b, HD 143361-b, HD 28185-b, HD 111232-b) are low-mass Jupiter-like planets in eccentric intermediate and long-period orbits. On the other hand, the companion to HD 43848 appears to be a long period brown dwarf in a very eccentric orbit.

Most stars are formed in a cluster or association, where the number density of stars can be high. This means that a large fraction of initially-single stars will undergo close encounters with other stars and/or exchange into binaries. We describe how such close encounters and exchange encounters can affect the properties of a planetary system around a single star. We define a singleton as a single star which has never suffered close encounters with other stars or spent time within a binary system. It may be that planetary systems similar to our own solar system can only survive around singletons. Close encounters or the presence of a stellar companion will perturb the planetary system, often leaving planets on tighter and more eccentric orbits. Thus planetary systems which initially resembled our own solar system may later more closely resemble some of the observed exoplanet systems.

Symplectic integrators are the tool of choice for many researchers studying dynamical systems because of their good long-term energy conservation properties. For systems with a dominant central mass, symplectic integrators are also highly efficient. In this chapter, I describe the theory of symplectic integrators in terms of Lie series. I show how conventional symplectic algorithms have been adapted for use in binary-star systems to study problems such as the dynamical stability of multi-planet systems and the accretion of planets from planetesimals. This is achieved by devising new coordinate systems for the wide-binary and close-binary cases separately. I show how the performance of these algorithms can be improved at little extra cost using symplectic correctors. Finally, I discuss drawbacks of these algorithms, in particular in dealing with close encounters with one or both members of the binary, and the prospects for overcoming these problems.

In this letter we consider the evolution of a planetary system around a star inside a wide binary. We simulate numerically the evolution of the planetary orbits for both co-planar and highly-inclined systems. We find that the Kozai mechanism operates in the latter case. This produces a highly eccentric outer planet whose orbit crosses those of some of the inner planets. Strong planet-planet interactions then follow resulting in the ejection of one or more planets. We note that planetary systems resembling our solar system, formed around single stars in stellar clusters may exchange into binaries and thus will be vulnerable to planet stripping. This process will reduce the number of solar-system like planetary systems, and may produce at least some of the observed extra-solar planets.

We calculate herein the late stages of terrestrial planet accumulation around a solar type star that has a binary companion with semimajor axis larger than the terrestrial planet region. We perform more than one hundred simulations to survey binary parameter space and to account for sensitive dependence on initial conditions in these dynamical systems. As expected, sufficiently wide binaries leave the planet formation process largely unaffected. As a rough approximation, binary stars with periastron $q_B > 10$ AU have minimal effect on terrestrial planet formation within $\sim 2$ AU of the primary, whereas binary stars with $q_B \la$ 5 AU restrict terrestrial planet formation to within $\sim$ 1 AU of the primary star. Given the observed distribution of binary orbital elements for solar type primaries, we estimate that about 40 -- 50 percent of the binary population is wide enough to allow terrestrial planet formation to take place unimpeded. The large number of simulations allows for us to determine the distribution of results -- the distribution of plausible terrestrial planet systems -- for effectively equivalent starting conditions. We present (rough) distributions for the number of planets, their masses, and their orbital elements.

In the core-accretion model, gas-giant planets form solid cores which then accrete gaseous envelopes. Tidal interactions with disk gas cause a core to undergo inward type-I migration in 10^4 to 10^5 years. Cores must form faster than this to survive. Giant planets clear a gap in the disk and undergo inward type-II migration in <10^6 years if observed disk accretion rates apply to the disk as a whole. Type-II migration times exceed typical disk lifetimes if viscous accretion occurs mainly in the surface layers of disks. Low turbulent viscosities near the midplane may allow planetesimals to form by coagulation of dust grains. The radius r of such planetesimals is unknown. If r<0.5 km, the core formation time is shorter than the type-I migration timescale and cores will survive. Migration is substantial in most cases, leading to a wide range of planetary orbits, consistent with the observed variety of extrasolar systems. When r is of order 100m and midplane alpha is of order 3 times 10^-5, giant planets similar to those in the Solar System can form.

The mass and period ratios of the two planets orbiting 47 UMa suggest a possible kinship to the Jupiter-Saturn pair in our solar system. We explore the current dynamical state of the 47 UMa system with numerical integrations, and compare the results with analytic secular theory. We find that the planets in the system are likely participating in a secular resonance in which the difference in the longitudes of pericenter librates around zero. Alternately, it is possible that the system is participating in the 7:3 mean motion resonance. We show that stability considerations restrict the mutual inclination between the two planets to 40 degrees or less, and that this result is relatively insensitive to the total mass of the two planets. We present hydrodynamical simulations which measure the torques exerted on the planets by a hypothesized external protoplanetary disk. We show that planetary migration in response to torques from the disk may have led to capture of the system into a 7:3 mean-motion resonance, although it is unclear how the eccentricities of the planets would have been damped after capture occured. We show that Earth-mass planets can survive for long periods in some regions of the habitable zone of the nominal co-planar system. A set of planetary accretion calculations, however, shows that it is unlikely that large terrestrial planets can form in the 47 UMa habitable zone.

We explore the possibility that extrasolar planets might be found in the 1:1 mean-motion resonance. There are a variety of stable co-orbtial configurations, and we specifically examine three different versions of the 1:1 resonance. These include tadpole and horseshoe type orbits, as well as a more exotic configuration which occurs when one planet has a highly eccentric orbit while the other planet moves on a nearly circular orbit. We show that pairs of planets in 1:1 resonance yield characteristic radial velocity signatures which are not prone to the sin(i) degeneracy. Indeed, Keplerian fits to the radial velocities cannot reveal the presence of two planets in the 1:1 resonance. We discuss a dynamical fitting method for such systems, and illustrate its use with a simulated data set. Finally, we argue that hydrodynamical simulations and torqued three-body calculations indicate that 1:1 resonant pairs might readily form and migrate within protostellar disks.

We show that short-term perturbations among massive planets in multiple planet systems can result in radial velocity variations of the central star which differ substantially from velocity variations derived assuming the planets are executing independent Keplerian motions. We discuss two alternate fitting methods which can lead to an improved dynamical description of multiple planet systems. In the first method, the osculating orbital elements are determined via a Levenberg-Marquardt minimization scheme driving an N-body integrator. The second method is an improved analytic model in which orbital elements are allowed to vary according to a simple model for resonant interactions between the planets. Both of these methods can determine the true masses for the planets by eliminating the sin(i) degeneracy inherent in fits that assume independent Keplerian motions. We apply our fitting methods to the GJ876 radial velocity data (Marcy et al. 2001), and argue that the mass factors for the two planets are likely in the 1.25-2.0 range

Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly with order in timestep, rendering higher-order methods impractical. However, symplectic integrators are often applied to systems in which perturbations between bodies are a small factor of the force due to a dominant central mass. In this case, it is possible to create optimized symplectic algorithms that require fewer substeps per timestep. This is achieved by only considering error terms of order epsilon, and neglecting those of order epsilon^2, epsilon^3 etc. Here we devise symplectic algorithms with 4 and 6 substeps per step which effectively behave as 4th and 6th-order integrators when epsilon is small. These algorithms are more efficient than the usual 2nd and 4th-order methods when applied to planetary systems.