The partition of irreversible heating between ions and electrons in compressively driven (but subsonic) collisionless turbulence is investigated by means of nonlinear gyrokinetic simulations. We derive a prescription for the ion-to-electron heating ratio $Q_{\text{i}}/Q_{\text{e}}$ as a function of the compressive-to-Alfv\'enic driving power ratio $P_{\text{compr}}/P_{\text{AW}}$, of the ratio of ion thermal pressure to magnetic pressure $\beta_{\text{i}}$, and of the ratio of ion-to-electron background temperatures $T_{\text{i}}/T_{\text{e}}$. It is shown that $Q_{\text{i}}/Q_{\text{e}}$ is an increasing function of $P_{\text{compr}}/P_{\text{AW}}$. When the compressive driving is sufficiently large, $Q_{\text{i}}/Q_{\text{e}}$ approaches $\simeq P_{\text{compr}}/P_{\text{AW}}$. This indicates that, in turbulence with large compressive fluctuations, the partition of heating is decided at the injection scales, rather than at kinetic scales. Analysis of phase-space spectra shows that the energy transfer from inertial-range compressive fluctuations to sub-Larmor-scale kinetic Alfv\'en waves is absent for both low and high $\beta_{\text{i}}$, meaning that the compressive driving is directly connected to the ion entropy fluctuations, which are converted into ion thermal energy. This result suggests that preferential electron heating is a very special case requiring low $\beta_{\text{i}}$ and no, or weak, compressive driving. Our heating prescription has wide-ranging applications, including to the solar wind and to hot accretion disks such as M87 and Sgr A*.

It is shown that in low-beta, weakly collisional plasmas, such as the solar corona, some instances of the solar wind, the aurora, inner regions of accretion discs, their coronae, and some laboratory plasmas, Alfv\'enic fluctuations produce no ion heating within the gyrokinetic approximation, i.e., as long as their amplitudes (at the Larmor scale) are small and their frequencies stay below the ion Larmor frequency (even as their spatial scales can be above or below the ion Larmor scale). Thus, all low-frequency ion heating in such plasmas is due to compressive fluctuations ("slow modes"). Because these fluctuations energetically decouple from the Alfv\'enic ones already in the inertial range, the above conclusion means that the energy partition between ions and electrons in low-beta plasmas is decided at the outer scale, where turbulence is launched, and can be determined from magnetohydrodynamic (MHD) models of the relevant astrophysical systems. Any additional ion heating must come from non-gyrokinetic mechanisms such as cyclotron heating or the stochastic heating owing to distortions of ions' Larmor orbits. An exception to these conclusions occurs in the Hall limit, i.e., when the ratio of the ion to electron temperatures is as low as the ion beta (equivalently, the electron beta is order unity). In this regime, slow modes couple to Alfv\'enic ones well above the Larmor scale (viz., at the ion inertial or ion sound scale), so the Alfv\'enic and compressive cascades join and then separate again into two cascades of fluctuations that linearly resemble kinetic Alfv\'en and ion cyclotron waves, with the former heating electrons and the latter ions. The two cascades are shown to decouple, scalings for them are derived, and it is argued physically that the two species will be heated by them at approximately equal rates.

Does overall thermal equilibrium exist between ions and electrons in a weakly collisional, magnetised, turbulent plasma---and, if not, how is thermal energy partitioned between ions and electrons? This is a fundamental question in plasma physics, the answer to which is also crucial for predicting the properties of far-distant astronomical objects such as accretion discs around black holes. In the context of discs, this question was posed nearly two decades ago and has since generated a sizeable literature. Here we provide the answer for the case in which energy is injected into the plasma via Alfv\'enic turbulence: collisionless turbulent heating typically acts to disequilibrate the ion and electron temperatures. Numerical simulations using a hybrid fluid-gyrokinetic model indicate that the ion-electron heating-rate ratio is an increasing function of the thermal-to-magnetic energy ratio, $\beta_\mathrm{i}$: it ranges from $\sim0.05$ at $\beta_\mathrm{i}=0.1$ to at least $30$ for $\beta_\mathrm{i} \gtrsim 10$. This energy partition is approximately insensitive to the ion-to-electron temperature ratio $T_\mathrm{i}/T_\mathrm{e}$. Thus, in the absence of other equilibrating mechanisms, a collisionless plasma system heated via Alfv\'enic turbulence will tend towards a nonequilibrium state in which one of the species is significantly hotter than the other, viz., hotter ions at high $\beta_\mathrm{i}$, hotter electrons at low $\beta_\mathrm{i}$. Spectra of electromagnetic fields and the ion distribution function in 5D phase space exhibit an interesting new magnetically dominated regime at high $\beta_i$ and a tendency for the ion heating to be mediated by nonlinear phase mixing ("entropy cascade") when $\beta_\mathrm{i}\lesssim1$ and by linear phase mixing (Landau damping) when $\beta_\mathrm{i}\gg1$

This study shows that a relativistic Hall effect significantly changes the properties of wave propagation by deriving a linear dispersion relation for relativistic Hall magnetohydrodynamics (HMHD). Whereas, in non- relativistic HMHD, the phase and group velocities of fast magnetosonic wave become anisotropic with an increasing Hall effect, the relativistic Hall effect brings upper bounds to the anisotropies. The Alfve\'en wave group velocity with strong Hall effect also becomes less anisotropic than non-relativistic case. Moreover, the group velocity surfaces of Alfv\'en and fast waves coalesce into a single surface in the direction other than near perpendicular to the ambient magnetic field. It is also remarkable that a characteristic scale length of the relativistic HMHD depends on ion temperature, magnetic field strength, and density while the non-relativistic HMHD scale length, i.e., ion skin depth, depends only on density. The modified characteristic scale length increases as the ion temperature increases and decreases as the magnetic field strength increases.

Two types of Eulerian action principles for relativistic extended magnetohydrodynamics (MHD) are formulated. With the first, the action is extremized under the constraints of density, entropy, and Lagrangian label conservation, which leads to a Clebsch representation for a generalized momentum and a generalized vector potential. The second action arises upon transformation to physical field variables, giving rise to a covariant bracket action principle, i.e., a variational principle in which constrained variations are generated by a degenerate Poisson bracket. Upon taking appropriate limits, the action principles lead to relativistic Hall MHD and well-known relativistic ideal MHD. For the first time, the Hamiltonian formulation of relativistic Hall MHD with electron thermal inertia (akin to [Comisso \textit{et al.}, Phys. Rev. Lett. {\bf 113}, 045001 (2014)] for the electron--positron plasma) is introduced. This thermal inertia effect allows for violation of the frozen-in magnetic flux condition in marked contrast to nonrelativistic Hall MHD that does satisfy the frozen-in condition. We also find violation of the frozen-in condition is accompanied by freezing-in of an alternative flux determined by a generalized vector potential. Finally, we derive a more general 3+1 Poisson bracket for nonrelativistic extended MHD, one that does not assume smallness of the electron ion mass ratio.

Charged particles in a magnetosphere are spontaneously attracted to a planet while increasing their kinetic energy via inward diffusion process. A constraint on particles' micro-scale adiabatic invariants restricts the class of motions available to the system, giving rise to a proper frame on which particle diffusion occurs. We investigate the inward diffusion process by numerical simulation of particles on constrained phase space. The results reveal the emergence of inhomogeneous density gradient and anisotropic heating, which is consistent with spacecraft observations, experimental observations, and the recently formulated diffusion model on the constrained phase space.