On the stability of a planet between Mars and the asteroid belt: Implications for the Planet V hypothesis

The stability of an additional planet between the orbit of Mars and the asteroid belt is examined in the context of the Planet V hypothesis. In this model, the Solar System initially contained a fifth terrestrial planet, “Planet V,” which was removed after ∼700 Myr, a possible trigger for the late heavy bombardment on the inner planets. The model is investigated using 96 N-body integrations of the 8 major planets with an additional body between Mars and the asteroid belt. In more than 1/4 of simulations, Planet V survives for 1000 Myr. In most other cases, Planet V collides with the Sun or hits another planet after several hundred Myr, leaving 4 surviving terrestrial planets. In 24/96 simulations, Planet V is lost by ejection or collision with the Sun while the other four terrestrial planets survive without undergoing a collision. In 18 cases, Planet V is removed at least 200 Myr after the beginning of the simulation. The endstate depends sensitively on the mass of Planet V. Collision with the Sun is likely when Planet V's mass is 0.25 Mars masses or less. When Planet V is more massive than this, collisions involving it and/or other terrestrial planets become commonplace. In unstable systems, the times of first encounter and first collision/ejection depend on the initial aphelion distance of Mars. Reducing Mars's aphelion distance increases these times and also increases the fraction of systems surviving for 1000 Myr. When Mars's current orbit is used, the stability of Planet V increases when these two planets are widely separated initially. Planet V's aphelion distance Q typically begins to cross the asteroid belt within a few tens to a few hundred Myr, and its orbit last leaves the belt several hundred Myr later in most cases. The total time spent with Q>2.1 AU is typically less than 200 Myr.     

How common are Earth-Moon planetary systems?

The Earth’s comparatively massive moon, formed via a giant impact on the proto-Earth, has played an important role in the development of life on our planet, both in the history and strength of the ocean tides and in stabilizing the chaotic spin of our planet. Here we show that massive moons orbiting terrestrial planets are not rare. A large set of simulations by Morishima et al. (Morishima, R., Stadel, J., Moore, B. [2010]. Icarus. 207, 517-535), where Earth-like planets in the habitable zone form, provides the raw simulation data for our study. We use limits on the collision parameters that may guarantee the formation of a circumplanetary disk after a protoplanet collision that could form a satellite and study the collision history and the long term evolution of the satellites qualitatively. In addition, we estimate and quantify the uncertainties in each step of our study. We find that giant impacts with the required energy and orbital parameters for producing a binary planetary system do occur with more than 1 in 12 terrestrial planets hosting a massive moon, with a low-end estimate of 1 in 45 and a high-end estimate of 1 in 4.     

Effects of impacts on the atmospheric evolution: Comparison between Mars, Earth, and Venus

Classified as a terrestrial planet, Venus, Mars, and Earth are similar in several aspects such as bulk composition and density. Their atmospheres on the other hand have significant differences. Venus has the densest atmosphere, composed of CO₂ mainly, with atmospheric pressure at the planet's surface 92 times that of the Earth, while Mars has the thinnest atmosphere, composed also essentially of CO₂, with only several millibars of atmospheric surface pressure. In the past, both Mars and Venus could have possessed Earth-like climate permitting the presence of surface liquid water reservoirs. Impacts by asteroids and comets could have played a significant role in the evolution of the early atmospheres of the Earth, Mars, and Venus, not only by causing atmospheric erosion but also by delivering material and volatiles to the planets. Here we investigate the atmospheric loss and the delivery of volatiles for the three terrestrial planets using a parameterized model that takes into account the impact simulation results and the flux of impactors given in the literature. We show that the dimensions of the planets, the initial atmospheric surface pressures and the volatiles contents of the impactors are of high importance for the impact delivery and erosion, and that they might be responsible for the differences in the atmospheric evolution of Mars, Earth and Venus.     

Environmental predictions for the Beagle 2 lander, based on GCM climate simulations

The Mars climate database (MCD) is a database of statistics based on output from physically consistent numerical model simulations which describe the climate and surface environment of Mars. It is used here to predict the meteorological environment of the Beagle 2 lander site. The database was constructed directly on the basis of output from multiannual integrations of two general circulation models, developed jointly at Laboratoire de Météorologie Dynamique du Centre National de la Recherche Scientifique, France, and the University of Oxford, UK. In an atmosphere with dust opacities similar to that observed by Mars Global Surveyor, the predicted surface temperature at the time of landing (L_s=322°, 13:00 local time), is , and varying between ∼186 and over the Martian day. The predicted air temperature at above the surface, the height of the fully extended Beagle 2 robot arm, is at the time of landing. The expected mean wind near the surface on landing is north-westerly in direction, becoming more southerly over the mission. An increase in mean surface pressure is expected during the mission. Heavy global dust storm predictions are discussed; conditions which may only occur in the extreme as the expected time of landing is around the end of the main dust storm period. Past observations show approximately a one in five chance of a large-scale dust storm in a whole Mars year over the landing region, Isidis Planitia. This statistic results from observations of global, encircling, regional and local dust storms but does not include any small-scale dust “events” such as dust devils.     

The formation and habitability of terrestrial planets in the presence of close-in giant planets

‘Hot jupiters,’ giant planets with orbits very close to their parent stars, are thought to form farther away and migrate inward via interactions with a massive gas disk. If a giant planet forms and migrates quickly, the planetesimal population has time to re-generate in the lifetime of the disk and terrestrial planets may form [P.J. Armitage, A reduced efficiency of terrestrial planet formation following giant planet migration, Astrophys. J. 582 (2003) L47-L50]. We present results of simulations of terrestrial planet formation in the presence of hot/warm jupiters, broadly defined as having orbital radii ?0.5 AU. We show that terrestrial planets similar to those in the Solar System can form around stars with hot/warm jupiters, and can have water contents equal to or higher than the Earth's. For small orbital radii of hot jupiters (e.g., 0.15, 0.25 AU) potentially habitable planets can form, but for semi-major axes of 0.5 AU or greater their formation is suppressed. We show that the presence of an outer giant planet such as Jupiter does not enhance the water content of the terrestrial planets, but rather decreases their formation and water delivery timescales. We speculate that asteroid belts may exist interior to the terrestrial planets in systems with close-in giant planets.     

Eccentricities of Planets in Binary Systems

The most puzzling property of the extrasolar planets discovered by recent radial velocity surveys is their high orbital eccentricities, which are very difficult to explain within our current theoretical paradigm for planet formation. Current data reveal that at least 25% of these planets, including some with particularly high eccentricities, are orbiting a component of a binary star system. The presence of a distant companion can cause significant secular perturbations in the orbit of a planet. At high relative inclinations, large-amplitude, periodic eccentricity perturbations can occur. These are known as “Kozai cycles” and their amplitude is purely dependent on the relative orbital inclination. Assuming that every planet host star also has a (possibly unseen, e.g., substellar) distant companion, with reasonable distributions of orbital parameters and masses, we determine the resulting eccentricity distribution of planets and compare it to observations? We find that perturbations from a binary companion always appear to produce an excess of planets with both very high (?0.6) and very low (e ? 0.1) eccentricities. The paucity of near-circular orbits in the observed sample implies that at least one additional mechanism must be increasing eccentricities. On the other hand, the overproduction of very high eccentricities observed in our models could be combined with plausible circularization mechanisms (e.g., friction from residual gas) to create more planets with intermediate eccentricities (e? 0.1–0.6).     

The evolution of the water distribution in a viscous protoplanetary disk

Astronomical observations have shown that protoplanetary disks are dynamic objects through which mass is transported and accreted by the central star. This transport causes the disks to decrease in mass and cool over time, and such evolution is expected to have occurred in our own solar nebula. Age dating of meteorite constituents shows that their creation, evolution, and accumulation occupied several Myr, and over this time disk properties would evolve significantly. Moreover, on this timescale, solid particles decouple from the gas in the disk and their evolution follows a different path. It is in this context that we must understand how our own solar nebula evolved and what effects this evolution had on the primitive materials contained within it. Here we present a model which tracks how the distribution of water changes in an evolving disk as the water-bearing species experience condensation, accretion, transport, collisional destruction, and vaporization. Because solids are transported in a disk at different rates depending on their sizes, the motions will lead to water being concentrated in some regions of a disk and depleted in others. These enhancements and depletions are consistent with the conditions needed to explain some aspects of the chemistry of chondritic meteorites and formation of giant planets. The levels of concentration and depletion, as well as their locations, depend strongly on the combined effects of the gaseous disk evolution, the formation of rapidly migrating rubble, and the growth of immobile planetesimals. Understanding how these processes operate simultaneously is critical to developing our models for meteorite parent body formation in the Solar System and giant planet formation throughout the galaxy. We present examples of evolution under a range of plausible assumptions and demonstrate how the chemical evolution of the inner region of a protoplanetary disk is intimately connected to the physical processes which occur in the outer regions.     

Dynamic evolution model of isothermal voids and shocks

We explore self-similar hydrodynamic evolution of central voids embedded in an isothermal gas of spherical symmetry under the self-gravity. More specifically, we study voids expanding at constant radial speeds in an isothermal gas and construct all types of possible void solutions without or with shocks in surrounding envelopes. We examine properties of void boundaries and outer envelopes. Voids without shocks are all bounded by overdense shells and either inflows or outflows in the outer envelope may occur. These solutions, referred to as type void solutions, are further divided into subtypes and according to their characteristic behaviours across the sonic critical line (SCL). Void solutions with shocks in envelopes are referred to as type voids and can have both dense and quasi-smooth edges. Asymptotically, outflows, breezes, inflows, accretions and static outer envelopes may all surround such type voids. Both cases of constant and varying temperatures across isothermal shock fronts are analyzed; they are referred to as types and void shock solutions. We apply the ‘phase net matching procedure’ to construct various self-similar void solutions. We also present analysis on void generation mechanisms and describe several astrophysical applications. By including self-gravity, gas pressure and shocks, our isothermal self-similar void (ISSV) model is adaptable to various astrophysical systems such as planetary nebulae, hot bubbles and superbubbles in the interstellar medium as well as supernova remnants.      

Orbital stability of systems of closely-spaced planets

An investigation of the stability of systems of 1 M_⊕ (Earth-mass) bodies orbiting a Sun-like star has been conducted for virtual times reaching 10 billion years. For the majority of the tests, a symplectic integrator with a fixed timestep of between 1 and 10 days was employed; however, smaller timesteps and a Bulirsch-Stoer integrator were also selectively utilized to increase confidence in the results. In most cases, the planets were started on initially coplanar, circular orbits, and the longitudinal initial positions of neighboring planets were widely separated. The ratio of the semimajor axes of consecutive planets in each system was approximately uniform (so the spacing between consecutive planets increased slowly in terms of distance from the star). The stability time for a system was taken to be the time at which the orbits of two or more planets crossed. Our results show that, for a given class of system (e.g., three 1 M_⊕ planets), orbit crossing times vary with planetary spacing approximately as a power law over a wide range of separation in semimajor axis. Chaos tests indicate that deviations from this power law persist for changed initial longitudes and also for small but non-trivial changes in orbital spacing. We find that the stability time increases more rapidly at large initial orbital separations than the power-law dependence predicted from moderate initial orbital separations. Systems of five planets are less stable than systems of three planets for a specified semimajor axis spacing. Furthermore, systems of less massive planets can be packed more closely, being about as stable as 1 M_⊕ planets when the radial separation between planets is scaled using the mutual Hill radius. Finally, systems with retrograde planets can be packed substantially more closely than prograde systems with equal numbers of planets.     

Survival of Trojan-type companions of Neptune during primordial planet migration

We investigate the survivability of Trojan-type companions of Neptune during primordial radial migration of the giant planets Jupiter, Saturn, Uranus, and Neptune. We adopt the usual planet migration model in which the migration speed decreases exponentially with a characteristic time scale τ (the e-folding time). We perform a series of numerical simulations, each involving the migrating giant planets plus ∼1000 test particle Neptune Trojans with initial distributions of orbital eccentricity, inclination, and libration amplitude similar to those of the known jovian Trojans asteroids. We analyze these simulations to measure the survivability of Neptune's Trojans as a function of migration rate. We find that orbital migration with the characteristic time scale τ=10⁶ years allows about 35% of preexisting Neptune Trojans to survive to 5τ, by which time the giant planets have essentially reached their final orbits. In contrast, slower migration with τ=10⁷ years yields only a ∼5% probability of Neptune Trojans surviving to a time of 5τ. Interestingly, we find that the loss of Neptune Trojans during planetary migration is not a random diffusion process. Rather, losses occur almost exclusively during discrete prolonged episodes when Trojan particles are swept by secondary resonances associated with mean-motion commensurabilities of Uranus with Neptune. These secondary resonances arise when the circulation frequencies, f, of critical arguments for Uranus-Neptune mean-motion near-resonances (e.g., f^UN_1:2, f^UN_4:7) are commensurate with harmonics of the libration frequency of the critical argument for the Neptune-Trojan 1:1 mean-motion resonance (f^NT_1:1). Trojans trapped in the secondary resonances typically have their libration amplitudes amplified until they escape the 1:1 resonance with Neptune. Trojans with large libration amplitudes are susceptible to loss during sweeping by numerous high-order secondary resonances (e.g., f^UN_1:2≈11f^NT_1:1). However, for the slower migration, with τ=10⁷ years, even tightly bound Neptune Trojans with libration amplitudes below 10° can be lost when they become trapped in 1:3 or 1:2 secondary resonances between f^UN_1:2 and f^NT_1:1. With τ=10⁷ years the 1:2 secondary resonance was responsible for the single greatest episode of loss, ejecting nearly 75% of existing Neptune Trojans. This episode occurred during the late stages of planetary migration when the remnant planetesimal disk would have been largely dissipated. We speculate that if the number of bodies liberated during this event was sufficiently high they could have caused a spike in the impact rate throughout the Solar System.