-
Magnetically Driven Turbulence in the Inner Regions of Protoplanetary Disks
Authors:
David G. Rea,
Jacob B. Simon,
Daniel Carrera,
Geoffroy Lesur,
Wladimir Lyra,
Debanjan Sengupta,
Chao-Chin Yang,
Andrew N. Youdin
Abstract:
Given the important role turbulence plays in the settling and growth of dust grains in protoplanetary disks, it is crucial that we determine whether these disks are turbulent and to what extent. Protoplanetary disks are weakly ionized near the mid-plane, which has led to a paradigm in which largely laminar magnetic field structures prevail deeper in the disk, with angular momentum being transporte…
▽ More
Given the important role turbulence plays in the settling and growth of dust grains in protoplanetary disks, it is crucial that we determine whether these disks are turbulent and to what extent. Protoplanetary disks are weakly ionized near the mid-plane, which has led to a paradigm in which largely laminar magnetic field structures prevail deeper in the disk, with angular momentum being transported via magnetically launched winds. Yet, there has been little exploration on the precise behavior of the gas within the bulk of the disk. We carry out 3D, local shearing box simulations that include all three low-ionization effects (Ohmic diffusion, ambipolar diffusion, and the Hall effect) to probe the nature of magnetically driven gas dynamics 1-30 AU from the central star. We find that gas turbulence can persist with a generous yet physically motivated ionization prescription (order unity Elsasser numbers). The gas velocity fluctuations range from 0.03-0.09 of the sound speed $c_s$ at the disk mid-plane to $\sim c_s$ near the disk surface, and are dependent on the initial magnetic field strength. However, the turbulent velocities do not appear to be strongly dependent on the field polarity, and thus appear to be insensitive to the Hall effect. The mid-plane turbulence has the potential to drive dust grains to collision velocities exceeding their fragmentation limit, and likely reduces the efficacy of particle clumping in the mid-plane, though it remains to be seen if this level of turbulence persists in disks with lower ionization levels.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
A solution for the density dichotomy problem of Kuiper Belt objects with multi-species streaming instability and pebble accretion
Authors:
Manuel H. Cañas,
Wladimir Lyra,
Daniel Carrera,
Leonardo Krapp,
Debanjan Sengupta,
Jacob B. Simon,
Orkan M. Umurhan,
Chao-Chin Yang,
Andrew Youdin
Abstract:
Kuiper belt objects show an unexpected trend, whereby large bodies have increasingly higher densities, up to five times greater than their smaller counterparts. Current explanations for this trend assume formation at constant composition, with the increasing density resulting from gravitational compaction. However, this scenario poses a timing problem to avoid early melting by decay of $^{26}$Al.…
▽ More
Kuiper belt objects show an unexpected trend, whereby large bodies have increasingly higher densities, up to five times greater than their smaller counterparts. Current explanations for this trend assume formation at constant composition, with the increasing density resulting from gravitational compaction. However, this scenario poses a timing problem to avoid early melting by decay of $^{26}$Al. We aim to explain the density trend in the context of streaming instability and pebble accretion. Small pebbles experience lofting into the atmosphere of the disk, being exposed to UV and partially losing their ice via desorption. Conversely, larger pebbles are shielded and remain more icy. We use a shearing box model including gas and solids, the latter split into ices and silicate pebbles. Self-gravity is included, allowing dense clumps to collapse into planetesimals. We find that the streaming instability leads to the formation of mostly icy planetesimals, albeit with an unexpected trend that the lighter ones are more silicate-rich than the heavier ones. We feed the resulting planetesimals into a pebble accretion integrator with a continuous size distribution, finding that they undergo drastic changes in composition as they preferentially accrete silicate pebbles. The density and masses of large KBOs are best reproduced if they form between 15 and 22\,AU. Our solution avoids the timing problem because the first planetesimals are primarily icy, and $^{26}$Al is mostly incorporated in the slow phase of silicate pebble accretion. Our results lend further credibility to the streaming instability and pebble accretion as formation and growth mechanisms.
△ Less
Submitted 8 January, 2024;
originally announced January 2024.
-
Streaming Instability and Turbulence: Conditions for Planetesimal Formation
Authors:
Jeonghoon Lim,
Jacob B. Simon,
Rixin Li,
Philip J. Armitage,
Daniel Carrera,
Wladimir Lyra,
David G. Rea,
Chao-Chin Yang,
Andrew N. Youdin
Abstract:
The streaming instability (SI) is a leading candidate for planetesimal formation, which can concentrate solids through two-way aerodynamic interactions with the gas. The resulting concentrations can become sufficiently dense to collapse under particle self-gravity, forming planetesimals. Previous studies have carried out large parameter surveys to establish the critical particle to gas surface den…
▽ More
The streaming instability (SI) is a leading candidate for planetesimal formation, which can concentrate solids through two-way aerodynamic interactions with the gas. The resulting concentrations can become sufficiently dense to collapse under particle self-gravity, forming planetesimals. Previous studies have carried out large parameter surveys to establish the critical particle to gas surface density ratio ($Z$), above which SI-induced concentration triggers planetesimal formation. The threshold $Z$ depends on the dimensionless stopping time ($τ_s$, a proxy for dust size). However, these studies neglected both particle self-gravity and external turbulence. Here, we perform 3D stratified shearing box simulations with both particle self-gravity and turbulent forcing, which we characterize via $α_D$ that measures turbulent diffusion. We find that forced turbulence, at amplitudes plausibly present in some protoplanetary disks, can increase the threshold $Z$ by up to an order of magnitude. For example, for $τ_s = 0.01$, planetesimal formation occurs when $Z \gtrsim 0.06$, $\gtrsim 0.1$, and $\gtrsim 0.2$ at $α_D = 10^{-4}$, $10^{-3.5}$, and $10^{-3}$, respectively. We provide a single fit to the critical $Z$ as a function of $α_D$ and $τ_s$ required for the SI to work (though limited to the range $τ_s = 0.01$--0.1). Our simulations also show that planetesimal formation requires a mid-plane particle-to-gas density ratio that exceeds unity, with the critical value being independent of $α_D$. Finally, we provide the estimation of particle scale height that accounts for both particle feedback and external turbulence.
△ Less
Submitted 19 December, 2023;
originally announced December 2023.
-
The Streaming Instability Cannot Form Planetesimals from mm-size Grains in Pressure Bumps
Authors:
Daniel Carrera,
Jacob B. Simon
Abstract:
We present evidence that it is unlikely that the streaming instability (SI) can form planetesimals from mm grains inside axisymmetric pressure bumps. We conducted the largest simulation of the SI so far (7 million CPU hours), consisting of a large slice of the disk with mm grains, a solar-like dust-to-gas ratio ($Z = 0.01$), and the largest pressure bump that does not cause gravitational instabili…
▽ More
We present evidence that it is unlikely that the streaming instability (SI) can form planetesimals from mm grains inside axisymmetric pressure bumps. We conducted the largest simulation of the SI so far (7 million CPU hours), consisting of a large slice of the disk with mm grains, a solar-like dust-to-gas ratio ($Z = 0.01$), and the largest pressure bump that does not cause gravitational instability (GI) in the particle layer. We used a high resolution of $1000/H$ to resolve as many SI unstable modes as possible. The simulation produced a long-lived particle over-density far exceeding the SI criteria (i.e., a critical solid abundance to headwind parameter ratio $Z/Π$) where strong clumping would occur if these conditions were present over an extended region of the disk; yet we observed none. The likely reason is that the time it takes particles to cross the high-$Z/Π$ region ($t_{\rm cross}$) is shorter than the growth timescale of the SI ($t_{\rm grow}$). We propose an added criterion for planetesimal formation by the SI -- that $t_{\rm cross} > t_{\rm grow}$. We show that any bump larger than the one in this run would form planetesimals by the GI instead of the SI. Our results significantly restrict the pathways to planet formation: Either protoplanetary disks regularly form grains larger than 1~mm, or planetesimals do not form by the SI in axisymmetric pressure bumps. Since bumps large enough to induce the GI are likely Rossby-wave unstable, we propose that mm grains may only form planetesimals in vortices.
△ Less
Submitted 23 May, 2022; v1 submitted 29 April, 2022;
originally announced April 2022.
-
Using Bayesian Deep Learning to infer Planet Mass from Gaps in Protoplanetary Disks
Authors:
Sayantan Auddy,
Ramit Dey,
Min-Kai Lin,
Daniel Carrera,
Jacob B. Simon
Abstract:
Planet induced sub-structures, like annular gaps, observed in dust emission from protoplanetary disks provide a unique probe to characterize unseen young planets. While deep learning based model has an edge in characterizing the planet's properties over traditional methods, like customized simulations and empirical relations, it lacks in its ability to quantify the uncertainty associated with its…
▽ More
Planet induced sub-structures, like annular gaps, observed in dust emission from protoplanetary disks provide a unique probe to characterize unseen young planets. While deep learning based model has an edge in characterizing the planet's properties over traditional methods, like customized simulations and empirical relations, it lacks in its ability to quantify the uncertainty associated with its predictions. In this paper, we introduce a Bayesian deep learning network "DPNNet-Bayesian" that can predict planet mass from disk gaps and provides uncertainties associated with the prediction. A unique feature of our approach is that it can distinguish between the uncertainty associated with the deep learning architecture and uncertainty inherent in the input data due to measurement noise. The model is trained on a data set generated from disk-planet simulations using the \textsc{fargo3d} hydrodynamics code with a newly implemented fixed grain size module and improved initial conditions. The Bayesian framework enables estimating a gauge/confidence interval over the validity of the prediction when applied to unknown observations. As a proof-of-concept, we apply DPNNet-Bayesian to dust gaps observed in HL Tau. The network predicts masses of $ 86.0 \pm 5.5 M_{\Earth} $, $ 43.8 \pm 3.3 M_{\Earth} $, and $ 92.2 \pm 5.1 M_{\Earth} $ respectively, which are comparable to other studies based on specialized simulations.
△ Less
Submitted 23 February, 2022;
originally announced February 2022.
-
Migration traps as the root cause of the Kepler dichotomy
Authors:
Brianna Zawadzki,
Daniel Carrera,
Eric B. Ford
Abstract:
It is often assumed that the "Kepler dichotomy" -- the apparent excess of planetary systems with a single detected transiting planet in the Kepler catalog -- reflects an intrinsic bimodality in the mutual inclinations of planetary orbits. After conducting 600 simulations of planet formation followed by simulated Kepler observations, we instead propose that the apparent dichotomy reflects a diverge…
▽ More
It is often assumed that the "Kepler dichotomy" -- the apparent excess of planetary systems with a single detected transiting planet in the Kepler catalog -- reflects an intrinsic bimodality in the mutual inclinations of planetary orbits. After conducting 600 simulations of planet formation followed by simulated Kepler observations, we instead propose that the apparent dichotomy reflects a divergence in the amount of migration and the separation of planetary semimajor axes into distinct "clusters". We find that our simulated high-mass systems migrate rapidly, bringing more planets into orbital periods of less than 200 days. The outer planets are often caught in a migration trap -- a range of planet masses and locations in which a dominant co-rotation torque prevents inward migration -- which splits the system into two clusters. If clusters are sufficiently separated, the inner cluster remains dynamically cold, leading to low mutual inclinations and a higher probability of detecting multiple transiting planets. Conversely, our simulated low-mass systems typically bring fewer planets inside 200 days, forming a single cluster that quickly becomes dynamically unstable, leading to collisions and high mutual inclinations. We propose an alternative explanation for the apparent Kepler dichotomy in which migration traps during formation lead to fewer planets inside the Kepler detection window, and where mutual inclinations play only a secondary role. If our scenario is correct, then Kepler's STIPs (Systems with Tightly-packed Inner Planets) are a sample of planets that escaped capture by co-rotation traps, and their sizes may be a valuable probe into the structure of protoplanetary discs.
△ Less
Submitted 17 August, 2022; v1 submitted 10 February, 2022;
originally announced February 2022.
-
Resilience of Planetesimal Formation in Weakly-Reinforced Pressure Bumps
Authors:
Daniel Carrera,
Andrew Thomas,
Jacob B. Simon,
Matthew A. Small,
Katherine A. Kretke,
Hubert Klahr
Abstract:
The discovery that axisymmetric dust rings are ubiquitous in protoplanetary disks has provoked a flurry of research on the role of pressure bumps in planet formation. High-resolution simulations by our group have shown that even a modest bump can collect enough dust to trigger planetesimal formation by the streaming instability. In this work, we probe the limits of planetesimal formation when the…
▽ More
The discovery that axisymmetric dust rings are ubiquitous in protoplanetary disks has provoked a flurry of research on the role of pressure bumps in planet formation. High-resolution simulations by our group have shown that even a modest bump can collect enough dust to trigger planetesimal formation by the streaming instability. In this work, we probe the limits of planetesimal formation when the external source of pressure bump reinforcement is extremely weak. We conduct simulations of radially elongated shearing boxes to capture the entire bump, which is generated and maintained over some timescale $t_{\rm reinf}$ by a Newtonian relaxation scheme. We find that planetesimal formation is extremely resilient for cm-sized grains. We reduced the strength of reinforcement by up to a factor of 100 and the location and initial masses of planetesimals were essentially unaffected. However, we do find that strong reinforcement causes much faster pebble drift compared to the the standard pebble drift rates. The resulting larger pebble flux enhances the planetesimal growth rate by pebble accretion. We hypothesize that to sustain the bump, our code has to extract angular momentum (the strength of this negative torque depends on $t_{\rm reinf}$), and some of this torque is transferred to the particles, causing them to drift faster for a stronger torque (i.e., smaller $t_{\rm reinf}$). Since any physical process that sustains a pressure bump must do so by torquing the gas, we conjecture that the effect on pebble drift is a real phenomenon, motivating further work with physically realistic sources to generate the bump.
△ Less
Submitted 1 February, 2022; v1 submitted 18 August, 2021;
originally announced August 2021.
-
Rapid Formation of Super-Earths Around Low-Mass Stars
Authors:
Brianna Zawadzki,
Daniel Carrera,
Eric Ford
Abstract:
NASA's TESS mission is expected to discover hundreds of M dwarf planets. However, few studies focus on how planets form around low-mass stars. We aim to better characterize the formation process of M dwarf planets to fill this gap and aid in the interpretation of TESS results. We use ten sets of N-body planet formation simulations which vary in whether a gas disc is present, initial range of embry…
▽ More
NASA's TESS mission is expected to discover hundreds of M dwarf planets. However, few studies focus on how planets form around low-mass stars. We aim to better characterize the formation process of M dwarf planets to fill this gap and aid in the interpretation of TESS results. We use ten sets of N-body planet formation simulations which vary in whether a gas disc is present, initial range of embryo semi-major axes, and initial solid surface density profile. Each simulation begins with 147 equal-mass embryos around a 0.2 solar mass star and runs for 100 Myr. We find that planets form rapidly, with most collisions occurring within the first 1 Myr. The presence of a gas disc reduces the final number of planets relative to a gas-free environment and causes planets to migrate inward. We find that roughly a quarter of planetary systems experience their final giant impact inside the gas disc, suggesting that some super-Earths may be able to reaccrete an extended gaseous envelope after their final giant impact, though these may be affected by additional process such as photoevaporation. In addition, we find that the final distribution of planets does not retain a memory of the slope of the initial surface density profile, regardless of whether or not a gas disc is present. Thus, our results suggest that present-day observations are unlikely to provide sufficient information to accurately reverse-engineer the initial distribution of solids.
△ Less
Submitted 1 March, 2021;
originally announced March 2021.
-
Protoplanetary Disk Rings as Sites for Planetesimal Formation
Authors:
Daniel Carrera,
Jacob B. Simon,
Rixin Li,
Katherine A. Kretke,
Hubert Klahr
Abstract:
Axisymmetric dust rings are a ubiquitous feature of young protoplanetary disks. These rings are likely caused by pressure bumps in the gas profile; a small bump can induce a traffic jam-like pattern in the dust density, while a large bump may halt radial dust drift entirely. The resulting increase in dust concentration may trigger planetesimal formation by the streaming instability (SI), as the SI…
▽ More
Axisymmetric dust rings are a ubiquitous feature of young protoplanetary disks. These rings are likely caused by pressure bumps in the gas profile; a small bump can induce a traffic jam-like pattern in the dust density, while a large bump may halt radial dust drift entirely. The resulting increase in dust concentration may trigger planetesimal formation by the streaming instability (SI), as the SI itself requires some initial concentration. Here we present the first 3D simulations of planetesimal formation in the presence of a pressure bump modeled specifically after those observed by ALMA. In particular, we place a pressure bump at the center of a large 3D shearing box, along with an initial solid-to-gas ratio of $Z = 0.01$, and we include both particle back-reaction and particle self-gravity. We consider both mm-sized and cm-sized particles separately. For simulations with cm-sized particles, we find that even a small pressure bump leads to the formation of planetesimals via the streaming instability; a pressure bump does {\it not} need to fully halt radial particle drift for the SI to become efficient. Furthermore, pure gravitational collapse via concentration in pressure bumps (such as would occur at sufficiently high concentrations and without the streaming instability) is not responsible for planetesimal formation. For mm-sized particles, we find tentative evidence that planetesimal formation does not occur. This result, if it holds up at higher resolution and for a broader range of parameters, could put strong constraints on where in protoplanetary disks planetesimals can form. Ultimately, however, our results suggest that for cm-sized particles, planetesimal formation in pressure bumps is an extremely robust process.
△ Less
Submitted 12 January, 2021; v1 submitted 4 August, 2020;
originally announced August 2020.
-
Architectures of Exoplanetary Systems. III: Eccentricity and Mutual Inclination Distributions of AMD-stable Planetary Systems
Authors:
Matthias Y. He,
Eric B. Ford,
Darin Ragozzine,
Daniel Carrera
Abstract:
The angular momentum deficit (AMD) of a planetary system is a measure of its orbital excitation and a predictor of long-term stability. We adopt the AMD-stability criteria to constrain the orbital architectures for exoplanetary systems. Previously, He, Ford, & Ragozzine (2019) (arXiv:1907.07773v2) showed through forward modelling (SysSim) that the observed multiplicity distribution can be well rep…
▽ More
The angular momentum deficit (AMD) of a planetary system is a measure of its orbital excitation and a predictor of long-term stability. We adopt the AMD-stability criteria to constrain the orbital architectures for exoplanetary systems. Previously, He, Ford, & Ragozzine (2019) (arXiv:1907.07773v2) showed through forward modelling (SysSim) that the observed multiplicity distribution can be well reproduced by two populations consisting of a low and a high mutual inclination component. Here, we show that a broad distribution of mutual inclinations arising from systems at the AMD-stability limit can also match the observed Kepler population. We show that distributing a planetary system's maximum AMD amongst its planets results in a multiplicity-dependent distribution of eccentricities and mutual inclinations. Systems with intrinsically more planets have lower median eccentricities and mutual inclinations, and this trend is well described by power-law functions of the intrinsic planet multiplicity ($n$): $\tildeμ_{e,n} \propto n^{-1.74_{-0.07}^{+0.11}}$ and $\tildeμ_{i,n} \propto n^{-1.73_{-0.08}^{+0.09}}$, where $\tildeμ_{e,n}$ and $\tildeμ_{i,n}$ are the medians of the eccentricity and inclination distributions. We also find that intrinsic single planets have higher eccentricities ($σ_{e,1} \sim 0.25$) than multi-planet systems, and that the trends with multiplicity appear in the observed distributions of period-normalized transit duration ratios. We show that the observed preferences for planet size orderings and uniform spacings are more extreme than what can be produced by the detection biases of the Kepler mission alone. Finally, we find that for systems with detected transiting planets between 5 and 10 days, there is another planet with a greater radial velocity signal $\simeq~53\%$ of the time.
△ Less
Submitted 9 December, 2020; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Increasing the Discovery Space in Astrophysics - A Collation of Six Submitted White Papers
Authors:
G. Fabbiano,
M. Elvis,
A. Accomazzi,
G. B. Berriman,
N. Brickhouse,
S. Bose,
D. Carrera,
I. Chilingarian,
F. Civano,
B. Czerny,
R. D'Abrusco,
B. Diemer,
J. Drake,
R. Emami Meibody,
J. R. Farah,
G. G. Fazio,
E. Feigelson,
F. Fornasini,
Jay Gallagher,
J. Grindlay,
L. Hernquist,
D. J. James,
M. Karovska,
V. Kashyap,
D. -W. Kim
, et al. (24 additional authors not shown)
Abstract:
We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the kn…
▽ More
We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.
△ Less
Submitted 18 March, 2019; v1 submitted 15 March, 2019;
originally announced March 2019.
-
Planet formation: The case for large efforts on the computational side
Authors:
Wladimir Lyra,
Thomas Haworth,
Bertram Bitsch,
Simon Casassus,
Nicolás Cuello,
Thayne Currie,
Andras Gáspár,
Hannah Jang-Condell,
Hubert Klahr,
Nathan Leigh,
Giuseppe Lodato,
Mordecai-Mark Mac Low,
Sarah Maddison,
George Mamatsashvili,
Colin McNally,
Andrea Isella,
Sebastián Pérez,
Luca Ricci,
Debanjan Sengupta,
Dimitris Stamatellos,
Judit Szulágyi,
Richard Teague,
Neal Turner,
Orkan Umurhan,
Jacob White
, et al. (32 additional authors not shown)
Abstract:
Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been a…
▽ More
Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.
△ Less
Submitted 11 March, 2019;
originally announced March 2019.
-
Planet-planet scattering as the source of the highest eccentricity exoplanets
Authors:
Daniel Carrera,
Sean R. Raymond,
Melvyn B. Davies
Abstract:
Most giant exoplanets discovered by radial velocity surveys have much higher eccentricities than those in the solar system. The planet--planet scattering mechanism has been shown to match the broad eccentricity distribution, but the highest-eccentricity planets are often attributed to Kozai-Lidov oscillations induced by a stellar companion. Here we investigate whether the highly eccentric exoplane…
▽ More
Most giant exoplanets discovered by radial velocity surveys have much higher eccentricities than those in the solar system. The planet--planet scattering mechanism has been shown to match the broad eccentricity distribution, but the highest-eccentricity planets are often attributed to Kozai-Lidov oscillations induced by a stellar companion. Here we investigate whether the highly eccentric exoplanet population can be produced entirely by scattering. We ran 500 N-body simulations of closely packed giant-planet systems that became unstable under their own mutual perturbations. We find that the surviving bound planets can have eccentricities up to $e > 0.99$, with a maximum of 0.999017 in our simulations. This suggests that there is no maximum eccentricity that can be produced by planet--planet scattering. Importantly, we find that extreme eccentricities are not extremely rare; the eccentricity distribution for all giant exoplanets with $e > 0.3$ is consistent with all planets concerned being generated by scattering. Our results show that the discovery of planets with extremely high eccentricities does not necessarily signal the action of the Kozai-Lidov mechanism.
△ Less
Submitted 9 September, 2019; v1 submitted 6 March, 2019;
originally announced March 2019.
-
Formation of short-period planets by disk migration
Authors:
Daniel Carrera,
Eric B. Ford,
Andre Izidoro
Abstract:
Protoplanetary disks are thought to be truncated at orbital periods of around 10 days. Therefore, origin of rocky short period planets with $P < 10$ days is a puzzle. We propose that many of these planets may form through the Type-I migration of planets locked into a chain of mutual mean motion resonances. We ran N-body simulations of planetary embryos embedded in a protoplanetary disk. The embryo…
▽ More
Protoplanetary disks are thought to be truncated at orbital periods of around 10 days. Therefore, origin of rocky short period planets with $P < 10$ days is a puzzle. We propose that many of these planets may form through the Type-I migration of planets locked into a chain of mutual mean motion resonances. We ran N-body simulations of planetary embryos embedded in a protoplanetary disk. The embryos experienced gravitational scatterings, collisions, disk torques, and dampening of orbital eccentricity and inclination. We then modelled Kepler observations of these planets using a forward model of both the transit probability and the detection efficiency of the Kepler pipeline. We found that planets become locked into long chains of mean motion resonances that migrate in unison. When the chain reaches the edge of the disk, the inner planets are pushed past the edge due to the disk torques acting on the planets farther out in the chain. Our simulated systems successfully reproduce the observed period distribution of short period Kepler planets between 1 and 2 $R_\oplus$. However, we obtain fewer closely packed short period planets than in the Kepler sample. Our results provide valuable insight into the planet formation process, and suggests that resonance locks, migration, and dynamical instabilities play important roles the the formation and evolution of close-in small exoplanets.
△ Less
Submitted 5 March, 2019;
originally announced March 2019.
-
Identifying inflated super-Earths and photo-evaporated cores
Authors:
Daniel Carrera,
Eric B. Ford,
Andre Izidoro,
Daniel Jontof-Hutter,
Sean N. Raymond,
Angie Wolfgang
Abstract:
We present empirical evidence, supported by a planet formation model, to show that the curve $R/R_\oplus = 1.05\,(F/F_\oplus)^{0.11}$ approximates the location of the so-called photo-evaporation valley. Planets below that curve are likely to have experienced complete photo-evaporation, and planets just above it appear to have inflated radii; thus we identify a new population of inflated super-Eart…
▽ More
We present empirical evidence, supported by a planet formation model, to show that the curve $R/R_\oplus = 1.05\,(F/F_\oplus)^{0.11}$ approximates the location of the so-called photo-evaporation valley. Planets below that curve are likely to have experienced complete photo-evaporation, and planets just above it appear to have inflated radii; thus we identify a new population of inflated super-Earths and mini-Neptunes. Our N-body simulations are set within an evolving protoplanetary disk and include prescriptions for orbital migration, gas accretion, and atmospheric loss due to giant impacts. Our simulated systems broadly match the sizes and periods of super-Earths in the Kepler catalog. They also reproduce the relative sizes of adjacent planets in the same system, with the exception of planet pairs that straddle the photo-evaporation valley. This latter group is populated by planet pairs with either very large or very small size ratios ($R_{\rm out} / R_{\rm in} \gg 1$ or $R_{\rm out} / R_{\rm in} \ll 1$) and a dearth of size ratios near unity. It appears that this feature could be reproduced if the planet outside the photo-evaporation valley (typically the outer planet, but some times not) has its atmosphere significantly expanded by stellar irradiation. This new population of planets may be ideal targets for future transit spectroscopy observations with the upcoming James Webb Space Telescope.
△ Less
Submitted 5 September, 2018; v1 submitted 13 April, 2018;
originally announced April 2018.
-
Toward an initial mass function for giant planets
Authors:
Daniel Carrera,
Melvyn B. Davies,
Anders Johansen
Abstract:
The distribution of exoplanet masses is not primordial. After the initial stage of planet formation is complete, the gravitational interactions between planets can lead to the physical collision of two planets, or the ejection of one or more planets from the system. When this occurs, the remaining planets are typically left in more eccentric orbits. Here we use present-day eccentricities of the ob…
▽ More
The distribution of exoplanet masses is not primordial. After the initial stage of planet formation is complete, the gravitational interactions between planets can lead to the physical collision of two planets, or the ejection of one or more planets from the system. When this occurs, the remaining planets are typically left in more eccentric orbits. Here we use present-day eccentricities of the observed exoplanet population to reconstruct the initial mass function of exoplanets before the onset of dynamical instability. We developed a Bayesian framework that combines data from N-body simulations with present-day observations to compute a probability distribution for the planets that were ejected or collided in the past. Integrating across the exoplanet population, we obtained an estimate of the initial mass function of exoplanets. We find that the ejected planets are primarily sub-Saturn type planets. While the present-day distribution appears to be bimodal, with peaks around $\sim 1 M_{\rm J}$ and $\sim 20 M_\oplus$, this bimodality does not seem to be primordial. Instead, planets around $\sim 60 M_\oplus$ appear to be preferentially removed by dynamical instabilities. Attempts to reproduce exoplanet populations using population synthesis codes should be mindful of the fact that the present population has been depleted of intermediate-mass planets. Future work should explore how the system architecture and multiplicity might alter our results.
△ Less
Submitted 25 March, 2017;
originally announced March 2017.
-
Planetesimal formation by the streaming instability in a photoevaporating disk
Authors:
Daniel Carrera,
Uma Gorti,
Anders Johansen,
Melvyn B. Davies
Abstract:
Recent years have seen growing interest in the streaming instability as a candidate mechanism to produce planetesimals. However, these investigations have been limited to small-scale simulations. We now present the results of a global protoplanetary disk evolution model that incorporates planetesimal formation by the streaming instability, along with viscous accretion, photoevaporation by EUV, FUV…
▽ More
Recent years have seen growing interest in the streaming instability as a candidate mechanism to produce planetesimals. However, these investigations have been limited to small-scale simulations. We now present the results of a global protoplanetary disk evolution model that incorporates planetesimal formation by the streaming instability, along with viscous accretion, photoevaporation by EUV, FUV, and X-ray photons, dust evolution, the water ice line, and stratified turbulence. Our simulations produce massive (60-130 $M_\oplus$) planetesimal belts beyond 100 au and up to $\sim 20 M_\oplus$ of planetesimals in the middle regions (3-100 au). Our most comprehensive model forms 8 $M_\oplus$ of planetesimals inside 3 au, where they can give rise to terrestrial planets. The planetesimal mass formed in the inner disk depends critically on the timing of the formation of an inner cavity in the disk by high-energy photons. Our results show that the combination of photoevaporation and the streaming instability are efficient at converting the solid component of protoplanetary disks into planetesimals. Our model, however, does not form enough early planetesimals in the inner and middle regions of the disk to give rise to giant planets and super-Earths with gaseous envelopes. Additional processes such as particle pileups and mass loss driven by MHD winds may be needed to drive the formation of early planetesimal generations in the planet forming regions of protoplanetary disks.
△ Less
Submitted 22 March, 2017;
originally announced March 2017.
-
Concentrating small particles in protoplanetary disks through the streaming instability
Authors:
Chao-Chin Yang,
Anders Johansen,
Daniel Carrera
Abstract:
Laboratory experiments indicate that direct growth of silicate grains via mutual collisions can only produce particles up to roughly millimeters in size. On the other hand, recent simulations of the streaming instability have shown that mm/cm-sized particles require an excessively high metallicity for dense filaments to emerge. Using a numerical algorithm for stiff mutual drag force, we perform si…
▽ More
Laboratory experiments indicate that direct growth of silicate grains via mutual collisions can only produce particles up to roughly millimeters in size. On the other hand, recent simulations of the streaming instability have shown that mm/cm-sized particles require an excessively high metallicity for dense filaments to emerge. Using a numerical algorithm for stiff mutual drag force, we perform simulations of small particles with significantly higher resolutions and longer simulation times than in previous investigations. We find that particles of dimensionless stopping time $τ_\mathrm{s} = 10^{-2}$ and $10^{-3}$ -- representing mm- and cm-sized particles interior of the water ice line -- concentrate themselves via the streaming instability at a solid abundance of a few percent. We thus revise a previously published critical solid abundance curve for the regime of $τ_\mathrm{s} \ll 1$. The solid density in the concentrated regions reaches values higher than the Roche density, indicating that direct collapse of particles down to mm sizes into planetesimals is possible. Our results hence bridge the gap in particle size between direct dust growth limited by bouncing and the streaming instability.
△ Less
Submitted 29 June, 2017; v1 submitted 21 November, 2016;
originally announced November 2016.
-
Survival of habitable planets in unstable planetary systems
Authors:
Daniel Carrera,
Melvyn B. Davies,
Anders Johansen
Abstract:
Many observed giant planets lie on eccentric orbits. Such orbits could be the result of strong scatterings with other giant planets. The same dynamical instability that produces these scatterings may also cause habitable planets in interior orbits to become ejected, destroyed, or be transported out of the habitable zone. We say that a habitable planet has resilient habitability if it is able to av…
▽ More
Many observed giant planets lie on eccentric orbits. Such orbits could be the result of strong scatterings with other giant planets. The same dynamical instability that produces these scatterings may also cause habitable planets in interior orbits to become ejected, destroyed, or be transported out of the habitable zone. We say that a habitable planet has resilient habitability if it is able to avoid ejections and collisions and its orbit remains inside the habitable zone. Here we model the orbital evolution of rocky planets in planetary systems where giant planets become dynamically unstable. We measure the resilience of habitable planets as a function of the observed, present-day masses and orbits of the giant planets. We find that the survival rate of habitable planets depends strongly on the giant planet architecture. Equal-mass planetary systems are far more destructive than systems with giant planets of unequal masses. We also establish a link with observation; we find that giant planets with present-day eccentricities higher than 0.4 almost never have a habitable interior planet. For a giant planet with an present-day eccentricity of 0.2 and semimajor axis of 5 AU orbiting a Sun-like star, 50% of the orbits in the habitable zone are resilient to the instability. As semimajor axis increases and eccentricity decreases, a higher fraction of habitable planets survive and remain habitable. However, if the habitable planet has rocky siblings, there is a significant risk of rocky planet collisions that would sterilize the planet.
△ Less
Submitted 29 August, 2016; v1 submitted 4 May, 2016;
originally announced May 2016.
-
How to form asteroids from mm-sized grains
Authors:
Daniel Carrera,
Anders Johansen,
Melvyn B. Davies
Abstract:
The size distribution of asteroids in the solar system suggests that they formed top-down, with 100-1000 km bodies forming from the gravitational collapse of dense clumps of small solid particles. We investigate the conditions under which solid particles can form dense clumps in a protoplanetary disc. We used a hydrodynamic code to model the solid-gas interaction in disc. We found that particles d…
▽ More
The size distribution of asteroids in the solar system suggests that they formed top-down, with 100-1000 km bodies forming from the gravitational collapse of dense clumps of small solid particles. We investigate the conditions under which solid particles can form dense clumps in a protoplanetary disc. We used a hydrodynamic code to model the solid-gas interaction in disc. We found that particles down to millimeter size can form dense clumps, but only in regions where solids make $\sim$ 8% of the local surface density. More generally, we mapped the range of particle sizes and concentrations that is consistent with the formation of particle clumps.
△ Less
Submitted 7 April, 2016;
originally announced April 2016.
-
How to form planetesimals from mm-sized chondrules and chondrule aggregates
Authors:
Daniel Carrera,
Anders Johansen,
Melvyn B. Davies
Abstract:
The size distribution of asteroids and Kuiper belt objects in the solar system is difficult to reconcile with a bottom-up formation scenario due to the observed scarcity of objects smaller than $\sim$100 km in size. Instead, planetesimals appear to form top-down, with large $100-1000$ km bodies forming from the rapid gravitational collapse of dense clumps of small solid particles. In this paper we…
▽ More
The size distribution of asteroids and Kuiper belt objects in the solar system is difficult to reconcile with a bottom-up formation scenario due to the observed scarcity of objects smaller than $\sim$100 km in size. Instead, planetesimals appear to form top-down, with large $100-1000$ km bodies forming from the rapid gravitational collapse of dense clumps of small solid particles. In this paper we investigate the conditions under which solid particles can form dense clumps in a protoplanetary disk. We use a hydrodynamic code to model the interaction between solid particles and the gas inside a shearing box inside the disk, considering particle sizes from sub-millimeter-sized chondrules to meter-sized rocks. We find that particles down to millimeter sizes can form dense particle clouds through the run-away convergence of radial drift known as the streaming instability. We make a map of the range of conditions (strength of turbulence, particle mass-loading, disk mass, and distance to the star) which are prone to producing dense particle clumps. Finally, we estimate the distribution of collision speeds between mm-sized particles. We calculate the rate of sticking collisions and obtain a robust upper limit on the particle growth timescale of $\sim$$10^5$ years. This means that mm-sized chondrule aggregates can grow on a timescale much smaller than the disk accretion timescale ($\sim$$10^6 - 10^7$ years). Our results suggest a pathway from the mm-sized grains found in primitive meteorites to fully formed asteroids. We speculate that asteroids may form from a positive feedback loop in which coagualation leads to particle clumping driven by the streaming instability. This clumping, in turn reduces collision speeds and enhances coagulation.} Future simulations should model coagulation and the streaming instability together to explore this feedback loop further.
△ Less
Submitted 24 April, 2015; v1 submitted 21 January, 2015;
originally announced January 2015.