Henrik N. Latter

Local outflows from turbulent accretion disks

Aims. The aim of this paper is to investigate the properties of accretion disks threaded by a weak vertical magnetic field, with a particular focus on the interplay between magnetohydrodynamic (MHD) turbulence driven by the magnetorotational instability (MRI) and outflows that might be launched from the disk. Methods. For that purpose, we use a set of numerical simulations performed with the MHD code RAMSES in the framework of the shearing box model. We concentrate on the case of a rather weak vertical magnetic field such that the initial ratio β0 of the thermal and magnetic pressures in the disk midplane equals 10 4 . Results. As reported recently, we find that MHD turbulence drives an efficient outflow out of the computational box. We demonstrate a strong sensitivity of that result to the box size: enlargements in the radial and vertical directions lead to a reduction of up to an order of magnitude in the mass-loss rate. Such a dependence prevents any realistic estimates of disk mass-loss rates being derived using shearing-box simulations. We find however that the flow morphology is robust and independent of the numerical details of the simulations. Its properties display some features and approximate invariants that are reminiscent of the Blandford & Payne launching mechanism, but differences exist. For the magnetic field strength considered in this paper, we also find that angular momentum transport is most likely dominated by MHD turbulence, the saturation of which scales with the magnetic Prandtl number, the ratio of viscosity and resistivity, in a way that is in good agreement with expectations based on unstratified simulations. Conclusions. This paper thus demonstrates for the first time that accretion disks can simultaneously exhibit MRI-driven MHD turbulence along with magneto-centrifugally accelerated outflows. However, in contradiction with previously published results, such outflows probably have little impact on the disk dynamics.

MRI channel flows and their parasites

Local simulations of the magnetorotational instability (MRI) in accretion discs can exhibit recurrent coherent structures called channel flows. The formation and destruction of these structures may play a role in the development and saturation of MRI-induced turbulence, and consequently help us understand the time-dependent accretion behaviour of certain astrophysical objects. Previous investigations have revealed that channel solutions are attacked by various parasitic modes, foremost of which is an analogue of the Kelvin-Helmholtz instability. We revisit these instabilities and show how they relate to the classical instabilities of plasma physics, the kink and pinch modes. However, we argue that in most cases channels emerge from developed turbulence and are eventually destroyed by turbulent mixing, not by the parasites. The exceptions are clean isolated channels, which appear in systems near criticality or which emerge from low amplitude initial conditions. These structures can achieve large amplitudes and are only then destroyed, giving rise to eruptive behaviour.

Differential rotation and convection in the Sun

We show that the differential rotation profile of the solar convection zone, apart from inner and outer boundary layers, can be reproduced with great accuracy if the isorotation contours correspond to characteristics of the thermal wind equation. This requires that there be a formal quantitative relationship involving the entropy and the angular velocity. Earlier work has suggested that this could arise from magnetohydrodynamic stability constraints; here, we argue that purely hydrodynamical processes could also lead to such a result. Of special importance to the hydrodynamical solution is the fact that the thermal wind equation is insensitive to radial entropy gradients. This allows a much more general class of solutions to fit the solar isorotation contours, beyond just those in which the entropy itself must be a function of the angular velocity. In particular, for this expanded class, the thermal wind solution of the solar rotation profile remains valid even when large radial entropy gradients are present. A clear and explicit example of this class of solution appears to be present in published numerical simulations of the solar convective zone. Though hydrodynamical in character, the theory is not sensitive to the presence of weak magnetic fields. Thus, the identification of solar isorotation contours with the characteristics of the thermal wind equation appears to be robust, accommodating, but by no means requiring, magnetic field dynamics.

MRI channel flows in vertically stratified models of accretion discs

Simulations of the magnetorotational instability (MRI) in ‘unstratified’ shearing boxes exhibit powerful coherent flows, whereby the fluid vertically splits into countermoving planar jets or ‘channels’. Channel flows correspond to certain axisymmetric linear MRI modes, and their preponderance follows from the remarkable fact that they are approximate non-linear solutions of the MHD equations in the limit of weak magnetic fields. We show in this paper, analytically and with one-dimensional numerical simulations, that this property is also shared by certain axisymmetric MRI modes in vertically stratified shearing boxes. These channel flows rapidly capture significant amounts of magnetic and kinetic energy, and thus are vulnerable to secondary shear instabilities. We examine these parasites in the vertically stratified context, and estimate the maximum amplitudes that channels attain before they are destroyed. These estimates suggest that a dominant channel flow will usually drive the discs magnetic field to thermal strengths. The prominence of these flows and their destruction place enormous demands on simulations, but channels in their initial stages also offer a useful check on numerical codes. These benchmarks are especially valuable given the increasing interest in the saturation of the stratified MRI. Lastly, we speculate on the potential connection between ‘run-away’ channel flows and outburst behaviour in protostellar and dwarf nova discs.

On the vertical-shear instability in astrophysical discs

We explore the linear stability of astrophysical discs exhibiting vertical shear, which arises when there is a radial variation in the temperature or entropy. Such discs are subject to a ‘vertical-shear instability’, which recent non-linear simulations have shown to drive hydrodynamic activity in the MRI-stable regions of protoplanetary discs. We first revisit locally isothermal discs using the quasi-global reduced model derived by Nelson et al. This analysis is then extended to global axisymmetric perturbations in a cylindrical domain. We also derive and study a reduced model describing discs with power-law radial entropy profiles (‘locally polytropic discs’), which are somewhat more realistic in that they possess physical (as opposed to numerical) surfaces. The fastest growing modes have very short wavelengths and are localized at the disc surfaces (if present), where the vertical shear is maximal. An additional class of modestly growing vertically global body modes is excited, corresponding to destabilized classical inertial waves (‘r modes’). We discuss the properties of both types of modes, and stress that those that grow fastest occur on the shortest available length-scales (determined either by the numerical grid or the physical viscous length). This ill-posedness makes simulations of the instability difficult to interpret. We end with some brief speculation on the non-linear saturation and resulting angular momentum transport.

Hydrodynamic instability in warped astrophysical discs

Warped astrophysical discs are usually treated as laminar viscous flows, which have anomalous properties when the disc is nearly Keplerian and the viscosity is small: fast horizontal shearing motions and large torques are generated, which cause the warp to evolve rapidly, in some cases at a rate that is inversely proportional to the viscosity. However, these flows are often subject to a linear hydrodynamic instability, which may produce small-scale turbulence and modify the large-scale dynamics of the disc. We use a warped shearing sheet to compute the oscillatory laminar flows in a warped disc and to analyse their linear stability by the Floquet method. We find widespread hydrodynamic instability deriving from the parametric resonance of inertial waves. Even very small, unobservable warps in nearly Keplerian discs of low viscosity can be expected to generate hydrodynamic turbulence, or at least wave activity, by this mechanism.

Viscous overstability and eccentricity evolution in three-dimensional gaseous discs

We investigate the growth or decay rate of the fundamental mode of even symmetry in a viscous accretion disc. This mode occurs in eccentric discs and is known to be potentially overstable. We determine the vertical structure of the disc and its modes, treating radiative energy transport in the diffusion approximation. In the limit of very long radial wavelength, an analytical criterion for viscous overstability is obtained, which involves the effective shear and bulk viscosity, the adiabatic exponent, and the opacity law of the disc. This differs from the prediction of a two-dimensional model. On shorter wavelengths (a few times the disc thickness), the criterion for overstability is more difficult to satisfy because of the different vertical structure of the mode. In a low-viscosity disc a third regime of intermediate wavelengths appears, in which the overstability is suppressed as the horizontal velocity perturbations develop significant vertical shear. We suggest that this effect determines the damping rate of eccentricity in protoplanetary discs, for which the long-wavelength analysis is inapplicable and overstability is unlikely to occur on any scale. In thinner accretion discs and in decretion discs around Be stars overstability may occur only on the longest wavelengths, leading to the preferential excitation of global eccentric modes.

Hysteresis and thermal limit cycles in MRI simulations of accretion discs

The recurrent outbursts that characterize low-mass binary systems reflect thermal state changes in their associated accretion discs. The observed outbursts are connected to the strong variation in disc opacity as hydrogen ionizes near 5000 K. This physics leads to accretion disc models that exhibit bistability and thermal limit cycles, whereby the disc jumps between a family of cool and low-accreting states and a family of hot and efficiently accreting states. Previous models have parametrized the disc turbulence via an alpha (or ‘eddy’) viscosity. In this paper we treat the turbulence more realistically via a suite of numerical simulations of the magnetorotational instability (MRI) in local geometry. Radiative cooling is included via a simple but physically motivated prescription. We show the existence of bistable equilibria and thus the prospect of thermal limit cycles, and in so doing demonstrate that MRI-induced turbulence is compatible with the classical theory. Our simulations also show that the turbulent stress and pressure perturbations are only weakly dependent on each other on orbital times; as a consequence, thermal instability connected to variations in turbulent heating (as opposed to radiative cooling) is unlikely to operate, in agreement with previous numerical results. Our work presents a first step towards unifying simulations of full magnetohydrodynamic turbulence with the correct thermal and radiative physics of the outbursting discs associated with dwarf novae, low-mass X-ray binaries and possibly young stellar objects.

Tidal disruption of satellites and formation of narrow rings

In this paper, we investigate the formation of narrow planetary rings such as those found around Uranus and Saturn through the tidal disruption of a weak, gravitationally bound satellite that migrates within its Roche limit. Using N-body simulations, we study the behaviour of rubble piles placed on circular orbits at different distances from a central planet. We consider both homogeneous satellites and differentiated bodies containing a denser core. We show that the Roche limit for a rubble pile is closer to the planet than for a fluid body of the same mean density. The Roche limit for a differentiated body is also closer to the planet than for a homogeneous satellite of the same mean density. Within its Roche limit, a homogeneous satellite totally disrupts and forms a narrow ring. The initial stages of the disruption are similar to the evolution of a viscous fluid ellipsoid, which can be computed semi-analytically. On the other hand, when a differentiated satellite is just within the Roche limit only the mantle is disrupted. This process is similar to Roche lobe overflow in interacting binary stars and produces two narrow rings on either side of a remnant satellite. We argue that the Uranian rings, and possibly their shepherd satellites, could have been formed through the tidal disruption of a number of protomoons that were formed inside the corotation radius of Uranus and migrated slowly inwards as a result of tidal interaction with the planet.

A comparison of local simulations and reduced models of MRI-induced turbulence

We run mean-field shearing-box numerical simulations with a temperature-dependent resistivity and compare them to a reduced dynamical model. Our simulations reveal the co-existence of two quasi-steady states, a ‘quiet’ state and an ‘active’ turbulent state, confirming the predictions of the reduced model. The initial conditions determine on which state the simulation ultimately settles. The active state is strongly influenced by the geometry of the computational box and the thermal properties of the gas. Cubic domains support permanent channel flows, bar-shaped domains exhibit eruptive behaviour, and horizontal slabs give rise to infrequent channels. Meanwhile, longer cooling time-scales lead to higher saturation amplitudes.