The vertical temperature structure of a protoplanetary disk bears on several processes relevant to planet formation, such as gas and dust grain chemistry, ice lines and convection. The temperature profile is controlled by irradiation from the central star and by any internal source of heat such as might arise from gas accretion. We investigate the heat and angular momentum transport generated by the resistive dissipation of magnetic fields in laminar disks. We use local one-dimensional simulations to obtain vertical temperature profiles for typical conditions in the inner disk (0.5 to 4 au). Using simple assumptions for the gas ionization and opacity, the heating and cooling rates are computed self-consistently in the framework of radiative non-ideal magnetohydrodynamics. We characterize steady solutions that are symmetric about the midplane and which may be associated with saturated Hall-shear unstable modes. We also examine the dissipation of electric currents driven by global accretion-ejection structures. In both cases we obtain significant heating for a sufficiently high opacity. Strong magnetic fields can induce an order-unity temperature increase in the disk midplane, a convectively unstable entropy profile, and a surface emissivity equivalent to a viscous heating of $\alpha \sim 10^{-2}$. These results show how magnetic fields may drive efficient accretion and heating in weakly ionized disks where turbulence might be inefficient, at least for a range of radii and ages of the disk.

Planet formation scenarios can be constrained by the ratio of the gaseous envelope mass relative to the solid core mass in the observed exoplanet populations. One-dimensional calculations find a critical (maximal) core mass for quasi-static envelopes to exist, suggesting that envelopes around more massive cores should collapse due to a `core-nucleated' instability. We study self-gravitating planetary envelopes via hydrodynamic simulations, progressively increasing the dimensionality of the problem. We characterize the core-nucleated instability and its non-linear evolution into runaway gas accretion in one-dimensional spherical envelopes. We show that rotationally-supported envelopes can enter a runaway accretion regime via polar shocks in a two-dimensional axisymmetric model. This picture remains valid for high-mass cores in three dimensions, where the gas gravity mainly adds up to the core gravity and enhances the mass accretion rate of the planet in time. We relate the core-nucleated instability to the absence of equilibrium connecting the planet to its parent disk and discuss its relevance for massive planet formation.

Massive planetary cores embedded in protoplanetary discs are believed to accrete extended atmospheres, providing a pathway to forming gas giants and gas-rich super-Earths. The properties of these atmospheres strongly depend on the nature of the coupling between the atmosphere and the surrounding disc. We examine the formation of gaseous envelopes around massive planetary cores via three-dimensional inviscid and isothermal hydrodynamic simulations. We focus the changes in the envelope properties as the core mass varies from low (sub-thermal) to high (super-thermal) values, a regime relevant to close-in super-Earths. We show that global envelope properties such as the amount of rotational support or turbulent mixing are mostly sensitive to the ratio of the Bondi radius of the core to its physical size. High-mass cores are fed by supersonic inflows arriving along the polar axis and shocking on the densest parts of the envelope, driving turbulence and mass accretion. Gas flows out of the core's Hill sphere in the equatorial plane, describing a global mass circulation through the envelope. The shell of shocked gas atop the core surface delimits regions of slow (inside) and fast (outside) material recycling by gas from the surrounding disc. While recycling hinders the runaway growth towards gas giants, the inner regions of protoplanetary atmospheres, more immune to mixing, may remain bound to the planet.

Measurements of exoplanetary masses and radii have revealed a population of massive super-Earths --- planets sufficiently large that, according to one dimensional models, they should have turned into gas giants. To better understand the origin of these objects, we carry out hydrodynamical simulations of planetary cores embedded in a nascent protoplanetary disk. In this first paper of a series, to gain intuition as well as to develop useful diagnostics, we focus on two-dimensional simulations of the flow around protoplanetary cores. We use the pluto code to study isothermal and adiabatic envelopes around cores of sub- to super-thermal masses, fully resolving the envelope properties down to the core surface. Owing to the conservation of vortensity, envelopes acquire a substantial degree of rotational support when the core mass increases beyond the thermal mass, suggesting a limited applicability of one-dimensional models for describing the envelope structure. The finite size of the core (relatively large for super-Earths) also controls the amount of rotational support in the entire envelope. Steady non-axisymmetric shocks develop in the supersonic envelopes of high-mass cores, triggering mass accretion and turbulent mixing in their interiors. We also examine the influence of the gas self-gravity on the envelope structure. Although it only weakly alters the properties of the envelopes, the gas gravity has significant effect on the properties of the density waves triggered by the core in the protoplanetary disk.

The planet-forming region of protoplanetary disks is cold, dense, and therefore weakly ionized. For this reason, magnetohydrodynamic (MHD) turbulence is thought to be mostly absent, and another mechanism has to be found to explain gas accretion. It has been proposed that magnetized winds, launched from the ionized disk surface, could drive accretion in the presence of a large-scale magnetic field. The efficiency and the impact of these surface winds on the disk structure is still highly uncertain. We present the first global simulations of a weakly ionized disk that exhibits large-scale magnetized winds. We also study the impact of self-organization, which was previously demonstrated only in non-stratified models. We perform numerical simulations of stratified disks with the PLUTO code. We compute the ionization fraction dynamically, and account for all three non-ideal MHD effects: ohmic and ambipolar diffusions, and the Hall drift. Simplified heating and cooling due to non-thermal radiation is also taken into account in the disk atmosphere. We find that disks can be accreting or not, depending on the configuration of the large-scale magnetic field. Magnetothermal winds, driven both by magnetic acceleration and heating of the atmosphere, are obtained in the accreting case. In some cases, these winds are asymmetric, ejecting predominantly on one side of the disk. The wind mass loss rate depends primarily on the average ratio of magnetic to thermal pressure in the disk midplane. The non-accreting case is characterized by a meridional circulation, with accretion layers at the disk surface and decretion in the midplane. Finally, we observe self-organization, resulting in axisymmetric rings of density and associated pressure "bumps". The underlying mechanism and its impact on observable structures are discussed.

Recent observations revealed organised structures in protoplanetary disks, such as axisymmetric rings or horseshoe concen- trations evocative of large-scale vortices. These structures are often interpreted as the result of planet-disc interactions. However, these disks are also known to be unstable to the magneto-rotational instability (MRI) which is believed to be one of the dominant angular momentum transport mechanism in these objects. It is therefore natural to ask if the MRI itself could produce these structures without invoking planets. The nonlinear evolution of the MRI is strongly affected by the low ionisation fraction in protoplanetary disks. The Hall effect in particular, which is dominant in dense and weakly ionised parts of these objects, has been shown to spontaneously drive self- organising flows in shearing box simulations. Here, we investigate the behaviour of global MRI-unstable disc models dominated by the Hall effect and characterise their dynamics. We perform 3D unstratified Hall-MHD simulations of keplerian disks for a broad range of Hall, ohmic and ambipolar Elsasser numbers. We confirm the transition from a turbulent to an organised state as the intensity of the Hall effect is increased. We observe the formation of zonal flows, their number depending on the available magnetic flux and on the intensity of the Hall effect. For intermediate Hall intensity, the flow self-organises into long-lived magnetised vortices. Neither the addition of a toroidal field nor ohmic or ambipolar diffusion drastically change this picture in the range of parameters we have explored. The ability of these structures to trap dust particles in this configuration is demonstrated. We conclude that Hall-MRI driven organisation is a plausible scenario which could explain some of the structures found in recent observations.