Evolution of giant gaseous protoplanets embedded in the primitive solar nebula

Evolutionary calculations are presented, in a spherically symmetric approximation, for a protoplanet of 1 Jovian mass with homogeneous solar composition during the early phase of quasi-static contraction prior to the dissociation of molecular hydrogen. In contrast to earlier calculations which assume that protoplanets are isolated, this study invokes a time-dependent surface boundary condition that simulates physical conditions in an evolving primitive solar nebula. In a first set of calculations the protoplanet is surrounded by a “thermal bath” whose temperature varies with time and whose pressure is small and constant in time. Under a wide range of parameters the result is evaporation and complete dispersal of the object. Conditions required for the protoplanet to survive are discussed. In a second set of calculations both the temperature and pressure at the surface vary with time according to models of the solar nebula. In this case the protoplanet is not dispersed, but the evolution is accelerated or retarded relative to that of an isolated protoplanet, depending upon whether the entropy in the nebula is lower than or higher than, respectively, the entropy in the outer layers of the protoplanet. Processes by which terrestrial planets can form in the cores of giant gaseous protoplanets are discussed.     

The migration and growth of protoplanets in protostellar discs

We investigate the gravitational interaction of a Jovian-mass protoplanet with a gaseous disc with aspect ratio and kinematic viscosity expected for the protoplanetary disc from which it formed. Different disc surface density distributions are investigated. We focus on the tidal interaction with the disc with the consequent gap formation and orbital migration of the protoplanet. Non-linear two-dimensional hydrodynamic simulations are employed using three independent numerical codes.
A principal result is that the direction of the orbital migration is always inwards and such that the protoplanet reaches the central star in a near-circular orbit after a characteristic viscous time‐scale of ∼10⁴ initial orbital periods. This is found to be independent of whether the protoplanet is allowed to accrete mass or not. Inward migration is helped by the disappearance of the inner disc, and therefore the positive torque it would exert, because of accretion on to the central star. Maximally accreting protoplanets reach about 4 Jovian masses on reaching the neighbourhood of the central star. Our results indicate that a realistic upper limit for the masses of closely orbiting giant planets is ∼5 Jupiter masses, if they originate in protoplanetary discs similar to the minimum-mass solar nebula. This is because of the reduced accretion rates obtained for planets of increasing mass.
Assuming that some process such as termination of the inner disc through a magnetospheric cavity stops the migration, the range of masses estimated for a number of close orbiting giant planets as well as their inward orbital migration can be accounted for by consideration of disc–protoplanet interactions during the late stages of giant planet formation.     

eta Carinae: Binarity Confirmed

Gas giant planets have been detected in orbit around an increasing number of nearby stars. Two theories have been advanced for the formation of such planets: core accretion and disk instability. Core accretion, the generally accepted mechanism, requires several million years or more to form a gas giant planet in a protoplanetary disk like the solar nebula. Disk instability, on the other hand, can form a gas giant protoplanet in a few hundred years. However, disk instability has previously been thought to be important only in relatively massive disks. New three-dimensional, "locally isothermal," hydrodynamical models without velocity damping show that a disk instability can form Jupiter-mass clumps, even in a disk with a mass (0.091 M middle dot in circle within 20 AU) low enough to be in the range inferred for the solar nebula. The clumps form with initially eccentric orbits, and their survival will depend on their ability to contract to higher densities before they can be tidally disrupted at successive periastrons. Because the disk mass in these models is comparable to that apparently required for the core accretion mechanism to operate, the models imply that disk instability could obviate the core accretion mechanism in the solar nebula and elsewhere.     

Structure and evolution of isolated giant gaseous protoplanets

The structure and evolution of isolated giant gaseous protoplanets in the mass range 0.3 to 4.5 Jovian masses is investigated. Under the assumptions of the calculations, the following properties are found: (1) The central region of protoplanets of mass less than about 1 Jovian mass is, at some evolutionary epoch, thermodynamically favorable to the liquification of major interstellar grain constituents. Grains in this region can grow and infall to form a planetary core in tens to hundreds of years. (2) All protoplanets studied are convective through-out most of their interior. This property is in contrast to Bodenheimer's fully radiative proto-Jupiter models. We attribute the difference to the use of improved opacities. The presence of convection has at least two important consequences. First, it can mix grains into the central regions during planetary core formation, possibly allowing a core of mass ～ 1 Earth mass to grow. Second, convection can transport angular momentum outward as the protoplanet quasi-statically contracts. (3) The thermal contraction time depends sensitively on the surface opacity (T < 200°K). This opacity is uncertain within a factor of 5. The contraction times imply that some protoplanets can remain stable against tidal disruption by the proto-Sun and solar nebula during core-forming stages.     

Orbital evolution and accretion of protoplanets tidally interacting with a gas disk: II. Solid surface density evolution with type-I migration

This paper investigates the surface density evolution of a planetesimal disk due to the effect of type-I migration by carrying out N-body simulation and through analytical method, focusing on terrestrial planet formation. The coagulation and the growth of the planetesimals take place in the abundant gas disk except for a final stage. A protoplanet excites density waves in the gas disk, which causes the torque on the protoplanet. The torque imbalance makes the protoplanet suffer radial migration, which is known as type-I migration. Type-I migration time scale derived by the linear theory may be too short for the terrestrial planets to survive, which is one of the major problems in the planet formation scenario. Although the linear theory assumes a protoplanet being in a gas disk alone, Kominami et al. [Kominami, J., Tanaka, H., Ida, S., 2005. Icarus 167, 231-243] showed that the effect of the interaction with the planetesimal disk and the neighboring protoplanets on type-I migration is negligible. The migration becomes pronounced before the planet's mass reaches the isolation mass, and decreases the solid component in the disk. Runaway protoplanets form again in the planetesimal disk with decreased surface density. In this paper, we present the analytical formulas that describe the evolution of the solid surface density of the disk as a function of gas-to-dust ratio, gas depletion time scale and semimajor axis, which agree well with our results of N-body simulations. In general, significant depletion of solid material is likely to take place in inner regions of disks. This might be responsible for the fact that there is no planet inside Mercury's orbit in our Solar System. Our most important result is that the final surface density of solid components (Σd) and mass of surviving planets depend on gas surface density (Σg) and its depletion time scale (τdep) but not on initial Σd; they decrease with increase in Σg and τdep. For a fixed gas-to-dust ratio and τdep, larger initial Σd results in smaller final Σd and smaller surviving planets, because of larger Σg. To retain a specific amount of Σd, the efficient disk condition is not an initially large Σd but the initial Σd as small as the specified final one and a smaller gas-to-dust ratio. To retain Σd comparable to that of the minimum mass solar nebula (MMSN), a disk must have the same Σd and a gas-to-dust ratio that is smaller than that of MMSN by a factor of 1.3×(τdep/1 Myr) at ∼1 AU. (Equivalently, type-I migration speed is slower than that predicted by the linear theory by the same factor.) The surviving planets are Mars-sized ones in this case; in order to form Earth-sized planets, their eccentricities must be pumped up to start orbit crossing and coagulation among them. At ∼5 AU, Σd of MMSN is retained under the same condition, but to form a core massive enough to start runaway gas accretion, a gas-to-dust ratio must be smaller than that of MMSN by a factor of 3×τdep/1 Myr.     

On the evolution of giant protoplanets forming in circumbinary discs

We present the results of hydrodynamic simulations of Jovian mass protoplanets that form in circumbinary discs. The simulations follow the orbital evolution of the binary plus protoplanet system acting under their mutual gravitational forces, and forces exerted by the viscous circumbinary disc. The evolution involves the clearing of the inner circumbinary disc initially, so that the binary plus protoplanet system orbits within a low density cavity. Continued interaction between disc and protoplanet causes inward migration of the planet towards the inner binary. Subsequent evolution can take three distinct paths: (i) the protoplanet enters the 4 : 1 mean motion resonance with the binary, but is gravitationally scattered through a close encounter with the secondary star; (ii) the protoplanet enters the 4 : 1 mean motion resonance, the resonance breaks, and the planet remains in a stable orbit just outside the resonance; (iii) when the binary has initial eccentricity   e _bin≥ 0.2  , the disc becomes eccentric, leading to a stalling of the planet migration, and the formation of a stable circumbinary planet.
These results have implications for a number of issues in the study of extrasolar planets. The ejection of protoplanets in close binary systems provides a source of 'free-floating planets', which have been discovered recently. The formation of a large, tidally truncated cavity may provide an observational signature of circumbinary planets during formation. The existence of protoplanets orbiting stably just outside a mean motion resonance (4 : 1) in the simulations indicate that such sites may harbour planets in binary star systems, and these could potentially be observed. Finally, the formation of stable circumbinary planets in eccentric binary systems indicates that circumbinary planets may not be uncommon.     

Growth rate of a migrating protoplanet

We investigated the collision rate (i.e., the growth rate) of a migrating protoplanet with planetesimals. The collision rate strongly depends on the orbital elements of planetesimals (e.g., their eccentricities and inclinations). Thus we calculated the orbital evolutions of 2000 planetesimals in the vicinity of the migrating protoplanet and obtained the collision rate by counting the number of collisions with the protoplanet. For slow migration, the protoplanet makes a gap around its orbit in the planetesimals disk. On the other hand, for rapid migration, the protoplanet cannot shepherd planetesimals and keeps catching planetesimals. The obtained collision rate becomes larger with an increase in the migration speed. The comparison of the obtained collision rates with that of the previous work with no migration shows that the rapid migration of a protoplanet can enhance the collision rate by more than the factor 10. Using the obtained collision rate, we examined the growth of a migrating protoplanet. Our results suggest that, due to the enhancement of the collision rate, planets can be formed before they fall to the sun.     

Planet formation process as a phase transition

We discuss an adiabatic phase transition between a nebular ring and a protoplanet, taking into account the gravitational corrections to the perfect gas entropy. We compare the results of this model with those found in a previous paper from different assumptions, and we try a qualitative comparison with the structure of the real solar system, using the data for the four giant planets. We discuss the relevance of some primordial phenomena in the formation process of these planets such as a mass depletion within the asteroidal belt, a mass loss from the external regions of the nebula and a decrease of Neptune's orbital radius due to ejection of cometary material. The first and the second process seem to be significantly supported by the results of the phase transition model.     

Planetesimal capture in the disk instability model

We follow the contraction and evolution of a typical Jupiter-mass clump created by the disk instability mechanism, and compute the rate of planetesimal capture during this evolution. We show that such a clump has a slow contraction phase lasting ∼3×¹⁰⁵ yr. By following the trajectories of planetesimals as they pass through the envelope of the protoplanet, we compute the cross-section for planetesimal capture at all stages of the protoplanet's evolution. We show that the protoplanet can capture a large fraction of the solid material in its feeding zone, which will lead to an enrichment of the protoplanet in heavy elements. The exact amount of this enrichment depends upon, but is not very sensitive to the size and random speed of the planetesimals.     

Radial diffusion rate of planetesimals due to gravitational encounters

We study the rate of radial diffusion of planetesimals due to mutual gravitational encounters under Hill’s approximations in the three-body problem. Planetesimals orbiting a central star radially migrate inward and outward as a result of mutual gravitational encounters and transfer angular momentum. We calculate the viscosity in a disk of equal-sized planetesimals due to their mutual gravitational encounters using three-body orbital integrations, and obtain a semianalytic expression that reproduces the numerical results. We find that the viscosity is independent of the velocity dispersion of planetesimals when the velocity dispersion is so small that Kepler shear dominates planetesimals’ relative velocities. On the other hand, in high-velocity cases where random velocities dominate the relative velocities, the viscosity is a decreasing function of the velocity dispersion, and is found to agree with previous estimates under the two-body approximation neglecting the solar gravity. We also calculate the rate of radial diffusion of planetesimals due to gravitational scattering by a massive protoplanet. Using these results, we discuss a condition for formation of nonuniform radial surface density distribution of planetesimals by gravitational perturbation of an embedded protoplanet.     

Mass distribution in the solar system

The present-day observed mass distribution in the solar system including the Sun is shown to be compatible with the idea of the splitting of a number of ring-shaped rotating clouds of particles in the equatorial plane of a single contracting nebula. The formation of such a nebula is discussed and it is inferred that during the course of contraction this nebula has remained a sphere of uniform density spinning with the Keplerian velocity of its surface layer. The mass of a planet is taken as the portion of this spherical solar nebula gained at the time of splitting by its gaseous ring of dimensions satisfying Roche and accretional limits.     

Circumplanetary disc properties obtained from radiation hydrodynamical simulations of gas accretion by protoplanets

We investigate the properties of circumplanetary discs formed in three-dimensional, self-gravitating radiation hydrodynamical models of gas accretion by protoplanets. We determine disc sizes, scaleheights, and density and temperature profiles for different protoplanet masses, in solar nebulae of differing grain opacities.
We find that the analytical prediction of circumplanetary disc radii in an evacuated gap  ( R _Hill/3)  from Quillen & Trilling yields a good estimate for discs formed by high-mass protoplanets. The radial density profiles of the circumplanetary discs may be described by power laws between   r ⁻²  and   r ^−3/2  . We find no evidence for the ring-like density enhancements that have been found in some previous models of circumplanetary discs. Temperature profiles follow a  ∼ r ^−7/10  power law regardless of protoplanet mass or nebula grain opacity. The discs invariably have large scaleheights  ( H / r > 0.2)  , making them thick in comparison with their encompassing circumstellar discs, and they show no flaring.     

Evaporation and condensation of giant interstellar clouds in a hot-gas environment

Gas phases of the interstellar medium (ISM) coexist locally, penetrate each other and mix by means of dynamical and plasmaphysical processes. E.g. heat conduction from the hot to the cooler gas leads to energy and mass exchange between the gas phases. Analytical solutions exist under which evaporation of cloudy material or condensation of hot gas onto the clouds' surface dominate. Since these results are derived for stationary and static conditions and under ideal assumptions, they do not necessarily hold for a dynamical ISM. On the other hand, the mass and energy exchange between the gas phases is of great importance for the energy budget of the ISM and by this influences the evolution of galaxies. This led us to investigate the evolution of interstellar clouds in a hot gas by means of numerical simulations. At first, we compare static models with the analytical results and found that interstellar clouds with parameters requiring analytically evaporation are, in contrast, accreting surrounding material if self-gravitation and cooling are implied. For the more realistic case, where clouds are embedded in a streaming hot gas, the models show that Kelvin-Helmholtz instability which leads to the disruption of the clouds is suppressed by heat conduction so that the clouds are stabilized to survive. This revised version was published online in September 2006 with corrections to the Cover Date.     

Protoplanetary core formation by rain-out of iron drops

Using the stochastic collection equation we find that the time scale for rain out of liquid iron in a Saturn mass protoplanet is rapid compared with other evolutionary time scales and hence iron protoplanetary core formation is inevitable. The survival of this core during subsequent protoplanetary evolution and the consequences of the rain-out on the evolution are also discussed.     

Calculations of the early evolution of Jupiter

The evolution of the protoplanet Jupiter is followed, using a hydrodynamic computer code with radiative energy transport. Jupiter is assumed to have formed as a subcondensation in the primitive solar nebula at a density just high enough for gravitational collapse to occur. The initial state has a density of 1.5 × 10^?11 g cm^?3 and a temperature of 43 K; the calculations are carried to an equilibrium state where the central density reaches 0.5 g cm^?3 and the central temperature reaches 2.5 × 10⁴ K. During the early part of the evolution the object contracts in quasi-hydrostatic equilibrium; later on hydrodynamic collapse occurs, induced by the dissociation of hydrogen molecules. After dissociation is complete, the planet regains hydrostatic equilibrium with a radius of a few times the present value. Further evolution beyond this point is not treated here; however the results are consistent with the existence of a high-luminosity phase shortly after the planet settles into its final quasistatic contraction.     

Disk Planet Interactions And Early Evolution in Young Planetary Systems

We study and review disk protoplanet interactions using local shearing box simulations. These suffer the disadvantage of having potential artefacts arising from periodic boundary conditions but the advantage, when compared to global simulations, of being able to capture much of the dynamics close to the protoplanet at high resolution for low computational cost. Cases with and without self sustained MHD turbulence are considered. The conditions for gap formation and the transition from type I migration are investigated and found to depend on whether the single parameter M _p R ³/(M_* H ³), with M _p, M_*, R, and H being the protoplanet mass, the central mass, the orbital radius and the disk semi-thickness, respectively, exceeds a number of order unity. We also investigate the coorbital torques experienced by a moving protoplanet in an inviscid disk. This is done by demonstrating the equivalence of the problem for a moving protoplanet to one where the protoplanet is in a fixed orbit which the disk material flows through radially as a result of the action of an appropriate external torque. For sustainable coorbital torques to be realized a quasi steady state must be realized in which the planet migrates through the disk without accreting significant mass. In that case, although there is sensitivity to computational parameters, in agreement with earlier work by Masset and Papaloizou [2003, ApJ, 588, 494] based on global simulations, the coorbital torques are proportional to the migration speed and result in a positive feedback on the migration, enhancing it and potentially leading to a runaway. This could lead to fast migration for protoplanets in the Saturn mass range in massive disks and may be relevant to the mass period correlation for extrasolar planets which gives a preponderance of sub Jovian masses at short orbital periods.     

Hydrodynamic instability of the solar nebula in the presence of a planetary core

When a planetary core composed of condensed matter is accumulated in the primitive solar nebula, the gas of the nebula becomes gravitationally concentrated as an envelope surrounding the planetary core. Models of such gaseous envelopes have been constructed subject to the assumption that the gas everywhere is on the same adiabat as that in the surrounding nebula. The gaseous envelope extends from the surface of the core to the distance at which the gravitational attraction of core plus envelope becomes equal to the gradient of the gravitational potential in the solar nebula; at this point the pressure and temperature of the gas in the envelope are required to attain the background values characteristics of the solar nebula. In general, as the mass of the condensed core increases, increasing amounts of gas became concentrated in the envelope, and these envelopes are stable against hydrodynamic instabilities. However, the core mass then goes through a maximum and starts to decrease. In most of the models tested, the envelopes were hydrodynamically unstable beyond the peak in the core mass. An unstable situation was always created if it was insisted that the core mass contain a larger amount of matter than given by these solutions. For an initial adiabat characterized by a temperature of 450°K and a pressure of 5 × 10^?6 atm, the maximum core mass at which instability occurs is approximately 115 earth masses; this value is rather insensitive to the position in the solar nebula or to the background pressure of the solar nebula. However, if the adiabat is lowered, then the core mass corresponding to instability is decreased. Since the core masses found by Podolak and Cameron for the giant planets are significantly less than the critical core mass corresponding to the initial solar nebula adiabat, we conclude that the giant planets obtained their large amounts of hydrogen and helium by a hydrodynamic collapse process in the solar nebula only after the nebula had been subjected to a considerable period of cooling.     

Did iron meteorites form in the asteroid belt?—Evidence from thermodynamic models

The major iron meteorite groups are defined essentially by their Ga, Ge, and Ni contents. It now seems clear that the differences between their abundances of Ga and Ge were produced by the process of condensation and accretion in the primordial solar nebula. The simplest interpretation of the Ni abundance, and its variations between the groups, is also that it was fixed during condensation and accretion; more particularly, it reflects the oxidation state of the nebula during condensation and accretion. The abundance patterns of 17 other trace elements have been examined and are consistent with this model. It is believed to be the simplest model published and most consistent with analogous calculations for the chondrites. If it is correct, then the iron meteorite groups formed over a very wide range of pressures, 10^?4 to 10^?8 atm. Such a range could only be found in a restricted region of the nebula, such as the asteroid belt, if a complex accretion sequence inside a protoplanet occurred. More likely, the iron meteorites were formed in widely dispersed regions of the nebula and only one group formed in the asteroid belt, probably group IIIAB. Groups IAB and IIAB formed nearer the Sun, and group IVA formed much further out, say, beyond the orbit of Jupiter.     

Gas capture and rare gas retention by accreting planets in the solar nebula

In this paper, the physico-chemical effects of the nebula gas on the planets are reviewed from a standpoint of planetary formation in the solar nebula.The proto-Earth growing in the nebula was surrounded by a primordial atmosphere with a solar chemical composition and solar isotopic composition. When the mass of the proto-Earth was greater than 0.3 times the present Earth mass, the surface was molten because of the blanketing effect of the atmosphere. Therefore, the primordial rare gasses contained in the primordial atmosphere dissolved into the molten Earth material without fractionation and in particular the dissolved neon is expected to be conserved in the present Earth material. Hence, if dissolved neon with a solar isotopic ratio is discovered in the Earth material, it will indicate that the Earth was formed in the nebula and that the dissolved rare gases were one of the sources which degassed to form the present atmosphere.     

3. Solar System Formation and Early Evolution: the First 100 Million Years

The solar system, as we know it today, is about 4.5 billion years old. It is widely believed that it was essentially completed 100 million years after the formation of the Sun, which itself took less than 1 million years, although the exact chronology remains highly uncertain. For instance: which, of the giant planets or the terrestrial planets, formed first, and how? How did they acquire their mass? What was the early evolution of the “primitive solar nebula” (solar nebula for short)? What is its relation with the circumstellar disks that are ubiquitous around young low-mass stars today? Is it possible to define a “time zero” (t ₀), the epoch of the formation of the solar system? Is the solar system exceptional or common? This astronomical chapter focuses on the early stages, which determine in large part the subsequent evolution of the proto-solar system. This evolution is logarithmic, being very fast initially, then gradually slowing down. The chapter is thus divided in three parts: (1) The first million years: the stellar era. The dominant phase is the formation of the Sun in a stellar cluster, via accretion of material from a circumstellar disk, itself fed by a progressively vanishing circumstellar envelope. (2) The first 10 million years: the disk era. The dominant phase is the evolution and progressive disappearance of circumstellar disks around evolved young stars; planets will start to form at this stage. Important constraints on the solar nebula and on planet formation are drawn from the most primitive objects in the solar system, i.e., meteorites. (3) The first 100 million years: the “telluric” era. This phase is dominated by terrestrial (rocky) planet formation and differentiation, and the appearance of oceans and atmospheres.