Aaron C. Boley

Clumps in the Outer Disk by Disk Instability: Why They are Initially Gas Giants and the Legacy of Disruption

Abstract We explore the initial conditions for fragments in the extended regions ( r ≳ 50 AU ) of gravitationally unstable disks. We combine analytic estimates for the fragmentation of spiral arms with 3D SPH simulations to show that initial fragment masses are in the gas giant regime. These initial fragments will have substantial angular momentum, and should form disks with radii of a few AU. We show that clumps will survive for multiple orbits before they undergo a second, rapid collapse due to H 2 dissociation and that it is possible to destroy bound clumps by transporting them into the inner disk. The consequences of disrupted clumps for planet formation, dust processing, and disk evolution are discussed. We argue that it is possible to produce Earth-mass cores in the outer disk during the earliest phases of disk evolution.

Constraining the Planetary System of Fomalhaut Using High-resolution ALMA Observations

The dynamical evolution of planetary systems leaves observable signatures in debris disks. Optical images trace micron-sized grains, which are strongly affected by stellar radiation and need not coincide with their parent body population. Observations of millimeter-sized grains accurately trace parent bodies, but previous images lack the resolution and sensitivity needed to characterize the rings morphology. Here we present ALMA 350 GHz observations of the Fomalhaut debris ring. These observations demonstrate that the parent body population is 13-19 AU wide with a sharp inner and outer boundary. We discuss three possible origins for the ring and suggest that debris confined by shepherd planets is the most consistent with the rings morphology.

CHONDRULE FORMATION IN BOW SHOCKS AROUND ECCENTRIC PLANETARY EMBRYOS

Recent isotopic studies of Martian meteorites by Dauphas and Pourmand have established that large ({approx}3000 km radius) planetary embryos existed in the solar nebula at the same time that chondrules-millimeter-sized igneous inclusions found in meteorites-were forming. We model the formation of chondrules by passage through bow shocks around such a planetary embryo on an eccentric orbit. We numerically model the hydrodynamics of the flow and find that such large bodies retain an atmosphere with Kelvin-Helmholtz instabilities allowing mixing of this atmosphere with the gas and particles flowing past the embryo. We calculate the trajectories of chondrules flowing past the body and find that they are not accreted by the protoplanet, but may instead flow through volatiles outgassed from the planets magma ocean. In contrast, chondrules are accreted onto smaller planetesimals. We calculate the thermal histories of chondrules passing through the bow shock. We find that peak temperatures and cooling rates are consistent with the formation of the dominant, porphyritic texture of most chondrules, assuming a modest enhancement above the likely solar nebula average value of chondrule densities (by a factor of 10), attributable to settling of chondrule precursors to the midplane of the disk or turbulent concentration. We calculate themorexa0» rate at which a planetary embryos eccentricity is damped and conclude that a single planetary embryo scattered into an eccentric orbit can, over {approx}10{sup 5} years, produce {approx}10{sup 24} g of chondrules. In principle, a small number (1-10) of eccentric planetary embryos can melt the observed mass of chondrules in a manner consistent with all known constraints.«xa0less

ON THE POSSIBILITY OF ENRICHMENT AND DIFFERENTIATION IN GAS GIANTS DURING BIRTH BY DISK INSTABILITY

We investigate the coupling between rock-size solids and gas during the formation of gas giant planets by disk fragmentation in the outer regions of massive disks. In this study, we use three-dimensional radiative hydrodynamic simulations and model solids as a spatial distribution of particles. We assume that half of the total solid fraction is in small grains and half in large solids. The former are perfectly entrained with the gas and set the opacity in the disk, while the latter are allowed to respond to gas drag forces, with the back reaction on the gas taken into account. To explore the maximum effects of gas-solid interactions, we first consider 10 cm size particles. We then compare these results to a simulation with 1 km size particles, which explores the low-drag regime. We show that (1) disk instability planets have the potential to form large cores due to aerodynamic capturing of rock-size solids in spiral arms before fragmentation; (2) temporary clumps can concentrate tens of M ⊕ of solids in very localized regions before clump disruption; (3) the formation of permanent clumps, even in the outer disk, is dependent on the grain-size distribution, i.e., the opacity; (4) nonaxisymmetric structure in the disk can create disk regions that have a solids-to-gas ratio greater than unity; (5) the solid distribution may affect the fragmentation process; (6) proto-gas giants and proto-brown dwarfs can start as differentiated objects prior to the H2 collapse phase; (7) spiral arms in a gravitationally unstable disk are able to stop the inward drift of rock-size solids, even redistributing them to larger radii; and (8) large solids can form spiral arms that are offset from the gaseous spiral arms. We conclude that planet embryo formation can be strongly affected by the growth of solids during the earliest stages of disk accretion.

MIGRATION OF GAS GIANT PLANETS IN GRAVITATIONALLY UNSTABLE DISKS

Characterization of migration in gravitationally unstable disks is necessary to understand the fate of protoplanets formed by disk instability. As part of a larger study, we are using a three-dimensional radiative hydrodynamics code to investigate how an embedded gas giant planet interacts with a gas disk that undergoes gravitational instabilities (GIs). This Letter presents results from simulations with a Jupiter-mass planet placed in orbit at 25 AU within a 0.14 M{sub sun} disk. The disk spans 5-40 AU around a 1 M{sub sun} star and is initially marginally unstable. In one simulation, the planet is inserted prior to the eruption of GIs; in another, it is inserted only after the disk has settled into a quasi-steady GI-active state, where heating by GIs roughly balances radiative cooling. When the planet is present from the beginning, its own wake stimulates growth of a particular global mode with which it strongly interacts, and the planet plunges inward 6 AU in about 10{sup 3} years. In both cases with embedded planets, there are times when the planets radial motion is slow and varies in direction. At other times, when the planet appears to be interacting with strong spiral modes, migration both inward and outwardmorexa0» can be relatively rapid, covering several AUs over hundreds of years. Migration in both cases appears to stall near the inner Lindblad resonance of a dominant low-order mode. Planet orbit eccentricities fluctuate rapidly between about 0.02 and 0.1 throughout the GI-active phases of the simulations.«xa0less

The collapse of protoplanetary clumps formed through disc instability: 3D simulations of the pre-dissociation phase

We present 3D smoothed particle hydrodynamics simulations of the collapse of clumps formed through gravitational instability in the outer part of a protoplanetary disc. The initial conditions are taken directly from a global disc simulation, and a realistic equation of state is used to follow the clumps as they contract over several orders of magnitude in density, approaching the molecular hydrogen dissociation stage. The effects of clump rotation, asymmetries and radiative cooling are studied. Rotation provides support against fast collapse, but non-axisymmetric modes develop and efficiently transport angular momentum outwards, forming a circumplanetary disc. This transport helps the clump reach the dynamical collapse phase, resulting from molecular hydrogen dissociation, on a thousand-year time-scale, which is smaller than time-scales predicted by some previous spherical 1D collapse models. Extrapolation to the threshold of the runaway hydrogen dissociation indicates that the collapse time-scales can be shorter than inward migration time-scales, suggesting that clumps could survive tidal disruption and deliver a protogas giant to distances of even a few au from the central star.

Interactions between Moderate- and Long-period Giant Planets: Scattering Experiments for Systems in Isolation and with Stellar Flybys

The chance that a planetary system will interact with another member of its host stars nascent cluster would be greatly increased if gas giant planets form in situ on wide orbits. In this paper, we explore the outcomes of planet-planet scattering for a distribution of multi-planet systems that all have one of the planets on an initial orbit of 100 AU. The scattering experiments are run with and without stellar flybys. We convolve the outcomes with distributions for protoplanetary disk and stellar cluster sizes to generalize the results where possible. We find that the frequencies of large mutual inclinations and high eccentricities are sensitive to the number of planets in a system, but not strongly to stellar flybys. However, flybys do play a role in changing the low and moderate portions of the mutual inclination distributions, and erase dynamically cold initial conditions on average. Wide-orbit planets can be mixed throughout the planetary system, and in some cases, can potentially become hot Jupiters, which we demonstrate using scattering experiments that include a tidal damping model. If planets form in situ on wide orbits, then there will be discernible differences in the proper-motion distributions of a sample of wide-orbit planets compared withmorexa0» a pure scattering formation mechanism. Stellar flybys can enhance the frequency of ejections in planetary systems, but autoionization is likely to remain the dominant source of free-floating planets.«xa0less

THE HEAVY-ELEMENT COMPOSITION OF DISK INSTABILITY PLANETS CAN RANGE FROM SUB- TO SUPER-NEBULAR

Transit surveys combined with Doppler data have revealed a class of gas giant planets that are massive and highly enriched in heavy-elements (e.g., HD 149026b, GJ436b, and HAT-P-20b). It is tempting to consider these planets as validation of core accretion plus gas capture because it is often assumed that disk instability planets should be of nebular composition. We show in this paper, to the contrary, that gas giants that form by disk instability can have a variety of heavy-element compositions, ranging from sub- to super-nebular values. High levels of enrichment can be achieved through one or multiple mechanisms, including enrichment at birth, planetesimal capture, and differentiation plus tidal stripping. As a result, the metallicity of an individual gas giant cannot be used to discriminate between gas giant formation modes.

Chemistry in a gravitationally unstable protoplanetary disc

Until now, axisymmetric, α-disc models have been adopted for calculations of the chemical composition of protoplanetary discs. While this approach is reasonable for many discs, it is not appropriate when self-gravity is important. In this case, spiral waves and shocks cause temperature and density variations that affect the chemistry. We have adopted a dynamical model of a solar-mass star surrounded by a massive (0.39xa0M⊙), self-gravitating disc, similar to those that may be found around Class 0 and early Class I protostars, in a study of disc chemistry. We find that for each of a number of species, e.g. H2O, adsorption and desorption dominate the changes in the gas-phase fractional abundance; because the desorption rates are very sensitive to temperature, maps of the emissions from such species should reveal the locations of shocks of varying strengths. The gas-phase fractional abundances of some other species, e.g. CS, are also affected by gas-phase reactions, particularly in warm shocked regions. We conclude that the dynamics of massive discs have a strong impact on how they appear when imaged in the emission lines of various molecular species.

HIGH-TEMPERATURE PROCESSING OF SOLIDS THROUGH SOLAR NEBULAR BOW SHOCKS: 3D RADIATION HYDRODYNAMICS SIMULATIONS WITH PARTICLES

A fundamental, unsolved problem in solar system formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks have been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals and planetary embryos on eccentric orbits can produce bow shocks as they move supersonically through the disk gas, and are one possible source of chondrule-melting shocks. We investigate chondrule formation in bow shocks around planetoids through three-dimensional radiation hydrodynamics simulations. A new radiation transport algorithm that combines elements of flux-limited diffusion and Monte Carlo methods is used to capture the complexity of radiative transport around bow shocks. An equation of state that includes the rotational, vibrational, and dissociation modes of H2 is also used. Solids are followed directly in the simulations and their thermal histories are recorded. Adiabatic expansion creates rapid cooling of the gas, and tail shocks behind the embryo can cause secondary heating events. Radiative transport is efficient, and bow shocks around planetoids can have luminosities ~few× 10–8 L ☉. While barred and radial chondrule textures could be produced in the radiative shocks explored here, porphyritic chondrules may only be possible in the adiabatic limit. We present a series of predicted cooling curves that merit investigation in laboratory experiments to determine whether the solids produced by bow shocks are represented in the meteoritic record by chondrules or other solids.