Large-scale vertical magnetic fields are believed to play a key role in the evolution of protoplanetary discs. Associated with non-ideal effects, such as ambipolar diffusion, they are known to launch a wind that could drive accretion in the outer part of the disc ($R> 1$ AU). They also potentially lead to self-organisation of the disc into large-scale axisymmetric structures, similar to the rings recently imaged by sub-millimetre or near-infrared instruments (ALMA and SPHERE). The aim of this paper is to investigate the mechanism behind the formation of these gaseous rings, but also to understand the dust dynamics and its emission in discs threaded by a large-scale magnetic field. To this end, we performed global magneto-hydrodynamics (MHD) axisymmetric simulations with ambipolar diffusion using a modified version of the PLUTO code. We explored different magnetisations with the midplane $\beta$ parameter ranging from $10^5$ to $10^3$ and included dust grains -- treated in the fluid approximation -- ranging from $100 \mu$m to 1 cm in size. We first show that the gaseous rings (associated with zonal flows) are tightly linked to the existence of MHD winds. Secondly, we find that millimetre-size dust is highly sedimented, with a typical scale height of 1 AU at $R=100$ AU for $\beta=10^4$, compatible with recent ALMA observations. We also show that these grains concentrate into pressure maxima associated with zonal flows, leading to the formation of dusty rings. Using the radiative transfer code MCFOST, we computed the dust emission and make predictions on the ring-gap contrast and the spectral index that one might observe with interferometers like ALMA.

Gravitational instability (GI) controls the dynamics of young massive protoplanetary discs. Apart from facilitating gas accretion on to the central protostar, it must also impact on the process of planet formation: directly through fragmentation, and indirectly through the turbulent concentration of small solids. To understand the latter process, it is essential to determine the dust dynamics in such a turbulent flow. For that purpose, we conduct a series of 3D shearing box simulations of coupled gas and dust, including the gas's self-gravity and scanning a range of Stokes numbers, from 0.001 to ~0.2. First, we show that the vertical settling of dust in the midplane is significantly impeded by gravitoturbulence, with the dust scale-height roughly 0.6 times the gas scale height for centimetre grains. This is a result of the strong vertical diffusion issuing from (a) small-scale inertial-wave turbulence feeding off the GI spiral waves and (b) the larger-scale vertical circulations that naturally accompany the spirals. Second, we show that at R=50 AU concentration events involving sub-metre particles and yielding order 1 dust to gas ratios are rare and last for less than an orbit. Moreover, dust concentration is less efficient in 3D than in 2D simulations. We conclude that GI is not especially prone to the turbulent accumulation of dust grains. Finally, the large dust scale-height measured in simulations could be, in the future, compared with that of edge-on discs seen by ALMA, thus aiding detection and characterisation of GI in real systems.

Rings and gaps have been observed in a wide range of protoplanetary discs, from young systems like HLTau to older discs like TW Hydra. Recent disc simulations have shown that magnetohydrodynamic (MHD) turbulence (in the ideal or non-ideal regime) can lead to the formation of rings and be an alternative to the embedded planets scenario. In this paper, we investigate how these ring form in this context and seek a generic formation process, taking into account the various dissipative regimes and magnetizations probed by the past simulations. We identify the existence of a linear and secular instability, driven by MHD winds, and giving birth to rings of gas having a width larger than the disc scale height. We show that the linear theory is able to make reliable predictions regarding the growth rates, ring/gap contrast and spacing, by comparing these predictions to a series of 2D (axisymmetric) and 3D MHD numerical simulations. In addition, we demonstrate that these rings can act as dust traps provided that the disc is sufficiently magnetised, with plasma beta lower than $10^4$. Given its robustness, the process identified in this paper could have important implications, not only for protoplanetary discs but also for a wide range of accreting systems threaded by large-scale magnetic fields.

The origin of large-scale and coherent magnetic fields in astrophysical discs is an important and long standing problem. Researchers commonly appeal to a turbulent dynamo, sustained by the magneto-rotational instability (MRI), to supply the large-scale field. But research over the last decade in particular has demonstrated that various non-ideal MHD effects can impede or extinguish the MRI, especially in protoplanetary disks. In this paper we propose a new scenario, by which the magnetic field is generated and sustained via the gravitational instability (GI). We use 3D stratified shearing box simulations to characterise the dynamo and find that it works at low magnetic Reynolds number (from unity to ~100) for a wide range of cooling times and boundary conditions. The process is kinematic, with a relatively fast growth rate ($< 0.1 \Omega$), and it shares some properties of mean field dynamos. The magnetic field is generated via the combination of differential rotation and spiral density waves, the latter providing compressible horizontal motions and large-scale vertical rolls. At greater magnetic Reynolds numbers the build up of large-scale field is diminished and instead small-scale structures emerge from the breakdown of twisted flux ropes. We propose that GI may be key to the dynamo engine not only in young protoplanetary discs but also in some AGN and galaxies.

MHD turbulence plays a crucial role in the dust dynamics of protoplanetary discs. It affects planet formation, vertical settling and is one possible origin of the large scale axisymmetric structures, such as rings, recently imaged by ALMA and SPHERE. Among the variety of MHD processes, the magnetorotational instability (MRI) has raised particular interest since it provides a source of turbulence and potentially organizes the flow into radial structures. However, the weak ionization of discs prevents the MRI from being excited beyond 1 AU. The strong sedimentation of millimetre dust measured in T-Tauri discs is also in contradiction with predictions based on ideal MRI turbulence. In this paper, we study the effects of non-ideal MHD and winds on the dynamics and sedimentation of dust grains. We consider a weakly ionized plasma subject to ambipolar diffusion characterizing the disc outer regions (>1 AU). For that, we perform numerical MHD simulations in the stratified shearing box, using the PLUTO code. Our simulations show that the mm-cm dust is contained vertically in a very thin layer, with typical heightscale ~0.4 AU at 30 AU, compatible with recent ALMA observations. Horizontally, the grains are trapped within pressure maxima induced by ambipolar diffusion, leading to the formation of dust rings. For micrometer grains, dust and gas scaleheights are similar. In this regime, the settling cannot be explained by a simple 1D diffusion theory but results from a large scale 2D circulation induced by both MHD winds and zonal flows. Overall, our results show that non-ideal MHD effects and their related winds play a major role in shaping the radial and vertical distribution of dust in protoplanetary discs. Leading to substantial accretion efficiency, non-ideal effects also a promising avenue to reconcile the low turbulent activity measured in discs with their relatively high accretion rates.

Spiral density waves dominate several facets of accretion disc dynamics --- planet-disc interactions and gravitational instability (GI) most prominently. Though they have been examined thoroughly in two-dimensional simulations, their vertical structures in the non-linear regime are somewhat unexplored. This neglect is unwarranted given that any strong vertical motions associated with these waves could profoundly impact dust dynamics, dust sedimentation, planet formation, and the emissivity of the disc surface. In this paper we combine linear calculations and shearing box simulations in order to investigate the vertical structure of spiral waves for various polytropic stratifications and wave amplitudes. For sub-adiabatic profiles we find that spiral waves develop a pair of counter-rotating poloidal rolls. Particularly strong in the nonlinear regime, these vortical structures issue from the baroclinicity supported by the background vertical entropy gradient. They are also intimately connected to the disk's g-modes which appear to interact nonlinearly with the density waves. Furthermore, we demonstrate that the poloidal rolls are ubiquitous in gravitoturbulence, emerging in the vicinity of GI spiral wakes, and potentially transporting grains off the disk midplane. Other than hindering sedimentation and planet formation, this phenomena may bear on observations of the disk's scattered infrared luminosity. The vortical features could also impact on the turbulent dynamo operating in young protoplanetary discs subject to GI, or possibly even galactic discs.

Though usually treated in isolation, the magnetorotational and gravitational instabilities (MRI and GI) may coincide at certain radii and evolutionary stages of protoplanetary discs and active galactic nuclei. Their mutual interactions could profoundly influence several important processes, such as accretion variability and outbursts, fragmentation and disc truncation, or large-scale magnetic field production. Direct numerical simulations of both instabilities are computationally challenging and remain relatively unexplored. In this paper, we aim to redress this neglect via a set of 3D vertically stratified shearing-box simulations, combining self-gravity and magnetic fields. We show that gravito-turbulence greatly weakens the zero-net-flux MRI. In the limit of efficient cooling (and thus enhanced GI), the MRI is completely suppressed, and yet strong magnetic fields are sustained by the gravitoturbulence. This turbulent `spiral wave' dynamo may have widespread application, especially in galactic discs. Finally, we present preliminary work showing that a strong net-vertical-flux revives the MRI and supports a magnetically dominated state, in which the GI is secondary.

Young protoplanetary discs and the outer radii of active galactic nucleii may be subject to gravitational instability and, as a consequence, fall into a `gravitoturbulent' state. While in this state, appreciable angular momentum can be transported. Alternatively, the gas may collapse into bound clumps, the progenitors of planets or stars. In this paper, we numerically characterize the properties of 3D gravitoturbulence, focussing especially on its dependence on numerical parameters (resolution, domain size) and its excitation of small-scale dynamics. Via a survey of vertically stratified shearing box simulations with PLUTO and RODEO, we find (a) evidence that certain gravitoturbulent properties are independent of horizontal box size only when the box is larger than $\simeq 40 H_0$, where $H_0$ is the height scale, (b) at high resolution, small-scale isotropic turbulence appears off the midplane around $z\simeq 0.5 -1 H_0$, and (c) this small-scale dynamics results from a parametric instability, involving the coupling of inertial waves with a large-scale axisymmetric epicyclic mode. This mode oscillates at a frequency close to $\Omega$ and is naturally excited by gravito-turbulence, via a nonlinear process to be determined. The small-scale turbulence we uncover has potential implications for a wide range of disc physics, e.g. turbulent saturation levels, fragmentation, turbulent mixing, and dust settling.

Jets, from the protostellar to the AGN context, have been extensively studied but their connection to the turbulent dynamics of the underlying accretion disc is poorly understood. Following a similar approach to Lesur et al. (2013), we examine the role of the magnetorotational instability (MRI) in the production and acceleration of outflows from discs. Via a suite of one-dimensional shearing-box simulations of stratified discs we show that magneto-centrifugal winds exhibit cyclic activity with a period of 10-20 Omega^{-1}, a few times the orbital period. The cycle seems to be more vigorous for strong vertical field; it is robust to the variation of relevant parameters and independent of numerical details. The convergence of these solutions (in particular the mass loss rate) with vertical box size is also studied. By considering a sequence of magnetohydrostatic equilibria and their stability, the periodic activity may be understood as the succession of the following phases: (a) a dominant MRI channel mode, (b) strong magnetic field generation, (c) consequent wind launching, and ultimately (d) vertical expulsion of the excess magnetic field by the expanding and accelerating gas associated with the wind. We discuss potential connections between this behaviour and observed time-variability in disk-jet systems.

In Keplerian accretion disks, turbulence and magnetic fields may be jointly excited through a subcritical dynamo process involving the magnetorotational instability (MRI). High-resolution simulations exhibit a tendency towards statistical self-organization of MRI dynamo turbulence into large-scale cyclic dynamics. Understanding the physical origin of these structures, and whether they can be sustained and transport angular momentum efficiently in astrophysical conditions, represents a significant theoretical challenge. The discovery of simple periodic nonlinear MRI dynamo solutions has recently proven useful in this respect, and has notably served to highlight the role of turbulent magnetic diffusion in the seeming decay of the dynamics at low magnetic Prandtl number Pm (magnetic diffusivity larger than viscosity), a common regime in accretion disks. The connection between these simple structures and the statistical organization reported in turbulent simulations remained elusive, though. Here, we report the numerical discovery in moderate aspect ratio Keplerian shearing boxes of new periodic, incompressible, three-dimensional nonlinear MRI dynamo solutions with a larger dynamical complexity reminiscent of such simulations. These "chimera" cycles are characterized by multiple MRI-unstable dynamical stages, but their basic physical principles of self-sustainment are nevertheless identical to those of simpler cycles found in azimuthally elongated boxes. In particular, we find that they are not sustained at low Pm either due to subcritical turbulent magnetic diffusion. These solutions offer a new perspective into the transition from laminar to turbulent instability-driven dynamos, and may prove useful to devise improved statistical models of turbulent accretion disk dynamos.

Gravitational instability (GI) features in several aspects of protostellar disk evolution, most notably in angular momentum transport, fragmentation, and the outbursts exemplified by FU Ori and EX Lupi systems. The outer regions of protostellar discs may also be coupled to magnetic fields, which could then modify the development of GI. To understand the basic elements of their interaction, we perform local 2D ideal and resistive MHD simulations with an imposed toroidal field. In the regime of moderate plasma beta, we find that the system supports a hot gravito-turbulent state, characterised by considerable magnetic energy and stress and a surprisingly large Toomre parameter $Q~10$. This result has potential implications for disk structure, vertical thickness, ionisation, etc. Our simulations also reveal the existence of long-lived and dense `magnetic islands' or plasmoids. Lastly, we find that the presence of a magnetic field has little impact on the fragmentation criterion of the disk. Though our focus is on protostellar disks, some of our results may be relevant for the outer radii of AGN.

The magnetorotational (MRI) dynamo has long been considered one of the possible drivers of turbulent angular momentum transport in astrophysical accretion disks. However, various numerical results suggest that this dynamo may be difficult to excite in the astrophysically relevant regime of magnetic Prandtl number (Pm) significantly smaller than unity, for reasons currently not well understood. The aim of this article is to present the first results of an ongoing numerical investigation of the role of both linear and nonlinear dissipative effects in this problem. Combining a parametric exploration and an energy analysis of incompressible nonlinear MRI dynamo cycles representative of the transitional dynamics in large aspect ratio shearing boxes, we find that turbulent magnetic diffusion makes the excitation and sustainment of this dynamo at moderate magnetic Reynolds number (Rm) increasingly difficult for decreasing Pm. This results in an increase in the critical Rm of the dynamo for increasing kinematic Reynolds number (Re), in agreement with earlier numerical results. Given its very generic nature, we argue that turbulent magnetic diffusion could be an important determinant of MRI dynamo excitation in disks, and may also limit the efficiency of angular momentum transport by MRI turbulence in low Pm regimes.

Magnetorotational dynamo action in Keplerian shear flow is a three-dimensional, nonlinear magnetohydrodynamic process whose study is relevant to the understanding of accretion and magnetic field generation in astrophysics. Transition to this form of dynamo is subcritical and shares many characteristics of transition to turbulence in non-rotating hydrodynamic shear flows. This suggests that these different fluid systems become active through similar generic bifurcation mechanisms, which in both cases have eluded detailed understanding so far. In this paper, we investigate numerically the bifurcation mechanisms at work in the incompressible Keplerian magnetorotational dynamo problem in the shearing box framework. Using numerical techniques imported from dynamical systems research, we show that the onset of chaotic dynamo action at magnetic Prandtl numbers larger than unity is primarily associated with global homoclinic and heteroclinic bifurcations of nonlinear magnetorotational dynamo cycles. These global bifurcations are supplemented by local bifurcations of cycles marking the beginning of period-doubling cascades. This suggests that nonlinear magnetorotational dynamo cycles provide the pathway to turbulent injection of both kinetic and magnetic energy in incompressible magnetohydrodynamic Keplerian shear flow in the absence of an externally imposed magnetic field. Studying the nonlinear physics and bifurcations of these cycles in different regimes and configurations may subsequently help to better understand the conditions of excitation of magnetohydrodynamic turbulence and instability-driven dynamos in various astrophysical systems and laboratory experiments. The detailed characterization of global bifurcations provided for this three-dimensional subcritical fluid dynamics problem may also prove useful for the problem of transition to turbulence in hydrodynamic shear flows.