Showing 1–50 of 68 results for author: Squire, J

Extended Cyclotron Resonant Heating of the Turbulent Solar Wind

Authors: Trevor A. Bowen, Ivan Y. Vasko, Stuart D. Bale, Benjamin D. G. Chandran, Alexandros Chasapis, Thierry Dudok de Wit, Alfred Mallet, Michael McManus, Romain Meyrand, Marc Pulupa, Jonathan Squire

Abstract: Circularly polarized, nearly parallel propagating waves are prevalent in the solar wind at ion-kinetic scales. At these scales, the spectrum of turbulent fluctuations in the solar wind steepens, often called the transition-range, before flattening at sub-ion scales. Circularly polarized waves have been proposed as a mechanism to couple electromagnetic fluctuations to ion gyromotion, enabling ion-s… ▽ More Circularly polarized, nearly parallel propagating waves are prevalent in the solar wind at ion-kinetic scales. At these scales, the spectrum of turbulent fluctuations in the solar wind steepens, often called the transition-range, before flattening at sub-ion scales. Circularly polarized waves have been proposed as a mechanism to couple electromagnetic fluctuations to ion gyromotion, enabling ion-scale dissipation that results in observed ion-scale steepening. Here, we study Parker Solar Probe observations of an extended stream of fast solar wind ranging from 15-55 solar radii. We demonstrate that, throughout the stream, transition-range steepening at ion-scales is associated with the presence of significant left handed ion-kinetic scale waves, which are thought to be ion-cyclotron waves. We implement quasilinear theory to compute the rate at which ions are heated via cyclotron resonance with the observed circularly polarized waves given the empirically measured proton velocity distribution functions. We apply the Von Karman decay law to estimate the turbulent decay of the large-scale fluctuations, which is equal to the turbulent energy cascade rate. We find that the ion-cyclotron heating rates are correlated with, and amount to a significant fraction of, the turbulent energy cascade rate, implying that cyclotron heating is an important dissipation mechanism in the solar wind. △ Less

Submitted 14 June, 2024; originally announced June 2024.

Microphysical Regulation of Non-Ideal MHD in Weakly-Ionized Systems: Does the Hall Effect Matter?

Authors: Philip F. Hopkins, Jonathan Squire, Raphael Skalidis, Nadine H. Soliman

Abstract: The magnetohydrodynamics (MHD) equations plus 'non-ideal' (Ohmic, Hall, ambipolar) resistivities are widely used to model weakly-ionized astrophysical systems. We show that if gradients in the magnetic field become too steep, the implied charge drift speeds become much faster than microphysical signal speeds, invalidating the assumptions used to derive both the resistivities and MHD equations them… ▽ More The magnetohydrodynamics (MHD) equations plus 'non-ideal' (Ohmic, Hall, ambipolar) resistivities are widely used to model weakly-ionized astrophysical systems. We show that if gradients in the magnetic field become too steep, the implied charge drift speeds become much faster than microphysical signal speeds, invalidating the assumptions used to derive both the resistivities and MHD equations themselves. Generically this situation will excite microscale instabilities that suppress the drift and current. We show this could be relevant at low ionization fractions especially if Hall terms appear significant, external forces induce supersonic motions, or dust grains become a dominant charge carrier. Considering well-established treatments of super-thermal drifts in laboratory, terrestrial, and Solar plasmas as well as conduction and viscosity models, we generalize a simple prescription to rectify these issues, where the resistivities are multiplied by a correction factor that depends only on already-known macroscopic quantities. This is generalized for multi-species and weakly-ionized systems, and leaves the equations unchanged in the drift limits for which they are derived, but restores physical behavior (driving the system back towards slow drift by diffusing away small-scale gradients in the magnetic field) if the limits are violated. This has important consequences: restoring intuitive behaviors such as the system becoming hydrodynamic in the limit of zero ionization; suppressing magnetic structure on scales below a critical length which can comparable to circumstellar disk sizes; limiting the maximum magnetic amplification; and suppressing the effects of the Hall term in particular. This likely implies that the Hall term does not become dynamically important under most conditions of interest in these systems. △ Less

Submitted 9 May, 2024; originally announced May 2024.

Comments: 16 pages, 4 figures. Submitted to the Open Journal of Astrophysics. Comments welcome

Self-organization in collisionless, high-$β$ turbulence

Authors: S. Majeski, M. W. Kunz, J. Squire

Abstract: The MHD equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only po… ▽ More The MHD equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, $β\gg 1$ turbulence, focusing on the role of a self-organizational effect known as `magneto-immutability'. We describe this self-organization analytically through a high-$β$, reduced ordering of the CGL-MHD equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses, and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small ($\sim 5\%$) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-versus-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 42 pages, 13 figures

The effects of finite electron inertia on helicity-barrier-mediated turbulence

Authors: T. Adkins, R. Meyrand, J. Squire

Abstract: Understanding the partitioning of turbulent energy between ions and electrons in weakly collisional plasmas is crucial for the accurate interpretation of observations and modelling of various astrophysical phenomena. Many such plasmas are "imbalanced", wherein the large-scale energy input is dominated by Alfvénic fluctuations propagating in a single direction. In this paper, we demonstrate that wh… ▽ More Understanding the partitioning of turbulent energy between ions and electrons in weakly collisional plasmas is crucial for the accurate interpretation of observations and modelling of various astrophysical phenomena. Many such plasmas are "imbalanced", wherein the large-scale energy input is dominated by Alfvénic fluctuations propagating in a single direction. In this paper, we demonstrate that when strongly-magnetised plasma turbulence is imbalanced, nonlinear conservation laws imply the existence of a critical value of the electron plasma beta that separates two dramatically different types of turbulence in parameter space. For betas below the critical value, the free energy injected on the largest scales is able to undergo a familiar Kolmogorov-type cascade to small scales where it is dissipated, heating electrons. For betas above the critical value, the system forms a "helicity barrier" that prevents the cascade from proceeding past the ion Larmor radius, causing the majority of the injected free energy to be deposited into ion heating. Physically, the helicity barrier results from the inability of the system to adjust to the disparity between the perpendicular-wavenumber scalings of the free energy and generalised helicity below the ion Larmor radius; restoring finite electron inertia can annul, or even reverse, this disparity, giving rise to the aforementioned critical beta. We relate this physics to the "dynamic phase alignment" mechanism, and characterise various other important features of the helicity barrier, including the nature of the nonlinear wavenumber-space fluxes, dissipation rates, and energy spectra. The existence of such a critical beta has important implications for heating, as it suggests that the dominant recipient of the turbulent energy -- ions or electrons -- can depend sensitively on the characteristics of the plasma at large scales. △ Less

Submitted 29 May, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: 33 pages, 10 figures

An Analytic Model For Magnetically-Dominated Accretion Disks

Authors: Philip F. Hopkins, Jonathan Squire, Eliot Quataert, Norman Murray, Kung-Yi Su, Ulrich P. Steinwandel, Kyle Kremer, Claude-Andre Faucher-Giguere, Sarah Wellons

Abstract: Recent numerical cosmological radiation-magnetohydrodynamic-thermochemical-star formation simulations have resolved the formation of quasar accretion disks with Eddington or super-Eddington accretion rates onto supermassive black holes (SMBHs) down to a few hundred gravitational radii. These 'flux-frozen' and hyper-magnetized disks appear to be qualitatively distinct from classical $α$ disks and m… ▽ More Recent numerical cosmological radiation-magnetohydrodynamic-thermochemical-star formation simulations have resolved the formation of quasar accretion disks with Eddington or super-Eddington accretion rates onto supermassive black holes (SMBHs) down to a few hundred gravitational radii. These 'flux-frozen' and hyper-magnetized disks appear to be qualitatively distinct from classical $α$ disks and magnetically-arrested disks: the midplane pressure is dominated by toroidal magnetic fields with plasma $β\ll 1$ powered by advection of magnetic flux from the interstellar medium (ISM), and they are super-sonically and trans-Alfvenically turbulent with cooling times short compared to dynamical times yet remain gravitationally stable owing to magnetic support. In this paper, we present a simple analytic similarity model for such disks. For reasonable assumptions, the model is entirely specified by the boundary conditions (inflow rate at the BH radius of influence [BHROI]). We show that the scalings from this model are robust to various detailed assumptions, agree remarkably well with the simulations (given their simplicity), and demonstrate the self-consistency and gravitational stability of such disks even in the outer accretion disk (approaching the BHROI) at hyper-Eddington accretion rates. △ Less

Submitted 12 March, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 12 pages, 1 figure. Replaced with accepted version. Companion paper to 'FORGE'D IN FIRE II' (arXiv:2310.04506, part of a series with arXiv:2309.13115 -- animations of the simulations referred to here can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html)

FORGE'd in FIRE II: The Formation of Magnetically-Dominated Quasar Accretion Disks from Cosmological Initial Conditions

Authors: Philip F. Hopkins, Jonathan Squire, Kung-Yi Su, Ulrich P. Steinwandel, Kyle Kremer, Yanlong Shi, Michael Y. Grudic, Sarah Wellons, Claude-Andre Faucher-Giguere, Daniel Angles-Alcazar, Norman Murray, Eliot Quataert

Abstract: In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc s… ▽ More In a companion paper, we reported the self-consistent formation of quasar accretion disks with inflow rates $\sim 10\,{\rm M_{\odot}\,yr^{-1}}$ down to <300 Schwarzschild radii from cosmological radiation-magneto-thermochemical-hydrodynamical galaxy and star formation simulations. We see the formation of a well-defined, steady-state accretion disk which is stable against star formation at sub-pc scales. The disks are optically thick, with radiative cooling balancing accretion, but with properties that are distinct from those assumed in most previous accretion disk models. The pressure is strongly dominated by (primarily toroidal) magnetic fields, with a plasma $β\sim 10^{-4}$ even in the disk midplane. They are qualitatively distinct from magnetically elevated or arrested disks. The disks are strongly turbulent, with trans-Alfvenic and highly super-sonic turbulence, and balance this via a cooling time that is short compared to the disk dynamical time, and can sustain highly super-Eddington accretion rates. Their surface and 3D densities at $\sim 10^{3}-10^{5}$ gravitational radii are much lower than in a Shakura-Sunyaev disk, with important implications for their thermo-chemistry and stability. We show how the magnetic field strengths and geometries arise from rapid advection of flux with the inflow from much weaker galaxy-scale fields in these 'flux-frozen' disks, and how this stabilizes the disk and gives rise to efficient torques. Re-simulating without magnetic fields produces catastrophic fragmentation with a vastly smaller, lower-$\dot{M}$ Shakura-Sunyaev-like disk. △ Less

Submitted 18 January, 2024; v1 submitted 6 October, 2023; originally announced October 2023.

Comments: 46 pages, 30 figures. Updated to match accepted version from The Open Journal of Astrophysics. Part of a series with arXiv:2309.13115 . Animations of the simulations can be viewed at http://www.tapir.caltech.edu/~phopkins/Site/animations/Movies_zoom.html

Electron-ion heating partition in imbalanced solar-wind turbulence

Authors: Jonathan Squire, Romain Meyrand, Matthew W. Kunz

Abstract: A likely candidate mechanism to heat the solar corona and solar wind is low-frequency "Alfvénic" turbulence sourced by magnetic fluctuations near the solar surface. Depending on its properties, such turbulence can heat different species via different mechanisms, and the comparison of theoretical predictions to observed temperatures, wind speeds, anisotropies, and their variation with heliocentric… ▽ More A likely candidate mechanism to heat the solar corona and solar wind is low-frequency "Alfvénic" turbulence sourced by magnetic fluctuations near the solar surface. Depending on its properties, such turbulence can heat different species via different mechanisms, and the comparison of theoretical predictions to observed temperatures, wind speeds, anisotropies, and their variation with heliocentric radius provides a sensitive test of this physics. Here we explore the importance of normalized cross helicity, or imbalance, for controlling solar-wind heating, since it is a key parameter of magnetized turbulence and varies systematically with wind speed and radius. Based on a hybrid-kinetic simulation in which the forcing's imbalance decreases with time -- a crude model for a plasma parcel entrained in the outflowing wind -- we demonstrate how significant changes to the turbulence and heating result from the "helicity barrier" effect. Its dissolution at low imbalance causes its characteristic features -- strong perpendicular ion heating with a steep "transition-range" drop in electromagnetic fluctuation spectra -- to disappear, driving more energy into electrons and parallel ion heat, and halting the emission of ion-scale waves. These predictions seem to agree with a diverse array of solar-wind observations, offering to explain a variety of complex correlations and features within a single theoretical framework. △ Less

Submitted 8 October, 2023; v1 submitted 24 August, 2023; originally announced August 2023.

Reflection-driven turbulence in the super-Alfvénic solar wind

Authors: Romain Meyrand, Jonathan Squire, Alfred Mallet, Benjamin D. G. Chandran

Abstract: In magnetized, stratified astrophysical environments such as the Sun's corona and solar wind, Alfvénic fluctuations ''reflect'' from background gradients, enabling nonlinear interactions and thus dissipation of their energy into heat. This process, termed ''reflection-driven turbulence,'' is thought to play a crucial role in coronal heating and solar-wind acceleration, explaining a range of observ… ▽ More In magnetized, stratified astrophysical environments such as the Sun's corona and solar wind, Alfvénic fluctuations ''reflect'' from background gradients, enabling nonlinear interactions and thus dissipation of their energy into heat. This process, termed ''reflection-driven turbulence,'' is thought to play a crucial role in coronal heating and solar-wind acceleration, explaining a range of observational correlations and constraints. Building on previous works focused on the inner heliosphere, here we study the basic physics of reflection-driven turbulence using reduced magnetohydrodynamics in an expanding box -- the simplest model that can capture the local turbulent plasma dynamics in the super-Alfvénic solar wind. Although idealized, our high-resolution simulations and simple theory reveal a rich phenomenology that is consistent with a diverse range of observations. Outwards-propagating fluctuations, which initially have high imbalance, decay nonlinearly to heat the plasma, becoming more balanced and magnetically dominated. Despite the high imbalance, the turbulence is strong because Elsässer collisions are suppressed by reflection, leading to ''anomalous coherence'' between the two Elsässer fields. This coherence, together with linear effects, causes the turbulence to anomalously grow the ''anastrophy'' (squared magnetic potential) as it decays, forcing the energy to rush to larger scales and forming a ''$1/f$-range'' energy spectrum as it does so. At late times, the expansion overcomes the nonlinear and Alfvénic physics, forming isolated, magnetically dominated ''Alfvén vortex'' structures that minimize their nonlinear dissipation. These results can plausibly explain the observed radial and wind-speed dependence of turbulence imbalance, residual energy, plasma heating, and fluctuation spectra, as well as making testable predictions for future observations. △ Less

Submitted 20 August, 2023; originally announced August 2023.

Galactic Cosmic-ray Scattering due to Intermittent Structures

Authors: Iryna S. Butsky, Philip F. Hopkins, Philipp Kempski, Sam B. Ponnada, Eliot Quataert, Jonathan Squire

Abstract: Cosmic rays (CRs) with energies $\ll$ TeV comprise a significant component of the interstellar medium (ISM). Major uncertainties in CR behavior on observable scales (much larger than CR gyroradii) stem from how magnetic fluctuations scatter CRs in pitch angle. Traditional first-principles models, which assume these magnetic fluctuations are weak and uniformly scatter CRs in a homogeneous ISM, stru… ▽ More Cosmic rays (CRs) with energies $\ll$ TeV comprise a significant component of the interstellar medium (ISM). Major uncertainties in CR behavior on observable scales (much larger than CR gyroradii) stem from how magnetic fluctuations scatter CRs in pitch angle. Traditional first-principles models, which assume these magnetic fluctuations are weak and uniformly scatter CRs in a homogeneous ISM, struggle to reproduce basic observables such as the dependence of CR residence times and scattering rates on rigidity. We therefore explore a new category of "patchy" CR scattering models, wherein CRs are predominantly scattered by intermittent strong scattering structures with small volume-filling factors. These models produce the observed rigidity dependence with a simple size distribution constraint, such that larger scattering structures are rarer but can scatter a wider range of CR energies. To reproduce the empirically-inferred CR scattering rates, the mean free path between scattering structures must be $\ell_{\rm mfp} \sim 10$ pc at GeV energies. We derive constraints on the sizes, internal properties, mass/volume-filling factors, and the number density any such structures would need to be both physically and observationally consistent. We consider a range of candidate structures, both large-scale (e.g. H II regions) and small-scale (e.g. intermittent turbulent structures, perhaps even associated with radio plasma scattering) and show that while many macroscopic candidates can be immediately ruled out as the primary CR scattering sites, many smaller structures remain viable and merit further theoretical study. We discuss future observational constraints that could test these models. △ Less

Submitted 23 January, 2024; v1 submitted 11 August, 2023; originally announced August 2023.

Comments: 9 pages, 3 figures, accepted to MNRAS

Mediation of Collisionless Turbulent Dissipation Through Cyclotron Resonance

Authors: Trevor A. Bowen, Stuart D. Bale, Benjamin D. G. Chandran, Alexandros Chasapis, Christopher H. K. Chen, Thierry Dudok de Wit, Alfred Mallet, Romain Meyrand, Jonathan Squire

Abstract: The dissipation of magnetized turbulence is fundamental to understanding energy transfer and heating in astrophysical systems. Collisionless interactions, such as resonant wave-particle process, are known to play a role in shaping turbulent astrophysical environments. Here, we present evidence for the mediation of turbulent dissipation in the solar wind by ion-cyclotron waves. Our results show tha… ▽ More The dissipation of magnetized turbulence is fundamental to understanding energy transfer and heating in astrophysical systems. Collisionless interactions, such as resonant wave-particle process, are known to play a role in shaping turbulent astrophysical environments. Here, we present evidence for the mediation of turbulent dissipation in the solar wind by ion-cyclotron waves. Our results show that ion-cyclotron waves interact strongly with magnetized turbulence, indicating that they serve as a major pathway for the dissipation of large-scale electromagnetic fluctuations. We further show that the presence of cyclotron waves significantly weakens observed signatures of intermittency in sub-ion-kinetic turbulence, which are known to be another pathway for dissipation. These observations results suggest that in the absence of cyclotron resonant waves, non-Gaussian, coherent structures are able to form at sub-ion-kinetic scales, and are likely responsible for turbulent heating. We further find that the cross helicity, i.e. the level of Alfvénicity of the fluctuations, correlates strongly with the presence of ion-scale waves, demonstrating that dissipation of collisionless plasma turbulence is not a universal process, but that the pathways to heating and dissipation at small scales are controlled by the properties of the large-scale turbulent fluctuations. We argue that these observations support the existence of a helicity barrier, in which highly Alfvénic, imbalanced, turbulence is prevented from cascading to sub-ion scales thus resulting in significant ion-cyclotron resonant heating. Our results may serve as a significant step in constraining the nature of turbulent heating in a wide variety of astrophysical systems. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Pressure anisotropy and viscous heating in weakly collisional plasma turbulence

Authors: Jonathan Squire, Matthew W Kunz, Lev Arzamasskiy, Zade Johnston, Eliot Quataert, Alexander A Schekochihin

Abstract: Pressure anisotropy can strongly influence the dynamics of weakly collisional, high-beta plasmas, but its effects are missed by standard magnetohydrodynamics (MHD). Small changes to the magnetic-field strength generate large pressure-anisotropy forces, heating the plasma, driving instabilities, and rearranging flows, even on scales far above the particles' gyroscales where kinetic effects are trad… ▽ More Pressure anisotropy can strongly influence the dynamics of weakly collisional, high-beta plasmas, but its effects are missed by standard magnetohydrodynamics (MHD). Small changes to the magnetic-field strength generate large pressure-anisotropy forces, heating the plasma, driving instabilities, and rearranging flows, even on scales far above the particles' gyroscales where kinetic effects are traditionally considered important. Here, we study the influence of pressure anisotropy on turbulent plasmas threaded by a mean magnetic field (Alfvénic turbulence). Extending previous results that were concerned with Braginskii MHD, we consider a wide range of regimes and parameters using a simplified fluid model based on drift kinetics with heat fluxes calculated using a Landau-fluid closure. We show that viscous (pressure-anisotropy) heating dissipates between a quarter and half of the turbulent cascade power injected at large scales; this does not depend strongly on either plasma beta or the ion-to-electron temperature ratio. This will in turn influence the plasma's thermodynamics by regulating energy partition between different dissipation channels (e.g., electron and ion heat). Due to the pressure anisotropy's rapid dynamical feedback onto the flows that create it -- an effect we term `magneto-immutability' -- the viscous heating is confined to a narrow range of scales near the forcing scale, supporting a nearly conservative, MHD-like inertial-range cascade, via which the rest of the energy is transferred to small scales. Despite the simplified model, our results -- including the viscous heating rate, distributions, and turbulent spectra -- compare favourably to recent hybrid-kinetic simulations. This is promising for the more general use of extended-fluid (or even MHD) approaches to model weakly collisional plasmas such as the intracluster medium, hot accretion flows, and the solar wind. △ Less

Submitted 6 July, 2023; v1 submitted 20 February, 2023; originally announced March 2023.

Comments: Accepted for publication in Journal of Plasma Physics

Parker Solar Probe: Four Years of Discoveries at Solar Cycle Minimum

Authors: N. E. Raouafi, L. Matteini, J. Squire, S. T. Badman, M. Velli, K. G. Klein, C. H. K. Chen, W. H. Matthaeus, A. Szabo, M. Linton, R. C. Allen, J. R. Szalay, R. Bruno, R. B. Decker, M. Akhavan-Tafti, O. V. Agapitov, S. D. Bale, R. Bandyopadhyay, K. Battams, L. Berčič, S. Bourouaine, T. Bowen, C. Cattell, B. D. G. Chandran, R. Chhiber , et al. (32 additional authors not shown)

Abstract: Launched on 12 Aug. 2018, NASA's Parker Solar Probe had completed 13 of its scheduled 24 orbits around the Sun by Nov. 2022. The mission's primary science goal is to determine the structure and dynamics of the Sun's coronal magnetic field, understand how the solar corona and wind are heated and accelerated, and determine what processes accelerate energetic particles. Parker Solar Probe returned a… ▽ More Launched on 12 Aug. 2018, NASA's Parker Solar Probe had completed 13 of its scheduled 24 orbits around the Sun by Nov. 2022. The mission's primary science goal is to determine the structure and dynamics of the Sun's coronal magnetic field, understand how the solar corona and wind are heated and accelerated, and determine what processes accelerate energetic particles. Parker Solar Probe returned a treasure trove of science data that far exceeded quality, significance, and quantity expectations, leading to a significant number of discoveries reported in nearly 700 peer-reviewed publications. The first four years of the 7-year primary mission duration have been mostly during solar minimum conditions with few major solar events. Starting with orbit 8 (i.e., 28 Apr. 2021), Parker flew through the magnetically dominated corona, i.e., sub-Alfvénic solar wind, which is one of the mission's primary objectives. In this paper, we present an overview of the scientific advances made mainly during the first four years of the Parker Solar Probe mission, which go well beyond the three science objectives that are: (1) Trace the flow of energy that heats and accelerates the solar corona and solar wind; (2) Determine the structure and dynamics of the plasma and magnetic fields at the sources of the solar wind; and (3) Explore mechanisms that accelerate and transport energetic particles. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: 157 pages, 65 figures

Journal ref: Space Science Reviews 2023

Microphysically modified magnetosonic modes in collisionless, high-$β$ plasmas

Authors: Stephen Majeski, Matthew W. Kunz, Jonathan Squire

Abstract: With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high-$β$ collisionless plasmas, with particular attention to their excitation of, and reaction to, kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpen… ▽ More With the support of hybrid-kinetic simulations and analytic theory, we describe the nonlinear behaviour of long-wavelength non-propagating (NP) modes and fast magnetosonic waves in high-$β$ collisionless plasmas, with particular attention to their excitation of, and reaction to, kinetic micro-instabilities. The perpendicularly pressure balanced polarization of NP modes produces an excess of perpendicular pressure over parallel pressure in regions where the plasma $β$ is increased. For mode amplitudes $δB/B_0 \gtrsim 0.3$, this excess excites the mirror instability. Particle scattering off these micro-scale mirrors frustrates the nonlinear saturation of transit-time damping, ensuring that large-amplitude NP modes continue their decay to small amplitudes. At asymptotically large wavelengths, we predict that the mirror-induced scattering will be large enough to interrupt transit-time damping entirely, isotropizing the pressure perturbations and morphing the collisionless NP mode into the magnetohydrodynamic (MHD) entropy mode. In fast waves, a fluctuating pressure anisotropy drives both mirror and firehose instabilities when the wave amplitude satisfies $δB/B_0 \gtrsim 2β^{-1}$. The induced particle scattering leads to delayed shock formation and MHD-like wave dynamics. Taken alongside prior work on self-interrupting Alfvén waves and self-sustaining ion-acoustic waves, our results establish a foundation for new theories of electromagnetic turbulence in low-collisionality, high-$β$ plasmas such as the intracluster medium, radiatively inefficient accretion flows, and the near-Earth solar wind. △ Less

Submitted 3 April, 2023; v1 submitted 5 January, 2023; originally announced January 2023.

Comments: 45 pages, 21 figures, accepted to Journal of Plasma Physics

A new buoyancy instability in galaxy clusters due to streaming cosmic rays

Authors: Philipp Kempski, Eliot Quataert, Jonathan Squire

Abstract: Active Galactic Nuclei (AGN) are believed to provide the energy that prevents runaway cooling of gas in the cores of galaxy clusters. However, how this energy is transported and thermalized throughout the Intracluster Medium (ICM) remains unclear. In recent work we showed that streaming cosmic rays (CRs) destabilise sound waves in dilute ICM plasmas. Here we show that CR streaming in the presence… ▽ More Active Galactic Nuclei (AGN) are believed to provide the energy that prevents runaway cooling of gas in the cores of galaxy clusters. However, how this energy is transported and thermalized throughout the Intracluster Medium (ICM) remains unclear. In recent work we showed that streaming cosmic rays (CRs) destabilise sound waves in dilute ICM plasmas. Here we show that CR streaming in the presence of gravity also destabilises a pressure-balanced wave. We term this new instability the CR buoyancy instability (CRBI). In stark contrast to standard results without CRs, the pressure-balanced mode is highly compressible at short wavelengths due to CR streaming. Maximal growth rates are of order $(p_c / p_g) β^{1/2} ω_{\rm ff}$, where $p_c/p_g$ is the ratio of CR pressure to thermal gas pressure, $β$ is the ratio of thermal to magnetic pressure and $ω_{\rm ff}$ is the free-fall frequency. The CRBI operates alongside buoyancy instabilities driven by background heat fluxes, i.e. the heat-flux-driven buoyancy instability (HBI) and the magneto-thermal instability (MTI). When the thermal mean free path $l_{\rm mfp}$ is $\ll$ the gas scale height $H$, the HBI/MTI set the growth rate on large scales, while the CRBI sets the growth rate on small scales. Conversely, when $l_{\rm mfp} \sim H$ and $(p_c/p_g) β^{1/2} \gtrsim 1$, CRBI growth rates exceed HBI/MTI growth rates even on large scales. Our results suggest that CR-driven instabilities may be partially responsible for the sound waves/weak shocks and turbulence observed in galaxy clusters. CR-driven instabilities generated near radio bubbles may also play an important role redistributing AGN energy throughout clusters. △ Less

Submitted 20 July, 2022; originally announced July 2022.

Comments: Submitted to MNRAS

Kinetic Turbulence in Collisionless High-Beta Plasmas

Authors: Lev Arzamasskiy, Matthew W. Kunz, Jonathan Squire, Eliot Quataert, Alexander A. Schekochihin

Abstract: We present results from three-dimensional hybrid-kinetic simulations of Alfvénic turbulence in a high-beta, collisionless plasma. The key feature of such turbulence is the interplay between local wave--wave interactions between the fluctuations in the cascade and the non-local wave-particle interactions associated with kinetic micro-instabilities driven by anisotropy in the thermal pressure (namel… ▽ More We present results from three-dimensional hybrid-kinetic simulations of Alfvénic turbulence in a high-beta, collisionless plasma. The key feature of such turbulence is the interplay between local wave--wave interactions between the fluctuations in the cascade and the non-local wave-particle interactions associated with kinetic micro-instabilities driven by anisotropy in the thermal pressure (namely, firehose, mirror, and ion-cyclotron). We present theoretical estimates for, and calculate directly from the simulations, the effective collisionality and plasma viscosity in pressure-anisotropic high-beta turbulence, demonstrating that, for strong Alfvénic turbulence, the effective parallel-viscous scale is comparable to the driving scale of the cascade. Below this scale, the kinetic-energy spectrum indicates an Alfvénic cascade with a slope steeper than $-5/3$ due to the anisotropic viscous stress. The magnetic-energy spectrum is shallower than $-5/3$ near the ion-Larmor scale due to fluctuations produced by the firehose instability. Most of the cascade energy (80-90%) is dissipated as ion heating through a combination of Landau damping and anisotropic viscous heating. Our results have implications for models of particle heating in low-luminosity accretion onto supermassive black holes, the effective viscosity of the intracluster medium, and the interpretation of near-Earth solar-wind observations. △ Less

Submitted 2 March, 2023; v1 submitted 11 July, 2022; originally announced July 2022.

Comments: 34 pages, 19 figures, accepted to PRX

On the construction of general large-amplitude spherically polarised Alfvén waves

Authors: Jonathan Squire, Alfred Mallet

Abstract: In a magnetised plasma on scales well above ion kinetic scales, any constant-magnitude magnetic field, accompanied by parallel Alfvénic flows, forms a nonlinear solution in an isobaric, constant-density background. These structures, which are also known as spherically polarised Alfvén waves, are observed ubiquitously in the solar wind, presumably created by the growth of small-amplitude fluctuatio… ▽ More In a magnetised plasma on scales well above ion kinetic scales, any constant-magnitude magnetic field, accompanied by parallel Alfvénic flows, forms a nonlinear solution in an isobaric, constant-density background. These structures, which are also known as spherically polarised Alfvén waves, are observed ubiquitously in the solar wind, presumably created by the growth of small-amplitude fluctuations as they propagate outwards in the corona. Here, we present a computational method to construct such solutions of arbitrary amplitude with general multi-dimensional structure, and explore some of their properties. The difficulty lies in computing a zero-divergence, constant-magnitude magnetic field, which leaves a single, quasi-free function to define the solution, while requiring strong constraints on any individual component of the field. Motivated by the physical process of wave growth in the solar wind, our method circumvents this issue by starting from low-amplitude Alfvénic fluctuations dominated by a strong mean field, then "growing" magnetic perturbations into the large-amplitude regime. We present example solutions with nontrivial structure in one, two, and three dimensions, demonstrating a clear tendency to form very sharp gradients or discontinuities, unless the solution is one dimensional. As well as being useful as an input for other calculations, particularly the study of parametric decay, our results provide a natural explanation for the extremely sharp field discontinuities observed across magnetic-field switchbacks in the low solar wind. △ Less

Submitted 30 August, 2022; v1 submitted 15 June, 2022; originally announced June 2022.

Comments: Accepted for publication in Journal of Plasma Physics Letters

On the properties of Alfvénic switchbacks in the expanding solar wind: the influence of the Parker spiral

Authors: Jonathan Squire, Zade Johnston, Alfred Mallet, Romain Meyrand

Abstract: Switchbacks -- rapid, large deflections of the solar wind's magnetic field -- have generated significant interest as possible signatures of the key mechanisms that heat the corona and accelerate the solar wind. In this context, an important task for theories of switchback formation and evolution is to understand their observable distinguishing features, allowing them to be assessed in detail using… ▽ More Switchbacks -- rapid, large deflections of the solar wind's magnetic field -- have generated significant interest as possible signatures of the key mechanisms that heat the corona and accelerate the solar wind. In this context, an important task for theories of switchback formation and evolution is to understand their observable distinguishing features, allowing them to be assessed in detail using spacecraft data. Here, we work towards this goal by studying the influence of the Parker spiral on the evolution of Alfvénic switchbacks in an expanding plasma. Using simple analytic arguments based on the physics of one-dimensional spherically polarized (constant-field-magnitude) Alfvén waves, we find that, by controlling the wave's obliquity, a Parker spiral strongly impacts switchback properties. Surprisingly, the Parker spiral can significantly enhance switchback formation, despite normalized wave amplitudes growing more slowly in its presence. In addition, switchbacks become strongly asymmetric: large switchbacks preferentially involve magnetic-field rotation in the plane of the Parker spiral (tangential deflections) rather than perpendicular (normal) rotations, and such deflections are strongly "tangentially skewed," meaning switchbacks always involve field rotations in the same direction (towards the positive-radial direction for an outwards mean field). In a companion paper, we show that these properties also occur in turbulent 3-D fields with switchbacks, with various caveats. These results demonstrate that substantial care is needed in assuming that specific features of switchbacks can be used to infer properties of the low corona; asymmetries and nontrivial correlations can develop as switchbacks propagate due to the interplay between expansion and spherically polarized, divergence-free magnetic fields. △ Less

Submitted 23 May, 2022; v1 submitted 19 May, 2022; originally announced May 2022.

On the properties of Alfvénic switchbacks in the expanding solar wind: three-dimensional numerical simulations

Authors: Zade Johnston, Jonathan Squire, Alfred Mallet, Romain Meyrand

Abstract: Switchbacks -- abrupt reversals of the magnetic field within the solar wind -- have been ubiquitously observed by Parker Solar Probe (PSP). Their origin, whether from processes near the solar surface or within the solar wind itself, remains under debate, and likely has key implications for solar wind heating and acceleration. Here, using three-dimensional expanding box simulations, we examine the… ▽ More Switchbacks -- abrupt reversals of the magnetic field within the solar wind -- have been ubiquitously observed by Parker Solar Probe (PSP). Their origin, whether from processes near the solar surface or within the solar wind itself, remains under debate, and likely has key implications for solar wind heating and acceleration. Here, using three-dimensional expanding box simulations, we examine the properties of switchbacks arising from the evolution of outwards-propagating Alfvén waves in the expanding solar wind in detail. Our goal is to provide testable predictions that can be used to differentiate between properties arising from solar surface processes and those from the in-situ evolution of Alfvén waves in switchback observations by PSP. We show how the inclusion of the Parker spiral causes magnetic field deflections within switchbacks to become asymmetric, preferentially deflecting in the plane of the Parker spiral and rotating in one direction towards the radial component of the mean field. The direction of the peak of the magnetic field distribution is also shown to be different from the mean field direction due to its highly skewed nature. Compressible properties of switchbacks are also explored, with magnetic-field-strength and density fluctuations being either correlated or anticorrelated depending on the value of $β$, agreeing with predictions from theory. We also measure dropouts in magnetic-field strength and density spikes at the boundaries of these synthetic switchbacks, both of which have been observed by PSP. The agreement of these properties with observations provide further support for the Alfvén wave model of switchbacks. △ Less

Submitted 2 December, 2022; v1 submitted 19 May, 2022; originally announced May 2022.

Comments: Revised argument in Section 2D, results unchanged; typos corrected

A Measurement of the Effective Mean-Free-Path of Solar Wind Protons

Authors: J. T. Coburn, C. H. K. Chen, J. Squire

Abstract: Weakly collisional plasmas are subject to nonlinear relaxation processes, which can operate at rates much faster than the particle collision frequencies. This causes the plasma to respond like a magnetised fluid despite having long particle mean-free-paths. In this Letter the effective collisional mechanisms are modelled in the plasma kinetic equation to produce density, pressure and magnetic fiel… ▽ More Weakly collisional plasmas are subject to nonlinear relaxation processes, which can operate at rates much faster than the particle collision frequencies. This causes the plasma to respond like a magnetised fluid despite having long particle mean-free-paths. In this Letter the effective collisional mechanisms are modelled in the plasma kinetic equation to produce density, pressure and magnetic field responses to compare with spacecraft measurements of the solar wind compressive fluctuations at 1 AU. This enables a measurement of the effective mean-free-path of the solar wind protons, found to be 4.35 $\times 10^5$ km, which is $\sim 10^3$ times shorter than the collisional mean-free-path. These measurements are shown to support the effective fluid behavior of the solar wind at scales above the proton gyroradius and demonstrate that effective collision processes alter the thermodynamics and transport of weakly collisional plasmas. △ Less

Submitted 24 March, 2022; originally announced March 2022.

Hydro-, Magnetohydro-, and Dust-Gas Dynamics of Protoplanetary Disks

Authors: G. Lesur, B. Ercolano, M. Flock, M. -K. Lin, C. -C. Yang, J. A. Barranco, P. Benitez-Llambay, J. Goodman, A. Johansen, H. Klahr, G. Laibe, W. Lyra, P. Marcus, R. P. Nelson, J. Squire, J. B. Simon, N. Turner, O. M. Umurhan, A. N. Youdin

Abstract: The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control t… ▽ More The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control the observational appearance and the secular evolution of these disks. This chapter is dedicated to the most recent developments in our understanding of the dynamics of gaseous and dusty disks, covering hydrodynamic and magnetohydrodynamic turbulence, gas-dust instabilities, dust clumping and disk winds. We show how these physical processes have been tested from observations and highlight standing questions that should be addressed in the future. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: 32 pages, 7 figures, to appear in Protostars and Planets VII, eds: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

On large-scale dynamos with stable stratification and the application to stellar radiative zones

Authors: Valentin Skoutnev, Jonathan Squire, Amitava Bhattacharjee

Abstract: Our understanding of large-scale magnetic fields in stellar radiative zones remains fragmented and incomplete. Such magnetic fields, which must be produced by some form of dynamo mechanism, are thought to dominate angular-momentum transport, making them crucial to stellar evolution. A major difficulty is the effect of stable stratification, which generally suppresses dynamo action. We explore the… ▽ More Our understanding of large-scale magnetic fields in stellar radiative zones remains fragmented and incomplete. Such magnetic fields, which must be produced by some form of dynamo mechanism, are thought to dominate angular-momentum transport, making them crucial to stellar evolution. A major difficulty is the effect of stable stratification, which generally suppresses dynamo action. We explore the effects of stable stratification on mean-field dynamo theory with a particular focus on a non-helical large-scale dynamo (LSD) mechanism known as the magnetic shear-current effect. We find that the mechanism is robust to increasing stable stratification as long as the original requirements for its operation are met: a source of shear and non-helical magnetic fluctuations (e.g. from a small-scale dynamo). Both are plausibly sourced in the presence of differential rotation. Our idealized direct numerical simulations, supported by mean-field theory, demonstrate the generation of near equipartition large-scale toroidal fields. Additionally, a scan over magnetic Reynolds number shows no change in the growth or saturation of the LSD, providing good numerical evidence of a dynamo mechanism resilient to catastrophic quenching, which has been an issue for helical dynamos. These properties -- the absence of catastrophic quenching and robustness to stable stratification -- make the mechanism a plausible candidate for generating in-situ large-scale magnetic fields in stellar radiative zones. △ Less

Submitted 16 September, 2022; v1 submitted 3 March, 2022; originally announced March 2022.

Standard Self-Confinement and Extrinsic Turbulence Models for Cosmic Ray Transport are Fundamentally Incompatible with Observations

Authors: Philip F. Hopkins, Jonathan Squire, Iryna S. Butsky, Suoqing Ji

Abstract: Models for cosmic ray (CR) dynamics fundamentally depend on the rate of CR scattering from magnetic fluctuations. In the ISM, for CRs with energies ~MeV-TeV, these fluctuations are usually attributed either to 'extrinsic turbulence' (ET) - a cascade from larger scales - or 'self-confinement' (SC) - self-generated fluctuations from CR streaming. Using simple analytic arguments and detailed live num… ▽ More Models for cosmic ray (CR) dynamics fundamentally depend on the rate of CR scattering from magnetic fluctuations. In the ISM, for CRs with energies ~MeV-TeV, these fluctuations are usually attributed either to 'extrinsic turbulence' (ET) - a cascade from larger scales - or 'self-confinement' (SC) - self-generated fluctuations from CR streaming. Using simple analytic arguments and detailed live numerical CR transport calculations in galaxy simulations, we show that both of these, in standard form, cannot explain even basic qualitative features of observed CR spectra. For ET, any spectrum that obeys critical balance or features realistic anisotropy, or any spectrum that accounts for finite damping below the dissipation scale, predicts qualitatively incorrect spectral shapes and scalings of B/C and other species. Even if somehow one ignored both anisotropy and damping, observationally-required scattering rates disagree with ET predictions by orders-of-magnitude. For SC, the dependence of driving on CR energy density means that it is nearly impossible to recover observed CR spectral shapes and scalings, and again there is an orders-of-magnitude normalization problem. But more severely, SC solutions with super-Alfvenic streaming are unstable. In live simulations, they revert to either arbitrarily-rapid CR escape with zero secondary production, or to bottleneck solutions with far-too-strong CR confinement and secondary production. Resolving these fundamental issues without discarding basic plasma processes requires invoking different drivers for scattering fluctuations. These must act on a broad range of scales with a power spectrum obeying several specific (but plausible) constraints. △ Less

Submitted 8 November, 2022; v1 submitted 3 December, 2021; originally announced December 2021.

Comments: 36 pages, 7 figures. Updated to match published version, added section discussing 'meso-scale' phenomenology

Numerical Study of Cosmic Ray Confinement through Dust Resonant Drag Instabilities

Authors: Suoqing Ji, Jonathan Squire, Philip F. Hopkins

Abstract: We investigate the possibility of cosmic ray (CR) confinement by charged dust grains through resonant drag instabilities (RDIs). We perform magnetohydrodynamic particle-in-cell simulations of magnetized gas mixed with charged dust and cosmic rays, with the gyro-radii of dust and GeV CRs on $\sim\mathrm{AU}$ scales fully resolved. As a first study, we focus on one type of RDI wherein charged grains… ▽ More We investigate the possibility of cosmic ray (CR) confinement by charged dust grains through resonant drag instabilities (RDIs). We perform magnetohydrodynamic particle-in-cell simulations of magnetized gas mixed with charged dust and cosmic rays, with the gyro-radii of dust and GeV CRs on $\sim\mathrm{AU}$ scales fully resolved. As a first study, we focus on one type of RDI wherein charged grains drift super-Alfv{é}nically, with Lorentz forces strongly dominating over drag forces. Dust grains are unstable to the RDIs and form concentrated columns and sheets, whose scale grows until saturating at the simulation box size. Initially perfectly-streaming CRs are strongly scattered by RDI-excited Alfv{é}n waves, with the growth rate of the CR perpendicular velocity components equaling the growth rate of magnetic field perturbations. These rates are well-predicted by analytic linear theory. CRs finally become isotropized and drift at least at $\sim v_\mathrm{A}$ by unidirectional Alfvén waves excited by the RDIs, with a uniform distribution of the pitch angle cosine $μ$ and a flat profile of the CR pitch angle diffusion coefficient $D_{μμ}$ around $μ= 0$, without the "$90$ degree pitch angle problem." With CR feedback on the gas included, $D_{μμ}$ decreases by a factor of a few, indicating a lower CR scattering rate, because the backreaction on the RDI from the CR pressure adds extra wave damping, leading to lower quasi-steady-state scattering rates. Our study demonstrates that the dust-induced CR confinement can be very important under certain conditions, e.g., the dusty circumgalactic medium around quasars or superluminous galaxies. △ Less

Submitted 22 April, 2022; v1 submitted 1 December, 2021; originally announced December 2021.

Comments: 14 pages, 10 figures, published on MNRAS

On the optical properties of Resonant Drag Instabilities: Variability of Asymptotic Giant Branch and R Coronae Borealis stars

Authors: Ulrich P. Steinwandel, Alexander A. Kaurov, Philip F. Hopkins, Jonathan Squire

Abstract: In dusty cool-star outflow or ejection events around AGB or RCB-like stars, dust is accelerated by radiation from the star and coupled to the gas via collisional drag forces. But it has recently been shown that such dust-gas mixtures are unstable to a super-class of instabilities called the resonant drag instabilities (RDIs), which promote dust clustering. We therefore consider idealized simulatio… ▽ More In dusty cool-star outflow or ejection events around AGB or RCB-like stars, dust is accelerated by radiation from the star and coupled to the gas via collisional drag forces. But it has recently been shown that such dust-gas mixtures are unstable to a super-class of instabilities called the resonant drag instabilities (RDIs), which promote dust clustering. We therefore consider idealized simulations of the RDIs operating on a spectrum of dust grain sizes subject to radiative acceleration (allowing for different grain optical properties), coupled to the gas with a realistic drag law, including or excluding the effects of magnetic fields and charged grains, and calculate for the first time how the RDIs could contribute to observed variability. We show that the RDIs naturally produce significant variations ($\sim 10-20\%$ $1σ$-level) in the extinction, corresponding to $\sim 0.1-1\,$mag level in the stellar types above, on timescales of order months to a year. The fluctuations are surprisingly robust to finite source-size effects as they are dominated by large-scale modes, which also means their spatial structure could be resolved in some nearby systems. We also quantify how this produces variations in the line-of-sight grain size-distribution. All of these variations are similar to those observed, suggesting that the RDIs may play a key role driving observed variability in dust extinction within dusty outflow/ejection events around cool stars. We further propose that the measured variations in grain sizes could directly be used to identify the presence of the RDIs in close by systems with observations. △ Less

Submitted 17 November, 2021; originally announced November 2021.

Comments: submitted to MNRAS, 13 pages, 12 Figures, comments welcome!

The In Situ Signature of Cyclotron Resonant Heating

Authors: Trevor A. Bowen, Benjmin D. G. Chandran, Jonathan Squire, Stuart D. Bale, Die Duan, Kristopher G. Klein, Davin Larson, Alfred Mallet, Michael D. McManus, Romain Meyrand, Jaye L. Verniero, Lloyd D. Woodham

Abstract: The dissipation of magnetized turbulence is an important paradigm for describing heating and energy transfer in astrophysical environments such as the solar corona and wind; however, the specific collisionless processes behind dissipation and heating remain relatively unconstrained by measurements. Remote sensing observations have suggested the presence of strong temperature anisotropy in the sola… ▽ More The dissipation of magnetized turbulence is an important paradigm for describing heating and energy transfer in astrophysical environments such as the solar corona and wind; however, the specific collisionless processes behind dissipation and heating remain relatively unconstrained by measurements. Remote sensing observations have suggested the presence of strong temperature anisotropy in the solar corona consistent with cyclotron resonant heating. In the solar wind, in situ magnetic field measurements reveal the presence of cyclotron waves, while measured ion velocity distribution functions have hinted at the active presence of cyclotron resonance. Here, we present Parker Solar Probe observations that connect the presence of ion-cyclotron waves directly to signatures of resonant damping in observed proton-velocity distributions. We show that the observed cyclotron wave population coincides with both flattening in the phase space distribution predicted by resonant quasilinear diffusion and steepening in the turbulent spectra at the ion-cyclotron resonant scale. In measured velocity distribution functions where cyclotron resonant flattening is weaker, the distributions are nearly uniformly subject to ion-cyclotron wave damping rather than emission, indicating that the distributions can damp the observed wave population. These results are consistent with active cyclotron heating in the solar wind. △ Less

Submitted 28 November, 2022; v1 submitted 9 November, 2021; originally announced November 2021.

Nonlinear Interactions in Spherically Polarized Alfvénic Turbulence

Authors: Trevor A. Bowen, Samuel T. Badman, Stuart D. Bale, Thierry Dudok de Wit, Timothy S. Horbury, Kristopher G. Klein, Davin Larson, Alfred Mallet, Lorenzo Matteini, Michael D. McManus, Jonathan Squire

Abstract: Turbulent magnetic field fluctuations observed in the solar wind often maintain a constant magnitude condition accompanied by spherically polarized velocity fluctuations; these signatures are characteristic of large-amplitude Alfvén waves. Nonlinear energy transfer in Alfvénic turbulence is typically considered in the small-amplitude limit where the constant magnitude condition may be neglected; i… ▽ More Turbulent magnetic field fluctuations observed in the solar wind often maintain a constant magnitude condition accompanied by spherically polarized velocity fluctuations; these signatures are characteristic of large-amplitude Alfvén waves. Nonlinear energy transfer in Alfvénic turbulence is typically considered in the small-amplitude limit where the constant magnitude condition may be neglected; in contrast, nonlinear energy transfer in the large-amplitude limit remains relatively unstudied. We develop a method to analyze finite-amplitude turbulence through studying fluctuations as constant magnitude rotations in the stationary wave (de Hoffmann-Teller) frame, which reveals that signatures of finite-amplitude effects exist deep into the MHD range. While the dominant fluctuations are consistent with spherically-polarized large-amplitude Alfvén waves, the subdominant mode is relatively compressible. Signatures of nonlinear interaction between the finite-amplitude spherically polarized mode with the subdominant population reveal highly aligned transverse components. In theoretical models of Alfvénic turbulence, alignment is thought to reduce nonlinearity; our observations require that alignment is sufficient to either reduce shear nonlinearity such that non-Alfvénic interactions may be responsible for energy transfer in spherically polarized states, or that counter-propagating fluctuations maintain anomalous coherence, which is a predicted signature of reflection-driven turbulence. △ Less

Submitted 21 October, 2021; originally announced October 2021.

The Acoustic Resonant Drag Instability with a Spectrum of Grain Sizes

Authors: Jonathan Squire, Stefania Moroianu, Philip F. Hopkins

Abstract: We study the linear growth and nonlinear saturation of the "acoustic Resonant Drag Instability" (RDI) when the dust grains, which drive the instability, have a wide, continuous spectrum of different sizes. This physics is generally applicable to dusty winds driven by radiation pressure, such as occurs around red-giant stars, star-forming regions, or active galactic nuclei. Depending on the physica… ▽ More We study the linear growth and nonlinear saturation of the "acoustic Resonant Drag Instability" (RDI) when the dust grains, which drive the instability, have a wide, continuous spectrum of different sizes. This physics is generally applicable to dusty winds driven by radiation pressure, such as occurs around red-giant stars, star-forming regions, or active galactic nuclei. Depending on the physical size of the grains compared to the wavelength of the radiation field that drives the wind, two qualitatively different regimes emerge. In the case of grains that are larger than the radiation's wavelength -- termed the constant-drift regime -- the grain's equilibrium drift velocity through the gas is approximately independent of grain size, leading to strong correlations between differently sized grains that persist well into the saturated nonlinear turbulence. For grains that are smaller than the radiation's wavelength -- termed the non-constant-drift regime -- the linear instability grows more slowly than the single-grain-size RDI and only the larger grains exhibit RDI-like behavior in the saturated state. A detailed study of grain clumping and grain-grain collisions shows that outflows in the constant-drift regime may be effective sites for grain growth through collisions, with large collision rates but low collision velocities. △ Less

Submitted 9 December, 2021; v1 submitted 21 October, 2021; originally announced October 2021.

Comments: Accepted for publication in MNRAS

High-frequency heating of the solar wind triggered by low-frequency turbulence

Authors: Jonathan Squire, Romain Meyrand, Matthew W. Kunz, Lev Arzamasskiy, Alexander A. Schekochihin, Eliot Quataert

Abstract: The fast solar wind's high speeds and nonthermal features require that significant heating occurs well above the Sun's surface. Two leading theories seem incompatible: low-frequency "Alfvénic" turbulence, which transports energy outwards and is observed ubiquitously by spacecraft but struggles to explain the observed dominance of ion over electron heating; and high-frequency ion-cyclotron waves (I… ▽ More The fast solar wind's high speeds and nonthermal features require that significant heating occurs well above the Sun's surface. Two leading theories seem incompatible: low-frequency "Alfvénic" turbulence, which transports energy outwards and is observed ubiquitously by spacecraft but struggles to explain the observed dominance of ion over electron heating; and high-frequency ion-cyclotron waves (ICWs), which explain the nonthermal heating of ions but lack an obvious source. Here, we argue that the recently proposed "helicity barrier" effect, which limits electron heating by inhibiting the turbulent cascade of energy to the smallest scales, can unify these two paradigms. Our six-dimensional simulations show how the helicity barrier causes the large-scale energy to grow in time, generating small parallel scales and high-frequency ICW heating from low-frequency turbulence. The resulting turbulence and ion distribution function also closely match in-situ measurements from Parker Solar Probe and other spacecraft, explaining, among other features, the decades-long puzzle of the steep "transition range" observed in magnetic fluctuation spectra. The theory predicts a causal link between plasma expansion and the ion-to-electron heating ratio. Given the observational association between wind speed and expansion, we argue that the helicity barrier could play a key role in regulating the bimodal speed distribution of the solar wind. △ Less

Submitted 23 February, 2022; v1 submitted 7 September, 2021; originally announced September 2021.

A solar source of Alfvénic magnetic field switchbacks: {\em in situ} remnants of magnetic funnels on supergranulation scales

Authors: S. D. Bale, T. S. Horbury, M. Velli, M. I. Desai, J. S. Halekas, M. D. McManus, O. Panasenco, S. T. Badman, T. A. Bowen, B. D. G. Chandran, J. F. Drake, J. C. Kasper, R. Laker, A. Mallet, L Matteini, T. D. Phan, N. E. Raouafi, J. Squire, L. D. Woodham, T. Wooley

Abstract: One of the striking observations from the Parker Solar Probe (PSP) spacecraft is the prevalence in the inner heliosphere of large amplitude, Alfvénic magnetic field reversals termed 'switchbacks'. These $δB_R/B \sim \mathcal{O}(1$) fluctuations occur on a range of timescales and in {\em patches} separated by intervals of quiet, radial magnetic field. We use measurements from PSP to demonstrate tha… ▽ More One of the striking observations from the Parker Solar Probe (PSP) spacecraft is the prevalence in the inner heliosphere of large amplitude, Alfvénic magnetic field reversals termed 'switchbacks'. These $δB_R/B \sim \mathcal{O}(1$) fluctuations occur on a range of timescales and in {\em patches} separated by intervals of quiet, radial magnetic field. We use measurements from PSP to demonstrate that patches of switchbacks are localized within the extensions of plasma structures originating at the base of the corona. These structures are characterized by an increase in alpha particle abundance, Mach number, plasma $β$ and pressure, and by depletions in the magnetic field magnitude and electron temperature. These intervals are in pressure-balance, implying stationary spatial structure, and the field depressions are consistent with overexpanded flux tubes. The structures are asymmetric in Carrington longitude with a steeper leading edge and a small ($\sim$1$^\circ$) edge of hotter plasma and enhanced magnetic field fluctuations. Some structures contain suprathermal ions to $\sim$85 keV that we argue are the energetic tail of the solar wind alpha population. The structures are separated in longitude by angular scales associated with supergranulation. This suggests that these switchbacks originate near the leading edge of the diverging magnetic field funnels associated with the network magnetic field - the primary wind sources. We propose an origin of the magnetic field switchbacks, hot plasma and suprathermals, alpha particles in interchange reconnection events just above the solar transition region and our measurements represent the extended regions of a turbulent outflow exhaust. △ Less

Submitted 2 September, 2021; originally announced September 2021.

Comments: 18 pages, 9 figures, submitted to Astrophys. J

Adaptive critical balance and firehose instability in an expanding, turbulent, collisionless plasma

Authors: A. F. A. Bott, L. Arzamasskiy, M. W. Kunz, E. Quataert, J. Squire

Abstract: Using hybrid-kinetic particle-in-cell simulation, we study the evolution of an expanding, collisionless, magnetized plasma in which strong Alfvénic turbulence is persistently driven. Temperature anisotropy generated adiabatically by the plasma expansion (and consequent decrease in the mean magnetic-field strength) gradually reduces the effective elasticity of the field lines, causing reductions in… ▽ More Using hybrid-kinetic particle-in-cell simulation, we study the evolution of an expanding, collisionless, magnetized plasma in which strong Alfvénic turbulence is persistently driven. Temperature anisotropy generated adiabatically by the plasma expansion (and consequent decrease in the mean magnetic-field strength) gradually reduces the effective elasticity of the field lines, causing reductions in the linear frequency and residual energy of the Alfvénic fluctuations. In response, these fluctuations modify their interactions and spatial anisotropy to maintain a scale-by-scale "critical balance" between their characteristic linear and nonlinear frequencies. Eventually the plasma becomes unstable to kinetic firehose instabilities, which excite rapidly growing magnetic fluctuations at ion-Larmor scales. The consequent pitch-angle scattering of particles maintains the temperature anisotropy near marginal stability, even as the turbulent plasma continues to expand. The resulting evolution of parallel and perpendicular temperatures does not satisfy double-adiabatic conservation laws, but is described accurately by a simple model that includes anomalous scattering. Our results have implications for understanding the complex interplay between macro- and micro-scale physics in various hot, dilute, astrophysical plasmas, and offer predictions concerning power spectra, residual energy, ion-Larmor-scale spectral breaks, and non-Maxwellian features in ion distribution functions that may be tested by measurements taken in high-beta regions of the solar wind. △ Less

Submitted 11 November, 2021; v1 submitted 3 August, 2021; originally announced August 2021.

Comments: ApJL, in press. 9 pages, 5 figures

Dust in the Wind with Resonant Drag Instabilities: I. The Dynamics of Dust-Driven Outflows in GMCs and HII Regions

Authors: Philip F. Hopkins, Anna L. Rosen, Jonathan Squire, Georgia V. Panopoulou, Nadine H. Soliman, Darryl Seligman, Ulrich P. Steinwandel

Abstract: Radiation-dust driven outflows, where radiation pressure on dust grains accelerates gas, occur in many astrophysical environments. Almost all previous numerical studies of these systems have assumed that the dust was perfectly-coupled to the gas. However, it has recently been shown that the dust in these systems is unstable to a large class of resonant drag instabilities (RDIs) which de-couple the… ▽ More Radiation-dust driven outflows, where radiation pressure on dust grains accelerates gas, occur in many astrophysical environments. Almost all previous numerical studies of these systems have assumed that the dust was perfectly-coupled to the gas. However, it has recently been shown that the dust in these systems is unstable to a large class of resonant drag instabilities (RDIs) which de-couple the dust and gas dynamics and could qualitatively change the nonlinear outcome of these outflows. We present the first simulations of radiation-dust driven outflows in stratified, inhomogeneous media, including explicit grain dynamics and a realistic spectrum of grain sizes and charge, magnetic fields and Lorentz forces on grains (which dramatically enhance the RDIs), Coulomb and Epstein drag forces, and explicit radiation transport allowing for different grain absorption and scattering properties. In this paper we consider conditions resembling giant molecular clouds (GMCs), HII regions, and distributed starbursts, where optical depths are modest, single-scattering effects dominate radiation-dust coupling, Lorentz forces dominate over drag on grains, and the fastest-growing RDIs are similar, such as magnetosonic and fast-gyro RDIs. These RDIs generically produce strong size-dependent dust clustering, growing nonlinear on timescales that are much shorter than the characteristic times of the outflow. The instabilities produce filamentary and plume-like or 'horsehead' nebular morphologies that are remarkably similar to observed dust structures in GMCs and HII regions. Additionally, in some cases they strongly alter the magnetic field structure and topology relative to filaments. Despite driving strong micro-scale dust clumping which leaves some gas behind, an order-unity fraction of the gas is always efficiently entrained by dust. △ Less

Submitted 8 November, 2022; v1 submitted 9 July, 2021; originally announced July 2021.

Comments: 25 pages, 22 figures. Updated to match published version

A Consistent Reduced-Speed-of-Light Formulation of Cosmic Ray Transport Valid in Weak and Strong-Scattering Regimes

Authors: Philip F. Hopkins, Jonathan Squire, Iryna S. Butsky

Abstract: We derive a consistent set of moments equations for CR-magnetohydrodynamics, assuming a gyrotropic distribution function (DF). Unlike previous efforts we derive a closure, akin to the M1 closure in radiation hydrodynamics (RHD), that is valid in both the nearly-isotropic-DF and/or strong-scattering regimes, and the arbitrarily-anisotropic DF or free-streaming regimes, as well as allowing for aniso… ▽ More We derive a consistent set of moments equations for CR-magnetohydrodynamics, assuming a gyrotropic distribution function (DF). Unlike previous efforts we derive a closure, akin to the M1 closure in radiation hydrodynamics (RHD), that is valid in both the nearly-isotropic-DF and/or strong-scattering regimes, and the arbitrarily-anisotropic DF or free-streaming regimes, as well as allowing for anisotropic scattering and transport/magnetic field structure. We present the appropriate two-moment closure and equations for various choices of evolved variables, including the CR phase space distribution function, number density, total energy, kinetic energy, and their fluxes or higher moments, and the appropriate coupling terms to the gas. We show that this naturally includes and generalizes a variety of terms including convection/fluid motion, anisotropic CR pressure, streaming, diffusion, gyro-resonant/streaming losses, and re-acceleration. We discuss how this extends previous treatments of CR transport including diffusion and moments methods and popular forms of the Fokker-Planck equation, as well as how this differs from the analogous M1-RHD equations. We also present two different methods for incorporating a reduced speed of light (RSOL) to reduce timestep limitations: in both we carefully address where the RSOL (versus true c) must appear for the correct behavior to be recovered in all interesting limits, and show how current implementations of CRs with a RSOL neglect some additional terms. △ Less

Submitted 16 March, 2022; v1 submitted 18 March, 2021; originally announced March 2021.

Comments: 17 pages, 2 figures. Readers not interested in the detailed derivations can skip to Appendix A, which simply summarizes the key equations derived. Updated to match published MNRAS version

Journal ref: MNRAS, 2022, 509, 3, 3779

The impact of astrophysical dust grains on the confinement of cosmic rays

Authors: Jonathan Squire, Philip F Hopkins, Eliot Quataert, Philipp Kempski

Abstract: We argue that charged dust grains could significantly impact the confinement and transport of galactic cosmic rays. For sub-GeV to ~1000GeV cosmic rays, small-scale parallel Alfvén waves, which isotropize cosmic rays through gyro-resonant interactions, are also gyro-resonant with charged grains. If the dust is nearly stationary, as in the bulk of the interstellar medium, Alfvén waves are damped by… ▽ More We argue that charged dust grains could significantly impact the confinement and transport of galactic cosmic rays. For sub-GeV to ~1000GeV cosmic rays, small-scale parallel Alfvén waves, which isotropize cosmic rays through gyro-resonant interactions, are also gyro-resonant with charged grains. If the dust is nearly stationary, as in the bulk of the interstellar medium, Alfvén waves are damped by dust. This will reduce the amplitude of Alfvén waves produced by the cosmic rays through the streaming instability, thus enhancing cosmic-ray transport. In well-ionized regions, the dust damping rate is larger by a factor of ~10 than other mechanisms that damp parallel Alfvén waves at the scales relevant for ~GeV cosmic rays, suggesting that dust could play a key role in regulating cosmic-ray transport. In astrophysical situations in which the dust moves through the gas with super-Alfvénic velocities, Alfvén waves are rendered unstable, which could directly scatter cosmic rays. This interaction has the potential to create a strong feedback mechanism where dust, driven through the gas by radiation pressure, then strongly enhances the confinement of cosmic rays, increasing their capacity to drive outflows. This mechanism may act in the circumgalactic medium around star-forming galaxies and active galactic nuclei. △ Less

Submitted 13 January, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Violation of the zeroth law of turbulence in space plasmas

Authors: Romain Meyrand, Jonathan Squire, Alexander A. Schekochihin, William Dorland

Abstract: The zeroth law of turbulence states that, for fixed energy input into large-scale motions, the statistical steady state of a turbulent system is independent of microphysical dissipation properties. The behavior, which is fundamental to nearly all fluid-like systems from industrial processes to galaxies, occurs because nonlinear processes generate smaller and smaller scales in the flow, until the d… ▽ More The zeroth law of turbulence states that, for fixed energy input into large-scale motions, the statistical steady state of a turbulent system is independent of microphysical dissipation properties. The behavior, which is fundamental to nearly all fluid-like systems from industrial processes to galaxies, occurs because nonlinear processes generate smaller and smaller scales in the flow, until the dissipation -- no matter how small -- can thermalize the energy input. Using direct numerical simulations and theoretical arguments, we show that in strongly magnetized plasma turbulence such as that recently observed by the Parker Solar Probe (PSP) spacecraft, the zeroth law is routinely violated. Namely, when such turbulence is "imbalanced" -- when the large-scale energy input is dominated by Alfvén waves propagating in one direction (the most common situation in space plasmas) -- nonlinear conservation laws imply the existence of a "barrier" at scales near the ion gyroradius. This causes energy to build up over time at large scales. The resulting magnetic-energy spectra bear a strong similarity to those observed in situ, exhibiting a sharp, steep kinetic transition range above and around the ion-Larmor scale, with flattening at yet smaller scales, thus resolving the decade-long puzzle of the position and variability of ion-kinetic spectral breaks in plasma turbulence. The "barrier" effect also suggests that how a plasma is forced at large scales (the imbalance) may have a crucial influence on thermodynamic properties such as the ion-to-electron heating ratio. △ Less

Submitted 10 June, 2021; v1 submitted 6 September, 2020; originally announced September 2020.

Small-Scale Dynamo in Stably Stratified Turbulence

Authors: Valentin Skoutnev, Jonathan Squire, Amitava Bhattacharjee

Abstract: We present numerical investigations into three principal properties of the small-scale dynamo in stably stratified turbulence: the onset criterion, the growth rate, and the nature of the magnetic field anisotropy in the kinematic regime. The results suggest that all three dynamo properties are controlled by the scale separation between the Ozmidov scale and the viscous or resistive scale. In addit… ▽ More We present numerical investigations into three principal properties of the small-scale dynamo in stably stratified turbulence: the onset criterion, the growth rate, and the nature of the magnetic field anisotropy in the kinematic regime. The results suggest that all three dynamo properties are controlled by the scale separation between the Ozmidov scale and the viscous or resistive scale. In addition to the critical magnetic Reynolds number, this allows for the definition of critical buoyancy and magnetic buoyancy Reynolds numbers for stratified small-scale dynamo onset in the high and low magnetic Prandtl number regimes, respectively. The presence of a small-scale dynamo in stellar radiative zones could affect dynamics through resulting Maxwell stresses and/or influence large-scale dynamo mechanisms in regions of differential rotation. Taking the solar radiative zone as a representative example and applying the onset criterion, we find that the stratification is strong enough to make the small-scale dynamo marginally active in the stably stratified turbulence of the solar tachocline. △ Less

Submitted 3 February, 2021; v1 submitted 3 August, 2020; originally announced August 2020.

Journal ref: ApJ 906 61 (2021)

Self-sustaining sound in collisionless, high-beta plasma

Authors: M. W. Kunz, J. Squire, A. A. Schekochihin, E. Quataert

Abstract: Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0 \gtrsim 2/β$ (where $β\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… ▽ More Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0 \gtrsim 2/β$ (where $β\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 $δn/n_0 \gtrsim 2/β$ 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. △ Less

Submitted 17 September, 2020; v1 submitted 16 June, 2020; originally announced June 2020.

Comments: 25 pages, 10 figures, Journal of Plasma Physics, in press

Effects of Different Cosmic Ray Transport Models on Galaxy Formation

Authors: Philip F. Hopkins, T. K. Chan, Jonathan Squire, Eliot Quataert, Suoqing Ji, Dusan Keres, Claude-Andre Faucher-Giguere

Abstract: 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… ▽ More 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 $γ$-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. △ Less

Submitted 16 March, 2022; v1 submitted 6 April, 2020; originally announced April 2020.

Comments: 6 pages, 2 figures. Updated to match published MNRAS version

Journal ref: MNRAS, 2021, 501, 3, 3663

Fluctuation dynamo in a weakly collisional plasma

Authors: D. A. St-Onge, M. W. Kunz, J. Squire, A. A. Schekochihin

Abstract: 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 t… ▽ More 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. △ Less

Submitted 13 July, 2020; v1 submitted 21 March, 2020; originally announced March 2020.

Comments: 62 pages, 25 figures, 1 table, accepted for publication to Journal of Plasma Physics

Journal ref: J. Plasma Phys. 86 (2020) 905860503

Physical models of streaming instabilities in protoplanetary disks

Authors: Jonathan Squire, Philip F. Hopkins

Abstract: 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 in… ▽ More 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. △ Less

Submitted 2 August, 2020; v1 submitted 3 March, 2020; originally announced March 2020.

Comments: Accepted for publication in MNRAS

Testing Physical Models for Cosmic Ray Transport Coefficients on Galactic Scales: Self-Confinement and Extrinsic Turbulence at GeV Energies

Authors: Philip F. Hopkins, Jonathan Squire, T. K. Chan, Eliot Quataert, Suoqing Ji, Dusan Keres, Claude-Andre Faucher-Giguere

Abstract: 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 prope… ▽ More 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 $γ$-rays ($L_γ$), 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_γ$), 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_γ$ 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. △ Less

Submitted 16 March, 2022; v1 submitted 14 February, 2020; originally announced February 2020.

Comments: 29 pages, 12 figures. Replaced with published MNRAS version

Journal ref: MNRAS, 2021, 501, 3

In-situ switchback formation in the expanding solar wind

Authors: Jonathan Squire, Benjamin D. G. Chandran, Romain Meyrand

Abstract: 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énic fluctuations advected by a radially expanding flow. Starting from simple superpositions of low-amplitude outward-propagating wave… ▽ More 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énic 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. △ Less

Submitted 6 February, 2020; v1 submitted 23 January, 2020; originally announced January 2020.

Sound-Wave Instabilities in Dilute Plasmas with Cosmic Rays: Implications for Cosmic-Ray Confinement and the Perseus X-ray Ripples

Authors: Philipp Kempski, Eliot Quataert, Jonathan Squire

Abstract: 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 Brag… ▽ More 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én speed, which in a weakly collisional plasma depends on the pressure anisotropy ($Δp$) of the thermal plasma. We find that this $Δ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 αβ^{-1/2}$, where $β$ is the ratio of thermal to magnetic pressure, and $α$, typically $\lesssim 1$, depends on other dimensionless parameters. In high-$β$ 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. △ Less

Submitted 20 February, 2020; v1 submitted 14 November, 2019; originally announced November 2019.

Comments: Accepted for publication in MNRAS; 14 Pages, 7 Figures

On the Survival of Cool Clouds in the Circum-Galactic Medium

Authors: Zhihui Li, Philip F. Hopkins, Jonathan Squire, Cameron Hummels

Abstract: 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,… ▽ More 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. △ Less

Submitted 18 December, 2019; v1 submitted 5 September, 2019; originally announced September 2019.

Comments: v1: 13 pages, 12 figures, submitted to MNRAS; v2: accepted version

Simulating the Diverse Instabilities of Dust in Magnetized Gas

Authors: Philip F. Hopkins, Jonathan Squire, Darryl Seligman

Abstract: 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 instabilit… ▽ More 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 $β$, 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. △ Less

Submitted 16 March, 2022; v1 submitted 23 April, 2019; originally announced April 2019.

Comments: 30 pages, 29 figures, Revised to match published MNRAS version. Animations and higher-resolution images available at http://www.tapir.caltech.edu/~phopkins/Site/animations/dust-and-gas-in-astrophysic/ . Simulation code available at http://www.tapir.caltech.edu/~phopkins/Site/GIZMO.html

Journal ref: MNRAS, 2020, 496, 2, 2123

[Plasma 2020 Decadal] Disentangling the Spatiotemporal Structure of Turbulence Using Multi-Spacecraft Data

Authors: J. M. TenBarge, O. Alexandrova, S. Boldyrev, F. Califano, S. S. Cerri, C. H. K. Chen, G. G. Howes, T. Horbury, P. A. Isenberg, H. Ji, K. G. Klein, C. Krafft, M. Kunz, N. F. Loureiro, A. Mallet, B. A. Maruca, W. H. Matthaeus, R. Meyrand, E. Quataert, J. C. Perez, O. W. Roberts, F. Sahraoui, C. S. Salem, A. A. Schekochihin, H. Spence , et al. (4 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 13 March, 2019; originally announced March 2019.

Comments: 6 pages, submitted for the Plasma 2020 Decadal Survey

[Plasma 2020 Decadal] The Material Properties of Weakly Collisional, High-Beta Plasmas

Authors: M. W. Kunz, J. Squire, S. A. Balbus, S. D. Bale, C. H. K. Chen, E. Churazov, S. C. Cowley, C. B. Forest, C. F. Gammie, E. Quataert, C. S. Reynolds, A. A. Schekochihin, L. Sironi, A. Spitkovsky, J. M. Stone, I. Zhuravleva, E. G. Zweibel

Abstract: 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:… ▽ More 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. △ Less

Submitted 10 March, 2019; originally announced March 2019.

Comments: 6 pages, submitted for the Plasma 2020 Decadal Survey

Shearing-box simulations of MRI-driven turbulence in weakly collisional accretion discs

Authors: Philipp Kempski, Eliot Quataert, Jonathan Squire, Matthew W. Kunz

Abstract: 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 ($Δp$) to lie within the stability bounds that would be otherwise imposed by kinetic microinstabilities. We explore a broad region of parameter space by co… ▽ More 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 ($Δ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 ($α_{\rm A} / α_{\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. △ Less

Submitted 16 April, 2019; v1 submitted 14 January, 2019; originally announced January 2019.

Comments: Accepted for publication in MNRAS

Magneto-immutable turbulence in weakly collisional plasmas

Authors: Jonathan Squire, Alexander A. Schekochihin, Eliot Quataert, Matthew W. Kunz

Abstract: 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, weakl… ▽ More 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-$β$) 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 $β\gtrsim1$ turbulent solar wind. △ Less

Submitted 14 January, 2019; v1 submitted 29 November, 2018; originally announced November 2018.

Comments: Accepted for publication in J. Plasma Phys

Nonlinear Evolution of the Resonant Drag Instability in Magnetized Gas

Authors: Darryl Seligman, Philip F. Hopkins, Jonathan Squire

Abstract: 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… ▽ More 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énic 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. △ Less

Submitted 22 October, 2018; originally announced October 2018.

Comments: 8 pages, 8 figures, submitted to MNRAS

Nonlinear Evolution of Instabilities Between Dust and Sound Waves

Authors: Eric R. Moseley, Jonathan Squire, Philip F. Hopkins

Abstract: 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 l… ▽ More 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. △ Less

Submitted 25 July, 2019; v1 submitted 18 October, 2018; originally announced October 2018.

Comments: 11 pages, 9 figures