Dynamics and Collisional Evolution of Closely Packed Planetary Systems

Dynamics and Collisional Evolution of Closely Packed Planetary Systems


Authors:

Hwang et al

Abstract:
High-multiplicity Kepler systems (referred to as Kepler multis) are often tightly packed and may be on the verge of instability. Many systems of this type could have experienced past instabilities, where the compact orbits and often low densities make physical collisions likely outcomes. We use numerical simulations to study the dynamical instabilities and planet-planet interactions in a synthetically generated sample of closely-packed, high-multiplicity systems. We focus specifically on systems resembling Kepler-11, a Kepler multi with six planets, and run a suite of dynamical integrations, sampling the initial orbital parameters around the nominal values reported in Lissauer et al. (2011a), finding that most of the realizations are unstable, resulting in orbit crossings and, eventually, collisions and mergers. We study in detail the dependence of stability on the orbital parameters of the planets and planet-pair characteristics to identify possible precursors to instability, compare the systems that emerge from dynamical instabilities to the observed Kepler sample (after applying observational corrections), and propose possible observable signatures of these instabilities. We examine the characteristics of each planet-planet collision, categorizing collisions by the degree of contact and collision energy, and find that grazing collisions are more common than direct impacts. Since the structure of many planets found in Kepler multis is such that the mass is dominated by a rocky core, but the volume is dominated by a low-density gaseous envelope, the sticky-sphere approximation may not be valid, and we present hydrodynamic calculations of planet-planet collisions clearly deviating from this approximation. Finally, we rerun a subset of our dynamical calculations using instead a modified prescription to handle collisions, finding, in general, higher multiplicity remnant systems.

Collisions between Sintered Icy Aggregates

Collisions between Sintered Icy Aggregates

Authors:

Sirono et al

Abstract:

Collisions between sintered icy dust aggregates are numerically simulated. If the temperature of an icy aggregate is sufficiently high, sintering promotes molecular transport and a neck between adjacent grains grows. This growth changes the mechanical responses of the neck. We included this effect in a simulation code, and conducted collisional simulations. For porous aggregates, the critical velocity for growth, below which the mass of an aggregate increases, decreased from 50 m s−1 for the non-sintered case to 20 m s−1. For compacted aggregates, the main collisional outcome is bouncing. These results come from the fact that the strength of the neck is increased by sintering. The numerical results suggest that the collisional growth of icy grain aggregates is strongly affected by sintering.

Exoplanets Around Binary Stars are Ejected, Don't Collide

Planet Scattering Around Binaries: Ejections, Not Collisions

Authors:

Smullen et al

Abstract:

Transiting circumbinary planets discovered by Kepler provide unique insight into binary star and planet formation. Several features of this new found population, for example the apparent pile-up of planets near the innermost stable orbit, may distinguish between formation theories. In this work, we determine how planet-planet scattering shapes planetary systems around binaries as compared to single stars. In particular, we look for signatures that arise due to differences in dynamical evolution in binary systems. We carry out a parameter study of N-body scattering simulations for four distinct planet populations around both binary and single stars. While binarity has little influence on the final system multiplicity or orbital distribution, the presence of a binary dramatically effects the means by which planets are lost from the system. Most circumbinary planets are lost due to ejections rather than planet-planet or planet-star collisions. The most massive planet in the system tends to control the evolution. Systems similar to the only observed multi-planet circumbinary system, Kepler-47, can arise from much more tightly packed, unstable systems. Only extreme initial conditions introduce differences in the final planet populations. Thus, we suggest that any intrinsic differences in the populations are imprinted by formation.

Constraints on Planetesimal Collision Models in Debris Disks

Constraints on Planetesimal Collision Models in Debris Disks

Authors:

MacGregor et al

Abstract:

Observations of debris disks offer a window into the physical and dynamical properties of planetesimals in extrasolar systems through the size distribution of dust grains. In particular, the millimeter spectral index of thermal dust emission encodes information on the grain size distribution. We have made new VLA observations of a sample of seven nearby debris disks at 9 mm, with 3" resolution and ∼5 μJy/beam rms. We combine these with archival ATCA observations of eight additional debris disks observed at 7 mm, together with up-to-date observations of all disks at (sub)millimeter wavelengths from the literature to place tight constraints on the millimeter spectral indices and thus grain size distributions. The analysis gives a weighted mean for the slope of the power law grain size distribution, n(a)∝a−q, of ⟨q⟩=3.36±0.02, with a possible trend of decreasing q for later spectral type stars. We compare our results to a range of theoretical models of collisional cascades, from the standard self-similar, steady-state size distribution (q=3.5) to solutions that incorporate more realistic physics such as alternative velocity distributions and material strengths, the possibility of a cutoff at small dust sizes from radiation pressure, as well as results from detailed dynamical calculations of specific disks. Such effects can lead to size distributions consistent with the data, and plausibly the observed scatter in spectral indices. For the AU Mic system, the VLA observations show clear evidence of a highly variable stellar emission component; this stellar activity obviates the need to invoke the presence of an asteroid belt to explain the previously reported compact millimeter source in this system.

Collision velocity of dust grains in self-gravitating protoplanetary discs

Collision velocity of dust grains in self-gravitating protoplanetary discs

Authors:

Booth et al

Abstract:

We have conducted the first comprehensive numerical investigation of the relative velocity distribution of dust particles in self-gravitating protoplanetary discs with a view to assessing the viability of planetesimal formation via direct collapse in such environments. The viability depends crucially on the large sizes that are preferentially collected in pressure maxima produced by transient spiral features (Stokes numbers, St∼1); growth to these size scales requires that collision velocities remain low enough that grain growth is not reversed by fragmentation. We show that, for a single sized dust population, velocity driving by the disc's gravitational perturbations is only effective for St>3, while coupling to the gas velocity dominates otherwise. We develop a criterion for understanding this result in terms of the stopping distance being of order the disc scale height. Nevertheless, the relative velocities induced by differential radial drift in multi-sized dust populations are too high to allow the growth of silicate dust particles beyond St∼10−2 or 10−1 (10cm to m sizes at 30au), such Stokes numbers being insufficient to allow concentration of solids in spiral features. However, for icy solids (which may survive collisions up to several 10ms−1), growth to St∼1 (10m size) may be possible beyond 30au from the star. Such objects would be concentrated in spiral features and could potentially produce larger icy planetesimals/comets by gravitational collapse. These planetesimals would acquire moderate eccentricities and remain unmodified over the remaining lifetime of the disc.

Massive Collision of Planetesimals Detected Around HD 61005

Massive collision of planetesimals in the asymmetric disk around HD61005

Authors:

Olofsson et al

Abstract:

Debris disks offer valuable insights into the latest stages of circumstellar disk evolution, and can possibly trace the outcomes of planetary formation processes. In the age range 10 to 100 Myr, most of the gas is expected to have been removed from the system, giant planets (if any) must have already been formed, and the formation of terrestrial planets may be on-going. Pluto-sized planetesimals, and their debris released in a collisional cascade, are under their mutual gravitational influence, which may result into non-axisymmetric structures in the debris disk. High angular resolution observations are required to investigate these effects and constrain the dynamical evolution of debris disks. Furthermore, multi-wavelength observations can provide information about the dust dynamics by probing different grain sizes. Here we present new VLT/SPHERE and ALMA observations of the debris disk around the 40 Myr old solar-type star HD 61005. We resolve the disk at unprecedented resolution both in the near-infrared (in scattered and polarized light) and at millimeter wavelengths. We perform a detailed modeling of these observations, including the spectral energy distribution. Thanks to the new observations, we propose a solution for both the radial and azimuthal distribution of the dust grains in the debris disk. We find that the disk has a moderate eccentricity (e∼0.1), and the dust density to be twice larger at the pericenter compared to the apocenter. We propose that the morphology of the disk could be the consequence of a massive collision between ∼ 1000 km-sized bodies at ∼ 61 AU, which could be the onset of in situ super-Earth formation. Numerical simulations have shown that such asymmetry can last for ∼1000 orbits (∼ 0.5 Myr at 61 AU).

Collisions of Solid ice in Planetesimal Formation

Collisions of solid ice in planetesimal formation

Authors:

Deckers et al

Abstract:

We present collision experiments of centimetre projectiles on to decimetre targets, both made up of solid ice, at velocities of 15–45 m s−1 at an average temperature of Tavg=255.8±0.7KTavg=255.8±0.7K. In these collisions, the centimetre body gets disrupted and part of it sticks to the target. This behaviour can be observed up to an upper threshold, that depends on the projectile size, beyond which there is no mass transfer. In collisions of small particles, as produced by the disruption of the centimetre projectiles, we also find mass transfer to the target. In this way, the larger body can gain mass, although the efficiency of the initial mass transfer is rather low. These collision results can be applied to planetesimal formation near the snowline, where evaporation and condensation is expected to produce solid ice. In free fall collisions at velocities up to about 7 m s−1, we investigated the threshold to fragmentation and coefficient of restitution of centimetre ice spheres.

Theia was not Unique: Uranus Formed From two Protoplanets Colliding

Stimulated Radiative Molecular Association in the Early Solar System. II. Orbital Radii of the Planets and Other Satellites of the Sun

Author:

Lombardi

Abstract:

In a previous investigation, the orbital radii of regular satellites of Uranus, Jupiter, Neptune, and Saturn are shown to be directly related to photon energies in the spectra of atomic and molecular hydrogen. To explain these observations a model was developed involving stimulated radiative molecular association (SRMA) reactions among photons and atoms in the protosatellite disks of the planets. In the present investigation, the previously developed model is applied to the planets and important satellites of the Sun. A key component of the model involves resonance associated with SRMA. Through this resonance, thermal energy is extracted from the protosun's protoplanetary disk at specific distances from the protosun wherever there is a match between the local thermal energy of the disk and the energy of photons impinging on the disk. Orbital radii of the planets and satellites are related to photon energies (EP values) in the spectrum of atomic hydrogen. An expression determined previously is used to relate EP values to temperatures in the disk. Results indicate the surface temperature of the protosun at the time when the evolution of the planets begins is higher than the surface temperature of a typical T Tauri star. The present investigation offers an explanation for the existence of the asteroid and classical Kuiper belts and predicts that a primordial belt once existed in the vicinity of Neptune. It also indicates that Uranus is formed from two protoplanets and is thus consistent with the theory that the large tilt of Uranus's axis of rotation was created by the collision of two bodies.

Resolution Dependence of Disruptive Collisions between Planetesimals in the Gravity Regime

Resolution Dependence of Disruptive Collisions between Planetesimals in the Gravity Regime

Authors:

Genda et al

Abstract:

Collisions are a fundamental process in planet formation. If colliding objects simply merge, a planetary object can grow. However, if the collision is disruptive, planetary growth is prevented. Therefore, the impact conditions under which collisions are destructive are important in understanding planet formation. So far, the critical specific impact energy for a disruptive collision Q_D^* has been investigated for various types of collisions between objects ranging in scale from centimeters to thousands of kilometers. Although the values of Q_D^* have been calculated numerically while taking into consideration various physical properties such as self-gravity, material strength, and porosity, the dependence of Q_D^* on numerical resolution has not been sufficiently investigated. In this paper, using the smoothed particle hydrodynamics (SPH) method, we performed numerical simulations of collisions between planetesimals at various numerical resolutions (from 2 x 10^4 to 5 x 10^6 SPH particles) and investigated the resulting variation in Q_D^*. The value of Q_D^* is shown to decrease as the number of SPH particles increases, and the difference between the Q_D^* values for the lowest and highest investigated resolutions is approximately a factor of two. Although the results for 5 x 10^6 SPH particles do not fully converge, higher-resolution simulations near the impact site show that the value of Q_D^* for the case with 5 x 10^6 SPH particles is close to the expected converged value. Although Q_D^* depends on impact parameters and material parameters, our results indicate that at least 5 x 10^6 SPH particles are required for numerical simulations in disruptive collisions to obtain the value of Q_D^* within 20% error.

Irregular Satellite Swarms Around Fomalhaut b

Collisional Cascade Caclulations for Irregular Satellite Swarms in Fomalhaut b

Authors:

Kenyon et al

Abstract:

We describe an extensive suite of numerical calculations for the collisional evolution of irregular satellite swarms around 1--300 M-earth planets orbiting at 120 AU in the Fomalhaut system. For 10--100 M-earth planets, swarms with initial masses of roughly 1% of the planet mass have cross-sectional areas comparable to the observed cross-sectional area of Fomalhaut b. Among 30--300 M-earth planets, our calculations yield optically thick swarms of satellites for ages of 1-10 Myr. Observations with HST and ground-based AO instruments can constrain the frequency of these systems around stars in the beta Pic moving group and possibly other nearby associations of young stars.

When Protoplanetary Disks are Gravitationally Unstable, Planetesimals can Collide, Grow Larger Than 1 Meter

ORBITAL SURVIVAL OF METER-SIZE AND LARGER BODIES DURING GRAVITATIONALLY UNSTABLE PHASES OF PROTOPLANETARY DISK EVOLUTION

Author:

Boss

Abstract:

A long-standing problem in the collisional accretion of terrestrial planets is the possible loss of m-size bodies through their inward migration onto the protostar as a result of gas drag forces. Such inward migration can be halted, and indeed even reversed, in a protoplanetary disk with local pressure maxima, such as marginally gravitationally unstable (MGU) phases of evolution, e.g., FU Orionis events. Results are presented for a suite of three-dimensional models of MGU disks extending from 1 to 10 AU and containing solid particles with sizes of 1 cm, 10 cm, 1 m, or 10 m, subject to disk gas drag and gravitational forces. These hydrodynamical models show that over disk evolution time scales of $\sim 6\times {10}^{3}$ years or longer, during which over half the gaseous disk mass is accreted by the protostar, very few 1 and 10 m bodies are lost through inward migration: most bodies survive and orbit stably in the outer disk. A greater fraction of 1 and 10 cm particles are lost to the central protostar during these time periods, as such particles are more closely tied to the disk gas accreting onto the protostar, but even in these cases, a significant fraction survive and undergo transport from the hot inner disk to the cold outer disk, perhaps explaining the presence of small refractory particles in Comet Wild 2. Evidently MGU disk phases offer a means to overcome the m-sized migration barrier to collisional accumulation.

NICE Model Predicts Comets are Fragmented Remnants of Larger Bodies

Comets as collisional fragments of a primordial planetesimal disk

Authors:

Morbidelli et al

Abstract:

The Rosetta mission and its exquisite measurements have revived the debate on whether comets are pristine planetesimals or collisionally evolved objects. We investigate the collisional evolution experienced by the precursors of current comet nuclei during the early stages of the Solar System, in the context of the so-called "Nice Model". We consider two environments for the collisional evolution: (1) the trans-planetary planetesimal disk, from the time of gas removal until the disk was dispersed by the migration of the ice giants, and (2) the dispersing disk during the time that the scattered disk was formed. Simulations have been performed, using different methods in the two cases, to find the number of destructive collisions typically experienced by a comet nucleus of 2km radius. In the widely accepted scenario, where the dispersal of the planetesimal disk occurred at the time of the Late Heavy Bombardment about 4Gy ago, comet-sized planetesimals have a very small chance to survive against destructive collisions in the disk. On the extreme assumption that the disk was dispersed directly upon gas removal, there is a chance for a significant fraction of the planetesimals to remain intact. However, these survivors would still bear the marks of many non-destructive impacts. Thus, the Nice Model of Solar System evolution predicts that typical km-sized comet nuclei are predominantly fragments resulting from collisions experienced by larger parent bodies. An important goal for further research is to investigate, whether the observed properties of comet nuclei are compatible with such a collisional origin.

Variations Detected in 6 Extreme Debris Disks Hint at Exoplanetary Collisions

Planetary Collisions outside the Solar System: Time Domain Characterization of Extreme Debris Disks

Authors:

Meng et al

Abstract:

Luminous debris disks of warm dust in the terrestrial planet zones around solar-like stars are recently found to vary, indicative of ongoing large-scale collisions of rocky objects. We use Spitzer 3.6 and 4.5 {\mu}m time-series observations in 2012 and 2013 (extended to 2014 in one case) to monitor 5 more debris disks with unusually high fractional luminosities ("extreme debris disk"), including P1121 in the open cluster M47 (80 Myr), HD 15407A in the AB Dor moving group (80 Myr), HD 23514 in the Pleiades (120 Myr), HD 145263 in the Upper Sco Association (10 Myr), and the field star BD+20 307 (>1 Gyr). Together with the published results for ID8 in NGC 2547 (35 Myr), this makes the first systematic time-domain investigation of planetary impacts outside the solar system. Significant variations with timescales shorter than a year are detected in five out of the six extreme debris disks we have monitored. However, different systems show diverse sets of characteristics in the time domain, including long-term decay or growth, disk temperature variations, and possible periodicity.

A Nemesis Class Jovian Detected in HD15115's Protoplanetary Disk

Near-Infrared Image of the Debris Disk around HD 15115

Authors:

Sai et al

Abstract:

We present a Subaru/IRCS H-band image of the edge-on debris disk around the F2V star HD 15115. We detected the debris disk, which has a bow shape and an asymmetric surface brightness, at a projected separation of 1--3" (~50--150 AU). The disk surface brightness is ~0.5--1.5 mag brighter on the western side than on the eastern side. We use an inclined annulus disk model to probe the disk geometry. The model fitting suggests that the disk has an inner hole with a radius of 86 AU and an eccentricity of 0.06. The disk model also indicates that the amount of dust on the western side is 2.2 times larger than that on the eastern side. A several Jupiter-mass planet may exist at ≳45 AU and capture grains at the Lagrangian points to open the eccentric gap. This scenario can explain both the eccentric gap and the difference in the amount of dust. In case of the stellar age of several 100 Myr, a dramatic planetesimal collision possibly causes the dust to increase in the western side. Interstellar medium interaction is also considered as a possible explanation of the asymmetric surface brightness, however, it hardly affect large grains in the vicinity of the inner hole.