Revising the Giant Planet Mass-Metallicity Relation: Deciphering the Formation Sequence of Giant Planets

The rate at which giant planets accumulate solids and gas is a critical component of planet formation models, yet it is extremely challenging to predict from first principles. Characterizing the heavy element (everything other than hydrogen and helium) content of giant planets provides important clues about their provenance. Using thermal evolution models with updated H-He EOS and atmospheric boundary condition that varies with envelope metallicity, we quantify the bulk heavy element content of 147 warm ($< 1000$ K) giant planets with well-measured masses and radii, more than tripling the sample size studied in Thorngren et al. 2016. These measurements reveal that the population's heavy element mass follows the relation $M_{\rm Z} = M_{\rm core} + f_Z (M_{\rm p} - M_{\rm core})$, with $M_{\rm core} = 14.7^{+1.8}_{-1.6}$ Earth masses (M$_\oplus$), $f_Z = 0.09 \pm 0.01$, and an astrophysical scatter of $0.66 \pm 0.08 \times M_Z$. The classical core-accretion scenario ($Z_{\rm p} = 1$ at 10 M$_\oplus$ and $Z_{\rm p} = 0.5$ at 20 M$_\oplus$) is inconsistent with the population. At low planet masses ($<< 150$ M$_\oplus$), $M_{\rm Z} \sim M_{\rm core}$ and as a result, $Z_{\rm p} = M_{\rm Z} / M_{\rm p}$ declines linearly with $M_{\rm p}$. However, bulk metallicity does not continue to decline with planet mass and instead flattens out at $f_Z \sim 0.09$ ($\sim 7 \times$ solar metallicity). When normalized by stellar metallicity, $Z_{\rm p} / Z_\star$ flattens out at $3.3 \pm 0.5$ at high planet masses. This explicitly shows that giant planets continue to accrete material enriched in heavy elements during the gas accretion phase.