Massive Protoplanetary Research Articles

Thermal Structure and Millimeter Emission from a Protoplanetary Disk with Embedded Protoplanets from Radiative Transfer Modeling

The discovery of protoplanets and circumplanetary disks provides a unique opportunity to characterize planet formation through observations. Massive protoplanets shape the physical and chemical structure of their host circumstellar disk by accretion, localized emission, and disk depletion. In this work, we study the thermal changes induced within the disk by protoplanet accretion and synthetic predictions through hydrodynamical simulations with postprocessed radiative transfer with an emphasis on radio millimeter emission. We explored distinct growth conditions and varied both planetary accretion rates and the local dust-to-gas mass ratios for a protoplanet at 1200 K. The radiative transfer models show that beyond the effect of disk gaps, in most cases, the circumplanetary disk (CPD) and the planet’s emission locally increase the disk temperature. Moreover, depending on the local dust-to-gas depletion and accretion rate, the presence of the CPD may have detectable signatures in millimeter emission. It also has the power to generate azimuthal asymmetries that are important for continuum subtraction. Thus, if other means of detection of protoplanets are proven, the lack of corresponding evidence at other wavelengths can set limits on their growth timescales through a combined analysis of the local dust-to-gas ratio and the accretion rate.

Forming giant planets around late-M dwarfs: Pebble accretion and planet–planet collision

We propose a pebble-driven core accretion scenario to explain the formation of giant planets around the late-M dwarfs of M★=0.1– 0.2 M⊙. In order to explore the optimal disk conditions for giant planet, we performed N-body simulations to investigate the growth and dynamical evolution of both single and multiple protoplanets in the disks with both inner viscously heated and outer stellar irradiated regions. The initial masses of the protoplanets are either assumed to be equal to 0.01 M⊕ or calculated based on the formula derived from streaming instability simulations. Our findings indicate that massive planets are more likely to form in disks with longer lifetimes, higher solid masses, moderate to high levels of disk turbulence, and larger initial masses of protoplanets. In the single protoplanet growth cases, the highest planet core mass that can be reached is generally lower than the threshold necessary to trigger rapid gas accretion, which impedes the formation of giant planets. Nonetheless, in multi-protoplanet cases, the cores can exceed the pebble isolation mass barrier aided by frequent planet–planet collisions. This consequently speeds their gas accretion up and promotes giant planet formation, making the optimal parameter space to grow giant planets substantially wider. Taken together, our results suggest that even around very-low-mass stellar hosts, the giant planets with orbital periods of ≲100 days are still likely to form when lunar-mass protoplanets first emerge from planetesimal accretion and then grow rapidly by a combination of pebble accretion and planet–planet collisions in disks with a high supply of a pebble reservoir >50 M⊕ and a turbulent level of αt ~ 10−3−10−2.

Accretion of ice giant planets from massive protoplanets whose migration is blocked by Jupiter and Saturn

ABSTRACT The growth and dynamical evolution of protoplanets beyond Saturn through collisions and type I migration typically result in a highly chaotic dynamics, producing a diversity of outcomes depending on the initial conditions. Here we present the results of N-bodies numerical simulations aiming to make a detailed exploration of different initial conditions and potential outcomes for this dynamics. We consider Jupiter and Saturn at the imminence of crossing the 3:2 mean motion resonance in two possible positions based on the final and initial conditions of the Grand Tack and the Nice models, respectively; four different gas disc lifetimes, and a range of population sizes of planetary embryos beyond Saturn with different masses and orbital configurations in a total of 72 different setups. We present statistical analyses of our outcomes including planets that ‘jump’ to interior orbits of Jupiter, the frequency of close encounters between the planetary embryos and Jupiter, the number of ejections and collisions, and the dynamical effects for stabilizing the final planetary system in chains of mean motion resonances. Results show that independently of the initial configuration of Jupiter and Saturn and the gas lifetime, the dynamical evolution goes through three main phases. A few per cent of the simulations successfully produce Uranus and Neptune analogues, which may have implications on the ice giants’ composition and obliquities, the material ejected from the Solar System, and the conditions for the giant planet instability.

SimAb: A simple, fast, and flexible model to assess the effects of planet formation on the atmospheric composition of gas giants

Context. The composition of exoplanet atmospheres provides us with vital insight into their formation scenario. Conversely, planet formation processes shape the composition of atmospheres and imprint their specific signatures. In this context, models of planet formation containing key formation processes help supply clues to how planets form. This includes constraints on the metallicity and carbon-to-oxygen ratio (C/O ratio) of the planetary atmospheres. Gas giants in particular are of great interest due to the amount of information we can obtain about their atmospheric composition from their spectra, and also due to their relative ease of observation. Aims. We present a basic, fast, and flexible planet formation model, called Simulating Abundances (SimAb), that forms giant planets and allows us to study their primary atmospheric composition soon after their formation. Methods. In SimAb we introduce parameters to simplify the assumptions about the complex physics involved in the formation of a planet. This approach allows us to trace and understand the influence of complex physical processes on the formed planets. In this study we focus on four different parameters and how they influence the composition of the planetary atmospheres: initial protoplanet mass, initial orbital distance of the protoplanet, planetesimal ratio in the disk, and dust grain fraction in the disk. Results. We focus on the C/O ratio and the metallicity of the planetary atmosphere as an indicator of their composition. We show that the initial protoplanet core mass does not influence the final composition of the planetary atmosphere in the context of our model. The initial orbital distance affects the C/O ratio due to the different C/O ratios in the gas phase and the solid phase at different orbital distances. Additionally, the initial orbital distance together with the amount of accreted planetesimals cause the planet to have subsolar or supersolar metallicity. Furthermore, the C/O ratio is affected by the dust grain fraction and the planetesimal ratio. Planets that accrete most of their heavy elements through dust grains will have a C/O ratio close to the solar C/O ratio, while planets that accrete most of their heavy elements from the planetesimals in the disk will end up with a C/O ratio closer to the C/O ratio in the solid phase of the disk. Conclusions. By using the C/O ratio and metallicity together we can put a lower and upper boundary on the initial orbital distance where supersolar metallicity planets are formed. We show that planetesimals are the main source for reaching supersolar metallicity planets. On the other hand, planets that mainly accrete dust grains will show a more solar composition. Supersolar metallicity planets that initiate their formation farther than the CO ice line have a C/O ratio closer to the solar value.

Explaining mercury via a single giant impact is highly unlikely

ABSTRACTThe classical scenario of terrestrial planet formation is characterized by a phase of giant impacts among Moon-to-Mars mass planetary embryos. While the classic model and its adaptations have produced adequate analogues of the outer three terrestrial planets, Mercury’s origin remains elusive. Mercury’s high-core mass fraction compared to the Earth’s is particularly outstanding. Among collisional hypotheses, this feature has been long interpreted as the outcome of an energetic giant impact among two massive protoplanets. Here, we revisit the classical scenario of terrestrial planet formation with focus on the outcome of giant impacts. We have performed a large number of N-body simulations considering different initial distributions of planetary embryos and planetesimals. Our simulations tested the effects of different giant planet configurations, from virtually circular to very eccentric configurations. We compare the giant impacts produced in our simulations with those that are more likely to account for the formation of Mercury and the Moon according to smoothed hydrodynamic simulations. Impact events that could lead to Moon’s formation are observed in all our simulations with up to ∼20 per cent of all giant impacts, consistent with the range of the expected Moon-forming event conditions. On the other hand, Mercury-forming events via a single giant impact are extremely rare, accounting for less than ∼1 per cent of all giant impacts. Our results suggest that producing Mercury as a remnant of a single giant impact that strips out the mantle of a differentiated planetary object with Earth-like iron-silicate ratio is challenging and alternative scenarios may be required (e.g. multiple collisions).

Giant Planet’s Formation Caught in Action

Astronomers took a direct image of a massive protoplanet embedded in a protoplanetary disk. The system provides strong evidence for an as-yet-unconfirmed theory of planet formation.

The Three-dimensional Flow Field around Planets on Eccentric Orbits

Abstract We investigate the properties of the hydrodynamic flow around eccentric protoplanets and compare them with the often assumed case of a circular orbit. To this end, we perform a set of 3D hydrodynamic simulations of protoplanets with small eccentricities (e ≤ 0.1). We adopt an isothermal equation of state and concentrate resolution on the protoplanet to investigate flows down to the scale of the protoplanet's circumplanetary disk (CPD). We find enhanced prograde rotation exterior to the CPD for low planet masses undergoing subsonic eccentric motion. If the eccentricity is made large enough to develop a bow shock, this trend reverses and rotation becomes increasingly retrograde. The instantaneous eccentric flow field is dramatically altered compared to circular orbits. Whereas the latter exhibits a generic pattern of polar inflow and midplane outflow, the flow geometry depends on the orbital phase in the eccentric case. For even the modest eccentricities tested here, the dominant source of inflow can come from the midplane instead of the poles. We find that the amount of inflow and outflow increases for higher e and lower protoplanet masses, thereby recycling more gas through the planet’s Bondi radius. These increased fluxes may increase the pebble accretion rate for eccentric planets up to several times that of the circular orbit rate. In response to the eccentric motion, the structure and rotation of the planet’s bound CPD remain unchanged. Because the CPD regulates the eventual accretion of gas onto the planet, we predict little change to the gas accretion rates between eccentric and circular planets.

In situ formation of hot Jupiters with companion super-Earths

ABSTRACT Observations have confirmed the existence of multiple-planet systems containing a hot Jupiter and smaller planetary companions. Examples include WASP-47, Kepler-730, and TOI-1130. We examine the plausibility of forming such systems in situ using N-body simulations that include a realistic treatment of collisions, an evolving protoplanetary disc, and eccentricity/inclination damping of planetary embryos. Initial conditions are constructed using two different models for the core of the giant planet: a ‘seed-model’ and an ‘equal-mass-model’. The former has a more massive protoplanet placed among multiple small embryos in a compact configuration. The latter consists only of equal-mass embryos. Simulations of the seed-model lead to the formation of systems containing a hot Jupiter and super-Earths. The evolution consistently follows four distinct phases: early giant impacts; runaway gas accretion on to the seed protoplanet; disc damping-dominated evolution of the embryos orbiting exterior to the giant; a late chaotic phase after dispersal of the gas disc. Approximately 1 per cent of the equal-mass simulations form a giant and follow the same four-phase evolution. Synthetic transit observations of the equal-mass simulations provide an occurrence rate of 0.26 per cent for systems containing a hot Jupiter and an inner super-Earth, similar to the 0.2 per cent occurrence rate from actual transit surveys, but simulated hot Jupiters are rarely detected as single transiting planets, in disagreement with observations. A subset of our simulations form two close-in giants, similar to the WASP-148 system. The scenario explored here provides a viable pathway for forming systems with unusual architectures, but does not apply to the majority of hot Jupiters.

Searching for wide-orbit gravitational instability protoplanets with ALMA in the dust continuum

ABSTRACT Searches for young gas giant planets at wide separations have so far focused on techniques appropriate for compact (Jupiter-sized) planets. Here, we point out that protoplanets born through gravitational instability (GI) may remain in an initial pre-collapse phase for as long as the first 105–107 yr after formation. These objects are hundreds of times larger than Jupiter and their atmospheres are too cold (T ∼ tens of K) to emit in the near-infrared or Hα via accretion shocks. However, it is possible that their dust emission can be detected with Atacama Large Millimeter/submillimeter Array (ALMA), even around Classes I and II protoplanetary discs. In this paper, we produce synthetic observations of these protoplanets. We find that making a detection in a disc at 140 pc would require a few hundred minutes of ALMA band 6 observation time. Protoplanets with masses of 3–5 MJ have the highest chance of being detected; less massive objects require unreasonably long observation times (1000 min), while more massive ones collapse into giant planets before 105 yr. We propose that high-resolution surveys of young (105–106 yr), massive and face on discs offer the best chance for observing protoplanets. Such a detection would help to place constraints on the protoplanet mass spectrum, explain the turnover in the occurrence frequency of gas giants with system metallicity and constrain the prevalence of GI as a planet formation mechanism. Consistent lack of detection would be evidence against GI as a common planet formation mechanism.

Accretion of Gas Giants Constrained by the Tidal Barrier

Abstract After protoplanets have acquired sufficient mass to open partial gaps in their natal protostellar disks, residual gas continues to diffuse onto horseshoe streamlines under the effect of viscous dissipation, and to meander in and out of the planets’ Hill sphere. Within the Hill sphere, the horseshoe streamlines intercept gas flow in circumplanetary disks. The host star’s tidal perturbation induces a barrier across the converging streamlines’ interface. Viscous transfer of angular momentum across this tidal barrier determines the rate of mass diffusion from the horseshoe streamlines onto the circumplanetary disks, and eventually the accretion rate onto the protoplanets. We carry out a series of numerical simulations to test the influence of this tidal barrier on superthermal planets. In weakly viscous disks, the protoplanets’ accretion rate steeply decreases with their masses above the thermal limit. As their growth timescale exceeds the gas depletion timescale, their masses reach asymptotic values comparable to that of Jupiter. In relatively thick and strongly viscous disks, protoplanets’ asymptotic masses exceed several times that of Jupiter. Two-dimensional numerical simulations show that such massive protoplanets strongly excite the eccentricity of nearby horseshoe streamlines, destabilize orderly flow, substantially enhance the diffusion rate across the tidal barrier, and elevate their growth rate until their natal disk is severely depleted. In contrast, eccentric streamlines remain stable in three-dimensional simulations. Based on the upper falloff in the observed mass distribution of known exoplanets, we suggest that their natal disks had relatively low viscosity (α ∼ 10−3) and modest thickness (H/R ∼ 0.03–0.05).

Can Chondrules Be Produced by the Interaction of Jupiter with the Protosolar Disk?

Abstract Chondrules are crystallized droplets of silicate melt formed by rapid heating to high temperatures (>1800 K) of solid precursors followed by hours or days of cooling. The time interval of chondrule formation is consistent with the formation timescale of Jupiter in the core-accretion model (1–4 Myr). Here we investigate if the shocks generated by a massive planet could generate flash heating episodes necessary to form chondrules using high-resolution 2D simulations with the multifluid code RoSSBi. We use different radiative cooling prescriptions, planet masses, orbits, and disk models. Temperatures reached during flash heating can be deduced from chondrule observations and are achieved in a Minimum Mass Solar Nebula (MMSN) for a massive protoplanet (>0.75 M ♃) but only in cases in which radiative cooling is low enough to lead to nearly adiabatic conditions. More realistic thermodynamics undershoot the temperatures required in shocks for chondrule formation. However, these temperatures are reached when considering more massive disks (e.g., five MMSN), but these conditions lead to fast planet migration and too low cooling rates compared to those deduced from chondrule textures. Thus, it seems unlikely that shocks from Jupiter can form chondrules in most cases. Independent of the nebular mass, the simulations demonstrate that a massive planet that forms a gap triggers vortices, which act as dust traps for chondrule precursors. These vortices also provide a high-pressure environment consistent with cosmochemical evidence from chondrules. They only lack the flash heating source for a potential chondrule formation environment.

Segregation of heavy elements in an initial protoplanet

In this paper, we present a calculation for settling times of grains falling under the action of the gravitational field inside a Jupiter mass protoplanet formed via disk instability, where the grains are assumed to be grown by the process of cold welding. As the energy equation, we have considered the conductive-radiative case of heat transport. In our calculation, we have not used any density model, rather, we have estimated the distribution of physical variables inside the protoplanet. With the obtained distribution of the physical variables, we have investigated the growth of the grains having different initial radii (5 × 10−3- 10−1cm) and hence calculated their sedimentation times. The results are found to be in good agreement with those obtained in other investigations.

Exploring the conditions for forming cold gas giants through planetesimal accretion

The formation of cold gas giants similar to Jupiter and Saturn in orbit and mass is a great challenge for planetesimal-driven core accretion models because the core growth rates far from the star are low. Here we model the growth and migration of single protoplanets that accrete planetesimals and gas. We integrated the core growth rate using fits in the literature to N-body simulations, which provide the efficiency of accreting the planetesimals that a protoplanet migrates through. We take into account three constraints from the solar system and from protoplanetary discs: (1) the masses of the terrestrial planets and the comet reservoirs in Neptune’s scattered disc and the Oort cloud are consistent with a primordial planetesimal population of a few Earth masses per AU, (2) evidence from the asteroid belt and the Kuiper belt indicates that the characteristic planetesimal diameter is 100 km, and (3) observations of protoplanetary discs indicate that the dust is stirred by weak turbulence; this gas turbulence also excites the inclinations of planetesimals. Our nominal model built on these constraints results in maximum protoplanet masses of 0.1 Earth masses. Ignoring constraint (1) above, we show that even a planetesimal population of 1000 Earth masses, corresponding to 50 Earth masses per AU, fails to produce cold gas giants (although it successfully forms hot and warm gas giants). We conclude that a massive planetesimal reservoir is in itself insufficient to produce cold gas giants. The formation of cold gas giants by planetesimal accretion additionally requires that planetesimals are small and that the turbulent stirring is very weak, thereby violating all three above constraints.

Growth after the streaming instability

Context. Streaming instability is a key mechanism in planet formation, clustering pebbles into planetesimals with the help of self-gravity. It is triggered at a particular disk location where the local volume density of solids exceeds that of the gas. After their formation, planetesimals can grow into protoplanets by feeding from other planetesimals in the birth ring as well as by accreting inwardly drifting pebbles from the outer disk. Aims. We aim to investigate the growth of planetesimals into protoplanets at a single location through streaming instability. For a solar-mass star, we test the conditions under which super-Earths are able to form within the lifetime of the gaseous disk. Methods. We modified the Mercury N-body code to trace the growth and dynamical evolution of a swarm of planetesimals at a distance of 2.7 AU from the star. The code simulates gravitational interactions and collisions among planetesimals, gas drag, type I torque, and pebble accretion. Three distributions of planetesimal sizes were investigated: (i) a mono-dispersed population of 400 km radius planetesimals, (ii) a poly-dispersed population of planetesimals from 200 km up to 1000 km, (iii) a bimodal distribution with a single runaway body and a swarm of smaller, 100 km size planetesimals. Results. The mono-dispersed population of 400 km size planetesimals cannot form protoplanets of a mass greater than that of the Earth. Their eccentricities and inclinations are quickly excited, which suppresses both planetesimal accretion and pebble accretion. Planets can form from the poly-dispersed and bimodal distributions. In these circumstances, it is the two-component nature that damps the random velocity of the large embryo through the dynamical friction of small planetesimals, allowing the embryo to accrete pebbles efficiently when it approaches 10−2 M⊕. Accounting for migration, close-in super-Earth planets form. Super-Earth planets are likely to form when the pebble mass flux is higher, the disk turbulence is lower, or the Stokes number of the pebbles is higher. Conclusions. For the single site planetesimal formation scenario, a two-component mass distribution with a large embryo and small planetesimals promotes planet growth, first by planetesimal accretion and then by pebble accretion of the most massive protoplanet. Planetesimal formation at single locations such as ice lines naturally leads to super-Earth planets by the combined mechanisms of planetesimal accretion and pebble accretion.

Multiple Spiral Arms in the Disk around Intermediate-mass Binary HD 34700A

We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearly see significant non-azimuthal polarization Uphi consistent with multiple scattering within a disk at an inferred inclination ~42deg. Radiative transfer modeling along with a new Gaia distance suggest HD 37400A is a young (~5 Myr) system consisting of two intermediate-mass (~2Msun) stars surrounded by a transitional disk and not a solar-mass binary with a debris disk as previously classified. Conventional assumptions of the dust-to-gas ratio would rule out a gravitational instability origin to the spirals while hydrodynamical models using the known external companion or a hypothetical massive protoplanet in the cavity both have trouble reproducing the relatively large spiral arm pitch angles (~30deg) without fine tuning of gas temperature. We explore the possibility that material surrounding a massive protoplanet could explain the rim discontinuity after also considering effects of shadowing by an inner disk. Analysis of archival Hubble Space Telescope data suggests the disk is rotating counterclockwise as expected from the spiral arm structure and revealed a new low-mass companion at 6.45" separation. We include an appendix which sets out clear definitions of Q, U, Qphi, Uphi, correcting some confusion and errors in the literature.

Dynamics of multiple protoplanets embedded in gas and pebble discs and its dependence on Σ and ν parameters

Protoplanets of super-Earth size may get trapped in convergence zones for planetary migration and form gas giants there. These growing planets undergo accretion heating, which triggers a hot-trail effect that can reverse migration directions, increase planetary eccentricities, and prevent resonant captures of migrating planets. In this work, we study populations of embryos that are accreting pebbles under different conditions, by changing the surface density, viscosity, pebble flux, mass, and the number of protoplanets. For modelling, we used the FARGO-THORIN two-dimensional (2D) hydrocode, which incorporates a pebble disc as a second pressure-less fluid, the coupling between the gas and pebbles, and the flux-limited diffusion approximation for radiative transfer. We find that massive embryos embedded in a disc with high surface density (Σ = 990 g cm−2 at 5.2 au) undergo numerous “unsuccessful” two-body encounters that do not lead to a merger. Only when a third protoplanet arrives in the convergence zone do three-body encounters lead to mergers. For a low-viscosity disc (ν = 5 × 1013 cm2 s−1), a massive co-orbital is a possible outcome, for which a pebble isolation develops and the co-orbital is further stabilised. For more massive protoplanets (5 M⊕), the convergence radius is located further out, in the ice-giant zone. After a series of encounters, there is an evolution driven by a dynamical torque of a tadpole region, which is systematically repeated several times until the co-orbital configuration is disrupted and planets merge. This may be a way to solve the problem that co-orbitals often form in simulations but they are not observed in nature. In contrast, the joint evolution of 120 low-mass protoplanets (0.1 M⊕) reveals completely different dynamics. The evolution is no longer smooth, but rather a random walk. This is because the spiral arms, developed in the gas disc due to Lindblad resonances, overlap with each other and affect not only a single protoplanet but several in the surrounding area. Our hydrodynamical simulations may have important implications for N-body simulations of planetary migration that use simplified torque prescriptions and are thus unable to capture protoplanet dynamics in its full glory.

Theoretical Model of Hydrogen Line Emission from Accreting Gas Giants

Abstract Progress in understanding of giant planet formation has been hampered by a lack of observational constraints to growing protoplanets. Recently, detection of an Hα-emission excess via direct imaging was reported for the protoplanet LkCa15b orbiting the pre-main-sequence star LkCa15. However, the physical mechanism for the Hα emission is poorly understood. According to recent high-resolution three-dimensional hydrodynamic simulations of the flow accreting onto protoplanets, the disk gas flows down almost vertically onto and collides with the surface of a circumplanetary disk at a supersonic velocity and thus passes through a strong shockwave. The shock-heated gas is hot enough to generate Hα emission. Here we develop a one-dimensional radiative hydrodynamic model of the flow after the shock by detailed calculations of chemical reactions and electron transitions in hydrogen atoms, and quantify the hydrogen line emission in the Lyman-, Balmer-, and Paschen-series from the accreting gas giant system. We then demonstrate that the Hα intensity is strong enough to be detected with current observational techniques. Comparing our theoretical Hα intensity with the observed one from LkCa15b, we constrain the protoplanet mass and the disk gas density. Observation of hydrogen line emission from protoplanets is highly encouraged to obtain direct constraints of accreting gas giants, which will be key in understanding the formation of gas giants.

Planet Formation: An Optimized Population-synthesis Approach

Abstract The physics of planet formation is investigated using a population synthesis approach. We develop a simple model for planetary growth including pebble and gas accretion, as well as orbital migration in an evolving protoplanetary disk. The model is run for a population of 2000 stars with a range of disk masses, disk radii, and initial protoplanet orbits. The resulting planetary distribution is compared with the observed population of extrasolar planets, and the model parameters are improved iteratively using a particle swarm optimization scheme. The characteristics of the final planetary systems are mainly controlled by the pebble isolation mass, which is the mass of a planet that perturbs nearby gas enough to halt the inward flux of drifting pebbles and stop growth. The pebble isolation mass increases with orbital distance such that giant planet cores can only form in the outer disk. Giants migrate inward, populating a wide range of final orbital distances. The best model fits have large initial protoplanet masses, short pebble drift timescales, low disk viscosities, and short atmospheric cooling times, all of which promote rapid growth. The model successfully reproduces the observed frequency and distribution of giant planets and brown dwarfs. The fit for super-Earths is poorer for single-planet systems, but improves steadily when more protoplanets are included. Although the study was designed to match the extrasolar planet distribution, analogs of the solar system form in 1–2% of systems that contain at least four protoplanets.

The diverse lives of massive protoplanets in self-gravitating discs

Gas giant planets may form early-on during the evolution of protostellar discs, while these are relatively massive. We study how Jupiter-mass planet-seeds (termed protoplanets) evolve in massive, but gravitationally stable (Q>1.5), discs using radiative hydrodynamic simulations. We find that the protoplanet initially migrates inwards rapidly, until it opens up a gap in the disc. Thereafter, it either continues to migrate inwards on a much longer timescale or starts migrating outwards. Outward migration occurs when the protoplanet resides within a gap with gravitationally unstable edges, as a high fraction of the accreted gas is high angular momentum gas from outside the protoplanet's orbit. The effect of radiative heating from the protoplanet is critical in determining the direction of the migration and the eccentricity of the protoplanet. Gap opening is facilitated by efficient cooling that may not be captured by the commonly used beta-cooling approximation. The protoplanet initially accretes at a high rate (1e-3Mj/yr), and its accretion luminosity could be a few tenths of the host star's luminosity, making the protoplanet easily observable (albeit only for a short time). Due to the high gas accretion rate, the protoplanet generally grows above the deuterium-burning mass-limit. Protoplanet radiative feedback reduces its mass growth so that its final mass is near the brown dwarf-planet boundary. The fate of a young planet-seed is diverse and could vary from a gas giant planet on a circular orbit at a few AU from the central star to a brown dwarf on an eccentric, wide orbit.

ON SHOCKS DRIVEN BY HIGH-MASS PLANETS IN RADIATIVELY INEFFICIENT DISKS. I. TWO-DIMENSIONAL GLOBAL DISK SIMULATIONS

Recent observations of gaps and non-axisymmetric features in the dust distributions of transition disks have been interpreted as evidence of embedded massive protoplanets. However, comparing the predictions of planet–disk interaction models to the observed features has shown far from perfect agreement. This may be due to the strong approximations used for the predictions. For example, spiral arm fitting typically uses results that are based on low-mass planets in an isothermal gas. In this work, we describe two-dimensional, global, hydrodynamical simulations of disks with embedded protoplanets, with and without the assumption of local isothermality, for a range of planet-to-star mass ratios 1–10 for a 1 star. We use the Pencil Code in polar coordinates for our models. We find that the inner and outer spiral wakes of massive protoplanets () produce significant shock heating that can trigger buoyant instabilities. These drive sustained turbulence throughout the disk when they occur. The strength of this effect depends strongly on the mass of the planet and the thermal relaxation timescale; for a planet embedded in a thin, purely adiabatic disk, the spirals, gaps, and vortices typically associated with planet–disk interactions are disrupted. We find that the effect is only weakly dependent on the initial radial temperature profile. The spirals that form in disks heated by the effects we have described may fit the spiral structures observed in transition disks better than the spirals predicted by linear isothermal theory.