The 3D dust and opacity distribution of protoplanets in multi-fluid global simulations

TL;DR

This paper presents a 3D simulation of dust distribution and opacities in protoplanet envelopes, revealing significant latitudinal gradients and their impact on planetary formation models.

Contribution

It introduces a comprehensive 3D model of dust dynamics in protoplanet envelopes, highlighting the importance of dust settling and asymmetries on opacity distributions.

Findings

Dust settling creates strong latitudinal opacity gradients.

Spiral wakes contribute to dust asymmetry and opacity variation.

Opacity can vary by over two orders of magnitude between mid-plane and poles.

Abstract

The abundance and distribution of solids inside the Hill sphere are central to our understanding of the giant planet dichotomy. Here, we present a three-dimensional characterization of the dust density, mass flux, and mean opacities in the envelope of sub-thermal and super-thermal mass planets. We simulate the dynamics of multiple dust species in a global protoplanetary disk model accounting for dust feedback. We find that the meridional flows do not effectively stir dust grains at scales of the Bondi sphere. Thus the dust-settling driven by the stellar gravitational potential sets the latitudinal dust density gradient within the planet envelope. Not only does the planet's potential enhance this gradient, but also the spiral wakes serve as another source of asymmetry. These asymmetries substantially alter the inferred mean Rosseland and Planck opacities. In cases with the moderate-to-strong dust settling, the opacity gradient can range from a few percent to more than two orders of magnitude between the mid-plane and the polar regions of the Bondi sphere. Finally, we show that this strong latitudinal opacity gradient can introduce a transition between optically thick and thin regimes at the scales of the planet envelope. We suggest that this transition is likely to occur when the equilibrium scale height of hundred-micron-sized particles is smaller than the Hill radius of the forming planet. This work calls into question the adoption of a constant opacity derived from well-mixed distributions and demonstrates the need for global radiation hydrodynamics models of giant planet formation which account for dust dynamics.