skip to main content

Title: Global 3D radiation hydrodynamic simulations of proto-Jupiter’s convective envelope
ABSTRACT The core accretion model of giant planet formation has been challenged by the discovery of recycling flows between the planetary envelope and the disc that can slow or stall envelope accretion. We carry out 3D radiation hydrodynamic simulations with an updated opacity compilation to model the proto-Jupiter’s envelope. To isolate the 3D effects of convection and recycling, we simulate both isolated spherical envelopes and envelopes embedded in discs. The envelopes are heated at given rates to achieve steady states, enabling comparisons with 1D models. We vary envelope properties to obtain both radiative and convective solutions. Using a passive scalar, we observe significant mass recycling on the orbital time-scale. For a radiative envelope, recycling can only penetrate from the disc surface until ∼0.1–0.2 planetary Hill radii, while for a convective envelope, the convective motion can ‘dredge up’ the deeper part of the envelope so that the entire convective envelope is recycled efficiently. This recycling, however, has only limited effects on the envelopes’ thermal structure. The radiative envelope embedded in the disc has identical structure as the isolated envelope. The convective envelope has a slightly higher density when it is embedded in the disc. We introduce a modified 1D approach which more » can fully reproduce our 3D simulations. With our updated opacity and 1D model, we recompute Jupiter’s envelope accretion with a 10 M⊕ core, and the time-scale to runaway accretion is shorter than the disc lifetime as in prior studies. Finally, we discuss the implications of the efficient recycling on the observed chemical abundances of the planetary atmosphere (especially for super-Earths and mini-Neptunes). « less
; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Page Range or eLocation-ID:
453 to 474
Sponsoring Org:
National Science Foundation
More Like this

    Despite many methods developed to find young massive planets in protoplanetary discs, it is challenging to directly detect low-mass planets that are embedded in discs. On the other hand, the core-accretion theory suggests that there could be a large population of embedded low-mass young planets at the Kelvin-Helmholtz (KH) contraction phase. We adopt both 1D models and 3D simulations to calculate the envelopes around low-mass cores (several to tens of M⊕) with different luminosities, and derive their thermal fluxes at radio wavelengths. We find that, when the background disc is optically thin at radio wavelengths, radio observations can see through the disc and probe the denser envelope within the planet’s Hill sphere. When the optically thin disc is observed with the resolution reaching one disc scale height, the radio thermal flux from the planetary envelope around a 10 M⊕ core is more than 10 per cent higher than the flux from the background disc. The emitting region can be extended and elongated. Finally, our model suggests that the au-scale clump at 52 au in the TW Hydrae disc revealed by ALMA is consistent with the envelope of an embedded 10–20 M⊕ planet, which can explain the detected flux, the spectral index dip,more »and the tentative spirals. The observation is also consistent with the planet undergoing pebble accretion. Future ALMA and ngVLA observations may directly reveal more such low-mass planets, enabling us to study core growth and even reconstruct the planet formation history using the embedded ‘protoplanet’ population.

    « less
  2. Abstract

    UsingAthena++, we perform 3D radiation-hydrodynamic calculations of the radiative breakout of the shock wave in the outer envelope of a red supergiant (RSG) that has suffered core collapse and will become a Type IIP supernova. The intrinsically 3D structure of the fully convective RSG envelope yields key differences in the brightness and duration of the shock breakout (SBO) from that predicted in a 1D stellar model. First, the lower-density “halo” of material outside of the traditional photosphere in 3D models leads to a shock breakout at lower densities than 1D models. This would prolong the duration of the shock breakout flash at any given location on the surface to ≈1–2 hr. However, we find that the even larger impact is the intrinsically 3D effect associated with large-scale fluctuations in density that cause the shock to break out at different radii at different times. This substantially prolongs the SBO duration to ≈3–6 hr and implies a diversity of radiative temperatures, as different patches across the stellar surface are at different stages of their radiative breakout and cooling at any given time. These predicted durations are in better agreement with existing observations of SBO. The longer durations lower the predicted luminositiesmore »by a factor of 3–10 (Lbol∼ 1044erg s−1), and we derive the new scalings of brightness and duration with explosion energies and stellar properties. These intrinsically 3D properties eliminate the possibility of using observed rise times to measure the stellar radius via light-travel time effects.

    « less
  3. ABSTRACT (Sub)millimetre dust opacities are required for converting the observable dust continuum emission to the mass, but their values have long been uncertain, especially in discs around young stellar objects. We propose a method to constrain the opacity κν in edge-on discs from a characteristic optical depth τ0,ν, the density ρ0, and radius R0 at the disc outer edge through κν = τ0,ν/(ρ0R0), where τ0,ν is inferred from the shape of the observed flux along the major axis, ρ0 from gravitational stability considerations, and R0 from direct imaging. We applied the 1D semi-analytical model to the embedded, Class 0, HH 212 disc, which has high-resolution data in Atacama Large Millimetre/submillimetre Array (ALMA) bands 9, 7, 6, and 3 and Very Large Array Ka band (λ = 0.43, 0.85, 1.3, 2.9, and 9.1 mm). The modelling is extended to 2D through RADMC-3D radiative transfer calculations. We find a dust opacity of κν ≈ 1.9 × 10−2, 1.3 × 10−2, and 4.9 × 10−3 cm2 g−1 of gas and dust for ALMA bands 7, 6, and 3, respectively, with uncertainties dependent on the adopted stellar mass. The inferred opacities lend support to the widely used prescription κλ = 2.3 × 10−2(1.3mm/λ) cm2 g−1 . We inferred a temperature of ∼45 K at the disc outer edge that increasesmore »radially inwards. It is well above the sublimation temperatures of ices such as CO and N2, which supports the notion that the disc chemistry cannot be completely inherited from the protostellar envelope.« less

    Volume complete sky surveys provide evidence for a binary origin for the formation of isolated white dwarfs with magnetic fields in excess of a MegaGauss. Interestingly, not a single high-field magnetic white dwarf has been found in a detached system, suggesting that if the progenitors are indeed binaries, the companion must be removed or merge during formation. An origin scenario consistent with observations involves the engulfment, inspiral, and subsequent tidal disruption of a low-mass companion in the interior of a giant star during a common envelope phase. Material from the shredded companion forms a cold accretion disc embedded in the hot ambient around the proto-white dwarf. Entrainment of hot material may evaporate the disc before it can sufficiently amplify the magnetic field, which typically requires at least a few orbits of the disc. Using three-dimensional hydrodynamic simulations of accretion discs with masses between 1 and 10 times the mass of Jupiter inside the core of an Asymptotic Giant Branch star, we find that the discs survive for at least 10 orbits (and likely for 100 orbits), sufficient for strong magnetic fields to develop.


    During the core collapse of massive stars that do not undergo a canonical energetic explosion, some of the hydrogen envelope of a red supergiant (RSG) progenitor may infall on to the newborn black hole (BH). Within the athena++ framework, we perform 3D, hydrodynamical simulations of idealized models of supergiant convection and collapse in order to assess whether the infall of the convective envelope can give rise to rotationally supported material, even if the star has zero angular momentum overall. Our dimension-less, polytropic models are applicable to the optically thick hydrogen envelope of non-rotating RSGs and cover a factor of 20 in stellar radius. At all radii, the specific angular momentum due to random convective flows implies associated circularization radii of 10–1500 times the innermost stable circular orbit of the BH. During collapse, the angular momentum vector of the convective flows is approximately conserved and is slowly varying on the time-scale relevant to forming discs at small radii. Our results indicate that otherwise failed explosions of RSGs lead to the formation of rotationally supported flows that are capable of driving outflows to large radii and powering observable transients. When the BH is able to accrete most of the hydrogen envelope, the finalmore »BH spin parameter is ∼ 0.5, even though the star is non-rotating. For fractional accretion of the envelope, the spin parameter is generally lower and never exceeds 0.8. We discuss the implications of our results for transients produced by RSG collapse to a black hole.

    « less