
We study the effects of rotation on the growth and saturation of the
doublediffusive fingering (thermohaline) instability at low Prandtl number.
Using direct numerical simulations, we estimate the compositional transport
rates as a function of the relevant nondimensional parameters  the Rossby
number, inversely proportional to the rotation rate, and the density ratio
which measures the relative thermal and compositional stratifications. Within
our explored range of parameters, we generally find rotation to have little
effect on vertical transport. However, we also present one exceptional case
where a cyclonic large scale vortex (LSV) is observed at low density ratio and
fairly low Rossby number. The LSV leads to significant enhancement in the
fingering transport rates by concentrating compositionally dense downflows at
its core. We argue that the formation of such LSVs could be relevant to solving
the missing mixing problem in RGB stars.

This paper continues the systematic investigation of diffusive shear
instabilities initiated in Part I of this series. In this work, we primarily
focus on quantifying the impact of nonlocal mixing, which is not taken into
account in Zahn's mixing model (Zahn 1992). To do that, we present the results
of direct numerical simulations using the PADDI code with a new model setup
designed to contain both laminar and turbulent shear layers. As in Part I, we
use the Low P\'eclet Number approximation of Ligni\`eres (1999) to model the
evolution of the perturbations. Our main findings are twofold, and have
significant impact on our understanding of the role of diffusive shear
instabilities in stars. The first is that turbulence is not necessarily
generated whenever $J{\rm Pr} < (J{\rm Pr})_c$ which Zahn's criterion assumes
(Zahn 1974). We have demonstrated that the presence or absence of turbulent
mixing in this limit hysteretically depends on the history of the shear layer.
The second finding is that Zahn's nonlinear instability criterion only
approximately locates the edge of the turbulent layer, since mixing beyond the
region where $J{\rm Pr} < (J{\rm Pr})_c$ can also take place in a manner
analogous to convective overshoot. We found that the turbulent kinetic energy
decays roughly exponentially beyond the edge of the shearunstable region, on a
lengthscale $\delta$ that is directly proportional to the scale of the
turbulent eddies, which are themselves of the order of the Zahn scale (see Part
I). Both of these results suggest that mixing by diffusive shear instabilities
should be modeled with more care than is currently standard in stellar
evolution codes.

We present Direct Numerical Simulations of the transport of heat and heavy
elements across a doublediffusive interface or a doublediffusive staircase,
in conditions that are close to those one may expect to find near the boundary
between the heavyelement rich core and the hydrogenhelium envelope of giant
planets such as Jupiter. We find that the nondimensional ratio of the buoyancy
flux associated with heavy element transport to the buoyancy flux associated
with heat transport lies roughly between 0.5 and 1, which is much larger than
previous estimates derived by analogy with geophysical doublediffusive
convection. Using these results in combination with a coreerosion model
proposed by Guillot et al. (2004), we find that the entire core of Jupiter
would be eroded within less than 1Myr assuming that the coreenvelope boundary
is composed of a single interface. We also propose an alternative model that is
more appropriate in the presence of a wellestablished doublediffusive
staircase, and find that in this limit a large fraction of the core could be
preserved. These findings are interesting in the context of Juno's recent
results, but call for further modeling efforts to better understand the process
of core erosion from first principles.

Shearinduced turbulence could play a significant role in mixing momentum and
chemical species in stellar radiation zones, as discussed by Zahn (1974). In
this paper we analyze the results of direct numerical simulations of stratified
plane Couette flows, in the limit of rapid thermal diffusion, to measure the
turbulent diffusivity and turbulent viscosity as a function of the local shear
and the local stratification. We find that the stability criterion proposed by
Zahn (1974), namely that the product of the gradient Richardson number and the
Prandtl number must be smaller than a critical values $(J\Pr)_c$ for
instability, adequately accounts for the transition to turbulence in the flow,
with $(J\Pr)_c \simeq 0.007$. This result recovers and confirms the prior
findings of Prat et al. (2016). Zahn's model for the turbulent diffusivity and
viscosity (Zahn 1992), namely that the mixing coefficient should be
proportional to the ratio of the thermal diffusivity to the gradient Richardson
number, does not satisfactorily match our numerical data when applied as is. It
fails (as expected) in the limit of large stratification where the Richardson
number exceeds the aforementioned threshold for instability, but it also fails
in the limit of low stratification where the turbulent eddy scale becomes
limited by the computational domain size. We propose a revised model for
turbulent mixing by diffusive stratified shear instabilities, that now properly
accounts for both limits, fits our data satisfactorily, and recovers Zahn's
1992 model in the limit of large Reynolds numbers.

When particles settle through a stable temperature or salinity gradient they
can drive an instability known as sedimentary fingering convection. This
phenomenon is thought to occur beneath sedimentrich river plumes in lakes and
oceans, in the context of marine snow where decaying organic materials serve as
the suspended particles, or in the atmosphere in the presence of aerosols or
volcanic ash. Laboratory experiments of Houk and Green (1973) and Green (1987)
have shown sedimentary fingering convection to be similar to the more commonly
known thermohaline fingering convection in many ways. Here, we study the
phenomenon using 3D direct numerical simulations. We find evidence for layer
formation in sedimentary fingering convection in regions of parameter space
where it does not occur for nonsedimentary systems. This is due to two
complementary effects. Sedimentation affects the turbulent fluxes and broadens
the region of parameter space unstable to the $\gamma$instability (Radko 2003)
to include systems at larger density ratios. It also gives rise to a new
layering instability that exists in $\gamma$stable regimes. The former is
likely quite ubiquitous in geophysical systems for sufficiently large settling
velocities, while the latter probably grows too slowly to be relevant, at least
in the context of sediments in water.

Fingering convection (or thermohaline convection) is a weak yet important
kind of mixing that occurs in stablystratified stellar radiation zones in the
presence of an inverse meanmolecularweight gradient. Brown et al. (2013)
recently proposed a new model for mixing by fingering convection, which
contains no free parameter, and was found to fit the results of direct
numerical simulations in almost all cases. Notably, however, they found that
mixing was substantially enhanced above their predicted values in the few cases
where largescale gravity waves, followed by thermocompositional layering,
grew spontaneously from the fingering convection. This effect is wellknown in
the oceanographic context, and is attributed to the excitation of the socalled
"collective instability". In this work, we build on the results of Brown et al.
(2013) and of Traxler et al. (2011b) to determine the conditions under which
the collective instability may be expected. We find that it is only relevant in
stellar regions which have a relatively large Prandtl number (the ratio of the
kinematic viscosity to the thermal diffusivity), $O(10^{3})$ or larger. This
implies that the collective instability cannot occur in main sequence stars,
where the Prandtl number is always much smaller than this (except in the outer
layers of surface convection zones where fingering is irrelevant anyway). It
could in principle be excited in regions of high electron degeneracy, during He
core flash, or in the interiors of white dwarfs. We discuss the implications of
our findings for these objects, both from a theoretical and from an
observational point of view.

Doublediffusive instabilities are often invoked to explain enhanced
transport in stablystratified fluids. The moststudied natural manifestation
of this process, fingering convection, commonly occurs in the ocean's
thermocline and typically increases diapycnal mixing by two orders of magnitude
over molecular diffusion. Fingering convection is also often associated with
structures on much larger scales, such as thermohaline intrusions, gravity
waves and thermohaline staircases. In this paper, we present an exhaustive
study of the phenomenon from small to large scales. We perform the first
threedimensional simulations of the process at realistic values of the heat
and salt diffusivities and provide accurate estimates of the induced turbulent
transport. Our results are consistent with oceanic field measurements of
diapycnal mixing in fingering regions. We then develop a generalized meanfield
theory to study the stability of fingering systems to largescale
perturbations, using our calculated turbulent fluxes to parameterize
smallscale transport. The theory recovers the intrusive instability, the
collective instability, and the gammainstability as limiting cases. We find
that the fastestgrowing largescale mode depends sensitively on the ratio of
the background gradients of temperature and salinity (the density ratio). While
only intrusive modes exist at high density ratios, the collective and
gammainstabilities dominate the system at the low density ratios where
staircases are typically observed. We conclude by discussing our findings in
the context of staircase formation theory.

Regions of the ocean's thermocline unstable to salt fingering are often
observed to host thermohaline staircases, stacks of deep wellmixed convective
layers separated by thin stablystratified interfaces. Decades after their
discovery, however, their origin remains controversial. In this paper we use 3D
direct numerical simulations to shed light on the problem. We study the
evolution of an analogous doublediffusive system, starting from an initial
statistically homogeneous fingering state and find that it spontaneously
transforms into a layered state. By analysing our results in the light of the
meanfield theory developed in Paper I, a clear picture of the sequence of
events resulting in the staircase formation emerges. A collective instability
of homogeneous fingering convection first excites a field of gravity waves,
with a welldefined vertical wavelength. However, the waves saturate early
through regular but localized breaking events, and are not directly responsible
for the formation of the staircase. Meanwhile, slowergrowing, horizontally
invariant but vertically quasiperiodic gammamodes are also excited and grow
according to the gammainstability mechanism. Our results suggest that the
nonlinear interaction between these various meanfield modes of instability
leads to the selection of one particular gammamode as the staircase
progenitor. Upon reaching a critical amplitude, this progenitor overturns into
a fullyformed staircase. We conclude by extending the results of our
simulations to real oceanic parameter values, and find that the progenitor
gammamode is expected to grow on a timescale of a few hours, and leads to the
formation of a thermohaline staircase in about one day with an initial spacing
of the order of one to two metres.

The solar convection zone exhibits a strong level of differential rotation,
whereby the rotation period of the polar regions is about 2530% longer than
the equatorial regions. The Coriolis force associated with these zonal flows
perpetually "pumps" the convection zone fluid, and maintains a quasisteady
circulation, poleward near the surface. What is the influence of this
meridional circulation on the underlying radiative zone, and in particular,
does it provide a significant source of mixing between the two regions? In
Paper I, we began to study this question by assuming a fixed meridional flow
pattern in the convection zone and calculating its penetration depth into the
radiative zone. We found that the amount of mixing caused depends very
sensitively on the assumed flow structure near the radiativeconvective
interface. We continue this study here by including a simple model for the
convection zone "pump", and calculating in a selfconsistent manner the
meridional flows generated in the whole Sun. We find that the global
circulation timescale depends in a crucial way on two factors: the overall
stratification of the radiative zone as measured by the Rossby number times the
square root of the Prandtl number, and, for weakly stratified systems, the
presence or absence of stresses within the radiative zone capable of breaking
the TaylorProudman constraint. We conclude by discussing the consequences of
our findings for the solar interior and argue that a potentially important
mechanism for mixing in Main Sequence stars has so far been neglected.

The rotation rate of the solar radiative zone is an important diagnostic for
angularmomentum transport in the tachocline and below. In this paper we study
the contribution of viscous and magnetic stresses to the global
angularmomentum balance. By considering a simple linearized toy model, we
discuss the effects of field geometry and applied boundary conditions on the
predicted rotation profile and rotation rate of the radiative interior. We
compare these analytical predictions with fully nonlinear simulations of the
dynamics of the radiative interior, as well as with observations. We discuss
the implications of these results as constraints on models of the solar
interior.

We examine a physical process that leads to the efficient formation of gas
giant planets around intermediate mass stars. In the gaseous protoplanetary
disks surrounding rapidlyaccreting intermediatemass stars we show that the
midplane temperature (heated primarily by turbulent dissipation) can reach >
1000 K out to 1 AU. Thermal ionization of this hot gas couples the disk to the
magnetic field, allowing the magnetorotational instability (MRI) to generate
turbulence and transport angular momentum. Further from the central star the
ionization fraction decreases, decoupling the disk from the magnetic field and
reducing the efficiency of angular momentum transport. As the disk evolves
towards a quasisteady state, a local maximum in the surface density and in the
midplane pressure both develop at the inner edge of the MRIdead zone, trapping
inwardly migrating solid bodies. Small particles accumulate and coagulate into
planetesimals which grow rapidly until they reach isolation mass. In contrast
to the situation around solar type stars, we show that the isolation mass for
cores at this critical radius around the more massive stars is large enough to
promote the accretion of significant amounts of gas prior to disk depletion.
Through this process, we anticipate a prolific production of gas giants at ~1
AU around intermediatemass stars.

Meridional flows with velocities of a few meters per second are observed in
the uppermost regions of the solar convection zone. The amplitude and pattern
of the flows deeper in the solar interior, in particular near the top of the
radiative region, are of crucial importance to a wide range of solar
magnetohydrodynamical processes. In this paper, we provide a systematic study
of the penetration of largescale meridional flows from the convection zone
into the radiative zone. In particular, we study the effects of the assumed
boundary conditions applied at the convectiveradiative interface on the deeper
flows. Using simplified analytical models in conjunction with more complete
numerical methods, we show that penetration of the convectivelydriven
meridional flows into the deeper interior is not necessarily limited to a
shallow Ekman depth but can penetrate much deeper, depending on how the
convectiveradiative interface flows are modeled.

This series of papers investigates the early stages of planet formation by
modeling the evolution of the gas and solid content of protostellar disks from
the early T Tauri phase until complete dispersal of the gas. In this first
paper, I present a new set of simplified equations modeling the growth and
migration of various species of grains in a gaseous protostellar disk evolving
as a result of the combined effects of viscous accretion and photoevaporation
from the central star. Using the assumption that the grain size distribution
function always maintains a powerlaw structure approximating the average
outcome of the exact coagulation/shattering equation, the model focuses on the
calculation of the growth rate of the largest grains only. The coupled
evolution equations for the maximum grain size, the surface density of the gas
and the surface density of solids are then presented and solved
selfconsistently using a standard 1+1 dimensional formalism. I show that the
global evolution of solids is controlled by a leaky reservoir of small grains
at large radii, and propose an empirically derived evolution equation for the
total mass of solids, which can be used to estimate the total heavy element
retention efficiency in the planet formation paradigm. Consistency with
observation of the total mass of solids in the Minimum Solar Nebula augmented
with the mass of the Oort cloud sets strong upper limit on the initial grain
size distribution, as well as on the turbulent parameter $\alphat$. Detailed
comparisons with SED observations are presented in a following paper.

We study the interaction between gas and dust particles in a protoplanetary
disk, comparing analytical and numerical results. We first calculate
analytically the trajectories of individual particles undergoing gas drag in
the disk, in the asymptotic cases of very small particles (Epstein regime) and
very large particles (Stokes regime). Using a Boltzmann averaging method, we
then infer their collective behavior. We compare the results of this analytical
formulation against numerical computations of a large number of particles.
Using successive moments of the Boltzmann equation, we derive the equivalent
fluid equations for the average motion of the particles; these are
intrinsically different in the Epstein and Stokes regimes. We are also able to
study analytically the temporal evolution of a collection of particles with a
given initial sizedistribution provided collisions are ignored.

A quasisteady state model of the consequences of rotation on the
hydrodynamical structure of a stellar radiative zone is derived, by studying in
particular the role of centrifugal and baroclinic driving of meridional motions
in angularmomentum transport. This nonlinear problem is solved numerically
assuming axisymmetry of the system, and within some limits, it is shown that
there exist simple analytical solutions. The limit of slow rotation recovers
EddingtonSweet theory, whereas it is shown that in the limit of rapid
rotation, the system settles into a geostrophic equilibrium. The behaviour of
the system is found to be controlled by one parameter only, linked to the
Prantl number, the stratification and the rotation rate of the star.

Gough & McIntyre have suggested that the dynamics of the solar tachocline are
dominated by the advectiondiffusion balance between the differential rotation,
a largescale primordial field and baroclinicly driven meridional motions. This
paper presents the first part of a study of the tachocline, in which a model of
the rotation profile below the convection zone is constructed along the lines
suggested by Gough & McIntyre and solved numerically. In this first part, a
reduced model of the tachocline is derived in which the effects of
compressibility and energy transport on the system are neglected; the
meridional motions are driven instead by EkmanHartmann pumping. It is shown
that there exists only a narrow range of magnetic field strengths for which the
system can achieve a nearly uniform rotation. The results are discussed with
respect to observations and to the limitations of this initial approach. A
following paper combines the effects of realistic baroclinic driving and
stratification with a model that follows closely the lines of work of Gough &
McIntyre.