BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jonathan Squire

Jonathan Squire appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.plasm-phastro-ph.HEastro-ph.SRastro-ph.GAphysics.flu-dynphysics.comp-phphysics.space-phastro-ph.COastro-ph.EPastro-ph.IMphysics.acc-ph

Catalog footprint

What is connected

24works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

24 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The magnetic shear-current effect: generation of large-scale magnetic fields by the small-scale dynamo

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Coherent nonhelical shear dynamos driven by magnetic fluctuations at low Reynolds numbers

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Electromotive force due to magnetohydrodynamic fluctuations in sheared rotating turbulence

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generation of large-scale magnetic fields by small-scale dynamo in shear flows

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gauge properties of the guiding center variational symplectic integrator

Variational symplectic algorithms have recently been developed for carrying out long-time simulation of charged particles in magnetic fields. As a direct consequence of their derivation from a discrete variational principle, these algorithms have very good long-time energy conservation, as well as exactly preserving discrete momenta. We present stability results for these algorithms, focusing on understanding how explicit variational integrators can be designed for this type of system. It is found that for explicit algorithms an instability arises because the discrete symplectic structure does not become the continuous structure in the $t \rightarrow 0$ limit. We examine how a generalized gauge transformation can be used to put the Lagrangian in the "antisymmetric discretization gauge," in which the discrete symplectic structure has the correct form, thus eliminating the numerical instability. Finally, it is noted that the variational guiding center algorithms are not electromagnetically gauge invariant. By designing a model discrete Lagrangian, we show that the algorithms are approximately gauge invariant as long as $\boldsymbol{A}$ and $ϕ$ are relatively smooth. A gauge invariant discrete Lagrangian is very important in a variational particle-in-cell algorithm where it ensures current continuity and preservation of Gauss's law

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Geometric integration of the Vlasov-Maxwell system with a variational particle-in-cell scheme

A fully variational, unstructured, electromagnetic particle-in-cell integrator is developed for integration of the Vlasov-Maxwell equations. Using the formalism of Discrete Exterior Calculus, the field solver, interpolation scheme and particle advance algorithm are derived through minimization of a single discrete field theory action. As a consequence of ensuring that the action is invariant under discrete electromagnetic gauge transformations, the integrator exactly conserves Gauss's law.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Magnetorotational instability: nonmodal growth and the relationship of global modes to the shearing box

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Nonmodal growth of the magnetorotational instability

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Statistical simulation of the magnetorotational dynamo

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

VEST: abstract vector calculus simplification in Mathematica

We present a new package, VEST (Vector Einstein Summation Tools), that performs abstract vector calculus computations in Mathematica. Through the use of index notation, VEST is able to reduce three-dimensional scalar and vector expressions of a very general type to a well defined standard form. In addition, utilizing properties of the Levi-Civita symbol, the program can derive types of multi-term vector identities that are not recognized by reduction, subsequently applying these to simplify large expressions. In a companion paper (Burby et al., 2013), we employ VEST in the automation of the calculation of high-order Lagrangians for the single particle guiding center system in plasma physics, a computation which illustrates its ability to handle very large expressions. VEST has been designed to be simple and intuitive to use, both for basic checking of work and more involved computations.

0 citations0 reviewsNo code linkedOpen access