Normalized to: Latter, H.

**57** article(s) in total. **36** co-authors, from **1** to **18** common article(s). Median position in authors list is **2,0**.

[1] oai:arXiv.org:2003.13263 [pdf] - 2090658

Electric heating and angular momentum transport in laminar models of
protoplanetary disks

Comments: 19 pages, 21 figures, published in MNRAS; the previous version used
an incorrect gas density to compute the opacity (2.5x lower than indicated)
and ionization fraction

Submitted: **2020-03-30**, last modified: **2020-05-06**

The vertical temperature structure of a protoplanetary disk bears on several
processes relevant to planet formation, such as gas and dust grain chemistry,
ice lines and convection. The temperature profile is controlled by irradiation
from the central star and by any internal source of heat such as might arise
from gas accretion. We investigate the heat and angular momentum transport
generated by the resistive dissipation of magnetic fields in laminar disks. We
use local one-dimensional simulations to obtain vertical temperature profiles
for typical conditions in the inner disk (0.5 to 4 au). Using simple
assumptions for the gas ionization and opacity, the heating and cooling rates
are computed self-consistently in the framework of radiative non-ideal
magnetohydrodynamics. We characterize steady solutions that are symmetric about
the midplane and which may be associated with saturated Hall-shear unstable
modes. We also examine the dissipation of electric currents driven by global
accretion-ejection structures. In both cases we obtain significant heating for
a sufficiently high opacity. Strong magnetic fields can induce an order-unity
temperature increase in the disk midplane, a convectively unstable entropy
profile, and a surface emissivity equivalent to a viscous heating of $\alpha
\sim 10^{-2}$. These results show how magnetic fields may drive efficient
accretion and heating in weakly ionized disks where turbulence might be
inefficient, at least for a range of radii and ages of the disk.

[2] oai:arXiv.org:2003.01167 [pdf] - 2061792

Dust dynamics and vertical settling in gravitoturbulent protoplanetary
discs

Comments: Accepted in MNRAS, 13 pages

Submitted: **2020-03-02**

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

[3] oai:arXiv.org:2001.08693 [pdf] - 2069013

Global simulations of self-gravitating magnetized protoplanetary disks

Comments: Happy Chinese new year

Submitted: **2020-01-23**

In the early stages of a protoplanetary disk, when its mass is a significant
fraction of its star's, turbulence generated by gravitational instability (GI)
should feature significantly in the disk's evolution. At the same time, the
disk may be sufficiently ionised for magnetic fields to play some role in the
dynamics. Though usually neglected, the impact of magnetism on the GI may be
critical, with consequences for several processes: the efficiency of accretion,
spiral structure formation, fragmentation, and the dynamics of solids. In this
paper, we report on global three-dimensional magnetohydrodynamical simulations
of a self-gravitating protoplanetary disk using the meshless finite mass (MFM)
Lagrangian technique. We confirm that GI spiral waves trigger a dynamo that
amplifies an initial magnetic field to nearly thermal amplitudes (plasma beta <
10), an order of magnitude greater than that generated by the
magneto-rotational instability alone. We also determine the dynamo's nonlinear
back reaction on the gravitoturbulent flow: the saturated state is
substantially hotter, with an associated larger Toomre parameter and weaker,
more 'flocculent' spirals. But perhaps of greater import is the dynamo's
boosting of accretion via a significant Maxwell stress; mass accretion is
enhanced by factors of several relative to either pure GI or pure MRI. Our
simulations use ideal MHD, an admittedly poor approximation in protoplanetary
disks, and thus future studies should explore the full gamut of non-ideal MHD.
In preparation for that, we exhibit a small number of Ohmic runs that reveal
that the dynamo, if anything, is stronger in a non-ideal environment. This work
confirms that magnetic fields are a potentially critical ingredient in
gravitoturbulent young disks, possibly controlling their evolution, especially
via their enhancement of (potentially episodic) accretion.

[4] oai:arXiv.org:1910.10696 [pdf] - 1985418

HFQPOs and discoseismic mode excitation in eccentric, relativistic
discs. I. Hydrodynamical simulations

Comments: 14 pages, 13 figures. Submitted to MNRAS on July 17, 2019 (no report
yet received). Comments welcome

Submitted: **2019-10-23**

High-frequency quasi-periodic oscillations (HFQPOs) observed in the emission
of black hole X-ray binary systems promise insight into strongly curved
spacetime. `Discoseismic' modes with frequencies set by the intrinsic
properties of the central black hole, in particular `trapped inertial waves'
(r-modes), offer an attractive explanation for HFQPOs. To produce an observable
signature, however, such oscillations must be excited to sufficiently large
amplitudes. Turbulence driven by the magnetorotational instability (MRI) does
not appear to provide the necessary amplification, but r-modes may still be
excited via interaction with accretion disc warps or eccentricities. We present
global, hydrodynamic simulations of relativistic accretion discs, which
demonstrate for the first time the excitation of trapped inertial waves by an
imposed eccentricity in the flow. While the r-modes' saturated state depends on
the vertical boundary conditions used in our unstratified, cylindrical
framework, their excitation is unambiguous in all runs with eccentricity >0.005
near the ISCO. These simulations provide a proof of concept, demonstrating the
robustness of the trapped inertial wave excitation mechanism in a
non-magnetized context. We explore the competition between this excitation, and
damping by MHD turbulence and radial inflow in a companion paper.

[5] oai:arXiv.org:1901.05190 [pdf] - 1871477

Local simulations of MRI turbulence with meshless methods

Comments: This version has been accepted by ApJS

Submitted: **2019-01-16**, last modified: **2019-02-21**

The magneto-rotational instability (MRI) is one of the most important
processes in sufficiently ionized astrophysical disks. Grid-based simulations,
especially those using the local shearing box approximation, provide a powerful
tool to study the ensuing nonlinear turbulence. On the other hand, while
meshless methods have been widely used in both cosmology, galactic dynamics,
and planet formation they have not been fully deployed on the MRI problem. We
present local unstratified and vertically stratified MRI simulations with two
meshless MHD schemes: a recent implementation of SPH MHD (Price2012), and a MFM
MHD scheme with a constrained gradient divergence cleaning scheme, as
implemented in the GIZMO code \citep{Hopkins2017}. Concerning variants of the
SPH hydro force formulation we consider both the "vanilla" SPH and the PSPH
variant included in GIZMO. We find, as expected, that the numerical noise
inherent in these schemes affects turbulence significantly. A high order
kernel, free of the pairing instability, is necessary. Both schemes can
adequately simulate MRI turbulence in unstratified shearing boxes with net
vertical flux. The turbulence, however, dies out in zero-net-flux unstratified
boxes, probably due to excessive and numerical dissipation. In zero-net-flux
vertically stratified simulations, MFM can reproduce the MRI dynamo and its
characteristic butterfly diagram for several tens of orbits before ultimately
decaying. In contrast, extremely strong toroidal fields, as opposed to
sustained turbulence, develop in equivalent simulations using SPH MHD. This
unphysical state in SPH MHD is likely caused by a combination of excessive
artificial viscosity, numerical resistivity, and the relatively large residual
errors in the divergence of the magnetic field remaining even after cleaning
procedures are applied.

[6] oai:arXiv.org:1811.03083 [pdf] - 1781655

Quasi-periodic oscillations, trapped inertial waves and strong toroidal
magnetic fields in relativistic accretion discs

Comments: 16 pages, 6 figures, MNRAS accepted

Submitted: **2018-11-07**

The excitation of trapped inertial waves (r-modes) by warps and
eccentricities in the inner regions of a black hole accretion disc may explain
the high-frequency quasi-periodic oscillations (HFQPOs) observed in the
emission of Galactic X-ray binaries. However, it has been suggested that strong
vertical magnetic fields push the oscillations' trapping region toward the
innermost stable circular orbit (ISCO), where conditions could be unfavourable
for their excitation. This paper explores the effects of large-scale magnetic
fields that exhibit \textit{both} toroidal and vertical components, through
local and global linear analyses. We find that a strong toroidal magnetic field
can reduce the detrimental effects of a vertical field: in fact, the isolation
of the trapping region from the ISCO may be restored by toroidal magnetic
fields approaching thermal strengths. The toroidal field couples the r-modes to
the disc's magneto-acoustic response and inflates the effective pressure within
the oscillations. As a consequence, the restoring force associated with the
vertical magnetic field's tension is reduced. Given the analytical and
numerical evidence that accretion discs threaded by poloidal magnetic field
lines develop a strong toroidal component, our result provides further evidence
that the detrimental effects of magnetic fields on trapped inertial modes are
not as great as previously thought.

[7] oai:arXiv.org:1810.06589 [pdf] - 1771739

Gravitoturbulent dynamos in astrophysical discs

Comments: accepted in MNRAS in October 2018, 21 pages

Submitted: **2018-10-15**

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

[8] oai:arXiv.org:1810.02773 [pdf] - 1767743

On Kelvin-Helmholtz and parametric instabilities driven by coronal waves

Comments: 12 pages, 4 figures, accepted for publication MNRAS

Submitted: **2018-10-05**

The Kelvin-Helmholtz instability has been proposed as a mechanism to extract
energy from magnetohydrodynamic (MHD) kink waves in flux tubes, and to drive
dissipation of this wave energy through turbulence. It is therefore a
potentially important process in heating the solar corona. However, it is
unclear how the instability is influenced by the oscillatory shear flow
associated with an MHD wave. We investigate the linear stability of a
discontinuous oscillatory shear flow in the presence of a horizontal magnetic
field within a Cartesian framework that captures the essential features of MHD
oscillations in flux tubes. We derive a Mathieu equation for the Lagrangian
displacement of the interface and analyse its properties, identifying two
different instabilities: a Kelvin-Helmholtz instability and a parametric
instability involving resonance between the oscillatory shear flow and two
surface Alfv\'{e}n waves. The latter occurs when the system is Kelvin-Helmholtz
stable, thus favouring modes that vary along the flux tube, and as a
consequence provides an important and additional mechanism to extract energy.
When applied to flows with the characteristic properties of kink waves in the
solar corona, both instabilities can grow, with the parametric instability
capable of generating smaller scale disturbances along the magnetic field than
possible via the Kelvin-Helmholtz instability. The characteristic time-scale
for these instabilities is $\sim 100$ s, for wavelengths of $200$ km. The
parametric instability is more likely to occur for smaller density contrasts
and larger velocity shears, making its development more likely on coronal loops
than on prominence threads.

[9] oai:arXiv.org:1808.00267 [pdf] - 1729217

Hydrodynamic convection in accretion discs

Comments: Accepted for publication in MNRAS (20 pages, 16 figures, 4 tables)

Submitted: **2018-08-01**

The prevalence and consequences of convection perpendicular to the plane of
accretion discs have been discussed for several decades. Recent simulations
combining convection and the magnetorotational instability have given fresh
impetus to the debate, as the interplay of the two processes can enhance
angular momentum transport, at least in the optically thick outburst stage of
dwarf novae. In this paper we seek to isolate and understand the most generic
features of disc convection, and so undertake its study in purely
hydrodynamical models. First, we investigate the linear phase of the
instability, obtaining estimates of the growth rates both semi-analytically,
using one-dimensional spectral computations, as well as analytically, using
WKBJ methods. Next we perform three-dimensional, vertically stratified,
shearing box simulations with the conservative, finite-volume code PLUTO, both
with and without explicit diffusion coefficients. We find that hydrodynamic
convection can, in general, drive outward angular momentum transport, a result
that we confirm with ATHENA, an alternative finite-volume code. Moreover, we
establish that the sign of the angular momentum flux is sensitive to the
diffusivity of the numerical scheme. Finally, in sustained convection, whereby
the system is continuously forced to an unstable state, we observe the
formation of various coherent structures, including large- scale and
oscillatory convective cells, zonal flows, and small vortices.

[10] oai:arXiv.org:1803.08757 [pdf] - 1663412

Dissipative structures in magnetorotational turbulence

Comments: MNRAS, accepted

Submitted: **2018-03-23**

Via the process of accretion, magnetorotational turbulence removes energy
from a disk's orbital motion and transforms it into heat. Turbulent heating is
far from uniform and is usually concentrated in small regions of intense
dissipation, characterised by abrupt magnetic reconnection and higher
temperatures. These regions are of interest because they might generate
non-thermal emission, in the form of flares and energetic particles, or
thermally process solids in protoplanetary disks. Moreover, the nature of the
dissipation bears on the fundamental dynamics of the magnetorotational
instability (MRI) itself: local simulations indicate that the large-scale
properties of the turbulence (e.g. saturation levels, the stress-pressure
relationship) depend on the short dissipative scales. In this paper we
undertake a numerical study of how the MRI dissipates and the small-scale
dissipative structures it employs to do so. We use the Godunov code RAMSES and
unstratified compressible shearing boxes. Our simulations reveal that
dissipation is concentrated in ribbons of strong magnetic reconnection that are
significantly elongated in azimuth, up to a scale height. Dissipative
structures are hence meso-scale objects, and potentially provide a route by
which large scales and small scales interact. We go on to show how these
ribbons evolve over time --- forming, merging, breaking apart, and
disappearing. Finally, we reveal important couplings between the large-scale
density waves generated by the MRI and the small-scale structures, which may
illuminate the stress-pressure relationship in MRI turbulence.

[11] oai:arXiv.org:1802.06620 [pdf] - 1652489

Spiral density waves and vertical circulation in protoplanetary discs

Comments: accepted in MNRAS, 13 pages, 10 figures

Submitted: **2018-02-19**

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

[12] oai:arXiv.org:1802.04792 [pdf] - 1781558

Quasi-periodic oscillations and the global modes of relativistic, MHD
accretion discs

Comments: 20 pages, 17 figures, MNRAS accepted

Submitted: **2018-02-13**

The high-frequency quasi-periodic oscillations (HFQPOs) that punctuate the
light curves of X-ray binary systems present a window onto the intrinsic
properties of stellar-mass black holes and hence a testbed for general
relativity. One explanation for these features is that relativistic distortion
of the accretion disc's differential rotation creates a trapping region in
which inertial waves (r-modes) might grow to observable amplitudes. Local
analyses, however, predict that large-scale magnetic fields push this trapping
region to the inner disc edge, where conditions may be unfavorable for r-mode
growth. We revisit this problem from a pseudo-Newtonian but fully global
perspective, deriving linearized equations describing a relativistic,
magnetized accretion flow, and calculating normal modes with and without
vertical density stratification. In an unstratified model, the choice of
vertical wavenumber determines the extent to which vertical magnetic fields
drive the r-modes toward the inner edge. In a global model fully incorporating
density stratification, we confirm that this susceptibility to magnetic fields
depends on disc thickness. Our calculations suggest that in thin discs, r-modes
may remain independent of the inner disc edge for vertical magnetic fields with
plasma betas as low as $\beta\approx 100-300$. We posit that the appearance of
r-modes in observations may be more determined by a competition between
excitation and damping mechanisms near the ISCO than the modification of the
trapping region by magnetic fields.

[13] oai:arXiv.org:1711.07436 [pdf] - 1591476

Vortices and the saturation of the vertical shear instability in
protoplanetary disks

Comments: 16 pages, MNRAS accepted

Submitted: **2017-11-20**

If sufficiently irradiated by its central star, a protoplanetary disks falls
into an equilibrium state exhibiting vertical shear. This state may be subject
to a hydrodynamical instability, the `vertical shear instability' (VSI), whose
breakdown into turbulence transports a moderate amount of angular momentum
while also facilitating planet formation, possibly via the production of
small-scale vortices. In this paper, we show that VSI modes (a) exhibit
arbitrary spatial profiles and (b) remain nonlinear solutions to the
incompressible local equations, no matter their amplitude. The modes are
themselves subject to parasitic Kelvin-Helmholtz instability, though the disk
rotation significantly impedes the parasites and permits the VSI to attain
large amplitudes (fluid velocities < 10% the sound speed). This `delay' in
saturation probably explains the prominence of the VSI linear modes in global
simulations. More generally, the parasites may set the amplitude of VSI
turbulence in strongly irradiated disks. They are also important in breaking
the axisymmetry of the flow, via the unavoidable formation of vortices. The
vortices, however, are not aligned with the orbital plane and thus express a
pronounced $z$-dependence. We also briefly demonstrate that the vertical shear
has little effect on the magnetorotational instability, whereas magnetic fields
easily quench the VSI, a potential issue in the ionised surface regions of the
disk and also at larger radii.

[14] oai:arXiv.org:1709.06845 [pdf] - 1608468

Magnetorotational instability and dynamo action in gravitoturbulent
astrophysical discs

Comments: 23 pages, 16 figures, accepted in MNRAS

Submitted: **2017-09-20**

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

[15] oai:arXiv.org:1708.09285 [pdf] - 1587687

Local models of astrophysical discs

Comments: MNRAS, accepted

Submitted: **2017-08-30**

Local models of gaseous accretion discs have been successfully employed for
decades to describe an assortment of small-scale phenomena, from instabilities
and turbulence, to dust dynamics and planet formation. For the most part, they
have been derived in a physically motivated but essentially ad hoc fashion,
with some of the mathematical assumptions never made explicit nor checked for
consistency. This approach is susceptible to error, and it is easy to derive
local models that support spurious instabilities or fail to conserve key
quantities. In this paper we present rigorous derivations, based on an
asympototic ordering, and formulate a hierarchy of local models
(incompressible, Boussinesq, and compressible), making clear which is best
suited for a particular flow or phenomenon while spelling out explicitly the
assumptions and approximations of each. We also discuss the merits of the
anelastic approximation, emphasising that anelastic systems struggle to
conserve energy unless strong restrictions are imposed on the flow. The
problems encountered by the anelastic approximation are exacerbated by the
disk's differential rotation, but also attend non-rotating systems such as
stellar interiors. We conclude with a defence of local models and their
continued utility in astrophysical research.

[16] oai:arXiv.org:1706.06537 [pdf] - 1584929

Gravito-turbulence and the excitation of small-scale parametric
instability in astrophysical discs

Comments: 22 pages, 19 figures, accepted for publication in MNRAS

Submitted: **2017-06-20**

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

[17] oai:arXiv.org:1704.05926 [pdf] - 1582694

Turbulent fluctuations and the excitation of Z Cam outbursts

Comments: Accepted for publication in MNRAS, 16 pages, 15 figures

Submitted: **2017-04-19**

Z Cam variables are a subclass of dwarf nova that lie near a global
bifurcation between outbursting ("limit cycle") and non-outbursting
("standstill") states. It is believed that variations in the secondary star's
mass-injection rate instigate transitions between the two regimes. In this
paper we explore an alternative trigger for these transitions: stochastic
fluctuations in the disk's turbulent viscosity. We employ simple one-zone and
global viscous models which, though inappropriate for detailed matching to
observed light curves, clearly indicate that turbulent disk fluctuations induce
outbursts when the system is sufficiently close to the global bifurcation
point. While the models easily produce the observed "outburst dip" pairs
exhibited by Z Cam and Nova-like variables, they struggle to generate long
trains of outbursts. We conclude that mass transfer variability is the dominant
physical process determining the overall Z Cam standstill outburst pattern, but
that viscous stochasticity provides an additional ingredient explaining some of
the secondary features observed.

[18] oai:arXiv.org:1703.00211 [pdf] - 1540203

MRI turbulence and thermal instability in accretion disks

Comments: 16 pages, 16 figures, MNRAS accepted

Submitted: **2017-03-01**

A long-standing puzzle in the study of black-hole accretion concerns the
presence or not of thermal instability. Classical theory predicts the
encircling accretion disk is unstable, as do some self-consistent MHD
simulations of the flow. Yet observations of strongly accreting sources
generally fail to exhibit cyclic or unstable dynamics on the expected
timescales. This paper checks whether turbulent fluctuations impede thermal
instability. It also asks if it makes sense to conduct linear stability
analyses on a turbulent background. These issues are explored with a set of MRI
simulations in thermally unstable local boxes in combination with stochastic
equations that approximate the disk energetics. These models show that the
disk's thermal behaviour deviates significantly from laminar theory, though
ultimately a thermal runaway does occur. We find that the disk temperature
evolves as a biased random walk, rather than increasing exponentially, and thus
generates a broad spread of outcomes, with instability often delayed for
several thermal times. We construct a statistical theory that describes some of
this behaviour, emphasising the importance of the `escape time' and its
associated probability distribution. In conclusion, turbulent fluctuations on
there own cannot stabilise a disk, but they can weaken and delay thermal
instability.

[19] oai:arXiv.org:1701.04312 [pdf] - 1767314

Planetary rings and other astrophysical disks

Comments: To appear in: Tiscareno & Murray Eds (2017), Planetary Ring Systems,
Cambridge University Press

Submitted: **2017-01-16**

This chapter explores the physics shared by planetary rings and the various
disks that populate the Universe. It begins with an observational overview,
ranging from protoplanetary disks to spiral galaxies, and then compares and
contrasts these astrophysical disks with the rings of the Solar System.
Emphasis is placed on fundamental physics and dynamics, and how research into
the two classes of object connects. Topics covered include disk formation,
accretion, collisional processes, waves, instabilities, and satellite-disk
interactions.

[20] oai:arXiv.org:1609.07930 [pdf] - 1513712

On dust-gas gravitational instabilities in protoplanetary discs

Comments: MNRAS accepted. 14 pages, 6 figures

Submitted: **2016-09-26**

In protoplanetary disks the aerodynamical friction between particles and gas
induces a variety of instabilities that facilitate planet formation. Of these
we examine the so-called `secular gravitational instability' (SGI) in the
two-fluid approximation, deriving analytical expressions for its stability
criteria and growth rates. Concurrently, we present a physical explanation of
the instability that shows how it manifests upon an intermediate range of
lengthscales exhibiting geostrophic balance in the gas component. In contrast
to a single-fluid treatment, the SGI is quenched within a critical disk radius,
as large as 10 AU and 30 AU for cm and mm sized particles respectively,
although establishing robust estimates is hampered by uncertainties in the
parameters (especially the strength of turbulence) and deficiencies in the
razor-thin disk model we employ. It is unlikely, however, that the SGI is
relevant for well-coupled dust. We conclude by applying these results to the
question of planetesimal formation and the provenance of large-scale dust
rings.

[21] oai:arXiv.org:1608.08557 [pdf] - 1475527

Magnetorotationally driven wind cycles in local disc models

Comments: 19 pages, 16 figures, accepted in MNRAS (August, 30th 2016)

Submitted: **2016-08-30**

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

[22] oai:arXiv.org:1606.04509 [pdf] - 1444082

The magnetorotational instability in debris-disc gas

Comments: 8 pages, 5 figures, accepted for publication in MNRAS

Submitted: **2016-06-14**

Debris discs are commonly swathed in gas which can be observed in UV, in fine
structure lines in FIR, and in resolved maps of CO emission. Carbon and oxygen
are overabundant in such gas, but it is severely depleted in hydrogen. As a
consequence, its ionisation fraction is remarkably high, suggesting
magnetohydrodynamic (MHD) processes may be important. In particular, the gas
may be subject to the magnetorotational instability (MRI), and indeed recent
modelling of $\beta$ Pictoris requires an anomalous viscosity to explain the
gas's observed radial structure. In this paper we explore the possibility that
the MRI is active in debris-disc gas and responsible for the observed mass
transport. We find that non-ideal MHD and dust-gas interactions play a
subdominant role, and that linear instability is viable at certain radii.
However, owing to low gas densities, the outer parts of the disc could be
stabilised by a weak ambient magnetic field, though it is difficult to
constrain such a field. Even if the MRI is stabilised by too strong a field, a
magnetocentrifugal wind may be launched in its place and this could lead to
equivalent (non-turbulent) transport. Numerical simulations of the vertically
stratified MRI in conditions appropriate to the debris disc gas should be able
to determine the nature of the characteristic behaviour at different radii, and
decide on the importance of the MRI (and MHD more generally) on the evolution
of these discs.

[23] oai:arXiv.org:1606.03012 [pdf] - 1483375

On the survival of zombie vortices in protoplanetary discs

Comments: 6 pages, 8 figures, 1 table. Submitted to MNRAS Letters

Submitted: **2016-06-09**

Recently it has been proposed that the zombie vortex instability (ZVI) could
precipitate hydrodynamical activity and angular momentum transport in
unmagnetised regions of protoplanetary discs, also known as "dead zones". In
this letter we scrutinise, with high resolution 3D spectral simulations, the
onset and survival of this instability in the presence of viscous and thermal
physics. First, we find that the ZVI is strongly dependent on the nature of the
viscous operator. Although the ZVI is easily obtained with hyper-diffusion, it
is difficult to sustain with physical (second order) diffusion operators up to
Reynolds numbers as high as 10^7. This sensitivity is probably due to the ZVI's
reliance on critical layers, whose characteristic lengthscale, structure, and
dynamics are controlled by viscous diffusion. Second, we observe that the ZVI
is sensitive to radiative processes, and indeed only operates when the Peclet
number is greater than a critical value ~10^4, or when the cooling time is
longer than ~10 Omega^{-1}. As a consequence, the ZVI struggles to appear at
R>0.3 AU in standard 0.01 M_sun T Tauri disc models, though younger more
massive disks provide a more hospitable environment. Together these results
question the prevalence of the ZVI in protoplanetary discs.

[24] oai:arXiv.org:1605.01873 [pdf] - 1427864

Gravitoturbulence in magnetised protostellar discs

Comments: MNRAS submitted, revised after reviewer comments

Submitted: **2016-05-06**

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

[25] oai:arXiv.org:1510.06247 [pdf] - 1319697

On the convective overstability in protoplanetary discs

Comments: MNRAS accepted

Submitted: **2015-10-21**

This paper explores the driving of low-level hydrodynamical activity in
protoplanetary-disc dead zones. A small adverse radial entropy gradient,
ordinarily stabilised by rotation, excites oscillatory convection (`convective
overstability') when thermal diffusion, or cooling, is neither too strong nor
too weak. I revisit the linear theory of the instability, discuss its
prevalence in protoplanetary discs, and show that unstable modes are exact
nonlinear solutions in the local Boussinesq limit. Overstable modes cannot grow
indefinitely, however, as they are subject to a secondary parametric
instability that limits their amplitudes to relatively low levels. If parasites
set the saturation level of the ensuing turbulence then the convective
overstability is probably too weak to drive significant angular momentum
transport or to generate vortices. But I also discuss an alternative, and far
more vigorous, saturation route that generates radial `layers' or `zonal flows'
(witnessed also in semiconvection). Numerical simulations are required to
determine which outcome is favoured in realistic discs, and consequently how
important the instability is for disc dynamics.

[26] oai:arXiv.org:1510.01214 [pdf] - 1308344

The stress-pressure relationship in simulations of MRI-induced
turbulence

Comments: Accepted for publication in MNRAS

Submitted: **2015-10-05**

We determine how MRI-turbulent stresses depend on gas pressure via a suite of
unstratified shearing box simulations. Earlier numerical work reported only a
very weak dependence at best, results that call into question the canonical
alpha-disk model and the thermal stability results that follow from it. Our
simulations, in contrast, exhibit a stronger relationship, and show that
previous work was box-size limited: turbulent `eddies' were artificially
restricted by the numerical domain rather than by the scale height.
Zero-net-flux runs without physical diffusion coefficients yield a stress
proportional to $P^{0.5}$, where P is pressure. The stresses are also
proportional to the grid length and hence remain numerically unconverged. The
same runs with physical diffusivities, however, give a result closer to an
alpha-disk: the stress is proportional to $P^{0.9}$. Net-flux simulations
without explicit diffusion exhibit stresses proportional to $P^{0.5}$, but
stronger imposed fields weaken this correlation. In summary, compressibility is
important for the saturation of the MRI, but the exact stress-pressure
relationship is difficult to ascertain in local simulations because of
numerical convergence issues and the influence of any imposed flux. As a
consequence, the interpretation of thermal stability behaviour in local
simulations is a problematic enterprise.

[27] oai:arXiv.org:1508.03012 [pdf] - 1292267

Local and global aspects of the linear MRI in accretion disks

Comments: 13 pages, 9 figures, MNRAS accepted

Submitted: **2015-08-12**

We revisit the linear MRI in a cylindrical model of an accretion disk and
uncover a number of attractive results overlooked in previous treatments. In
particular, we elucidate the connection between local axisymmetric modes and
global modes, and show that a local channel flow corresponds to the evanescent
part of a global mode. In addition, we find that the global problem reproduces
the local dispersion relation without approximation, a result that helps
explain the success the local analysis enjoys in predicting global growth
rates. MRI channel flows are nonlinear solutions to the governing equations in
the local shearing box. However, only a small subset of MRI modes share the
same property in global disk models, providing further evidence that the
prominence of channels in local boxes is artificial. Finally, we verify our
results via direct numerical simulations with the Godunov code RAMSES.

[28] oai:arXiv.org:1503.06953 [pdf] - 1232099

On the vertical-shear instability in astrophysical discs

Comments: 18 pages, 15 figures, accepted for publication in MNRAS (2015 March
22)

Submitted: **2015-03-24**

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.

[29] oai:arXiv.org:1411.3236 [pdf] - 1223135

Vortex cycles at the inner edges of dead zones in protoplanetary disks

Comments: 13 pages, 15 figures

Submitted: **2014-11-12**

In protoplanetary disks, the inner boundary between the turbulent and laminar
regions is a promising site for planet formation because solids may become
trapped at the interface itself or in vortices generated by the Rossby wave
instability. The disk thermodynamics and the turbulent dynamics at that
location are entwined because of the importance of turbulent dissipation on
thermal ionization and, conversely, of thermal ionisation on the turbulence.
However, most previous work has neglected this dynamical coupling and have thus
missed a key element of the physics in this region. In this paper, we aim to
determine how the the interplay between ionization and turbulence impacts on
the formation and evolution of vortices at the interface between the active and
the dead zones. Using the Godunov code RAMSES, we have performed a 3D
magnetohydrodynamic global numerical simulation of a cylindrical model of an
MRI--turbulent protoplanetary disk, including thermodynamical effects as well
as a temperature-dependant resistivity. The comparison with an analogous 2D
viscous simulation has been extensively used to help identify the relevant
physical processes and the disk's long-term evolution. We find that a vortex
formed at the interface, due to Rossby wave instability, migrates inward and
penetrates the active zone where it is destroyed by turbulent motions.
Subsequently, a new vortex emerges a few tens of orbits later at the interface,
and the new vortex migrates inward too. The sequence repeats itself, resulting
in cycles of vortex formation, migration, and disruption. This behavior is
successfully reproduced using two different codes. In this paper, we
characterize this vortex life cycle and discuss its implications for planet
formation at the dead/active interface. Our simulations highlight the
importance of thermodynamical processes for the vortex evolution at the dead
zone inner edge.

[30] oai:arXiv.org:1404.4247 [pdf] - 1292112

The ballistic transport instability in Saturn's rings III: numerical
simulations

Comments: 11 pages, 9 figures, MNRAS accepted

Submitted: **2014-04-16**

Saturn's inner B-ring and its C-ring support wavetrains of contrasting
amplitudes but with similar length scales, 100-1000 km. In addition, the inner
B-ring is punctuated by two intriguing `flat' regions between radii 93,000 km
and 98,000 km in which the waves die out, whereas the C-ring waves coexist with
a forest of plateaus, narrow ringlets, and gaps. In both regions the waves are
probably generated by a large-scale linear instability whose origin lies in the
meteoritic bombardment of the rings: the ballistic transport instability. In
this paper, the third in a series, we numerically simulate the long-term
nonlinear evolution of this instability in a convenient local model. Our C-ring
simulations confirm that the unstable system forms low-amplitude wavetrains
possessing a preferred band of wavelengths. B-ring simulations, on the other
hand, exhibit localised nonlinear wave `packets' separated by linearly stable
flat zones. Wave packets travel slowly while spreading in time, a result that
suggests the observed flat regions in Saturn's B-ring are shrinking. Finally,
we present exploratory runs of the inner B-ring edge which reproduce earlier
numerical results: ballistic transport can maintain the sharpness of a
spreading edge while building a `ramp' structure at its base. Moreover, the
ballistic transport instability can afflict the ramp region, but only in
low-viscosity runs.

[31] oai:arXiv.org:1404.4245 [pdf] - 1292111

The ballistic transport instability in Saturn's rings II: nonlinear wave
dynamics

Comments: 14 pages, 13 figures, MNRAS accepted

Submitted: **2014-04-16**

The ejecta discharged by impacting meteorites can redistribute a planetary
ring's mass and angular momentum. This `ballistic transport' of ring properties
instigates a linear instability that could generate the 100--1000-km
undulations observed in Saturn's inner B-ring and in its C-ring. We present
semi-analytic results demonstrating how the instability sustains steadily
travelling nonlinear wavetrains. At low optical depths, the instability
produces approximately sinusoidal waves of low amplitude, which we identify
with those observed between radii 77,000 and 86,000 km in the C-ring. On the
other hand, optical depths of 1 or more exhibit hysteresis, whereby the ring
falls into multiple stable states: the homogeneous background equilibrium or
large-amplitude wave states. Possibly the `flat zones' and `wave zones' between
radii 93,000 and 98,000 km in the B-ring correspond to the stable homogeneous
and wave states, respectively. In addition, we test the linear stability of the
wavetrains and show that only a small subset are stable. In particular, stable
solutions all possess wavelengths greater than the lengthscale of fastest
linear growth. We supplement our calculations with a weakly nonlinear analysis
that suggests the C-ring reproduces some of the dynamics of the complex
Ginzburg--Landau equation. In the third paper in the series, these results will
be tested and extended with numerical simulations.

[32] oai:arXiv.org:1402.6525 [pdf] - 1203584

Thermodynamics of the dead-zone inner edge in protoplanetary disks

Comments: 16 pages, 16 figures. Accepted in Astronomy and Astrophysics

Submitted: **2014-02-26**

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

[33] oai:arXiv.org:1303.0264 [pdf] - 1164954

Hydrodynamic instability in warped astrophysical discs

Comments: 17 pages, 7 figures, revised version, to be published in MNRAS

Submitted: **2013-03-01**, last modified: **2013-07-01**

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

[34] oai:arXiv.org:1303.0263 [pdf] - 1164953

Local and global dynamics of warped astrophysical discs

Comments: 17 pages, 10 figures, revised version, to be published in MNRAS

Submitted: **2013-03-01**, last modified: **2013-07-01**

Astrophysical discs are warped whenever a misalignment is present in the
system, or when a flat disc is made unstable by external forces. The evolution
of the shape and mass distribution of a warped disc is driven not only by
external influences but also by an internal torque, which transports angular
momentum through the disc. This torque depends on internal flows driven by the
oscillating pressure gradient associated with the warp, and on physical
processes operating on smaller scales, which may include instability and
turbulence. We introduce a local model for the detailed study of warped discs.
Starting from the shearing sheet of Goldreich & Lynden-Bell, we impose the
oscillating geometry of the orbital plane by means of a coordinate
transformation. This warped shearing sheet (or box) is suitable for analytical
and computational treatments of fluid dynamics, magnetohydrodynamics, etc., and
it can be used to compute the internal torque that drives the large-scale
evolution of the disc. The simplest hydrodynamic states in the local model are
horizontally uniform laminar flows that oscillate at the orbital frequency.
These correspond to the nonlinear solutions for warped discs found in previous
work by Ogilvie, and we present an alternative derivation and generalization of
that theory. In a companion paper we show that these laminar flows are often
linearly unstable, especially if the disc is nearly Keplerian and of low
viscosity. The local model can be used in future work to determine the
nonlinear outcome of the hydrodynamic instability of warped discs, and its
interaction with others such as the magnetorotational instability.

[35] oai:arXiv.org:1210.6664 [pdf] - 1152426

Local outflows from turbulent accretion disks

Comments: 15 pages, 14 figures, changes according to referee report. Accepted
in Astronomy and Astrophysics

Submitted: **2012-10-24**, last modified: **2013-02-14**

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

[36] oai:arXiv.org:1210.3358 [pdf] - 1152060

Large-scale N-body simulations of the viscous overstability in Saturn's
rings

Comments: 15 pages, 12 figures, accepted for publication by MNRAS

Submitted: **2012-10-11**, last modified: **2013-01-24**

We present results from large-scale particle simulations of the viscous
overstability in Saturn's rings. The overstability generates a variety of
structure on scales covering a few hundred metres to several kilometres,
including axisymmetric wavetrains and larger-scale modulations. Such patterns
have been observed in Saturn's rings by the Cassini spacecraft. Our simulations
model the collisional evolution of particles in a co-rotating patch of the
disk. These are the largest N-body simulations of the viscous overstability yet
performed. The radial box size is five orders of magnitude larger than a
typical particle radius, and so describes a 20-50 km radial portion of the
rings. Its evolution is tracked for more than 10,000 orbits. In agreement with
hydrodynamics, our N-body simulations reveal that the viscous overstability
exhibits a rich set of dynamics characterised by nonlinear travelling waves
with wavelengths of a few hundred meters. In addition, wave defects, such as
sources and shocks, punctuate this bed of waves and break them up into
large-scale divisions of radial width ~5 km. We find that the wavelength of the
travelling waves is positively correlated with the mean optical depth. In order
to assess the role of the numerical boundary conditions and also background
ring structure, we include simulations of broad spreading rings and simulations
with a gradient in the background surface density. Overall, our numerical
results and approach provide a tool with which to interpret Cassini occultation
observations of microstructure in Saturn's rings. We present an example of such
a synthetic occultation observation and discuss what features to expect. We
also make the entire source code freely available.

[37] oai:arXiv.org:1209.2943 [pdf] - 1151389

The ballistic transport instability in Saturn's rings I: formalism and
linear theory

Comments: 14 pages, 7 figures, MNRAS accepted

Submitted: **2012-09-13**

Planetary rings sustain a continual bombardment of hypervelocity meteoroids
that erode the surfaces of ring particles on time scales of 10^5 - 10^7 years.
The debris ejected from such impacts re-accretes on to the ring, though often
at a slightly different orbital radius from the point of emission. This
`ballistic transport' leads to a rearrangement of the disk's mass and angular
momentum, and gives rise to a linear instability that generates structure on
relatively large scales. It is likely that the 100-km undulations in Saturn's
inner B-ring and the plateaus and 1000-km waves in Saturn's C-ring are
connected to the nonlinear saturation of the instability. In this paper the
physical problem is reformulated so as to apply to a local patch of disk (the
shearing sheet). This new streamlined model helps facilitate our physical
understanding of the instability, and also makes more tractable the analysis of
its nonlinear dynamics. We concentrate on the linear theory in this paper,
showing that the instability is restricted to a preferred range of intermediate
wavenumbers and optical depths. We subsequently apply these general results to
the inner B-ring and the C-ring and find that in both regions the ballistic
transport instability should be near marginality, a fact that may have
important consequences for its prevalence and nonlinear development. Owing to
damping via self-gravity wakes, the instability should not be present in the
A-ring. A following paper will explore the instability's nonlinear saturation
and how it connects to the observed large-scale structure.

[38] oai:arXiv.org:1207.4727 [pdf] - 1124961

Hysteresis and thermal limit cycles in MRI simulations of accretion
discs

Comments: 15 pages, 8 figures, 4 tables; MNRAS, accepted; v2 minor corrections
and clarifications

Submitted: **2012-07-19**, last modified: **2012-08-21**

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

[39] oai:arXiv.org:1203.6572 [pdf] - 1117641

The dynamics of inner dead-zone boundaries in protoplanetary disks

Comments: 16 pages, 10 figures; MNRAS accepted; v3 final corrections

Submitted: **2012-03-29**, last modified: **2012-08-21**

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.

[40] oai:arXiv.org:1206.2677 [pdf] - 1124078

Tidal disruption of satellites and formation of narrow rings

Comments: Accepted for publication in MNRAS. Some figures have been compressed
to fit into astro-ph size guidelines. Please contact authors if full
resolution images are required

Submitted: **2012-06-12**

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

[41] oai:arXiv.org:1202.3440 [pdf] - 1116676

The HBI in a quasi-global model of the intracluster medium

Comments: MNRAS, in press; minor modifications from v1

Submitted: **2012-02-15**, last modified: **2012-04-27**

In this paper we investigate how convective instabilities influence heat
conduction in the intracluster medium (ICM) of cool-core galaxy clusters. The
ICM is a high-beta, weakly collisional plasma in which the transport of
momentum and heat is aligned with the magnetic field. The anisotropy of heat
conduction, in particular, gives rise to instabilities that can access energy
stored in a temperature gradient of either sign. We focus on the heat-flux
buoyancy-driven instability (HBI), which feeds on the outwardly increasing
temperature profile of cluster cool cores. Our aim is to elucidate how the
global structure of a cluster impacts on the growth and morphology of the
linear HBI modes when in the presence of Braginskii viscosity, and ultimately
on the ability of the HBI to thermally insulate cores. We employ an idealised
quasi-global model, the plane-parallel atmosphere, which captures the essential
physics -- e.g. the global radial profile of the cluster -- while letting the
problem remain analytically tractable. Our main result is that the dominant HBI
modes are localised to the the innermost (~<20%) regions of cool cores. It is
then probable that, in the nonlinear regime, appreciable field-line insulation
will be similarly localised. Thus, while radio-mode feedback appears necessary
in the central few tens of kpc, heat conduction may be capable of offsetting
radiative losses throughout most of a cool core over a significant fraction of
the Hubble time. Finally, our linear solutions provide a convenient numerical
test for the nonlinear codes that tackle the saturation of such convective
instabilities in the presence of anisotropic transport.

[42] oai:arXiv.org:1203.4192 [pdf] - 1117392

The gravitational instability of a stream of co-orbital particles

Comments: 11 pages, 7 figures, 1 table. MNRAS, accepted

Submitted: **2012-03-19**

We describe the dynamics of a stream of equally spaced macroscopic particles
in orbit around a central body (e.g. a planet or star). A co-orbital
configuration of small bodies may be subject to gravitational instability,
which takes the system to a spreading, disordered and collisional state. We
detail the linear instability's mathematical and physical features using the
shearing sheet model and subsequently track its nonlinear evolution with local
N-body simulations. This model provides a convenient tool with which to
understand the gravitational and collisional dynamics of narrow belts, such as
Saturn's F-ring and the streams of material wrenched from tidally disrupted
bodies. In particular, we study the tendency of these systems to form
long-lived particle aggregates. Finally, we uncover an unexpected connection
between the linear dynamics of the gravitational instability and the
magnetorotational instability.

[43] oai:arXiv.org:1111.3809 [pdf] - 1091700

Global model of differential rotation in the Sun

Comments: 26 pages, 7 figures, accepted to MNRAS

Submitted: **2011-11-16**

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.

[44] oai:arXiv.org:1104.5396 [pdf] - 1076285

On linear dust-gas streaming instabilities in protoplanetary discs

Comments: 9 pages, 1 figure. Accepted to MNRAS

Submitted: **2011-04-28**

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.

[45] oai:arXiv.org:1006.2308 [pdf] - 1033059

Spherical single-roll dynamos at large magnetic Reynolds numbers

Comments: 36 pages, 8 figures. V2: incorrect labels in Fig3 corrected. The
article appears in Physics of Fluids, 22, 066601, and may be found at
http://pof.aip.org/phfle6/v22/i6/p066601_s1 . (Copyright 2010 American
Institute of Physics. This article may be downloaded for personal use only.
Any other use requires prior permission of the author and the American
Institute of Physics)

Submitted: **2010-06-11**, last modified: **2010-06-30**

This paper concerns kinematic helical dynamos in a spherical fluid body
surrounded by an insulator. In particular, we examine their behaviour in the
regime of large magnetic Reynolds number $\Rm$, for which dynamo action is
usually concentrated upon a simple resonant stream-surface. The dynamo
eigensolutions are computed numerically for two representative single-roll
flows using a compact spherical harmonic decomposition and fourth-order
finite-differences in radius. These solutions are then compared with the growth
rates and eigenfunctions of the Gilbert and Ponty (2000) large $\Rm$ asymptotic
theory. We find good agreement between the growth rates when $\Rm>10^4$, and
between the eigenfunctions when $\Rm>10^5$.

[46] oai:arXiv.org:1006.3419 [pdf] - 1033152

Hydrodynamical simulations of viscous overstability in Saturn's rings

Comments: 30 pages, 13 figures. Icarus, accepted

Submitted: **2010-06-17**

We perform axisymmetric hydrodynamical simulations that describe the
nonlinear outcome of the viscous overstability in dense planetary rings. These
simulations are particularly relevant for Cassini observations of fine-scale
structure in Saturn's A and B-ring, which take the form of periodic
microstructure on the 0.1 km scale, and irregular larger-scale structure on
1-10 km. Nonlinear wavetrains dominate all the simulations, and we associate
them with the observed periodic microstructure. The waves can undergo small
chaotic fluctuations in their phase and amplitude, and may be punctuated by
more formidable `wave defects' distributed on longer scales. It is unclear,
however, whether the defects are connected to the irregular larger-scale
variations observed by Cassini. The long-term behaviour of the simulations is
dominated by the imposed boundary conditions, and more generally by the
limitations of the local model we use: the shearing box. When periodic boundary
conditions are imposed, the system eventually settles on a uniform travelling
wave of a predictable wavelength, while reflecting boundaries, and boundaries
with buffer zones, maintain a disordered state. The simulations omit
self-gravity, though we examine its influence in future work.

[47] oai:arXiv.org:1005.5100 [pdf] - 1032796

The effect of the tachocline on differential rotation in the Sun

Comments: 20 Pages, 4 figures, to appear in MNRAS

Submitted: **2010-05-27**

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.

[48] oai:arXiv.org:1004.0109 [pdf] - 1026063

MRI channel flows in vertically-stratified models of accretion disks

Comments: 17 pages, 12 figures. MNRAS, accepted.

Submitted: **2010-04-01**

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

[49] oai:arXiv.org:1002.3750 [pdf] - 1025290

Resistive double-diffusive instability in the dead-zones of protostellar
disks

Comments: 16 pages, 5 figures. MNRAS Accepted. V2: cosmetic changes to bring in
line with MNRAS version

Submitted: **2010-02-19**, last modified: **2010-03-23**

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.

[50] oai:arXiv.org:0907.0641 [pdf] - 250747

Inertial waves near corotation in 3D hydrodynamical disks

Comments: 17 pages, 5 figures. MNRAS accepted. V2 - Section 5.2 moved to
appendix and errors removed

Submitted: **2009-07-03**, last modified: **2009-10-20**

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.

[51] oai:arXiv.org:0907.5075 [pdf] - 1003161

On Differential Rotation and Convection in the Sun

Comments: 16 pages, 2 figures. Accepted for publication in MNRAS

Submitted: **2009-07-29**

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.

[52] oai:arXiv.org:0904.3190 [pdf] - 1001979

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

Comments: MNRAS accepted, 9 pages, 11 figures

Submitted: **2009-04-21**

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.

[53] oai:arXiv.org:0904.0143 [pdf] - 1001832

The viscous overstability, nonlinear wavetrains, and finescale structure
in dense planetary rings

Comments: 52 pages, 8 figures. Accepted in Icarus

Submitted: **2009-04-01**

This paper addresses the fine-scale axisymmetric structure exhibited in
Saturn's A and B-rings. We aim to explain both the periodic microstructure on
150-220m, revealed by the Cassini UVIS and RSS instruments, and the irregular
variations in brightness on 1-10km, reported by the Cassini ISS. We propose
that the former structures correspond to the peaks and troughs of the nonlinear
wavetrains that form naturally in a viscously overstable disk. The latter
variations on longer scales may correspond to modulations and defects in the
wavetrains' amplitudes and wavelength. We explore these ideas using a simple
hydrodynamical model which captures the correct qualitative behaviour of a disk
of inelastically colliding particles, while also permitting us to make progress
with analytic and semi-analytic techniques. Specifically, we calculate a family
of travelling nonlinear density waves and determine their stability properties.
Detailed numerical simulations that confirm our basic results will appear in a
following paper.

[54] oai:arXiv.org:0901.4712 [pdf] - 250691

MRI channel flows and their parasites

Comments: 17 pages, 17 figures, accepted by MNRAS

Submitted: **2009-01-29**

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.

[55] oai:arXiv.org:0803.4123 [pdf] - 11294

Dense planetary rings and the viscous overstability

Comments: Accepted in Icarus; 58 pages; 15 figures

Submitted: **2008-03-28**

This paper examines the onset of the viscous overstability in dense
particulate rings. First, we formulate a dense gas kinetic theory that is
applicable to the Saturnian system. Our model is essentially that of Araki and
Tremaine (1986), which we show can be both simplified and generalised. Second,
we put this model to work computing the equilibrium properties of dense
planetary rings, which we subsequently compare with the results of N-body
simulations, namely those of Salo (1991). Finally, we present the linear
stability analyses of these equilibrium states, and derive criteria for the
onset of viscous overstability in the self-gravitating and non-self-gravitating
cases. These are framed in terms of particle size, orbital frequency, optical
depth, and the parameters of the collision law. Our results compare favourably
with the simulations of Salo et al. (2001). The accuracy and practicality of
the continuum model we develop encourages its general use in future
investigations of nonlinear phenomena.

[56] oai:arXiv.org:astro-ph/0609164 [pdf] - 84756

The linear stability of dilute particulate rings

Comments:

Submitted: **2006-09-06**

Irregular structure in planetary rings is often attributed to the intrinsic
instabilities of a homogeneous state undergoing Keplerian shear. Previously
these have been analysed with simple hydrodynamic models. We instead employ a
kinetic theory, in which we solve the linearised moment equations derived in
Shu and Stewart 1985 for a dilute ring. This facilitates an examination of
velocity anisotropy and non-Newtonian stress, and their effects on the viscous
and viscous/gravitational instabilities thought to occur in Saturn's rings.
Because we adopt a dilute gas model, the applicability of our results to the
actual dense rings of Saturn are significantly curtailled. Nevertheless this
study is a necessary preliminary before an attack on the difficult problem of
dense ring dynamics. We find the Shu and Stewart formalism admits analytic
stability criteria for the viscous overstability, viscous instability, and
thermal instability. These criteria are compared with those of a hydrodynamic
model incorporating the effective viscosity and cooling function computed from
the kinetic steady state. We find the two agree in the `hydrodynamic limit'
(i.e. many collisions per orbit) but disagree when collisions are less
frequent, when we expect the viscous stress to be increasingly non-Newtonian
and the velocity distribution increasingly anisotropic. In particular,
hydrodynamics predicts viscous overstability for a larger portion of parameter
space. We also numerically solve the linearised equations of the more accurate
Goldreich and Tremaine 1978 kinetic model and discover its linear stability to
be qualitatively the same as that of Shu and Stewart's. Thus the simple
collision operator adopted in the latter would appear to be an adequate
approximation for dilute rings, at least in the linear regime.

[57] oai:arXiv.org:astro-ph/0609158 [pdf] - 84750

Viscous overstability and eccentricity evolution in three-dimensional
gaseous discs

Comments: 11 pages, 8 figures

Submitted: **2006-09-06**

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