Search

Whistler waves propagating nearly parallel to the ambient magnetic field experience a nonlinear instability due to transverse currents when the background plasma has a population of sufficiently low energy electrons. Intriguingly, this nonlinear process may generate oblique electrostatic waves, including whistlers near the resonance cone with properties resembling oblique chorus waves in the Earth’s magnetosphere. Focusing on the generation of oblique whistlers, earlier analysis of the instability is extended here to the case where low-energy background plasma consists of both a “cold” population with energy of a few eV and a “warm” electron component with energy of the order of 100 eV. This is motivated by spacecraft observations in the Earth’s magnetosphere where oblique chorus waves were shown to interact resonantly with the warm electrons. The main new results are: 1) the instability producing oblique electrostatic waves is sensitive to the shape of the electron distribution at low energies. In the whistler range of frequencies, two distinct peaks in the growth rate are typically present for the model considered: a peak associated with the warm electron population at relatively low wavenumbers and a peak associated with the cold electron population at relatively high wavenumbers; 2) overall, the instability producing oblique whistler waves near the resonance cone persists (with a reduced growth rate) even in the cases where the temperature of the cold population is relatively high, including cases where cold population is absent and only the warm population is included; 3) particle-in-cell simulations show that the instability leads to heating of the background plasma and formation of characteristic plateau and beam features in the parallel electron distribution function in the range of energies resonant with the instability. The plateau/beam features have been previously detected in spacecraft observations of oblique chorus waves. However, they have been attributed to external sources and have been proposed to be the mechanism generating oblique chorus. In the present scenario, the causality link is reversed and the instability generating oblique whistler waves is shown to be a possible mechanism for formation of the plateau and beam features.

Understanding plasma turbulence requires a synthesis of experiments, observations, theory, and simulations. In the case of kinetic plasmas such as the solar wind, the lack of collisions renders the fluid closures such as viscosity meaningless and one needs to resort to higher-order fluid models or kinetic models. Typically, the computational expense in such models is managed by simulating artificial values of certain parameters such as the ratio of the Alfvén speed to the speed of light ( v _A / c ) or the relative mass ratio of ions and electrons ( m _i / m _e ). Although, typically care is taken to use values as close as possible to realistic values within the computational constraints, these artificial values could potentially introduce unphysical effects. These unphysical effects could be significant at sub-ion scales, where kinetic effects are the most important. In this paper, we use the 10-moment fluid model in the Gkeyll framework to perform controlled numerical experiments, systematically varying the ion–electron mass ratio from a small value down to the realistic proton–electron mass ratio. We show that the unphysical mass ratio has a significant effect on the kinetic range dynamics as well as the heating of both plasma species. The dissipative process for both ions and electrons becomes more compressive in nature, although the ions remain nearly incompressible in all cases. The electrons move from being dominated by incompressive viscous-like heating/dissipation to very compressive heating/dissipation dominated by compressions/rarefactions. While the heating change is significant for the electrons, a mass ratio of m _i / m _e ∼ 250 captures the asymptotic behavior of electron heating.

The interaction of β -particles with the weakly ionized plasma background is an important mechanism for powering the kilonova (KN) transient signal from neutron star mergers. For this purpose, we present an implementation of the approximate fast-particle collision kernel, described by Inokuti following the seminal formulation of Bethe, in a spectral solver of the Vlasov–Maxwell–Boltzmann equation. In particular, we expand the fast-particle plane-wave atomic excitation kernel into coefficients of the Hermite basis, and derive the relevant discrete spectral system. In this fast-particle limit, the approach permits the direct use of atomic data, including optical oscillator strengths, normally applied to photon–matter interaction. The resulting spectral matrix is implemented in the MASS-APP spectral solver framework, in a way that avoids full matrix storage per spatial zone. We numerically verify aspects of the matrix construction, and present a proof-of-principle 3D simulation of a 2D axisymmetric KN ejecta snapshot. Our preliminary numerical results indicate that a reasonable choice of Hermite basis parameters for β -particles in the KN is a bulk velocity parameter u = 0, a thermal velocity parameter α = 0.5 c , and a 9 × 9 × 9 mode velocity basis set (Hermite orders of 0–8 in each dimension). For interior-ejecta sample zones, we estimate that the ratio of thermalization from large-angle (≳2.°5) bound excitation scattering to total thermalization is ∼0.002–0.003.

In the Earth’s dipole magnetosphere finite-gyroradius effects produce a shift of the atmospheric loss cone away from the direction of the magnetic field. This loss-cone shift is theoretically described by the “Mozer transform” [Mozer, F. S. (1966). Proton trajectories in the radiation belt. J. Geophys. Res. 71:2701], which is based upon the curvature drift of particles crossing the equatorial plane. For positive ions the northern and southern loss cones both shift westward and for electrons the northern and southern loss cones both shift eastward. This loss-cone shift is part of a coordinate-system transform, with the transformed coordinates better organizing the behavior of particle orbits in the dipole magnetic field (e.g. first adiabatic invariants, mirror heights, and bounce times). In this report it is demonstrated that the transformed coordinate system also properly organizes pitch-angle diffusion. This improved organization of the diffusion is true whether the angular scattering is produced by plasma-wave scattering or by field-line-curvature (FLC) scattering. It is shown that FLC scattering and the loss cone shift are linked, so that if FLC scattering is occurring, there is a loss cone shifted away from the magnetic-field direction and the Mozer-transformed coordinates are needed.

This perspective advocates for the need of a combined system science approach to global magnetospheric models and to spacecraft magnetospheric data to answer the question “Do simulations behave in the same manner as the magnetosphere does?” (instead of the standard validation question “How well do simulations reproduce spacecraft data?”). This approach will 1) validate global magnetospheric models statistically, without the need for a direct comparison against spacecraft data, 2) expose the deficiencies of the models, and 3) provide physics support to the system analysis performed on the magnetospheric system.

When the Earth’s moon is in the supersonic solar wind, the darkside of the Moon and the lunar plasma wake can be very dangerous charging environments. In the absence of photoelectron emission (dark) and in the absence of cool plasma (wake), the emission or collection of charge to reduce electrical potentials is difficult. Unique extreme charging events may occur during impulsive solar-energetic-electron (SEE) events when the lunar wake is dominated by relativistic electrons, with the potential to charge and differentially charge objects on and above the lunar surface to very-high negative electrical potentials. In this report the geometry of the magnetic connections from the Sun to the lunar nightside are explored; these magnetic connections are the pathways for SEEs from the Sun. Rudimentary charging calculations for objects in the relativistic-electron environment of the lunar wake are performed. To enable these charging calculations, secondary-electron yields for impacts by relativistic electrons are derived. Needed lunar electrical-grounding precautions for SEE events are discussed. Calls are made 1) for future dynamic simulations of the plasma wake in the presence of time-varying SEE-event relativistic electrons and time-varying solar-wind magnetic-field orientations and 2) for future charging calculations in the relativistic-electron wake environment and on the darkside lunar surface.

Insufficiently accurate magnetic-field-line mapping between the aurora and the equatorial magnetosphere prevents us from determining the cause of many types of aurora. An important example is the longstanding question of how the magnetosphere drives low-latitude (growth-phase) auroral arcs: a large number of diverse generator mechanisms have been hypothesized but equatorial magnetospheric measurements cannot be unambiguously connected to arcs in the ionosphere, preventing the community from identifying the correct generator mechanisms. Here a mission concept is described to solve the magnetic-connection problem. From an equatorial instrumented spacecraft, a powerful energetic-electron beam is fired into the atmospheric loss cone resulting in an optical beam spot in the upper atmosphere that can be optically imaged from the ground, putting the magnetic connection of the equatorial spacecraft’s measurements into the context of the aurora. Multiple technical challenges that must be overcome for this mission concept are discussed: these include spacecraft charging, beam dynamics, beam stability, detection of the beam spot in the presence of aurora, and the safety of nearby spacecraft.

The 1970s and 1980s were heydays for using active electron beam experiments to probe some of the fundamental physical processes that occur throughout the heliosphere and in astrophysical contexts. Electron beam experiments were used to study spacecraft charging and spacecraft-plasma coupling; beam-plasma interaction physics; magnetic bounce and drift physics; auroral physics; wave generation; and military applications. While these experiments were enormously successful, they were also limited by the technologies that were available at that time. New advances in space instrumentation, data collection, and accelerator technologies enable a revolutionary new generation of active experiments using electron beams in space. In this paper we discuss such an experiment, the Beam Plasma Interactions Experiment (Beam PIE), a sounding rocket experiment designed to (a) advance high-electron mobility transistor-based radio frequency (RF) linear accelerator electron technology for space applications and (b) study the production of whistler and X-mode waves by modulated electron beams.

Planned active space experiments and ideas for future active space experiments are reviewed. Three active experiments being readied are DSX (Demonstration and Space eXperiments), SMART (Space Measurement of Rocket-released Turbulence), and BeamPIE (Beam Plasma Interaction Experiment). Ideas for future experiments include relativistic-electron-beam experiments for magnetic-field-line tracing, relativistic-electron-beam experiments to probe the middle atmosphere, plasma-wave launching using superparamagnetic-nanoparticle amplification of magnetic fields, the heavy-ion mass loading of collisionless magnetic-field-line reconnection, the use of electrostatically charged tethers to pitch-angle scatter radiation-belt particles, cold plasma releases to modify magnetospheric plasma physics, and neutral-gas releases to enhance neutral-particle imaging of the magnetosphere. Technologies that are being developed to enable future space active experiments are reviewed: this includes the development of compact relativistic accelerators, superparamagnetic particle amplified antennae, CubeSats, and a new understanding of how to control dynamic spacecraft charging. New capabilities to use laboratory facilities to design space active experiments as well as new computer-simulation capabilities to design and understand space active experiments are reviewed.

Energetic electron beams have been proposed for tracing magnetic field lines from the magnetosphere down to the ionosphere, in active experiments aimed at diagnosing mechanisms at play in the coupling between magnetosphere and ionosphere. It is recognized however that in the absence of an efficient mitigation technique, this approach would lead to unacceptably large spacecraft charging and positive potential buildup, which would result in environmental hazard for the spacecraft. This problem would be particularly acute in low density regions of the magnetosphere of interest in the study of magnetic field reconnection and substorm dynamics. A solution to this predicament could consist of creating a plasma contactor whereby a gas puff would be ionized, leading to the evacuation of positive charges and collection of cold electrons, thus compensating for the charges lost in the electron beam. A possible alternative is presented here, which consists of attaching a large passive conducting surface to the spacecraft, a “tethered capacitor”, from which negative charges would be drawn to compensate for those lost from the beam. This capacitor would then charge to a large positive potential, leaving the spacecraft and electron gun at a lower, acceptable positive potential. The tethered capacitor could have a relatively small mass; consisting only of a thin conducting surface that would be “inflated” as a result of repulsive electrostatic forces. This charge mitigation concept, as applied to active electron beam experiments, is explored using three dimensional particle-in-cell (PIC) simulations from which scaling laws can be inferred for the spacecraft and tethered capacitor potentials under proposed electron beam operations.

Spectral (transform) methods for solution of Vlasov-Maxwell system have shown significant promise as numerical methods capable of efficiently treating fluid-kinetic coupling in magnetized plasmas. We discuss SpectralPlasmaSolver (SPS), an implementation of three-dimensional, fully electromagnetic algorithm based on a decomposition of the plasma distribution function in Hermite modes in velocity space and Fourier modes in physical space. A fully-implicit time discretization is adopted for numerical stability and to ensure exact conservation laws for total mass, momentum and energy. The SPS code is parallelized using Message Passing Interface for distributed memory architectures. Application of the method to analysis of kinetic range of scales in plasma turbulence under conditions typical of the solar wind is demonstrated. With only 4 Hermite modes per velocity dimension, the algorithm yields damping rates of kinetic Alfvén waves with accuracy of 50% or better, which is sufficient to obtain a model of kinetic scales capable of reproducing many of the expected statistical properties of turbulent fluctuations. With increasing number of Hermite modes, progressively more accurate values for collisionless damping rates are obtained. Fully nonlinear simulations of decaying turbulence are presented and successfully compared with similar simulations performed using Particle-In-Cell method.

Whistler waves propagating nearly parallel to the ambient magnetic field experience a nonlinear instability due to transverse currents when the background plasma has a population of sufficiently low energy electrons. Intriguingly, this nonlinear process may generate oblique electrostatic waves, including whistlers near the resonance cone with properties resembling oblique chorus waves in the Earth's magnetosphere. Focusing on the generation of oblique whistlers, earlier analysis of the instability is extended here to the case where low-energy background plasma consists of both a "cold" population with energy of a few eV and a "warm" electron component with energy of the order of 100 eV. This is motivated by spacecraft observations in the Earth's magnetosphere where oblique chorus waves were shown to interact resonantly with the warm electrons. The main new results are: 1) the instability producing oblique electrostatic waves is sensitive to the shape of the electron distribution at low energies. In the whistler range of frequencies, two distinct peaks in the growth rate are typically present for the model considered: a peak associated with the warm electron population at relatively low wavenumbers and a peak associated with the cold electron population at relatively high wavenumbers; 2) overall, the instability producing oblique whistler waves near the resonance cone persists (with a reduced growth rate) even in the cases where the temperature of the cold population is relatively high, including cases where cold population is absent and only the warm population is included; 3) particle-in-cell simulations show that the instability leads to heating of the background plasma and formation of characteristic plateau and beam features in the parallel electron distribution function in the range of energies resonant with the instability. The plateau/beam features have been previously detected in spacecraft observations of oblique chorus waves. However, they have been attributed to external sources and have been proposed to be the mechanism generating oblique chorus. In the present scenario, the causality link is reversed and the instability generating oblique whistler waves is shown to be a possible mechanism for formation of the plateau and beam features. [ABSTRACT FROM AUTHOR]

Recent observational and analytical studies suggested that a new regime of kinetic turbulence may exist in plasma environments with low electron beta (Chen & Boldyrev 2017). Such a regime, termed inertial kinetic-Alfvén turbulence (iKAW), is relevant for the solar corona, Earth’s magnetosheath, and other astrophysical systems where the electron and ion plasma beta parameters satisfy the condition βe ≪ βi ≲ 1. In this paper we present kinetic numerical simulations that confirm the existence of the iKAW regime. Specifically, the simulations demonstrate a transition at scales below electron inertial length de when βe ≪ βi ≲ 1. Spectral slopes and other statistical properties of turbulence at sub-de scales are consistent with the phenomenological theory of inertial kinetic-Alfvén turbulence proposed by Chen & Boldyrev (2017) and with the recent observations in the Earth’s magnetosheath. [ABSTRACT FROM AUTHOR]

Boundaries between space plasmas occur in numerous contexts and scales, from astrophysical jets to planetary magnetospheres. Mass and momentum transport across boundaries poses a fundamental problem in magnetospheric physics. Kelvin-Helmholtz instability (KHI) is a promising mechanism to facilitate transport. Although previous studies have suggested KHI occurrence in various space plasmas, theory predicts that compressibility prevents KHI excitation at boundaries with large density gradients because of previously considered boundary structures where density varies with velocity. Based on the observations of a large density gradient boundary by MAVEN at Mars, where we can observe an extreme case, in this study, we show that it is the entropy, instead of the previously considered density, that varies with the velocity in the real velocity-sheared boundary. The entropy-based boundary structure places the velocity shear in a lower-density region than the traditional density-based structure and weakens the compressibility effect. This new boundary structure thus enables KHI excitation even at large density gradient boundaries, such as at the ionopause of unmagnetized planets and the plasmapause of magnetized planets. The result suggests the ubiquitous occurrence of KHI in the plasma universe and emphasizes its important role in planetary cold plasma escape from unmagnetized planets. [ABSTRACT FROM AUTHOR]

Hollow cathode plasma contactors have an extensive history as spacecraft neutralizers, but questions remain on whether the technology can neutralize significant positive charging in tenuous space plasmas. Simulations and a representative semianalytical model (which we term the ion emission model) predict effective neutralization in these scenarios as ion current is emitted from the surface of the quasi-neutral contactor plasma according to the space-charge limit. This article focuses on experimental plasma measurements of a hollow cathode biased to several positive voltages with respect to a surrounding vacuum chamber in order to validate the ion emission model. Particular attention is paid to the emitted ions (those reaching the chamber wall) and to the ion-rich sheath region where the ion emission process takes place. Retarding potential analyzer (RPA) measurements were performed to understand ion flow velocity and direction as these parameters relate directly to the current predicted by the ion emission model. Planar probe measurements were performed at the chamber wall to compare local emission currents to model predictions. Evidence of collisions within the plasma (particularly charge exchange collisions) and simple models predicting ion energetics are presented. These experiments suggest that ion emission from the contactor plasma acts according to the ion emission model and adds to the physical understanding of spacecraft neutralization using hollow cathode plasma contactors as it may occur in tenuous space plasmas. [ABSTRACT FROM AUTHOR]

Ganguli and Crabtree have written a comment about a recent article by the authors listed above on radiation-belt remediation. They have objected to our evaluation of the Naval Research Laboratory’s chemical release concept which states that this concept may be impractical due to an apparently low overall efficiency. In their comment, they provide a scientific argument and refer to the published literature to counter our statement. Here, we provide more details on our numerical calculations and experimental results which led to this evaluation. [ABSTRACT FROM AUTHOR]

Energetic electrons can be trapped in Earth’s magnetic field, forming the radiation belts (also known as the Van Allen Belts). These electrons, which can originate from the solar wind or a high-altitude nuclear explosion (HANE), have the potential to damage satellites in low-Earth orbit (LEO). For example, in 1962, the U.S. detonated a nuclear device at an altitude of about 400 km in the Starfish experiment. The resulting enhancement of the radiation belts disabled several satellites within a few months and energetic electrons remained in the radiation belts for up to several years. In order to address this potential vulnerability, schemes have been proposed to drain electrons from the radiation belts, with the most promising approaches based on using high-power very-low-frequency (VLF) waves to scatter the electrons into more field-aligned trajectories, forcing them to precipitate into Earth’s atmosphere. This paper will provide an overview of enhanced electron distributions in the radiation belts as well as approaches to VLF wave belt remediation including the use of either antennas or relativistic electrons beams in space to generate the VLF waves. [ABSTRACT FROM AUTHOR]

We describe a new electrostatic particle-in-cell (PIC) code in curvilinear geometry called curvilinear PIC (CPIC). The code models the microscopic (kinetic) evolution of a plasma with the PIC method, coupled with an adaptive computational grid that can conform to arbitrarily shaped domains. CPIC is particularly suited for multiscale problems associated with the interaction of complex objects with plasmas. A map is introduced between the physical space and the logical space, where the grid is uniform and Cartesian. In CPIC, most of the operations are performed in logical space. CPIC was designed following criteria of versatility, robustness, and performance. Its main features are the use of structured meshes, a scalable field solver based on the black box multigrid algorithm and a hybrid mover, where particles' position is in logical space while the velocity is in physical space. Test examples involving the interaction of a plasma with material boundaries are presented. [ABSTRACT FROM PUBLISHER]

We present simulations of the nonlinear dynamics of a nonneutral plasma in a Penning trap. We use KANDINSKY, an advanced cylindrical particle-in-cell (PlC) code based on a new method that avoids the traditional problems of cylindrical PlC codes. The principal novelty resides in defining a particle volume in addition to the usual particle charge. Results from the two-dimensional drift-Poisson model and a recently developed model with compression effects are presented. [ABSTRACT FROM AUTHOR]