IC project: Exploring Planet-disk Interaction in Dusty Protoplanetary Disks and Their Observational Implications
Abstract not provided.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Abstract not provided.
Aspects of star formation are considered, taking into account those observations and models which have a relation to the formation of the solar system. Attention is given to the galactic structure, the stellar properties, planetary indications, interstellar clouds, observations of forming stars, the conditions for cloud collapse, models of star formation, the modeling of secondary features, and aspects of dynamical evolution. Aspects of cosmochemistry are also explored. Questions of chronology are discussed along with nuclide variations, chemical variations among chondritic meteorites, and chemical variations among differentiated objects. A description of planet formation is provided, taking into consideration, gas-dust interaction, sun and accretion disk dynamics, planetesimal swarms, planet-disk interactions, planetary systems, terrestrial planet formation, and gaseous protoplanet contraction.
Under the support of NASA Origins grant, we studied the formation of gaps in protoplanetary disks due the tidal interaction between a fully grown protoplanet and protostellar disk. The result of this study is published in the Astrophysical Journal, (vol 514, 344-367, 1999) and in several conference proceedings. The main focus of this work is to analyze planet-disk interaction during the final stages of protoplanetary formation.
The proposal achieved many of its objectives. The main area of investigation was the interaction of young planets with surrounding protostellar disks. The topics of interest include: 1) Simulations of Planet-Disk Interactions; 2) Secular Interactions Between Inclined Planets and a Gaseous Disk; 3) On the Tilting of Protostellar Disks by Resonant Tidal Effects; 4) Three-Dimensional Waves in Thermally Stratified Disks; and 5) Predictions of the Distribution of Planets. A list of publications resulting from this grant is also presented.
When a planet orbits inside a debris disk like the disk around Vega or Beta Pictoris, the planet may be invisible, but the patterns it creates in the disk may give it away. Observing and decoding these patterns may be the only way we can detect exo-Neptunes orbiting more than 20 AU from their stars, and the only way we can spot planets in systems undergoing the late stages of planet formation. Fortunately, every few months, a new image of a debris disk appears with curious structures begging for explanation. I'll describe some new ideas in the theory of these planet-disk interactions and provide a buyers guide to the latest models (and the planets they predict).
Following the optical imaging of exoplanet candidate Fomalhaut b (Fom b), we present a numerical model of how Fomalhaut's debris disk is gravitationally shaped by a single interior planet. The model is simple, adaptable to other debris disks, and can be extended to accommodate multiple planets. If Fom b is the dominant perturber of the belt, then to produce the observed disk morphology it must have a mass M(sub pl) < 3M(sub J), an orbital semimajor axis a(sub pl) > 101.5AU, and an orbital eccentricity e(sub pl) = 0.11 - 0.13. These conclusions are independent of Fom b's photometry. To not disrupt the disk, a greater mass for Fom b demands a smaller orbit farther removed from the disk; thus, future astrometric measurement of Fom b's orbit, combined with our model of planet-disk interaction, can be used to determine the mass more precisely. The inner edge of the debris disk at a approximately equals 133AU lies at the periphery of Fom b's chaotic zone, and the mean disk eccentricity of e approximately equals 0.11 is secularly forced by the planet, supporting predictions made prior to the discovery of Fom b. However, previous mass constraints based on disk morphology rely on several oversimplifications. We explain why our constraint is more reliable. It is based on a global model of the disk that is not restricted to the planet's chaotic zone boundary. Moreover, we screen disk parent bodies for dynamical stability over the system age of approximately 100 Myr, and model them separately from their dust grain progeny; the latter's orbits are strongly affected by radiation pressure and their lifetimes are limited to approximately 0.1 Myr by destructive grain-grain collisions. The single planet model predicts that planet and disk orbits be apsidally aligned. Fomalhaut b's nominal space velocity does not bear this out, but the astrometric uncertainties are difficult to quantify. Even if the apsidal misalignment proves real, our calculated upper mass limit of 3 M(sub J) still holds. Parent bodies are evacuated from mean-motion resonances with Fom b; these empty resonances are akin to the Kirkwood gaps opened by Jupiter. The belt contains at least 3M(sub Earth) of solids that are grinding down to dust, their velocity dispersions stirred so strongly by Fom b that collisions are destructive. Such a large mass in solids is consistent with Fom b having formed in situ.
Jupiter’s core is generally assumed to have formed beyond the snow line. Here we consider an alternative scenario, that Jupiter’s core may have accumulated in the innermost parts of the protoplanetary disk. A growing body of research suggests that small particles (“pebbles”) continually drift inward through the disk. If a fraction of drifting pebbles is trapped at the inner edge of the disk a several Earth-mass core can quickly grow. Subsequently, the core may migrate outward beyond the snow line via planet-disk interactions. Of course, to reach the outer Solar System Jupiter’s core must traverse the terrestrial planet-forming region. We use N-body simulations including synthetic forces from an underlying gaseous disk to study how the outward migration of Jupiter’s core sculpts the terrestrial zone. If the outward migration is fast (τ(sub mig) ∼ 10(exp 4) years), the core simply migrates past resident planetesimals and planetary embryos. However, if its migration is slower (τ(sub mig) ∼ 10(exp 5) years) the core removes solids from the inner disk by shepherding objects in mean motion resonances. In many cases the disk interior to 0.5-1 AU is cleared of embryos and most planetesimals. By generating a mass deficit close to the Sun, the outward migration of Jupiter’s core may thus explain the absence of terrestrial planets closer than Mercury. Jupiter’s migrating core often stimulates the growth of another large (∼ Earth-mass) core – that may provide a seed for Saturn’s core – trapped in exterior resonance. The migrating core also may transport a fraction of terrestrial planetesimals, such as the putative parent bodies of iron meteorites, to the asteroid belt.
Context : Super-Earths can form at large orbital radii and migrate inward due to tidal interactions with the circumstellar disk. In this scenario, convergent migration may occur and lead to the formation of resonant pairs of planets. Aim : We explore the conditions under which convergent migration and resonance capture take place, and what dynamical consequences can be expected on the dust distribution surrounding the resonant pair. Methods : We combine hydrodynamic planet–disk interaction models with dust evolution calculations to investigate the signatures produced in the dust distribution by a pair of planets in mean-motion resonances. Results : We find that convergent migration takes place when the outer planet is the more massive. However, convergent migration also depends on the local properties of the disk, and divergent migration may result as well. For similar disk parameters, the capture in low degree resonances (e.g., 2:1 or 3:2) is preferred close to the star where the resonance strength can more easily overcome the tidal torques exerted by the gaseous disk. Farther away from the star, convergent migration may result in capture in high degree resonances. The dust distribution shows potentially observable features typically when the planets are trapped in a 2:1 resonance. In other cases, with higher degree resonances (e.g., 5:4 or 6:5) dust features may not be sufficiently pronounced to be easily observable. Conclusions : The degree of resonance established by a pair of super-Earths may be indicative of the location in the disk where capture occurred. There can be significant differences in the dust distribution around a single super-Earth and a pair of super-Earths in resonance.
Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.
Binary systems can be born surrounded by circumbinary discs. The gaseous disc around either of the two stellar companions can have its life extended by the supply of mass arriving from the circumbinary disc. The objective of this study is to investigate the gravitational interactions exerted by a compact and eccentric binary system on the circumbinary and circumprimary discs, and the resulting transport of gas and solids between the disc components. We assume that the gas in the system behaves as a fluid, and we model its evolution by means of high-resolution hydrodynamical simulations. Dust grains are modelled as Lagrangian particles that interact with the gas and the stars. Our models indicate that significant fluxes of gas and dust proceed from the circumbinary disc towards the circumprimary disc. For the applied system parameters, grains of certain sizes are segregated outside the tidal gap generated by the stars. Consequently, the size distribution of the transported dust is not continuous, but presents a gap in the millimetre size range. In close binaries, the lifetime of an isolated circumprimary disc is found to be short, ~10 5 years, because of its low mass. However, because of the influx of gas from beyond the tidal gap, the disc around the primary star can survive much longer, ~10 6 years, as long as gas accretion from the circumbinary disc continues. The supply of solids and the extended lifetime of a circumbinary disc also aids in the possible formation of giant planets. Compared to close binary systems without a circumbinary disc, we expect a higher frequency of singleplanet or multiple-planet systems. Additionally, a planetesimal or debris belt can form in the proximity of the truncation radius of the circumprimary disc and/or around the location of the exterior edge of the tidal gap.
Context. Giant planets are observed orbiting the primary stars of close binary systems. Such planets may have formed in compact circumprimary disks, which once surrounded these stars, under conditions much different than those encountered around single stars. Aims. In order to quantify the effects of the strong gravitational perturbations exerted on circumprimary disk material, the three-dimensional (3D) dynamics of gas and dust in orbit around the primary star of a compact and eccentric binary system was modeled by applying the stellar and orbital parameters of γ Cephei, a well-known system that can be representative of a class of close binaries. Methods. Circumprimary gas was approximated as an Eulerian viscous and compressible fluid and modeled by means of 3D hydrodynamical simulations, assuming locally isothermal conditions in the medium around the primary star. Dust grains were modeled as Lagrangean particles, subjected to gravity and aerodynamic drag forces. Models that include a giant planet were also considered. Results. Models indicate that spiral density waves excited around pericenter passage propagate toward the inner boundary of the disk, through at least a few pressure scale-heights from the mid-plane, inducing radial and vertical mixing in the gas. However, perturbations imparted to gas, both in terms of eccentricity and precession, are far weaker than previously estimated by two-dimensional (2D) simulations. Models predict small eccentricities, ≲0.03, and slow retrograde precession. The addition of a giant planet does not change the low eccentricity state of the disk. The parameters applied to the disk would lead to the formation of a massive planet, many times the mass of Jupiter, in agreement with some observations. Micron to mm-size dust grains are well coupled to the gas, resulting in similar dynamics and statistically similar distributions of orbital elements. The planet only affects the dust distributions locally. In agreement with outcomes of recent 2D models, the lifetime of an isolated circumprimary disk would be brief, ~10 5 years, because of its compact nature, requiring a long-term external supply of mass to allow for the in situ formation of a giant planet.
In this paper, we present H band spectroscopic and Hα photometric observations of HD 100546 obtained with the Gemini Planet Imager and the Magellan Visible AO camera. We detect H band emission at the location of the protoplanet HD 100546 b, but show that the choice of data processing parameters strongly affects the morphology of this source. It appears point-like in some aggressive reductions, but rejoins an extended disk structure in the majority of the others. Furthermore, we demonstrate that this emission appears stationary on a timescale of 4.6 years, inconsistent at the 2σ level with a Keplerian clockwise orbit at 59 au in the disk plane. The H band spectrum of the emission is inconsistent with any type of low effective temperature object or accreting protoplanetary disk. It strongly suggests a scattered-light origin, as this is consistent with the spectrum of the star and the spectra extracted at other locations in the disk. A non-detection at the 5σ level of HD 100546 b in differential Hα imaging places an upper limit, assuming the protoplanet lies in a gap free of extinction, on the accretion luminosity of 1.7 × 10 -4 L ⊙ and $M\dot{M}\lt 6.3\times {10}^{-7}\,{M}_{\mathrm{Jup}}^{2}\,{\mathrm{yr}}^{-1}$ for 1 R Jup . These limits are comparable to the accretion luminosity and accretion rate of T-Tauri stars or LkCa 15 b. Taken together, these lines of evidence suggest that the H band source at the location of HD 100546 b is not emitted by a planetary photosphere or an accreting circumplanetary disk but is a disk feature enhanced by the point-spread function subtraction process. Lastly, this non-detection is consistent with the non-detection in the K band reported in an earlier study but does not exclude the possibility that HD 100546 b is deeply embedded.
Dust growth is often neglected when building models of protoplanetary disks due to its complexity and computational expense. However, it does play a major role in shaping the evolution of protoplanetary dust and planet formation. In this paper, we present a numerical model coupling 2D hydrodynamic evolution of a protoplanetary disk, including a Jupiter-mass planet, and dust coagulation. This is obtained by including multiple dust fluids in a single grid-based hydrodynamic simulation and solving the Smoluchowski equation for dust coagulation on top of solving for the hydrodynamic evolution. We find that fragmentation of dust aggregates trapped in a pressure bump outside of the planetary gap leads to an enhancement in the density of small grains. Additionally, we compare the results obtained from the full-coagulation treatment to the commonly used, fixed-dust-size approach and to previously applied, less computationally intensive methods for including dust coagulation. We find that the full-coagulation results cannot be reproduced using the fixed-size treatment, but some can be mimicked using a relatively simple method for estimating the characteristic dust size in every grid cell.
The omnipresence of super-Earths suggests that they are able to be retained in natal disks around low-mass stars, whereas exoplanets' mass distributions indicate that some cores have transformed into gas giants through runaway gas accretion at ≳1 au from solar-type stars. In this paper, we show that transition to runaway gas accretion by cores may be self-impeded by an increase of the grain opacity in their envelope after they have acquired sufficient mass (typically ~ 10M ⊕ ) to enter a pebble-isolation phase. The accumulation of approximately millimeter- to meter-sized pebbles in their migration barriers enhances their local fragmentation rates. The freshly produced submillimeter grains pass through the barrier, elevate the effective dust opacity, and reduce the radiative flux in the core's envelope. These effects alone are adequate to suppress the transition to runaway accretion and preserve super-Earths in the stellar proximity (~0.1 au), albeit entropy advection between the envelope and the disk can further reduce the accretion rate. At an intermediate distance (~1 au) from their host stars, the escalation in the dust opacity dominates over entropy advection in stalling the transition to runaway accretion for marginally pebble-isolated cores. Beyond a few astronomical unit, the transformation of more-massive cores to gas giants is reachable before severe depletion of disk gas. This requirement can be satisfied either in extended disks with large scale height via orderly accretion of migrating pebbles or through the mergers of oligarchic protoplanetary embryos, and can account for the correlated occurrence of long-period gas giants and close-in super-Earths.
During their formation, emerging protoplanets tidally interact with their natal disks. Proto–gas giant planets, with Hill radii larger than the disk thickness, open gaps and quench gas flow in the vicinity of their orbits. It is usually assumed that their type II migration is coupled to the viscous evolution of the disk. Although this hypothesis provides an explanation for the origin of close-in planets, it also encounters a predicament on the retention of long-period orbits for most gas giant planets. Moreover, numerical simulations indicate that the planets' migrations are not solely determined by the viscous diffusion of their natal disk. In this work, we carry out a series of hydrodynamic simulations combined with analytic studies to examine the transition between different paradigms of type II migration. We find a range of planetary mass for which gas continues to flow through a severely depleted gap so that the surface density distribution in the disk region beyond the gap is maintained in a quasi-steady state. The associated gap profile modifies the location of corotation and Lindblad resonances. In the proximity of the planet's orbit, high-order Lindblad and corotation torque are weakened by the gas depletion in the gap, while low-order Lindblad torques near the gap walls preserve their magnitude. Consequently, the intrinsic surface density distribution of the disk delicately determines both the pace and direction of the planets' type II migration. We show that this effect might stall the inward migration of giant planets and preserve them in disk regions where the surface density is steep.
The tidal perturbation of embedded protoplanets on their natal disks has been widely attributed to be the cause of gap-ring structures in submillimeter images of protoplanetary disks around T Tauri stars. Numerical simulations of this process have been used to propose scaling of characteristic dust-gap width/gap-ring distance with respect to planet mass. Applying such scaling to analyze observed gap samples yields a continuous mass distribution for a rich population of hypothetical planets in the range of several Earth to Jupiter masses. In contrast, the conventional core-accretion scenario of planet formation predicts a bimodal mass function due to (1) the onset of runaway gas accretion above ~20 Earth masses and (2) suppression of accretion induced by gap opening. Here, we examine the dust disk response to the tidal perturbation of eccentric planets as a possible resolution of this paradox. Based on simulated gas and dust distributions, we show the gap-ring separation of Neptune-mass planets with small eccentricities might become comparable to that induced by Saturn-mass planets on circular orbits. Furthermore, this degeneracy may obliterate the discrepancy between the theoretical bimodal mass distribution and the observed continuous gap width distribution. Despite damping due to planet–disk interaction, modest eccentricity may be sustained either in the outer regions of relatively thick disks or through resonant excitation among multiple super Earths. Moreover, the ring-like dust distribution induced by planets with small eccentricities is axisymmetric even in low viscosity environments, consistent with the paucity of vortices in Atacama Large Millimeter/submillimeter Array images.
Here, we present the discovery of a brown dwarf companion to the debris disk host star HR 2562. This object, discovered with the Gemini Planet Imager (GPI), has a projected separation of 20.3 ± 0.3 au ($0\buildrel{\prime\prime}\over{.} 618\pm 0\buildrel{\prime\prime}\over{.} 004$) from the star. With the high astrometric precision afforded by GPI, we have confirmed, to more than 5σ, the common proper motion of HR 2562B with the star, with only a month-long time baseline between observations. Spectral data in the J-, H-, and K-bands show a morphological similarity to L/T transition objects. We assign a spectral type of L7 ± 3 to HR 2562B and derive a luminosity of log(L ${}_{\mathrm{bol}}$/${L}_{\odot })=-4.62\pm 0.12$, corresponding to a mass of 30 ± 15 ${M}_{\mathrm{Jup}}$ from evolutionary models at an estimated age of the system of 300–900 Myr. Although the uncertainty in the age of the host star is significant, the spectra and photometry exhibit several indications of youth for HR 2562B. The source has a position angle that is consistent with an orbit in the same plane as the debris disk recently resolved with Herschel. Additionally, it appears to be interior to the debris disk. Though the extent of the inner hole is currently too uncertain to place limits on the mass of HR 2562B, future observations of the disk with higher spatial resolution may be able to provide mass constraints. This is the first brown-dwarf-mass object found to reside in the inner hole of a debris disk, offering the opportunity to search for evidence of formation above the deuterium burning limit in a circumstellar disk.
In this work, we study how the dust coagulation/fragmentation will influence the evolution and observational appearances of vortices induced by a massive planet embedded in a low-viscosity disk by performing global 2D high-resolution hydrodynamical simulations. Within the vortex, due to its higher gas surface density and steeper pressure gradients, dust coagulation, fragmentation, and drift (to the vortex center) are all quite efficient, producing dust particles ranging from 1 μm to ~1.0 cm, as well as an overall high dust-to-gas ratio (above unity). In addition, the dust size distribution is quite nonuniform inside the vortex, with the mass-weighted average dust size at the vortex center (~4.0 mm) being a factor of ~10 larger than other vortex regions. Both large (~millimeter) and small (tens of microns) particles contribute strongly to affect the gas motion within the vortex. As such, we find that the inclusion of dust coagulation has a significant impact on the vortex lifetime and the typical vortex lifetime is about 1000 orbits. Furthermore, after the initial gaseous vortex is destroyed, the dust spreads into a ring with a few remaining smaller gaseous vortices with a high dust concentration and a large maximum size (~millimeter). At late time, the synthetic dust continuum images for the coagulation case show as a ring inlaid with several hot spots at the 1.33 mm band, while only distinct hot spots remain at 7.0 mm.