We apply field-particle correlations -- a technique that tracks the time-averaged velocity-space structure of the energy density transfer rate between electromagnetic fields and plasma particles -- to data drawn from a hybrid Vlasov-Maxwell simulation of Alfv\'en Ion-Cyclotron turbulence. Energy transfer in this system is expected to include both Landau and cyclotron wave-particle resonances, unlike previous systems to which the field-particle correlation technique has been applied. In this simulation, the energy transfer rate mediated by the parallel electric field $E_\parallel$ comprises approximately $60\%$ of the total rate, with the remainder mediated by the perpendicular electric field $E_\perp$. The parallel electric field resonantly couples to protons, with the canonical bipolar velocity-space signature of Landau damping identified at many points throughout the simulation. The energy transfer mediated by $E_\perp$ preferentially couples to particles with $v_{tp} \lesssim v_\perp \lesssim 3 v_{tp}$ in agreement with the expected formation of a cyclotron diffusion plateau. Our results demonstrate clearly that the field-particle correlation technique can distinguish distinct channels of energy transfer using single-point measurements, even at points in which multiple channels act simultaneously, and can be used to determine quantitatively the rates of particle energization in each channel.

Data from the first two orbits of the Sun by Parker Solar Probe reveal that the solar wind sunward of 50 solar radii is replete with plasma waves and instabilities. One of the most prominent plasma wave power enhancements in this region appears near the electron cyclotron frequency (f_ce). Most of this wave power is concentrated in electric field fluctuations near 0.7 f_ce and f_ce, with strong harmonics of both frequencies extending above f_ce. At least two distinct, often concurrent, wave modes are observed, preliminarily identified as electrostatic whistler-mode waves and electron Bernstein waves. Wave intervals range in duration from a few seconds to hours. Both the amplitudes and number of detections of these near-f_ce waves increase significantly with decreasing distance to the Sun, suggesting that they play an important role in the evolution of electron populations in the near-Sun solar wind. Correlations are found between the detection of these waves and properties of solar wind electron populations, including electron core drift, implying that these waves play a role in regulating the heat flux carried by solar wind electrons. Observation of these near-f_ce waves is found to be strongly correlated with near-radial solar wind magnetic field configurations with low levels of magnetic turbulence. A scenario for the growth of these waves is presented which implies that regions of low-turbulence near-radial magnetic field are a prominent feature of solar wind structure near the Sun.

Understanding the physical processes in the solar wind and corona which actively contribute to heating, acceleration, and dissipation is a primary objective of NASA's Parker Solar Probe (PSP) mission. Observations of coherent electromagnetic waves at ion scales suggests that linear cyclotron resonance and non-linear processes are dynamically relevant in the inner heliosphere. A wavelet-based statistical study of coherent waves in the first perihelion encounter of PSP demonstrates the presence of transverse electromagnetic waves at ion resonant scales which are observed in 30-50\% of radial field intervals. Average wave amplitudes of approximately 4 nT are measured, while the mean duration of wave events is of order 20 seconds; however long duration wave events can exist without interruption on hour-long timescales. Though ion scale waves are preferentially observed during intervals with a radial mean magnetic field, we show that measurement constraints, associated with single spacecraft sampling of quasi-parallel waves superposed with anisotropic turbulence, render the measured quasi-parallel ion-wave spectrum unobservable when the mean magnetic field is oblique to the solar wind flow; these results imply that the occurrence of coherent ion-scale waves is not limited to a radial field configuration. The lack of strong radial scaling of characteristic wave amplitudes and duration suggests that the waves are generated {\em{in-situ}} through plasma instabilities. Additionally, observations of proton distribution functions indicate that temperature anisotropy may drive the observed ion-scale waves.

Turbulence is ubiquitously observed in nearly collisionless heliospheric plasmas, including the solar wind and corona and the Earth's magnetosphere. Understanding the collisionless mechanisms responsible for the energy transfer from the turbulent fluctuations to the particles is a frontier in kinetic turbulence research. Collisionless energy transfer from the turbulence to the particles can take place reversibly, resulting in non-thermal energy in the particle velocity distribution functions (VDFs) before eventual collisional thermalization is realized. Exploiting the information contained in the fluctuations in the VDFs is valuable. Here we apply a recently developed method based on VDFs, the field-particle correlation technique, to a $\beta=1$, solar-wind-like, low-frequency Alfv\'enic turbulence simulation with well resolved phase space to identify the field-particle energy transfer in velocity space. The field-particle correlations reveal that the energy transfer, mediated by the parallel electric field, results in significant structuring of the ion and electron VDFs in the direction parallel to the magnetic field. Fourier modes representing the length scales between the ion and electron gyroradii show that energy transfer is resonant in nature, localized in velocity space to the Landau resonances for each Fourier mode. The energy transfer closely follows the Landau resonant velocities with varying perpendicular wavenumber $k_\perp$ and plasma $\beta$. This resonant signature, consistent with Landau damping, is observed in all diagnosed Fourier modes that cover the dissipation range of the simulation.

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.

How turbulent energy is dissipated in weakly collisional space and astrophysical plasmas is a major open question. Here, we present the application of a field-particle correlation technique to directly measure the transfer of energy between the turbulent electromagnetic field and electrons in the Earth's magnetosheath, the region of solar wind downstream of the Earth's bow shock. The measurement of the secular energy transfer from the parallel electric field as a function of electron velocity shows a signature consistent with Landau damping. This signature is coherent over time, close to the predicted resonant velocity, similar to that seen in kinetic Alfv\'en turbulence simulations, and disappears under phase randomisation. This suggests that electron Landau damping could play a significant role in turbulent plasma heating, and that the technique is a valuable tool for determining the particle energisation processes operating in space and astrophysical plasmas.

Laboratory experiments provide a valuable complement to explore the fundamental physics of space plasmas without the limitations inherent to spacecraft measurements. Specifically, experiments overcome the restriction that spacecraft measurements are made at only one (or a few) points in space, enable greater control of the plasma conditions and applied perturbations, can be reproducible, and are orders of magnitude less expensive than launching spacecraft. Here I highlight key open questions about the physics of space plasmas and identify the aspects of these problems that can potentially be tackled in laboratory experiments. Several past successes in laboratory space physics provide concrete examples of how complementary experiments can contribute to our understanding of physical processes at play in the solar corona, solar wind, planetary magnetospheres, and outer boundary of the heliosphere. I present developments on the horizon of laboratory space physics, identifying velocity space as a key new frontier, highlighting new and enhanced experimental facilities, and showcasing anticipated developments to produce improved diagnostics and innovative analysis methods. A strategy for future laboratory space physics investigations will be outlined, with explicit connections to specific fundamental plasma phenomena of interest.

Understanding the removal of energy from turbulent fluctuations in a magnetized plasma and the consequent energization of the constituent plasma particles is a major goal of heliophysics and astrophysics. Previous work has shown that nonlinear interactions among counterpropagating Alfven waves---or Alfven wave collisions---are the fundamental building block of astrophysical plasma turbulence and naturally generate current sheets in the strongly nonlinear limit. A nonlinear gyrokinetic simulation of a strong Alfven wave collision is used to examine the damping of the electromagnetic fluctuations and the associated energization of particles that occurs in self-consistently generated current sheets. A simple model explains the flow of energy due to the collisionless damping and the associated particle energization, as well as the subsequent thermalization of the particle energy by collisions. The net particle energization by the parallel electric field is shown to be spatially intermittent, and the nonlinear evolution is essential in enabling that spatial non-uniformity. Using the recently developed field-particle correlation technique, we show that particles resonant with the Alfven waves in the simulation dominate the energy transfer, demonstrating conclusively that Landau damping plays a key role in the spatially intermittent damping of the electromagnetic fluctuations and consequent energization of the particles in this strongly nonlinear simulation.

Plasma turbulence occurs ubiquitously in space and astrophysical plasmas, mediating the nonlinear transfer of energy from large-scale electromagnetic fields and plasma flows to small scales at which the energy may be ultimately converted to plasma heat. But plasma turbulence also generically leads to a tangling of the magnetic field that threads through the plasma. The resulting wander of the magnetic field lines may significantly impact a number of important physical processes, including the propagation of cosmic rays and energetic particles, confinement in magnetic fusion devices, and the fundamental processes of turbulence, magnetic reconnection, and particle acceleration. The various potential impacts of magnetic field line wander are reviewed in detail, and a number of important theoretical considerations are identified that may influence the development and saturation of magnetic field line wander in astrophysical plasma turbulence. The results of nonlinear gyrokinetic simulations of kinetic Alfven wave turbulence of sub-ion length scales are evaluated to understand the development and saturation of the turbulent magnetic energy spectrum and of the magnetic field line wander. It is found that turbulent space and astrophysical plasmas are generally expected to contain a stochastic magnetic field due to the tangling of the field by strong plasma turbulence. Future work will explore how the saturated magnetic field line wander varies as a function of the amplitude of the plasma turbulence and the ratio of the thermal to magnetic pressure, known as the plasma beta.

Under the low density and high temperature conditions typical of heliospheric plasmas, the macroscopic evolution of the heliosphere is strongly affected by the kinetic plasma physics governing fundamental microphysical mechanisms. Kinetic turbulence, collision less magnetic reconnection, particle acceleration, and kinetic instabilities are four poorly understood, grand-challenge problems that lie at the new frontier of kinetic heliophysics. The increasing availability of high cadence and high phase-space resolution measurements of particle velocity distributions by current and upcoming spacecraft missions and of massively parallel nonlinear kinetic simulations of weakly collisional heliospheric plasmas provides the opportunity to transform our understanding of these kinetic mechanisms through the full utilization of the information contained in the particle velocity distributions. Several major considerations for future investigations of kinetic heliophysics are examined. Turbulent dissipation followed by particle heating is highlighted as an inherently two-step process in weakly collisional plasmas, distinct from the more familiar case in fluid theory. Concerted efforts must be made to tackle the big-data challenge of visualizing the high-dimensional (3D-3V) phase space of kinetic plasma theory through physics-based reductions. Furthermore, the development of innovative analysis methods that utilize full velocity-space measurements, such as the field-particle correlation technique, will enable us to gain deeper insight into these four grand-challenge problems of kinetic heliophysics. A systems approach to tackle the multi-scale problem of heliophysics through a rigorous connection between the kinetic physics at microscales and the self-consistent evolution of the heliosphere at macro scales will propel the field of kinetic heliophysics into the future.

Determining the physical mechanisms that extract energy from turbulent fluctuations in weakly collisional magnetized plasmas is necessary for a more complete characterization of the behavior of a variety of space and astrophysical plasmas. Such a determination is complicated by the complex nature of the turbulence as well as observational constraints, chiefly that in situ measurements of such plasmas are typically only available at a single point in space. Recent work has shown that correlations between electric fields and particle velocity distributions constructed from single-point measurements produce a velocity-dependent signature of the collisionless damping mechanism. We extend this work by constructing field-particle correlations using data sets drawn from single points in strongly driven, turbulent, electromagnetic gyrokinetic simulations to demonstrate that this technique can identify the collisionless mechanisms operating in such systems. The correlation's velocity-space structure agrees with expectations of resonant mechanisms transferring energy collisionlessly in turbulent systems. This work motivates the eventual application of field-particle correlations to spacecraft measurements in the solar wind, with the ultimate goal to determine the physical mechanisms that dissipate magnetized plasma turbulence.

Plasma turbulence is ubiquitous in space and astrophysical plasmas, playing an important role in plasma energization, but the physical mechanisms leading to dissipation of the turbulent energy remain to be definitively identified. Kinetic simulations in two dimensions (2D) have been extensively used to study the dissipation process. How the limitation to 2D affects energy dissipation remains unclear. This work provides a model of comparison between two- and three-dimensional (3D) plasma turbulence using gyrokinetic simulations; it also explores the dynamics of distribution functions during the dissipation process. It is found that both 2D and 3D nonlinear gyrokinetic simulations of a low-beta plasma generate electron velocity-space structures with the same characteristics as that of linear Landau damping of Alfv\'en waves in a 3D linear simulation. The continual occurrence of the velocity-space structures throughout the turbulence simulations suggests that the action of Landau damping may be responsible for the turbulent energy transfer to electrons in both 2D and 3D, and makes possible the subsequent irreversible heating of the plasma through collisional smoothing of the velocity-space fluctuations. Although, in the 2D case where variation along the equilibrium magnetic field is absent, it may be expected that Landau damping is not possible, a common trigonometric factor appears in the 2D resonant denominator, leaving the resonance condition unchanged from the 3D case. The evolution of the 2D and 3D cases is qualitatively similar. However, quantitatively the nonlinear energy cascade and subsequent dissipation is significantly slower in the 2D case.

Turbulence profoundly affects particle transport and plasma heating in many astrophysical plasma environments, from galaxy clusters to the solar corona and solar wind to Earth's magnetosphere. Both fluid and kinetic simulations of plasma turbulence ubiquitously generate coherent structures, in the form of current sheets, at small scales, and the locations of these current sheets appear to be associated with enhanced rates of dissipation of the turbulent energy. Therefore, illuminating the origin and nature of these current sheets is critical to identifying the dominant physical mechanisms of dissipation, a primary aim at the forefront of plasma turbulence research. Here we present evidence from nonlinear gyrokinetic simulations that strong nonlinear interactions between counterpropagating Alfven waves, or strong Alfven wave collisions, are a natural mechanism for the generation of current sheets in plasma turbulence. Furthermore, we conceptually explain this current sheet development in terms of the nonlinear dynamics of Alfven wave collisions, showing that these current sheets arise through constructive interference among the initial Alfven waves and nonlinearly generated modes. The properties of current sheets generated by a strong Alfven wave collisions are compared to published observations of current sheets in the Earth's magnetosheath and the solar wind, and the nature of these current sheets leads to the expectation that Landau damping of the constituent Alfven waves plays a dominant role in the damping of turbulently generated current sheets.

An innovative field-particle correlation technique is proposed that uses single-point measurements of the electromagnetic fields and particle velocity distribution functions to investigate the net transfer of energy from fields to particles associated with the collisionless damping of turbulent fluctuations in weakly collisional plasmas, such as the solar wind. In addition to providing a direct estimate of the local rate of energy transfer between fields and particles, it provides vital new information about the distribution of that energy transfer in velocity space. This velocity-space signature can potentially be used to identify the dominant collisionless mechanism responsible for the damping of turbulent fluctuations in the solar wind. The application of this novel field-particle correlation technique is illustrated using the simplified case of the Landau damping of Langmuir waves in an electrostatic 1D-1V Vlasov-Poisson plasma, showing that the procedure both estimates the local rate of energy transfer from the electrostatic field to the electrons and indicates the resonant nature of this interaction. Modifications of the technique to enable single-point spacecraft measurements of fields and particles to diagnose the collisionless damping of turbulent fluctuations in the solar wind are discussed, yielding a method with the potential to transform our ability to maximize the scientific return from current and upcoming spacecraft missions, such as the Magnetospheric Multiscale (MMS) and Solar Probe Plus missions.

Particle velocity distributions measured in the weakly collisional solar wind are frequently found to be non-Maxwellian, but how these non-Maxwellian distributions impact the physics of plasma turbulence in the solar wind remains unanswered. Using numerical solutions of the linear dispersion relation for a collisionless plasma with a bi-Maxwellian proton velocity distribution, we present a unified framework for the four proton temperature anisotropy instabilities, identifying the associated stable eigenmodes, highlighting the unstable region of wavevector space, and presenting the properties of the growing eigenfunctions. Based on physical intuition gained from this framework, we address how the proton temperature anisotropy impacts the nonlinear dynamics of the \Alfvenic fluctuations underlying the dominant cascade of energy from large to small scales and how the fluctuations driven by proton temperature anisotropy instabilities interact nonlinearly with each other and with the fluctuations of the large-scale cascade. We find that the nonlinear dynamics of the large-scale cascade is insensitive to the proton temperature anisotropy, and that the instability-driven fluctuations are unlikely to cause significant nonlinear evolution of either the instability-driven fluctuations or the turbulent fluctuations of the large-scale cascade.

A dynamical approach, rather than the usual statistical approach, is taken to explore the physical mechanisms underlying the nonlinear transfer of energy, the damping of the turbulent fluctuations, and the development of coherent structures in kinetic plasma turbulence. It is argued that the linear and nonlinear dynamics of Alfven waves are responsible, at a very fundamental level, for some of the key qualitative features of plasma turbulence that distinguish it from hydrodynamic turbulence, including the anisotropic cascade of energy and the development of current sheets at small scales. The first dynamical model of kinetic turbulence in the weakly collisional solar wind plasma that combines self-consistently the physics of Alfven waves with the development of small-scale current sheets is presented and its physical implications are discussed. This model leads to a simplified perspective on the nature of turbulence in a weakly collisional plasma: the nonlinear interactions responsible for the turbulent cascade of energy and the formation of current sheets are essentially fluid in nature, while the collisionless damping of the turbulent fluctuations and the energy injection by kinetic instabilities are essentially kinetic in nature.

The weak collisionality typical of turbulence in many diffuse astrophysical plasmas invalidates an MHD description of the turbulent dynamics, motivating the development of a more comprehensive theory of kinetic turbulence. In particular, a kinetic approach is essential for the investigation of the physical mechanisms responsible for the dissipation of astrophysical turbulence and the resulting heating of the plasma. This chapter reviews the limitations of MHD turbulence theory and explains how kinetic considerations may be incorporated to obtain a kinetic theory for astrophysical plasma turbulence. Key questions about the nature of kinetic turbulence that drive current research efforts are identified. A comprehensive model of the kinetic turbulent cascade is presented, with a detailed discussion of each component of the model and a review of supporting and conflicting theoretical, numerical, and observational evidence.

Motivated by the upcoming Solar Orbiter and Solar Probe Plus missions, qualitative and quantitative predictions are made for the effects of the violation of the Taylor hypothesis on the magnetic energy frequency spectrum measured in the near-Sun environment. The synthetic spacecraft data method is used to predict observational signatures of the violation for critically balanced Alfv\'enic turbulence or parallel fast/whistler turbulence. The violation of the Taylor hypothesis can occur in the slow flow regime, leading to a shift of the entire spectrum to higher frequencies, or in the dispersive regime, in which the dissipation range spectrum flattens at high frequencies. It is found that Alfv\'enic turbulence will not significantly violate the Taylor hypothesis, but whistler turbulence will. The flattening of the frequency spectrum is therefore a key observational signature for fast/whistler turbulence.

Turbulence in astrophysical and space plasmas is dominated by the nonlinear interaction of counterpropagating Alfven waves. Most Alfven wave turbulence theories have been based on ideal plasma models, such as incompressible MHD, for Alfven waves at large scales. However, in the inertial Alfven wave regime (vA > vthe), relevant to magnetospheric plasmas, how the turbulent nonlinear interactions are modified by the dispersive nature of the waves remains to be explored. Here we present the first laboratory evidence of the nonlinear interaction in the inertial regime. A comparison is made with the theory for MHD Alfven waves.

It is often asserted or implicitly assumed, without justification, that the results of two-dimensional investigations of plasma turbulence are applicable to the three-dimensional plasma environments of interest. A projection method is applied to derive two scalar equations that govern the nonlinear evolution of the Alfvenic and pseudo-Alfvenic components of ideal incompressible magnetohydrodynamic (MHD) plasma turbulence. The mathematical form of these equations makes clear the inherently three-dimensional nature of plasma turbulence, enabling an analysis of the nonlinear properties of two-dimensional limits often used to study plasma turbulence. In the anisotropic limit k_perp >>k_parallel that naturally arises in magnetized plasma systems, the perpendicular 2D limit retains the dominant nonlinearities that are mediated only by the Alfvenic fluctuations but lacks the wave physics associated with the linear term that is necessary to capture the anisotropic cascade of turbulent energy. In the in-plane 2D limit, the nonlinear energy transfer is controlled instead by the pseudo-Alfven waves, with the Alfven waves relegated to a passive role. In the oblique 2D limit, an unavoidable azimuthal dependence connecting the wavevector components will likely cause artificial azimuthal asymmetries in the resulting turbulent dynamics. Therefore, none of these 2D limits is sufficient to capture fully the rich three-dimensional nonlinear dynamics critical to the evolution of plasma turbulence.

The interpretation of single-point spacecraft measurements of solar wind turbulence is complicated by the fact that the measurements are made in a frame of reference in relative motion with respect to the turbulent plasma. The Taylor hypothesis---that temporal fluctuations measured by a stationary probe in a rapidly flowing fluid are dominated by the advection of spatial structures in the fluid rest frame---is often assumed to simplify the analysis. But measurements of turbulence in upcoming missions, such as Solar Probe Plus, threaten to violate the Taylor hypothesis, either due to slow flow of the plasma with respect to the spacecraft or to the dispersive nature of the plasma fluctuations at small scales. Assuming that the frequency of the turbulent fluctuations is characterized by the frequency of the linear waves supported by the plasma, we evaluate the validity of the Taylor hypothesis for the linear kinetic wave modes in the weakly collisional solar wind. The analysis predicts that a dissipation range of solar wind turbulence supported by whistler waves is likely to violate the Taylor hypothesis, while one supported by kinetic Alfven waves is not.

The quasilinear premise is a hypothesis for the modeling of plasma turbulence in which the turbulent fluctuations are represented by a superposition of randomly-phased linear wave modes, and energy is transferred among these wave modes via nonlinear interactions. We define specifically what constitutes the quasilinear premise, and present a range of theoretical arguments in support of the relevance of linear wave properties even in a strongly turbulent plasma. We review evidence both in support of and in conflict with the quasilinear premise from numerical simulations and measurements of plasma turbulence in the solar wind. Although the question of the validity of the quasilinear premise remains to be settled, we suggest that the evidence largely supports the value of the quasilinear premise in modeling plasma turbulence and that its usefulness may also be judged by the insights gained from such an approach, with the ultimate goal to develop the capability to predict the evolution of any turbulent plasma system, including the spectrum of turbulent fluctuations, their dissipation, and the resulting plasma heating.

Motivated by recent observations of distinct parallel and perpendicular signatures in magnetic helicity measurements segregated by wave period and angle between the local magnetic field and the solar wind velocity, this paper undertakes a comparison of three intervals of \emph{Ulysses} data with synthetic time series generated from a physically motivated turbulence model. From these comparisons, it is hypothesized that the observed signatures result from a perpendicular cascade of \Alfvenic fluctuations and a local, non-turbulent population of ion cyclotron or whistler waves generated by temperature anisotropy instabilities. By constraining the model's free parameters through comparison to \emph{in situ} data, it is found that, on average, $ \sim 95\%$ of the power near dissipative scales is contained in a perpendicular \Alfvenic cascade and that the parallel fluctuations are propagating nearly unidirectionally. The effects of aliasing on magnetic helicity measurements are considered and shown to be significant near the Nyquist frequency.

Results of the first validation of large guide field, $B_g / \delta B_0 \gg 1$, gyrokinetic simulations of magnetic reconnection at a fusion and solar corona relevant $\beta_i = 0.01$ and solar wind relevant $\beta_i = 1$ are presented, where $\delta B_0$ is the reconnecting field. Particle-in-cell (PIC) simulations scan a wide range of guide magnetic field strength to test for convergence to the gyrokinetic limit. The gyrokinetic simulations display a high degree of morphological symmetry, to which the PIC simulations converge when $\beta_i B_g / \delta B_0 \gtrsim 1$ and $B_g / \delta B_0 \gg 1$. In the regime of convergence, the reconnection rate, relative energy conversion, and overall magnitudes are found to match well between the PIC and gyrokinetic simulations, implying that gyrokinetics is capable of making accurate predictions well outside its regime of formal applicability. These results imply that in the large guide field limit many quantities resulting from the nonlinear evolution of reconnection scale linearly with the guide field.

A unique method of driving Alfvenic turbulence via an oscillating Langevin antenna is presented. This method of driving is motivated by a desire to inject energy into a finite domain numerical simulation in a manner that models the nonlinear transfer of energy from fluctuations in the turbulent cascade at scales larger than the simulation domain.. The oscillating Langevin antenna is shown to capture the essential features of the larger scale turbulence and efficiently couple to the plasma, generating steady-state turbulence within one characteristic turnaround time. The antenna is also sufficiently flexible to explore both strong and weak regimes of Alfvenic plasma turbulence.

Turbulence in space and astrophysical plasmas is governed by the nonlinear interactions between counterpropagating Alfven waves. Here we present the theoretical considerations behind the design of the first laboratory measurement of an Alfven wave collision, the fundamental interaction underlying Alfvenic turbulence. By interacting a relatively large-amplitude, low-frequency Alfven wave with a counterpropagating, smaller-amplitude, higher-frequency Alfven wave, the experiment accomplishes the secular nonlinear transfer of energy to a propagating daughter Alfven wave. The predicted properties of the nonlinearly generated daughter Alfven wave are outlined, providing a suite of tests that can be used to confirm the successful measurement of the nonlinear interaction between counterpropagating Alfven waves in the laboratory.

We present the first study of the formation and dissipation of current sheets at electron scales in a wave-driven, weakly collisional, 3D kinetic turbulence simulation. We investigate the relative importance of dissipation associated with collisionless damping via resonant wave-particle interactions versus dissipation in small-scale current sheets in weakly collisional plasma turbulence. Current sheets form self-consistently from the wave-driven turbulence, and their filling fraction is well correlated to the electron heating rate. However, the weakly collisional nature of the simulation necessarily implies that the current sheets are not significantly dissipated via Ohmic dissipation. Rather, collisionless damping via the Landau resonance with the electrons is sufficient to account for the measured heating as a function of scale in the simulation, without the need for significant Ohmic dissipation. This finding suggests the possibility that the dissipation of the current sheets is governed by resonant wave-particle interactions and that the locations of current sheets correspond spatially to regions of enhanced heating.

The nonlinear interaction between counterpropagating Alfven waves is the physical mechanism underlying the cascade of energy to small scales in astrophysical plasma turbulence. Beginning with the equations for incompressible MHD, an asymptotic analytical solution for the nonlinear evolution of these Alfven wave collisions is derived in the weakly nonlinear limit. The resulting qualitative picture of nonlinear energy transfer due to this mechanism involves two steps: first, the primary counterpropagating Alfven waves interact to generate an inherently nonlinear, purely magnetic secondary fluctuation with no parallel variation; second, the two primary waves each interact with this secondary fluctuation to transfer energy secularly to two tertiary Alfven waves. These tertiary modes are linear Alfven waves with the same parallel wavenumber as the primary waves, indicating the lack of a parallel cascade. The amplitude of these tertiary modes increases linearly with time due to the coherent nature of the resonant four-wave interaction responsible for the nonlinear energy transfer. The implications of this analytical solution for turbulence in astrophysical plasmas is discussed. The solution presented here provides valuable intuition about the nonlinear interactions underlying magnetized plasma turbulence, in support of an experimental program to verify in the laboratory the nature of this fundamental building block of astrophysical plasma turbulence.

This paper presents the numerical verification of an asymptotic analytical solution for the nonlinear interaction between counterpropagating Alfven waves, the fundamental building block of astrophysical plasma turbulence. The analytical solution, derived in the weak turbulence limit using the equations of incompressible MHD, is compared to a nonlinear gyrokinetic simulation of an Alfven wave collision. The agreement between these methods signifies that the incompressible solution satisfactorily describes the essential dynamics of the nonlinear energy transfer, even under the weakly collisional plasma conditions relevant to many astrophysical environments.

Turbulence is a phenomenon found throughout space and astrophysical plasmas. It plays an important role in solar coronal heating, acceleration of the solar wind, and heating of the interstellar medium. Turbulence in these regimes is dominated by Alfven waves. Most turbulence theories have been established using ideal plasma models, such as incompressible MHD. However, there has been no experimental evidence to support the use of such models for weakly to moderately collisional plasmas which are relevant to various space and astrophysical plasma environments. We present the first experiment to measure the nonlinear interaction between two counterpropagating Alfven waves, which is the building block for astrophysical turbulence theories. We present here four distinct tests that demonstrate conclusively that we have indeed measured the daughter Alfven wave generated nonlinearly by a collision between counterpropagating Alfven waves.

Turbulence is a ubiquitous phenomenon in space and astrophysical plasmas, driving a cascade of energy from large to small scales and strongly influencing the plasma heating resulting from the dissipation of the turbulence. Modern theories of plasma turbulence are based on the fundamental concept that the turbulent cascade of energy is caused by the nonlinear interaction between counterpropagating Alfven waves, yet this interaction has never been observationally or experimentally verified. We present here the first experimental measurement in a laboratory plasma of the nonlinear interaction between counterpropagating Alfven waves, the fundamental building block of astrophysical plasma turbulence. This measurement establishes a firm basis for the application of theoretical ideas developed in idealized models to turbulence in realistic space and astrophysical plasma systems.

We motivate the importance of studying kinetic scale turbulence for understanding the macroscopic properties of the heliosphere, such as the heating of the solar wind. We then discuss the technique by which kinetic scale density fluctuations can be measured using the spacecraft potential, including a calculation of the timescale for the spacecraft potential to react to the density changes. Finally, we compare the shape of the density spectrum at ion scales to theoretical predictions based on a cascade model for kinetic turbulence. We conclude that the shape of the spectrum, including the ion scale flattening, can be captured by the sum of passive density fluctuations at large scales and kinetic Alfven wave turbulence at small scales.

Kinetic plasma theory is used to generate synthetic spacecraft data to analyze and interpret the compressible fluctuations in the inertial range of solar wind turbulence. The kinetic counterparts of the three familiar linear MHD wave modes---the fast, Alfven, and slow waves---are identified and the properties of the density-parallel magnetic field correlation for these kinetic wave modes is presented. The construction of synthetic spacecraft data, based on the quasi-linear premise---that some characteristics of magnetized plasma turbulence can be usefully modeled as a collection of randomly phased, linear wave modes---is described in detail. Theoretical predictions of the density-parallel magnetic field correlation based on MHD and Vlasov-Maxwell linear eigenfunctions are presented and compared to the observational determination of this correlation based on 10 years of Wind spacecraft data. It is demonstrated that MHD theory is inadequate to describe the compressible turbulent fluctuations and that the observed density-parallel magnetic field correlation is consistent with a statistically negligible kinetic fast wave energy contribution for the large sample used in this study. A model of the solar wind inertial range fluctuations is proposed comprised of a mixture of a critically balanced distribution of incompressible Alfvenic fluctuations and a critically balanced or more anisotropic than critical balance distribution of compressible slow wave fluctuations. These results imply that there is little or no transfer of large scale turbulent energy through the inertial range down to whistler waves at small scales.

The magnetic variance anisotropy ($\mathcal{A}_m$) of the solar wind has been used widely as a method to identify the nature of solar wind turbulent fluctuations; however, a thorough discussion of the meaning and interpretation of the $\mathcal{A}_m$ has not appeared in the literature. This paper explores the implications and limitations of using the $\mathcal{A}_m$ as a method for constraining the solar wind fluctuation mode composition and presents a more informative method for interpreting spacecraft data. The paper also compares predictions of the $\mathcal{A}_m$ from linear theory to nonlinear turbulence simulations and solar wind measurements. In both cases, linear theory compares well and suggests the solar wind for the interval studied is dominantly Alfv\'{e}nic in the inertial and dissipation ranges to scales $k \rho_i \simeq 5$.

A numerical simulation of kinetic plasma turbulence is performed to assess the applicability of critical balance to kinetic, dissipation scale turbulence. The analysis is performed in the frequency domain to obviate complications inherent in performing a local analysis of turbulence. A theoretical model of dissipation scale critical balance is constructed and compared to simulation results, and excellent agreement is found. This result constitutes the first evidence of critical balance in a kinetic turbulence simulation and provides evidence of an anisotropic turbulence cascade extending into the dissipation range. We also perform an Eulerian frequency analysis of the simulation data and compare it to the results of a previous study of magnetohydrodynamic turbulence simulations.

A refined cascade model for kinetic turbulence in weakly collisional astrophysical plasmas is presented that includes both the transition between weak and strong turbulence and the effect of nonlocal interactions on the nonlinear transfer of energy. The model describes the transition between weak and strong MHD turbulence and the complementary transition from strong kinetic Alfven wave (KAW) turbulence to weak dissipating KAW turbulence, a new regime of weak turbulence in which the effects of shearing by large scale motions and kinetic dissipation play an important role. The inclusion of the effect of nonlocal motions on the nonlinear energy cascade rate in the dissipation range, specifically the shearing by large-scale motions, is proposed to explain the nearly power-law energy spectra observed in the dissipation range of both kinetic numerical simulations and solar wind observations.

We use a large, statistical set of measurements from the Wind spacecraft at 1 AU, and supporting synthetic spacecraft data based on kinetic plasma theory, to show that the compressible component of inertial range solar wind turbulence is primarily in the kinetic slow mode. The zero-lag cross correlation C(delta n, delta B_parallel) between proton density fluctuations delta n and the field-aligned (compressible) component of the magnetic field delta B_parallel is negative and close to -1. The typical dependence of C(delta n,delta B_parallel) on the ion plasma beta_i is consistent with a spectrum of compressible wave energy that is almost entirely in the kinetic slow mode. This has important implications for both the nature of the density fluctuation spectrum and for the cascade of kinetic turbulence to short wavelengths, favoring evolution to the kinetic Alfven wave mode rather than the (fast) whistler mode.

This paper employs a recent turbulent heating prescription to predict the ratio of proton-to-total heating due to the kinetic dissipation of Alfvenic turbulence as a function of heliocentric distance. Comparing to a recent empirical estimate for this turbulent heating ratio in the high-speed solar wind, the prediction shows good agreement with the empirical estimate for R >~ 0.8 AU, but predicts less ion heating than the empirical estimate at smaller heliocentric radii. At these smaller radii, the turbulent heating prescription, calculated in the gyrokinetic limit, fails because the turbulent cascade is predicted to reach the proton cyclotron frequency before Landau damping terminates the cascade. These findings suggest that the turbulent cascade can reach the proton cyclotron frequency at R <~ 0.8 AU, leading to a higher level of proton heating than predicted by the turbulent heating prescription in the gyrokinetic limit. At larger heliocentric radii, R >~ 0.8 AU, this turbulent heating prescription contains all of the necessary physical mechanisms needed to reproduce the empirically estimated proton-to-total heating ratio.

The first three-dimensional, nonlinear gyrokinetic simulation of plasma turbulence resolving scales from the ion to electron gyroradius with a realistic mass ratio is presented, where all damping is provided by resolved physical mechanisms. The resulting energy spectra are quantitatively consistent with a magnetic power spectrum scaling of $k^{-2.8}$ as observed in \emph{in situ} spacecraft measurements of the "dissipation range" of solar wind turbulence. Despite the strongly nonlinear nature of the turbulence, the linear kinetic \Alfven wave mode quantitatively describes the polarization of the turbulent fluctuations. The collisional ion heating is measured at sub-ion-Larmor radius scales, which provides the first evidence of the ion entropy cascade in an electromagnetic turbulence simulation.

The gyrokinetic simulation code AstroGK is developed to study fundamental aspects of kinetic plasmas and for applications mainly to astrophysical problems. AstroGK is an Eulerian slab code that solves the electromagnetic Gyrokinetic-Maxwell equations in five-dimensional phase space, and is derived from the existing gyrokinetics code GS2 by removing magnetic geometry effects. Algorithms used in the code are described. The code is benchmarked using linear and nonlinear problems. Serial and parallel performance scalings are also presented.

The ratio of ion to electron heating due to the dissipation of Alfvenic turbulence in astrophysical plasmas is calculated based on a cascade model for turbulence in weakly collisional plasmas. Conditions for validity of this model are discussed, a prescription for the turbulent heating is presented, and it is applied to predict turbulent heating in accretion disks and the interstellar medium.

Two-dimensional electrostatic turbulence in magnetized weakly-collisional plasmas exhibits a cascade of entropy in phase space [Phys. Rev. Lett. 103, 015003 (2009)]. At scales smaller than the gyroradius, this cascade is characterized by the dimensionless ratio D of the collision time to the eddy turnover time measured at the scale of the thermal Larmor radius. When D >> 1, a broad spectrum of fluctuations at sub-Larmor scales is found in both position and velocity space. The distribution function develops structure as a function of v_{perp}, the velocity coordinate perpendicular to the local magnetic field. The cascade shows a local-scale nonlinear interaction in both position and velocity spaces, and Kolmogorov's scaling theory can be extended into phase space.

Measurements of small-scale turbulent fluctuations in the solar wind find a non-zero right-handed magnetic helicity. This has been interpreted as evidence for ion cyclotron damping. However, theoretical and empirical evidence suggests that the majority of the energy in solar wind turbulence resides in low frequency anisotropic kinetic Alfven wave fluctuations that are not subject to ion cyclotron damping. We demonstrate that a dissipation range comprised of kinetic Alfven waves also produces a net right-handed fluctuating magnetic helicity signature consistent with observations. Thus, the observed magnetic helicity signature does not necessarily imply that ion cyclotron damping is energetically important in the solar wind.

The proton temperature anisotropy in the solar wind is known to be constrained by the theoretical thresholds for pressure anisotropy-driven instabilities. Here we use approximately 1 million independent measurements of gyroscale magnetic fluctuations in the solar wind to show for the first time that these fluctuations are enhanced along the temperature anisotropy thresholds of the mirror, proton oblique firehose, and ion cyclotron instabilities. In addition, the measured magnetic compressibility is enhanced at high plasma beta ($\beta_\parallel \gtrsim 1$) along the mirror instability threshold but small elsewhere, consistent with expectations of the mirror mode. The power in this frequency (the 'dissipation') range is often considered to be driven by the solar wind turbulent cascade, an interpretation which should be qualified in light of the present results. In particular, we show that the short wavelength magnetic fluctuation power is a strong function of collisionality, which relaxes the temperature anisotropy away from the instability conditions and reduces correspondingly the fluctuation power.

One proposed mechanism for heating the solar wind, from close to the sun to beyond 10 AU, invokes low-frequency, oblique, Alfven-wave turbulence. Because small-scale oblique Alfven waves (kinetic Alfven waves) are compressive, the measured density fluctuations in the solar wind place an upper limit on the amplitude of kinetic Alfven waves and hence an upper limit on the rate at which the solar wind can be heated by low-frequency, Alfvenic turbulence. We evaluate this upper limit for both coronal holes at 5 solar radii and in the near-Earth solar wind. At both radii, the upper limit we find is consistent with models in which the solar wind is heated by low-frequency Alfvenic turbulence. At 1 AU, the upper limit on the turbulent heating rate derived from the measured density fluctuations is within a factor of 2 of the measured solar wind heating rate. Thus if low-frequency Alfvenic turbulence contributes to heating the near-Earth solar wind, kinetic Alfven waves must be one of the dominant sources of solar wind density fluctuations at frequencies of order 1 Hz. We also present a simple argument for why density fluctuation measurements do appear to rule out models in which the solar wind is heated by non-turbulent high-frequency waves ``sweeping'' through the ion-cyclotron resonance, but are compatible with heating by low-frequency Alfvenic turbulence.

Different results for the cascade power in strong, incompressible MHD turbulence with nonzero cross helicity appear in the literature. In this paper, we discuss the conditions under which these different results are valid. We define z+ to be the rms amplitude of Alfven waves propagating parallel to the background magnetic field, and z- to be the rms amplitude of Alfven waves propagating anti-parallel to the background magnetic field. Nonzero cross helicity implies that z+ and z- differ, and we take z- to be less than z+. We find that the mechanism that generates the z- fluctuations strongly affects the cascade power, because it controls the coherence time for interactions between oppositely directed wave packets at the outer scale. In particular, for fixed values of z+ and z-, the cascade power is in many cases larger when the z- fluctuations are generated by the reflection of z+ fluctuations than when the z- fluctuations are generated by forcing that is only weakly correlated with the z+ fluctuations.

We present a theoretical framework for plasma turbulence in astrophysical plasmas (solar wind, interstellar medium, galaxy clusters, accretion disks). The key assumptions are that the turbulence is anisotropic with respect to the mean magnetic field and frequencies are low compared to the ion cyclotron frequency. The energy injected at the outer scale scale has to be converted into heat, which ultimately cannot be done without collisions. A KINETIC CASCADE develops that brings the energy to collisional scales both in space and velocity. Its nature depends on the physics of plasma fluctuations. In each of the physically distinct scale ranges, the kinetic problem is systematically reduced to a more tractable set of equations. In the "inertial range" above the ion gyroscale, the kinetic cascade splits into a cascade of Alfvenic fluctuations, which are governed by the RMHD equations at both the collisional and collisionless scales, and a passive cascade of compressive fluctuations, which obey a linear kinetic equation along the moving field lines associated with the Alfvenic component. In the "dissipation range" between the ion and electron gyroscales, there are again two cascades: the kinetic-Alfven-wave (KAW) cascade governed by two fluid-like Electron RMHD equations and a passive phase-space cascade of ion entropy fluctuations. The latter cascade brings the energy of the inertial-range fluctuations that was damped by collisionless wave-particle interaction at the ion gyroscale to collisional scales in the phase space and leads to ion heating. The KAW energy is similarly damped at the electron gyroscale and converted into electron heat. Kolmogorov-style scaling relations are derived for these cascades. Astrophysical and space-physical applications are discussed in detail.

The limitations of Hall MHD as a model for turbulence in weakly collisional plasmas are explored using quantitative comparisons to Vlasov-Maxwell kinetic theory over a wide range of parameter space. The validity of Hall MHD in the cold ion limit is shown, but spurious undamped wave modes exist in Hall MHD when the ion temperature is finite. It is argued that turbulence in the dissipation range of the solar wind must be one, or a mixture, of three electromagnetic wave modes: the parallel whistler, oblique whistler, or kinetic Alfven waves. These modes are generally well described by Hall MHD. Determining the applicability of linear kinetic damping rates in turbulent plasmas requires a suite of fluid and kinetic nonlinear numerical simulations. Contrasting fluid and kinetic simulations will also shed light on whether the presence of spurious wave modes alters the nonlinear couplings inherent in turbulence and will illuminate the turbulent dynamics and energy transfer in the regime of the characteristic ion kinetic scales.

Howes et al. Reply to Comment on "Kinetic Simulations of Magnetized Turbulence in Astrophysical Plasmas" arXiv:0711.4355

This paper studies the turbulent cascade of magnetic energy in weakly collisional magnetized plasmas. A cascade model is presented, based on the assumptions of local nonlinear energy transfer in wavenumber space, critical balance between linear propagation and nonlinear interaction times, and the applicability of linear dissipation rates for the nonlinearly turbulent plasma. The model follows the nonlinear cascade of energy from the driving scale in the MHD regime, through the transition at the ion Larmor radius into the kinetic Alfven wave regime, in which the turbulence is dissipated by kinetic processes. The turbulent fluctuations remain at frequencies below the ion cyclotron frequency due to the strong anisotropy of the turbulent fluctuations, k_parallel << k_perp (implied by critical balance). In this limit, the turbulence is optimally described by gyrokinetics; it is shown that the gyrokinetic approximation is well satisfied for typical slow solar wind parameters. Wave phase velocity measurements are consistent with a kinetic Alfven wave cascade and not the onset of ion cyclotron damping. The conditions under which the gyrokinetic cascade reaches the ion cyclotron frequency are established. Cascade model solutions imply that collisionless damping provides a natural explanation for the observed range of spectral indices in the dissipation range of the solar wind. The dissipation range spectrum is predicted to be an exponential fall off; the power-law behavior apparent in observations may be an artifact of limited instrumental sensitivity. The cascade model is motivated by a programme of gyrokinetic simulations of turbulence and particle heating in the solar wind.

This paper describes a conceptual framework for understanding kinetic plasma turbulence as a generalized form of energy cascade in phase space. It is emphasized that conversion of turbulent energy into thermodynamic heat is only achievable in the presence of some (however small) degree of collisionality. The smallness of the collision rate is compensated by the emergence of small-scale structure in the velocity space. For gyrokinetic turbulence, a nonlinear perpendicular phase mixing mechanism is identified and described as a turbulent cascade of entropy fluctuations simultaneously occurring at spatial scales smaller than the ion gyroscale and in velocity space. Scaling relations for the resulting fluctuation spectra are derived. An estimate for the collisional cutoff is provided. The importance of adequately modeling and resolving collisions in gyrokinetic simulations is biefly discussed, as well as the relevance of these results to understanding the dissipation-range turbulence in the solar wind and the electrostatic microturbulence in fusion plasmas.

This letter presents the first ab initio, fully electromagnetic, kinetic simulations of magnetized turbulence in a homogeneous, weakly collisional plasma at the scale of the ion Larmor radius (ion gyroscale). Magnetic and electric-field energy spectra show a break at the ion gyroscale; the spectral slopes are consistent with scaling predictions for critically balanced turbulence of Alfven waves above the ion gyroscale (spectral index -5/3) and of kinetic Alfven waves below the ion gyroscale (spectral indices of -7/3 for magnetic and -1/3 for electric fluctuations). This behavior is also qualitatively consistent with in situ measurements of turbulence in the solar wind. Our findings support the hypothesis that the frequencies of turbulent fluctuations in the solar wind remain well below the ion cyclotron frequency both above and below the ion gyroscale.

The transfer of turbulent energy through an inertial range from the driving scale to dissipative scales in a kinetic plasma followed by the conversion of this energy into heat is a fundamental plasma physics process. A theoretical foundation for the study of this process is constructed, but the details of the kinetic cascade are not well understood. Several important properties are identified: (a) the conservation of a generalized energy by the cascade; (b) the need for collisions to increase entropy and realize irreversible plasma heating; and (c) the key role played by the entropy cascade--a dual cascade of energy to small scales in both physical and velocity space--to convert ultimately the turbulent energy into heat. A strategy for nonlinear numerical simulations of kinetic turbulence is outlined. Initial numerical results are consistent with the operation of the entropy cascade. Inertial range turbulence arises in a broad range of space and astrophysical plasmas and may play an important role in the thermalization of fusion energy in burning plasmas.

We present a cascade model for turbulence in weakly collisional plasmas that follows the nonlinear cascade of energy from the large scales of driving in the MHD regime to the small scales of the kinetic Alfven wave regime where the turbulence is dissipated by kinetic processes. Steady-state solutions of the model for the slow solar wind yield three conclusions: (1) beyond the observed break in the magnetic energy spectrum, one expects an exponential cut-off; (2) the widely held interpretation that this dissipation range obeys power-law behavior is an artifact of instrumental sensitivity limitations; and, (3) over the range of parameters relevant to the solar wind, the observed variation of dissipation range spectral indices from -2 to -4 is naturally explained by the varying effectiveness of Landau damping, from an undamped prediction of -7/3 to a strongly damped index around -4.

Magnetohydrodynamic (MHD) turbulence is encountered in a wide variety of astrophysical plasmas, including accretion disks, the solar wind, and the interstellar and intracluster medium. On small scales, this turbulence is often expected to consist of highly anisotropic fluctuations with frequencies small compared to the ion cyclotron frequency. For a number of applications, the small scales are also collisionless, so a kinetic treatment of the turbulence is necessary. We show that this anisotropic turbulence is well described by a low frequency expansion of the kinetic theory called gyrokinetics. This paper is the first in a series to examine turbulent astrophysical plasmas in the gyrokinetic limit. We derive and explain the nonlinear gyrokinetic equations and explore the linear properties of gyrokinetics as a prelude to nonlinear simulations. The linear dispersion relation for gyrokinetics is obtained and its solutions are compared to those of hot-plasma kinetic theory. These results are used to validate the performance of the gyrokinetic simulation code {\tt GS2} in the parameter regimes relevant for astrophysical plasmas. New results on global energy conservation in gyrokinetics are also derived. We briefly outline several of the problems to be addressed by future nonlinear simulations, including particle heating by turbulence in hot accretion flows and in the solar wind, the magnetic and electric field power spectra in the solar wind, and the origin of small-scale density fluctuations in the interstellar medium.

We introduce Gradient Particle Magnetohydrodynamics (GPM), a new Lagrangian method for magnetohydrodynamics based on gradients corrected for the locally disordered particle distribution. The development of a numerical code for MHD simulation using the GPM algorithm is outlined. Validation tests simulating linear and nonlinear sound waves, linear MHD waves, advection of magnetic fields in a magnetized vortex, hydrodynamical shocks, and three-dimensional collapse are presented, demonstrating the viability of an MHD code using GPM. The characteristics of a GPM code are discussed and possible avenues for further development and refinement are mentioned. We conclude with a view of how GPM may complement other methods currently in development for the next generation of computational astrophysics.

The effect of magnetic shear and shear flow on local gravitationally induced instabilities is investigated. A simple model is constructed allowing for an arbitrary entropy gradient and a shear plasma flow in the Boussinesq approximation. A transformation to shearing magnetic coordinates achieves a model with plasma flow along the magnetic field lines where the coordinate lines are coincident with the field lines. The solution for the normal modes of the system depends on two parameters: the Alfven Mach number of the plasma flow and the entropy gradient. The behavior of the unstable normal modes of this system is summarized by a stability diagram. Important characteristics of this stability diagram are the following: magnetic shear is stabilizing and the entropy gradient must exceed a threshold value for unstable mode growth to occur; flow acts to suppress mode growth in a substantially unstable regime as expected, yet near marginal stability it can lessen the stabilizing effect of magnetic shear and enhance the growth rates of the instability; and, as the Alfven Mach number approaches one, the instability is completely stabilized. Analytical work is presented supporting the characteristics of the stability diagram and illuminating the physical mechanisms controlling the behavior of the model. The implications of this work for astrophysical and fusion applications and the potential for future research extending the results to include compressibility are discussed.