Runaway Accretion Research Articles

The Not-so Dramatic Effect of Advective Flows on Gas Accretion

Super-Earths and mini-Neptunes are the most common types of exoplanets discovered, yet the physics of their formation are still debated. Standard core accretion models in gas-rich environments find that typical mini-Neptune mass planets would blow up into Jupiters before the underlying disk gas dissipates away. The injection of entropy from the protoplanetary disk into forming gaseous envelopes has recently been put forward as a mechanism to delay this runaway accretion, specifically at short orbital distances. Here, we reevaluate this line of reasoning by incorporating recycling flows of gas into a numerical one-dimensional thermodynamic model with a more realistic equation of state and opacities and the thermal state of the advective flow. At 0.1 au, we find that advective flows are only able to produce mini-Neptunes if they can penetrate below ∼0.25 of the planet’s gravitational sphere of influence. Otherwise, the gas-to-core mass ratio (GCR) reaches above ∼10%, which is too large to explain the measured properties of mini-Neptunes, necessitating other gas-limiting processes such as late-time core assembly. The effect of entropy advection on gas accretion weakens even further beyond 0.1 au. We present an updated scaling relation between GCR and the penetration depth of the advective flows, which varies nontrivially with orbital distances, core masses, and dusty versus dust-free opacity. We further demonstrate how measurements of planet mass distribution beyond ∼1 au using future instruments such as the Nancy Grace Roman Space Telescope could be used to disambiguate between different formation conditions of gas-poor planets.

Supermassive black holes from runaway mergers and accretion in nuclear star clusters

Abstract Rapid formation of supermassive black holes occurs in dense nuclear star clusters that are initially gas-dominated. Stellar-mass black hole remnants of the most massive cluster stars sink into the core, where a massive runaway black hole forms as a consequence of combined effects of repeated mergers and Eddington-limited gas accretion. The associated gravitational-wave signals of high-redshift extreme mass-ratio inspirals are a unique signature of the nuclear star cluster scenario.

A Mini-Neptune Orbiting the Metal-poor K Dwarf BD+29 2654

We report the discovery and Doppler mass measurement of a 7.4 days 2.3 R ⊕ mini-Neptune around a metal-poor K dwarf BD+29 2654 (TOI-2018). Based on a high-resolution Keck/HIRES spectrum, the Gaia parallax, and multiwavelength photometry from the UV to the mid-infrared, we found that the host star has K, , [Fe/H] = − 0.58 ± 0.18, M * = 0.57 ± 0.02 M ⊙, and R * = 0.62 ± 0.01 R ⊙. Precise Doppler measurements with Keck/HIRES revealed a planetary mass of M p = 9.2 ± 2.1 M ⊕ for TOI-2018 b. TOI-2018 b has a mass and radius that are consistent with an Earthlike core, with a ∼1%-by-mass hydrogen/helium envelope or an ice–rock mixture. The mass of TOI-2018 b is close to the threshold for runaway accretion and hence giant planet formation. Such a threshold is predicted to be around 10M ⊕ or lower for a low-metallicity (low-opacity) environment. If TOI-2018 b is a planetary core that failed to undergo runaway accretion, it may underline the reason why giant planets are rare around low-metallicity host stars (one possibility is their shorter disk lifetimes). With a K-band magnitude of 7.1, TOI-2018 b may be a suitable target for transmission spectroscopy with the James Webb Space Telescope. The system is also amenable to metastable Helium observation; the detection of a Helium exosphere would help distinguish between a H/He-enveloped planet and a water world.

GRB 211211A: A Neutron Star–White Dwarf Merger?

The gamma-ray burst GRB 211211A and its associated kilonova-like emission were reported recently. A significant difference between this association event and GRB 170817A/AT 2017gfo is that GRB 211211A has a very long duration. In this Letter, we show that this association event may arise from a neutron star–white dwarf (NS–WD) merger if the central engine leaves a magnetar behind. Within the NS–WD merger, the main burst of GRB 211211A could be produced by magnetic bubble eruptions from toroidal magnetic field amplification of the premerger NS. This toroidal field amplification can be induced by the runaway accretion from the WD debris disk if the disk is in low initial entropy and efficient wind, while the extended emission of GRB 211211A is likely involved with magnetic propelling. The observed energetics and duration of the prompt emission of GRB 211211A can be fulfilled in comparison with those of accretion in a hydrodynamical thermonuclear simulation, as long as the WD has a mass ≳1M ⊙. Moreover, if the X-ray plateau in GRB afterglows is due to the magnetar spin-down radiation, GRB optical afterglows and kilonova-like emission can be jointly well modeled combining the standard forward shock with the radioactive decay power of 56Ni adding a rotational power input from the postmerger magnetar.

Modelling the evolution of silicate/volatile accretion discs around white dwarfs

ABSTRACT A growing number of debris discs have been detected around metal-polluted white dwarfs. They are thought to be originated from tidally disrupted exoplanetary bodies and responsible for metal accretion on to host WDs. To explain (1) the observationally inferred accretion rate higher than that induced by Poynting–Robertson drag, $\dot{M}_{\rm PR}$, and (2) refractory-rich photosphere composition indicating the accretion of terrestrial rocky materials, previous studies proposed runaway accretion of silicate particles due to gas drag by the increasing silicate vapour produced by the sublimation of the particles. Because re-condensation of the vapour diffused beyond the sublimation line was neglected, we revisit this problem by one-dimensional advection/diffusion simulation that consistently incorporates silicate sublimation/condensation and back-reaction to particle drift due to gas drag in the solid-rich disc. We find that the silicate vapour density in the region overlapping the solid particles follows the saturating vapour pressure and that no runaway accretion occurs if the re-condensation is included. This always limits the accretion rate from mono-compositional silicate discs to $\dot{M}_{\rm PR}$ in the equilibrium state. Alternatively, by performing additional simulations that couple the volatile gas (e.g. water vapour), we demonstrate that the volatile gas enhances the silicate accretion to $\gt \dot{M}_{\rm PR}$ through gas drag. The refractory-rich accretion is simultaneously reproduced when the initial volatile fraction of disc is $\lesssim 10$ wt per cent because of the suppression of volatile accretion due to the efficient back-reaction of solid to gas. The discs originating from C-type asteroid analogues might be a possible clue to the high-$\dot{M}$ puzzle.

TOI-3757 b: A Low-density Gas Giant Orbiting a Solar-metallicity M Dwarf

We present the discovery of a new Jovian-sized planet, TOI-3757 b, the lowest-density transiting planet known to orbit an M dwarf (M0V). This planet was discovered around a solar-metallicity M dwarf, using Transiting Exoplanet Survey Satellite photometry and confirmed with precise radial velocities from the Habitable-zone Planet Finder (HPF) and NEID. With a planetary radius of 12.0 R ⊕ and mass of 85.3 M ⊕, not only does this object add to the small sample of gas giants (∼10) around M dwarfs, but also its low density ( g cm−3) provides an opportunity to test theories of planet formation. We present two hypotheses to explain its low density; first, we posit that the low metallicity of its stellar host (∼0.3 dex lower than the median metallicity of M dwarfs hosting gas giants) could have played a role in the delayed formation of a solid core massive enough to initiate runaway accretion. Second, using the eccentricity estimate of 0.14 ± 0.06, we determine it is also plausible for tidal heating to at least partially be responsible for inflating the radius of TOI-3757b b. The low density and large scale height of TOI-3757 b makes it an excellent target for transmission spectroscopy studies of atmospheric escape and composition (transmission spectroscopy measurement of ∼ 190). We use HPF to perform transmission spectroscopy of TOI-3757 b using the helium 10830 Å line. Doing this, we place an upper limit of 6.9% (with 90% confidence) on the maximum depth of the absorption from the metastable transition of He at ∼10830 Å, which can help constraint the atmospheric mass-loss rate in this energy-limited regime.

Giant planet formation in Palatini gravity

Some parts of the accretion model of the jovian planets' formation are studied in the context of Palatini gravity. We mainly focus on the critical core mass, that is, a mass for which there is no hydrostatic equilibrium solution for the planet's envelope, which is a starting point of the runaway accretion. We also discuss the conditions needed to be satisfied by a planet such that it can posses a massive gaseous envelope around a solid core.

The Critical Core Mass of Rotating Planets

Abstract The gravitational harmonics measured from the Juno and Cassini spacecraft help us specify the internal structure and chemical elements of Jupiter and Saturn, respectively. However, we still do not know much about the impact of rotation on the planetary internal structure as well as on their formation. The centrifugal force induced by the rotation deforms the planetary shape and partially counteracts the gravitational force. Thus, rotation will affect the critical core mass of the exoplanet. Once the atmospheric mass becomes comparable to the critical core mass, the planet will enter the runaway accretion phase and become a gas giant. We have confirmed that the critical core masses of rotating planets depend on the stiffness of the polytrope, the outer boundary conditions, and the thickness of the isothermal layer. The critical core mass with the Bondi boundary condition is determined by the surface properties. The critical core mass of a rotating planet will increase with the core gravity (i.e., the innermost density). For the Hill boundary condition, the soft polytrope shares the same properties as planets with the Bondi boundary condition. Because the total mass for planets with the Hill boundary condition increases with the decrease of the polytropic index, a higher core gravity is required for rotating planets. As a result, the critical core mass in the stiff Hill model sharply increases. The rotational effects become more important when the radiative and convective regions coexist. Further, the critical core mass of planets with the Hill (Bondi) boundary increases noticeably as the radiative layer becomes thinner (thicker).

In Situ Formation of Super-Earth/Sub-Neptune Driven by the Planetary Rotation

Abstract Kepler’s observation shows that many of the detected planets are super-Earths. They are inside a range of critical masses overlapping the core masses (2–20 M ⊕), which would trigger the runaway accretion and develop the gas giants. Thus, super-Earths/sub-Neptunes can be formed by restraining runaway growth of gaseous envelopes. We assess the effect of planetary rotation in delaying the mass growth. The centrifugal force, induced by spin, will offset a part of the gravitational force and deform the planet. Tracking the change in structure, we find that the temperature at the radiative–convective boundary (RCB) is approximate to the boundary temperature. Since rotation reduces the radiation energy densities in the convective and radiative layers, RCB will penetrate deeper. The cooling luminosity would decrease. Under this condition, the evolutionary timescale can exceed the disk lifetime (10 Myr), and a super-Earth/sub-Neptune could be formed after undergoing additional mass-loss processes. In the dusty atmosphere, even a lower angular velocity can also promote a super-Earth/sub-Neptune forming. Therefore, we conclude that rotation can slow down the planet’s cooling and then promote a super-Earth/sub-Neptune forming.

The effect of late giant collisions on the atmospheres of protoplanets and the formation of cold sub-Saturns

ABSTRACT We investigate the origins of cold sub-Saturns (CSSs), an exoplanetary population inferred from microlensing surveys. If confirmed, these planets would rebut a theorized gap in planets’ mass distribution between those of Neptune and Jupiter caused by the rapid runaway accretion of supercritical cores. In an attempt to resolve this theoretical-observational disparity, we examine the outcomes of giant collisions between sub-critical protoplanets. Due to the secular interaction among protoplanets, these events may occur in rapidly depleting discs. We show that impactors ∼5 per cent the mass of near-runaway envelopes around massive cores can efficiently remove these envelopes entirely via a thermally driven super-Eddington wind emanating from the core itself, in contrast with the stellar Parker winds usually considered. After a brief cooling phase, the merged cores resume accretion. But, the evolution time-scale of transitional discs is too brief for the cores to acquire sufficiently massive envelopes to undergo runaway accretion despite their large combined masses. Consequently, these events lead to the emergence of CSS without their transformation into gas giants. We show that these results are robust for a wide range of disc densities, grain opacities, and silicate abundance in the envelope. Our fiducial case reproduces CSS with heavy (≳30 M⊕) cores and less massive (a few M⊕) sub-critical envelopes. We also investigate the other limiting cases, where continuous mergers of comparable-mass cores yield CSS with wider ranges of core-to-envelope mass ratios and envelope opacities. Our results indicate that it is possible for CSS and Uranus and Neptune to emerge within the framework of well-studied processes, and they may be more common than previously postulated.

TKS X: Confirmation of TOI-1444b and a Comparative Analysis of the Ultra-short-period Planets with Hot Neptunes

Abstract We report the discovery of TOI-1444b, a 1.4 R ⊕ super-Earth on a 0.47 day orbit around a Sun-like star discovered by TESS. Precise radial velocities from Keck/HIRES confirmed the planet and constrained the mass to be 3.87 ± 0.71M ⊕. The RV data set also indicates a possible nontransiting, 16 day planet (11.8 ± 2.9M ⊕). We report a tentative detection of phase-curve variation and a secondary eclipse of TOI-1444b in the TESS bandpass. TOI-1444b joins the growing sample of 17 ultra-short-period planets (USPs) with well-measured masses and sizes, most of which are compatible with an Earth-like composition. We take this opportunity to examine the expanding sample of ultra-short-period planets (<2R ⊕) and contrast them with the newly discovered sub-day ultrahot Neptunes (>3R ⊕, >2000F ⊕ TOI-849 b, LTT9779 b, and K2-100). We find that (1) USPs have predominately Earth-like compositions with inferred iron core mass fractions of 0.32 ± 0.04 and have masses below the threshold of runaway accretion (∼10M ⊕), while ultrahot Neptunes are above the threshold and have H/He or other volatile envelopes. (2) USPs are almost always found in multi-planet systems consistent with a secular interaction formation scenario; ultrahot Neptunes (P orb ≲1 day) tend to be “lonely,” similar to longer-period hot Neptunes (P orb1–10 days) and hot Jupiters. (3) USPs occur around solar-metallicity stars while hot Neptunes prefer higher metallicity hosts. (4) In all these respects, ultrahot Neptunes show more resemblance to hot Jupiters than the smaller USP planets, although ultrahot Neptunes are rarer than both USPs and hot Jupiters by 1–2 orders of magnitude.

New mid-infrared imaging constraints on companions and protoplanetary disks around six young stars

Context. Mid-infrared (mid-IR) imaging traces the sub-micron and micron-sized dust grains in protoplanetary disks and it offers constraints on the geometrical properties of the disks and potential companions, particularly if those companions have circumplanetary disks. Aims. We use the VISIR instrument and its upgrade NEAR on the VLT to take new mid-IR images of five (pre-)transition disks and one circumstellar disk with proposed planets and obtain the deepest resolved mid-IR observations to date in order to put new constraints on the sizes of the emitting regions of the disks and the presence of possible companions. Methods. We derotated and stacked the data to find the disk properties. Where available, we compare the data to PRODIMO (Protoplanetary Disk Model) radiation thermo-chemical models to achieve a deeper understanding of the underlying physical processes within the disks. We applied the circularised point spread function subtraction method to find upper limits on the fluxes of possible companions and model companions with circumplanetary disks. Results. We resolved three of the six disks and calculated position angles, inclinations, and (upper limits to) sizes of emission regions in the disks, improving upper limits on two of the unresolved disks. In all cases the majority of the mid-IR emission comes from small inner disks or the hot inner rims of outer disks. We refined the existing PRODIMO HD 100546 model spectral energy distribution (SED) fit in the mid-IR by increasing the PAH abundance relative to the ISM, adopting coronene as the representative PAH, and increasing the outer cavity radius to 22.3 AU. We produced flux estimates for putative planetary-mass companions and circumplanetary disks, ruling out the presence of planetary-mass companions with L > 0.0028 L⊙ for a > 180 AU in the HD 100546 system. Upper limits of 0.5–30 mJy are obtained at 8–12 μm for potential companions in the different disks. We rule out companions with L > 10−2 L⊙ for a > 60 AU in TW Hydra, a > 110 AU in HD 169142, a > 150 AU in HD 163296, and a > 160 AU in HD 36112. Conclusions. The mid-IR emission comes from the central regions and traces the inner areas of the disks, including inner disks and inner rims of outer disks. Planets with mid-IR luminosities corresponding to a runaway accretion phase can be excluded from the HD 100546, HD 169142, TW Hydra, and HD 36112 systems at separations >1′′. We calculated an upper limit to the occurrence rate of wide-orbit massive planets with circumplanetary disks of 6.2% (68% confidence). Future observations with METIS on the ELT will be able to achieve a factor of 10 better sensitivity with a factor of five better spatial resolution. MIRI on JWST will be able to achieve 250 times better sensitivity. Both will possibly detect the known companions to all six targets.

STARFORGE: the effects of protostellar outflows on the IMF

ABSTRACT The initial mass function (IMF) of stars is a key quantity affecting almost every field of astrophysics, yet it remains unclear what physical mechanisms determine it. We present the first runs of the STAR FORmation in Gaseous Environments project, using a new numerical framework to follow the formation of individual stars in giant molecular clouds (GMCs) using the gizmo code. Our suite includes runs with increasingly complex physics, starting with isothermal ideal magnetohydrodynamics (MHD) and then adding non-isothermal thermodynamics and protostellar outflows. We show that without protostellar outflows the resulting stellar masses are an order of magnitude too high, similar to the result in the base isothermal MHD run. Outflows disrupt the accretion flow around the protostar, allowing gas to fragment and additional stars to form, thereby lowering the mean stellar mass to a value similar to that observed. The effect of jets upon global cloud evolution is most pronounced for lower mass GMCs and dense clumps, so while jets can disrupt low-mass clouds, they are unable to regulate star formation in massive GMCs, as they would turn an order unity fraction of the mass into stars before unbinding the cloud. Jets are also unable to stop the runaway accretion of massive stars, which could ultimately lead to the formation of stars with masses ${\gt}500\, \mathrm{M}_{\rm \odot }$. Although we find that the mass scale set by jets is insensitive to most cloud parameters (i.e. surface density, virial parameter), it is strongly dependent on the momentum loading of the jets (which is poorly constrained by observations) as well as the temperature of the parent cloud, which predicts slightly larger IMF variations than observed. We conclude that protostellar jets play a vital role in setting the mass scale of stars, but additional physics are necessary to reproduce the observed IMF.

The Preservation of Super-Earths and the Emergence of Gas Giants after Their Progenitor Cores Have Entered the Pebble-isolation Phase

Abstract 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.

The imprint of the protoplanetary disc in the accretion of super-Earth envelopes

ABSTRACT Super-Earths are by far the most dominant type of exoplanet, yet their formation is still not well understood. In particular, planet formation models predict that many of them should have accreted enough gas to become gas giants. Here, we examine the role of the protoplanetary disc in the cooling and contraction of the protoplanetary envelope. In particular, we investigate the effects of (1) the thermal state of the disc as set by the relative size of heating by accretion or irradiation, and whether its energy is transported by radiation or convection, and (2) advection of entropy into the outer envelope by disc flows that penetrate the Hill sphere, as found in 3D global simulations. We find that, at 5 and 1 au, this flow at the level reported in the non-isothermal simulations where it penetrates only to ∼0.3 times the Hill radius has little effect on the cooling rate since most of the envelope mass is concentrated close to the core, and far from the flow. On the other hand, at 0.1 au, the envelope quickly becomes fully radiative, nearly isothermal, and thus cannot cool down, stalling gas accretion. This effect is significantly more pronounced in convective discs, leading to envelope mass orders of magnitude lower. Entropy advection at 0.1 au in either radiative or convective discs could therefore explain why super-Earths failed to undergo runaway accretion. These results highlight the importance of the conditions and energy transport in the protoplanetary disc for the accretion of planetary envelopes.

Deflating Super-puffs: Impact of Photochemical Hazes on the Observed Mass–Radius Relationship of Low-mass Planets

Abstract The observed mass–radius relationship of low-mass planets informs our understanding of their composition and evolution. Recent discoveries of low-mass, large-radius objects (“super-puffs”) have challenged theories of planet formation and atmospheric loss, as their high inferred gas masses make them vulnerable to runaway accretion and hydrodynamic escape. Here we propose that high-altitude photochemical hazes could enhance the observed radii of low-mass planets and explain the nature of super-puffs. We construct model atmospheres in radiative-convective equilibrium and compute rates of atmospheric escape and haze distributions, taking into account haze coagulation, sedimentation, diffusion, and advection by an outflow wind. We develop mass–radius diagrams that include atmospheric lifetimes and haze opacity, which is enhanced by the outflow, such that young (∼0.1–1 Gyr), warm (T eq ≥ 500 K), low-mass objects (M c < 4 M ⊕) should experience the most apparent radius enhancement due to hazes, reaching factors of three. This reconciles the densities and ages of the most extreme super-puffs. For Kepler-51b, the inclusion of hazes reduces its inferred gas mass fraction to <10%, similar to that of planets on the large-radius side of the sub-Neptune radius gap. This suggests that Kepler-51b may be evolving toward that population and that some warm sub-Neptunes may have evolved from super-puffs. Hazes also render transmission spectra of super-puffs and sub-Neptunes featureless, consistent with recent measurements. Our hypothesis can be tested by future observations of super-puffs’ transmission spectra at mid-infrared wavelengths, where we predict that the planet radius will be half of that observed in the near-infrared.

How planets grow by pebble accretion

Context. Proto-planets embedded in their natal disks acquire hot envelopes as they grow and accrete solids. This ensures that the material they accrete – pebbles, as well as (small) planetesimals – will vaporize to enrich their atmospheres. Enrichment modifies an envelope’s structure and significantly alters its further evolution. Aims. Our aim is to describe the formation of planets with polluted envelopes from the moment that impactors begin to sublimate to beyond the disk’s eventual dissipation. Methods. We constructed an analytical interior structure model, characterized by a hot and uniformly mixed high-Z vapor layer surrounding the core, located below the usual unpolluted radiative-convective regions. Our model assumes an ideal equation of state and focuses on identifying trends rather than precise calculations. The expressions we derived are applicable to all single-species pollutants, but we used SiO2 to visualize our results. Results. The evolution of planets with uniformly mixed polluted envelopes follows four potential phases. Initially, the central core grows directly through impacts and rainout until the envelope becomes hot enough to vaporize and absorb all incoming solids. We find that a planet reaches runaway accretion when the sum of its core and vapor mass exceeds a value that we refer to as the critical metal mass – a criterion that supersedes the traditional critical core mass. The critical metal mass scales positively with both the pollutant’s evaporation temperature and with the planet’s core mass. Hence, planets at shorter orbital separations require the accretion of more solids to reach runaway as they accrete less volatile materials. If the solids accretion rate dries up, we identify the decline of the mean molecular weight – dilution – as a mechanism to limit gas accretion during a polluted planet’s embedded cooling phase. When the disk ultimately dissipates, the envelope’s inner temperature declines and its vapor eventually rains out, augmenting the mass of the core. The energy release that accompanies this does not result in significant mass-loss, as it only occurs after the planet has substantially contracted.

K2-19b and c are in a 3:2 Commensurability but out of Resonance: A Challenge to Planet Assembly by Convergent Migration

Abstract K2-19b and c were among the first planets discovered by NASA’s K2 mission and together stand in stark contrast with the physical and orbital properties of the solar system planets. The planets are between the size of Uranus and Saturn at 7.0 ± 0.2 and 4.1 ± 0.2 , respectively, and reside a mere 0.1% outside the nominal 3:2 mean-motion resonance. They represent a different outcome of the planet formation process than the solar system, as well as the vast majority of known exoplanets. We measured the physical and orbital properties of these planets using photometry from K2, Spitzer, and ground-based telescopes, along with radial velocities from Keck/HIRES. Through a joint photodynamical model, we found that the planets have moderate eccentricities of e ≈ 0.20 and well-aligned apsides Δϖ ≈ 0°. The planets occupy a strictly nonresonant configuration: the resonant angles circulate rather than librate. This defies the predictions of standard formation pathways that invoke convergent or divergent migration, both of which predict Δϖ ≈ 180° and eccentricities of a few percent or less. We measured masses of M p,b = 32.4 ± 1.7 and M p,c = 10.8 ± 0.6 . Our measurements, with 5% fractional uncertainties, are among the most precise of any sub-Jovian exoplanet. Mass and size reflect a planet’s core/envelope structure. Despite having a relatively massive core of ≈ 15 , K2-19b is envelope-rich, with an envelope mass fraction of roughly 50%. This planet poses a challenge to standard models of core-nucleated accretion, which predict that cores ≳10 will quickly accrete gas and trigger runaway accretion when the envelope mass exceeds that of the core.

Global 3D radiation-hydrodynamic simulations of gas accretion: Opacity-dependent growth of Saturn-mass planets

The full spatial structure and temporal evolution of the accretion flow into the envelopes of growing gas giants in their nascent discs is only accessible in simulations. Such simulations are constrained in their approach of computing the formation of gas giants by dimensionality, resolution, consideration of self-gravity, energy treatment and the adopted opacity law. Our study explores how a number of these parameters affect the measured accretion rate of a Saturn-mass planet. We present a global 3D radiative hydrodynamics framework using the FARGOCA-code. The planet is represented by a gravitational potential with a smoothing length at the location of the planet. No mass or energy sink is used; instead luminosity and gas accretion rates are self-consistently computed. We find that the gravitational smoothing length must be resolved by at least ten grid cells to obtain converged measurements of the gas accretion rates. Secondly, we find gas accretion rates into planetary envelopes that are compatible with previous studies, and continue to explain those via the structure of our planetary envelopes and their luminosities. Our measured gas accretion rates are formally in the stage of Kelvin–Helmholtz contraction due to the modest entropy loss that can be obtained over the simulation timescale, but our accretion rates are compatible with those expected during late run-away accretion. Our detailed simulations of the gas flow into the envelope of a Saturn-mass planet provide a framework for understanding the general problem of gas accretion during planet formation and highlight circulation features that develop inside the planetary envelopes. Those circulation features feedback into the envelope energetics and can have further implications for transporting dust into the inner regions of the envelope.

Self-gravitating planetary envelopes and the core-nucleated instability

Abstract Planet formation scenarios can be constrained by the ratio of the gaseous envelope mass relative to the solid core mass in the observed exoplanet populations. One-dimensional calculations find a critical (maximal) core mass for quasi-static envelopes to exist, suggesting that envelopes around more massive cores should collapse due to a ‘core-nucleated’ instability. We study self-gravitating planetary envelopes via hydrodynamic simulations, progressively increasing the dimensionality of the problem. We characterize the core-nucleated instability and its non-linear evolution into runaway gas accretion in one-dimensional spherical envelopes. We show that rotationally supported envelopes can enter a runaway accretion regime via polar shocks in a two-dimensional axisymmetric model. This picture remains valid for high-mass cores in three dimensions, where the gas gravity mainly adds up to the core gravity and enhances the mass accretion rate of the planet in time. We relate the core-nucleated instability to the absence of equilibrium connecting the planet to its parent disc and discuss its relevance for massive planet formation.