Another aspect is that [1]} calculate the extinguishing column density from a distance \(r=0.03\) –\(0.1~{R_{\textnormal {Hill}}}\)  outwards from the planet, due to the smoothing.
Given the rescaling, for PDS 70 b at 22 au this corresponds to \(r=80\) , 340, and \(\approx 200~R_{\mathrm {J}}\)  as a starting radius for the integration for the \(M_{\mathrm {p}}=0.48\) , 0.95, and 2.48 \(M_{\mathrm {J}}\)  simulations. Thus, their column densities are smaller than what higher-resolution simulations would yield.
Nevertheless,
the qualitative result
of a decreasing extinction with increasing mass should be robust.
