A Systematic Study of Planetary Envelope Growth with 3D Radiation-Hydrodynamics Simulations

Abstract: In the core accretion model of planet formation, envelope cooling regulates the accretion of material and ultimately sets the timescale to form a giant planet. Given the diversity of planet-forming environments, opacity uncertainties, and the advective transport of energy by 3-dimensional recycling flows, it is unclear whether 1D models can adequately describe envelope structure and accretion in all regimes. Even in 3D models, it is unclear whether approximate radiative transfer methods sufficiently model envelope cooling particularly at the planetary photosphere. To address these uncertainties, we present a suite of 3D radiation hydrodynamics simulations employing methods that directly solve the transfer equation. We perform a parameter space study, formulated in terms of dimensionless parameters, for a variety of envelope optical depths and cooling times. We find that the thermodynamic structure of the envelope ranges from adiabatic to isothermal based on the cooling time and by extension, the background disk temperature and density. Our models show general agreement with 1D static calculations, suggesting a limited role of recycling flows in determining envelope structure. By adopting a dimensionless framework, these models can be applied to a wide range of formation conditions and assumed opacities. In particular, we dimensionalize them to the case of a super-Earth and proto-Jupiter and place upper limits on the 3D mass accretion rates prior to runaway growth. Finally, we evaluate the fidelity of approximate radiative transfer methods and find that even in the most challenging cases, more approximate methods are sufficiently accurate and worth their savings in computational cost.