A magnetohydrodynamic (MHD) shock is a discontinuity through which plasma flows, while the mass, momentum and energy are conserved, and the entropy increases. In a collisionless system, the entropy is strictly constant along the phase-space trajectory. The MHD required entropy increase is attributed to dissipative processes causing relaxation to local thermodynamic equilibrium. The search for such dissipative processes in the shock observations usually explores variations of the entropy density or entropy per particle across the shock front. We show that both quantities may increase or decrease across the shock front even when the entropy is strictly conserved. Entropy production should be assessed using the entropy flux density.

The solar wind is observed to undergo substantial heating as it expands through the heliosphere, with measured temperature profiles exceeding those expected from adiabatic cooling. A plausible source of this heating is reflection-driven turbulence (RDT), in which gradients in the background Alfvén speed partially reflect outward-propagating Alfvén waves, seeding counter-propagating fluctuations that interact and dissipate via turbulence. Previous RDT models assume a radial-background magnetic field, but at larger radii the interplanetary field is known to be twisted into the Parker spiral (PS). Here, we generalise RDT phenomenology to include a PS, using three-dimensional expanding-box magnetohydrodynamic simulations to test the ideas and compare the resulting turbulence with the radial-background-field case. We argue that the underlying RDT dynamics remains broadly similar with a PS, but the controlling scales change: as the azimuthal field grows it ‘cuts across’ perpendicularly stretched, pancake-like eddies, producing outer scales perpendicular to the magnetic field that are much smaller than in the radial-background case. Consequently, the outer-scale nonlinear turnover time increases more slowly with heliocentric distance in PS geometry, weakening the tendency (seen in radial-background models) for the cascade to ‘freeze’ into quasi-static, magnetically dominated structures. This allows the system to dissipate a larger fraction of the fluctuation energy as heat, also implying that the turbulence remains strongly imbalanced (with high normalised cross-helicity) out to larger heliocentric distances. We complement our heating results with a detailed characterisation of the turbulence (e.g. spectra, switchbacks and compressive fractions), providing a set of concrete predictions for comparison with spacecraft observations.

This study investigates ion kinetic effects during the parametric decay instability (PDI) of parallel-propagating Alfvén waves under plasma conditions characteristic of the Earth’s ionosphere. By using a series of hybrid particle-in-cell simulations, we examine the evolution of ion velocity distribution functions (VDFs) in ultra-low-beta plasmas. Our numerical campaign systematically explores the dependence on key parameters (plasma beta, pump-wave amplitude and polarisation, and ion composition). To emphasise the role of kinetic effects, we choose to trigger the PDI with a dispersive mother wave with wavelength comparable to the ion characteristic inertial length. Our results reveal pronounced non-thermal VDF modifications, including parallel heating and the formation of secondary ion beams, linked to the nonlinear evolution of parametric decay instability. By varying the plasma beta and the pump-wave amplitude, we identify a critical regime where rapid and complete broadening of the velocity distribution function is observed, triggering bidirectional ion acceleration. Notably, simulations modelling realistic ionospheric conditions demonstrate that even low-amplitude Alfvénic perturbations can induce significant VDF spreading and ion beam generation, with hydrogen ions exhibiting stronger effects than oxygen. These non-thermal microscopic processes offer a plausible mechanism for particle precipitation in space weather events. This work represents the first comprehensive study with hybrid simulations of PDI-driven ion kinetics in ultra-low-beta plasmas, providing quantitative estimates for the time delay between electromagnetic wave impact and ion VDF modification, and new insights into wave–particle interactions that may contribute to ion acceleration, precipitation processes and space plasma dynamics.

Inertial Alfvén waves are thought to accelerate electrons to auroral energies via their parallel electric field in the Earth’s magnetosphere. During active geomagnetic times, it is estimated that a significant percentage of electron precipitation energy into the Earth’s ionosphere can be attributed to these waves. However, self-consistent wave/particle interactions of inertial Alfvén waves with the accelerated electron population are not well understood. We show that recent self-consistent models have a strong nonlinearity in them. A reduced set of equations which describe this nonlinear steepening is derived and shown to agree with drift-kinetic simulations and other published studies. From this reduced set of equations, many properties of the nonlinearity are derived and shown to agree with simulations. This includes the time and length scales and connecting the speed of the wave to the perturbation maximum value.

We report on a first-principles numerical study of magnetic reconnection in plasmas with different initial ion-to-electron temperature ratios. In cases where this ratio is significantly below unity, we observe intense wave activity in the diffusion region, driven by the ion-acoustic instability. Our analysis shows that the dominant macroscopic effect of this instability is to drive substantial ion heating. In contrast to earlier studies reporting significant anomalous resistivity, we find that anomalous contributions due to the ion-acoustic instability are minimal. These results shed light on the dynamical impact of this instability on reconnection processes, offering new insights into the fundamental physics governing collisionless reconnection.

We study the interaction of an ion with a fluctuation in the electromagnetic fields that is localised in both space and time. We study the scale dependence of the interaction in both space and time, deriving a generic form for the ion’s energy change, which involves an exponential cutoff based on the characteristic time scale of the electromagnetic fluctuation. This leads to diffusion in energy in both $v_\perp$ and $v_\parallel$. We show how to apply our results to general plasma physics phenomena, and specifically to Alfvénic turbulence and to reconnection. Our theory can be viewed as a unification of previous models of stochastic ion heating, cyclotron heating and reconnection heating in a single theoretical framework.

Two-fluid simulations using local Landau-fluid closures derived from linear theory provide an efficient computational framework for plasma modelling, since they bridge the gap between computationally intensive kinetic simulations and fluid descriptions. Their accuracy in representing kinetic effects depends critically on the validity of the linear approximation used in the derivation: the plasma should not be too far from local thermodynamic equilibrium (LTE). However, many of the problems where these models are of particular interest (such as plasma turbulence and instabilities) are in fact quite far from LTE. The question then arises as to whether kinetic-scale processes are still sufficiently well captured outside of the theoretical regime of applicability of the closure. In this paper, we show that two-fluid simulations with Landau-fluid closures can effectively reproduce the energy spectra obtained with fully kinetic Vlasov simulations, used as references, as long as the local closure parameter is appropriately chosen. Our findings validate the usage of two-fluid simulations with a Landau-fluid closure as a possible alternative to fully kinetic simulations of turbulence, in cases where being able to simulate extremely large domains is of particular interest.

We present a simple, analytically solvable magnetohydrodynamics model of current sheet formation through X-point collapse under optically thin radiative cooling. Our results show that cooling accelerates the collapse of the X-point along the inflows, but strong cooling can arrest or even reverse the current sheet elongation in the outflow direction. Hence, we detail a modification to the radiatively cooled Sweet–Parker model developed by Uzdensky & McKinney (Phys. Plasmas, 1962, vol. 18, issue 4, p. 042105) to allow for varying current sheet length. The steady-state solution shows that, when radiative cooling dominates compressional heating, the current sheet length is shorter than the system size, with an increased reconnection rate compared with the classical Sweet–Parker rate. The model and subsequent results lay out the groundwork for a more complete theoretical understanding of magnetic reconnection in regimes dominated by optically thin radiative cooling.

We develop a new scaling theory for the resistive tearing mode instability of a current sheet with a strong shear flow across the layer. The growth rate decreases with increasing flow shear and is completely stabilized as the shear flow becomes Alfvénic: both in the constant-$\varPsi$ regime, as in previous results, but we also show that the growth rate is in fact suppressed more strongly in the nonconstant-$\varPsi$ regime. As a consequence, for sufficiently large flow shear, the maximum of the growth rate is always affected by the shear suppression, and the wavenumber at which this maximum growth rate is attained is an increasing function of the strength of the flow shear. These results may be important for the onset of reconnection in imbalanced MHD turbulence.

We study the evolution of collisionless plasmas that, due to their macroscopic evolution, are susceptible to the firehose instability, using both analytic theory and hybrid-kinetic particle-in-cell simulations. We establish that, depending on the relative magnitude of the plasma $\beta$, the characteristic time scale of macroscopic evolution and the ion-Larmor frequency, the saturation of the firehose instability in high-$\beta$ plasmas can result in three qualitatively distinct thermodynamic (and electromagnetic) states. By contrast with the previously identified ‘ultra-high-beta’ and ‘Alfvén-inhibiting’ states, the newly identified ‘Alfvén-enabling’ state, which is realised when the macroscopic evolution time $\tau$ exceeds the ion-Larmor frequency by a $\beta$-dependent critical parameter, can support linear Alfvén waves and Alfvénic turbulence because the magnetic tension associated with the plasma’s macroscopic magnetic field is never completely negated by anisotropic pressure forces. We characterise these states in detail, including their saturated magnetic-energy spectra. The effective collision operator associated with the firehose fluctuations is also described; we find it to be well approximated in the Alfvén-enabling state by a simple quasi-linear pitch-angle scattering operator. The box-averaged collision frequency is $\nu _{\textrm {eff}} \sim \beta /\tau$, in agreement with previous results, but certain subpopulations of particles scatter at a much larger (or smaller) rate depending on their velocity in the direction parallel to the magnetic field. Our findings are essential for understanding low-collisionality astrophysical plasmas including the solar wind, the intracluster medium of galaxy clusters and black hole accretion flows. We show that all three of these plasmas are in the Alfvén-enabling regime of firehose saturation and discuss the implications of this result.

We describe a new model for the study of weakly collisional, magnetised plasmas derived from exploiting the separation of the dynamics parallel and perpendicular to the magnetic field. This unique system of equations retains the particle dynamics parallel to the magnetic field while approximating the perpendicular dynamics through a spectral expansion in the perpendicular degrees of freedom, analogous to moment-based fluid approaches. In so doing, a hybrid approach is obtained that is computationally efficient enough to allow for larger-scale modelling of plasma systems while eliminating a source of difficulty in deriving fluid equations applicable to magnetised plasmas. We connect this system of equations to historical asymptotic models and discuss advantages and disadvantages of this approach, including the extension of this parallel-kinetic-perpendicular moment beyond the typical region of validity of these more traditional asymptotic models. This paper forms the first of a multi-part series on this new model, covering the theory and derivation, alongside demonstration benchmarks of this approach that include shocks and magnetic reconnection.

In one of the leading theories for the origin of the solar wind, photospheric motions launch Alfvén waves (AWs) that propagate along open magnetic-field lines through the solar atmosphere and into the solar wind. The radial variation in the Alfvén speed causes some of the AWs to reflect, and counter-propagating AWs subsequently interact to produce Alfveńic turbulence, in which AW energy cascades from long wavelengths to short wavelengths and dissipates, heating the plasma. In this paper we develop a one-dimensional two-fluid solar-wind model that includes Alfvénic turbulence, proton temperature anisotropy and a novel method for apportioning the turbulent heating rate between parallel proton heating, perpendicular proton heating and electron heating. We employ a turbulence model that accounts for recent observations from NASA’s Parker Solar Probe, which find that AW fluctuations in the near-Sun solar wind are intermittent and less anisotropic than in previous models of anisotropic magnetohydrodynamic turbulence. Our solar-wind model reproduces a wide range of remote observations of the corona and in-situ measurements of the solar wind, and our turbulent heating model consists of analytic equations that could be usefully incorporated into other solar-wind models and numerical models of more distant astrophysical plasmas.

From the near-Earth solar wind to the intracluster medium of galaxy clusters, collisionless, high-beta, magnetized plasmas pervade our universe. Energy and momentum transport from large-scale fields and flows to small-scale motions of plasma particles is ubiquitous in these systems, but a full picture of the underlying physical mechanisms remains elusive. The transfer is often mediated by a turbulent cascade of Alfvénic fluctuations as well as a variety of kinetic instabilities; these processes tend to be multi-scale and/or multi-dimensional, which makes them difficult to study using spacecraft missions and numerical simulations alone. Meanwhile, existing laboratory devices struggle to produce the collisionless, high ion beta ($\beta _i \gtrsim 1$), magnetized plasmas across the range of scales necessary to address these problems. As envisioned in recent community planning documents, it is therefore important to build a next generation laboratory facility to create a $\beta _i \gtrsim 1$, collisionless, magnetized plasma in the laboratory for the first time. A working group has been formed and is actively defining the necessary technical requirements to move the facility towards a construction-ready state. Recent progress includes the development of target parameters and diagnostic requirements as well as the identification of a need for source-target device geometry. As the working group is already leading to new synergies across the community, we anticipate a broad community of users funded by a variety of federal agencies (including National Aeronautics and Space Administration, Department of Energy and National Science Foundation) to make copious use of the future facility.

The Vlasov–Maxwell equations provide kinetic simulations of collisionless plasmas, but numerically solving them on classical computers is often impractical. This is due to the computational resource constraints imposed by the time evolution in the six-dimensional phase space, which requires broad spatial and temporal scales. The novelty of this study is to implement a quantum–classical hybrid Vlasov–Maxwell solver and the rigorous numerical scheme evaluation by numerical simulations. Specifically, the Vlasov solver implements the Hamiltonian simulation based on quantum singular value transformation, coupled with a classical Maxwell solver. We perform numerical simulation of a one-dimensional advection test and a one-spatial-dimension, one-velocity-dimension two-stream instability test on the Qiskit-Aer-GPU quantum circuit emulator with an A100 GPU. The computational complexity of our quantum algorithm can potentially be reduced from the classical $\mathcal{O}(N^6T^2/\epsilon )$ to $\mathcal{O}\left (\text{poly}(\log {N})\left (NT+T\log \left (T/\epsilon \right )\right )\right )$ for the $N$ grid system, simulation time $T$ and error tolerance $\epsilon$ in the limit where the number of queries is large enough and the error is small enough. Furthermore, the numerical analysis reveals that our quantum algorithm is robust under larger time steps compared with classical algorithms with the constraint of Courant–Friedrichs–Lewy condition.

Electron-only reconnection (E-REC) is a process recently observed in the Earth’s magnetosheath, where magnetic reconnection occurs at electron kinetic scales, and ions do not couple to the reconnection process. Electron-only reconnection is likely to have a significant impact on the energy conversion and dissipation of turbulence cascades at kinetic scales in some settings. This paper investigates E-REC under different intensities of strong guide fields (the ratio between the guide field and the in-plane asymptotic field strength is 5, 10 and 20, respectively) via two-dimensional fully kinetic particle-in-cell simulations, focusing on electron heating. The simulations are initialized with a force-free current sheet equilibrium under various intensities of strong guide fields. Similarly to previous experimental studies, electron temperature anisotropy along separatrices is observed, which is found to be mainly caused by the variations of parallel temperature. Both regions of anisotropy and parallel temperature increase/decrease along separatrices become thinner with increasing guide fields. Besides, we find a transition from a quadrupolar to a hexapolar (six-polar) to an octopolar (eight-polar) structure in temperature anisotropy and parallel temperature as the guide field intensifies. Non-Maxwellian electron velocity distribution functions (EVDFs) at different locations in the three simulations are observed. Our results show that parallel electron velocity varies notably with different guide field intensities and finite parallel electron heat flux density is observed. The three simulations exhibit features of the Chew–Goldberger–Low theory, with the level of consistency increasing as the guide field strength increases. This explains the electron parallel temperature variations and the shape of the EVDFs observed along the separatrices. This work may provide insights into the understanding of electron heating and parallel heat flux density in E-REC observed in the turbulent magnetosheath.

Collisionless shocks are frequently analysed using the magnetohydrodynamic (MHD) formalism, even though the required collisionality hypothesis is not fulfilled. In a previous work (Bret & Narayan, 2018 J. Plasma Phys. vol. 84, p. 905840604), we presented a model of collisionless shock displaying an important departure from the expected MHD behaviour, in the case of a strong flow aligned magnetic field. This model was non-relativistic. Here, it is extended to the relativistic regime, considering zero upstream pressure and upstream Lorentz factor $\gg 1$. The result agrees satisfactorily with Particle-in-Cell simulations and shows a similar, and important, departure from the MHD prediction. In the strong field regime, the density jump $r$, seen in the downstream frame, behaves like $r \sim 2 + 1/\gamma _{\mathrm{up}}$, while MHD predicts 4 ($\gamma _{\mathrm{up}}$ is the Lorentz factor of the upstream measured in the downstream frame). Only pair plasmas are considered.

The presence of multi-component protons with their distinct features is confirmed by various space missions in the Earth’s outer magnetosphere regions. Isotropic cold protons and anisotropic hot protons significantly influence/modify the dispersion behaviour of various modes and instabilities and regulate the magnetospheric dynamics effectively. Our present study pays attention to the left-hand-polarised proton cyclotron mode, which gets unstable in the large proton temperature anisotropy condition, i.e. $T_{\perp p}\gt T_{\parallel p}$. Such favourable thermal conditions for protons are extensively observed during the compression of the solar wind against the Earth’s magnetic field. To reveal the wave dynamics in more detail, i.e. time-scale variations in the cold and hot proton temperatures and resulting wave-energy density, we further allow the time evolution of our model bi-Maxwellian distribution function in response to the proton cyclotron instability. Based on velocity-moment techniques, we formulated a set of equations comprising an instantaneous dispersion relation, dynamical perpendicular and parallel temperature relations and a wave-energy density equation. For the graphical illustrations of our mathematical results, we choose initial conditions that are relevant to magnetospheric space environments and reported in various experimental studies. Our exact numerical analysis shows the notable impact of hot proton temperature anisotropy and relative density on the real frequency, growth rate, evolution of initial distributions and wave-energy density of the proton cyclotron instability. Such detailed outcomes will be quite helpful for global/local magnetospheric experimental and simulation studies.

In this study, we conducted an electrical analysis of the effects of cold plasma on the properties of distilled water, using a corona discharge in a tip–plane configuration. The discharge was initiated by applying a voltage of 7.17 kV with a 2 mm gap between the tip and the water surface. We investigated the impact of plasma treatment on the total dissolved solids (TDS) and conductivity of 20 mL of distilled water, with exposure times ranging from 2 to 12 min. The results show that plasma treatment leads to a significant increase in conductivity and TDS, with a proportional increase relative to the exposure time. In addition to these measurements, we performed a detailed electrical analysis to evaluate the energy efficiency of the plasma treatment. This analysis involved calculating the useful power and energy efficiency using an equivalent electrical model of the corona discharge reactor, as direct measurement of these parameters is challenging in this context. The model allowed us to calculate energy consumption and analyse the electrical behaviour of the system throughout the treatment process. This study also enables us to monitor, control and optimize the energy during plasma treatment, providing insights into the energy dynamics involved. The findings have potential applications in improving energy efficiency in industrial and environmental processes.

Electron energisation by magnetic reconnection has historically been studied in the Lagrangian guiding-centre framework. Insights from such studies include that Fermi acceleration in magnetic islands can accelerate electrons to high energies. An alternative Eulerian fluid formulation of electron energisation was recently used to study electron energisation during magnetic reconnection in the absence of magnetic islands. Here, we use particle-in-cell simulations to compare the Eulerian and Lagrangian models of electron energisation in a set-up where reconnection leads to magnetic island formation. We find the largest energisation at the edges of magnetic islands. There, energisation related to the diamagnetic drift dominates in the Eulerian model, while the Fermi related term dominates in the Lagrangian model. The models predict significantly different energisation rates locally. A better agreement is found after integrating over the simulation domain. We show that strong magnetic curvature can break the magnetic moment conservation assumed by the Lagrangian model, leading to erroneous results. The Eulerian fluid model is a complete fluid description and accurately models bulk energisation. However, local measurements of its constituent energisation terms need not reflect locations where plasma is heated or accelerated. The Lagrangian guiding centre model can accurately describe the energisation of particles, but it cannot describe the evolution of the fluid energy. We conclude that while both models can be valid, they describe two fundamentally different quantities, and care should be taken when choosing which model to use.

We derive a set of simplified equations that can be used for numerical studies of reduced magnetohydrodynamic turbulence within a small patch of the radially expanding solar wind. We allow the box to be either stationary in the Sun’s frame or to be moving at an arbitrary velocity along the background magnetic-field lines, which we take to be approximately radial. We focus in particular on the case in which the box moves at the same speed as outward-propagating Alfvén waves. To aid in the design and optimization of future numerical simulations, we express the equations in terms of scalar potentials and Clebsch coordinates. The equations we derive will be particularly useful for conducting high-resolution numerical simulations of reflection-driven magnetohydrodynamic turbulence in the solar wind, and may also be useful for studying turbulence within other astrophysical outflows.