The kinetic stability of collisionless, sloshing beam-ion ( $45^\circ$ pitch angle) plasma is studied in a three-dimensional (3-D) simple magnetic mirror, mimicking the Wisconsin high-temperature superconductor axisymmetric mirror experiment. The collisional Fokker–Planck code CQL3D-m provides a slowing-down beam-ion distribution to initialize the kinetic-ion/fluid-electron code Hybrid-VPIC, which then simulates free plasma decay without external heating or fuelling. Over $1$ – $10\;\mathrm{\unicode{x03BC} s}$ , drift-cyclotron loss-cone (DCLC) modes grow and saturate in amplitude. The DCLC scatters ions to a marginally stable distribution with gas-dynamic rather than classical-mirror confinement. Sloshing ions can trap cool (low-energy) ions in an electrostatic potential well to stabilize DCLC, but DCLC itself does not scatter sloshing beam-ions into the said well. Instead, cool ions must come from external sources such as charge-exchange collisions with a low-density neutral population. Manually adding cool $\mathord {\sim } 1\;\mathrm{keV}$ ions improves beam-ion confinement several-fold in Hybrid-VPIC simulations, which qualitatively corroborates prior measurements from real mirror devices with sloshing ions.

The resonant right-hand instability (RHI) is often the dominant mode driven by reflected ions upstream of Earth’s quasi-parallel bow shock. In the tradition of Peter Gary, this paper further explores the right-hand instability using numerical solutions of the plasma dispersion relation and non-linear kinetic simulations, with parameters inspired by observations from NASA’s Magnetospheric Multiscale (MMS) mission. Agreement is found between the ion distributions in the particle-in-cell simulations and Magnetospheric Multiscale spacecraft data, which show the gyrophase bunching characteristic of the instability. The non-linear structures created by right-hand instability tend to be stronger when the plasma beta is lower. These structures have sizes of around 100 to 200 ion inertial lengths perpendicular to the magnetic field, presenting planet-sized disturbances to the magnetosphere. 2d and 3D hybrid particle-in-cell simulations show that modes with a range of propagation angles oblique to the magnetic field are excited, providing a ground to understand previous statistical studies of observed foreshock waves.

Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.

Abstract Magnetospheric Multiscale observed an extended current layer concurrent with a strong parallel electric field. The current layer is embedded in the reconnection exhaust and is consistent with an extended magnetized jet expected in roughly symmetric moderate guide field reconnection. The strong parallel electric field observed at the jet balances a strong gradient in the parallel electron pressure caused by a transition between magnetized electrons with Te‖ > Te⊥ on one side of the jet and demagnetized, roughly isotropic electrons across the thin current layer. Simulation results show how this transition can occur and how the associated electron pressure gradients in Ohm's law balance the parallel electric field.

made in both visible light (2000-6000 Angstroms and x-radiation (10-100 keV). Peak current densities reached approx. 3kAmp/sq. cm. at approx. 20 cm from the diode where the beam appeared to consist of two separate coaxial components both of near uniform current density. The central component contained approx. 65% of the total beam current and was initially within a radius of approx. 1.4 cm. The outer component of the beam started with a radius of approx. 5 cm (defined by the anode foil) and appeared to expand steadily. The visible emissions belonged to the species N2 and N2+ and were synchronous with the electron beam current. The absolute emission in the band N2+ 1N(0-0) agreed with the calculated fluorescence intensity within experimental uncertainty, but the emission in the band N(2)2P(0-0) was approx 6 times larger than suggested by simple fluorescence calculations. A simple technique was devised to select a fraction of the electron beam emerging from the diode to produce a pinched matched beam.

Various systems can be modeled as a point-like explosion of ionized debris into a magnetized, collisionless background plasma – including astrophysical examples, active experiments in space, and laser-driven laboratory experiments. Debris streaming from the explosion parallel to the magnetic field may drive multiple resonant and non-resonant ion-ion beam instabilities, some of which can efficiently couple the debris energy to the background and may even support the formation of shocks. We present a large-scale hybrid (kinetic ions + fluid electrons) particle-in-cell (PIC) simulation, extending hundreds of ion inertial lengths from a 3-D explosion, that resolves these instabilities. We show that the character of these instabilities differs notably from the 1-D equivalent by the presence of unique transverse structure. Additional 2-D simulations explore how the debris beam length, width, density, and speed affect debris-background coupling, with implications for the generation of quasi-parallel shocks.

Anisotropic electron heating predictions of the equations of state of Le et al. (2009); https://doi.org/10.1103/PhysRevLett.102.085001 (“the EoS”) are applied to reconnection in the Earth's magnetosheath. The model is applicable to open systems where electrons are streaming along magnetic field lines into the reconnection region, sourced by fixed external reservoirs of plasma ambient to the reconnection region. While spacecraft observations have previously shown the EoS to hold in the region near the X ‐line, we find that for an event observed far downstream (∼100 d i ) from the X ‐line both inflows and the exhaust follow the predictions of the EoS. Furthermore, the model underlying the EoS is extended to include additional perpendicular heating terms relevant to the considered event undergoing active magnetic compression with local electron trapping. The agreement between the spacecraft observations and the EoS emphasizes the roles of the ambient plasma sources and the dynamics of trapped and passing electrons in setting and controlling the anisotropic electron heating at large scale in naturally occurring plasma configurations.

We explore the dynamic magnetosphere of Mercury by employing a three‐dimensional hybrid particle-in-cell (particle ions and massless fluid electrons) code – hybrid-VPIC. The newly developed hybrid-VPIC code (based on the high-performance fully kinetic Vector Particle-In-Cell, VPIC code) incorporates ion kinetics (beam and anisotropy driven instabilities) that are critical for foreshock and magnetosheath physics, as well as the Hall effect which is important for collisionless magnetic reconnection; therefore, it is particularly well suited for investigating the kinetic physics of Mercury's dynamic magnetosphere. The simulation results are in good agreement with MESSENGER’s magnetic field measurements during its second Mercury flyby. We will investigate collisionless magnetic reconnection (including flux transfer events or FTEs and ion velocity distribution functions) and foreshock physics (including plasma turbulence and particle acceleration) in this study.

The coupling between a rapidly expanding cloud of ionized debris and an ambient magnetized plasma is revisited with a hybrid (kinetic ion/fluid electron) simulation code that allows a study over a wide range of plasma parameters. Over a specified range of hypothetical conditions, simple scaling laws in terms of the total debris mass and explosion speed are derived and verified for the maximal size of the debris cloud and the fraction of debris that free-streams from the burst along the magnetic field. The amount of debris that escapes from the burst with minimal coupling to the background magnetic field increases with the debris gyroradius. Test cases with two different debris species--including a heavy minority species with a relatively large gyroradius--highlight how the collisionless coupling of the debris depends on the single-particle trajectories as well as the overall conservation of energy and momentum.

We present a drift kinetic model for the free expansion of a thermal plasma out of a magnetic nozzle. This problem relates to plasma space propulsion systems, natural environments such as the solar wind, and end losses from the expander region of mirror magnetically confined fusion concepts such as the gas dynamic trap. The model incorporates trapped and passing orbit types encountered in the mirror expander geometry and maps to an upstream thermal distribution. This boundary condition and quasineutrality require the generation of an ambipolar potential drop of ∼5Te/e, forming a thermal barrier for the electrons. The model for the electron and ion velocity distributions and fluid moments is confirmed with data from a fully kinetic simulation. Finally, the model is extended to account for a population of fast sloshing ions arising from neutral beam heating within a magnetic mirror, again resulting in good agreement with a corresponding kinetic simulation.

In this paper, we derive a new drift-kinetic method for estimating gradients in the plasma properties through a velocity space distribution at a single point. The gradients are intrinsically related to agyrotropic features of the distribution function. This method predicts the gradients in the magnetized distribution function and can predict gradients of arbitrary moments of the gyrotropic background distribution function. The method allows for estimates on density and pressure gradients on the scale of a Larmor radius, proving to resolve smaller scales than any method currently available to spacecraft. The model is verified with a set of fully kinetic VPIC particle-in-cell simulations.

A valid fluid model for electrons in collisionless space plasmas is desirable for understanding the structure and evolution of magnetic reconnection geometries. Additionally, such a fluid model would be useful for the simulation of systems too large to be tractable in a fully kinetic model. Using Magnetospheric Multiscale spacecraft observations, we provide direct confirmation of the Lê et al., 2009 (https://doi.org/10.1103/PhysRevLett.102.085001) equations of state for the electron pressure tensor during guide field reconnection and demonstrate how the closure can be applied in efficient numerical simulation, yielding new physical insight to the electron heating problem. Furthermore, we use the Lê et al., 2009 equations of state to predict a scaling of electron heating in the exhaust comparable to the available observational data.

Aided by fully kinetic simulations, spacecraft observations of magnetic reconnection in Earth's magnetotail are analyzed. The structure of the electron diffusion region is in quantitative agreement with the numerical model. Of special interest, the spacecraft data reveal how reconnection is mediated by off-diagonal stress in the electron pressure tensor breaking the frozen-in law of the electron fluid.

Fully kinetic simulations of asymmetric magnetic reconnection reveal the presence of magnetic-field-aligned beams of electrons flowing toward the topological magnetic x line. Within the ∼6d_{e} electron-diffusion region, the beams become oblique to the local magnetic field, providing a unique signature of the electron-diffusion region where the electron frozen-in law is broken. The numerical predictions are confirmed by in situ Magnetospheric Multiscale spacecraft observations during asymmetric reconnection at Earth's dayside magnetopause.

The recent study of signal interference circuits, which find its origins in earlier work related to active channelized filters, has introduced new methods for shaping the frequency response of electrical systems. This paper seeks to extend this thread of research by investigating the frequency response shaping of electromechanical resonators which are embedded in feedforward, signal interference control architectures. In particular, mathematical models are developed to explore the behavior of linear resonators that are embedded in two prototypical signal interference control topologies, which can exhibit a variety of qualitatively distinct frequency domain behaviors with component-level tuning. Experimental approaches are then used to demonstrate the proposed designs' utility.

Previous work has shown that Fermi acceleration can be an effective heating mechanism during magnetic island coalescence, where electrons may undergo repeated reflections as the magnetic field lines contract. This energization has the potential to account for the power-law distributions of particle energy inferred from observations of solar flares. Here, we develop a generalized framework for the analysis of Fermi acceleration that can incorporate the effects of compressibility and non-uniformity along field lines, which have commonly been neglected in previous treatments of the problem. Applying this framework to the simplified case of the uniform flux tube allows us to find both the power-law scaling of the distribution function and the rate at which the power-law behavior develops. We find that a guide magnetic field of order unity effectively suppresses the development of power-law distributions.

This work describes an analytical framework suitable for the analysis of large-scale arrays of coupled resonators, including those which feature amplitude and phase dynamics, inherent element-level parameter variation, nonlinearity, and/or noise. In particular, this analysis allows for the consideration of coupled systems in which the number of individual resonators is large, extending as far as the continuum limit corresponding to an infinite number of resonators. Moreover, this framework permits analytical predictions for the amplitude and phase dynamics of such systems. The utility of this analytical methodology is explored through the analysis of a system of N non-identical resonators with global coupling, including both reactive and dissipative components, physically motivated by an electromagnetically-transduced microresonator array. In addition to the amplitude and phase dynamics, the behavior of the system as the number of resonators varies is investigated and the convergence of the discrete system to the infinite-N limit is characterized.