Researcher profile

Jonathan Squire

Jonathan Squire contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2022arXiv

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

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.

preprint2022arXiv

Effects of Different Cosmic Ray Transport Models on Galaxy Formation

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Numerical Study of Cosmic Ray Confinement through Dust Resonant Drag Instabilities

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Simulating the Diverse Instabilities of Dust in Magnetized Gas

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.

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2021arXiv

Small-Scale Dynamo in Stably Stratified Turbulence

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.

preprint2021arXiv

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

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.

preprint2020arXiv

In-situ switchback formation in the expanding solar wind

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.

preprint2020arXiv

Physical models of streaming instabilities in protoplanetary disks

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.

preprint2020arXiv

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

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.

preprint2019arXiv

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

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.