Circumplanetary discs (CPDs) form around young gas giants and are thought to be the sites of moon formation as well as an intermediate reservoir of gas that feeds the growth of the gas giant. How the physical properties of such CPDs are affected by the planetary mass and the overall opacity is relatively poorly understood. In order to clarify this, we use the global radiation hydrodynamics code FARGOCA, with a grid structure that allows resolving the planetary gravitational potential sufficiently well for a CPD to form. We then study the gas flows and density/temperature structures that emerge as a function of planet mass, opacity and potential depth. Our results indicate interesting structure formation for Jupiter-mass planets at low opacities, which we subsequently analyse in detail. Using an opacity level that is 100 times lower than that of ISM dust, our Jupiter-mass protoplanet features an envelope that is sufficiently cold for a CPD to form, and a free-fall region separating the CPD and the circumstellar disc emerges. Interestingly, this free-fall region appears to be a result of supersonic erosion of outer envelope material, as opposed to the static structure formation that one would expect at low opacities. Our analysis reveals that the planetary spiral arms seem to pose a significant pressure barrier that needs to be overcome through radiative cooling in order for gas to free-fall onto the CPD. The circulation inside the CPD is near-keplerian and modified by the presence of CPD spiral arms. For high opacities we recover results from the literature, finding an essentially featureless hot envelope. With this work, we demonstrate the first simulation and analysis of a complete detachment process of a protoplanet from its parent disc in a 3D radiation hydrodynamics setting.

The full spatial structure and temporal evolution of the accretion flow into the envelopes of growing gas giants in their nascent discs is only accessible in simulations. Such simulations are constrained in their approach of computing the formation of gas giants by dimensionality, resolution, consideration of self-gravity, energy treatment and the adopted opacity law. Our study explores how a number of these parameters affect that measured accretion rate of a Saturn-mass planet. We present a global 3D radiative hydrodynamics framework using the FARGOCA-code. The planet is represented by a gravitational potential with a smoothing length at the location of the planet. No mass or energy sink is used, instead luminosity and gas accretion rates are self-consistently computed. We find that the gravitational smoothing length must be resolved by at least 10 grid cells to obtain converged measurements of the gas accretion rates. Secondly, we find gas accretion rates into planetary envelopes compatible with previous studies, and continue to explain those via the structure of our planetary envelopes and their luminosities. Our measured gas accretion rates are formally in the stage of Kelvin-Helmholtz contraction due to the modest entropy loss that can be obtained over the simulation time-scale, but our accretion rates are compatible with those expected during late run-away accretion. Our detailed simulations of the gas flow into the envelope of a Saturn-mass planet provides a framework for understanding the general problem of gas accretion during planet formation and highlights circulation features that develop inside the planetary envelopes. Those circulation features feedback into the envelope energetics and can have further implications for transporting dust into the inner regions of the envelope.