Astrophysical Disks Research Articles

Prospects for observing the magnetorotational instability in the plasma Couette experiment

Many astrophysical disks, such as protoplanetary disks, are in a regime where non-ideal, plasma-specific magnetohydrodynamic (MHD) effects can significantly influence the behaviour of the magnetorotational instability (MRI). The possibility of studying these effects in the plasma Couette experiment (PCX) is discussed. An incompressible, dissipative global stability analysis is developed to include plasma-specific two-fluid effects and neutral collisions, which are inherently absent in analyses of Taylor–Couette flows (TCFs) in liquid metal experiments. It is shown that with boundary driven flows, a ion-neutral collision drag body force significantly affects the azimuthal velocity profile, thus limiting the flows to regime where the MRI is not present. Electrically driven flow (EDF) is proposed as an alternative body force flow drive in which the MRI can destabilize at more easily achievable plasma parameters. Scenarios for reaching MRI relevant parameter space and necessary hardware upgrades are described.

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 nonlinear 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. (2013). 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. In all cases the fastest growing modes have very short wavelengths and are localised 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 destabilised 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 nonlinear saturation and resulting angular momentum transport.

ON THE ANISOTROPIC NATURE OF MRI-DRIVEN TURBULENCE IN ASTROPHYSICAL DISKS

The magnetorotational instability is thought to play an important role in enabling accretion in sufficiently ionized astrophysical disks. The rate at which MRI-driven turbulence transports angular momentum is related to both the strength of the amplitudes of the fluctuations on various scales and the degree of anisotropy of the underlying turbulence. This has motivated several studies of the distribution of turbulent power in spectral space. In this paper, we investigate the anisotropic nature of MRI-driven turbulence using a pseudo-spectral code and introduce novel ways to robustly characterize the underlying turbulence. We show that the general flow properties vary in a quasi-periodic way on timescales comparable to 10 inverse angular frequencies motivating the temporal analysis of its anisotropy. We introduce a 3D tensor invariant analysis to quantify and classify the evolution of the anisotropic turbulent flow. This analysis shows a continuous high level of anisotropy, with brief sporadic transitions towards two- and three-component isotropic turbulent flow. This temporal-dependent anisotropy renders standard shell-average, especially when used simultaneously with long temporal averages, inadequate for characterizing MRI-driven turbulence. We propose an alternative way to extract spectral information from the turbulent magnetized flow, whose anisotropic character depends strongly on time. This consists of stacking 1D Fourier spectra along three orthogonal directions that exhibit maximum anisotropy in Fourier space. The resulting averaged spectra show that the power along each of the three independent directions differs by several orders of magnitude over most scales, except the largest ones. Our results suggest that a first-principles theory to describe fully developed MRI-driven turbulence will likely have to consider the anisotropic nature of the flow at a fundamental level.

One-armed spirals in locally isothermal, radially structured self-gravitating discs

We describe a new mechanism that leads to the destabilisation of non-axisymmetric waves in astrophysical discs with an imposed radial temperature gradient. This might apply, for example, to the outer parts of protoplanetary discs. We use linear density wave theory to show that non-axisymmetric perturbations generally do not conserve their angular momentum in the presence of a forced temperature gradient. This implies an exchange of angular momentum between linear perturbations and the background disc. In particular, when the disturbance is a low-frequency trailing wave and the disc temperature decreases outwards, this interaction is unstable and leads to the growth of the wave. We demonstrate this phenomenon through numerical hydrodynamic simulations of locally isothermal discs in 2D using the FARGO code and in 3D with the ZEUS-MP and PLUTO codes. We consider radially structured discs with a self-gravitating region which remains stable in the absence of a temperature gradient. However, when a temperature gradient is imposed we observe exponential growth of a one-armed spiral mode (azimuthal wavenumber $m=1$) with co-rotation radius outside the bulk of the spiral arm, resulting in a nearly-stationary one-armed spiral pattern. The development of this one-armed spiral does not require the movement of the central star, as found in previous studies. Because destabilisation by a forced temperature gradient does not explicitly require disc self-gravity, we suggest this mechanism may also affect low-frequency one-armed oscillations in non-self-gravitating discs.

A viscous instability in axially symmetric laminar shear flows

A viscous instability in shearing laminar axisymmetric hydrodynamic flows around a gravitating center is described. In the linearized hydrodynamic equations written in the Boussinesq approximation with microscopic molecular transport coefficients, the instability arises when the viscous dissipation is taken into account in the energy equation. Using the local WKB approximation, we derive a third-order algebraic dispersion equation with two modes representing the modified Rayleigh modes R+ and R-, and the third X-mode. We show that in thin accretion flows the viscosity destabilizes one of the Rayleigh modes in a wide range of wavenumbers, while the X-mode always remains stable. In Keplerian flows, the instability increment is found to be a few Keplerian rotational periods at wavelengths with $kr\sim 10-50$. This instability may cause turbulence in astrophysical accretion discs even in the absence of magnetic field.

Dissipative effects on the sustainment of a magnetorotational dynamo in Keplerian shear flow

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.

Giant Planet Magnetodiscs and Aurorae—An Introduction

This volume contains the reports discussed during the Workshop “Giant Planet Magnetodiscs and Aurorae” held 26–30 November 2012, at the International Space Science Institute, organised together with the Europlanet project, supported by FP7 (Grant No. 228319). Magnetodiscs are large current sheets surrounding Jupiter and Saturn (also Uranus and Neptune) that are filled with plasma principally originating in the natural satellites of these worlds. They are also solar system analogues for astrophysical discs. Magnetodiscs are special features of the fast rotating giant planets, a special feature of rotationally driven magnetospheres. Their structure is modified by variability in their plasma sources and by the solar wind. Auroral signatures in the optical and radio wavebands allow a diagnostic of these dynamical processes and enable the visualisation of these large plasma and field structures.

Транзиентная динамика возмущений в астрофизических дисках

Транзиентная динамика возмущений в астрофизических дисках

Sensitivity of the magnetorotational instability to the shear parameter in stratified simulations

The magnetorotational instability (MRI) is a shear instability and thus its sensitivity to the shear parameter q = −d ln Ω/d ln r is of interest to investigate. Motivated by astrophysical discs, most (but not all) previous MRI studies have focused on the Keplerian value of q = 1.5. Using simulation with eight vertical density scaleheights, we contribute to the subset of studies addressing the effect of varying q in stratified numerical simulations. We discuss why shearing boxes cannot easily be used to study q > 2 and thus focus on q < 2. As per previous simulations, which were either unstratified or stratified with a smaller vertical domain, we find that the q dependence of stress for the stratified case is not linear, contrary to the Shakura–Sunyaev model. We find that the scaling agrees with Abramowicz, Brandenburg & Lasota who found it to be proportional to the shear to vorticity ratio q/(2 − q). We also find however that the shape of the magnetic and kinetic energy spectra are relatively insensitive to q and that the ratio of Maxwell stress to magnetic energy ratio also remains nearly independent of q. This is consistent with a theoretical argument in which the rate of amplification of the azimuthal field depends linearly on q and the turbulent correlation time τ depends inversely on q. As such, we measure the correlation time of the turbulence and find that indeed it is inversely proportional to q.

Scattering line polarization in rotating, optically thick disks

To interpret observations of astrophysical disks it is essential to understand the formation process of the emitted light. If the disk is optically thick, scattering dominated and permeated by a Keplerian velocity field, Non-Local Thermodynamic Equilibrium radiative transfer modeling must be done to compute the emergent spectrum from a given disk model. We investigate Non-local thermodynamic equilibrium polarized line formation in different simple disk models and aim to demonstrate the importance of both radiative transfer effects and scattering as well as the effects of velocity fields. We self-consistently solve the coupled equations of radiative transfer and statistical equilibrium for a two level atom model by means of Jacobi iteration. We compute scattering polarization, that is Q/I and U/I line profiles. The degree of scattering polarization is significantly influenced by the inclination of the disk with respect to observer, but also by the optical thickness of the disk and the presence of rotation. Stokes U shows double-lobed profiles with amplitude which increases with the disk rotation. Our results suggest that the line profiles, especially the polarized ones, emerging from gaseous disks differ significantly from the profiles predicted by simple approximations. The profiles are diverse in shape, but typically symmetric in Stokes Q and antisymmetric in Stokes U. A clear indicator of disk rotation is the presence of Stokes U, which might prove to be a useful diagnostic tool. We also demonstrate that, for moderate rotational velocities, an approximate treatment can be used, where non-local thermodynamic equilibrium radiative transfer is done in the velocity field-free approximation and Doppler shift is applied in the process of spatial integration over the whole emitting surface.

Hydrodynamic instability in eccentric astrophysical discs

Eccentric Keplerian discs are believed to be unstable to three-dimensional hydrodynamical instabilities driven by the time-dependence of fluid properties around an orbit. These instabilities could lead to small-scale turbulence, and ultimately modify the global disc properties. We use a local model of an eccentric disc, derived in a companion paper, to compute the nonlinear vertical ("breathing mode") oscillations of the disc. We then analyse their linear stability to locally axisymmetric disturbances for any disc eccentricity and eccentricity gradient using a numerical Floquet method. In the limit of small departures from a circular reference orbit, the instability of an isothermal disc is explained analytically. We also study analytically the small-scale instability of an eccentric neutrally stratified polytropic disc with any polytropic index using a WKB approximation. We find that eccentric discs are generically unstable to the parametric excitation of small-scale inertial waves. The nonlinear evolution of these instabilities should be studied in numerical simulations, where we expect them to lead to a decay of the disc eccentricity and eccentricity gradient as well as to induce additional transport and mixing. Our results highlight that it is essential to consider the three-dimensional structure of eccentric discs, and their resulting vertical oscillatory flows, in order to correctly capture their evolution.

Local and global dynamics of eccentric astrophysical discs

We formulate a local dynamical model of an eccentric disc in which the dominant motion consists of elliptical Keplerian orbits. The model is a generalization of the well known shearing sheet, and is suitable for both analytical and computational studies of the local dynamics of eccentric discs. It is spatially homogeneous in the horizontal dimensions but has a time-dependent geometry that oscillates at the orbital frequency. We show how certain averages of the stress tensor in the local model determine the large-scale evolution of the shape and mass distribution of the disc. The simplest solutions of the local model are laminar flows consisting of a (generally nonlinear) vertical oscillation of the disc. Eccentric discs lack vertical hydrostatic equilibrium because of the variation of the vertical gravitational acceleration around the eccentric orbit, and in some cases because of the divergence of the orbital velocity field associated with an eccentricity gradient. We discuss the properties of the laminar solutions, showing that they can exhibit extreme compressional behaviour for eccentricities greater than about $0.5$, especially in discs that behave isothermally. We also derive the linear evolutionary equations for an eccentric disc that follow from the laminar flows in the absence of a shear viscosity. In a companion paper we show that these solutions are linearly unstable and we determine the associated growth rates and unstable modes.

ACCRETION DISK RADIATION DYNAMICS AND THE COSMIC BATTERY

We investigate the dynamics of radiation in an astrophysical accretion disk around a Kerr black hole. The source of the radiation is the accretion disk itself, and not a central object as in previous studies of the Poynting-Robertson effect. We generate numerical sky maps from photon trajectories that originate on the surface of the disk as seen from the inner edge of the disk at the position of the innermost stable circular orbit (ISCO). We investigate several accretion disk morphologies with a Shakura-Sunyaev surface temperature distribution. Finally, we calculate the electromotive source of the Cosmic Battery mechanism around the inner edge of the accretion disk and obtain characteristic timescales for the generation of astrophysical magnetic fields.

Modification of the jeans instability criterion for fractal-structure astrophysical objects in the framework of nonextensive statistics

Unlike classical studies in which the gravitational instability criterion for astrophysical disks is derived in the framework of traditional kinetics or hydrodynamics, we propose to consider the totality of fluffy dust clusters of various astrophysical objects, in particular, protoplanetary subdisks, as a special type of continuous medium, i.e., fractal medium for which there are points and areas not filled with its components. Within the deformed Tsallis statistics formalism, which is intended to describe the behavior of anomalous systems with strong gravitational interaction and fractal nature of phase space, we derive, on the basis of the modified kinetic equation (with the collision integral in the Bhatnagar-Gross-Krook form), the generalized hydrodynamic Euler equations for a medium with the fractal mass dimension. Considering the linearization of the q-hydrodynamics equations, we investigate the instability of an infinitely homogeneous medium to obtain a simplified version of the modified gravitational instability criterion for an astrophysical disk with fractal structure.

Global simulations of magnetorotational turbulence - III. Influence of field configuration and mass injection

The stresses produced by magnetorotational turbulence can provide effective angular momentum transport in accretion disks. However, questions remain about the ability of simulated disks to reproduce observationally inferred stress-to-gas-pressure ratios. In this paper we present a set of high resolution global magnetohydrodynamic disk simulations which are initialised with different field configurations: purely toroidal, vertical field lines, and nested poloidal loops. A mass source term is included which allows the total disk mass to equilibrate in simulations with long run times, and also enables the impact of rapid mass injection to be explored. Notably different levels of angular momentum transport are observed during the early-time transient disk evolution. However, given sufficient time to relax, the different models evolve to a statistically similar quasi-steady state with a stress-to-gas-pressure ratio, $\alpha \sim 0.032-0.036$. The indication from our results is that {\it steady, isolated} disks may be unable to maintain a large-scale magnetic field or produce values for the stress-to-gas-pressure ratio implied by some observations. Supplementary simulations exploring the influence of trapping magnetic field, injecting vertical field, and rapidly injecting additional mass into the disk show that large stresses ($\alpha \sim 0.1-0.25$) can be induced by these mechanisms. The simulations highlight the common late-time evolution and characteristics of turbulent disks for which the magnetic field is allowed to evolve freely. If the boundaries of the disk, the rate of injection of magnetic field, or the rate of mass replenishment are modified to mimic astrophysical disks, markedly different disk evolution occurs.

Turbulence decay towards the linearly stable regime of Taylor–Couette flow

AbstractTaylor–Couette (TC) flow is used to probe the hydrodynamical (HD) stability of astrophysical accretion disks. Experimental data on the subcritical stability of TC flow are in conflict about the existence of turbulence (cf. Ji et al. (Nature, vol. 444, 2006, pp. 343–346) and Paoletti et al. (Astron. Astroph., vol. 547, 2012, A64)), with discrepancies attributed to end-plate effects. In this paper we numerically simulate TC flow with axially periodic boundary conditions to explore the existence of subcritical transitions to turbulence when no end plates are present. We start the simulations with a fully turbulent state in the unstable regime and enter the linearly stable regime by suddenly starting a (stabilizing) outer cylinder rotation. The shear Reynolds number of the turbulent initial state is up to $Re_s \lesssim 10^5$ and the radius ratio is $\eta =0.714$. The stabilization causes the system to behave as a damped oscillator and, correspondingly, the turbulence decays. The evolution of the torque and turbulent kinetic energy is analysed and the periodicity and damping of the oscillations are quantified and explained as a function of shear Reynolds number. Though the initially turbulent flow state decays, surprisingly, the system is found to absorb energy during this decay.

Asymmetric evolution of magnetic reconnection in collisionless accretion disk

An evolution of a magnetic reconnection in a collisionless accretion disk is investigated using a 2.5 dimensional hybrid code simulation. In astrophysical disks, magnetorotational instability (MRI) is considered to play an important role by generating turbulence in the disk and contributes to an effective angular momentum transport through a turbulent viscosity. Magnetic reconnection, on the other hand, also plays an important role on the evolution of the disk through a dissipation of a magnetic field enhanced by a dynamo effect of MRI. In this study, we developed a hybrid code to calculate an evolution of a differentially rotating system. With this code, we first confirmed a linear growth of MRI. We also investigated a behavior of a particular structure of a current sheet, which would exist in the turbulence in the disk. From the calculation of the magnetic reconnection, we found an asymmetric structure in the out-of-plane magnetic field during the evolution of reconnection, which can be understood by a coupling of the Hall effect and the differential rotation. We also found a migration of X-point whose direction is determined only by an initial sign of J0×Ω0, where J0 is the initial current density in the neutral sheet and Ω0 is the rotational vector of the background Keplerian rotation. Associated with the migration of X-point, we also found a significant enhancement of the perpendicular magnetic field compared to an ordinary MRI. MRI-Magnetic reconnection coupling and the resulting magnetic field enhancement can be an effective process to sustain a strong turbulence in the accretion disk and to a transport of angular momentum.

EFFECTIVE Q CRITERION FOR DISK STABILITY IN AN EXTERNAL POTENTIAL

standard Q criterion (with Q > 1) describes the stability against local, axisymmetric perturbations in a disk supported by rotation and random motion. Most astrophysical disks, however, are under the influence of an external gravitational potential, which can significantly affect their stability. A typical example is a galactic disk embedded in a dark matter halo. Here, we do a linear perturbation analysis for a disk in an external potential and obtain a generalized dispersion relation and the effective stability criterion. An external potential, such as that due to the dark matter halo concentric with the disk, contributes to the unperturbed rotational field and significantly increases its stability. We obtain the values for the effective Q parameter for the Milky Way and for a low surface brightness galaxy, UGC 7321. We find that in each case the stellar disk by itself is barely stable and it is the dark matter halo that stabilizes the disk against local, axisymmetric gravitational instabilities. Thus, the dark matter halo is necessary to ensure local disk stability. This result has been largely missed so far because in practice the Q parameter for a galactic disk is obtained using the observed rotational field that already includes the effect of the halo.

Rossby wave instability in astrophysical discs

A brief review is given of the Rossby wave instability in astrophysical discs. In non-self-gravitating discs, around for example a newly forming stars, the instability can be triggered by an axisymmetric bump at some radius r0 in the disc surface mass-density. It gives rise to exponentially growing non-axisymmetric perturbation (, m = 1,2,…) in the vicinity of r0 consisting of anticyclonic vortices. These vortices are regions of high pressure and consequently act to trap dust particles which in turn can facilitate planetesimal growth in proto-planetary discs. The Rossby vortices in the discs around stars and black holes may cause the observed quasi-periodic modulations of the disc's thermal emission.

Accretion in giant planet circumplanetary discs

During the final growth phase of giant planets, accretion is thought to be controlled by a surrounding circumplanetary disk. Current astrophysical accretion disk models rely on hydromagnetic turbulence or gravitoturbulence as the source of effective viscosity within the disk. However, the magnetically-coupled accreting region in these models is so limited that the disk may not support inflow at all radii, or at the required rate. Here, we examine the conditions needed for self-consistent accretion, in which the disk is susceptible to accretion driven by magnetic fields or gravitational instability. We model the disk as a Shakura-Sunyaev $\alpha$ disk and calculate the level of ionisation, the strength of coupling between the field and disk using Ohmic, Hall and Ambipolar diffusevities for both an MRI and vertical field, and the strength of gravitational instability. We find that the standard constant-$\alpha$ disk is only coupled to the field by thermal ionisation within $30\,R_J$ with strong magnetic diffusivity prohibiting accretion through the bulk of the midplane. In light of the failure of the constant-$\alpha$ disk to produce accretion consistent with its viscosity we drop the assumption of constant-$\alpha$ and present an alternate model in which $\alpha$ varies radially according to the level magnetic turbulence or gravitoturbulence. We find that a vertical field may drive accretion across the entire disk, whereas MRI can drive accretion out to $\sim200\,R_J$, beyond which Toomre's $Q=1$ and gravitoturbulence dominates. The disks are relatively hot ($T\gtrsim800\,$K), and consequently massive ($M_{\text{disk}}\sim0.5\,M_J$).