Gap-opening Planets Research Articles

Observational diagnostics of elongated planet-induced vortices with realistic planet formation time-scales

Gap-opening planets can generate dust-trapping vortices that may explain some of the latest discoveries of high-contrast crescent-shaped dust asymmetries in transition discs. While planet-induced vortices were previously thought to have concentrated shapes, recent computational work has shown that these features naturally become much more elongated in the gas when simulations account for the relatively long timescale over which planets accrete their mass. In this work, we conduct two-fluid hydrodynamical simulations of vortices induced by slowly-growing Jupiter-mass planets in discs with very low viscosity ($\alpha = 3 \times 10^{-5}$). We simulate the dust dynamics for four particle sizes spanning 0.3 mm to 1 cm in order to produce synthetic ALMA images. In our simulations, we find that an elongated vortex still traps dust, but not directly at its center. With a flatter pressure bump and disruptions from the planet's overlapping spiral density waves, the dust instead circulates around the vortex. This motion (1) typically carries the peak off-center, (2) spreads the dust out over a wider azimuthal extent $\geq 180^{\circ}$, (3) skews the azimuthal profile towards the front of the vortex, and (4) can also create double peaks in newly-formed vortices. In particular, we expect that the most defining observational signature, a peak offset of more than $30^{\circ}$, should be detectable $>30\%$ of the time in observations with a beam diameter of at most the planet's separation from its star.

Slowing Down Type II Migration of Gas Giants to Match Observational Data

Abstract The mass and semimajor axis distribution of gas giants in exoplanetary systems obtained by radial velocity surveys shows that super-Jupiter-mass planets are piled up at ≳1 au, while Jupiter/sub-Jupiter-mass planets are broadly distributed from ∼0.03 au to beyond 1 au. This feature has not been explained by theoretical predictions. In order to reconcile this inconsistency, we investigate evolution of gas giants with a new type II migration formula by Kanagawa et al., by comparing the migration, growth timescales of gas giants, and disk lifetime, and by population synthesis simulation. While the classical migration model assumes that a gas giant opens up a clear gap in the protoplanetary disk and the planet migration is tied to the disk gas accretion, recent high-resolution simulations show that the migration of gap-opening planets is decoupled from the disk gas accretion and Kanagawa et al. proposed that type II migration speed is nothing other than type I migration speed with the reduced disk gas surface density in the gap. We show that with this new formula, type II migration is significantly reduced for super-Jupiter-mass planets, if the disk accretion is driven by the disk wind as suggested by recent magnetohydrodynamic simulations. Population synthesis simulations show that super-Jupiter-mass planets remain at ≳1 au without any additional ingredient such as disk photoevaporation. Therefore, the mystery of the pile-up of gas giants at ≳1 au will be theoretically solved if the new formula is confirmed and wind-driven disk accretion dominates.

Radial Migration of Gap-opening Planets in Protoplanetary Disks. I. The Case of a Single Planet

Abstract A large planet orbiting a star in a protoplanetary disk opens a density gap along its orbit due to the strong disk–planet interaction and migrates with the gap in the disk. It is expected that in the ideal case, a gap-opening planet migrates at the viscous drift speed, which is referred to as type II migration. However, recent hydrodynamic simulations have shown that, in general, the gap-opening planet is not locked to the viscous disk evolution. A new physical model is required to explain the migration speed of such a planet. For this reason, we re-examined the migration of a planet in the disk, by carrying out the two-dimensional hydrodynamic simulations in a wide parameter range. We have found that the torque exerted on the gap-opening planet depends on the surface density at the bottom of the gap. The planet migration slows down as the surface density of the bottom of the gap decreases. Using the gap model developed in our previous studies, we have constructed an empirical formula of the migration speed of the gap-opening planets, which is consistent with the results given by the hydrodynamic simulations performed by us and other researchers. Our model easily explains why the migration speed of the gap-opening planets can be faster than the viscous gas drift speed. It can also predict the planet mass at which the type I migration is no longer adequate due to the gap development in the disk, providing a gap formation criterion based on planetary migration.

Characterizing the Variable Dust Permeability of Planet-induced Gaps

Abstract Aerodynamic theory predicts that dust grains in protoplanetary disks will drift radially inward on comparatively short timescales. In this context, it has long been known that the presence of a gap opened by a planet can significantly alter the dust dynamics. In this paper, we carry out a systematic study employing long-term numerical simulations aimed at characterizing the critical particle size for retention outside a gap as a function of particle size, as well as various key parameters defining the protoplanetary disk model. To this end, we perform multifluid hydrodynamical simulations in two dimensions, including different dust species, which we treat as pressureless fluids. We initialize the dust outside of the planet’s orbit and study under which conditions dust grains are able to cross the gap carved by the planet. In agreement with previous work, we find that the permeability of the gap depends both on dust dynamical properties and the gas disk structure: while small dust follows the viscously accreting gas through the gap, dust grains approaching a critical size are progressively filtered out. Moreover, we introduce and compute a depletion factor that enables us to quantify the way in which higher viscosity, smaller planet mass, or a more massive disk can shift this critical size to larger values. Our results indicate that gap-opening planets may act to deplete the inner reaches of protoplanetary disks of large dust grains—potentially limiting the accretion of solids onto forming terrestrial planets.

Secular evolution of eccentricity in protoplanetary discs with gap-opening planets

We explore the evolution of the eccentricity of an accretion disc perturbed by an embedded planet whose mass is sufficient to open a large gap in the disc. Various methods for representing the orbit-averaged motion of an eccentric disc are discussed. We characterize the linear instability which leads to the growth of eccentricity by means of hydrodynamical simulations. We numerically recover the known result that eccentricity growth in the disc is possible when the planet-to-star mass ratio exceeds 0.003. For mass ratios larger than this threshold, the precession rates and growth rates derived from simulations, as well as the shape of the eccentric mode, compare well with the predictions of a linear theory of eccentric discs. We study mechanisms by which the eccentricity growth eventually saturates into a non-linear regime.

What is the Mass of a Gap-opening Planet?

Abstract High-contrast imaging instruments such as GPI and SPHERE are discovering gap structures in protoplanetary disks at an ever faster pace. Some of these gaps may be opened by planets forming in the disks. In order to constrain planet formation models using disk observations, it is crucial to find a robust way to quantitatively back out the properties of the gap-opening planets, in particular their masses, from the observed gap properties, such as their depths and widths. Combining 2D and 3D hydrodynamics simulations with 3D radiative transfer simulations, we investigate the morphology of planet-opened gaps in near-infrared scattered-light images. Quantitatively, we obtain correlations that directly link intrinsic gap depths and widths in the gas surface density to observed depths and widths in images of disks at modest inclinations under finite angular resolution. Subsequently, the properties of the surface density gaps enable us to derive the disk scale height at the location of the gap h, and to constrain the quantity M p 2/α, where M p is the mass of the gap-opening planet and α characterizes the viscosity in the gap. As examples, we examine the gaps recently imaged by VLT/SPHERE, Gemini/GPI, and Subaru/HiCIAO in HD 97048, TW Hya, HD 169142, LkCa 15, and RX J1615.3-3255. Scale heights of the disks and possible masses of the gap-opening planets are derived assuming each gap is opened by a single planet. Assuming α = 10−3, the derived planet masses in all cases are roughly between 0.1 and 1 M J.

Slowly-growing gap-opening planets trigger weaker vortices

The presence of a giant planet in a low-viscosity disc can create a gap edge in the disc's radial density profile sharp enough to excite the Rossby Wave Instability. This instability may evolve into dust-trapping vortices that might explain the "banana-shaped" features in recently observed asymmetric transition discs with inner cavities. Previous hydrodynamical simulations of planet-induced vortices have neglected the timescale of hundreds to thousands of orbits to grow a massive planet to Jupiter-size. In this work, we study the effect of a giant planet's runaway growth timescale on the lifetime and characteristics of the resulting vortex. For two different planet masses (1 and 5 Jupiter masses) and two different disc viscosities ($\alpha$=3$\times 10^{-4}$ and 3$\times10^{-5}$), we compare the vortices induced by planets with several different growth timescales between 10 and 4000 planet orbits. In general, we find that slowly-growing planets create significantly weaker vortices with lifetimes and surface densities reduced by more than $50\%$. For the higher disc viscosity, the longest growth timescales in our study inhibit vortex formation altogether. Additionally, slowly-growing planets produce vortices that are up to twice as elongated, with azimuthal extents well above $180^{\circ}$ in some cases. These unique, elongated vortices likely create a distinct signature in the dust observations that differentiates them from the more concentrated vortices that correspond to planets with faster growth timescales. Lastly, we find that the low viscosities necessary for vortex formation likely prevent planets from growing quickly enough to trigger the instability in self-consistent models.

Stability of resonant configurations during the migration of planets and constraints on disk-planet interactions

We study the stability of mean-motion resonances (MMR) between two planets during their migration in a protoplanetary disk. We use an analytical model of resonances, and describe the effect of the disk by a migration timescale (T_{m,i}) and an eccentricity damping timescale (T_{e,i}) for each planet (i=1,2 respectively for the inner and outer planet). We show that the resonant configuration is stable if T_{e,1}/T_{e,2} > (e_1/e_2)^2. This general result can be used to put constraints on specific models of disk-planet interactions. For instance, using classical prescriptions for type I migration, we show that when the angular momentum deficit (AMD) of the inner orbit is larger than the outer's orbit AMD, resonant systems must have a locally inverted disk density profile to stay locked in resonance during the migration. This inversion is very untypical of type I migration and our criterion can thus provide an evidence against classical type I migration. That is indeed the case for the Jupiter-mass resonant systems HD 60532b, c (3:1 MMR), GJ 876b, c (2:1 MMR), and HD 45364b, c (3:2 MMR). This result may be an evidence for type II migration (gap opening planets), which is compatible with the large masses of these planets.

THE MIGRATION OF GAP-OPENING PLANETS IS NOT LOCKED TO VISCOUS DISK EVOLUTION

Most standard descriptions of Type II migration state that massive, gap-opening planets must migrate at the viscous drift rate. This is based on the idea that the disk is separated into an inner and outer region and gas is considered unable to cross the gap. In fact, gas easily crosses the gap on horseshoe orbits, nullifying this necessary premise which would set the migration rate. In this work, it is demonstrated using highly accurate numerical calculations that the actual migration rate is dependent on disk and planet parameters, and can be significantly larger or smaller than the viscous drift rate. In the limiting case of a disk much more massive than the secondary, the migration rate saturates to a constant which is sensitive to disk parameters and is not necessarily of order viscous rate. In the opposite limit of a low-mass disk, the migration rate decreases linearly with disk mass. Steady-state solutions in the low disk mass limit show no pile-up outside the secondary's orbit, and no corresponding drainage of the inner disk.