Dust settling in local simulations of turbulent protoplanetary disks

In this paper, we study the effect of MHD turbulence on the dynamics of dust particles in protoplanetary disks. We vary the size of the particles and relate the dust evolution to the turbulent velocity fluctuations. We performed numerical simulations using two Eulerian MHD codes, both based on finite difference techniques: ZEUS--3D and NIRVANA. These were local shearing box simulations incorporating vertical stratification. Both ideal and non ideal MHD simulations with midplane dead zones were carried out. The codes were extended to incorporate different models for the dust as an additional fluid component. Good agreement between results obtained using the different approaches was obtained. The simulations show that a thin layer of very small dust particles is diffusively spread over the full vertical extent of the disk. We show that a simple description obtained using the diffusion equation with a diffusion coefficient simply expressed in terms of the velocity correlations accurately matches the results. Dust settling starts to become apparent for particle sizes of the order of 1 to 10 centimeters for which the gas begins to decouple in a standard solar nebula model at 5.2 AU. However, for particles which are 10 centimeters in size, complete settling toward a very thin midplane layer is prevented by turbulent motions within the disk, even in the presence of a midplane dead zone of significant size. These results indicate that, when present, MHD turbulence affects dust dynamics in protoplanetary disks. We find that the evolution and settling of the dust can be accurately modelled using an advection diffusion equation that incorporates vertical settling. The value of the diffusion coefficient can be calculated from the turbulent velocity field when that is known for a time of several local orbits.Comment: 15 pages, 16 figures, accepted in Astronomy & Astrophysic

MRI-driven angular momentum transport in protoplanetary disks

Angular momentum transport in accretion disk has been the focus of intense research in theoretical astrophysics for many decades. In the past twenty years, MHD turbulence driven by the magnetorotational instability has emerged as an efficient mechanism to achieve that goal. Yet, many questions and uncertainties remain, among which the saturation level of the turbulence. The consequences of the magnetorotational instability for planet formation models are still being investigated. This lecture, given in September 2012 at the school "Role and mechanisms of angular momentum transport in the formation and early evolution of stars" in Aussois (France), aims at introducing the historical developments, current status and outstanding questions related to the magnetorotational instability that are currently at the forefront of academic research.Comment: 51 pages, 16 figures, to appear in the proceedings of the Evry Schatzman School 2012 of PNPS and CNRS/INSU on the "Role and mechanisms of angular momentum transport during the formation and early evolution of stars", Eds. P.Hennebelle & C.Charbonne

Thermodynamics of the dead-zone inner edge in protoplanetary disks

In protoplanetary disks, the inner boundary between the turbulent and laminar regions could be a promising site for planet formation, thanks to the trapping of solids at the boundary itself or in vortices generated by the Rossby wave instability. At the interface, the disk thermodynamics and the turbulent dynamics are entwined because of the importance of turbulent dissipation and thermal ionization. Numerical models of the boundary, however, have neglected the thermodynamics, and thus miss a part of the physics. The aim of this paper is to numerically investigate the interplay between thermodynamics and dynamics in the inner regions of protoplanetary disks by properly accounting for turbulent heating and the dependence of the resistivity on the local temperature. Using the Godunov code RAMSES, we performed a series of 3D global numerical simulations of protoplanetary disks in the cylindrical limit, including turbulent heating and a simple prescription for radiative cooling. We find that waves excited by the turbulence significantly heat the dead zone, and we subsequently provide a simple theoretical framework for estimating the wave heating and consequent temperature profile. In addition, our simulations reveal that the dead-zone inner edge can propagate outward into the dead zone, before staling at a critical radius that can be estimated from a mean-field model. The engine driving the propagation is in fact density wave heating close to the interface. A pressure maximum appears at the interface in all simulations, and we note the emergence of the Rossby wave instability in simulations with extended azimuth. Our simulations illustrate the complex interplay between thermodynamics and turbulent dynamics in the inner regions of protoplanetary disks. They also reveal how important activity at the dead-zone interface can be for the dead-zone thermodynamic structure.Comment: 16 pages, 16 figures. Accepted in Astronomy and Astrophysic

Dissipative and nonaxisymmetric standard-MRI in Kepler disks

Deviations from axial symmetry are necessary to maintain self-sustained MRI-turbulence. We define the parameters region where nonaxisymmetric MRI is excited and study dependence of the unstable modes structure and growth rates on the relevant parameters. We solve numerically the linear eigenvalue problem for global axisymmetric and nonaxisymmetric modes of standard-MRI in Keplerian disks with finite diffusion. For small magnetic Prandtl number the microscopic viscosity completely drops out from the analysis so that the stability maps and the growth rates expressed in terms of the magnetic Reynolds number Rm and the Lundquist number S do not depend on the magnetic Prandtl number Pm. The minimum magnetic field for onset of nonaxisymmetric MRI grows with Rm. For given S all nonaxisymmetric modes disappear for sufficiently high Rm. This behavior is a consequence of the radial fine-structure of the nonaxisymmetric modes resulting from the winding effect of differential rotation. It is this fine-structure which presents severe resolution problems for the numerical simulation of MRI at large Rm. For weak supercritical magnetic fields only axisymmetric modes are unstable. Nonaxisymmetric modes need stronger fields and not too fast rotation. If Pm is small its real value does not play any role in MRI.Comment: 4 pages, 6 figures, A&A Lette

Spiral-driven accretion in protoplanetary discs - I. 2D models

We numerically investigate the dynamics of a 2D non-magnetised protoplanetary disc surrounded by an inflow coming from an external envelope. We find that the accretion shock between the disc and the inflow is unstable, leading to the generation of large-amplitude spiral density waves. These spiral waves propagate over long distances, down to radii at least ten times smaller than the accretion shock radius. We measure spiral-driven outward angular momentum transport with 1e-4 1e-8 Msun/yr. We conclude that the interaction of the disc with its envelope leads to long-lived spiral density waves and radial angular momentum transport with rates that cannot be neglected in young non-magnetised protostellar discs.Comment: 4 pages, 4 figures, accepted in A&A Letter

The influence of turbulence during magnetized core collapse and its consequences on low-mass star formation

[Abridged] Theoretical and numerical studies of star formation have shown that magnetic field (B) has a strong influence on both disk formation and fragmentation; even a relatively low B can prevent these processes. However, very few studies investigated the combined effects of B and turbulence. We study the effects of turbulence in magnetized core collapse, focusing on the magnetic diffusion, the orientation of the angular momentum (J) of the protostellar core, and on its consequences on disk formation, fragmentation and outflows. We perform 3D, AMR, MHD simulations of magnetically supercritical collapsing dense cores of 5 Msun using the MHD code RAMSES. A turbulent velocity field is imposed as initial conditions, characterised by a Kolmogorov power spectrum. Different levels of turbulence and magnetization are investigated, as well as 3 realisations for the turbulent velocity field. Magnetic diffusion, orientation of the rotation axis with respect to B, transport of J, disk formation, fragmentation and outflows formation are studied. The turbulent velocity field imposed as initial conditions contains a non-zero J, responsible for a misalignment of the rotation axis. Turbulence is also responsible for an effective turbulent diffusivity in the vicinity of the core. Both effects are responsible for a significant decrease of the magnetic braking, and facilitate the formation of early massive disks for not too high magnetization. Fragmentation can occur even with mu ~ 5 at late time in contrast with 1 Msun cores for which fragmentation is prevented for such values of mu. Slow asymmetric outflows are launched. They carry a mass which is comparable to the mass within the core. Because of misalignment and turbulent diffusion, massive disk formation is possible though their mass and size are still reduced compared to the hydrodynamical case. We find that for mu >= 5, fragmentation can happen.Comment: 15 pages, 21 figures, submitted in A&

MHD simulations of the magnetorotational instability in a shearing box with zero net flux. II. The effect of transport coefficients

We study the influence of the choice of transport coefficients (viscosity and resistivity) on MHD turbulence driven by the magnetorotational instability (MRI) in accretion disks. We follow the methodology described in paper I: we adopt an unstratified shearing box model and focus on the case where the net vertical magnetic flux threading the box vanishes. For the most part we use the finite difference code ZEUS, including explicit transport coefficients in the calculations. However, we also compare our results with those obtained using other algorithms (NIRVANA, the PENCIL code and a spectral code) to demonstrate both the convergence of our results and their independence of the numerical scheme. We find that small scale dissipation affects the saturated state of MHD turbulence. In agreement with recent similar numerical simulations done in the presence of a net vertical magnetic flux, we find that turbulent activity (measured by the rate of angular momentum transport) is an increasing function of the magnetic Prandtl number Pm for all values of the Reynolds number Re that we investigated. We also found that turbulence disappears when the Prandtl number falls below a critical value Pm_c that is apparently a decreasing function of Re. For the limited region of parameter space that can be probed with current computational resources, we always obtained Pm_c>1. We conclude that the magnitudes of the transport coefficients are important in determining the properties of MHD turbulence in numerical simulations in the shearing box with zero net flux, at least for Reynolds numbers and magnetic Prandtl numbers that are such that transport is not dominated by numerical effects and thus can be probed using current computational resources.Comment: 10 pages, 13 figures, accepted in A&A. Numerical results improved, minor changes in the tex