2 Electromagnetic linear dispersion relation

In a magnetised plasma with various species $s$, the linear dispersion relation for electromagnetic wave modes of (complex) frequency $\omega$ and wave mode $k$ that propagate parallel to a constant background magnetic field (here taken to be ${{\boldsymbol B}}_0 = B_0 \hat {\boldsymbol {x}}$) is given by (Schlickeiser Reference Schlickeiser2002)

(2.1)\begin{align} D^{{\pm}} & = \omega^2 - k^2c^2 + \sum_s \chi^{{\pm}}_s = 0, \quad {\rm where} \end{align}

(2.2)\begin{align} \chi^{{\pm}}_s & = \frac{ q^2_s }{m_s } \int {\rm d}^3u \frac{f_{s,0} (u_{{\parallel}} , u_{{\perp}}) }{\gamma} \left[ \frac{\omega -k \upsilon_{{\parallel}}}{k \upsilon_{{\parallel}}-\omega \pm \varOmega_s } - \frac{\upsilon_{{\perp}}^2 c^{{-}2}(k^2c^2-\omega^2 )} {2 (k \upsilon_{{\parallel}} - \omega \pm \varOmega_s )^2}\right]. \end{align}

Here, $\chi _s^{\pm }$ is the linear response for species $s$, which is characterised by charge $q_s$, mass $m_s$ and the equilibrium gyrotropic phase-space distribution, $f_{s,0}$. The light speed in vacuum is denoted by $c$ and the spatial part of the 4-velocity is ${{\boldsymbol u}} =\gamma \boldsymbol{\upsilon}$, where $u_{\parallel }$ and $u_{\perp }$ are the parallel and perpendicular velocities, respectively, which are defined with respect to the direction of the uniform background magnetic field ${{\boldsymbol B}}_0$. The magnitude of the velocity is defined such that $u^2 = u_{\parallel }^2 + u_{\perp }^2$ and $\gamma =\sqrt {1+u^2/c^2}$. The non-relativistic cyclotron frequency of species $s$ is $\varOmega _{s,0} = q_s B_0/ m_s$, and the relativistic one is $\varOmega _s = \varOmega _{s,0} / \gamma$.

For simplicity, we take $f_{s,0}$ to be a gyrotropic ring distribution with fixed parallel and perpendicular velocities for all particles of species $s$, i.e. $f_{s,0}(u_{\parallel },u_{\perp }) = n_s \delta (u_{\parallel } - \gamma \upsilon_{\rm dr}) \delta (u_{\perp } - \gamma \upsilon_{\perp })/ (2 \pi u_{\perp })$, where $n_s$ is a uniform density of the plasma species $s$, and all of its particles are drifting along ${{\boldsymbol B}}_0$ with constant speed $\upsilon_{{\rm dr},s}$ and have the same perpendicular velocity $\upsilon_{\perp,s}$. In this case, the linear response reduces to $\zeta _s^{\pm }$ given by (Holcomb & Spitkovsky Reference Holcomb and Spitkovsky2019; Weidl, Winske & Niemann Reference Weidl, Winske and Niemann2019; Shalaby et al. Reference Shalaby, Thomas and Pfrommer2021)

(2.3)\begin{equation} \chi_s^{{\pm}} \rightarrow \zeta_s^{{\pm}} ( \upsilon_{{\rm dr},s},\upsilon_{{\perp},s} ,n_s) =\frac{ \omega^2_s }{ \gamma} \left[\frac{\omega -k \upsilon_{{\rm dr},s} }{k \upsilon_{{\rm dr},s} -\omega \pm \varOmega_s } -\frac{\upsilon_{{\perp},s} c^{{-}2}(k^2c^2-\omega^2 )} {2 (k \upsilon_{{\rm dr},s} - \omega \pm \varOmega_s)^2}\right].\end{equation}

Here, $\omega _{s} = \sqrt { n_{ s} q^2_{s}/(m_{ s} \epsilon _0) }$ is the plasma frequency of species $s$ and $\epsilon _0$ is the permittivity of free space. The solutions of $D^{\pm }=0$ correspond to right/left polarisation states, and the corresponding first-order perturbed electric and magnetic fields obey $E_{y,1}(k,\omega ) \pm i E_{z,1}(k,\omega ) =0$ and $B_{y,1}(k,\omega ) \pm i B_{z,1}(k,\omega ) =0$, respectively. In the analysis below, the velocity distributions for both the background plasma and CR species are graphically illustrated in figure 1; a cold distribution for background electrons and ions and a gyrotropic distribution for CRs.

Figure 1.

Schematic representation of the velocity distributions adopted for our analytical calculation in this paper. We chose a cold distribution for background electrons and ions, and a gyrotropic distribution for CR electrons and ions. While a gyrotropic CR distribution is adopted here for analytical tractability, the nature of resonances that lead to the instabilities discussed in this work is independent of such a choice.

2.1 Normal modes of electron–ion plasma

The normal modes of an electron–ion plasma are well known and extensively studied in the literature (see, e.g. Boyd & Sanderson Reference Boyd and Sanderson2003). In the cold limit (also called the cold-hydrodynamic (HD) limit), one can obtain the plasma rest-frame wave modes by solving the following dispersion relation:

(2.4)\begin{equation} \omega^2- k^2 c^2 + \zeta_{\rm e}^{-} (0,0,n_0) + \zeta_{\rm i}^{-} (0,0,n_0) = 0, \end{equation}

where $\zeta _{\rm i(e)}$ denotes the contribution of the cold stationary ion (electron) species with a fixed uniform number density $n_0$. Therefore, for $k>0$, the real part of the solutions, $\omega _{\rm r}={\rm Re}(\omega )>0 ~ (<0)$, indicate that the sense of rotation of the magnetic eigenmodes is the same as that of electrons (ions) around the constant background magnetic field ${{\boldsymbol B}}_0$. The normal modes for such a case include light characteristics ($\omega _{\rm r} = k c$) in the short-wavelength regime, which turn into $\omega _{\rm r} \sim \omega _{\rm p}$ for small values of $k$, where $\omega _{\rm p} = \sqrt {\omega _{\rm i}^2+\omega _{\rm e}^2}$ is the total plasma frequency.

The other set of wave modes is the electron- and ion-cyclotron waves in the short-wavelength regime, $k d_{\rm i} \ll 1$, where $d_{\rm i}=c/\omega _{\rm i}$ is the ion skin depth. More precisely, the electron-cyclotron branch includes the forward-propagating Alfvén waves ($\omega _{\rm r} = k \upsilon_\mathrm {A}$) for $k d_{\rm i} \ll 1$, where $\upsilon_\mathrm {A} = B_0/\sqrt { \mu _0 n_{\rm i} m_{\rm i}} = \varOmega _{\rm i,0} c/\omega _{\rm i}$ is the Alfvén speed for ions. This branch turns into whistler waves for shorter wavelengths at around $k d_{\rm i} > 1$, which become electron-cyclotron waves with $\omega _{\rm r} = |\varOmega _{\rm e}|$ for $k d_{\rm e} = k d_{\rm i}/\sqrt {m_{\rm r}} >1$, where $d_{\rm e} = c /\omega _{\rm e}$ is the electron skin depth. The ion-cyclotron branch is conceptually simpler: it includes the backward-propagating Alfvén waves for $k d_{\rm i} \ll 1$ and the ion-cyclotron wave modes with $\omega _{\rm r} = - \varOmega _{\rm i}$ for $k d_{\rm i} \geq 1$. An example for the electron- and ion-cyclotron branches is depicted by the black solid lines in the top panel of figure 2.

Figure 2.

Solutions of the dispersion relation with low-density drifting CRs in the rest frame of the background plasma ((2.6)). We use the following parameters: $\upsilon_\mathrm {A} = 10^{-4} c$, $\alpha = n_{\rm cr}/n_0 = 10^{-6}$, $m_{\rm i}/m_{\rm e} = 36$, $v_{\perp, {\rm e}}=0$ and $\upsilon_{\perp,{\rm i}} = \upsilon_\mathrm {A}$. The solutions are, in general, eight complex values of $\omega$ at each wave mode $k$. In the top panel, we show the real values, $\omega _{\rm r}$, of three solutions: the CR ion cyclotron (coloured lines with different line styles, for three different parameter choices) as well as the background ion- and electron-cyclotron waves (black solid lines). We choose three values for the CR drift speed along ${{\boldsymbol B}}_0$ and indicate each case with a different colour and line style. In the bottom panel, we show the growth rate of the CR ion-cyclotron solution, which is also the fastest growth rate, $\varGamma$, at each $k$, i.e. it has the largest imaginary part of $\omega$ of all solutions. In the top panel, $\omega _{\rm r}$ represents the rate by which the electric and magnetic field perturbation vectors rotate. Here, $\omega _{\rm r}>0$ ($\omega _{\rm r}<0$) indicates the same sense of rotation as the gyro-motion of electrons (ions) around the large-scale background magnetic field ${{\boldsymbol B}}_0$. Vertical dashed lines are solutions of (2.7) and indicate values of the wave mode for which the CR ion-cyclotron wave in the background plasma rest frame is in exact resonance with the background cyclotron waves. These points exactly correspond to the locations of the peaks in the growth rate.

2.2 Resonances in the presence of drifting low-density CRs

Whether a drifting electron–ion plasma population (CR population) can excite resonant plasma instabilities is most easily seen by identifying the intersection points of background and CR wave modes in the $\omega$-$k$ plane: this yields solutions of wave modes for which the rotation rate of both background and CRs are in resonance. Thus, it is for these wave modes that we expect the largest energy exchange between CRs and the background plasma populations. In the following, we take a closer look at these resonant interactions in the case of CRs with a low number density, $n_{\rm cr}$, in comparison with the density of the background plasma $n_0$. We characterise this by the CR-to-background number density ratioFootnote ¹ $\alpha \equiv n_{\rm cr}/n_0 \ll 1$ and adopt a mean relative drift speed $v_{\rm dr}$ of the CR population along the large-scale background magnetic field. In the background frame, one can derive the expected CR modes by solving the dispersion relation

(2.5)\begin{equation} \omega^2- k^2 c^2 + \zeta_{\rm e}^{-} (v_{\rm dr},0, \alpha n_0) + \zeta_{\rm i}^{-} (v_{\rm dr},0, \alpha n_0) = 0. \end{equation}

As anticipated, one solution of (2.5) is the CR ion-cyclotron branch, which is comparable to that of the background, albeit with a Doppler-shifted rotation frequency, that is $\omega \rightarrow \omega + k \upsilon_{\rm dr}$. Additionally, in the ion-cyclotron branch, the backward-propagating Alfvén waves start to rotate at the Doppler-shifted ion-cyclotron frequency at $k d_{\rm cr,i} = k d_{\rm i}/ \sqrt {\alpha }\sim 1$. In other words, within the ion-cyclotron branch, wave modes rotate with the Doppler-shifted ion-cyclotron frequency at $k d_{\rm i} \gtrsim \sqrt {\alpha }$, which is much less than unity. This means that, as observed from the background, this branch mainly produces waves rotating with Doppler-shifted ion-cyclotron frequency $\omega _{\rm r} = - \varOmega _{\rm i} + k \upsilon_{\rm dr}$. We refer to this mode as the CR ion-cyclotron wave mode.

To verify that the individual wave modes of CRs and those of the background plasma are indeed the sum of the background and CR wave modes as seen in the background frame, we solve the full dispersion relation

(2.6)\begin{equation} \omega^2- k^2 c^2 + \zeta_{\rm e}^{-} (0,0, n_0) + \zeta_{\rm i}^{-} (0,0,n_0) + \zeta_{\rm e}^{-} (\upsilon_{\rm dr},\upsilon_{{\perp}, {\rm e}}, \alpha n_0) + \zeta_{\rm i}^{-} (\upsilon_{\rm dr},\upsilon_{{\perp},{\rm i}}, \alpha n_0) = 0. \end{equation}

Note that, in our set-up, the zeroth-order current due to CR ions is compensated by CR electrons drifting at the same speed. Here, we have added the possibility for CR populations to have a non-zero perpendicular velocity. That is, CRs are distributed uniformly on a ring in the perpendicular velocity space, characterised by a radius of $\upsilon_{\perp }$. Solutions for the full dispersion relation are shown in figure 2, where we use a reduced mass ratio $m_{\rm r} = m_{\rm i}/m_{\rm e} = 36$ for visual purposes to reduce the separation between electron- and ion-cyclotron frequencies; $\sqrt {m_{\rm r}}/2 = 3$. We also adopt an ion Alfvén speed $\upsilon_\mathrm {A}$ at $10^{-4}c$, the density ratio $\alpha = n_{\rm cr}/n_0 = 10^{-6}$, $\upsilon_{\perp, {\rm e}}=0$ and $\upsilon_{\perp,{\rm i}} = \upsilon_\mathrm {A}$. We use different values of $\upsilon_{\rm dr}/\upsilon_\mathrm {A} \approx \{3.7, 3, 2.7\}$, such that various interesting types of resonances are expected in the solutions of the dispersion relation. We note that this choice of parameters enables us to show the physics of resonance on a linear scale in the background plasma frame. In § 2.3, we adopt more physically motivated parameters and show the resulting solutions of the dispersion relation in the CR drift frame.

For all values of $\upsilon_{\rm dr}$, we obtain the same exact background wave modes as described in the previous section, which are shown as black lines in the top panel of figure 2. Differing values of $\upsilon_{\rm dr}$ result in various CR ion-cyclotron wave modes, which are shown as different colours and which exactly agree with $\omega _{\rm r} = - \varOmega _{\rm i} + k \upsilon_{\rm dr}$ expected for the various cases. In the bottom panel of figure 2, we show the fastest growth rate obtained by solving the full dispersion relation in (2.6). The real part of the most unstable modes are the CR ion-cyclotron waves that we show with different colours and line styles in the top panel of figure 2. That is, in the unstable regions, only the rotating and propagating electromagnetic waves that are supported by the gyrating and drifting CRs are exponentially growing. We note that in regions where the fastest growth rate is zero, any branch can be chosen, and we select the CR ion-cyclotron wave mode in this case.

To find the locations where the CR ion-cyclotron wave modes are in resonance with the background wave modes, we solve

(2.7)\begin{equation} [\omega^2- k^2 c^2 + \zeta_{\rm e}^{-} (0,0,n_0) + \zeta_{\rm i}^{-} (0,0,n_0)]_{\omega ={-} \varOmega_{\rm i} + k \upsilon_{\rm dr}} = 0, \end{equation}

for wave mode $k$. That is, the solutions are those wave modes for which the CR ion-cyclotron branch intersects the background wave modes in the cold-HD limit. This gives a fourth-order polynomial which, in general, has four solutions for $k$ for any value $\upsilon_{\rm dr}/\upsilon_\mathrm {A}$, albeit these solutions may be degenerate.

For the case of $\upsilon_{\rm dr}/\upsilon_\mathrm {A} = 2.7 < \sqrt {m_{\rm r}}/2$ (solid-red curves in figure 2), there exist four locations where the CR ion-cyclotron wave modes are in resonance with the background wave modes, i.e. (2.7) has four distinct roots. These expected locations for resonances are indicated with red-dashed vertical lines in the top and bottom panels of figure 2: there are two resonances for $k d_{\rm i} <1$ and two resonances for $k d_{\rm i} >1$. The two solutions for $k d_{\rm i} <1$ represent the intersections of the CR ion-cyclotron wave mode with both forward- and backward-propagating Alfvén waves, and thus resonance occurs approximately at $k d_{\rm i} = \upsilon_\mathrm {A}/(\upsilon_{\rm dr} \pm \upsilon_\mathrm {A})$. This corresponds to the two resonant peaks in the growth rate due to the gyro-resonant (streaming) instability (Kulsrud & Pearce Reference Kulsrud and Pearce1969). The other two resonances at smaller scales at $k d_{\rm i} >1$ correspond to the two peaks in the growth rate due to the intermediate-scale instability (Shalaby et al. Reference Shalaby, Thomas and Pfrommer2021). In the next section, we show the growth rates when using realistic values of the mass ratio, for which the fastest growth rate due to the intermediate-scale instability is more than an order of magnitude greater than that at the gyroscales.

In figure 2, the purple dash-dotted curves represent the case of $\upsilon_{\rm dr}/\upsilon_\mathrm {A} \sim \sqrt {m_{\rm r}}/2$. We see that the two resonant peaks of the intermediate-scale instability at $k d_{\rm i} >1$ merge into a single resonance. As in the previous case, the two peaks of the gyroscale instability for $k d_{\rm i} <1$ remain, but move closer together. As $\upsilon_{\rm dr}/\upsilon_\mathrm {A}$ grows larger than $\sqrt {m_{\rm r}}/2$, the resonances at $k d_{\rm i} >1$ disappear, and thus, the intermediate-scale instability is no longer able to drive wave modes with $k d_{\rm i} > 1$ unstable (blue dashed curves in figure 2). That is, $\upsilon_{\rm dr}/\upsilon_\mathrm {A} \leq \sqrt {m_{\rm r}}/2$ is a condition for the possibility of resonance between CR ion-cyclotron wave modes and background wave modes at $k d_{\rm i} >1$, which is also the condition for the intermediate-scale instability.Footnote ² The numerical identification of this instability condition was previously conducted by Shalaby et al. (Reference Shalaby, Thomas and Pfrommer2021).

It is important to acknowledge that, in a realistic environment, background plasmas possess a finite temperature. Consequently, when incorporating this into the dispersion relation, it can influence wave modes occurring on scales smaller than the ion Debye length, which is much smaller in comparison with the ion skin depth for non-relativistic plasmas (Reville et al. Reference Reville, Kirk, Duffy and O'Sullivan2008; Zweibel & Everett Reference Zweibel and Everett2010). Thus, it does not have any impact on the instabilities discussed in this paper. Furthermore, discussions and analysis of simulations in Shalaby et al. (Reference Shalaby, Thomas and Pfrommer2021) and Lemmerz et al. (Reference Lemmerz, Shalaby, Thomas and Pfrommer2023) reveal that ion-cyclotron thermal damping has negligible effects on the driven wave modes.

2.3 Instabilities in the rest frame of the CRs

In this section, we demonstrate that the intermediate-scale instability indeed drives CR comoving ion-cyclotron wave modes as found in Shalaby et al. (Reference Shalaby, Thomas and Pfrommer2021). To this end, we present the solution of the dispersion relation in the CR rest frame, where the background plasma is drifting with velocity $-\upsilon_{\rm dr}$ and hence drifting anti-parallel with respect to the background magnetic field ${{\boldsymbol B}}_0$. The dispersion relation that we solve is given by

(2.8)\begin{equation} \omega^2\!- k^2 c^2 \!+ \zeta_{\rm e}^{-} (-\upsilon_{\rm dr},0, n_0) + \zeta_{\rm i}^{-} (-\upsilon_{\rm dr},0,n_0) + \zeta_{\rm e}^{-} (0,\upsilon_{{\perp}, {\rm e}}, \alpha n_0) + \zeta_{\rm i}^{-} (0,\upsilon_{{\perp}, {\rm i}}, \alpha n_0) = 0. \end{equation}

Figure 3 shows some solutions of (2.8) for $\upsilon_\mathrm {A} = 10^{-4} c$, $\alpha = n_{\rm cr}/n_0 = 10^{-6}$,Footnote ³ $m_{\rm r} = m_{\rm i}/m_{\rm e} = 1836$, $v_{\perp, {\rm e}}=0$ and $\upsilon_{\perp, {\rm i}} = \upsilon_\mathrm {A}$, along with various background plasma drift speeds such that $\upsilon_{\rm dr}/\upsilon_\mathrm {A} \sim \{22.28, 21.45 , 19.281 \} \sim \{1.05, 1, 0.9 \} \sqrt {m_{\rm r}}/2$. In the top panel of that figure, the CR ion-cyclotron wave modes are shown as a black line for different values of $\upsilon_{\rm dr}$. Meanwhile, the background cyclotron wave modes are rotated differently depending on the drift speed, leading to intersections at both long ($kd_{\rm i} <1$) and short ($kd_{\rm i} >1$) wavelengths for $\upsilon_{\rm dr}/\upsilon_\mathrm {A} \leq \sqrt {m_{\rm r}}/2$. These intersections lead to instabilities at these wave modes. The bottom panel of figure 3 shows the growth rates of the ion-cyclotron wave modes, which are also the fastest-growing modes. This demonstrates that the driven modes are ion-cyclotron waves at all wavelengths in the rest frame of CRs. If the condition for the intermediate-scale instability is not fulfilled ($\upsilon_{\rm dr}/\upsilon_\mathrm {A} > \sqrt {m_{\rm r}}/2$), the instability no longer operates and short-wavelengths modes are stable.

Figure 3.

Solutions of the dispersion relation of low-density CRs in the rest frame of the CRs, i.e. using (2.8), with the following parameters: $\upsilon_\mathrm {A} = 10^{-4} c$, $\alpha = n_{\rm cr}/n_0 = 10^{-6}$, $m_{\rm i}/m_{\rm e} = 1836$, $v_{\perp,{\rm e}}=0$ and $\upsilon_{\perp, {\rm i}} = \upsilon_\mathrm {A}$. We vary the relative drift speed (the background plasma is drifting with $\upsilon_{\rm dr}$ anti-parallel to background magnetic field ${{\boldsymbol B}}_0$) to obtain similar cases as in figure 2 but with a realistic ion-to-electron mass ratio. In the top panel, we show the rotation rate of wave modes, with the solid-black line representing the CR ion-cyclotron wave mode which is the same for all cases. The background wave modes are the same in the background rest frame, however, as seen in the rest frame of the CRs, the rotation and growth rates of the background waves are different and are thus indicated with various colours and line styles. In the bottom panel, we show the growth rates of the CR ion-cyclotron waves in all cases, and we find that these are, also, the fastest-growing rates for the instability. This shows that, in the case of realistic $m_{\rm r}$, the dominance of intermediate-scale growth rate compared with that at the gyroscale ($kd_{\rm i}<1$) is much more pronounced and that the growth rate for the forward-propagating Alfvén waves is also larger in comparison with the backward-propagating wave.

The solutions in figure 3 show that, for realistic values of $m_{\rm r}$, the intermediate-scale instability growth rates significantly dominate over gyroscale growth rates. This dominance is even more pronounced at larger pitch angles because the growth rate at the peaks of the intermediate-scale instability is proportional to $(\upsilon_{\perp,{\rm i}}/\upsilon_\mathrm {A})^{2/3}$ (Shalaby et al. Reference Shalaby, Thomas and Pfrommer2021). Additionally, it is important to note that the growth rate is higher for the resonance with forward-propagating Alfvén waves at the gyroscale ($kd_{\rm i}<1$) compared with that with the backward-propagating waves.

3 Instabilities with approximate background plasma descriptions

In this section, we revisit the dispersion relation in the background plasma rest frame to investigate the effect of approximating the background plasma description on the nature of the emerging instabilities. One commonly used approximation for the background plasma dispersion assumes the MHD dispersion relation, where it is assumed that $\omega \ll \varOmega _{\rm i} \ll |\varOmega _{\rm e}|$. In this limit, the third and fourth terms in (2.6) are reduced to $\omega ^2 c^2/\upsilon_\mathrm {A}^2$ and the first term ($\omega ^2$) is neglected, due to the fact that $\upsilon_\mathrm {A}^2 \ll c^2 \Rightarrow \omega ^2 \ll \omega ^2 c^2/\upsilon_\mathrm {A}^2$ . The dispersion relation is for

(3.1)\begin{equation} {\rm MHD:}\ \frac{\omega^2 c^2}{\upsilon_\mathrm{A}^2}- k^2 c^2 + \zeta_{\rm e}^{-} (\upsilon_{\rm dr},\upsilon_{{\perp},{\rm e}}, \alpha n_0) + \zeta_{\rm i}^{-} (\upsilon_{\rm dr},\upsilon_{{\perp}, {\rm i}}, \alpha n_0) = 0. \end{equation}

This approximation is typical in works aiming at computing various types of CR-driven instabilities (see, e.g. Zweibel Reference Zweibel2003; Bell Reference Bell2004; Amato & Blasi Reference Amato and Blasi2009; Bai et al. Reference Bai, Ostriker, Plotnikov and Stone2019). Another possible approximation involves considering the impact of finite but small values of $k d_{\rm i}$, which are wave modes comparable to the ion skin depth, as is typically done in the Hall-MHD approximation. In this case, the dispersion relation approximates the behaviour of Alfvén and whistler waves found in the cold HD dispersion relation. Within the Hall-MHD approximation, the background plasma contributions (see, e.g. § 14.4.4 of Goedbloed, Keppens & Poedts Reference Goedbloed, Keppens and Poedts2010) are reduced to $\omega ^2 c^2/[\upsilon_\mathrm {A}^2 (1+ k d_{\rm i})^2]$, leading to the following modified dispersion relation:

(3.2)\begin{equation} \text{Hall-MHD:} \ \frac{\omega^2 c^2}{\upsilon_\mathrm{A}^2 (1+ k d_{\rm i})^2 }- k^2 c^2 + \zeta_{\rm e}^{-} (\upsilon_{\rm dr},\upsilon_{{\perp}, {\rm e}}, \alpha n_0) + \zeta_{\rm i}^{-} (\upsilon_{\rm dr},\upsilon_{{\perp},{\rm i}}, \alpha n_0) = 0.\end{equation}

In figure 4, we show the solutions of the dispersion relations with various background plasma assumptions, i.e. the solutions for ((2.6), (3.1) and (3.2)) using the following parameters: $\upsilon_\mathrm {A} = 10^{-4} c$, $\alpha = n_{\rm cr}/n_0 = 10^{-6}$, $m_{\rm i}/m_{\rm e} = 1836$, $v_{\perp, {\rm e}}=0$, $\upsilon_{\perp, {\rm i}} = \upsilon_\mathrm {A}$ and $\upsilon_{\rm dr}/\upsilon_\mathrm {A} = 0.98 \sqrt {m_{\rm r}}/2 \approx 20.99$. In all cases, we fix the CR drift speed, resulting in the same CR ion-cyclotron wave mode (shown as a black line in the top panel) that follows the expected dispersion relation $\omega _{\rm r} = k \upsilon_\mathrm {A} - \varOmega _{\rm i}$. When using the MHD background approximation, the background wave modes follow the expected dispersion relation $\omega _{\rm r} = \pm k \upsilon_\mathrm {A}$ (shown as orange dash-dotted curves in the top panel). In the case of Hall-MHD, the dispersion of background wave modes, as expected, follows $\omega _{\rm r} = \pm k \upsilon_\mathrm {A} (1+k d_{\rm i})$; shown as green dashed curves in the top panel.

Figure 4.

Solutions of the dispersion relation for various approximations of the background plasma using the following parameters: $\upsilon_\mathrm {A} = 10^{-4} c$, $\alpha = n_{\rm cr}/n_0 = 10^{-6}$, $m_{\rm i}/m_{\rm e\ }= 1836$, $v_{\perp, {\rm e}}=0$, $\upsilon_{\perp,{\rm i}} = \upsilon_\mathrm {A}$ and $\upsilon_{\rm dr}/\upsilon_\mathrm {A} = 0.98 \sqrt {m_{\rm r}}/2 \sim 20.99$. The top panel shows the rotation rate of various wave modes, while the bottom panel shows the growth rate of the corresponding unstable CR ion-cyclotron wave mode, which has the fastest growth rate. The black line in the top panel represents the CR ion-cyclotron wave mode, which follows the same dispersion in all cases, namely $\omega _{\rm r} = k \upsilon_{\rm dr} - \varOmega _{\rm i}$. In the MHD case, the background wave modes follow a dispersion relation of $\omega = \pm k \upsilon_\mathrm {A}$ (orange curves in the top panel), while in the Hall-MHD case, the background waves follow a dispersion relation of $\omega = \pm k \upsilon_\mathrm {A} (1+k d_{\rm i})$. The bottom panel reveals that the growth rates of the CR ion-cyclotron wave mode at the gyroscale ($k d_{\rm i}<1$) are almost identical for different approximations of the background plasma. However, at intermediate scales where $k d_{\rm i}>1$, the MHD background approximation fails to capture the fastest-growing instability, namely the intermediate-scale instability with growth rates shown by the red curve. In contrast, the Hall-MHD approximation reproduces the first peak of the dominant instability growth rate, although at a reduced rate and at a larger wavelength. For all wavelengths shorter than this peak, the use Hall-MHD approximation leads to wrong growth rates. Vertical lines indicate the predicted intersection points of CR ion-cyclotron waves and background wave modes for the different approximations.

In the bottom panel of figure 4, we show the growth rates of the CR ion-cyclotron waves (which are the fastest growth rates) for different approximations of the background plasma. At large (gyro) scales, all approximations produce similar growth rates peaking at the same wave modes, indicating that using any of these approximations for the background plasma leads to correctly capturing the gyroscale instability. However, at smaller scales, i.e. scales where the intermediate-scale instability operates ($kd_{\rm i}>1$), the MHD approximation of the background plasma dispersion wrongly predicts complete stability at these scales. In the case of Hall-MHD, the growth of the first peak of the intermediate-scale instability is approximately reproduced but at a longer wavelength. For shorter wavelengths, the use of the Hall-MHD approximation predicts completely wrong growth rates compared with those found when using the cold-HD dispersion for the background plasma. An additional worrisome implication arising from the use of the Hall-MHD approximation is its incorrect prediction that shorter-wavelength modes are universally unstable, even in cases where the condition for driving the intermediate-scale instability is not satisfied.

4 Nature of CR scattering at different resonant scales

So far, we have focused on the influence of gyrotropic CR populations in driving perpendicular electromagnetic perturbations. In the absence of such perturbations, particles follow trajectories characterised by a constant drift speed along the direction of the magnetic field vector, while simultaneously gyrating at their gyrofrequency. In this section, our goal is to evaluate the effects of these perturbations on the particle trajectories, specifically how particles are scattered by the induced parallel electromagnetic perturbations. To accomplish this, we examine the Lorentz force acting on a particle with velocity $\boldsymbol{\upsilon} =\{ \upsilon_{x} , \upsilon_{y}, \upsilon_{z} \}$, charge $q_\mathrm {s}$ and mass $m_\mathrm {s}$ caused by these perturbations. In Fourier space, along the particle trajectory, the momentum equation is

(4.1)\begin{equation} \frac{\mathrm{d} \gamma \boldsymbol{\upsilon}}{\mathrm{d}t} = \frac{q_\mathrm{s}}{m_\mathrm{s}} ( \delta {{{\boldsymbol E}}}^k + \boldsymbol{\upsilon} \boldsymbol{\times} \delta {{{\boldsymbol B}}}^k ) =\frac{q_\mathrm{s}}{m_\mathrm{s}} \begin{bmatrix} \upsilon_{y} \delta B^{k}_{z} - \upsilon_{z} \delta B^{k}_{y}\\ -(\upsilon_{x}-\upsilon_{\rm ph}) \delta B^{k}_{z}\\ (\upsilon_{x}-\upsilon_{\rm ph})\delta B^{k}_{y} \end{bmatrix}. \end{equation}

Here, we use the fact that ${{\boldsymbol k}} \parallel {{\boldsymbol B}}_0$ and both vectors are aligned with $\hat {\boldsymbol {x}}$. Furthermore, we utilise the relation ${{\boldsymbol k}} \boldsymbol {\times } \delta {{\boldsymbol E}}^k = \omega \delta {{\boldsymbol B}}^k$ to find $\delta {{\boldsymbol E}}^k = - \boldsymbol{\upsilon} _{\rm ph} \boldsymbol {\times } \delta {{\boldsymbol B}}^k$. This indicates that particle scattering in the parallel direction is significantly influenced by the phase difference between the particle's perpendicular velocities and the magnetic field perturbations. On the other hand, scattering of particles in the perpendicular direction crucially depends on the difference between the relative drift speed of the particles and the waves.

In previous sections, we found that the propagation of CRs destabilise waves at different wavelengths. The real frequency of the unstable wave modes, in the rest frame of the background, is always given by $\omega _{\rm r} = k \upsilon_{\rm dr} - \varOmega _{\rm i}$. Consequently, the phase velocity of the unstable waves can be expressed as $\upsilon_{\rm ph} = \upsilon_{\rm dr} - \upsilon_\mathrm {A}/(k d_{\rm i})$. At the gyroscale ($k d_{\rm i} < 1$), wave growth peaks at wave modes where $k d_{\rm i} \approx \upsilon_\mathrm {A}/(\upsilon_{\rm dr} \mp \upsilon_\mathrm {A})$. This leads to a phase velocity of approximately $\upsilon_{\rm ph} \approx \pm \upsilon_\mathrm {A}$ for forward ($+$) and backward ($-$) propagating Alfvén waves. On the other hand, the fastest growth due to the intermediate-scale instability occurs for $k d_{\rm i} > 1$. Since typically $\upsilon_{\rm dr} \gg \upsilon_\mathrm {A}$, the phase velocity of the driven unstable modes is approximately $\upsilon_{\rm ph} \approx \upsilon_{\rm dr}$. To summarise, the phase speed of the growing wave modes can be expressed as follows:

(4.2)\begin{equation} \upsilon_{\rm ph} = \left\{ \begin{array}{ll} \pm \upsilon_\mathrm{A}, & kd_{\rm i} \ll 1, \\ \upsilon_{\rm dr}, & kd_{\rm i} > 1. \end{array} \right. \end{equation}

Therefore, although both the gyro and intermediate-scale instabilities are resonant instabilities, they lead to perturbations with distinct phase speeds. Due to the significant disparity in phase velocity between these two scales, the scattering of CR ions by these electromagnetic perturbations exhibits notable differences, as can be seen from (4.1). The wave modes driven by the intermediate-scale instability scatter particles that drive them only in the direction parallel to ${{\boldsymbol B}}_0$, since $\upsilon_{x} = \upsilon_{\rm dr} \sim \upsilon_{\rm ph}$. However, particles with different parallel drift speeds can be scattered both in parallel and perpendicular directions due to wave modes driven by the intermediate-scale instability. On the other hand, scattering occurs in both parallel and perpendicular directions at the gyroscale. It can be demonstrated from (4.1) that the gyroscale waves result in energy-conserving scattering in the frame of the driven Alfvén waves (Shalaby et al. Reference Shalaby, Thomas and Pfrommer2021).

5 Summary

This paper examines the physics of resonant instabilities driven by CR ions with a gyrotropic momentum distribution. These instabilities occur most rapidly when a resonance between the Doppler-shifted background and CR wave modes occurs. The Doppler shift arises from the relative drift between the background and CR plasma. The relative drift sets the location of resonances and hence the most unstable wavelengths of various instabilities. It leads to two important resonant instabilities: the gyroscale instability at large scales (Kulsrud & Pearce Reference Kulsrud and Pearce1969) and the recently found intermediate-scale instability (Shalaby et al. Reference Shalaby, Thomas and Pfrommer2021).

The gyroscale instability has peak growth rates when the CR ion-cyclotron wave mode resonates with both, forward- and backward-propagating Alfvén waves of the background plasma. In the background frame, the resonances occur when $- \varOmega _{\rm i} + k \upsilon_{\rm dr} \rightarrow \pm k \upsilon_\mathrm {A}$, which results in the most unstable wave modes at $k_{\rm g,i}^{\pm }d_{\rm i} = \upsilon_\mathrm {A} /(\upsilon_{\rm dr} \mp \upsilon_\mathrm {A})$. From figure 3, it is evident that resonance with the forward-propagating Alfvén wave leads to faster growth at larger wave modes $k_{\rm g,i}^+$, thus it is called the ion gyroscale. Since $\upsilon_{\rm dr} \gg \upsilon_\mathrm {A}$, both growth peaks at the gyroscales occur for wavelengths larger than the ion skin depth, i.e. $k_{\rm g,i}^{\pm }d_{\rm i} <1$. Therefore, in the linear regime, the gyroscale instability can be accurately described when approximating the electron–ion background plasma using either MHD or Hall-MHD approximations (see figure 4).

For an electron–ion background plasma, additional resonances occur at shorter wavelengths, where $k d_{\rm i} > 1$. These resonances give rise to the intermediate-scale instability that destabilise wave modes between the ion gyroscale, $k_{\rm g,i}^+$ and the electron gyroscale, $k_{\rm g,e} = m_{\rm r} k_{\rm g,i}^+$, where $m_{\rm r}$ represents the ion-to-electron mass ratio. At these resonant scales, the peak growth rates of the intermediate-scale modes are significantly larger compared with those at the gyroscales. This establishes the intermediate-scale instability as the fastest instability in the linear regime of the resonant CR-driven instabilities. Moreover, the dominance of the peak growth of the intermediate-scale instability is further amplified when CR ions possess larger pitch angles, corresponding to a larger perpendicular velocity ($\upsilon_{\perp }$) (see table 1 of Shalaby et al. Reference Shalaby, Thomas and Pfrommer2021).

As depicted in figure 4, the utilisation of the MHD approximation for describing the linear response of the background plasma hinders the occurrence of resonances at short wavelengths. Consequently, such an approximation suppresses the dominant, intermediate-scale instability. Conversely, employing the Hall-MHD approximation for the background plasma captures a resonance at $k d_{\rm i} > 1$, and, while the fastest growth rate associated with the intermediate-scale instability is approximately captured, it is associated with a wrong wavenumber. Additionally, this approximation erroneously predicts the presence of intermediate-scale instability even when the conditions for the instability are not met.

We argue in § 1 that background temperatures are unlikely to impact the growth of the instabilities addressed in this paper. Supporting evidence for this assertion can be found in the simulations conducted by Shalaby et al. (Reference Shalaby, Thomas and Pfrommer2021), where the background plasma was characterised by high temperatures, yet exhibited excellent agreement with the growth rates predicted from the dispersion relation assuming a cold background plasma. On the other hand, it is natural to contemplate whether these instabilities, particularly the newly discovered intermediate-scale instability, persist under different and potentially more realistic velocity distributions for CR ions. The fundamental explanation for the origin of these instabilities lies in the resonance between the ion-cyclotron wave modes of CR ions and the background wave modes. Consequently, any velocity distribution that supports CR ion-cyclotron modes while adhering to the instability conditions will excite this instability. This argument is supported by circumstantial evidence from particle-in-cell simulations of shocks conducted by Shalaby et al. (Reference Shalaby, Lemmerz, Thomas and Pfrommer2022), where the instability is clearly driven by CRs with a thermal velocity distribution.

Thus, while our choice of the CR ion distribution ensures the ease of repeatability in our analytical calculations, it does not imply that the instability is exclusively associated with such a choice of CR ion velocity distribution. The choice of a particular CR ion velocity distribution could impact the growth rates of different instabilities. However, qualitative considerations show that it is highly likely that the intermediate-scale instabilities will remain dominant even in these cases. An analytical demonstration of this is, however, deferred to future studies.

6 Outlook

The presence of the intermediate-scale instability is vital for the efficiency of electron acceleration in parallel electron–ion non-relativistic shocks as seen in fully kinetic particle-in-cell simulations (Shalaby et al. Reference Shalaby, Broderick, Chang, Pfrommer, Lamberts and Puchwein2017a, Reference Shalaby, Thomas and Pfrommer2021). Moreover, when the condition for the instability is not met in simulations, a notable decrease in the efficiency of electron acceleration is observed (Shalaby et al. Reference Shalaby, Lemmerz, Thomas and Pfrommer2022). That is, the significance of this instability extends beyond its potential to regulate the transport of CRs in various astrophysical scenarios. It could substantially influence particle injection and acceleration processes occurring at shocks, the escape of CRs from their sources into the interstellar medium and CR-driven galactic winds. This novel understanding of the fundamental physics underlying CR-driven resonant plasma instabilities will contribute to elucidate the critical role played by CRs in many astrophysical environments.