A model of global magnetic reconnection rate in relativistic collisionless plasmas

A model of global magnetic reconnection rate in relativistic collisionless plasmas is developed and validated by the fully kinetic simulation. Through considering the force balance at the upstream and downstream of the diffusion region, we show that the global rate is bounded by a value $sim 0.3$ even when the local rate goes up to $sim O(1)$ and the local inflow speed approaches the speed of light in strongly magnetized plasmas. The derived model is general and can be applied to magnetic reconnection under widely different circumstances.

Pressure-anisotropy-driven microturbulence and magnetic-field evolution in shearing, collisionless plasma

The nonlinear state of a high-beta collisionless plasma is investigated when an imposed linear shear amplifies or diminishes a uniform magnetic field, driving pressure anisotropies and hence firehose/mirror instabilities. The evolution of the resulting microscale turbulence is considered when the shear is switched off or reversed after one shear time (mimicking local behaviour of a macroscopic flow), so a new macroscale configuration is superimposed on the microscale state left behind by the previous one. There is a threshold value of plasma beta: when $betallOmega/S$ (ion cyclotron frequency/shear rate), the emergence of firehose/mirror fluctuations driven unstable by shear and their disappearance when the shear is removed/reversed are quasi-instantaneous compared to the shear time, viz., the decay time of these fluctuations is $simbeta/Omega ll 1/S$ (this result follows from the free decay of the fluctuations being constrained by the same marginal-stability thresholds as their growth). In contrast, when $betagtrsimOmega/S$ (ultra-high beta), the old microscale state can only be removed on the shear timescale. In this regime, driven firehose fluctuations grow secularly to order-unity amplitudes, compensating for the decay of the mean field and so pinning pressure anisotropy at marginal stability with no appreciable scattering of particles---which is unlike what happens at moderate $beta$. When the shear reverses, the shearing away of this firehose turbulence compensates for the increase in the mean field and thus prevents growth of the pressure anisotropy, stopping the system from going mirror-unstable. Therefore, at ultra-high beta, the system stays close to the firehose threshold, the mirror instability is largely suppressed, while the mean magnetic energy barely changes at all. Implications for plasma dynamo and thus the origin of cosmic magnetism are discussed.

Inertial-Range Kinetic Turbulence in Pressure-Anisotropic Astrophysical Plasmas

A theoretical framework for low-frequency electromagnetic (drift-)kinetic turbulence in a collisionless, multi-species plasma is presented. The result generalises reduced magnetohydrodynamics (RMHD) and kinetic RMHD (Schekochihin et al. 2009) for pressure-anisotropic plasmas, allowing for species drifts---a situation routinely encountered in the solar wind and presumably ubiquitous in hot dilute astrophysical plasmas (e.g. intracluster medium). Two main objectives are achieved. First, in a non-Maxwellian plasma, the relationships between fluctuating fields (e.g., the Alfven ratio) are order-unity modified compared to the more commonly considered Maxwellian case, and so a quantitative theory is developed to support quantitative measurements now possible in the solar wind. The main physical feature of low-frequency plasma turbulence survives the generalisation to non-Maxwellian distributions: Alfvenic and compressive fluctuations are energetically decoupled, with the latter passively advected by the former; the Alfvenic cascade is fluid, satisfying RMHD equations (with the Alfven speed modified by pressure anisotropy and species drifts), whereas the compressive cascade is kinetic and subject to collisionless damping. Secondly, the organising principle of this turbulence is elucidated in the form of a generalised kinetic free-energy invariant. It is shown that non-Maxwellian features in the distribution function reduce the rate of phase mixing and the efficacy of magnetic stresses; these changes influence the partitioning of free energy amongst the various cascade channels. As the firehose or mirror instability thresholds are approached, the dynamics of the plasma are modified so as to reduce the energetic cost of bending magnetic-field lines or of compressing/rarefying them. Finally, it is shown that this theory can be derived as a long-wavelength limit of non-Maxwellian slab gyrokinetics.

Modeling imbalanced collisionless Alfven wave turbulence with nonlinear diffusion equations

A pair of nonlinear diffusion equations in Fourier space} is used to study the dynamics of strong Alfven-wave turbulence, from MHD to electron scales. Special attention is paid to the regime of imbalance between the energies of counter-propagating waves commonly observed in the solar wind (SW), especially in regions relatively close to the Sun. In the collisionless regime where dispersive effects arise at scales comparable to or larger than those where dissipation becomes effective, the imbalance produced by a given injection rate of generalized cross-helicity (GCH), which is an invariant, is much larger than in the corresponding collisional regime described by the usual (or reduced) magnetohydrodynamics. The combined effect of high imbalance and ion Landau damping induces a steep energy spectrum for the transverse magnetic field at sub-ion scales. This spectrum is consistent with observations in highly Alfvenic regions of the SW, such as trailing edges, but does not take the form of a transition range continued at smaller scales by a shallower spectrum. This suggests that the observed spectra displaying such a transition result from the superposition of contributions originating from various streams with different degrees of imbalance. Furthermore, when imbalanced energy injection is supplemented at small scales in an already fully developed turbulence, for example under the effect of magnetic reconnection, a significant enhancement of the imbalance at all scales is observed.

Nonrelativistic collisionless shocks in weakly magnetized electron--ion plasmas: two-dimensional particle-in-cell simulation of perpendicular shock

A two-dimensional particle-in-cell simulation is performed to investigate weakly magnetized perpendicular shocks with a magnetization parameter of 6 x 10^-5, which is equivalent to a high Alfven Mach number M_A of ~130. It is shown that current filaments form in the foot region of the shock due to the ion-beam--Weibel instability (or the ion filamentation instability) and that they generate a strong magnetic field there. In the downstream region, these current filaments also generate a tangled magnetic field that is typically 15 times stronger than the upstream magnetic field. The thermal energies of electrons and ions in the downstream region are not in equipartition and their temperature ratio is T_e / T_i ~ 0.3 - 0.4. Efficient electron acceleration was not observed in our simulation, although a fraction of the ions are accelerated slightly on reflection at the shock. The simulation results agree very well with the Rankine-Hugoniot relations. It is also shown that electrons and ions are heated in the foot region by the Buneman instability (for electrons) and the ion-acoustic instability (for both electrons and ions). However, the growth rate of the Buneman instability is significantly reduced due to the relatively high temperature of the reflected ions. For the same reason, ion-ion streaming instability does not grow in the foot region.