Simulations of the Fomalhaut System Within Its Local Galactic Environment

Fomalhaut A is among the most well-studied nearby stars and has been discovered to possess a putative planetary object as well as a remarkable eccentric dust belt. This eccentric dust belt has often been interpreted as the dynamical signature of one or more planets that elude direct detection. However, the system also contains two other stellar companions residing ~100,000 AU from Fomalhaut A. We have designed a new symplectic integration algorithm to model the evolution of Fomalhaut A's planetary dust belt in concert with the dynamical evolution of its stellar companions to determine if these companions are likely to have generated the dust belt's morphology. Using our numerical simulations, we find that close encounters between Fomalhaut A and B are expected, with a ~25% probability that the two stars have passed within at least 400 AU of each other at some point. Although the outcomes of such encounter histories are extremely varied, these close encounters nearly always excite the eccentricity of Fomalhaut A's dust belt and occasionally yield morphologies very similar to the observed belt. With these results, we argue that close encounters with Fomalhaut A's stellar companions should be considered a plausible mechanism to explain its eccentric belt, especially in the absence of detected planets capable of sculpting the belt's morphology. More broadly, we can also conclude from this work that very wide binary stars may often generate asymmetries in the stellar debris disks they host.Comment: Accepted to MNRAS, 22 pages, 15 figures, 2 appendice

Formation of short-period planets by disk migration

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.Comment: 12 pages, 11 figures, submitted to MNRA

Oort cloud (exo)planets

Dynamical instabilities among giant planets are thought to be nearly ubiquitous, and culminate in the ejection of one or more planets into interstellar space. Here we perform N-body simulations of dynamical instabilities while accounting for torques from the galactic tidal field. We find that a fraction of planets that would otherwise have been ejected are instead trapped on very wide orbits analogous to those of Oort cloud comets. The fraction of ejected planets that are trapped ranges from 1-10%, depending on the initial planetary mass distribution. The local galactic density has a modest effect on the trapping efficiency and the orbital radii of trapped planets. The majority of Oort cloud planets survive for Gyr timescales. Taking into account the demographics of exoplanets, we estimate that one in every 200-3000 stars could host an Oort cloud planet. This value is likely an overestimate, as we do not account for instabilities that take place at early enough times to be affected by their host stars' birth cluster, or planet stripping from passing stars. If the Solar System's dynamical instability happened after birth cluster dissolution, there is a ~7% chance that an ice giant was captured in the Sun's Oort cloud.Comment: MNRAS Letters, in press. Blog post about paper at https://planetplanet.net/2023/06/21/oort-cloud-exoplanets

Implantation of asteroids from the terrestrial planet region: The effect of the timing of the giant planet instability

The dynamical architecture and compositional diversity of the asteroid belt strongly constrain planet formation models. Recent Solar System formation models have shown that the asteroid belt may have been born empty and later filled with objects from the inner ($5 au) of the solar system. In this work, we focus on the implantation of inner solar system planetesimals into the asteroid belt - envisioned to represent S and/or E- type asteroids - during the late-stage accretion of the terrestrial planets. It is widely accepted that the solar system's giant planets formed in a more compact orbital configuration and evolved to their current dynamical state due to a planetary dynamical instability. In this work, we explore how the implantation efficiency of asteroids from the terrestrial region correlates with the timing of the giant planet instability, which has proven challenging to constrain. We carried out a suite of numerical simulations of the accretion of terrestrial planets considering different initial distributions of planetesimals in the terrestrial region and dynamical instability times. Our simulations show that a giant planet dynamical instability occurring at$t\gtrapprox5$ Myr -- relative to the time of the sun's natal disk dispersal -- is broadly consistent with the current asteroid belt, allowing the total mass carried out by S-complex type asteroids to be implanted into the belt from the terrestrial region. Finally, we conclude that an instability that occurs coincident with the gas disk dispersal is either inconsistent with the empty asteroid belt scenario, or may require that the gas disk in the inner solar system have dissipated at least a few Myr earlier than the gas in the outer disk (beyond Jupiter's orbit).Comment: Under review in Icaru

Formation of Super-Earths

Super-Earths are the most abundant planets known to date and are characterized by having sizes between that of Earth and Neptune, typical orbital periods of less than 100 days and gaseous envelopes that are often massive enough to significantly contribute to the planet's overall radius. Furthermore, super-Earths regularly appear in tightly-packed multiple-planet systems, but resonant configurations in such systems are rare. This chapters summarizes current super-Earth formation theories. It starts from the formation of rocky cores and subsequent accretion of gaseous envelopes. We follow the thermal evolution of newly formed super-Earths and discuss their atmospheric mass loss due to disk dispersal, photoevaporation, core-cooling and collisions. We conclude with a comparison of observations and theoretical predictions, highlighting that even super-Earths that appear as barren rocky cores today likely formed with primordial hydrogen and helium envelopes and discuss some paths forward for the future.Comment: Invited review accepted for publication in the 'Handbook of Exoplanets,' Planet Formation section, Springer Reference Works, Juan Antonio Belmonte and Hans Deeg, Ed

The link between Athor and EL meteorites does not constrain the timing of the giant planet instability

The asteroid Athor, residing today in the inner main asteroid belt, has been recently associated as the source of EL enstatite meteorites to Earth. It has been argued that Athor formed in the terrestrial region -- as indicated by similarity in isotopic compositions between Earth and EL meteorites -- and was implanted in the belt $\gtrsim$60 Myr after the formation of the solar system. A recently published study modelling Athor's implantation in the belt (Avdellidou et al 2024) further concluded, using an idealized set of numerical simulations, that Athor cannot have been scattered from the terrestrial region and implanted at its current location unless the giant planet dynamical instability occurred {\em after} Athor's implantation ($\gtrsim$60~Myr). In this work, we revisit this problem with a comprehensive suite of dynamical simulations of the implantation of asteroids into the belt during the terrestrial planet accretion. We find that Athor-like objects can in fact be implanted into the belt long after the giant planets' dynamical instability. The probability of implanting Athor analogs when the instability occurs at $\lesssim15$~Myr is at most a factor of $\sim$2 lower than that of an instability occurring at $\sim100$~Myr after the solar system formation. Moreover, Athor's implantation can occur up to $\gtrsim$100 Myr after the giant planet instability. We conclude that Athor's link to EL meteorites does not constrain the timing of the solar system's dynamical instability.Comment: Submitted. Comments welcom

The eccentricity distribution of giant planets and their relation to super-Earths in the pebble accretion scenario

Observations of the population of cold Jupiter planets ($r>$1 AU) show that nearly all of these planets orbit their host star on eccentric orbits. For planets up to a few Jupiter masses, eccentric orbits are thought to be the outcome of planet-planet scattering events taking place after gas dispersal. We simulate the growth of planets via pebble and gas accretion as well as the migration of multiple planetary embryos in their gas disc. We then follow the long-term dynamical evolution of our formed planetary system up to 100 Myr after gas disc dispersal. We investigate the importance of the initial number of protoplanetary embryos and different damping rates of eccentricity and inclination during the gas phase for the final configuration of our planetary systems. We constrain our model by comparing the final dynamical structure of our simulated planetary systems to that of observed exoplanet systems. Our results show that the initial number of planetary embryos has only a minor impact on the final orbital eccentricity distribution of the giant planets, as long as damping of eccentricity and inclination is efficient. If damping is inefficient (slow), systems with a larger initial number of embryos harbor larger average eccentricities. In addition, for slow damping rates, we observe that scattering events already during the gas disc phase are common and that the giant planets formed in these simulations match the observed giant planet eccentricity distribution best. These simulations also show that massive giant planets (above Jupiter mass) on eccentric orbits are less likely to host inner super-Earths as these get lost during the scattering phase, while systems with less massive giant planets on nearly circular orbits should harbor systems of inner super-Earths. Finally, our simulations predict that giant planets are on average not single, but live in multi-planet systems.Comment: 25 pages, 18 figures, accepted by A&A, with language correction

Implications of Jupiter Inward Gas-Driven Migration for the Inner Solar System

The migration history of Jupiter in the sun's natal disk remains poorly constrained. Here we consider how Jupiter's migration affects small-body reservoirs and how this constrains its original orbital distance from the Sun. We study the implications of large-scale and inward radial migration of Jupiter for the inner solar system while considering the effects of collisional evolution of planetesimals. We use analytical prescriptions to simulate the growth and migration of Jupiter in the gas disk. We assume the existence of a planetesimal disk inside Jupiter's initial orbit. This planetesimal disk received an initial total mass and size-frequency distribution (SFD). Planetesimals feel the effects of aerodynamic gas drag and collide with one another, mostly while shepherded by the migrating Jupiter. Our main goal is to measure the amount of mass in planetesimals implanted into the main asteroid belt (MAB) and the SFD of the implanted population. We also monitor the amount of dust produced during planetesimal collisions. We find that the SFD of the planetesimal population implanted into the MAB tends to resemble that of the original planetesimal population interior to Jupiter. We also find that unless very little or no mass existed between 5 au and Jupiter's original orbit, it would be difficult to reconcile the current low mass of the MAB with the possibility that Jupiter migrated from distances beyond 15 au. This is because the fraction of the original disk mass that gets implanted into the MAB is very large. Finally, we discuss the implications of our results in terms of dust production to the so-called NC-CC isotopic dichotomy.Comment: Accepted for publication in The Astrophysical Journal Letters; In pres

Survival and dynamics of rings of co-orbital planets under perturbations

In co-orbital planetary systems, two or more planets share the same orbit around their star. Here we test the dynamical stability of co-orbital rings of planets perturbed by outside forces. We test two setups: i) 'stationary' rings of planets that, when unperturbed, remain equally-spaced along their orbit; and ii) horseshoe constellation systems, in which planets are continually undergoing horseshoe librations with their immediate neighbors. We show that a single rogue planet crossing the planets' orbit more massive than a few lunar masses (0.01-0.04 Earth masses) systematically disrupts a co-orbital ring of 6, 9, 18, or 42 Earth-mass planets located at 1 au. Stationary rings are more resistant to perturbations than horseshoe constellations, yet when perturbed they can transform into stable horseshoe constellation systems. Given sufficient time, any co-orbital ring system will be perturbed into either becoming a horseshoe constellation or complete destabilization.Comment: 5 pages, 4 figures. Re-submitted to MNRAS. Blog post about co-orbital constellations here: https://planetplanet.net/2023/04/20/constellations-of-co-orbital-planets

Constellations of co-orbital planets: horseshoe dynamics, long-term stability, transit timing variations, and potential as SETI beacons

Co-orbital systems contain two or more bodies sharing the same orbit around a planet or star. The best-known flavors of co-orbital systems are tadpoles (in which two bodies' angular separations oscillate about the L4/L5 Lagrange points $60^\circ$ apart) and horseshoes (with two bodies periodically exchanging orbital energy to trace out a horseshoe shape in a co-rotating frame). Here, we use N-body simulations to explore the parameter space of many-planet horseshoe systems. We show that up to 24 equal-mass, Earth-mass planets can share the same orbit at 1 au, following a complex pattern in which neighboring planets undergo horseshoe oscillations. We explore the dynamics of horseshoe constellations, and show that they can remain stable for billions of years and even persist through their stars' post-main sequence evolution. With sufficient observations, they can be identified through their large-amplitude, correlated transit timing variations. Given their longevity and exotic orbital architectures, horseshoe constellations may represent potential SETI beacons.Comment: 10 pages, 10 figures. Published in MNRAS. YouTube playlist with animations of horseshoe constellation systems here: https://www.youtube.com/playlist?list=PLelMZVM3ka3F335LGLxkxrD1ieiLJYQ5N . Blog post here: https://planetplanet.net/2023/04/20/constellations-of-co-orbital-planets