NBIA logo


Horizon 2020 - ERC Starting Grant

"First-principles global MHD disc simulations:
Defining planet-forming environments
in early solar systems"

The aim of my research project is to create the most realistic computer simulations of protoplanetary discs of gas and dust, thus defining the environment that shapes the early development of planetary systems.

Project overview

In their role as planet nurseries, protoplanetary discs are of key interest to planet formation theory. Their dynamical, radiative and thermodynamic properties critically define the environment for embedded solids: dust grains, pebbles and planetesimals. In short, the building blocks of planet formation. The discs' dynamics and structure in turn depend critically on the influence of magnetic fields that couple to tenuously ionised and low-density regions. Being comparatively cold and dense, the ionisation state of the disc plasma is dominated by external far-UV, X-Ray, and cosmic-ray radiation, leading to a layered vertical structure - with turbulent surface layers and a magnetically-decoupled midplane. This classic dead-zone picture is now turned upside-down by previously ignored microphysical effects.

back to top

Highlight: Magnetocentrifugal winds

Protoplanetary disks accrete onto their central T Tauri star via magnetic stresses. When the effect of ambipolar diffusion is included, and in the presence of a vertical magnetic field, the disk remains laminar in its planet-forming inner region, and a magnetocentrifugal disk wind forms that provides an important mechanism for removing angular momentum. We have performed global MHD simulations where the time-dependent gas-phase electron and ion fractions are computed under FUV and X-ray ionization with a simplified recombination chemistry.

Radiative-magnetic wind model: Detail from a proof-of-concept global MHD simulation of a PPD with ambipolar diffusion, radiative transfer and stellar irradiation, showing the magnetic field strength (colour), velocity arrows (black) and magnetic field lines (white), and iso-contours of the radiation temperature (grey).

To investigate whether the mass loading of the wind is potentially affected by the limited vertical extent of our existing simulations, we attempt to develop a model of a realistic disk atmosphere. To this end, by accounting for stellar irradiation and diffuse reprocessing of radiation, we aim at improving our models towards more realistic thermodynamic properties.

related publications

back to top

Highlight: Jupiter acting as a dust barrier

Aerodynamic theory predicts that dust grains in protoplanetary disks will drift radially inward on comparatively short timescales. In this context, it has long been known that the presence of a gap opened by a planet can alter the dust dynamics significantly. Lead by PhD student Philipp Weber, we have carried out a systematic study employing long-term numerical simulations aimed at characterizing the critical particle-size for retention outside a gap as a function of particle size and for various key parameters defining the protoplanetary disk model.

Dust stream around an embedded planet: Surface density of ten micron-sized dust that was initialized in a single azimuthal grid cell outside the planet's location. The panels show the distribution of dust after 8, and 35 orbits, respectively, in a frame that is co-rotating with the planet, and the circle illustrates the planet's sphere of influence.

In agreement with previous work, we find that the permeability of the gap depends both on dust dynamical properties and the gas disk structure: while small dust follows the viscously accreting gas through the gap, dust grains approaching a critical size are progressively filtered out.

Gap permeability as a function of particle size: The figure shows radial profiles of the dust surface density, Σd normalized to the initial gas surface density, Σg,0. Particles larger in size than about a tenth of a centimeter begin to aerodynamically decouple from the gas flow, and, as a result, become significantly depleted in the inner system.

We introduce and compute a depletion factor (see figure above) that enables us to quantify the way in which higher viscosity, smaller planet mass, or a more massive disk can shift this critical size to larger values. Our results indicate that gap-opening planets may act to deplete the inner reaches of protoplanetary disks of large dust grains, potentially limiting the accretion of solids onto forming terrestrial planets.

related publications

back to top

Highlight: Self-organization with Hall-MHD

Imaging of dust continuum emitted from disks around nearby protostars reveals diverse substructure. In recent years, theoretical efforts have been intensified to investigate in how far the intrinsic dynamics of protoplanetary disks (PPDs) can lead to such features. Turbulence in the realm of non-ideal magnetohydrodynamics (MHD) is one candidate for explaining the generation of zonal flows which can lead to local dust enhancements.

Spontaneous self-organization of magnetic flux: In this snapshot from a Hall-MHD simulation, the turbulent fluctuations created by the magneto-rotational instability have undergone an extended phase of self-organization leading to zonal bands and vortices. In our case, four bands of vertical magnetic flux are obtained.

In an effort led by PhD student Leonardo Krapp, we consider combinations of vertical and azimuthal initial net flux and perform 3D non-ideal MHD simulations aimed at studying self-organization induced by the Hall effect in turbulent PPDs. We include dust grains, treated in the fluid approximation, in order to study their evolution subject to the emerging zonal flows. In the regime of a dominant Hall effect, we robustly obtain large-scale organized concentrations in the vertical field that remain stable for many orbits.

Accumulation of dust into rings: Grains of various sizes (top three panels) concentrate at locations in the disk, where the flow deviates from its sub-Keplerian equilibrium profile (4th panel), related to the accumulation of magnetic flux (5th panel). The aerodynamic coupling of the particles is characterized by their so-called Stokes number, St, describing their relative coupling to the gas.

Our research has, moreover, found that including a moderately strong net-azimuthal magnetic flux can significantly alter the dynamics, partially preventing the self-organization of zonal flows. However, in all cases, we demonstrate that the emerging features are capable of accumulating dust grains for a range of Stokes numbers.

related publications

back to top

This project has received funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No 638596).