The fundamental difference between incompressible ideal magnetohydrodynamics and the dynamics of a non-conducting fluid is that magnetic fields exert a tension force that opposes their bending; magnetic fields behave like elastic strings threading the fluid. It is natural, therefore, to expect that a magnetic field tangled at small length scales should resist a large-scale shear in an elastic way, much as a ball of tangled elastic strings responds elastically to an impulse. Furthermore, a tangled field should support the propagation of `magnetoelastic waves', the isotropic analogue of Alfv\'en waves on a straight magnetic field. Here, we study magnetoelasticity in the idealised context of an equilibrium tangled field configuration. In contrast to previous treatments, we explicitly account for intermittency of the Maxwell stress, and show that this intermittency necessarily decreases the frequency of magnetoelastic waves in a stable field configuration. We develop a mean-field formalism to describe magnetoelastic behaviour, retaining leading-order corrections due to the coupling of large- and small-scale motions, and solve the initial-value problem for viscous fluids subjected to a large-scale shear, showing that the development of small-scale motions results in anomalous viscous damping of large-scale waves. Finally, we test these analytic predictions using numerical simulations of standing waves on tangled, linear force-free magnetic-field equilibria.

We model the light curves of the novel and extremely luminous transient ASASSN-15lh at nine different frequencies, from infrared to ultra-violet photon energies, as an evolving relativistic disc produced in the aftermath of a tidal disruption event. Good fits to all nine light curves are simultaneously obtained when $M_{\rm acc} \simeq 0.07 M_{\odot}$ is accreted onto a black hole of mass $M \simeq 10^9 M_{\odot}$, and near maximal rotation $a/r_g = 0.99$. The best fit black hole mass is consistent with a number of existing estimates from galactic scaling relationships. If confirmed, our results represent the detection of one of the most massive rapidly spinning black holes to date, and are strong evidence for a TDE origin for ASASSN-15lh. This would be the first TDE to be observed in the disc dominated state at optical and infrared frequencies.

About 10% of stars more massive than ${\approx}\,1.5\,\mathrm{M}_\odot$ have strong, large-scale surface magnetic fields and are being discussed as progenitors of highly-magnetic white dwarfs and magnetars. The origin of these fields remains uncertain. Recent 3D magnetohydrodynamical simulations have shown that strong magnetic fields can be generated in the merger of two massive stars. Here, we follow the long-term evolution of such a 3D merger product in a 1D stellar evolution code. During a thermal relaxation phase after the coalescence, the merger product reaches critical surface rotation, sheds mass and then spins down primarily because of internal mass readjustments. The spin of the merger product after thermal relaxation is mainly set by the co-evolution of the star-torus structure left after coalescence. This evolution is still uncertain, so we also consider magnetic braking and other angular-momentum-gain and -loss mechanisms that may influence the final spin of the merged star. Because of core compression and mixing of carbon and nitrogen in the merger, enhanced nuclear burning drives a transient convective core that greatly contributes to the rejuvenation of the star. Once the merger product relaxed back to the main sequence, it continues its evolution similar to that of a genuine single star of comparable mass. It is a slow rotator that matches the magnetic blue straggler $\tau$ Sco. Our results show that merging is a promising mechanism to explain some magnetic massive stars and it may also be key to understand the origin of the strong magnetic fields of highly-magnetic white dwarfs and magnetars.

We perform a detailed numerical and analytical study of the properties of observed light curves from relativistic thin discs, focussing on observational bands most appropriate for comparison with tidal disruption events (TDEs). We make use of asymptotic expansion techniques applied to the spectral emission integral, using time dependent disc temperature profiles appropriate for solutions of the relativistic thin disc equation. Rather than a power law associated with bolometric disc emission $L \sim t^{-n}$, the observed X-ray flux from disc-dominated TDEs will typically have the form of a power law multiplied by an exponential (see eq. 91). While precise details are somewhat dependent on the nature of the ISCO stress and disc-observer orientational angle, the general form of the time-dependent flux is robust and insensitive to the exact disc temperature profile. We present numerical fits to the UV and X-ray light curves of ASASSN-14li, a particularly well observed TDE. This modelling incorporates strong gravity optics. The full 900 days of ASASSN-14li X-ray observations are very well fit by a simple relativistic disc model, significantly improving upon previous work. The same underlying model also fits the final 1000 days of ASASSN-14li observations in three different UV bandpasses. Finally, we demonstrate that the analytic formulae reproduce the properties of full numerical modelling at both UV and X-ray wavelengths with great fidelity.

Magnetic fields are ubiquitous in the Universe. The Sun's magnetic field drives the solar wind and causes solar flares and other energetic surface phenomena that profoundly affect space weather here on Earth. The first magnetic field in a star other than the Sun was detected in 1947 in the peculiar A-type star 78 Vir. It is now known that the magnetic fields of the Sun and other low-mass stars (<1.5 solar masses) are generated in-situ by a dynamo process in their turbulent, convective envelopes. Unlike such stars, intermediate-mass and high-mass stars (>1.5 solar masses; referred to as "massive" stars here) have relatively quiet, radiative envelopes where a solar-like dynamo cannot operate. However, about 10% of them, including 78 Vir, have strong, large-scale surface magnetic fields whose origin has remained a major mystery. The massive star $\tau$ Sco is a prominent member of this group and appears to be surprisingly young compared to other presumably coeval members of the Upper Scorpius association. Here, we present the first 3D magneto-hydrodynamical simulations of the coalescence of two massive main-sequence stars and 1D stellar evolution computations of the subsequent evolution of the merger product that can explain $\tau$ Sco's magnetic field, apparent youth and other observed characteristics. We argue that field amplification in stellar mergers is a general mechanism to form strongly-magnetised massive stars. These stars are promising progenitors of those neutron stars that host the strongest magnetic fields in the Universe, so-called magnetars, and that may give rise to some of the enigmatic fast radio bursts. Strong magnetic fields affect the explosions of core-collapse supernovae and, moreover, those magnetic stars that have rapidly-rotating cores at the end of their lives might provide the right conditions to power long-duration gamma-ray bursts and super-luminous supernovae.

We present solutions to the relativistic thin disc evolutionary equation using an ${\alpha}$-model for the turbulent stress tensor. Solutions with a finite stress at the innermost stable circular orbit (ISCO) give rise to bolometric light curves with a shallow power law time dependence, in good agreement with those observed in tidal disruption events. A self-similar model based on electron scattering opacity, for example, yields a power law index of $-11/14$, as opposed to $-19/16$ for the case of zero ISCO stress. These solutions correspond to an extended period of relaxation of the evolving disc which, like the light curves they produce, is not sustainable indefinitely. Cumulative departures from the approximation of exact circular orbits cause the power law index to evolve slowly with time, leading eventually to the steeper fall-off associated with traditional zero ISCO stress models. These modified solutions are discussed in detail in a companion paper.

We present solutions to the relativistic thin disc evolutionary equation using a modified description of the mean fluid flow within the disc. The model takes into account the effects of sub-circular velocities in the innermost disc regions, and resolves otherwise unsustainable behaviour present in simple finite ISCO stress disc models. We show that the behaviour of a relativistic thin disc evolving with a finite ISCO stress is comprised of three distinct stages which join the ordinarily distinct finite and vanishing ISCO stress solutions into a fully continuous model parameterisation. The most important prediction of our model is the existence of an intermediate stage of "stalled accretion", controlled by a single dimensionless parameter. The hallmarks of this evolutionary phase appear to have been seen in GRMHD simulations as well as in the late time X-ray observations of tidal disruption events, but dedicated simulations and extended observations are needed for a deeper understanding.

This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed.

An encounter between a passing star and a massive black hole at the centre of a galaxy, a so-called tidal disruption event or TDE, may leave a debris disc that subsequently accretes onto the hole. We solve for the time evolution of such a TDE disc, making use of an evolutionary equation valid for both the Newtonian and Kerr regimes. The late time luminosity emergent from such a disc is of interest as a model diagnostic, as it tends to follow a power law decline. The original simple ballistic fallback model, with equal mass in equal energy intervals, produces a -5/3 power law, while standard viscous disc descriptions yield a somewhat more shallow decline, with an index closer to -1.2. Of four recent, well-observed tidal disruption event candidates however, all had fall-off power law indices smaller than 1 in magnitude. In this work, we revisit the problem of thin disc evolution, solving this reduced problem in full general relativity. Our solutions produce power law indices that are in much better accord with observations. The late time observational data from many TDEs are generally supportive, not only of disc accretion models, but of finite stress persisting down to the innermost stable circular orbit.

In the classical theory of thin disc accretion discs, the constraints of mass and angular momentum conservation lead to a diffusion-like equation for the turbulent evolution of the surface density. Here, we revisit this problem, extending the Newtonian analysis to the regime of Kerr geometry relevant to black holes. A diffusion-like equation once again emerges, but now with a singularity at the radius at which the effective angular momentum gradient passes through zero. The equation may be analysed using a combination of WKB, local techniques, and matched asymptotic expansions. It is shown that imposing the boundary condition of a vanishing stress tensor (more precisely the radial-azimuthal component thereof) allows smooth stable modes to exist external to the angular momentum singularity, the innermost stable circular orbit, while smoothly vanishing inside this location. The extension of the disc diffusion equation to the domain of general relativity introduces a new tool for numerical and phenomenolgical studies of accretion discs, and may prove to be a useful technique for understanding black hole X-ray transients.

In the standard thin disc formalism the dimensionless $\alpha$ parameter is usually assumed to be constant. However, there are good theoretical reasons for believing, as well as evidence from simulations, that $\alpha$ is dependent on intrinsic disc properties. In this paper we analyse the conditions for the stability of a thin accretion disc in which $\alpha$ is a function of the magnetic Prandtl number, the ratio of collisional viscosity to resistivity. In the inner disc, where the free electron opacity and radiation viscosity dominate, the disc is unstable if $\alpha$ is proportional to the magnetic Prandtl number with an exponent, $n$, and $6/13<n<10/3$. This is within the range of values for the power-law index found in MHD simulations with simple energetics. We calculate the evolution of the unstable disc within the $\alpha$ formalism and show that the physically accessible solutions form a limit cycle, analogous to the behaviour seen in recurrent dwarf novae. It is noteworthy that the time-dependent global behaviour of the instability results in cyclic heating of the inner section of the disc, when parameters appropriate for an X-ray binary system are used. We calculate a model spectrum of the disc in the flaring and quiescent states and show that the behaviour is compatible with X-ray observations of the thermal accretion disc in flaring X-ray binary systems.

Understanding what determines the strength of MHD turbulence in accretion discs is a question of fundamental theoretical and observational importance. In this work we investigate whether the dependence of the turbulent accretion disc stress ($\alpha$) on the magnetic Prandtl number (Pm) is sufficiently sensitive to induce thermal-viscous instability using 3D MHD simulations. We first investigate whether the $\alpha$-Pm dependence, found by many previous authors, has a physical or numerical origin by conducting a suite of local shearing-box simulations. We find that a definite $\alpha$-Pm dependence persists when simultaneously increasing numerical resolution and decreasing the absolute values of both the viscous and resistive dissipation coefficients. This points to a physical origin of the $\alpha$-Pm dependence. Using a further set of simulations which include realistic turbulent heating and radiative cooling, and by giving Pm a realistic physical dependence on the plasma temperature and density, we demonstrate that the $\alpha$-Pm dependence is sufficiently strong to lead to a local instability. We confirm that the instability manifests itself as an unstable limit cycle by mapping the local thermal-equilibrium curve of the disc. This is the first self-consistent MHD simulation demonstrating the Pm instability from first principles. This result is important because a physical Pm instability would lead to the global propagation of heating and cooling fronts and a transition between disc states on timescales compatible with the observed hard/soft state transitions in black hole binaries.

A conserved stress energy tensor for weak field gravitational waves propagating in vacuum is derived directly from the linearized wave equation alone, for an arbitrary gauge using standard general relativity. In any harmonic gauge, the form of the tensor leads directly to the classical expression for the outgoing wave energy. The method described here, however, is a much simpler, shorter, and more physically motivated approach than is the customary procedure, which involves a lengthy and cumbersome second-order (in wave-amplitude) calculation starting with the Einstein tensor. Our method has the added advantage of exhibiting the direct coupling between the outgoing energy flux in gravitational waves and the work done by the gravitational field on the sources. For nonharmonic gauges, the directly derived wave stress tensor has an apparent index asymmetry. This coordinate artefact may be removed, and the symmetrized (still gauge-invariant) tensor then takes on its widely-used form. Angular momentum conservation follows immediately. For any harmonic gauge, however, the stress tensor is manifestly symmetric from the start and its derivation depends, in its entirety, on the structure of the linearized wave equation.

The Goldreich-Schubert-Fricke (GSF) instability is a rotational instability that is thought to contribute to the transfer of angular momentum in differentially rotating stars. It has been included in recent codes of stellar evolution in a diffusion-like approximation, under the assumption that the kinematic viscosity $\nu$ is unimportant for the development of the instability. As noted previously by other authors, for most stellar applications this may not be a valid approximation. We discuss this issue in detail, solving the dispersion relation of the perturbed modes for realistic values of $\nu$ in the bulk of the radiative zone of the Sun and of three red giant stars at various evolutionary stages. We find that the GSF instability is triggered only in layers of extremely strong shear. In a simple case study, we also investigate the effect of a small deviation from axisymmetry or a small background magnetic field. We find that, like the viscosity, these have a stabilising effect. We conclude that this instability is probably far less efficient in transporting angular momentum than is often assumed, and may not even be present.

Much of astrophysics consists of the study of ionised gas under the influence of gravitational and magnetic fields. Thus, it is not possible to understand the astrophysical universe without a detailed knowledge of the dynamics of magnetised fluids. Fluid dynamics is, however, a notoriously tricky subject, in which it is all too easy for one's a priori intuition to go astray. In this review, we seek to guide the reader through a series of illuminating yet deceptive problems, all with an enlightening twist. We cover a broad range of topics including the instabilities acting in accretion discs, the hydrodynamics governing the convective zone of the Sun, the magnetic shielding of a cooling galaxy cluster, and the behaviour of thermal instabilities and evaporating clouds. The aim of this review is to surprise and intrigue even veteran astrophysical theorists with an idiosynchratic choice of problems and counterintuitive results. At the same time, we endeavour to bring forth the fundamental ideas, to set out important assumptions, and to describe carefully whatever novel techniques may be appropriate to the problem at hand. By beginning at the beginning, and analysing a wide variety of astrophysical settings, we seek not only to make this review suitable for fluid dynamic veterans, but to engage novice recruits as well with what we hope will be an unusual and instructive introduction to the subject.

Barotropic rotation and radiative equilibrium are mutually incompatible in stars. The issue is often addressed by allowing for a meridional circulation, but this is not devoid of theoretical complications. Models of rotation in the Sun which maintain strict radiative equilibrium, making use of the observation that the Sun is not in a state of barotropic rotation, have recently been suggested. To investigate the dynamical behaviour of these solutions, we study the local stability of stratified, weakly magnetized, differentially rotating fluids to non-axisymmetric perturbations. Finite heat conductivity, kinematic viscosity, and resistivity are present. The evolution of local embedded perturbations is governed by a set of coupled, ordinary differential equations with time-dependent coefficients. Two baroclinic models of rotation for the upper radiative zone and tachocline are studied: (i) an interpolation based on helioseismology data, (ii) a theoretical solution directly compatible with radiative equilibrium. The growth of the local Goldreich-Schubert-Fricke instability appears to be suppressed, largely because of the viscosity. An extensive exploration of wavenumber space is carried out, with and without a magnetic field. Although we easily find classical local instabilities when they ought formally to be present, for the Sun the analysis reveals neither unstable solutions, nor even solutions featuring a large transient growth. We have not ruled out larger scale or nonlinear instabilities, nor have we rigorously proven local stability. But rotational configurations in close agreement with observations, generally thought to be vulnerable to the classic local Goldreich-Schubert-Fricke instability, do appear to be locally stable under rather general circumstances.

Central jetted active galactic nuclei (AGN) appear to heat the core regions of the intracluster medium (ICM) in cooling-core galaxy clusters and groups, thereby preventing a cooling catastrophe. However, the physical mechanism(s) by which the directed flow of kinetic energy is thermalized throughout the ICM core remains unclear. We examine one widely discussed mechanism whereby the AGN induces subsonic turbulence in the ambient medium, the dissipation of which provides the ICM heat source. Through controlled inviscid 3-d hydrodynamic simulations, we verify that explosive AGN-like events can launch gravity waves (g-modes) into the ambient ICM which in turn decay to volume-filling turbulence. In our model, however, this process is found to be inefficient, with less than 1% of the energy injected by the AGN activity actually ending up in the turbulence of the ambient ICM. This efficiency is an order of magnitude or more too small to explain the observations of AGN-feedback in galaxy clusters and groups with short central cooling times. Atmospheres in which the g-modes are strongly trapped/confined have an even lower efficiency since, in these models, excitation of turbulence relies on the g-modes' ability to escape from the center of the cluster into the bulk ICM. Our results suggest that, if AGN-induced turbulence is indeed the mechanism by which the AGN heats the ICM core, its driving may rely on physics beyond that captured in our ideal hydrodynamic model.

It is well known that the combination of barotropic rotation and radiative equilibrium are mutually incompatible in stars. The Sun's internal rotation is far from barotropic, however, which allows at least the theoretical possibility that the Sun's thermal balance is one of radiative equilibrium in the region of the tachocline near the outer boundary of the radiative zone. We show here that (i) the constraint of radiative equilibrium leads to a straightforward ordinary differential equation for the Sun's rotation profile, and (ii) solutions of this equation can be found that, to within current levels of accuracy, closely resemble the rotation profile deduced from helioseismology. More generally, we calculate how large a baroclinic deviation from uniform rotation is required to maintain radiative equilibrium without meridional circulation throughout the bulk of the radiative zone. Very little deviation is required, well below detectability. The feasibility of radiative equilibrium for the tachocline suggests that the issue of a spreading boundary layer may be less severe than previously thought.

We derive simple and explicit expressions for the high frequency spectrum of a classical accretion disc. Both stress-free and finite stress inner boundaries are considered. A classical accretion disc spectrum with a stress-free inner boundary departs from a Wien spectrum at large $\nu$, scaling as $\nu^{2.5}$ (as opposed to $\nu^3$) times the usual exponential cut-off. If there is finite stress at the inner disc boundary, the maximum disc temperature generally occurs at this edge, even at relatively modest values of the stress. In this case, the high frequency spectrum is proportional to $\nu^2$ times the exponential cut-off. If the temperature maximum is a local hot spot, instead of an axisymmetric ring, then an interior maximum produces a $\nu^2$ prefactor while an edge maximum yields $\nu^{1.5}$. Because of beaming effects, these latter findings should pertain to a classical relativistic disc. The asymptotics are in general robust and independent of the detailed temperature profile, provided only that the liberated free energy of differential rotation is dissipated locally, and may prove useful beyond the strict domain of classical disc theory. As observations continue to improve with time, our findings suggest the possibility of using the high energy spectral component of black hole candidates as a signature prediction of classical theory, as well as an diagnostic of the stress at the inner regions of an accretion disc.

The nearly equal lunar and solar angular sizes as subtended at the Earth is generally regarded as a coincidence. This is, however, an incidental consequence of the tidal forces from these bodies being comparable. Comparable magnitudes implies strong temporal modulation, as the forcing frequencies are nearly but not precisely equal. We suggest that on the basis of paleogeographic reconstructions, in the Devonian period, when the first tetrapods appeared on land, a large tidal range would accompany these modulated tides. This would have been conducive to the formation of a network of isolated tidal pools, lending support to A.S. Romer's classic idea that the evaporation of shallow pools was an evolutionary impetus for the development of chiridian limbs in aquatic tetrapodomorphs. Romer saw this as the reason for the existence of limbs, but strong selection pressure for terrestrial navigation would have been present even if the limbs were aquatic in origin. Since even a modest difference in the Moon's angular size relative to the Sun's would lead to a qualitatively different tidal modulation, the fact that we live on a planet with a Sun and Moon of close apparent size is not entirely coincidental: it may have an anthropic basis.

In protoplanetary disks, the inner radial boundary between the MRI turbulent (`active') and MRI quiescent (`dead') zones plays an important role in models of the disk evolution and in some planet formation scenarios. In reality, this boundary is not well-defined: thermal heating from the star in a passive disk yields a transition radius close to the star (<0.1 au), whereas if the disk is already MRI active, it can self-consistently maintain the requisite temperatures out to a transition radius of roughly 1 au. Moreover, the interface may not be static; it may be highly fluctuating or else unstable. In this paper, we study a reduced model of the dynamics of the active/dead zone interface that mimics several important aspects of a real disk system. We find that MRI-transition fronts propagate inward (a `dead front' suppressing the MRI) if they are initially at the larger transition radius, or propagate outward (an `active front' igniting the MRI) if starting from the smaller transition radius. In both cases, the front stalls at a well-defined intermediate radius, where it remains in a quasi-static equilibrium. We propose that it is this new, intermediate stalling radius that functions as the true boundary between the active and dead zones in protoplanetary disks. These dynamics are likely implicated in observations of variable accretion, such as FU Ori outbursts, as well as in those planet formation theories that require the accumulation of solid material at the dead/active interface.

We examine the linear behavior of three-dimensional Lagrangian displacements in a stratified, shearing background. The isentropic and iso-rotation surfaces of the equilibrium flow are assumed to be axisymmetric, but otherwise fully two-dimensional. Three-dimensional magnetic fields are included in the perturbation equations; however the equilibrium is assumed to be well-described by purely hydrodynamic forces. The model, in principle very general, is used to study the behavior of fluid displacements in an environment resembling the solar convection zone. Some very suggestive results emerge. All but high-latitude displacements align themselves with the observed surfaces of constant angular velocity. The tendency for the angular velocity to remain constant with depth in the bulk of the convective zone, together with other critical features of the rotation profile, emerge from little more than a visual inspection of the governing equation. In the absence of a background axial angular velocity gradient, displacements exhibit no poleward bias, suggesting that solar convection "plays-off" of prexisting shear rather than creates it. We argue that baroclinic vorticity of precisely the right order is generated at the radiative/convective zone boundary due to centrifugal distortion of equipotential surfaces that is not precisely followed by isothermal surfaces. If so, many features of the Sun's internal rotation become more clear, including: i) the general appearance of the tachocline; ii) the extension of differential rotation well into the radiative zone; iii) the abrupt change of morphology of convective zone isorotation surfaces; and iv) the inability of current numerical simulations to reproduce the solar rotation profile without imposed entropy boundary conditions.

The interaction between settling of dust grains and magnetorotational instability (MRI) turbulence in protoplanetary disks is analyzed. We use a reduced system of coupled ordinary differential equations to represent the interaction between the diffusion of grains and the inhibition of the MRI. The coupled equations are styled on a Landau equation for the turbulence and a Fokker-Planck equation for the diffusion. The turbulence-grain interaction is probably most relevant near the outer edge of the disk's quiescent, or "dead" zone. Settling is most pronounced near the midplane, where a high dust concentration can self-consistently suppress the MRI. Under certain conditions, however, grains can reach high altitudes, a result of some observational interest. Finally, we show that the equilibrium solutions are linearly stable.

The isorotation contours of the solar convective zone (SCZ) show three distinct morphologies, corresponding to two boundary layers (inner and outer), and the bulk of the interior. Previous work has shown that the thermal wind equation together with informal arguments on the nature of convection in a rotating fluid could be used to deduce the shape of the isorotation surfaces in the bulk of the SCZ with great fidelity, and that the tachocline contours could also be described by relatively simple phenomenology. In this paper, we show that the form of these surfaces can be understood more broadly as a mathematical consequence of the thermal wind equation and a narrow convective shell. The analysis does not yield the angular velocity function directly, an additional surface boundary condition is required. But much can already be deduced without constructing the entire rotation profile. The mathematics may be combined with dynamical arguments put forth in previous works to the mutual benefit of each. An important element of our approach is to regard the constant angular velocity surfaces as an independent coordinate variable for what is termed the "residual entropy," a quantity that plays a key role in the equation of thermal wind balance. The difference between the dynamics of the bulk of the SCZ and the tachocline is due to a different functional form of the residual entropy in each region. We develop a unified theory for the rotational behavior of both the SCZ and the tachocline, using the solutions for the characteristics of the thermal wind equation. These characteristics are identical to the isorotation contours in the bulk of the SCZ, but the two deviate in the tachocline. The outer layer may be treated, at least descriptively, by similar mathematical techniques, but this region probably does not obey thermal wind balance.

We revisit, via a very simplified set of equations, a linear streaming instability (technically an overstability), which is present in, and potentially important for, dusty protoplanetary disks (Youdin & Goodman 2005). The goal is a better understanding of the physical origin of such instabilities, which are notoriously subtle. Rotational dynamics seem to be essential to this type of instability, which cannot be captured by one-dimensional Cartesian models. Dust `pileups' in moving pressure maxima are an important triggering mechanism of the instability, and drag feedback of dust upon the gas allows these maxima to be strengthened. Coriolis forces and the background drift counteract the effects of the pressure force.

We analyze the linear stability of a dilute, hot plasma, taking into account the effects of stratification and anisotropic thermal conduction. The work is motivated by attempts to understand the dynamics of the intracluster medium in galaxy clusters. We show that magnetic field configurations that nominally stabilize either the heat-flux driven buoyancy instability (associated with a positive thermal gradient) or the magnetothermal instability (negative thermal gradient) can lead to previously unrecognized g-mode overstabilities. The driving source of the overstability is either radiative cooling (positive temperature gradient) or the heat flux itself (negative temperature gradient). While the implications of these overstabilities have yet to be explored, we speculate that the cold fronts observed in many relaxed galaxy clusters may be related to their non-linear evolution.

In this paper, we present a model for the effects of the tachocline on the differential rotation in the solar convection zone. The mathematical technique relies on the assumption that entropy is nearly constant ("well-mixed") in isorotation surfaces both outside and within the tachocline. The resulting solutions exhibit nontrivial features that strikingly resemble the true tachocline isorotation contours in unexpected detail. This strengthens the mathematical premises of the theory. The observed rotation pattern in the tachocline shows strong quadrupolar structure, an important feature that is explicitly used in constructing our solutions. The tachocline is treated locally as an interior boundary layer of small but finite thickness, and an explicit global solution is then constructed. A dynamical link can thus be established between the internal jump in the angular velocity at the tachocline and the spread of angular velocities observed near the solar surface. In general, our results suggest that the bulk of the solar convection zone is in thermal wind balance, and that simple quadrupolar stresses, local in radius, mediate the tachocline transition from differential rotation to uniform rotation in the radiative interior.

We outline a novel linear instability that may arise in the dead-zones of protostellar disks, and possibly the fluid interiors of planets and protoplanets. In essence it is an axisymmetric buoyancy instability, but one that would not be present in a purely hydrodynamical gas. The necessary ingredients for growth include a negative radial entropy gradient (of any magnitude), weak magnetic fields, and efficient resistive diffusion (in comparison with thermal diffusion). The character of the instability is local, axisymmetric, and double-diffusive, and it attacks lengths much shorter than the resistive scale. Like the axisymmetric convective instability, it draws its energy from the negative radial entropy gradient; but by utilising the diffusing magnetic field, it can negate the stabilising influence of rotation. Its nonlinear saturated state, while not transporting appreciable angular momentum, could drive radial and vertical mixing, which may influence the temperature structure of the disk, dust dynamics and, potentially, planet formation.

Under the assumption of thermal wind balance and effective entropy mixing in constant rotation surfaces, the isorotational contours of the solar convective zone may be reproduced with great fidelity. Even at this early stage of development, this helioseismology fit may be used to put a lower bound on the midlatitude {\em radial} solar entropy gradient, which in good accord with standard mixing length theory. In this paper, we generalize this solar calculation to fully convective stars (and potentially planets), retaining the assumptions of thermal wind balance and effective entropy mixing in isorotational surfaces. It is found that each isorotation contour is of the form $R^2 = A+B\Phi(r)$, where $R$ is the radius from the rotation axis, $\Phi(r)$ is the (assumed spherical) gravitational potential, and $A$ and $B$ are constant along the contour. This result is applied to simple models of fully convective stars. Both solar-like surface rotation profiles (angular velocity decreasing toward the poles) as well as "antisolar" profiles (angular velocity increasing toward the poles) are modeled; the latter bear some suggestive resemblance to numerical simulations. We also perform exploratory studies of zonal surface flows similar to those seen in Jupiter and Saturn. In addition to providing a practical framework for understanding the results of large scale numerical simulations, our findings may also prove useful in dynamical calculations for which a simple but viable model for the background rotation profile in a convecting fluid is needed. Finally, our work bears directly on an important goal of the CoRoT program: to elucidate the internal structure of rotating, convecting stars.

This paper concerns the interaction between non-axisymmetric inertial waves and their corotation resonances in a hydrodynamical disk. Inertial waves are of interest because they can localise in resonant cavities circumscribed by Lindblad radii, and as a consequence exhibit discrete oscillation frequencies that may be observed. It is often hypothesised that these trapped eigenmodes are affiliated with the poorly understood QPO phenomenon. We demonstrate that a large class of non-axisymmetric 3D inertial waves cannot manifest as trapped normal modes. This class includes any inertial wave whose resonant cavity contains a corotation singularity. Instead, these `singular' modes constitute a continuous spectrum and, as an ensemble, are convected with the flow, giving rise to shearing waves. Lastly, we present a simple demonstration of how the corotation singularity stabilizes three-dimensional perturbations in a slender torus.

We perform a suite of simulations of cooling cores in clusters of galaxies in order to investigate the effect of the recently discovered heat flux buoyancy instability (HBI) on the evolution of cores. Our models follow the 3-dimensional magnetohydrodynamics (MHD) of cooling cluster cores and capture the effects of anisotropic heat conduction along the lines of magnetic field, but do not account for the cosmological setting of clusters or the presence of AGN. Our model clusters can be divided into three groups according to their final thermodynamical state: catastrophically collapsing cores, isothermal cores, and an intermediate group whose final state is determined by the initial configuration of magnetic field. Modeled cores that are reminiscent of real cluster cores show evolution towards thermal collapse on a time scale which is prolonged by a factor of ~2-10 compared with the zero-conduction cases. The principal effect of the HBI is to re-orient field lines to be perpendicular to the temperature gradient. Once the field has been wrapped up onto spherical surfaces surrounding the core, the core is insulated from further conductive heating (with the effective thermal conduction suppressed to less than 1/100th of the Spitzer value) and proceeds to collapse. We speculate that, in real clusters, the central AGN and possibly mergers play the role of "stirrers," periodically disrupting the azimuthal field structure and allowing thermal conduction to sporadically heat the core.

We show that the differential rotation profile of the solar convection zone, apart from inner and outer boundary layers, can be reproduced with great accu- racy 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.

The magnetohydrodynamical behavior (MHD) of accretion disks is reviewed. A detailed presentation of the fundamental MHD equations appropriate for protostellar disks is given. The combination of a weak (subthermal) magnetic field and Keplerian rotation is unstable to the magnetorotational instability (MRI), if the degree of ionization in the disk is sufficiently high. The MRI produces enhanced angular momentum and leads to a breakdown of laminar flow into turbulence. If the turbulent energy is dissipated locally, standard "$\alpha$" modeling should give a reasonable estimate for the disk structure. Because away from the central star the ionization fraction of protostellar disks is small, they are generally not in the regime of near perfect conductivity. Nonideal MHD effects are important. Of these, Ohmic dissipation and Hall electromotive forces are the most important. The presence of dust is also critical, as small interstellar scale grains absorb free charges that are needed for good magnetic coupling. On scales of AU's there may be a region near the disk midplane that is magnetically decoupled, a so-called {\em dead zone.} But the growth and settling of the grains as time evolves reduces their efficiency to absorb charge. With ionization provided by coronal X-rays from the central star (and possibly also cosmic rays), protostellar disks may be sufficiently magnetized throughout most of their lives to be MRI active, especially away from the disk midplane.

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.

This is a white paper submitted to the Stars and Stellar Evolution (SSE) Science Frontier Panel (SFP) of the NRC's Astronomy and Astrophysics 2010 Decadal Survey. The white paper is endorsed by the American Physical Society's (APS) Topical Group on Plasma Astrophysics (GPAP).

Local simulations of the magnetorotational instability (MRI) in accretion disks 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 the clean isolated channels which appear in systems near criticality or which emerge from low amplitude initial conditions. These structures inevitably achieve large amplitudes and are only then destroyed, giving rise to eruptive behaviour.

We show that a simple, modified version of the Magnetorotational Instability (MRI) can develop in the outer liquid core of the Earth, in the presence of a background shear. It requires either thermal wind, or a primary instability, such as convection, to drive a weak differential rotation within the core. The force balance in the Earth's core is very unlike classical astrophysical applications of the MRI (such as gaseous disks around stars). Here, the weak differential rotation in the Earth core yields an instability by its constructive interaction with the planet's much larger rotation rate. The resulting destabilising mechanism is just strong enough to counteract stabilizing resistive effects, and produce growth on geophysically interesting timescales. We give a simple physical explanation of the instability, and show that it relies on a force balance appropriate to the Earth's core, known as magnetostrophic balance.

The solar convective zone, or SCZ, is nearly adiabatic and marginally convectively unstable. But the SCZ is also in a state of differential rotation, and its dynamical stability properties are those of a weakly magnetized gas. This renders it far more prone to rapidly growing rotational baroclinic instabilities than a hydrodynamical system would be. These instabilities should be treated on the same footing as convective instabilites. If isentropic and isorotational surfaces coincide in the SCZ, the gas is marginally (un)stable to {\em both} convective and rotational disturbances. This is a plausible resolution for the instabilities associated with these more general rotating convective systems. This motivates an analysis of the thermal wind equation in which isentropes and isorotational surfaces are identical. The characteristics of this partial differential equation correspond to isorotation contours, and their form may be deduced even without precise knowledge of how the entropy and rotation are functionally related. Although the exact solution of the global SCZ problem in principle requires this knowledge, even the simplest models produce striking results in broad agreement with helioseismology data. This includes horizontal (i.e. quasi-spherical) isorotational contours at the poles, axial contours at the equator, and approximately radial contours at midlatitudes. The theory does not apply directly to the tachocline, where a simple thermal wind balance is not expected to be valid. The work presented here is subject to tests of self-consistency, among them the prediction that there should be good agreement between isentropes and isorotational contours in sufficiently well-resolved large scale numerical MHD simulations.

The role of thermal conduction in regulating the thermal behavior of cooling flows in galaxy clusters is reexamined. Recent investigations have shown that the anisotropic Coulomb heat flux caused by a magnetic field in a dilute plasma drives a dynamical instability. A long standing problem of cooling flow theory has been to understand how thermal conduction can offset radiative core losses without completely preventing them. In this Letter we propose that magnetohydrodynamic turbulence driven by the heat flux instability regulates field-line insulation and drives a reverse convective thermal flux, both of which may mediate the stabilization of the cooling cores of hot clusters. This model suggests that turbulent mixing should accompany strong thermal gradients in cooling flows. This prediction seems to be supported by the spatial distribution of metals in the central galaxies of clusters, which shows a much stronger correlation with the ambient hot gas temperature gradient than with the parent stellar population.

We investigate the behavior of the magnetic Prandtl number (ratio of microscopic viscosity to resistivity) for accretion sources. Generally this number is very small in standard accretion disk models, but can become larger than unity within $\sim 50$ Schwarzschild radii of the central mass. Recent numerical investigations suggest a marked dependence of the level of MHD turbulence on the value of the Prandtl number. Hence, black hole and neutron star accretors, i.e. compact X-ray sources, are affected. The astrophysical consequences of this could be significant, including a possible route to understanding the mysterious state changes that have long characterized these sources.

Axisymmetric incompressible modes of the magneto-rotational instability (MRI) with a vertical wavenumber are exact solutions of the non-linear local equations of motion for a disk (shearing box). They are referred to as "channel solutions". Here, we generalize a class of these solutions to include energy losses, viscous, and resistive effects. In the limit of zero shear, we recover the result that torsional Alfv\'en waves are exact solutions of the non-linear equations. Our method allows the extension of these solutions into the dissipative regime. These new solutions serve as benchmarks for simulations including dissipation and energy loss, and to calibrate numerical viscosity and resistivity in the Zeus3D code. We quantify the anisotropy of numerical dissipation and compute its scaling with time and space resolution. We find a strong dependence of the dissipation on the mean magnetic field that may affect the saturation state of the MRI as computed with Zeus3D. It is also shown that elongated grid cells generally preclude isotropic dissipation and that a Courant time step smaller than that which is commonly used should be taken to avoid spurious anti-diffusion of magnetic field.

We present an exact three-dimensional wave solution to the shearing sheet equations of motion. The existence of this solution argues against transient amplification as a route to turbulence in unmagnetized disks. Moreover, because the solution covers an extensive dynamical range in wavenumber space, it is an excellent test of the dissipative properties of numerical codes.

In dilute astrophysical plasmas, the collisional mean free path of a particle can exceed its Larmor radius. Under these conditions, the thermal conductivity and viscosity of the plasma can be dramatically altered. This alteration allows outwardly decreasing angular velocity or temperature gradients to become strongly destabilizing. This paper generalizes an earlier, simple analysis of the viscosity instability, by including the dynamical effects of magnetic field line tension. Such effects lower the growth rates found in the absence of such tension, but still allow growth rates in excess of the maximum of the standard magnetorotational instability. We find very good quantitative agreement with more complex kinetic treatments of the same process. The combination of large growth rates and large magnetic Prandtl number suggest that protogalactic disks are powerful dynamos.

Numerical MHD codes have become extraordinarily powerful tools with which to study accretion turbulence. They have been used primarily to extract values for the classical $\alpha$ parameter, and to follow complex evolutionary development. Energy transport, which is at the heart of classical disk theory, has yet to be explored in any detail. Further topics that should be explored by simulation include nonideal MHD, radiation physics, and outburst behavior related to the temperature sensitivity of the resistivity.

In this paper and a companion work, we report on the first global numerical simulations of self-gravitating magnetized tori, subject in particular to the influence of the magnetorotational instability (MRI). In this work, paper I, we restrict our calculations to the study of the axisymmetric evolution of such tori. Our goals are twofold: (1) to investigate how self-gravity influences the global structure and evolution of the disks; and (2) to determine whether turbulent density inhomogeneities can be enhanced by self-gravity in this regime. As in non self-gravitating models, the linear growth of the MRI is followed by a turbulent phase during which angular momentum is transported outward. As a result, self-gravitating tori quickly develop a dual structure composed of an inner thin Keplerian disk fed by a thicker self-gravitating disk, whose rotation profile is close to a Mestel disk. Our results show that the effects of self-gravity enhance density fluctuations much less than they smooth the disk, and giving it more coherence. We discuss the expected changes that will occur in 3D simulations, the results of which are presented in a companion paper.

We present 3D magnetohydrodynamic (MHD) numerical simulations of the evolution of self--gravitating and weakly magnetized disks with an adiabatic equation of state. Such disks are subject to the development of both the magnetorotational and gravitational instabilities, which transport angular momentum outward. As in previous studies, our hydrodynamical simulations show the growth of strong m=2 spiral structure. This spiral disturbance drives matter toward the central object and disappears when the Toomre parameter Q has increased well above unity. When a weak magnetic field is present as well, the magnetorotational instability grows and leads to turbulence. In that case, the strength of the gravitational stress tensor is lowered by a factor of about~2 compared to the hydrodynamical run and oscillates periodically, reaching very small values at its minimum. We attribute this behavior to the presence of a second spiral mode with higher pattern speed than the one which dominates in the hydrodynamical simulations. It is apparently excited by the high frequency motions associated with MHD turbulence. The nonlinear coupling between these two spiral modes gives rise to a stress tensor that oscillates with a frequency which is a combination of the frequencies of each of the modes. This interaction between MHD turbulence and gravitational instabilities therefore results in a smaller mass accretion rate onto the central object.

An exact nonlinear scaling transformation is presented for the local three-dimensional dynamical equations of motion for differentially rotating disks. The result is relevant to arguments that have been put forth claiming that numerical simulations lack the necessary numerical resolution to resolve nonlinear instabilities that are supposedly present. We show here that any time dependent velocity field satisfying the local equations of motion and existing on small length scales, has an exact rescaled counterpart that exists on arbitrarily larger scales as well. Large scale flows serve as a microscope to view small scale behavior. The absence of any large scale instabilities in local numerical simulations of Keplerian disks suggests that the equations in this form have no instabilities at any scale, and that finite Reynolds number suppression is not the reason for the exhibited stable behavior. While this argument does not rule out the possibility of global hydrodynamical instability, it does imply that differential rotation {\em per se} is not unstable in a manner analogous to shear layers or high Reynolds number Poiseuille flow. Analogies between the stability behavior of accretion disks and these flows are specious.

When the ion mean free path much exceeds the Larmor radius in a plasma, the viscous stress tensor is altered dramatically, and depends only upon quantities measured along the field lines. This regime corresponds to typical interstellar medium conditions in galaxies and protogalaxies, even if the magnetic field is extremely weak, with a negligible Lorentz force on all scales of interest. In this work, the only role of the magnetic field is to channel angular momentum transport along its lines of force. We show that differential rotation in such a gas is highly unstable, with a maximum growth rate exceeding that of the magnetorotational instability. The regime of interest has been treated previously by plasma kinetic methods. Where there is overlap, our work appears to be in agreement with the kinetic results. The nonlinear outcome of this instability is likely to be a turbulent process, significantly augmenting the magnetorotational instability, and important to the initial phases of the amplification of small galactic magnetic fields.

We present global 2D and 3D simulations of self-gravitating magnetized tori. We used the 2D calculations to demonstrate that the properties of the MRI are not affected by the presence of self-gravity: MHD turbulence and enhanced angular momentum transport follow the linear growth of the instability. In 3D, we have studied the interaction between an $m=2$ gravitational instability and MHD turbulence. We found its strength to be significantly decreased by the presence of the latter, showing that both instabilities strongly interact in their non-linear phases. We discuss the consequences of these results.

We study the local stability of stratified, differentially-rotating fluids to axisymmetric perturbations in the presence of a weak magnetic field and of finite resistivity, viscosity and heat conductivity. This is a generalization of the Goldreich-Schubert-Fricke (GSF) double-diffusive analysis to the magnetized and resistive, triple-diffusive case. Our fifth-order dispersion relation admits a novel branch which describes a magnetized version of multi-diffusive modes. We derive necessary conditions for axisymmetric stability in the inviscid and perfect-conductor (double-diffusive) limits. In each case, rotation must be constant on cylinders and angular velocity must not decrease with distance from the rotation axis for stability, irrespective of the relative strength of viscous, resistive and heat diffusion. Therefore, in both double-diffusive limits, solid body rotation marginally satisfies our stability criteria. The role of weak magnetic fields is essential to reach these conclusions. The triple-diffusive situation is more complex, and its stability criteria are not easily stated. Numerical analysis of our general dispersion relation confirms our analytic double-diffusive criteria, but also shows that an unstable double-diffusive situation can be significantly stabilized by the addition of a third, ostensibly weaker, diffusion process. We describe a numerical application to the Sun's upper radiative zone and establish that it would be subject to unstable multi-diffusive modes if moderate or strong radial gradients of angular velocity were present.

We present the first global simulations of self-gravitating magnetized tori. The simulations are performed with Zeus-2D and GLOBAL. We find the magnetorotational instability (MRI) to behave similarly in a self-gravitating environment as in previous simulations of non self-gravitating systems: enhancement of turbulent angular momentum transport follows the linear phase. The torus quickly develops a two component structure composed of an inner thick disk in Keplerian rotation and an outer massive disk. We compare this result with zero mass global simulations in 2D, and also present preliminary results of 3D simulations.

The effects of ambipolar diffusion on the linear stability of weakly ionised accretion discs are examined. Earlier work on this topic has focused on axial magnetic fields and perturbation wavenumbers. We consider here more general field and wavenumber geometries, and find that qualitatively new results are obtained. Provided a radial wavenumber and azimuthal field are present along with their axial counterparts, ambipolar diffusion will always be destabilising, with unstable local modes appearing at well-defined wavenumber bands. The wavenumber corresponding to the maximum growth rate need not, in general, lie along the vertical axis. Growth rates become small relative to the local angular velocity when the ion-neutral collision time exceeds the orbital time. In common with Hall electromotive forces, ambipolar diffusion destabilises both positive and negative angular velocity gradients. In at least some cases, therefore, uniformly rotating molecular cloud cores may reflect the marginally stable state of the ambipolar magnetorotational instability.

Just as correlations between fluctuating radial and azimuthal velocities produce a coherent stress contributing to the angular momentum transport in turbulent accretion disks, correlations in the velocity and temperature fluctuations produce a coherent energy flux. This nonadvective energy flux is always of secondary importance in thin radiative disks, but cannot be neglected in nonradiative flows, in which it completes the mean field description of turbulence. It is, nevertheless, generally ignored in accretion flow theory, with the exception of models explicitly driven by thermal convection, where it is modeled phenomenologically. This flux embodies both turbulent thermal convection as well as wave transport, and its presence is essential for a proper formulation of energy conservation, whether convection is present or not. The sign of the thermal flux is likely to be outward in real systems, but the restrictive assumptions used in numerical simulations may lead to inward thermal transport, in which case qualitatively new effects may be exhibited. We find, for example, that a static solution would require inward, not outward, thermal transport. Even if it were present, thermal convection would be unlikely to stifle accretion, but would simply add to the outward rotational energy flux that must already be present.

The status of our current understanding of angular momentum transport in accretion disks is reviewed. The last decade has seen a dramatic increase both in the recognition of key physical processes and in our ability to carry through direct numerical simulations of turbulent flow. Magnetic fields have at once powerful and subtle influences on the behavior of (sufficiently) ionized gas, rendering them directly unstable to free energy gradients. Outwardly decreasing angular velocity profiles are unstable. The breakdown of Keplerian rotation into MHD turbulence may be studied in some numerical detail, and key transport coefficients may be evaluated. Chandra observations of the Galactic Center support the existence of low luminosity accretion, which may ultimately prove amenable to global three-dimensional numerical simulation.

The processes of planet formation and migration depend intimately on the interaction between planetesimals and the gaseous disks in which they form. The formation of gaps in the disk can severely limit the mass of the planet and its migration toward the protostar. We investigate the process of gap formation through magnetohydrodynamic simulations in which internal stress arises self-consistently from turbulence generated by the magnetorotational instability. The simulations investigate three different planetary masses and two disk temperatures to bracket the tidal (thermal) and viscous gap opening conditions. The results are in general qualitative agreement with previous simulations of gap formation, but show significant differences. In the presence of MHD turbulence, the gaps produced are shallower and asymmetrically wider than those produced with pure hydrodynamics. The rate of gap formation is also slowed, with accretion occurring across the developing gap. Viscous hydrodynamics does not adequately describe the evolution, however, because planets capable of producing gaps also may be capable of affecting the level MHD turbulence in different regions of the disk.

Chaotic flow is studied in a series of numerical magnetohydrodynamical simulations that use the shearing box formalism. This mimics important features of local accretion disk dynamics. The magnetorotational instability gives rise to flow turbulence, and quantitative chaos parameters, such as the largest Lyapunov exponent, can be measured. Linear growth rates appear in these exponents even when the flow is fully turbulent. The extreme sensitivity to initial conditions associated with chaotic flows has practical implications, the most important of which is that hundreds of orbital times are needed to extract a meaningful average for the stress. If the evolution time in a disk is less than this, the classical $\alpha$ formalism will break down.

The authors respond to the commentary listed in the title.

A simple variational calculation is presented showing that a uniformly rotating barotropic fluid in an external potential attains a true energy minimum if and only if the rotation profile is everywhere subsonic. If regions of supersonic rotation are present, fluid variations exist that could take the sytem to states of lower energy. In any given system, these states may or may not be dynamically accessible, but their existence is important. It means that extending the degrees of freedom available to the fluid (say by weak magnetic fields) may open a path to fluid instabilities. Whether astrophysical gaseous nebula tend toward states of uniform rotation or toward more Keplerian core-disk systems appears to be largely a matter of whether the rotation profile is transonic or not. The suggestion is made that the length scale associated with coherent molecular cloud cores is related to the requirement that the cores be stable and rotate subsonically.

We review numerical simulations of MHD turbulence. The last decade has witnessed fundamental advances both in the technical capabilities of direct numerical simulation, and in our understanding of key physical processes. Magnetic fields tap directly into the free energy sources in a sufficiently ionized gas. The result is that adverse angular velocity and adverse temperature gradients, not the classical angular momentum and entropy gradients, destabilize laminar and stratified flow. This has profound consequences for astrophysical accretion flows, and has opened the door to a new era of numerical simulation experiments.}

We analyze three-dimensional magnetohydrodynamic (MHD) simulations of a nonradiative accretion flow around a black hole using a pseudo-Newtonian potential. The flow originates from a torus initially centered at 100 gravitational (Schwarzschild) radii. Accretion is driven by turbulent stresses generated self-consistently by the magnetorotational instability. The resulting flow has three well-defined dynamical components: a hot, thick, rotationally-dominated Keplerian disk; a surrounding magnetized corona with vigorous circulation and outflow; and a magnetically-confined jet along the centrifugal funnel wall. Inside of 10 gravitational radii, the disk becomes very hot, more toroidal, and highly intermittent. These results contrast sharply with quasi-spherical, self-similar viscous models. There are no significant dynamical differences between simulations that include resistive heating and those that do not. We conclude by deducing some simple radiative properties of our solutions, and apply the results to the accretion-powered Galactic center source Sgr A*.

The principles underlying a proposed class of black hole accretion models are examined. The flows are generally referred to as ``convection-dominated,'' and are characterized by inward transport of angular momentum by thermal convection and outward viscous transport, vanishing mass accretion, and vanishing local energy dissipation. We examine the viability of these ideas by explicitly calculating the leading order angular momentum transport of axisymmetric modes in magnetized, differentially rotating, stratified flows. The modes are destabilized by the generalized magnetorotational instability. There are no inward transporting modes at all unless the magnitude of the (imaginary) Brunt frequency is comparable to the epicyclic frequency, a condition requiring substantial levels of dissipation. When inward transporting modes do exist, they appear at long wavelengths, unencumbered by magnetic tension. Moreover, very general thermodynamic principles prohibit the complete recovery of irreversible dissipative energy losses, a central feature of convection-dominated models. Dissipation of the free energy of differential rotation inevitably requires outward angular momentum transport. Our results are in good agreement with global MHD simulations, which find significant levels of outward transport and energy dissipation, whether or not destabilizing entropy gradients are present.

We calculate the ionization fraction of protostellar alpha disks, taking into account vertical temperature structure, and the possible presence of trace metal atoms. Both thermal and X-ray ionization are considered. Previous investigations of layered disks used radial power-law models with isothermal vertical structure. But alpha models are used to model accretion, and the present work is a step towards a self-consistent treatment. The extent of the magnetically uncoupled (``dead'') zone depends sensitively on alpha, on the assumed accretion rate, and on the critical magnetic Reynolds number, below which MHD turbulence cannot be self-sustained. Its extent is extremely model-dependent. It is also shown that a tiny fraction of the cosmic abundance of metal atoms can dramatically affect the ionization balance. Gravitational instabilities are an unpromising source of transport, except in the early stages of disk formation.

The stability of a dilute plasma to local convective and rotational disturbances is examined. A subthermal magnetic field and finite thermal conductivity along the field lines are included in the analysis. Stability criteria similar in form to the classical H{\o}iland inequalities are found, but with angular velocity gradients replacing angular momentum gradients, and temperature gradients replacing entropy gradients. These criteria are indifferent to the properties of the magnetic field and to the magnitude of the thermal conductivity. Angular velocity gradients and temperature gradients are both free energy sources; it is not surprising that they are directly relevant to the stability of the gas. Magnetic fields and thermal conductivity provide the means by which these sources can be tapped. Previous studies have generally been based upon the classical H{\o}iland criteria, which are inappropriate for magnetized, dilute astrophysical plasmas. In sharp contrast to recent claims in the literature, the new stability criteria demonstrate that marginal flow stability is not a fundamental property of accreting plasmas thought to be associated with low luminosity X-ray sources.

We present a global magnetohydrodynamic (MHD) three dimensional simulation of a nonradiative accretion flow originating in a pressure supported torus. The evolution is controlled by the magnetorotational instability which produces turbulence. The flow forms a nearly Keplerian disk. The total pressure scale height in this disk is comparable to the vertical size of the initial torus. Gas pressure dominates only near the equator; magnetic pressure is more important in the surrounding atmosphere. A magnetically dominated bound outflow is driven from the disk. The accretion rate through the disk exceeds the final rate into the hole, and a hot torus forms inside 10 r_g. Hot gas, pushed up against the centrifugal barrier and confined by magnetic pressure, is ejected in a narrow, unbound, conical outflow. The dynamics are controlled by magnetic turbulence, not thermal convection, and a hydrodynamic alpha model is inadequate to describe the flow. The limitations of two dimensional MHD simulations are also discussed.

The effects of Hall electromotive forces (HEMFs) on the linear stability of protostellar disks are examined. Earlier work on this topic focused on axial field and perturbation wavenumbers. Here we treat the problem more generally. Both axisymmetric and nonaxisymmetric cases are investigated. Though seldom explicitly included in calculations, HEMFs appear to be important whenever Ohmic dissipation is. They allow for the appearance of electron whistler waves, and since these have right-handed polarization, a helicity factor is also introduced into the stability problem. This factor is the product of the components of the angular velocity and magnetic field along the perturbation wavenumber, and it is destabilizing when negative. Unless the field and angular velocity are exactly aligned, it is always possible to find destabilizing wavenumbers. HEMFs can destabilize any differential rotation law, even those with angular velocity increasing outward. Regardless of the sign of the angular velocity gradient, the maximum growth rate is always given in magnitude by the local Oort A value of the disk, as in the standard magnetorotational instability. The role of Hall EMFs may prove crucial to understanding how turbulence is maintained in the ``low state'' of eruptive disk systems.

The stratification of entropy and the stratification of angular momentum are closely analogous. Of particular interest is the behavior of disks in which angular momentum transport is controlled by convection, and heat transport by dynamical turbulence. In both instances we argue that the transport must proceed ``backwards'' relative to the sense one would expect from a simple enhanced diffusion approach. Reversed angular momentum transport has already been seen in numerical simulations; contra-gradient thermal diffusion should be amenable to numerical verification as well. These arguments also bear on the observed nonlinear local stability of isolated Keplerian disks. We also describe a diffusive instability that is the entropy analogue to the magnetorotational instability. It affects thermally stratified layers when Coulomb conduction and a weak magnetic field are present. The criterion for convective instability goes from one of upwardly decreasing entropy to one of upwardly decreasing temperature. The maximum growth rate is of order the inverse sound crossing time, independent of the thermal conductivity. The indifference of the growth rate to the conduction coefficient, its simple dynamical scaling, and the replacement in the stability criterion of a conserved quantity (entropy) gradient by a free energy (temperature) gradient are properties similar to those exhibited by the magnetorotational instability.

The dynamical state of the solar nebula depends critically upon whether or not the gas is magnetically coupled. The presence of a subthermal field will cause laminar flow to break down into turbulence. Magnetic coupling, in turn, depends upon the ionization fraction of the gas. The inner most region of the nebula ($\lesssim 0.1$ AU) is magnetically well-coupled, as is the outermost region ($\gtrsim 10$ AU). The magnetic status of intermediate scales ($\sim 1$ AU) is less certain. It is plausible that there is a zone adjacent to the inner disk in which turbulent heating self-consistently maintains the requisite ionization levels. But the region adjacent to the active outer disk is likely to be magnetically ``dead.'' Hall currents play a significant role in nebular magnetohydrodynamics. Though still occasionally argued in the literature, there is simply no evidence to support the once standard claim that differential rotation in a Keplerian disk is prone to break down into shear turbulence by nonlinear instabilities. There is abundant evidence---numerical, experimental, and analytic---in support of the stabilizing role of Coriolis forces. Hydrodynamical turbulence is almost certainly not a source of enhanced turbulence in the solar nebula, or in any other astrophysical accretion disk.

The dynamical foundations of $\alpha$ disk models are described. At the heart of the viscous formalism of accretion disk models are correlations in the fluctuating components of the disk velocity, magnetic field, and gravitational potential. We relate these correlations to the large scale mean flow dynamics used in phenomenological viscous disk models. MHD turbulence readily lends itself to the $\alpha$ formalsim, but transport by self-gravity does not. Nonlocal transport is an intrinsic property of turbulent self-gravitating disks, which in general cannot be captured by an $\alpha$ model. Local energy dissipation and $\alpha$-like behavior can be re-established if the pattern speeds associated with the amplitudes of an azimuthal Fourier decomposition of the turbulence are everywhere close to the local rotation frequency. In this situation, global wave transport must be absent. Shearing box simulations, which employ boundary conditions forcing local behavior, are probably not an adequate tool for modeling the behavior of self-gravitating disks. As a matter of principle, it is possible that disks which hover near the edge of gravitational stability may behave in accord with a local $\alpha$ model, but global simulations performed to date suggest matters are not this simple.

We employ a variety of numerical simulations in the local shearing box system to investigate in greater depth the local hydrodynamic stability of Keplerian differential rotation. In particular we explore the relationship of Keplerian shear to the nonlinear instabilities known to exist in simple Cartesian shear. The Coriolis force is the source of linear stabilization in differential rotation. We exploit the formal equivalence of constant angular momentum flows and simple Cartesian shear to examine the transition from stability to nonlinear instability. The manifestation of nonlinear instability in shear flows is known to be sensitive to initial perturbation and to the amount of viscosity; marginally (linearly) stable differentially rotating flows exhibit this same sensitivity. Keplerian systems, however, are completely stable; the strength of the stabilizing Coriolis force easily overwhelms any destabilizing nonlinear effects. In fact, nonlinear effects speed the decay of applied turbulence by producing a rapid cascade of energy to high wavenumbers where dissipation occurs. Our conclusions are tested with grid resolution experiments and by comparison with results from a code that employs an alternative numerical algorithm. The properties of hydrodynamic differential rotation are contrasted with magnetohydrodynamic differential rotation. The kinetic stress couples to the vorticity which limits turbulence, while the magnetic stress couples to the shear which promotes turbulence. Thus magnetohydrodynamic turbulence is uniquely capable of acting as a turbulent angular momentum transport mechanism in disks.

A new class of disk MHD equilibrium solutions is described, which is valid within the standard local (``shearing sheet'') approximation scheme. These solutions have the following remarkable property: velocity streamlines and magnetic lines of force rotate rigidly, even in the presence of differential rotation. This situation comes about because the Lorentz forces acting upon modified epicycles compel fluid elements to follow magnetic lines of force. Field line (and streamline) configurations may be elliptical or hyperbolic, prograde or retrograde. These structures have previously known hydrodynamical analogs: the ``planet'' solutions described by Goodman, Narayan, & Goldreich. The primary focus of this investigation is configurations in the disk plane. A related family of solutions lying in a vertical plane is briefly discussed; other families of solutions may exist. Whether these MHD structures are stable is not yet known, but could readily be determined by three-dimensional simulations. If stable or quasi-stable, these simple structures may find important applications in both accretion and galactic disks.

When present, extended disks of neutral hydrogen around spiral galaxies show a remarkably uniform velocity dispersion of approx 6 km/s. Since stellar winds and supernovae are largely absent in such regions, neither the magnitude nor the constancy of this number can be accounted for in the classical picture in which interstellar turbulence is driven by stellar energy sources. Here we suggest that magnetic fields with strengths of a few microgauss in these extended disks allow energy to be extracted from galactic differential rotation through MHD driven turbulence. The magnitude and constancy of the observed velocity dispersion may be understood if its value is Alfvenic. Moreover, by providing a simple explanation for a lower bound to the gaseous velocity fluctuations, MHD processes may account for the sharp outer edge to star formation in galaxy disks.