In magnetized, stratified environments such as the Sun's corona and solar wind, Alfvénic fluctuations ‘reflect’ from background gradients, enabling nonlinear interactions that allow their energy to dissipate into heat. This process, termed ‘reflection-driven turbulence’, likely plays a key role in coronal heating and solar-wind acceleration, explaining a range of detailed observational correlations and constraints. Building on previous works focused on the inner heliosphere, here we study the basic physics of reflection-driven turbulence using reduced magnetohydrodynamics in an expanding box – the simplest model that can capture local turbulent plasma dynamics in the super-Alfvénic solar wind. Although idealized, our high-resolution simulations and simple theory reveal a rich phenomenology that is consistent with a diverse range of observations. Outwards-propagating fluctuations, which initially have high imbalance (high cross-helicity), decay nonlinearly to heat the plasma, becoming more balanced and magnetically dominated. Despite the high imbalance, the turbulence is strong because Elsässer collisions are suppressed by reflection, leading to ‘anomalous coherence’ between the two Elsässer fields. This coherence, together with linear effects, causes the growth of ‘anastrophy’ (squared magnetic potential) as the turbulence decays, forcing the energy to rush to larger scales and forming a ‘$1/f$-range’ energy spectrum in the process. Eventually, expansion overcomes the nonlinear and Alfvénic physics, forming isolated, magnetically dominated ‘Alfvén vortices’ with minimal nonlinear dissipation. These results can plausibly explain the observed radial and wind-speed dependence of turbulence imbalance (cross-helicity), residual energy, fluctuation amplitudes, plasma heating and fluctuation spectra, as well as making a variety of testable predictions for future observations.

This investigation explores the potential formation of a relaxed equilibrium state, specifically the quadruple Beltrami state, in a three-component dusty plasma consisting of electrons, ions and negatively charged dust particles. This equilibrium state is derived by employing momentum-balanced equations along with Ampere's law. The quadruple Beltrami state is a composite of four Beltrami states, each associated with four distinct eigenvalues. Using the variational principle, we obtained the same relaxed state based on the system's constraints, which include magnetofluid energy, and the helicity of electrons, ions and dust particles. The unified flow is also derived. Dynamo action is investigated in two configurations: a rectangular geometry and a rectangular geometry with an internal conductor. Small-scale turbulent dynamo behaviour is observed in the former, while large-scale turbulent dynamo effects are noted in the latter. The magnitude of the magnetic field is found to be greater in the configuration with an internal conductor. Additionally, flow profiles are plotted as functions of Beltrami parameters and density variations of plasma species. This study contributes to the understanding of relaxation theory and the underlying physics of systems with an internal conductor, such as Saturn (planetary rings around a magnetosphere) and Jupiter magnetosphere, Uranus, Neptune, etc.

Using the ONEDFEL code we perform free electron laser simulations in the astrophysically important guide-field dominated regime. For wigglers’ (Alfvén waves) wavelengths of tens of kilometres and beam Lorentz factor ${\sim }10^3$, the resulting coherently emitted waves are in the centimetre range. Our simulations show a growth of the wave intensity over fourteen orders of magnitude, over the astrophysically relevant scale of approximately a few kilometres. The signal grows from noise (unseeded). The resulting spectrum shows fine spectral substructures, reminiscent of those observed in fast radio bursts.

The regulation of electron heat transport in high-$\beta$, weakly collisional, magnetized plasma is investigated. A temperature gradient oriented along a mean magnetic field can induce a kinetic heat-flux-driven whistler instability (HWI), which back-reacts on the transport by scattering electrons and impeding their flow. Previous analytical and numerical studies have shown that the heat flux for the saturated HWI scales as $\beta _e^{-1}$. These numerical studies, however, had limited scale separation and consequently large fluctuation amplitudes, which calls into question their relevance at astrophysical scales. To this end, we perform a series of particle-in-cell simulations of the HWI across a range of $\beta _e$ and temperature-gradient length scales under two different physical set-ups. The saturated heat flux in all of our simulations follows the expected $\beta _e^{-1}$ scaling, supporting the robustness of the result. We also use our simulation results to develop and implement several methods to construct an effective collision operator for whistler turbulence. The results point to an issue with the standard quasi-linear explanation of HWI saturation, which is analogous to the well-known $90^{\circ }$ scattering problem in the cosmic-ray community. Despite this limitation, the methods developed here can serve as a blueprint for future work seeking to characterize the effective collisionality caused by kinetic instabilities.

The decay of a turbulent magnetic field is slower with helicity than without. Furthermore, the magnetic correlation length grows faster for a helical than a non-helical field. Both helical and non-helical decay laws involve conserved quantities: the mean magnetic helicity density and the Hosking integral. Using direct numerical simulations in a triply periodic domain, we show quantitatively that in the fractionally helical case the mean magnetic energy density and correlation length are approximately given by the maximum of the values for the purely helical and purely non-helical cases. The time of switchover from one to the other decay law can be obtained on dimensional grounds and is approximately given by $I_{H}^{1/2}I_{M}^{-3/2}$, where $I_{H}$ is the Hosking integral and $I_{M}$ is the mean magnetic helicity density. An earlier approach based on the decay time is found to agree with our new result and suggests that the Hosking integral exceeds naive estimates by the square of the same resistivity-dependent factor by which also the turbulent decay time exceeds the Alfvén time. In the presence of an applied magnetic field, the mean magnetic helicity density is known to be not conserved, and we show that then also the Hosking integral is not conserved.

We study the tearing instability of a current sheet in a relativistic pair plasma with a power-law distribution function. We first estimate the growth rate analytically and then confirm the analytical results by solving numerically the dispersion equation, taking into account all exact particle trajectories within the reconnecting layer. We found that the instability is suppressed when the particle spectrum becomes harder.

Young stellar objects (YSOs) are protostars that exhibit bipolar outflows fed by accretion disks. Theories of the transition between disk and outflow often involve a complex magnetic field structure thought to be created by the disk coiling field lines at the jet base; however, due to limited resolution, these theories cannot be confirmed with observation and thus may benefit from laboratory astrophysics studies. We create a dynamically similar laboratory system by driving a $\sim$1 MA current pulse with a 200 ns rise through a $\approx$2 mm-tall Al cylindrical wire array mounted to a three-dimensional (3-D)-printed, stainless steel scaffolding. This system creates a plasma that converges on the centre axis and ejects cm-scale bipolar outflows. Depending on the chosen 3-D-printed load path, the system may be designed to push the ablated plasma flow radially inwards or off-axis to make rotation. In this paper, we present results from the simplest iteration of the load which generates radially converging streams that launch non-rotating jets. The temperature, velocity and density of the radial inflows and axial outflows are characterized using interferometry, gated optical and ultraviolet imaging, and Thomson scattering diagnostics. We show that experimental measurements of the Reynolds number and sonic Mach number in three different stages of the experiment scale favourably to the observed properties of YSO jets with $Re\sim 10^5\unicode{x2013}10^9$ and $M\sim 1\unicode{x2013}10$, while our magnetic Reynolds number of $Re_M\sim 1\unicode{x2013}15$ indicates that the magnetic field diffuses out of our plasma over multiple hydrodynamical time scales. We compare our results with 3-D numerical simulations in the PERSEUS extended magnetohydrodynamics code.

A collisionless shock structure results from the nonlinear interaction between charged particles and electromagnetic fields. Yet, a collisionless shock is globally governed by the mass, momentum and energy conservation requirements. A stable shock structure must ensure that the fluxes of the conserved quantities are constant on average, and, therefore, is determined by this necessity. Here, we study an observed high upstream temperature high Mach number shock and show that the conservation laws cannot be fulfilled unless the shock is spatially inhomogeneous along the shock front and time-dependent.

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, $\beta \gg 1$, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-$\beta$, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (${\sim }5\,\%$) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

The response tensor is derived for a relativistically streaming, strongly magnetized, one-dimensional Jüttner distribution of electrons and positrons, referred to as a pulsar plasma. This is used to produce a general treatment of wave dispersion in a pulsar plasma. Specifically, relativistic streaming, the spread in Lorentz factors in a pulsar rest frame and cyclotron resonances are taken into account. Approximations to the response tensor are derived by making approximations to relativistic plasma dispersion functions appearing in the general form of the response tensor. The cold-plasma limit, the highly relativistic limit and limits related to cyclotron resonances are considered. The theory developed in this paper has applications to generalized Faraday rotation in pulsars and magnetars.

Electron heating and acceleration in collisionless shocks is a long-standing problem. Rapid isotropization of heated electrons cannot be explained solely by the cross-shock potential. In addition, the macroscopic cross-shock potential prevents efficient reflection and injection into the diffusive acceleration regime. Recent observations have shown that small-scale electric fields are present in the shock front, together with the large-scale cross-shock potential. These small-scale fields have been found also in the upstream and downstream regions. Electron heating in shocks is produced by the combined action of the large- and small-scale fields. The large-scale potential determines the energy transferred to the electrons. The small-scale electrostatic fields scatter electrons. Here we study the scattering of electrons on the typical waveforms, namely solitary bipolar spikes and wavepackets. The main effect is the generation of backstreaming electrons with large pitch angles. It is found that wavepackets are more efficient in electron reflection in the interaction of electrons both with a single spike and with multiple spikes.

A specific set of dimensionless plasma and turbulence parameters is introduced to characterize the nature of turbulence and its dissipation in weakly collisional space and astrophysical plasmas. Key considerations are discussed for the development of predictive models of the turbulent plasma heating that characterize the partitioning of dissipated turbulent energy between the ion and electron species and between the perpendicular and parallel degrees of freedom for each species. Identifying the kinetic physical mechanisms that govern the damping of the turbulent fluctuations is a critical first step in constructing such turbulent heating models. A set of ten general plasma and turbulence parameters are defined, and reasonable approximations along with the exploitation of existing scaling theories for magnetohydrodynamic turbulence are used to reduce this general set of ten parameters to just three parameters in the isotropic temperature case: the ion plasma beta, the ion-to-electron temperature ratio and the isotropic driving wavenumber. A critical step forward in this study is to identify the dependence of all of the proposed kinetic mechanisms for turbulent damping in terms of the same set of fundamental plasma and turbulence parameters. Analytical estimations of the scaling of each damping mechanism on these fundamental parameters are presented. The power of this approach is illustrated in the development of the first phase diagram for the turbulent damping mechanisms as a function of the ion plasma beta and isotropic driving wavenumber for unity ion-to-electron temperature ratio, showing the regions of this two-dimensional parameter space in which ion Landau and transit-time damping, electron Landau and transit-time damping, ion cyclotron damping, ion stochastic heating, collisionless magnetic reconnection and kinetic ‘viscous’ heating play a role in the damping of the turbulent fluctuations.

To expand on recent work, we introduce collisional terms in the analysis of the warm ion–electron, two-fluid equations for a homogeneous plasma at rest. Consequently, the plasma is now described by six variables: the magnetisation, the ratio of masses over charges, the electron and ion sound speeds, the angle between the wave vector and the magnetic field and a new parameter describing the electron–ion collision frequency. This additional parameter does not introduce new wave modes compared with the collisionless case, but does result in complex mode frequencies. Both for the backward and forward propagating modes the imaginary components are negative and thus quantify collisional damping. We provide convenient (polynomial) expressions to quantify frequencies and damping rates in all short- and long-wavelength limits, including the cutoff and resonance limits, whilst the one-fluid magnetohydrodynamic limit is retained with the familiar undamped slow, Alfvén and fast waves. As collisions only introduce a damping, the previously introduced labelling of the wave modes S, A, F, M, O and X can be kept and assigned based on their long- and short-wavelength behaviour. The obtained damping at cutoff and resonance limits is parametrised with the collision frequency, and can be tailored to match known kinetic damping expressions. It is demonstrated that varying the angle can introduce crossings between the wave modes, as was already present in the ideal ion–electron case, but also a collision frequency exceeding a critical collision frequency can lead to crossings at angles where previously only avoided crossings were found.

Transit-time damping (TTD) is a process in which the magnetic mirror force – induced by the parallel gradient of magnetic field strength – interacts with resonant plasma particles in a time-varying magnetic field, leading to the collisionless damping of electromagnetic waves and the resulting energization of those particles through the perpendicular component of the electric field, $E_\perp$. In this study, we utilize the recently developed field–particle correlation technique to analyse gyrokinetic simulation data. This method enables the identification of the velocity-space structure of the TTD energy transfer rate between waves and particles during the damping of plasma turbulence. Our analysis reveals a unique bipolar pattern of energy transfer in the velocity-space characteristic of TTD. By identifying this pattern, we provide clear evidence of TTD's significant role in the damping of strong plasma turbulence. Additionally, we compare the TTD signature with that of Landau damping (LD). Although they both produce a bipolar pattern of phase-space energy density loss and gain about the parallel resonant velocity of the Alfvénic waves, they are mediated by different forces and exhibit different behaviours as the perpendicular velocity $v_\perp \to 0$. We also explore how the dominant damping mechanism varies with ion plasma beta $\beta _i$, showing that TTD dominates over LD for $\beta _i > 1$. This work deepens our understanding of the role of TTD in the damping of weakly collisional plasma turbulence and paves the way to seek the signature of TTD using in situ spacecraft observations of turbulence in space plasmas.

Understanding the partitioning of turbulent energy between ions and electrons in weakly collisional plasmas is crucial for the accurate interpretation of observations and modelling of various astrophysical phenomena. Many such plasmas are ‘imbalanced’, wherein the large-scale energy input is dominated by Alfvénic fluctuations propagating in a single direction. In this paper, we demonstrate that when strongly-magnetised plasma turbulence is imbalanced, nonlinear conservation laws imply the existence of a critical value of the electron plasma beta (the ratio of the thermal to magnetic pressures) that separates two dramatically different types of turbulence in parameter space. For betas below the critical value, the free energy injected on the largest scales is able to undergo a familiar Kolmogorov-type cascade to small scales where it is dissipated, heating electrons. For betas above the critical value, the system forms a ‘helicity barrier’ that prevents the cascade from proceeding past the ion Larmor radius, causing the majority of the injected free energy to be deposited into ion heating. Physically, the helicity barrier results from the inability of the system to adjust to the disparity between the perpendicular-wavenumber scalings of the free energy and generalised helicity below the ion Larmor radius; restoring finite electron inertia can annul, or even reverse, this disparity, giving rise to the aforementioned critical beta. We relate this physics to the ‘dynamic phase alignment’ mechanism (that operates under yet lower beta conditions and in pair plasmas), and characterise various other important features of the helicity barrier, including the nature of the nonlinear wavenumber-space fluxes, dissipation rates, and energy spectra. The existence of such a critical beta has important implications for heating, as it suggests that the dominant recipient of the turbulent energy, ions or electrons, can depend sensitively on the characteristics of the plasma at large scales.

In this paper, we investigate the kinetic stability of classical, collisional plasma – that is, plasma in which the mean-free-path $\lambda$ of constituent particles is short compared with the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared with the evolution time. Fluid equations are typically used to describe such plasmas, since their distribution functions are close to being Maxwellian. The small deviations from the Maxwellian distribution are calculated via the Chapman–Enskog (CE) expansion in $\lambda /L \ll 1$, and determine macroscopic momentum and heat fluxes in the plasma. Such a calculation is only valid if the underlying CE distribution function is stable at collisionless length scales and/or time scales. We find that at sufficiently high plasma $\beta$, the CE distribution function can be subject to numerous microinstabilities across a wide range of scales. For a particular form of the CE distribution function arising in strongly magnetised plasma (viz. plasma in which the Larmor periods of particles are much smaller than collision times), we provide a detailed analytic characterisation of all significant microinstabilities, including peak growth rates and their associated wavenumbers. Of specific note is the discovery of several new microinstabilities, including one at sub-electron-Larmor scales (the ‘whisper instability’) whose growth rate in certain parameter regimes is large compared with other instabilities. Our approach enables us to construct the kinetic stability maps of classical, two-species collisional plasma in terms of $\lambda$, the electron inertial scale $d_e$ and the plasma $\beta$. This work is of general consequence in emphasising the fact that high-$\beta$ collisional plasmas can be kinetically unstable; for strongly magnetised CE plasmas, the condition for instability is $\beta \gtrsim L/\lambda$. In this situation, the determination of transport coefficients via the standard CE approach is not valid.

Linearized gravity around a rotating black hole or compact object introduces the concept of a gravitomagnetic field, which originates from the matter–current in the rotating object. Plasma in proximity to this object is subsequently subjected to motion guided by this gravitomagnetic term (where mass serves as the effective charge) in addition to the conventional magnetic field. Such an interplay of fields complicates the accessible plasma waves of the system and thus merits exploration to delineate this interplay, to identify any observable signatures of the gravitomagnetic field in weakly magnetized systems, and to motivate future numerical work in a fully relativistic setting where the effects may be stronger. In this work we analyse the dispersion of the upper and lower hybrid electrostatic waves in a plasma immersed in both magnetic and gravitomagnetic fields. In particular, we discuss the effective augmentation or cancellation of the two fields under the right conditions for the upper hybrid wave. In contrast, the lower hybrid wave experiences a frequency up-shift from the gravitomagnetic field regardless of whether it is parallel or antiparallel to the magnetic field for the studied field strengths.

Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and $P_{1/3}$ multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number $M_A \approx 1.5$), supersonic (Sonic Mach number $M_S \approx 4-5$) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio $L/\delta \approx 100$, Lundquist number $S_L \approx 400$). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate = $\tau _{\text {cool}}^{-1}/\tau _{H}^{-1}\approx 100$); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Finally, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

Collisionless shocks are frequently analysed using the magnetohydrodynamics (MHD) formalism, even though MHD assumes a small mean free path. Yet, isotropy of pressure, the fruit of binary collisions and assumed in MHD, may not apply in collisionless shocks. This is especially true within a magnetized plasma, where the field can stabilize an anisotropy. In a previous article (Bret & Narayan, J. Plasma Phys., vol. 88, no. 6, 2022b, p. 905880615), a model was presented capable of dealing with the anisotropies that may arise at the front crossing. It was solved for any orientation of the field with respect to the shock front. Yet, for some values of the upstream parameters, several downstream solutions were found. Here, we complete the work started in Bret & Narayan (J. Plasma Phys., vol. 88, no. 6, 2022b, p. 905880615) by showing how to pick the physical solution out of the ones offered by the algebra. This is achieved by 2 means: (i) selecting the solution that has the downstream field obliquity closest to the upstream one. This criterion is exemplified on the parallel case and backed up by particle-in-cell simulations. (ii) Filtering out solutions which do not satisfy a criteria already invoked to trim multiple solutions in MHD: the evolutionarity criterion, that we assume valid in the collisionless case. The end result is a model in which a given upstream configuration results in a unique, or no downstream configuration (as in MHD). The largest departure from MHD is found for the case of a parallel shock.