Vlasov methods in space physics and astrophysics

Magnetosheath Jet Formation Influenced by Parameters in Solar Wind Structures

Magnetosheath jets are dynamic pressure enhancements observed in the terrestrial magnetosheath. Their generation mechanisms are currently debated but the majority of jets can be linked to foreshock processes. Recent results showed that jets are less numerous when coronal mass ejections (CMEs) cross the magnetosheath and more numerous when stream interaction regions (SIRs) cross it. Here, we show for the first time how the pronounced substructures of CMEs and SIRs are related to jet production. We distinguish between compression and magnetic ejecta (ME) regions for the CME as well as compression region associated with the stream interface and high-speed streams (HSSs) for the SIR. Based on THEMIS and OMNI data covering 2008-2021, we show the 2D probability distribution of jet occurrence using the cone angle and Alfvén Mach number. We compare this distribution with the values within each solar wind (SW) structure. We find that both high cone angles and low Alfvén Mach numbers within CME-MEs are unfavorable for jet production as they may inhibit a well-defined foreshock region. 1D histograms of all parameters show, which SW parameters govern jet occurrence in each SW structure. In terms of the considered parameters the most favorable conditions for jet generation are found for HSSs due to their associated low cone angles, low densities, and low magnetic field strengths.

Transmission of foreshock waves through Earth's bow shock

The Earth's magnetosphere and its bow shock, which is formed by the interaction of the supersonic solar wind with the terrestrial magnetic field, constitute a rich natural laboratory enabling in situ investigations of universal plasma processes. Under suitable interplanetary magnetic field conditions, a foreshock with intense wave activity forms upstream of the bow shock. So-called 30 s waves, named after their typical period at Earth, are the dominant wave mode in the foreshock and play an important role in modulating the shape of the shock front and affect particle reflection at the shock. These waves are also observed inside the magnetosphere and down to the Earth's surface, but how they are transmitted through the bow shock remains unknown. By combining state-of-the-art global numerical simulations and spacecraft observations, we demonstrate that the interaction of foreshock waves with the shock generates earthward-propagating, fast-mode waves, which reach the magnetosphere. These findings give crucial insight into the interaction of waves with collisionless shocks in general and their impact on the downstream medium.

Energy Flux Through the Magnetopause During Flux Transfer Events in Hybrid-Vlasov 2D Simulations

Solar wind-magnetosphere coupling drives magnetospheric dynamic phenomena by enabling energy exchange between magnetospheric and solar wind plasmas. In this study, we examine two-dimensional noon-midnight meridional plane simulation runs of the global hybrid-Vlasov code Vlasiator with southward interplanetary magnetic field driving. We compute the energy flux, which consists of the Poynting flux and hydrodynamic energy flux components, through the Earth's magnetopause during flux transfer events (FTEs). The results demonstrate the spatiotemporal variations of the energy flux along the magnetopause during an FTE, associating the FTE leading (trailing) edge with an energy injection into (escape from) the magnetosphere on the dayside. Furthermore, FTEs traveling along the magnetopause transport energy to the nightside magnetosphere. We identify the tail lobes as a primary entry region for solar wind energy into the magnetosphere, consistent with results from global magnetohydrodynamic simulations and observations.

Quantum-inspired method for solving the Vlasov-Poisson equations

Kinetic simulations of collisionless (or weakly collisional) plasmas using the Vlasov equation are often infeasible due to high-resolution requirements and the exponential scaling of computational cost with respect to dimension. Recently, it has been proposed that matrix product state (MPS) methods, a quantum-inspired but classical algorithm, can be used to solve partial differential equations with exponential speed-up, provided that the solution can be compressed and efficiently represented as a MPS within some tolerable error threshold. In this work, we explore the practicality of MPS methods for solving the Vlasov-Poisson equations for systems with one coordinate in space and one coordinate in velocity, and find that important features of linear and nonlinear dynamics, such as damping or growth rates and saturation amplitudes, can be captured while compressing the solution significantly. Furthermore, by comparing the performance of different mappings of the distribution functions onto the MPS, we develop an intuition of the MPS representation and its behavior in the context of solving the Vlasov-Poisson equations, which will be useful for extending these methods to higher-dimensional problems.

Electron Signatures of Reconnection in a Global eVlasiator Simulation

Geospace plasma simulations have progressed toward more realistic descriptions of the solar wind-magnetosphere interaction from magnetohydrodynamic to hybrid ion-kinetic, such as the state-of-the-art Vlasiator model. Despite computational advances, electron scales have been out of reach in a global setting. eVlasiator, a novel Vlasiator submodule, shows for the first time how electromagnetic fields driven by global hybrid-ion kinetics influence electrons, resulting in kinetic signatures. We analyze simulated electron distributions associated with reconnection sites and compare them with Magnetospheric Multiscale (MMS) spacecraft observations. Comparison with MMS shows that key electron features, such as reconnection inflows, heated outflows, flat-top distributions, and bidirectional streaming, are in remarkable agreement. Thus, we show that many reconnection-related features can be reproduced despite strongly truncated electron physics and an ion-scale spatial resolution. Ion-scale dynamics and ion-driven magnetic fields are shown to be significantly responsible for the environment that produces electron dynamics observed by spacecraft in near-Earth plasmas.

Quasi-Parallel Shock Reformation Seen by Magnetospheric Multiscale and Ion-Kinetic Simulations

Shock waves in collisionless plasmas are among the most efficient particle accelerators in space. Shock reformation is a process important to plasma heating and acceleration, but direct observations of reformation at quasi-parallel shocks have been lacking. Here, we investigate Earth's quasi-parallel bow shock with observations by the four Magnetospheric Multiscale spacecraft. The multi-spacecraft observations provide evidence of short large-amplitude magnetic structures (SLAMS) causing reformation of the quasi-parallel shock. We perform an ion-kinetic Vlasiator simulation of the bow shock and show that SLAMS reforming the bow shock recreates the multi-spacecraft measurements. This provides a method for identifying shock reformation in the future.

PIC methods in astrophysics: simulations of relativistic jets and kinetic physics in astrophysical systems

The Particle-In-Cell (PIC) method has been developed by Oscar Buneman, Charles Birdsall, Roger W. Hockney, and John Dawson in the 1950s and, with the advances of computing power, has been further developed for several fields such as astrophysical, magnetospheric as well as solar plasmas and recently also for atmospheric and laser-plasma physics. Currently more than 15 semi-public PIC codes are available which we discuss in this review. Its applications have grown extensively with increasing computing power available on high performance computing facilities around the world. These systems allow the study of various topics of astrophysical plasmas, such as magnetic reconnection, pulsars and black hole magnetosphere, non-relativistic and relativistic shocks, relativistic jets, and laser-plasma physics. We review a plethora of astrophysical phenomena such as relativistic jets, instabilities, magnetic reconnection, pulsars, as well as PIC simulations of laser-plasma physics (until 2021) emphasizing the physics involved in the simulations. Finally, we give an outlook of the future simulations of jets associated to neutron stars, black holes and their merging and discuss the future of PIC simulations in the light of petascale and exascale computing.

Ballooning-Interchange Instability in the Near-Earth Plasma Sheet and Auroral Beads: Global Magnetospheric Modeling at the Limit of the MHD Approximation

Explosive magnetotail activity has long been understood in the context of its auroral manifestations. While global models have been used to interpret and understand many magnetospheric processes, the temporal and spatial scales of some auroral forms have been inaccessible to global modeling creating a gulf between observational and theoretical studies of these phenomena. We present here an important step toward bridging this gulf using a newly developed global magnetosphere-ionosphere model with resolution capturing ≲ 30 km azimuthal scales in the auroral zone. In a global magnetohydrodynamic (MHD) simulation of the growth phase of a synthetic substorm, we find the self-consistent formation and destabilization of localized magnetic field minima in the near-Earth magnetotail. We demonstrate that this destabilization is due to ballooning-interchange instability which drives earthward entropy bubbles with embedded magnetic fronts. Finally, we show that these bubbles create localized field-aligned current structures that manifest in the ionosphere with properties matching observed auroral beads.

Comparative Analysis of the Vlasiator Simulations and MMS Observations of Multiple X-Line Reconnection and Flux Transfer Events

The Vlasiator hybrid-Vlasov code was developed to investigate global magnetospheric dynamics at ion-kinetic scales. Here we focus on the role of magnetic reconnection in the formation and evolution of magnetic islands at the low-latitude magnetopause, under southward interplanetary magnetic field conditions. The simulation results indicate that (1) the magnetic reconnection ion kinetics, including the Earthward pointing Larmor electric field on the magnetospheric side of an X-point and anisotropic ion distributions, are well-captured by Vlasiator, thus enabling the study of reconnection-driven magnetic island evolution processes, (2) magnetic islands evolve due to continuous reconnection at adjacent X-points, "coalescence" which refers to the merging of neighboring islands to create a larger island, "erosion" during which an island loses magnetic flux due to reconnection, and "division" which involves the splitting of an island into smaller islands, and (3) continuous reconnection at adjacent X-points is the dominant source of magnetic flux and plasma to the outer layers of magnetic islands resulting in cross-sectional growth rates up to + 0.3 RE 2/min. The simulation results are compared to the Magnetospheric Multiscale (MMS) measurements of a chain of ion-scale flux transfer events (FTEs) sandwiched between two dominant X-lines. The MMS measurements similarly reveal (1) anisotropic ion populations and (2) normalized reconnection rate ~0.18, in agreement with theory and the Vlasiator predictions. Based on the simulation results and the MMS measurements, it is estimated that the observed ion-scale FTEs may grow Earth-sized within ~10 min, which is comparable to the average transport time for FTEs formed in the subsolar region to the high-latitude magnetopause. Future simulations shall revisit reconnection-driven island evolution processes with improved spatial resolutions.