On the equivalence of the planet-satellite formation processes
Quantitative equivalence model in planet-satellite formation within Jupiter and Uranus systems
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Quantitative equivalence model in planet-satellite formation within Jupiter and Uranus systems
Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.
Formation of the planets in the Kepler-36 system is modeled by detailed numerical simulations according to the core-nucleated accretion scenario. The standard model is updated to include the dissolution of accreting rocky planetesimals in the gaseous envelope of the planet, leading to substantial enrichment of the envelope mass in heavy elements and a non-uniform composition with depth. For Kepler-36 c, models involving in situ formation and models involving orbital migration are considered. The results are compared with standard formation models. The calculations include the formation (accretion) phase, as well as the subsequent cooling phase, up to the age of Kepler-36 (7 Gyr). During the latter phase, mass loss induced by stellar XUV radiation is included. In all cases, the results fit the measured mass, 7.84 M, and radius, 3.68 R, of Kepler-36 c. Two parameters are varied to obtain these fits: the disk solid surface density at the formation location, and the "efficiency" factor in the XUV mass loss rate. The updated models are hotter and therefore less dense in the silicate portion of the planet and in the overlying layers of H/He, as compared with standard models. The lower densities mean that only about half as much H/He is needed to be accreted to fit the present-day mass and radius constraints. For Kepler-36 b, an updated in situ calculation shows that the entire H/He envelope is lost, early in the cooling phase, in agreement with observation.
ABSTRACT Many stars harbour multiplanet systems. As these stars expand late in their evolutions, the innermost planet may be engulfed, leading to a common envelope (CE) event. Even if this is insufficient to eject the envelope, it may expand the star further, causing additional CE events, with the last one unbinding what remains of the envelope. This multiplanet CE scenario may have broad implications for stellar and planetary evolution across a range of systems. We develop a simplified version and show that it may be able to explain the recently observed planet WD 1856 b.
A mechanism for the formation of the outer planets is proposed, the basis of which is the idea that the giant planets contain an excess of chemically condensable materials over solar composition. Planetary cores were formed by the clumping together of chemically condensed bodies forming a thin disk in the solar nebula. Gas surrounding a core becomes unstable against collapse onto the core. In the case of Jupiter and Saturn, much of the collapsing gas goes into orbit about the formed planet, forming a relatively thin circumplanetary disks with differential rotation in the prograde sense. For Uranus and Neptune, the dynamical collapse mechanism is unlikely. A disk of gas around Uranus may have been formed during a collision of the protoplanet with a large body. The circumplanetary disks then form the basis for formation of satellite systems, in which the Goldreich-Ward instability mechanism plays a role.
Given our presently inadequate understanding of the turbulent state of the solar and planetary nebulae, we believe the way to make progress in satellite formation is to consider two end member models that avoid over-reliance on specific choices of the turbulence (alpha), which is essentially a free parameter. The first end member model postulates turbulence decay once giant planet accretion ends. If so, Keplerian disks must eventually pass through the quiescent phases, so that the survival of satellites (and planets) ultimately hinges on gap-opening. In this scenario, the criterion for gap-opening itself sets the value for the gas surface density of the satellite disk.
In order to create a coherent scenario of satellite formation. the source of the solids (rock-metal and ice) that will eventually make up the satellites must be considered. While it is customary to use a solar composition mixture with a gas/solid mass ratio of about 100, at the tail end of the formation of the giant planet (when satellite formation is thought to have taken place) the fraction of solids entrained in the gas (particles with sizes lower than the decoupling size about 1 m for typical nebula parameters) is likely to be significantly lower than cosmic. In particular, in the core accretion model of giant planet formation one expects low dust and rubble content at late times due to particle coagulation leading to a collisional distribution of particle sizes with most of the mas residing in objects 1 km or larger, which are not coupled to the gas and whose dynamics must be followed independently. As a result, flow of gas into circumplanetary orbits is not sufficient to constrain the mas available to form satellites.
We report kinematic and thermal signatures associated with the directly imaged protoplanet candidate in the Elias 2–24 disc. Using the DSHARP (Disc Substructures at High Angular Resolution Project) ALMA (Atacama Large Millimetre/submillimetre Array) observations of the 12 CO J = 2 − 1 line, we show that the disc kinematics are perturbed, with a detached CO emission spot at the location of the planet candidate and traces of spiral wakes, and also that the observed CO emission intensities require local heating. While the foreground extinction hides the velocity channels associated with the planet, preventing a planet mass estimate, the level of gas heating implied by the CO emission indicates the presence of a warm, embedded giant planet. Comparison with models shows that this could either be a ≳5 M Jup or a lower mass (≳2 M Jup ) but accreting protoplanet.
Research on a variety of dynamical processes relevant to the formation of planets, satellites and ring systems is discussed. The main focus is on studies of accretionary formation of early protoplanets using a numerical model, structures and evolution of ring systems and individual bodies within planetary rings, and theories of lunar origin.
The origin of the outer planets and their satellite systems is discussed within the framework of astrophysical theories of the origin of stars and the nature of viscous accretion disks. The nature of the outer planet systems is reviewed, and the gas-instability and core-instability models of the origin of the outer planets are described. Critical tests of the models are discussed, showing that the core-instability model is favored. The main phases of evolution that the giant planets underwent during and after their formation are outlined, and the source of their current excess luminosity is considered in detail. The origin of the regular and irregular satellites of the outer planets is discussed, focusing on their formation in viscous accretion disks or by collision, and the possibility of capture by planets.
The origin of the rings of Uranus and Neptune is considered by performing two types of stochastic simulations of the collisional history of small moons: Monte Carlo simulations in which only the largest surviving fragments from each disruption is followed, and a Markov chain approach which makes it possible to follow the size distribution from each disruption to arbitrarily small sizes. Results indicate that the population of small satellites around Uranus and Neptune have evolved through catastrophic fragmentation since the end of planet and satellite formation 3 to 4 billion years ago.
Current models of planet and satellite formation are marred by our lack of understanding regarding the turbulent state of accretion disks. According to the Rayleigh criterion, Keplerian disks are hydrodynamically stable. Indeed, it has been argued that a carefully designed Taylor Couette experiment shows stability in the case of positive radial gradients in specific angular momentum even for high Reynolds numbers [1], in agreement with numerical simulations which consistently show turbulence decay [2]. Other possible sources of turbulence may fail due to low ionization, may decay as the optical depth decreases due to dust coagulation, may involve unrealistic boundary conditions, or result in limited transport. The difficulty stems not only from the degree of turbulence, but also from the kind of turbulence, and whether it may be characterized by an alpha parameter.
We present the results of hydrodynamical simulations of the orbital evolution of planets undergoing runaway gas accretion in radiative discs. We consider accreting disc models with constant mass flux through the disc, and where radiative cooling balances the effect of viscous heating and stellar irradiation. We assume that 20–30 M(sub ⨁) giant planet cores are formed in the region where viscous heating dominates and migrate outward under the action of a strong entropy-related corotation torque. In the case where gas accretion is neglected and for an α viscous stress parameter α = 2 × 10(exp -3), we find evidence for strong dynamical torques in accreting discs with accretion rates M ̇ ≳ 7 × 10(exsp -8)M(sub ⨀) yr(exp -1). Their main effect is to increase outward migration rates by a factor of ∼2 typically. In the presence of gas accretion, however, runaway outward migration is observed with the planet passing through the zero-torque radius and the transition between the viscous heating and stellar heating dominated regimes. The ability for an accreting planet to enter a fast migration regime is found to depend strongly on the planet growth rate but can occur for values of the mass flux through the disc of M ̇ ≳ 5 × 10(exp -8)M(sub ⨀) yr(exp -1).We find that an episode of runaway outward migration can cause an accreting planet formed in the 5–10 au region to temporarily orbit at star–planet separations as large as∼60–70 au. However, increase in the amplitude of the Lindblad torque associated with planet growth plus change in the streamline topology near the planet systematically cause the direction of migration to be reversed. Subsequent evolution corresponds to the planet migrating inward rapidly until it becomes massive enough to open a gap in the disc and migrate in the type II regime. Our results indicate that a planet can reach large orbital distances under the combined effect of dynamical torques and gas accretion, but an alternative mechanism is required to explain the presence of massive planets on wide orbits.
We model the subnebulae of Jupiter and Saturn wherein satellite accretion took place. We expect a giant planet subnebula to be composed of an optically thick (given gaseous opacity) inner region inside of the planet's centrifugal radius (located at r(sub c, sup J) = l5R(sub J) for Jupiter and r(sub c, sup S) = 22R(sub S) for Saturn), and an optically thin, extended outer disk out to a fraction of the planet's Roche lobe, which we choose to be R(sub roche)/5 (located at approximately 150R(sub J) near the inner irregular satellites for Jupiter, and approximately 200R(sub S) near Phoebe for Saturn). This places Titan and Ganymede in the inner disk, Callisto and Iapetus in the outer disk, and Hyperion in the transition region. The inner disk is the leftover of the gas accreted by the protoplanet. The outer disk results from the solar torque on nebula gas flowing into the protoplanet during the time of giant planet gap opening. For the sake of specificity, we use a cosmic mixture 'minimum mass' model to constrain the gas densities of the inner disks of Jupiter and Saturn (and also Uranus). For the total mass of the outer disk we use the simple scaling M(sub disk) = M(sub P)tau(sub gap)/tau(sub acc), where M(sub P) is the mass of the giant planet, tau(sub gap) is the gap opening timescale, and tau(sub acc) is the giant planet accretion time. This gives a total outer disk mass of approximately 100M(sub Callisto) for Jupiter and possibly approximately 200M(sub Iapetus) for Saturn (which contain enough condensables to form Callisto and Iapetus respectively). Our model has Ganymede at a subnebula temperature of approximately 250 K and Titan at approximately 100 K. The outer disks of Jupiter and Saturn have constant temperatures of 130 K and 90 K respectively.
Using an optically thick inner disk and an extended, optically thin outer disk as described in Mosqueira and Estrada, we compute the torque as a function of position in the subnebula, and show that although the torque exerted on the satellite is generally negative, which leads to inward migration as expected, there are regions of the disk where the torque is positive. For our model these regions of positive torque correspond roughly to the locations of Callisto and Iapetus. Though the outer location of zero torque depends on the (unknown) size of the transition region between the inner and outer disks, the result that Saturn's is found much farther out (at approximately 3r(sub c, sup S) where r(sub c, sup S) is Saturn's centrifugal radius) than Jupiter's (at approximately 2r(sub c, sup J), where r(sub c, sup J) is Jupiter's centrifugal radius) is mostly due to Saturn's less massive outer disk, and larger Hill radius. For a satellite to survive in the disk the timescale of satellite migration must be longer than the timescale for gas dissipation. For large satellites (approximately 1000 km) migration is dominated by the gas torque. We consider the possibility that the feedback reaction of the gas disk caused by the redistribution of gas surface density around satellites with masses larger than the inertial mass causes a large drop in the drift velocity of such objects, thus improving the likelihood that they will be left stranded following gas dissipation. We adapt the inviscid inertial mass criterion to include gas drag, and m-dependent non-local deposition of angular momentum.
Aspects of satellite origin and evolution which bear on the formation of the solar system are addressed, with emphasis on the chemical evidence contained in satellites concerning the modification of the interstellar material which went to form the solar system. New results on the tidal evolution of outer-planet satellite systems are presented. The gross physical and chemical properties of the satellite systems of the outer planets are outlined. Models for regular satellite formation are reviewed. Attention is given to the evolution of the satellites, with emphasis on tidal evolution, including a discussion of the dynamical origin of Triton. The volatile budgets of Titan, Triton, and Pluto/Charon are discussed, and models are presented for the origin of these budgets which tie them to the larger set of processes associated with solar system formation.
Solar system origin and planetary formation are discussed with emphasis placed on accretion disk dynamics, disk instabilities, giant gaseous protoplanets, condensation, sedimentation, coagulation, planetesimal swarm evolution, giant planet formation, and implications for satellites. Disk formation and the dynamics of a protosatellite disk are considered as well as satellite accretion, impact disruption and ablation, and satellite capture. Possible explanations for each of the satellite systems are offered. It is concluded that satellite formation involves a variety of processes.
The streaming instability is a promising mechanism to induce the formation of planetesimals. Nonetheless, this process has been found in previous studies to require either a dust-to-gas surface density ratio or a dust size that is enhanced compared to observed values. Employing two-dimensional global simulations of protoplanetary disks, we show that the vertical shear instability and the streaming instability in concert can cause dust concentration that is sufficient for planetesimal formation for lower surface density ratios and smaller dust sizes than the streaming instability in isolation, and in particular under conditions that are consistent with observational constraints. This is because dust overdensities forming in pressure bumps induced by the vertical shear instability act as seeds for the streaming instability and are enhanced by it. While our two-dimensional model does not include self-gravity, we find that strong dust clumping and the formation (and dissolution) of gravitationally unstable overdensities can be robustly inferred from the evolution of the maximum or the mean dust-to-gas volume density ratio. The vertical shear instability puffs up the dust layer to an average mid-plane dust-to-gas density ratio that is significantly below unity. We therefore find that reaching a mid-plane density ratio of one is not necessary to trigger planetesimal formation via the streaming instability when it acts in unison with the vertical shear instability.