Normalized to: Lesur, G.

**50** article(s) in total. **94** co-authors, from **1** to **12** common article(s). Median position in authors list is **2,0**.

[1] oai:arXiv.org:2006.01194 [pdf] - 2105821

Ring formation and dust dynamics in wind-driven protoplanetary discs:
global simulations

Comments: 17 pages, accepted in A&A

Submitted: **2020-06-01**

Large-scale vertical magnetic fields are believed to play a key role in the
evolution of protoplanetary discs. Associated with non-ideal effects, such as
ambipolar diffusion, they are known to launch a wind that could drive accretion
in the outer part of the disc ($R> 1$ AU). They also potentially lead to
self-organisation of the disc into large-scale axisymmetric structures, similar
to the rings recently imaged by sub-millimetre or near-infrared instruments
(ALMA and SPHERE). The aim of this paper is to investigate the mechanism behind
the formation of these gaseous rings, but also to understand the dust dynamics
and its emission in discs threaded by a large-scale magnetic field. To this
end, we performed global magneto-hydrodynamics (MHD) axisymmetric simulations
with ambipolar diffusion using a modified version of the PLUTO code. We
explored different magnetisations with the midplane $\beta$ parameter ranging
from $10^5$ to $10^3$ and included dust grains -- treated in the fluid
approximation -- ranging from $100 \mu$m to 1 cm in size. We first show that
the gaseous rings (associated with zonal flows) are tightly linked to the
existence of MHD winds. Secondly, we find that millimetre-size dust is highly
sedimented, with a typical scale height of 1 AU at $R=100$ AU for $\beta=10^4$,
compatible with recent ALMA observations. We also show that these grains
concentrate into pressure maxima associated with zonal flows, leading to the
formation of dusty rings. Using the radiative transfer code MCFOST, we computed
the dust emission and make predictions on the ring-gap contrast and the
spectral index that one might observe with interferometers like ALMA.

[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:1909.12258 [pdf] - 1979782

Magnetically-driven jets and winds from weakly magnetized accretion
disks

Comments: Accepted 23/09/2019

Submitted: **2019-09-26**, last modified: **2019-10-14**

Semi-analytical models of disk outflows have successfully described
magnetically-driven, self-confined super-Alfv\'enic jets from near Keplerian
accretion disks. These Jet Emitting Disks are possible for high levels of disk
magnetization $\mu$ defined as $\mu=2/\beta$ where beta is the usual plasma
parameter. In near-equipartition JEDs, accretion is supersonic and jets carry
away most of the disk angular momentum. However, these solutions prove
difficult to compare with cutting edge numerical simulations, for the reason
that numerical simulations show wind-like outflows but in the domain of small
magnetization. In this work, we present for the first time self-similar
isothermal solutions for accretion-ejection structures at small magnetization
levels. We elucidate the role of MRI-like structures in the acceleration
processes that drive this new class of solutions. The disk magnetization $\mu$
is the main control parameter: massive outflows driven by the pressure of the
toroidal magnetic field are obtained up to $\mu \sim 10^{-2}$, while more
tenuous centrifugally-driven outflows are obtained at larger $\mu$ values. The
generalized parameter space and the astrophysical consequences are discussed.
We believe that these new solutions could be a stepping stone in understanding
the way astrophysical disks drive either winds or jets. Defining jets as
self-confined outflows and winds as uncollimated outflows, we propose a simple
analytical criterion based on the initial energy content of the outflow, to
discriminate jets from winds. We show that jet solution are achieved at all
magnetization level, while winds could be obtained only in weakly magnetized
disks that feature heating.

[4] oai:arXiv.org:1812.02076 [pdf] - 1912713

Magnetic wind-driven accretion in dwarf novae

Comments: 7 pages, 6 figures, accepted for publication in A&A

Submitted: **2018-12-05**, last modified: **2019-05-14**

Dwarf novae (DNe) and X-ray binaries exhibit outbursts thought to be due to a
thermal-viscous instability in the accretion disk. The disk instability model
(DIM) assumes that accretion is driven by turbulent transport, customarily
attributed to the magneto-rotational instability (MRI). Recent results point
out that MRI turbulence alone fails to reproduce the light curves of DNe. We
aim to study the impact of wind-driven accretion on the light curves of DNe.
Local and global simulations show that magneto-hydrodynamic winds are present
when a magnetic field threads the disk, even for relatively high ratios of
thermal pressure to magnetic pressure ($\beta \approx 10^{5}$). These winds are
very efficient in removing angular momentum but do not heat the disk; they do
not behave as MRI-driven turbulence. We add wind-driven transport in the
angular momentum equation of the DIM, assuming a fixed magnetic configuration:
dipolar or constant with radius. We use prescriptions for the wind torque and
the turbulent torque derived from shearing box simulations. The wind torque
enhances the accretion of matter, resulting in light curves that look like DNe
outbursts when assuming a dipolar field with a moment
$\mu\approx10^{30}\,\mathrm{G\,cm^{3}}$. In the region where the wind
dominates, the disk is cold, optically thin and the accretion speed is sonic.
This acts as if the inner disk was truncated, leading to higher quiescent X-ray
luminosities from the white dwarf boundary layer than expected with the
standard DIM. The disk is stabilized if the wind-dominated region is large
enough, potentially leading to `dark' disks emitting little radiation.
Wind-driven accretion can play a key role in shaping the light curves of DNe
and X-ray binaries. Future studies will need to include the time evolution of
the magnetic field threading the disk to fully assess its impact on the
dynamics of the accretion flow.

[5] oai:arXiv.org:1904.07910 [pdf] - 1890477

Spontaneous ring formation in wind-emitting accretion discs

Comments: 17 pages, accepted in A&A

Submitted: **2019-04-16**

Rings and gaps have been observed in a wide range of protoplanetary discs,
from young systems like HLTau to older discs like TW Hydra. Recent disc
simulations have shown that magnetohydrodynamic (MHD) turbulence (in the ideal
or non-ideal regime) can lead to the formation of rings and be an alternative
to the embedded planets scenario. In this paper, we investigate how these ring
form in this context and seek a generic formation process, taking into account
the various dissipative regimes and magnetizations probed by the past
simulations. We identify the existence of a linear and secular instability,
driven by MHD winds, and giving birth to rings of gas having a width larger
than the disc scale height. We show that the linear theory is able to make
reliable predictions regarding the growth rates, ring/gap contrast and spacing,
by comparing these predictions to a series of 2D (axisymmetric) and 3D MHD
numerical simulations. In addition, we demonstrate that these rings can act as
dust traps provided that the disc is sufficiently magnetised, with plasma beta
lower than $10^4$. Given its robustness, the process identified in this paper
could have important implications, not only for protoplanetary discs but also
for a wide range of accreting systems threaded by large-scale magnetic fields.

[6] oai:arXiv.org:1811.08663 [pdf] - 2025343

Radiation thermo-chemical models of protoplanetary disks. Grain and
polycyclic aromatic hydrocarbon charging

Comments: accepted to A&A. 29 pages

Submitted: **2018-11-21**

Context. Disks around pre-main-sequence stars evolve over time by turbulent
viscous spreading. The main contender to explain the strength of the turbulence
is the Magneto-Rotational-Instability (MRI) model, whose efficiency depends on
the disk ionization fraction. Aims. We aim at computing self-consistently the
chemistry including PAH charge chemistry, the grain charging and an estimate of
an effective value of the turbulence alpha parameter in order to find
observational signatures of disk turbulence. Methods. We introduced PAH and
grain charging physics and their interplay with other gas-phase reactions in
the physico-chemical code ProDiMo. Non-ideal magnetohydrodynamics effects such
as Ohmic and ambipolar diffusion are parametrized to derive an effective value
for the turbulent parameter alpha_eff . We explored the effects of turbulence
heating and line broadening on CO isotopologue sub-millimeter lines. Results.
The spatial distribution of alpha_eff depends on various unconstrained disk
parameters such as the magnetic parameter beta_mag or the cosmic ray density
distribution inside the protoplanetary disks. The inner disk midplane shows the
presence of the so-called dead-zone where the turbulence is quasi-inexistent.
The disk is heated mostly by thermal accommodation on dust grains in the
dead-zone, by viscous heating outside the dead-zone up to a few hundred
astronomical units, and by chemical heating in the outer disk. The CO
rotational lines probe the warm molecular disk layers where the turbulence is
at its maximum. However, the effect of turbulence on the CO line profiles is
minimal and difficult to distinguish from the thermal broadening. Conclusions.
Viscous heating of the gas in the disk midplane outside the dead-zone is
efficient. The determination of alpha from CO rotational line observations
alone is challenging.

[7] oai:arXiv.org:1809.09131 [pdf] - 1791059

Turbulent and wind-driven accretion in dwarf novae threaded by a large
scale magnetic field

Comments: 13 pages, 7 figures, accepted in A&A. Typo in Figure 3 in the first
version. Revised now

Submitted: **2018-09-24**, last modified: **2018-11-21**

Dwarf novae (DNe) are accreting white dwarfs that show eruptions due to a
thermal-viscous instability in the accretion disk. The outburst timescales
constrain $\alpha$, the ratio of the viscous stress to the thermal pressure,
and so the mechanism of angular momentum transport. The eruptive state has
$\alpha\approx0.1$ while the quiescent state has $\alpha\approx0.03$. Turbulent
transport due to the magneto-rotational instability (MRI) is generally
considered to be the source of angular momentum transport in DNe. Here, we
perform 3D local magnetohydrodynamic (MHD) shearing box simulations including
vertical stratification, radiative transfer and a net constant vertical
magnetic flux to investigate how transport changes between the outburst and
quiescent states of DNe. We find that a constant $B_z$ provides a higher
$\alpha$ in quiescence than in outburst, in opposition to what is expected.
Including resistivity quenches MRI turbulence in quiescence, suppressing
transport, unless the magnetic field is high enough, which again leads to
$\alpha\approx0.1$. A major difference between simulations with a net poloidal
flux and simulations without is that angular momentum transport in the former
is shared between turbulent and wind-driven transport. We find that wind-driven
transport dominates in quiescence even for low magnetic fields $\sim 1$ G. This
can have a major impact on observational signatures since wind-driven transport
does not heat the disk. Furthermore, wind transport cannot be reduced to an
$\alpha$ prescription. We provide fits for $\alpha$ and the wind torque with
$\beta$, ratio of thermal to magnetic pressure. We conclude that the evolution
of the thermal-viscous instability, and its consequences on the outburst cycles
of CVs, needs to be seriously revised to take into account that most of the
accretion energy may be carried away by a wind instead of being locally
dissipated.

[8] oai:arXiv.org:1805.00458 [pdf] - 1759289

Dust settling and rings in the outer regions of protoplanetary discs
subject to ambipolar diffusion

Comments: 22 pages, 21 figures, accepted in A&A

Submitted: **2018-05-01**

MHD turbulence plays a crucial role in the dust dynamics of protoplanetary
discs. It affects planet formation, vertical settling and is one possible
origin of the large scale axisymmetric structures, such as rings, recently
imaged by ALMA and SPHERE. Among the variety of MHD processes, the
magnetorotational instability (MRI) has raised particular interest since it
provides a source of turbulence and potentially organizes the flow into radial
structures. However, the weak ionization of discs prevents the MRI from being
excited beyond 1 AU. The strong sedimentation of millimetre dust measured in
T-Tauri discs is also in contradiction with predictions based on ideal MRI
turbulence. In this paper, we study the effects of non-ideal MHD and winds on
the dynamics and sedimentation of dust grains. We consider a weakly ionized
plasma subject to ambipolar diffusion characterizing the disc outer regions (>1
AU). For that, we perform numerical MHD simulations in the stratified shearing
box, using the PLUTO code. Our simulations show that the mm-cm dust is
contained vertically in a very thin layer, with typical heightscale ~0.4 AU at
30 AU, compatible with recent ALMA observations. Horizontally, the grains are
trapped within pressure maxima induced by ambipolar diffusion, leading to the
formation of dust rings. For micrometer grains, dust and gas scaleheights are
similar. In this regime, the settling cannot be explained by a simple 1D
diffusion theory but results from a large scale 2D circulation induced by both
MHD winds and zonal flows. Overall, our results show that non-ideal MHD effects
and their related winds play a major role in shaping the radial and vertical
distribution of dust in protoplanetary discs. Leading to substantial accretion
efficiency, non-ideal effects also a promising avenue to reconcile the low
turbulent activity measured in discs with their relatively high accretion
rates.

[9] oai:arXiv.org:1710.05872 [pdf] - 1618336

Impact of convection and resistivity on angular momentum transport in
dwarf novae

Comments: 14 pages, 13 figures, 1 table; accepted for publication in A&A

Submitted: **2017-10-16**, last modified: **2017-10-17**

The eruptive cycles of dwarf novae (DN) are thought to be due to a
thermal-viscous instability in the accretion disk surrounding the white dwarf
(WD). This model has long been known to imply a stress to pressure ratio \alpha
~0.1 in outburst compared to \alpha ~ 0.01 in quiescence. Such an enhancement
in $\alpha$ has recently been observed in simulations of turbulent transport
driven by the magneto-rotational instability (MRI) when convection is present,
without requiring a net magnetic flux. We independently recover this result by
carrying out PLUTO MHD simulations of vertically stratified, radiative,
shearing boxes with the thermodynamics and opacities appropriate to DN. The
results are robust against the choice of vertical boundary conditions. The
thermal equilibrium solutions found by the simulations trace the well-known
S-curve in the density-temperature plane. We confirm that the high values of
\alpha ~ 0.1 occur near the tip of the hot branch of the S-curve, where
convection is active. However, we also present thermally-stable simulations at
lower temperatures that have standard values of \alpha ~ 0.03 despite the
presence of vigorous convection. We find no simple relationship between \alpha
and the strength of the convection, as measured by the ratio of convective to
radiative flux. The cold branch is only very weakly ionized so, in the second
part of this work, we studied the impact of non-ideal MHD effects on transport.
We include resistivity in the simulations and find that the MRI-driven
transport is quenched (\alpha ~ 0) below the critical density at which the
magnetic Reynolds number R_m \leq 10^4. This is problematic as X-ray emission
observed in quiescent systems requires ongoing accretion onto the WD. We verify
that these X-rays cannot self-sustain MRI-driven turbulence by photo-ionizing
the disk and discuss possible solutions to the issue of accretion in
quiescence.

[10] oai:arXiv.org:1709.09173 [pdf] - 1588920

Instability in strongly magnetized accretion discs: A global perspective

Comments: 23 pages, 7 Figures and 1 Table. Accepted for publication in MNRAS

Submitted: **2017-09-26**

We examine the properties of strongly magnetized accretion discs in a global
framework, with particular focus on the evolution of magnetohydrodynamic
instabilities such as the magnetorotational instability (MRI). Work by Pessah
and Psaltis showed that MRI is stabilized beyond a critical toroidal field in
compressible, differentially rotating flows and, also, reported the appearance
of two new instabilities beyond this field. Their results stemmed from
considering geometric curvature effects due to the suprathermal background
toroidal field, which had been previously ignored in weak-field studies.
However, their calculations were performed under the local approximation, which
poses the danger of introducing spurious behavior due to the introduction of
global geometric terms in an otherwise local framework. In order to avoid this,
we perform a global eigenvalue analysis of the linearized MHD equations in
cylindrical geometry. We confirm that MRI indeed tends to be highly suppressed
when the background toroidal field attains the Pessah-Psaltis limit. We also
observe the appearance of two new instabilities that emerge in the presence of
highly suprathermal toroidal fields. These results were additionally verified
using numerical simulations in PLUTO. There are, however, certain differences
between the the local and global results, especially in the vertical wavenumber
occupancies of the various instabilities, which we discuss in detail. We also
study the global eigenfunctions of the most unstable modes in the suprathermal
regime, which are inaccessible in the local analysis. Overall, our findings
emphasize the necessity of a global treatment for accurately modeling strongly
magnetized accretion discs.

[11] oai:arXiv.org:1705.03319 [pdf] - 1583253

Angular momentum transport in accretion disks: a hydrodynamical
perspective

Comments: 23 pages, 4 figures, proceeding for the Astrofluid conference in
honor of Jean-Paul Zahn (Paris, June 2016), references added

Submitted: **2017-05-09**, last modified: **2017-06-16**

The radial transport of angular momentum in accretion disk is a fundamental
process in the universe. It governs the dynamical evolution of accretion disks
and has implications for various issues ranging from the formation of planets
to the growth of supermassive black holes. While the importance of magnetic
fields for this problem has long been demonstrated, the existence of a source
of transport solely hydrodynamical in nature has proven more difficult to
establish and to quantify. In recent years, a combination of results coming
from experiments, theoretical work and numerical simulations has dramatically
improved our understanding of hydrodynamically mediated angular momentum
transport in accretion disk. Here, based on these recent developments, we
review the hydrodynamical processes that might contribute to transporting
angular momentum radially in accretion disks and highlight the many questions
that are still to be answered.

[12] oai:arXiv.org:1702.02857 [pdf] - 1550486

Universal small-scale structure in turbulence driven by
magnetorotational instability

Comments: 9 pages, 11 figures, to appear in Monthly Notices of the Royal
Astronomical Society

Submitted: **2017-02-09**

The intermittent small-scale structure of turbulence governs energy
dissipation in many astrophysical plasmas and is often believed to have
universal properties for sufficiently large systems. In this work, we argue
that small-scale turbulence in accretion disks is universal in the sense that
it is insensitive to the magnetorotational instability (MRI) and background
shear, and therefore indistinguishable from standard homogeneous
magnetohydrodynamic (MHD) turbulence at small scales. We investigate the
intermittency of current density, vorticity, and energy dissipation in
numerical simulations of incompressible MHD turbulence driven by the MRI in a
shearing box. We find that the simulations exhibit a similar degree of
intermittency as in standard MHD turbulence. We perform a statistical analysis
of intermittent dissipative structures and find that energy dissipation is
concentrated in thin sheet-like structures that span a wide range of scales up
to the box size. We show that these structures exhibit strikingly similar
statistical properties to those in standard MHD turbulence. Additionally, the
structures are oriented in the toroidal direction with a characteristic tilt of
approximately 17.5 degrees, implying an effective guide field in that
direction.

[13] oai:arXiv.org:1612.00883 [pdf] - 1557774

Global simulations of protoplanetary disks with net magnetic flux: I.
Non-ideal MHD case

Comments: 23 pages, 33 figures, accepted for publication in A&A

Submitted: **2016-12-02**, last modified: **2017-02-02**

The planet-forming region of protoplanetary disks is cold, dense, and
therefore weakly ionized. For this reason, magnetohydrodynamic (MHD) turbulence
is thought to be mostly absent, and another mechanism has to be found to
explain gas accretion. It has been proposed that magnetized winds, launched
from the ionized disk surface, could drive accretion in the presence of a
large-scale magnetic field. The efficiency and the impact of these surface
winds on the disk structure is still highly uncertain. We present the first
global simulations of a weakly ionized disk that exhibits large-scale
magnetized winds. We also study the impact of self-organization, which was
previously demonstrated only in non-stratified models. We perform numerical
simulations of stratified disks with the PLUTO code. We compute the ionization
fraction dynamically, and account for all three non-ideal MHD effects: ohmic
and ambipolar diffusions, and the Hall drift. Simplified heating and cooling
due to non-thermal radiation is also taken into account in the disk atmosphere.
We find that disks can be accreting or not, depending on the configuration of
the large-scale magnetic field. Magnetothermal winds, driven both by magnetic
acceleration and heating of the atmosphere, are obtained in the accreting case.
In some cases, these winds are asymmetric, ejecting predominantly on one side
of the disk. The wind mass loss rate depends primarily on the average ratio of
magnetic to thermal pressure in the disk midplane. The non-accreting case is
characterized by a meridional circulation, with accretion layers at the disk
surface and decretion in the midplane. Finally, we observe self-organization,
resulting in axisymmetric rings of density and associated pressure "bumps". The
underlying mechanism and its impact on observable structures are discussed.

[14] oai:arXiv.org:1611.07376 [pdf] - 1546302

Spiral-driven accretion in protoplanetary discs - III tri-dimensional
simulations

Comments: accepted for publication in A&A

Submitted: **2016-11-22**

Understanding how accretion proceeds in proto-planetary discs and more
generally their dynamics is a crucial issue for explaining the conditions in
which planets form. The role that accretion of gas from the surrounding
molecular cloud onto the disc may have on its structure needs to be quantified.
We perform tri-dimensional simulations using the Cartesian AMR code RAMSES of
an accretion disc subject to infalling material. For the aspect ratio of $H/R
\simeq 0.15$ and disk mass $M_d \simeq 10^{-2}$ M$_\odot$ used in our study, we
find that for typical accretion rates on the order of a few 10$^{-7}$ M$_\odot$
yr$^{-1}$, values of the $\alpha$ parameter as high as a few 10$^{-3}$ are
inferred. The mass that is accreted in the inner part of the disc is typically
at least $50\%$ of the total mass that has been accreted onto the disc. Our
results suggest that external accretion of gas at moderate values, onto
circumstellar discs may trigger prominent spiral arms, reminiscent of recent
observations made with various instruments, and lead to significant transport
through the disc. If confirmed from observational studies, such accretion may
therefore influence disc evolution.

[15] oai:arXiv.org:1610.10089 [pdf] - 1532346

Shadows and spirals in the protoplanetary disk HD 100453

Comments: A&A, accepted

Submitted: **2016-10-31**, last modified: **2016-11-03**

Understanding the diversity of planets requires to study the morphology and
the physical conditions in the protoplanetary disks in which they form. We
observed and spatially resolved the disk around the ~10 Myr old protoplanetary
disk HD 100453 in polarized scattered light with SPHERE/VLT at optical and
near-infrared wavelengths, reaching an angular resolution of ~0.02", and an
inner working angle of ~0.09". We detect polarized scattered light up to ~0.42"
(~48 au) and detect a cavity, a rim with azimuthal brightness variations at an
inclination of 38 degrees, two shadows and two symmetric spiral arms. The
spiral arms originate near the location of the shadows, close to the semi major
axis. We detect a faint spiral-like feature in the SW that can be interpreted
as the scattering surface of the bottom side of the disk, if the disk is
tidally truncated by the M-dwarf companion currently seen at a projected
distance of ~119 au. We construct a radiative transfer model that accounts for
the main characteristics of the features with an inner and outer disk
misaligned by ~72 degrees. The azimuthal brightness variations along the rim
are well reproduced with the scattering phase function of the model. While
spirals can be triggered by the tidal interaction with the companion, the close
proximity of the spirals to the shadows suggests that the shadows could also
play a role. The change in stellar illumination along the rim, induces an
azimuthal variation of the scale height that can contribute to the brightness
variations. Dark regions in polarized images of transition disks are now
detected in a handful of disks and often interpreted as shadows due to a
misaligned inner disk. The origin of such a misalignment in HD 100453, and of
the spirals, is unclear, and might be due to a yet-undetected massive companion
inside the cavity, and on an inclined orbit.

[16] oai:arXiv.org:1607.02903 [pdf] - 1530995

Magnetorotational dynamo chimeras. The missing link to turbulent
accretion disk dynamo models?

Comments: 12 pages, 8 figures, submitted to A&A

Submitted: **2016-07-11**

In Keplerian accretion disks, turbulence and magnetic fields may be jointly
excited through a subcritical dynamo process involving the magnetorotational
instability (MRI). High-resolution simulations exhibit a tendency towards
statistical self-organization of MRI dynamo turbulence into large-scale cyclic
dynamics. Understanding the physical origin of these structures, and whether
they can be sustained and transport angular momentum efficiently in
astrophysical conditions, represents a significant theoretical challenge. The
discovery of simple periodic nonlinear MRI dynamo solutions has recently proven
useful in this respect, and has notably served to highlight the role of
turbulent magnetic diffusion in the seeming decay of the dynamics at low
magnetic Prandtl number Pm (magnetic diffusivity larger than viscosity), a
common regime in accretion disks. The connection between these simple
structures and the statistical organization reported in turbulent simulations
remained elusive, though. Here, we report the numerical discovery in moderate
aspect ratio Keplerian shearing boxes of new periodic, incompressible,
three-dimensional nonlinear MRI dynamo solutions with a larger dynamical
complexity reminiscent of such simulations. These "chimera" cycles are
characterized by multiple MRI-unstable dynamical stages, but their basic
physical principles of self-sustainment are nevertheless identical to those of
simpler cycles found in azimuthally elongated boxes. In particular, we find
that they are not sustained at low Pm either due to subcritical turbulent
magnetic diffusion. These solutions offer a new perspective into the transition
from laminar to turbulent instability-driven dynamos, and may prove useful to
devise improved statistical models of turbulent accretion disk dynamos.

[17] 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.

[18] oai:arXiv.org:1603.02475 [pdf] - 1392915

Self-organisation in protoplanetary disks: global, non-stratified
Hall-MHD simulations

Comments: 16 pages, 21 figures

Submitted: **2016-03-08**

Recent observations revealed organised structures in protoplanetary disks,
such as axisymmetric rings or horseshoe concen- trations evocative of
large-scale vortices. These structures are often interpreted as the result of
planet-disc interactions. However, these disks are also known to be unstable to
the magneto-rotational instability (MRI) which is believed to be one of the
dominant angular momentum transport mechanism in these objects. It is therefore
natural to ask if the MRI itself could produce these structures without
invoking planets. The nonlinear evolution of the MRI is strongly affected by
the low ionisation fraction in protoplanetary disks. The Hall effect in
particular, which is dominant in dense and weakly ionised parts of these
objects, has been shown to spontaneously drive self- organising flows in
shearing box simulations. Here, we investigate the behaviour of global
MRI-unstable disc models dominated by the Hall effect and characterise their
dynamics. We perform 3D unstratified Hall-MHD simulations of keplerian disks
for a broad range of Hall, ohmic and ambipolar Elsasser numbers. We confirm the
transition from a turbulent to an organised state as the intensity of the Hall
effect is increased. We observe the formation of zonal flows, their number
depending on the available magnetic flux and on the intensity of the Hall
effect. For intermediate Hall intensity, the flow self-organises into
long-lived magnetised vortices. Neither the addition of a toroidal field nor
ohmic or ambipolar diffusion drastically change this picture in the range of
parameters we have explored. The ability of these structures to trap dust
particles in this configuration is demonstrated. We conclude that Hall-MRI
driven organisation is a plausible scenario which could explain some of the
structures found in recent observations.

[19] oai:arXiv.org:1602.01721 [pdf] - 1400446

Spiral-driven accretion in protoplanetary discs - II Self-similar
solutions

Comments: accepted for publication in A&A

Submitted: **2016-02-04**

Accretion discs are ubiquitous in the universe and it is a crucial issue to
understand how angular momentum and mass are being radially transported in
these objects. Here, we study the role played by non-linear spiral patterns
within hydrodynamical and non self-gravitating accretion disc assuming that
external disturbances such as infall onto the disc may trigger them. To do so,
we computed self-similar solutions that describe discs in which a spiral wave
propagates. Such solutions present both shocks and critical sonic points that
we carefully analyze. For all allowed temperatures and for several spiral
shocks, we calculated the wave structure. In particular we inferred the angle
of the spiral patern, the stress it exerts on the disc as well as the
associated flux of mass and angular momentum as a function of temperature. We
quantified the rate of angular momentum transport by means of the dimensionless
$\alpha$ parameter. For the thickest disc we considered (corresponding to $h/r$
values of about 1/3), we found values of $\alpha$ as high as $0.1$, and scaling
with the temperature $T$ such that $\alpha \propto T^{3/2} \propto (h/r)^3$.
The spiral angle scales with the temperature as $\arctan(r/h)$. The existence
of these solutions suggests that perturbations occurring at disc outer
boundaries, such as for example perturbations due to infall motions, can
propagate deep inside the disc and therefore should not be ignored, even when
considering small radii.

[20] oai:arXiv.org:1512.03739 [pdf] - 1326086

On the Nature of Magnetic Turbulence in Rotating, Shearing Flows

Comments: 5 pages, 6 figures, 1 table

Submitted: **2015-12-11**

The local properties of turbulence driven by the magnetorotational
instability (MRI) in rotating, shearing flows are studied in the framework of a
shearing-box model. Based on numerical simulations, we propose that the
MRI-driven turbulence comprises two components: the large-scale shear-aligned
strong magnetic field and the small-scale fluctuations resembling
magnetohydrodynamic (MHD) turbulence. The energy spectrum of the large-scale
component is close to $k^{-2}$, whereas the spectrum of the small-scale
component agrees with the spectrum of strong MHD turbulence $k^{-3/2}$. While
the spectrum of the fluctuations is universal, the outer-scale characteristics
of the turbulence are not; they depend on the parameters of the system, such as
the net magnetic flux. However, there is remarkable universality among the
allowed turbulent states -- their intensity $v_0$ and their outer scale
$\lambda_0$ satisfy the balance condition $v_0/\lambda_0\sim \mathrm
d\Omega/\mathrm d\ln r$, where $\mathrm d\Omega/\mathrm d\ln r$ is the local
orbital shearing rate of the flow. Finally, we find no sustained dynamo action
in the $\mathrm{Pm}=1$ zero net-flux case for Reynolds numbers as high as
$45\,000$, casting doubts on the existence of an MRI dynamo in the
$\mathrm{Pm}\leq 1$ regime.

[21] oai:arXiv.org:1509.04859 [pdf] - 1292361

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

Comments: 4 pages, 4 figures, accepted in A&A Letters

Submitted: **2015-09-16**

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

[22] oai:arXiv.org:1508.00904 [pdf] - 1288444

Magnetically driven accretion in protoplanetary discs

Comments: accepted to MNRAS after very minor revisions

Submitted: **2015-08-04**, last modified: **2015-09-04**

We characterize magnetically driven accretion at radii between 1 au and 100
au in protoplanetary discs, using a series of local non-ideal
magnetohydrodynamic (MHD) simulations. The simulations assume a Minimum Mass
Solar Nebula (MMSN) disc that is threaded by a net vertical magnetic field of
specified strength. Confirming previous results, we find that the Hall effect
has only a modest impact on accretion at 30 au, and essentially none at 100 au.
At 1-10 au the Hall effect introduces a pronounced bi-modality in the accretion
process, with vertical magnetic fields aligned to the disc rotation supporting
a strong laminar Maxwell stress that is absent if the field is anti-aligned. In
the anti-aligned case, we instead find evidence for bursts of turbulent stress
at 5-10 au, which we tentatively identify with the non-axisymmetric Hall-shear
instability. The presence or absence of these bursts depends upon the details
of the adopted chemical model, which suggests that appreciable regions of
actual protoplanetary discs might lie close to the borderline between laminar
and turbulent behaviour. Given the number of important control parameters that
have already been identified in MHD models, quantitative predictions for disc
structure in terms of only radius and accretion rate appear to be difficult.
Instead, we identify robust qualitative tests of magnetically driven accretion.
These include the presence of turbulence in the outer disc, independent of the
orientation of the vertical magnetic fields, and a Hall-mediated bi-modality in
turbulent properties extending from the region of thermal ionization to 10 au.

[23] oai:arXiv.org:1505.05661 [pdf] - 1245796

Angular momentum transport and large eddy simulations in
magnetorotational turbulence: the small Pm limit

Comments: Accepted for publication in A&A

Submitted: **2015-05-21**

Angular momentum transport in accretion discs is often believed to be due to
magnetohydrodynamic turbulence mediated by the magnetorotational instability.
Despite an abundant literature on the MRI, the parameters governing the
saturation amplitude of the turbulence are poorly understood and the existence
of an asymptotic behavior in the Ohmic diffusion regime is not clearly
established. We investigate the properties of the turbulent state in the small
magnetic Prandtl number limit. Since this is extremely computationally
expensive, we also study the relevance and range of applicability of the most
common subgrid scale models for this problem. Unstratified shearing boxes
simulations are performed both in the compressible and incompressible limits,
with a resolution up to 800 cells per disc scale height. The latter constitutes
the largest resolution ever attained for a simulation of MRI turbulence. In the
presence of a mean magnetic field threading the domain, angular momentum
transport converges to a finite value in the small Pm limit. When the mean
vertical field amplitude is such that {\beta}, the ratio between the thermal
and magnetic pressure, equals 1000, we find {\alpha}~0.032 when Pm approaches
zero. In the case of a mean toroidal field for which {\beta}=100, we find
{\alpha}~0.018 in the same limit. Both implicit LES and Chollet-Lesieur closure
model reproduces these results for the {\alpha} parameter and the power
spectra. A reduction in computational cost of a factor at least 16 (and up to
256) is achieved when using such methods. MRI turbulence operates efficiently
in the small Pm limit provided there is a mean magnetic field. Implicit LES
offers a practical and efficient mean of investigation of this regime but
should be used with care, particularly in the case of a vertical field.
Chollet-Lesieur closure model is perfectly suited for simulations done with a
spectral code.

[24] oai:arXiv.org:1410.3406 [pdf] - 934329

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

Comments: 7 pages, 6 figures

Submitted: **2014-10-13**

The magnetorotational (MRI) dynamo has long been considered one of the
possible drivers of turbulent angular momentum transport in astrophysical
accretion disks. However, various numerical results suggest that this dynamo
may be difficult to excite in the astrophysically relevant regime of magnetic
Prandtl number (Pm) significantly smaller than unity, for reasons currently not
well understood. The aim of this article is to present the first results of an
ongoing numerical investigation of the role of both linear and nonlinear
dissipative effects in this problem. Combining a parametric exploration and an
energy analysis of incompressible nonlinear MRI dynamo cycles representative of
the transitional dynamics in large aspect ratio shearing boxes, we find that
turbulent magnetic diffusion makes the excitation and sustainment of this
dynamo at moderate magnetic Reynolds number (Rm) increasingly difficult for
decreasing Pm. This results in an increase in the critical Rm of the dynamo for
increasing kinematic Reynolds number (Re), in agreement with earlier numerical
results. Given its very generic nature, we argue that turbulent magnetic
diffusion could be an important determinant of MRI dynamo excitation in disks,
and may also limit the efficiency of angular momentum transport by MRI
turbulence in low Pm regimes.

[25] oai:arXiv.org:1402.4133 [pdf] - 844794

Thanatology in Protoplanetary Discs: the combined influence of Ohmic,
Hall, and ambipolar diffusion on dead zones

Comments: 16 pages, 18 figures, 3 tables, accepted for publication in Astronomy
& Astrophysics

Submitted: **2014-02-17**, last modified: **2014-04-18**

Protoplanetary discs are poorly ionised due to their low temperatures and
high column densities, and are therefore subject to three "non-ideal"
magnetohydrodynamic effects: Ohmic dissipation, ambipolar diffusion, and the
Hall effect. The existence of magnetically driven turbulence in these discs has
been a central question since the discovery of the magnetorotational
instability. Early models considered Ohmic diffusion only and led to a scenario
of layered accretion, in which a magnetically "dead" zone in the disc midplane
is embedded within magnetically "active" surface layers at distances ~1-10 au
from the central protostellar object. Recent work has suggested that a
combination of Ohmic dissipation and ambipolar diffusion can render both the
midplane and surface layers of the disc inactive and that torques due to
magnetically driven outflows are required to explain the observed accretion
rates. We reassess this picture by performing three-dimensional numerical
simulations that include, for the first time, all three non-ideal MHD effects.
We find that the Hall effect can generically "revive" dead zones by producing a
dominant azimuthal magnetic field and a large-scale Maxwell stress throughout
the midplane, provided the angular velocity and magnetic field satisfy Omega.B
> 0. The attendant large magnetic pressure modifies the vertical density
profile and substantially increases the disc scale height beyond its
hydrostatic value. Outflows are produced, but are not necessary to explain
accretion rates <10^{-7} Msun/yr. The flow in the disc midplane is essentially
laminar, suggesting that dust sedimentation may be efficient. These results
demonstrate that, if the MRI is relevant for driving mass accretion in
protoplanetary discs, one must include the Hall effect to obtain even
qualitatively correct results.

[26] oai:arXiv.org:1401.7306 [pdf] - 1203034

Transport and Accretion in Planet-Forming Disks

Comments: Refereed review chapter accepted for publication in Protostars and
Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen,
C. Dullemond & Th. Henning. 24 pages, 7 figures. The talk presented at the
Protostars and Planets VI conference is available at
http://youtube.com/watch?v=tEgw0PXwkGE

Submitted: **2014-01-28**

Planets appear to form in environments shaped by the gas flowing through
protostellar disks to the central young stars. The flows in turn are governed
by orbital angular momentum transfer. In this chapter we summarize current
understanding of the transfer processes best able to account for the flows,
including magneto-rotational turbulence, magnetically-launched winds,
self-gravitational instability and vortices driven by hydrodynamical
instabilities. For each in turn we outline the major achievements of the past
few years and the outstanding questions. We underscore the requirements for
operation, especially ionization for the magnetic processes and heating and
cooling for the others. We describe the distribution and strength of the
resulting flows and compare with the long-used phenomenological
$\alpha$-picture, highlighting issues where the fuller physical picture yields
substantially different answers. We also discuss the links between magnetized
turbulence and magnetically-launched outflows, and between magnetized
turbulence and hydrodynamical vortices. We end with a summary of the status of
efforts to detect specific signatures of the flows.

[27] oai:arXiv.org:1306.5887 [pdf] - 1172278

Magnetic self-organisation in Hall-dominated magnetorotational
turbulence

Comments: 18 pages, 16 figures, accepted for publication in MNRAS

Submitted: **2013-06-25**

The magnetorotational instability (MRI) is the most promising mechanism by
which angular momentum is efficiently transported outwards in astrophysical
discs. However, its application to protoplanetary discs remains problematic.
These discs are so poorly ionised that they may not support magnetorotational
turbulence in regions referred to as `dead zones'. It has recently been
suggested that the Hall effect, a non-ideal magnetohydrodynamic (MHD) effect,
could revive these dead zones by enhancing the magnetically active column
density by an order of magnitude or more. We investigate this idea by
performing local, three-dimensional, resistive Hall-MHD simulations of the MRI
in situations where the Hall effect dominates over Ohmic dissipation. As
expected from linear stability analysis, we find an exponentially growing
instability in regimes otherwise linearly stable in resistive MHD. However,
instead of vigorous and sustained magnetorotational turbulence, we find that
the MRI saturates by producing large-scale, long-lived, axisymmetric structures
in the magnetic and velocity fields. We refer to these structures as zonal
fields and zonal flows, respectively. Their emergence causes a steep reduction
in turbulent transport by at least two orders of magnitude from extrapolations
based upon resistive MHD, a result that calls into question contemporary models
of layered accretion. We construct a rigorous mean-field theory to explain this
new behaviour and to predict when it should occur. Implications for
protoplanetary disc structure and evolution, as well as for theories of planet
formation, are briefly discussed.

[28] oai:arXiv.org:1306.3374 [pdf] - 1172039

Global bifurcations to subcritical magnetorotational dynamo action in
Keplerian shear flow

Comments: 43 pages, 17 figures, accepted for publication in J. Fluid Mech

Submitted: **2013-06-14**

Magnetorotational dynamo action in Keplerian shear flow is a
three-dimensional, nonlinear magnetohydrodynamic process whose study is
relevant to the understanding of accretion and magnetic field generation in
astrophysics. Transition to this form of dynamo is subcritical and shares many
characteristics of transition to turbulence in non-rotating hydrodynamic shear
flows. This suggests that these different fluid systems become active through
similar generic bifurcation mechanisms, which in both cases have eluded
detailed understanding so far. In this paper, we investigate numerically the
bifurcation mechanisms at work in the incompressible Keplerian
magnetorotational dynamo problem in the shearing box framework. Using numerical
techniques imported from dynamical systems research, we show that the onset of
chaotic dynamo action at magnetic Prandtl numbers larger than unity is
primarily associated with global homoclinic and heteroclinic bifurcations of
nonlinear magnetorotational dynamo cycles. These global bifurcations are
supplemented by local bifurcations of cycles marking the beginning of
period-doubling cascades. This suggests that nonlinear magnetorotational dynamo
cycles provide the pathway to turbulent injection of both kinetic and magnetic
energy in incompressible magnetohydrodynamic Keplerian shear flow in the
absence of an externally imposed magnetic field. Studying the nonlinear physics
and bifurcations of these cycles in different regimes and configurations may
subsequently help to better understand the conditions of excitation of
magnetohydrodynamic turbulence and instability-driven dynamos in various
astrophysical systems and laboratory experiments. The detailed characterization
of global bifurcations provided for this three-dimensional subcritical fluid
dynamics problem may also prove useful for the problem of transition to
turbulence in hydrodynamic shear flows.

[29] 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.

[30] oai:arXiv.org:1212.4045 [pdf] - 603678

Modelling colliding wind binaries with RAMSES, extension to special
relativity

Comments: Proceeding of Astronum 2012 Conference

Submitted: **2012-12-17**

We present high resolution simulations with RAMSES of supersonic colliding
stellar winds. The collision results in a double shock structure which is
subject to different instabilities. The Kelvin-Helmholtz instability (KHI)
introduces some mixing and variability. For isothermal winds, the Non-linear
Thin Shell Instability violently affects the interaction region. Properly
modelling these instabilities requires a high enough resolution and an adapted
numerical method, especially when one of the winds strongly dominates the other
one. At large scale, orbital motion is expected to turn the shocked zone into a
spiral but we find that in some configurations the KHI may disrupt the spiral.
A colliding wind structure is also expected in gamma-ray binaries composed of a
massive star and a young pulsar which emits a highly relativistic wind.
Numerical simulations are necessary to understand the geometry of such systems
and should take into account the relativistic nature of the pulsar wind. We
implemented a second order Godunov method to solve the equations of
relativistic hydrodynamics with RAMSES, including the possibility for Adaptive
Mesh Refinement. After a brief overview of our numerical implementation, we
will present preliminary simulation that contrast the structure of the flow in
gamma-ray and stellar binaries.

[31] oai:arXiv.org:1212.4142 [pdf] - 1158558

Colliding wind binaries and gamma-ray binaries : relativistic version of
the RAMSES code

Comments: Proceeding of the 5th International Symposium on High-Energy
Gamma-Ray Astronomy (Gamma2012). arXiv admin note: text overlap with
arXiv:1212.4045

Submitted: **2012-12-17**

Gamma-ray binaries are colliding wind binaries (CWB) composed of a massive
star a non-accreting pulsar with a highly relativistic wind. Particle
acceleration at the shocks results in emission going from extended radio
emission to the gamma-ray band. The interaction region is expected to show
common features with stellar CWB. Performing numerical simulations with the
hydrodynamical code RAMSES, we focus on their structure and stability and find
that the Kelvin-Helmholtz instability (KHI) can lead to important mixing
between the winds and destroy the large scale spiral structure. To investigate
the impact of the relativistic nature of the pulsar wind, we extend RAMSES to
relativistic hydrodynamics (RHD). Preliminary simulations of the interaction
between a pulsar wind and a stellar wind show important similarities with
stellar colliding winds with small relativistic corrections.

[32] oai:arXiv.org:1210.6660 [pdf] - 1152424

The magnetorotational instability as a jet launching mechanism

Comments: 14 pages, 17 figures, submitted to A&A

Submitted: **2012-10-24**

Magnetorotational turbulence and magnetically driven disc winds are often
considered as separate processes. However, realistic astrophysical discs are
expected to be subject to both effects, although possibly at different times
and locations. We investigate here the potential link between these two
phenomena using a mixed numerical and analytical approach. We show in
particular that large-scale MRI modes which dominate strongly magnetised discs
(plasma beta~10) naturally produce magnetically driven outflows in the
nonlinear regime. We show that these outflows share many similarities with
local and global disc wind solutions found in the literature. We also
investigate the 3D stability of these outflows and show that they are unstable
on dynamical timescales. The implications of these results for the transition
between a jet-emitting disc and a standard "viscous" disc are discussed.

[33] oai:arXiv.org:1202.2060 [pdf] - 1116515

Impact of orbital motion on the structure and stability of adiabatic
shocks in colliding wind binaries

Comments: Submitted to Astronomy and Astophysics

Submitted: **2012-02-09**

The collision of winds from massive stars in binaries results in the
formation of a double-shock structure with observed signatures from radio to
X-rays. We study the structure and stability of the colliding wind region as it
turns into a spiral due to orbital motion. We focus on adiabatic winds, where
mixing between the two winds is expected to be restricted to the
Kelvin-Helmholtz instability (KHI). Mixing of the Wolf-Rayet wind with
hydrogen-rich material is important for dust formation in pinwheel nebulae such
as WR 104, where the spiral structure has been resolved in infrared. We use the
hydrodynamical code RAMSES with an adaptive grid. A wide range of binary
systems with different wind velocities and mass loss rates are studied with 2D
simulations. A specific 3D simulation is performed to model WR 104. Orbital
motion leads to the formation of two distinct spiral arms where the KHI
develops differently. We find that the spiral structure is destroyed when there
is a large velocity gradient between the winds, unless the collimated wind is
much faster. We argue that the KHI plays a major role in maintaining or not the
structure. We discuss the consequences for various colliding wind binaries.
When stable, there is no straightforward relationship between the spatial step
of the spiral, the wind velocities, and the orbital period. Our 3D simulation
of WR 104 indicates that the colder, well-mixed trailing arm has more
favourable conditions for dust formation than the leading arm. The single-arm
infrared spiral follows more closely the mixing map than the density map,
suggesting the dust-to-gas ratio may vary between the leading and trailing
density spirals. However, the density is much lower than what dust formation
models require. Including radiative cooling would lead to higher densities, and
also to thin shell instabilities whose impact on the large structure remains
unknown.

[34] oai:arXiv.org:1201.5020 [pdf] - 1093205

On the interaction between tides and convection

Comments: 14 pages, 5 figures, to be published in MNRAS

Submitted: **2012-01-24**

We study the interaction between tides and convection in astrophysical bodies
by analysing the effect of a homogeneous oscillatory shear on a fluid flow.
This model can be taken to represent the interaction between a large-scale
periodic tidal deformation and a smaller-scale convective motion. We first
consider analytically the limit in which the shear is of low amplitude and the
oscillation period is short compared to the timescales of the unperturbed flow.
In this limit there is a viscoelastic response and we obtain expressions for
the effective elastic modulus and viscosity coefficient. The effective
viscosity is inversely proportional to the square of the oscillation frequency,
with a coefficient that can be positive, negative or zero depending on the
properties of the unperturbed flow. We also carry out direct numerical
simulations of Boussinesq convection in an oscillatory shearing box and measure
the time-dependent Reynolds stress. The results indicate that the effective
viscosity of turbulent convection falls rapidly as the oscillation frequency is
increased, attaining small negative values in the cases we have examined,
although significant uncertainties remain because of the turbulent noise. We
discuss the implications of this analysis for astrophysical tides.

[35] oai:arXiv.org:1109.1811 [pdf] - 1083929

Periodic magnetorotational dynamo action as a prototype of nonlinear
magnetic field generation in shear flows

Comments: 10 pages, 6 figures, accepted for publication in Physical Review E

Submitted: **2011-09-08**

The nature of dynamo action in shear flows prone to magnetohydrodynamic
instabilities is investigated using the magnetorotational dynamo in Keplerian
shear flow as a prototype problem. Using direct numerical simulations and
Newton's method, we compute an exact time-periodic magnetorotational dynamo
solution to the three-dimensional dissipative incompressible
magnetohydrodynamic equations with rotation and shear. We discuss the physical
mechanism behind the cycle and show that it results from a combination of
linear and nonlinear interactions between a large-scale axisymmetric toroidal
magnetic field and non-axisymmetric perturbations amplified by the
magnetorotational instability. We demonstrate that this large scale dynamo
mechanism is overall intrinsically nonlinear and not reducible to the standard
mean-field dynamo formalism. Our results therefore provide clear evidence for a
generic nonlinear generation mechanism of time-dependent coherent large-scale
magnetic fields in shear flows and call for new theoretical dynamo models.
These findings may offer important clues to understand the transitional and
statistical properties of subcritical magnetorotational turbulence.

[36] oai:arXiv.org:1012.2690 [pdf] - 1042522

Nonlinear energy transfers in accretion discs MRI turbulence. I-Net
vertical field case

Comments: 10 pages, 9 figures, accepted for publication in Astronomy &
Astrophysics

Submitted: **2010-12-13**

The magnetorotational instability (MRI) is believed to be responsible for
most of the angular momentum transport in accretion discs. However, molecular
dissipation processes may drastically change the efficiency of MRI turbulence
in realistic astrophysical situations. The physical origin of this dependency
is still poorly understood as linear and quasi linear theories fail to explain
it. In this paper, we look for the link between molecular dissipation processes
and MRI transport of angular momentum in non stratified shearing box
simulations including a mean vertical field. We show that magnetic helicity is
unimportant in the model we consider. We perform a spectral analysis on the
simulations tracking energy exchanges in spectral space when turbulence is
fully developed. We find that the energy exchanges are essentially direct (from
large to small scale) whereas some non linear interactions appear to be non
local in spectral space. We speculate that these non local interactions are
responsible for the correlation between turbulent transport and molecular
dissipation. We argue that this correlation should then disappear when a
significant scale separation is achieved and we discuss several methods by
which one can test this hypothesis.

[37] oai:arXiv.org:1009.5197 [pdf] - 1040852

Vortex migration in protoplanetary disks

Comments: 13 pages, 13 figures, accepted for publication in ApJ

Submitted: **2010-09-27**

We consider the radial migration of vortices in two-dimensional isothermal
gaseous disks. We find that a vortex core, orbiting at the local gas velocity,
induces velocity perturbations that propagate away from the vortex as density
waves. The resulting spiral wave pattern is reminiscent of an embedded planet.
There are two main causes for asymmetries in these wakes: geometrical effects
tend to favor the outer wave, while a radial vortensity gradient leads to an
asymmetric vortex core, which favors the wave at the side that has the lowest
density. In the case of asymmetric waves, which we always find except for a
disk of constant pressure, there is a net exchange of angular momentum between
the vortex and the surrounding disk, which leads to orbital migration of the
vortex. Numerical hydrodynamical simulations show that this migration can be
very rapid, on a time scale of a few thousand orbits, for vortices with a size
comparable to the scale height of the disk. We discuss the possible effects of
vortex migration on planet formation scenarios.

[38] oai:arXiv.org:1006.5669 [pdf] - 1033393

Supertransient magnetohydrodynamic turbulence in Keplerian shear flows

Comments: 4 pages, 5 figures, accepted for publication in Physical Review
Letters

Submitted: **2010-06-29**

A subcritical transition to turbulence in magnetized Keplerian shear flows is
investigated using a statistical approach. Three-dimensional numerical
simulations of the shearing box equations with zero net magnetic flux are
employed to determine the transition from decaying to sustained turbulence as a
function of the magnetic Reynolds number Rm. The results reveal no clear
transition to sustained turbulence as the average lifetime of the transients
grows as an exponential function of Rm, in accordance with a type-II
supertransient law.

[39] oai:arXiv.org:1005.2283 [pdf] - 1032441

Buoyancy-induced time delays in Babcock-Leighton flux-transport dynamo
models

Comments: 13 pages, 10 figures, accepted for publication in A&A

Submitted: **2010-05-13**

The Sun is a magnetic star whose cyclic activity is thought to be linked to
internal dynamo mechanisms. A combination of numerical modelling with various
levels of complexity is an efficient and accurate tool to investigate such
intricate dynamical processes. We investigate the role of the magnetic buoyancy
process in 2D Babcock-Leighton dynamo models, by modelling more accurately the
surface source term for poloidal field. Methods. To do so, we reintroduce in
mean-field models the results of full 3D MHD calculations of the non-linear
evolution of a rising flux tube in a convective shell. More specifically, the
Babcock-Leighton source term is modified to take into account the delay
introduced by the rise time of the toroidal structures from the base of the
convection zone to the solar surface. We find that the time delays introduced
in the equations produce large temporal modulation of the cycle amplitude even
when strong and thus rapidly rising flux tubes are considered. Aperiodic
modulations of the solar cycle appear after a sequence of period doubling
bifurcations typical of non-linear systems. The strong effects introduced even
by small delays is found to be due to the dependence of the delays on the
magnetic field strength at the base of the convection zone, the modulation
being much less when time delays remain constant. We do not find any
significant influence on the cycle period except when the delays are made
artificially strong. A possible new origin of the solar cycle variability is
here revealed. This modulated activity and the resulting butterfly diagram are
then more compatible with observations than what the standard Babcock-Leighton
model produces.

[40] oai:arXiv.org:1004.1384 [pdf] - 1026192

MRI-driven turbulent transport: the role of dissipation, channel modes
and their parasites

Comments: 11 pages, 9 figures. Accepted for publication in A&A

Submitted: **2010-04-08**

In the recent years, MRI-driven turbulent transport has been found to depend
in a significant way on fluid viscosity $\nu$ and resistivity $\eta$ through
the magnetic Prandtl number $Pm=\nu/\eta$. In particular, the transport
decreases with decreasing $Pm$; if persistent at very large Reynolds numbers,
this trend may lead to question the role of MRI-turbulence in YSO disks, whose
Prandtl number is usually very small. In this context, the principle objective
of the present investigation is to characterize in a refined way the role of
dissipation. Another objective is to characterize the effect of linear (channel
modes) and quasi-linear (parasitic modes) physics in the behavior of the
transport. These objectives are addressed with the help of a number of
incompressible numerical simulations. The horizontal extent of the box size has
been increased in order to capture all relevant (fastest growing) linear and
secondary parasitic unstable modes. The major results are the following: i- The
increased accuracy in the computation of transport averages shows that the
dependence of transport on physical dissipation exhibits two different regimes:
for $Pm \lesssim 1$, the transport has a power-law dependence on the magnetic
Reynolds number rather than on the Prandtl number; for $Pm > 1$, the data are
consistent with a primary dependence on $Pm$ for large enough ($\sim 10^3$)
Reynolds numbers. ii- The transport-dissipation correlation is not clearly or
simply related to variations of the linear modes growth rates. iii- The
existence of the transport-dissipation correlation depends neither on the
number of linear modes captured in the simulations, nor on the effect of the
parasitic modes on the saturation of the linear modes growth. iv- The transport
is usually not dominated by axisymmetric (channel) modes

[41] oai:arXiv.org:1002.4621 [pdf] - 1025371

On the angular momentum transport due to vertical convection in
accretion discs

Comments: 6 pages, 5 figures, accepted in MNRAS

Submitted: **2010-02-24**

The mechanism of angular momentum transport in accretion discs has long been
debated. Although the magnetorotational instability appears to be a promising
process, poorly ionized regions of accretion discs may not undergo this
instability. In this letter, we revisit the possibility of transporting angular
momentum by turbulent thermal convection. Using high-resolution spectral
methods, we show that strongly turbulent convection can drive outward angular
momentum transport at a rate that is, under certain conditions, compatible with
observations of discs. We find however that the angular momentum transport is
always much weaker than the vertical heat transport. These results indicate
that convection might be another way to explain global disc evolution, provided
that a sufficiently unstable vertical temperature profile can be maintained.

[42] oai:arXiv.org:0911.0663 [pdf] - 1018198

The subcritical baroclinic instability in local accretion disc models

Comments: 12 pages, 17 figures, Accepted in A&A

Submitted: **2009-11-03**, last modified: **2010-01-14**

(abridged) Aims: We present new results exhibiting a subcritical baroclinic
instability (SBI) in local shearing box models. We describe the 2D and 3D
behaviour of this instability using numerical simulations and we present a
simple analytical model describing the underlying physical process.
Results: A subcritical baroclinic instability is observed in flows stable for
the Solberg-Hoiland criterion using local simulations. This instability is
found to be a nonlinear (or subcritical) instability, which cannot be described
by ordinary linear approaches. It requires a radial entropy gradient weakly
unstable for the Schwartzchild criterion and a strong thermal diffusivity (or
equivalently a short cooling time). In compressible simulations, the
instability produces density waves which transport angular momentum outward
with typically alpha<3e-3, the exact value depending on the background
temperature profile. Finally, the instability survives in 3D, vortex cores
becoming turbulent due to parametric instabilities.
Conclusions: The subcritical baroclinic instability is a robust phenomenon,
which can be captured using local simulations. The instability survives in 3D
thanks to a balance between the 2D SBI and 3D parametric instabilities.
Finally, this instability can lead to a weak outward transport of angular
momentum, due to the generation of density waves by the vortices.

[43] oai:arXiv.org:1001.0109 [pdf] - 130353

The Validity of the Super-Particle Approximation during Planetesimal
Formation

Comments: 10 pages, 7 figures, accepted for publication in Astronomy and
Astrophysics

Submitted: **2009-12-31**, last modified: **2010-01-08**

The formation mechanism of planetesimals in protoplanetary discs is hotly
debated. Currently, the favoured model involves the accumulation of meter-sized
objects within a turbulent disc, followed by a phase of gravitational
instability. At best one can simulate a few million particles numerically as
opposed to the several trillion meter-sized particles expected in a real
protoplanetary disc. Therefore, single particles are often used as
super-particles to represent a distribution of many smaller particles. It is
assumed that small scale phenomena do not play a role and particle collisions
are not modeled. The super-particle approximation can only be valid in a
collisionless or strongly collisional system, however, in many recent numerical
simulations this is not the case.
In this work we present new results from numerical simulations of
planetesimal formation via gravitational instability. A scaled system is
studied that does not require the use of super-particles. We find that the
scaled particles can be used to model the initial phases of clumping if the
properties of the scaled particles are chosen such that all important
timescales in the system are equivalent to what is expected in a real
protoplanetary disc. Constraints are given for the number of particles needed
in order to achieve numerical convergence.
We compare this new method to the standard super-particle approach. We find
that the super-particle approach produces unreliable results that depend on
artifacts such as the gravitational softening in both the requirement for
gravitational collapse and the resulting clump statistics. Our results show
that short range interactions (collisions) have to be modelled properly.

[44] oai:arXiv.org:0907.1393 [pdf] - 1527982

Turbulent resistivity evaluation in MRI generated turbulence

Comments: 11 pages, 8 figures, accepted for publication in Astronomy and
Astrophysics

Submitted: **2009-07-08**

(abriged) MRI turbulence is a leading mechanism for the generation of an
efficient turbulent transport of angular momentum in an accretion disk through
a turbulent viscosity effect. It is believed that the same process could also
transport large-scale magnetic fields in disks, reshaping the magnetic
structures in these objects. This process, known as turbulent resistivity, has
been suggested and used in several accretion-ejection models and simulations to
produce jets. Still, the efficiency of MRI-driven turbulence to transport
large-scale magnetic fields is largely unknown.
We investigate this problem both analytically and numerically. We introduce a
linear calculation of the MRI in the presence of a spatially inhomogeneous mean
magnetic field. We show that, in this configuration, MRI modes lead to an
efficient magnetic field transport, on the order of the angular momentum
transport. We next use fully non linear simulations of MRI turbulence to
compute the turbulent resistivity in several magnetic configurations.
We find that the turbulent resistivity is on the order of the turbulent
viscosity in all our simulations, although somewhat lower. The turbulent
resistivity tensor is found to be highly anisotropic with a diffusion
coefficient 3 times greater in the radial direction than in the vertical
direction.
These results support the possibility of driving jets from turbulent disks;
the resulting jets may not be steady.

[45] oai:arXiv.org:0903.1720 [pdf] - 22225

On the Stability of Elliptical Vortices in Accretion Discs

Comments: 13 pages, 14 figures, accepted for publication by A&A. Animations
available on author's website

Submitted: **2009-03-10**

(Abriged) The existence of large-scale and long-lived 2D vortices in
accretion discs has been debated for more than a decade. They appear
spontaneously in several 2D disc simulations and they are known to accelerate
planetesimal formation through a dust trapping process. However, the issue of
the stability of these structures to the imposition of 3D disturbances is still
not fully understood, and it casts doubts on their long term survival.
Aim: We present new results on the 3D stability of elliptical vortices
embedded in accretion discs, based on a linear analysis and several non-linear
simulations.
Methods: We derive the linearised equations governing the 3D perturbations in
the core of an elliptical vortex, and we show that they can be reduced to a
Floquet problem. We solve this problem numerically in the astrophysical regime
and we present several analytical limits for which the mechanism responsible
for the instability can be explained. Finally, we compare the results of the
linear analysis to some high resolution simulations.
Results: We show that most anticyclonic vortices are unstable due to a
resonance between the turnover time and the local epicyclic oscillation period.
In addition, we demonstrate that a strong vertical stratification does not
create any additional stable domain of aspect ratio, but it significantly
reduces growth rates for relatively weak (and therefore elongated) vortices.
Conclusions: Elliptical vortices are always unstable, whatever the horizontal
or vertical aspect-ratio is. The instability can however be weak and is often
found at small scales, making it difficult to detect in low-order
finite-difference simulations.

[46] oai:arXiv.org:0809.3871 [pdf] - 16612

Localized magnetorotational instability and its role in the accretion
disc dynamo

Comments: 15 pages, 5 figures

Submitted: **2008-09-23**

(Abriged) The magnetorotational instability (MRI) is believed to be an
efficient way to transport angular momentum in accretion discs. It has also
been suggested as a way to amplify magnetic fields in discs, the instability
acting as a nonlinear dynamo. Recent numerical work has shown that a
large-scale magnetic field, which is predominantly azimuthal, can be sustained
by motions driven by the MRI of this same field. Following this idea, we
present an analytical calculation of the MRI in the presence of an azimuthal
field with a non-trivial vertical structure. We find that the mean radial EMF
associated to MRI modes tends to reduce the magnetic energy, acting like a
turbulent resistivity by mixing the non-uniform azimuthal field. Meanwhile, the
azimuthal EMF generates a radial field that, in combination with the Keplerian
shear, tends to amplify the azimuthal field and can therefore assist in the
dynamo process. This effect, however, is reversed for sufficiently strong
azimuthal fields, naturally leading to a saturation of the dynamo and possibly
to a cyclic behaviour of the magnetic field, as found in previous numerical
works.

[47] oai:arXiv.org:0807.1703 [pdf] - 14400

On Self-Sustained Dynamo Cycles in Accretion Discs

Comments: 11 pages, 11 figures, Accepted in A&A

Submitted: **2008-07-10**

(abridged) MHD turbulence is known to exist in shearing boxes with either
zero or nonzero net magnetic flux. However, the way turbulence survives in the
zero-net-flux case is not explained by linear theory and appears as a purely
numerical result. Aims: We look for a nonlinear mechanism able to explain the
persistence of MHD turbulence in shearing boxes with zero net magnetic flux,
and potentially leading to large-scale dynamo action. Method: Spectral
nonlinear simulations of the magnetorotational instability are shown to exhibit
a large-scale axisymmetric magnetic field, maintained for a few orbits. The
generation process of this field is investigated using the results of the
simulations and an inhomogeneous linear approach. Results: The mechanism by
which turbulence is sustained in zero-net-flux shearing boxes is shown to be
related to the existence of a large-scale azimuthal field, surviving for
several orbits. In particular, it is shown that MHD turbulence in shearing
boxes can be seen as a dynamo process coupled to a magnetorotational-type
instability.

[48] oai:arXiv.org:0705.3622 [pdf] - 1585

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

Comments: 10 pages, 13 figures, accepted in A&A. Numerical results improved,
minor changes in the text

Submitted: **2007-05-24**, last modified: **2007-10-04**

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

[49] oai:arXiv.org:0704.2943 [pdf] - 652

Impact of dimensionless numbers on the efficiency of MRI-induced
turbulent transport

Comments: 11 pages, 11 figures, accepted to MNRAS

Submitted: **2007-04-23**

The magneto-rotational instability is presently the most promising source of
turbulent transport in accretion disks. However, some important issues still
need to be addressed to quantify the role of MRI in disks; in particular no
systematic investigation of the role of the physical dimensionless parameters
of the problem on the dimensionless transport has been undertaken yet. First,
we complete existing investigations on the field strength dependence by showing
that the transport in high magnetic pressure disks close to marginal stability
is highly time-dependent and surprisingly efficient. Second, we bring to light
a significant dependence of the global transport on the magnetic Prandtl
number, with $\alpha\propto Pm^\delta$ for the explored range: $0.12<Pm<8$ and
$200<Re<6400$ ($\delta$ being in the range 0.25 to 0.5). We show that the
dimensionless transport is not correlated to the dimensionless linear growth
rate, contrarily to a largely held expectation. More generally, these results
stress the need to control dissipation processes in astrophysical simulations.

[50] oai:arXiv.org:astro-ph/0509541 [pdf] - 76122

On the relevance of subcritical hydrodynamic turbulence to accretion
disk transport

Comments: 21 pages, 17 figures, accepted in Astronomy and Astrophysics

Submitted: **2005-09-19**

Hydrodynamic unstratified keplerian flows are known to be linearly stable at
all Reynolds numbers, but may nevertheless become turbulent through nonlinear
mechanisms. However, in the last ten years, conflicting points of view have
appeared on this issue. We have revisited the problem through numerical
simulations in the shearing sheet limit. It turns out that the effect of the
Coriolis force in stabilizing the flow depends on whether the flow is cyclonic
(cooperating shear and rotation vorticities) or anticyclonic (competing shear
and rotation vorticities); keplerian flows are anticyclonic. We have obtained
the following results: i/ The Coriolis force does not quench turbulence in
subcritical flows; ii/ The resolution demand, when moving away from the
marginal stability boundary, is much more severe for anticyclonic flows than
for cyclonic ones. Presently available computer resources do not allow
numerical codes to reach the keplerian regime. iii/ The efficiency of turbulent
transport is directly correlated to the Reynolds number of transition to
turbulence $Rg$, in such a way that the Shakura-Sunyaev parameter $\alpha\sim
1/Rg$. iv/ Even the most optimistic extrapolations of our numerical data show
that subcritical turbulent transport would be too inefficient in keplerian
flows by several orders of magnitude for astrophysical purposes. v/ Our results
suggest that the data obtained for keplerian-like flows in a Taylor-Couette
settings are largely affected by secondary flows, such as Ekman circulation.