The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle--particle collisions allows departures from local thermodynamic equilibrium. These departures exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present a numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart. If instead the parallel viscous stress is left unabated -- a situation relevant to recent kinetic simulations of the fluctuation dynamo and to the early stages of the dynamo in a magnetized ICM -- the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics of the saturated state of the large-Prandtl-number (${\rm Pm}\gtrsim{1}$) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches magnetic-energy spectra. A prediction of this model, confirmed by our simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. In all cases that do exhibit a dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $\delta n/n_0 \gtrsim 2/\beta$ (where $\beta\gg{1}$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves' otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more `fluid-like' than kinetic. A nonlinear fluctuation-dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels, and weakly collisional thermodynamics when IAWs with $\delta n/n_0 \gtrsim 2/\beta$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

We develop simple, physically motivated models for drag-induced dust-gas streaming instabilities, which are thought to be crucial for clumping grains to form planetesimals in protoplanetary disks. The models explain, based on the physics of gaseous epicyclic motion and dust-gas drag forces, the most important features of the streaming instability and its simple generalisation, the disk settling instability. Some of the key properties explained by our models include the sudden change in the growth rate of the streaming instability when the dust-to-gas-mass ratio surpasses one, the slow growth rate of the streaming instability compared to the settling instability for smaller grains, and the main physical processes underlying the growth of the most unstable modes in different regimes. As well as providing helpful simplified pictures for understanding the operation of an interesting and fundamental astrophysical fluid instability, our models may prove useful for analysing simulations and developing nonlinear theories of planetesimal growth in disks.

Cosmic rays (CRs) with ~GeV energies can contribute significantly to the energy and pressure budget in the interstellar, circumgalactic, and intergalactic medium (ISM, CGM, IGM). Recent cosmological simulations have begun to explore these effects, but almost all studies have been restricted to simplified models with constant CR diffusivity and/or streaming speeds. Physical models of CR propagation/scattering via extrinsic turbulence and self-excited waves predict transport coefficients which are complicated functions of local plasma properties. In a companion paper, we consider a wide range of observational constraints to identify proposed physically-motivated cosmic-ray propagation scalings which satisfy both detailed Milky Way (MW) and extra-galactic $\gamma$-ray constraints. Here, we compare the effects of these models relative to simpler 'diffusion+streaming' models on galaxy and CGM properties at dwarf through MW mass scales. The physical models predict large local variations in CR diffusivity, with median diffusivity increasing with galacto-centric radii and decreasing with galaxy mass and redshift. These effects lead to a more rapid dropoff of CR energy density in the CGM (compared to simpler models), in turn producing weaker effects of CRs on galaxy star formation rates (SFRs), CGM absorption profiles and galactic outflows. The predictions of the more physical CR models tend to lie 'in between' models which ignore CRs entirely and models which treat CRs with constant diffusivity.

Weakly collisional, magnetised plasmas characterised by anisotropic viscosity and conduction are ubiquitous in galaxies, halos and the intracluster medium (ICM). Cosmic rays (CRs) play an important role in these environments as well, by providing additional pressure and heating to the thermal plasma. We carry out a linear stability analysis of weakly collisional plasmas with cosmic rays using Braginskii MHD for the thermal gas. We assume that the CRs stream at the Alfv\'en speed, which in a weakly collisional plasma depends on the pressure anisotropy ($\Delta p$) of the thermal plasma. We find that this $\Delta p$-dependence introduces a phase shift between the CR-pressure and gas-density fluctuations. This drives a fast-growing acoustic instability: CRs offset the damping of acoustic waves by anisotropic viscosity and give rise to wave growth when the ratio of CR pressure to gas pressure is $\gtrsim \alpha \beta^{-1/2}$, where $\beta$ is the ratio of thermal to magnetic pressure, and $\alpha$, typically $\lesssim 1$, depends on other dimensionless parameters. In high-$\beta$ environments like the ICM, this condition is satisfied for small CR pressures. We speculate that the instability studied here may contribute to the scattering of high-energy CRs and to the excitation of sound waves in galaxy-halo, group and cluster plasmas, including the long-wavelength X-ray fluctuations in \textit{Chandra} observations of the Perseus cluster. It may also be important in the vicinity of shocks in dilute plasmas (e.g., cluster virial shocks or galactic wind termination shocks), where the CR pressure is locally enhanced.

The microphysics of ~GeV cosmic ray (CR) transport on galactic scales remain deeply uncertain, with almost all studies adopting simple prescriptions (e.g. constant-diffusivity). We explore different physically-motivated, anisotropic, dynamical CR transport scalings in high-resolution cosmological FIRE simulations of dwarf and ~$L_{\ast}$ galaxies where scattering rates vary with local plasma properties motivated by extrinsic turbulence (ET) or self-confinement (SC) scenarios, with varying assumptions about e.g. turbulent power spectra on un-resolved scales, Alfven-wave damping, etc. We self-consistently predict observables including $\gamma$-rays ($L_{\gamma}$), grammage, residence times, and CR energy densities to constrain the models. We demonstrate many non-linear dynamical effects (not captured in simpler models) tend to enhance confinement. For example, in multi-phase media, even allowing arbitrary fast transport in neutral gas does not substantially reduce CR residence times (or $L_{\gamma}$), as transport is rate-limited by the ionized WIM and 'inner CGM' gaseous halo ($10^{4}-10^{6}$ K gas within 10-30 kpc), and $L_{\gamma}$ can be dominated by trapping in small 'patches.' Most physical ET models contribute negligible scattering of ~1-10 GeV CRs, but it is crucial to account for anisotropy and damping (especially of fast modes) or else scattering rates would violate observations. We show that the most widely-assumed scalings for SC models produce excessive confinement by factors >100 in the WIM and inner CGM, where turbulent and Landau damping dominate. This suggests either a breakdown of quasi-linear theory used to derive the CR transport parameters in SC, or that other novel damping mechanisms dominate in intermediate-density ionized gas.

Recent near-sun solar-wind observations from Parker Solar Probe have found a highly dynamic magnetic environment, permeated by abrupt radial-field reversals, or "switchbacks." We show that many features of the observed turbulence are reproduced by a spectrum of Alfv\'enic fluctuations advected by a radially expanding flow. Starting from simple superpositions of low-amplitude outward-propagating waves, our expanding-box compressible MHD simulations naturally develop switchbacks because (i) the normalized amplitude of waves grows due to expansion and (ii) fluctuations evolve towards spherical polarization (i.e., nearly constant field strength). These results suggest that switchbacks form in-situ in the expanding solar wind and are not indicative of impulsive processes in the chromosphere or corona.

We explore the survival of cool clouds in multi-phase circum-galactic media. We revisit the "cloud crushing problem" in a large survey of simulations including radiative cooling, self-shielding, self-gravity, magnetic fields, and anisotropic Braginskii conduction and viscosity (with saturation). We explore a wide range of parameters including cloud size, velocity, ambient temperature and density, as well as a variety of magnetic field configurations and cloud turbulence. We find that realistic magnetic fields and turbulence have weaker effects on cloud survival; the most important physics is radiative cooling and conduction. Self-gravity and self-shielding are important for clouds which are initially Jeans-unstable, but largely irrelevant otherwise. Non-self-gravitating, realistically magnetized clouds separate into four regimes: (1) At low column densities, clouds evaporate rapidly via conduction. (2) A "failed pressure confinement" regime, where the ambient hot gas cools too rapidly to provide pressure confinement for the cloud. (3) An "infinitely long-lived" regime, in which the cloud lifetime becomes longer than the cooling time of gas swept up in the leading bow shock, so the cloud begins to accrete and grow. (4) A "classical cloud destruction" regime, where clouds are eventually destroyed by instabilities. In the final regime, the cloud lifetime can exceed the naive cloud-crushing time owing to conduction-induced compression. However, small and/or slow-moving clouds can also evaporate more rapidly than the cloud-crushing time. We develop simple analytic models that explain the simulated cloud destruction times in this regime.

We study the non-linear evolution of the acoustic 'Resonant Drag Instability' (RDI) using numerical simulations. The acoustic RDI is excited in a dust-gas mixture when dust grains stream through gas, interacting with sound waves to cause a linear instability. We study this process in a periodic box by accelerating neutral dust with an external driving force. The instability grows as predicted by linear theory, eventually breaking into turbulence and saturating. As in linear theory, the non-linear behavior is characterized by three regimes - high, intermediate, and low wavenumbers - the boundary between which is determined by the dust-gas coupling strength and the dust-to-gas mass ratio. The high and intermediate wavenumber regimes behave similarly to one another, with large dust-to-gas ratio fluctuations while the gas remains largely incompressible. The saturated state is highly anisotropic: dust is concentrated in filaments, jets, or plumes along the direction of acceleration, with turbulent vortex-like structures rapidly forming and dissipating in the perpendicular directions. The low-wavenumber regime exhibits large fluctuations in gas and dust density, but the dust and gas remain more strongly coupled in coherent 'fronts' perpendicular to the acceleration. These behaviors are qualitatively different from those of dust 'passively' driven by external hydrodynamic turbulence, with no back-reaction force from dust onto gas. The virulent nature of these instabilities has interesting implications for dust-driven winds in a variety of astrophysical systems, including around cool-stars, in dusty torii around active-galactic-nuclei, and in and around giant molecular clouds.

Recently Squire & Hopkins showed that charged dust grains moving through magnetized gas under the influence of any external force (e.g. radiation pressure, gravity) are subject to a spectrum of instabilities. Qualitatively distinct instability families are associated with different Alfvenic or magnetosonic waves and drift or gyro motion. We present a suite of simulations exploring these instabilities, for grains in a homogeneous medium subject to an external acceleration. We vary parameters such as the ratio of Lorentz-to-drag forces on dust, plasma $\beta$, size scale, and acceleration. All regimes studied drive turbulent motions and dust-to-gas fluctuations in the saturated state, can rapidly amplify magnetic fields into equipartition with velocity fluctuations, and produce instabilities that persist indefinitely (despite random grain motions). Different parameters produce diverse morphologies and qualitatively different features in dust, but the saturated gas state can be broadly characterized as anisotropic magnetosonic or Alfvenic turbulence. Quasi-linear theory can qualitatively predict the gas turbulent properties. Turbulence grows from small to large scales, and larger-scale modes usually drive more vigorous gas turbulence, but dust velocity and density fluctuations are more complicated. In many regimes, dust forms structures (clumps, filaments, sheets) that reach extreme over-densities (up to $\gg 10^{9}$ times mean), and exhibit substantial sub-structure even in nearly-incompressible gas. These can be even more prominent at lower dust-to-gas ratios. In other regimes, dust self-excites scattering via magnetic fluctuations that isotropize and amplify dust velocities, producing fast, diffusive dust motions.

We present a systematic shearing-box investigation of MRI-driven turbulence in a weakly collisional plasma by including the effects of an anisotropic pressure stress, i.e. anisotropic (Braginskii) viscosity. We constrain the pressure anisotropy ($\Delta p$) to lie within the stability bounds that would be otherwise imposed by kinetic microinstabilities. We explore a broad region of parameter space by considering different Reynolds numbers and magnetic-field configurations, including net vertical flux, net toroidal-vertical flux and zero net flux. Remarkably, we find that the level of turbulence and angular-momentum transport are not greatly affected by large anisotropic viscosities: the Maxwell and Reynolds stresses do not differ much from the MHD result. Angular-momentum transport in Braginskii MHD still depends strongly on isotropic dissipation, e.g., the isotropic magnetic Prandtl number, even when the anisotropic viscosity is orders of magnitude larger than the isotropic diffusivities. Braginskii viscosity nevertheless changes the flow structure, rearranging the turbulence to largely counter the parallel rate of strain from the background shear. We also show that the volume-averaged pressure anisotropy and anisotropic viscous transport decrease with increasing isotropic Reynolds number (${\rm Re}$); e.g., in simulations with net vertical field, the ratio of anisotropic to Maxwell stress ($\alpha_{\rm A} / \alpha_{\rm M}$) decreases from $\sim 0.5$ to $\sim 0.1$ as we move from ${\rm Re} \sim 10^3$ to ${\rm Re} \sim 10^4$. Anisotropic transport may thus become negligible at high ${\rm Re}$. Anisotropic viscosity nevertheless becomes the dominant source of heating at large ${\rm Re}$, accounting for $\gtrsim 50 \%$ of the plasma heating. We conclude by briefly discussing the implications of our results for RIAFs onto black holes.

This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed.

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

We propose that pressure anisotropy causes weakly collisional turbulent plasmas to self-organize so as to resist changes in magnetic-field strength. We term this effect "magneto-immutability" by analogy with incompressibility (resistance to changes in pressure). The effect is important when the pressure anisotropy becomes comparable to the magnetic pressure, suggesting that in collisionless, weakly magnetized (high-$\beta$) plasmas its dynamical relevance is similar to that of incompressibility. Simulations of magnetized turbulence using the weakly collisional Braginskii model show that magneto-immutable turbulence is surprisingly similar, in most statistical measures, to critically balanced MHD turbulence. However, in order to minimize magnetic-field variation, the flow direction becomes more constrained than in MHD, and the turbulence is more strongly dominated by magnetic energy (a nonzero "residual energy"). These effects represent key differences between pressure-anisotropic and fluid turbulence, and should be observable in the $\beta\gtrsim1$ turbulent solar wind.

Recently, Squire & Hopkins (2017) showed any coupled dust-gas mixture is subject to a class of linear 'resonant drag instabilities' (RDI). These can drive large dust-to-gas ratio fluctuations even at arbitrarily small dust-to-gas mass ratios. Here, we identify and study both resonant and new non-resonant instabilities, in the simple case where the gas satisfies neutral hydrodynamics and supports acoustic waves ($\omega^{2}=c_{s}^{2}\,k^{2}$). The gas and dust are coupled via an arbitrary drag law and subject to external accelerations (e.g. gravity, radiation pressure). If there is any dust drift velocity, the system is unstable. The instabilities exist for {\em all} dust-to-gas ratios $\mu$ and their growth rates depend only weakly on $\mu$ around resonance, as $\sim\mu^{1/3}$ or $\sim \mu^{1/2}$ (depending on wavenumber). The behavior changes depending on whether the drift velocity is larger or smaller than the sound speed $c_{s}$. In the supersonic regime a 'resonant' instability appears with growth rate increasing without limit with wavenumber, even for vanishingly small $\mu$ and values of the coupling strength ('stopping time'). In the subsonic regime non-resonant instabilities always exist, but their growth rates no longer increase indefinitely towards small wavelengths. The dimensional scalings and qualitative behavior of the instability do not depend sensitively on the drag law or equation-of-state of the gas. The instabilities directly drive exponentially growing dust-to-gas-ratio fluctuations, which can be large even when the modes are otherwise weak. We discuss physical implications for cool-star winds, AGN-driven winds and torii, and starburst winds: the instabilities alter the character of these outflows and could drive clumping and/or turbulence in the dust and gas.

Squire & Hopkins (2017) showed that coupled dust-gas mixtures are generically subject to 'resonant drag instabilities' (RDIs), which drive violently-growing fluctuations in both. But the role of magnetic fields and charged dust has not yet been studied. We therefore explore the RDI in gas which obeys ideal MHD and is coupled to dust via both Lorentz forces and drag, with an external acceleration (e.g., gravity, radiation) driving dust drift through gas. We show this is always unstable, at all wavelengths and non-zero values of dust-to-gas ratio, drift velocity, dust charge, 'stopping time' or drag coefficient (for any drag law), or field strength; moreover growth rates depend only weakly (sub-linearly) on these parameters. Dust charge and magnetic fields do not suppress instabilities, but give rise to a large number of new instability 'families,' each with distinct behavior. The 'MHD-wave' (magnetosonic or Alfven) RDIs exhibit maximal growth along 'resonant' angles where the modes have a phase velocity matching the corresponding MHD wave, and growth rates increase without limit with wavenumber. The 'gyro' RDIs are driven by resonances between drift and Larmor frequencies, giving growth rates sharply peaked at specific wavelengths. Other instabilities include 'acoustic' and 'pressure-free' modes (previously studied), and a family akin to cosmic ray instabilities which appear when Lorentz forces are strong and dust streams super-Alfvenically along field lines. We discuss astrophysical applications in the warm ISM, CGM/IGM, HII regions, SNe ejecta/remnants, Solar corona, cool-star winds, GMCs, and AGN.

We investigate, for the first time, the nonlinear evolution of the magnetized "resonant drag instabilities" (RDIs). We explore magnetohydrodynamic (MHD) simulations of gas mixed with (uniform) dust grains subject to Lorentz and drag forces, using the GIZMO code. The magnetized RDIs exhibit fundamentally different behaviour than the purely acoustic RDIs. The dust organizes into coherent structures and the system exhibits strong dust-gas separation. In the linear and early nonlinear regime, the growth rates agree with linear theory and the dust self-organizes into two-dimensional planes or "sheets." Eventually the gas develops fully nonlinear, saturated Alfv\'enic and compressible fast-mode turbulence, which fills the under-dense regions with a small amount of dust, and drives a dynamo which saturates at equipartition of kinetic and magnetic energy. The dust density fluctuations exhibit significant non-Gaussianity, and the power spectrum is strongly weighted towards the largest (box-scale) modes. The saturation level can be understood via quasi-linear theory, as the forcing and energy input via the instabilities becomes comparable to saturated tension forces and dissipation in turbulence. The magnetized simulation presented here is just one case; it is likely that the magnetic RDIs can take many forms in different parts of parameter space.

We identify and study a number of new, rapidly growing instabilities of dust grains in protoplanetary disks, which may be important for planetesimal formation. The study is based on the recognition that dust-gas mixtures are generically unstable to a Resonant Drag Instability (RDI), whenever the gas, absent dust, supports undamped linear modes. We show that the "streaming instability" is an RDI associated with epicyclic oscillations; this provides simple interpretations for its mechanisms and accurate analytic expressions for its growth rates and fastest-growing wavelengths. We extend this analysis to more general dust streaming motions and other waves, including buoyancy and magnetohydrodynamic oscillations, finding various new instabilities. Most importantly, we identify the disk "settling instability," which occurs as dust settles vertically into the midplane of a rotating disk. For small grains, this instability grows many orders of magnitude faster than the standard streaming instability, with a growth rate that is independent of grain size. Growth timescales for realistic dust-to-gas ratios are comparable to the disk orbital period, and the characteristic wavelengths are more than an order of magnitude larger than the streaming instability (allowing the instability to concentrate larger masses). This suggests that in the process of settling, dust will band into rings then filaments or clumps, potentially seeding dust traps, high-metallicity regions that in turn seed the streaming instability, or even overdensities that coagulate or directly collapse to planetesimals.

We show that grains streaming through a fluid are generically unstable if their velocity, projected along some direction, matches the phase velocity of a fluid wave (linear oscillation). This can occur whenever grains stream faster than any fluid wave. The wave itself can be quite general--sound waves, magnetosonic waves, epicyclic oscillations, and Brunt-V\"ais\"al\"a oscillations each generate instabilities, for example. We derive a simple expression for the growth rates of these "resonant drag instabilities" (RDI). This expression (i) illustrates why such instabilities are so virulent and generic, and (ii) allows for simple analytic computation of RDI growth rates and properties for different fluids. As examples, we introduce several new instabilities, which could see application across a variety of physical systems from atmospheres to protoplanetary disks, the interstellar medium, and galactic outflows. The matrix-based resonance formalism we introduce can also be applied more generally in other (nonfluid) contexts, providing a simple means for calculating and understanding the stability properties of interacting systems.

Using two-dimensional hybrid-kinetic simulations, we explore the nonlinear "interruption" of standing and traveling shear-Alfv\'en waves in collisionless plasmas. Interruption involves a self-generated pressure anisotropy removing the restoring force of a linearly polarized Alfv\'enic perturbation, and occurs for wave amplitudes $\delta B_{\perp}/B_{0}\gtrsim \beta^{\,-1/2}$ (where $\beta$ is the ratio of thermal to magnetic pressure). We use highly elongated domains to obtain maximal scale separation between the wave and the ion gyroscale. For standing waves above the amplitude limit, we find that the large-scale magnetic field of the wave decays rapidly. The dynamics are strongly affected by the excitation of oblique firehose modes, which transition into long-lived parallel fluctuations at the ion gyroscale and cause significant particle scattering. Traveling waves are damped more slowly, but are also influenced by small-scale parallel fluctuations created by the decay of firehose modes. Our results demonstrate that collisionless plasmas cannot support linearly polarized Alfv\'en waves above $\delta B_{\perp}/B_{0}\sim \beta^{\,-1/2}$. They also provide a vivid illustration of two key aspects of low-collisionality plasma dynamics: (i) the importance of velocity-space instabilities in regulating plasma dynamics at high $\beta$, and (ii) how nonlinear collisionless processes can transfer mechanical energy directly from the largest scales into thermal energy and microscale fluctuations, without the need for a scale-by-scale turbulent cascade.

In collisionless and weakly collisional plasmas, such as hot accretion flows onto compact objects, the magnetorotational instability (MRI) can differ significantly from the standard (collisional) MRI. In particular, pressure anisotropy with respect to the local magnetic-field direction can both change the linear MRI dispersion relation and cause nonlinear modifications to the mode structure and growth rate, even when the field and flow perturbations are small. This work studies these pressure-anisotropy-induced nonlinearities in the weakly nonlinear, high-ion-beta regime, before the MRI saturates into strong turbulence. Our goal is to better understand how the saturation of the MRI in a low collisionality plasma might differ from that in the collisional regime. We focus on two key effects: (i) the direct impact of self-induced pressure-anisotropy nonlinearities on the evolution of an MRI mode, and (ii) the influence of pressure anisotropy on the "parasitic instabilities" that are suspected to cause the mode to break up into turbulence. Our main conclusions are: (i) The mirror instability regulates the pressure anisotropy in such a way that the linear MRI in a collisionless plasma is an approximate nonlinear solution once the mode amplitude becomes larger than the background field (just as in MHD). This implies that differences between the collisionless and collisional MRI become unimportant at large amplitudes. (ii) The break up of large amplitude MRI modes into turbulence via parasitic instabilities is similar in collisionless and collisional plasmas. Together, these conclusions suggest that the route to magnetorotational turbulence in a collisionless plasma may well be similar to that in a collisional plasma, as suggested by recent kinetic simulations. As a supplement to these findings, we offer guidance for the design of future kinetic simulations of magnetorotational turbulence.

We study the dynamics of large, charged dust grains in turbulent giant molecular clouds (GMCs). Massive dust grains behave as aerodynamic particles in primarily neutral dense gas, and thus are able to produce dramatic small-scale fluctuations in the dust-to-gas ratio. Hopkins & Lee (2016) directly simulated the dynamics of neutral dust grains in super-sonic MHD turbulence, typical of GMCs, and showed that the dust-to-gas fluctuations can exceed factor ~1000 on small scales, with important implications for star formation, stellar abundances, and dust behavior and growth. However, even in primarily neutral gas in GMCs, dust grains are negatively charged and Lorentz forces are non-negligible. Therefore, we extend our previous study by including the effects of Lorentz forces on charged grains (in addition to drag). For small charged grains (sizes <0.1 micron), Lorentz forces suppress dust-to-gas ratio fluctuations, while for large grains (sizes >1 micron), Lorentz forces have essentially no effect, trends that are well explained with a simple theory of dust magnetization. In some special intermediate cases, Lorentz forces can enhance dust-gas segregation. Regardless, for the physically expected scaling of dust charge with grain size, we find the most important effects depend on grain size (via the drag equation) with Lorentz forces/charge as a second-order correction. We show that the dynamics we consider are determined by three dimensionless numbers in the limit of weak background magnetic fields: the turbulent Mach number, a dust drag parameter (proportional to grain size) and a dust Lorentz parameter (proportional to grain charge); these allow us to generalize our simulations to a wide range of conditions.

We propose a model for the density statistics in supersonic turbulence, which play a crucial role in star-formation and the physics of the interstellar medium (ISM). Motivated by [Hopkins, MNRAS, 430, 1880 (2013)], the model considers the density to be arranged into a collection of strong shocks of width $\sim\! \mathcal{M}^{-2}$, where $\mathcal{M}$ is the turbulent Mach number. With two physically motivated parameters, the model predicts all density statistics for $\mathcal{M}>1$ turbulence: the density probability distribution and its intermittency (deviation from log-normality), the density variance-Mach number relation, power spectra, and structure functions. For the proposed model parameters, reasonable agreement is seen between model predictions and numerical simulations, albeit within the large uncertainties associated with current simulation results. More generally, the model could provide a useful framework for more detailed analysis of future simulations and observational data. Due to the simple physical motivations for the model in terms of shocks, it is straightforward to generalize to more complex physical processes, which will be helpful in future more detailed applications to the ISM. We see good qualitative agreement between such extensions and recent simulations of non-isothermal turbulence.

Passive scalar mixing (metals, molecules, etc.) in the turbulent interstellar medium (ISM) is critical for abundance patterns of stars and clusters, galaxy and star formation, and cooling from the circumgalactic medium. However, the fundamental scaling laws remain poorly understood in the highly supersonic, magnetized, shearing regime relevant for the ISM. We therefore study the full scaling laws governing passive-scalar transport in idealized simulations of supersonic turbulence. Using simple phenomenological arguments for the variation of diffusivity with scale based on Richardson diffusion, we propose a simple fractional diffusion equation to describe the turbulent advection of an initial passive scalar distribution. These predictions agree well with the measurements from simulations, and vary with turbulent Mach number in the expected manner, remaining valid even in the presence of a large-scale shear flow (e.g. rotation in a galactic disk). The evolution of the scalar distribution is not the same as obtained using simple, constant "effective diffusivity" as in Smagorinsky models, because the scale-dependence of turbulent transport means an initially Gaussian distribution quickly develops highly non-Gaussian tails. We also emphasize that these are mean scalings that only apply to ensemble behaviors (assuming many different, random scalar injection sites): individual Lagrangian "patches" remain coherent (poorly-mixed) and simply advect for a large number of turbulent flow-crossing times.

This work, which extends Squire et al. [ApJL, 830 L25 (2016)], explores the effect of self-generated pressure anisotropy on linearly polarized shear-Alfv\'en fluctuations in low-collisionality plasmas. Such anisotropies lead to stringent limits on the amplitude of magnetic perturbations in high-beta plasmas, above which a fluctuation can destabilize itself through the parallel firehose instability. This causes the wave frequency to approach zero, "interrupting" the wave and stopping its oscillation. These effects are explored in detail in the collisionless and weakly collisional "Braginskii" regime, for both standing and traveling waves. The focus is on simplified models in one dimension, on scales much larger than the ion gyroradius. The effect has interesting implications for the physics of magnetized turbulence in the high-beta conditions that are prevalent in many astrophysical plasmas.

It is shown that low-collisionality plasmas cannot support linearly polarized shear-Alfv\'en fluctuations above a critical amplitude $\delta B_{\perp}/B_{0} \sim \beta^{\,-1/2}$, where $\beta$ is the ratio of thermal to magnetic pressure. Above this cutoff, a developing fluctuation will generate a pressure anisotropy that is sufficient to destabilize itself through the parallel firehose instability. This causes the wave frequency to approach zero, interrupting the fluctuation before any oscillation. The magnetic field lines rapidly relax into a sequence of angular zig-zag structures. Such a restrictive bound on shear-Alfv\'en-wave amplitudes has far-reaching implications for the physics of magnetized turbulence in the high-$\beta$ conditions prevalent in many astrophysical plasmas, as well as for the solar wind at $\sim 1 \mathrm{AU}$ where $\beta \gtrsim 1$.

A novel large-scale dynamo mechanism, the magnetic shear-current effect, is discussed and explored. The effect relies on the interaction of magnetic fluctuations with a mean shear flow, meaning the saturated state of the small-scale dynamo can drive a large-scale dynamo -- in some sense the inverse of dynamo quenching. The dynamo is nonhelical, with the mean-field $\alpha$ coefficient zero, and is caused by the interaction between an off-diagonal component of the turbulent resistivity and the stretching of the large-scale field by shear flow. Following up on previous numerical and analytic work, this paper presents further details of the numerical evidence for the effect, as well as an heuristic description of how magnetic fluctuations can interact with shear flow to produce the required electromotive force. The pressure response of the fluid is fundamental to this mechanism, which helps explain why the magnetic effect is stronger than its kinematic cousin, and the basic idea is related to the well-known lack of turbulent resistivity quenching by magnetic fluctuations. As well as being interesting for its applications to general high Reynolds number astrophysical turbulence, where strong small-scale magnetic fluctuations are expected to be prevalent, the magnetic shear-current effect is a likely candidate for large-scale dynamo in the unstratified regions of ionized accretion disks. Evidence for this is discussed, as well as future research directions and the challenges involved with understanding details of the effect in astrophysically relevant regimes.

We propose a new mechanism for turbulent mean-field dynamo in which the magnetic fluctuations resulting from a small-scale dynamo drive the generation of large-scale magnetic fields. This is in stark contrast to the common idea that small-scale magnetic fields should be harmful to large-scale dynamo action. These dynamos occur in the presence of large-scale velocity shear and do not require net helicity, resulting from off-diagonal components of the turbulent resistivity tensor as the magnetic analogue of the "shear-current" effect. Given the inevitable existence of non-helical small-scale magnetic fields in turbulent plasmas, as well as the generic nature of velocity shear, the suggested mechanism may help to explain generation of large-scale magnetic fields across a wide range of astrophysical objects.

This article presents a calculation of the mean electromotive force arising from general small-scale magnetohydrodynamical turbulence, within the framework of the second-order correlation approximation. With the goal of improving understanding of the accretion disk dynamo, effects arising through small-scale magnetic fluctuations, velocity gradients, density and turbulence stratification, and rotation, are included. The primary result, which supplements numerical findings, is that an off-diagonal turbulent resistivity due to magnetic fluctuations can produce large-scale dynamo action -- the magnetic analogue of the "shear-current" effect. In addition, consideration of $\alpha$ effects in the stratified regions of disks gives the puzzling result that there is no strong prediction for a sign of $\alpha$, since the effects due to kinetic and magnetic fluctuations, as well as those due to shear and rotation, are each of opposing signs and tend to cancel each other.

Nonhelical shear dynamos are studied with a particular focus on the possibility of coherent dynamo action. The primary results -- serving as a follow up to the results of Squire & Bhattacharjee [arXiv:1506.04109 (2015)] -- pertain to the "magnetic shear-current effect" as a viable mechanism to drive large-scale magnetic field generation. This effect raises the interesting possibility that the saturated state of the small-scale dynamo could drive large-scale dynamo action, and is likely to be important in the unstratified regions of accretion disk turbulence. In this paper, the effect is studied at low Reynolds numbers, removing the complications of small-scale dynamo excitation and aiding analysis by enabling the use of quasi-linear statistical simulation methods. In addition to the magnetically driven dynamo, new results on the kinematic nonhelical shear dynamo are presented. These illustrate the relationship between coherent and incoherent driving in such dynamos, demonstrating the importance of rotation in determining the relative dominance of each mechanism.

Turbulence and dynamo induced by the magnetorotational instability (MRI) are analyzed using quasi-linear statistical simulation methods. It is found that homogenous turbulence is unstable to a large scale dynamo instability, which saturates to an inhomogenous equilibrium with a strong dependence on the magnetic Prandtl number (Pm). Despite its enormously reduced nonlinearity, the dependence of the angular momentum transport on Pm in the quasi-linear model is qualitatively similar to that of nonlinear MRI turbulence. This indicates that recent convergence problems may be related to large scale dynamo and suggests how dramatically simplified models may be used to gain insight into the astrophysically relevant regimes of very low or high Pm.

We study the magnetorotational instability (MRI) using nonmodal stability techniques. Despite the spectral instability of many forms of the MRI, this proves to be a natural method of analysis that is well-suited to deal with the non-self-adjoint nature of the linear MRI equations. We find that the fastest growing linear MRI structures on both local and global domains can look very different to the eigenmodes, invariably resembling waves shearing with the background flow (shear waves). In addition, such structures can grow many times faster than the least stable eigenmode over long time periods, and be localized in a completely different region of space. These ideas lead -- for both axisymmetric and non-axisymmetric modes -- to a natural connection between the global MRI and the local shearing box approximation. By illustrating that the fastest growing global structure is well described by the ordinary differential equations (ODEs) governing a single shear wave, we find that the shearing box is a very sensible approximation for the linear MRI, contrary to many previous claims. Since the shear wave ODEs are most naturally understood using nonmodal analysis techniques, we conclude by analyzing local MRI growth over finite time-scales using these methods. The strong growth over a wide range of wave-numbers suggests that nonmodal linear physics could be of fundamental importance in MRI turbulence.

We analyze the linear growth of the magnetorotational instability (MRI) in the short time limit using nonmodal methods. Our findings are quite different to standard results, illustrating that shearing wave energy can grow at the maximum MRI rate, $-d\Omega/d \ln r,$ for any choice of azimuthal and vertical wavelengths. In addition, by comparing the growth of shearing waves with static structures, we show that over short time-scales shearing waves will always be dynamically more important than static structures in the ideal limit. By demonstrating that fast linear growth is possible at all wavelengths, these results suggest that nonmodal linear physics could play a fundamental role in MRI turbulence.