2. Pressure anisotropy in the limits of fast transit time and no pitch angle scattering

In the present analysis we consider electrons sufficiently energetic that compared with $e\boldsymbol {v}\times \boldsymbol {B}$ the force of the electric field becomes unimportant to their single-pass guiding centre orbit motion. This restricts the study to superthermal electrons for which ${\mathcal {E}} \gg e \varPhi$, where $\varPhi \simeq T_e/e$ is the electrostatic potential. Given the low electron mass, the electron thermal motion is often sufficiently fast that each electron transit occurs on a time scale that is much faster than the time scale characterizing the evolution of the magnetic field. Considering this fast transit-time limit, the first and second adiabatic invariants, $\mu = m{v_\perp }^2/(2B)$ and $J=\oint v_{\parallel } \,\textrm {d}l$, respectively, are conserved quantities. We also assume a two-dimensional periodic system, without gradients in the $y$ direction, i.e. $\partial / \partial y=0$. In the drift-kinetic limit (assuming a low electron mass) the electrons then strictly follow the contours of constant magnetic flux, $\varPsi = A_y(x,z)$. Here $A_y(x,z)$ is the $y$ component of the magnetic vector potential, and assuming $B_y=0$, the magnetic field is given by ${\boldsymbol {B}}= \boldsymbol {\nabla } A_y \times {\boldsymbol {e}}_y$, where ${\boldsymbol {e}}_y$ is the unit vector in the $y$ direction.

The framework is readily generalized to other geometries as long as only a single local maximum is present in $B$ within the periodic domain considered. As shown by Egedal et al. (Reference Egedal, Schroeder and Lichko2021), multiple local maxima in $B$ are associated with separate locally trapped populations and in a dynamically changing configuration orbit transitions can lead to energization through $v_{\parallel }$ mixing. This effect of energization by mixing of separate trapped populations is not addressed here.

To illustrate how anisotropic features in $f(\boldsymbol {x},\boldsymbol {v})$ are generated during magnetic perturbations, we consider the simple compressional perturbation shown in figure 1, where $x$ represents the position along the flux tube. We introduce ${\tilde {B}}_m(t)$ as the maximum field strength observed along the flux tube, normalized by the background field $B_0$. Thus, for sinusoidal and periodic perturbations of amplitude $\Delta B$ and angular frequency $\omega$, we have

(2.1)\begin{equation} {\tilde{B}}_m(t)=1 + (\Delta B/B_0) |\sin(\omega t)|. \end{equation}

The spatial and temporal evolution of the normalized magnetic field is then described by

(2.2)\begin{equation} {\tilde{B}}(x,t)= 1-({\tilde{B}}_m(t)-1) \cos(k_{{\parallel}} x). \end{equation}

As indicated in figure 1, we consider the section of the flux tube with $-L/2 \leq x\leq L/2$, where $L=2{\rm \pi} /k_{\parallel }$, and with values of ${\tilde {B}}$ characterized by the range $2 - {\tilde {B}}_m \leq \tilde {B} \leq {\tilde {B}}_m$.

Figure 1. A trapped (green) and a passing (black) electron orbit within a compressional magnetic perturbation.

The limit of negligible pitch angle scattering is conveniently explored using the adiabatic invariants $\mu$ and $J$. For the present geometry the orbits can be categorized as either trapped or passing as illustrated in figure 1, and we define the action integral as

(2.3)\begin{equation} J= 4\int_0^{L_b} v_{{\parallel}} \,{\textrm{d} x}', \end{equation}

where for trapped particles $L_b$ is the location of the particle turning point, while $L_b=L/2$ for passing electrons. For trapped orbits the factor of 4 in front of the integral corresponds to the four sections of the full orbits with alternating signs of $xv_{\parallel }$. For the case of passing electrons the factor of 4 results when counting contributions from orbits with both signs of $v_{\parallel }$ (but otherwise identical values of $\mu$ and ${\mathcal {E}}$), making the definition of $J$ in (2.3) continuous across the trapped–passing boundary.

We first consider the limit of negligible pitch angle scattering, and assume that the initial distribution function $f=f_{0}(v^2)$ is isotropic along the uniform flux tube ${\tilde {B}}=1$ observed at $t=0$. It is now useful to introduce the normalized action integral and the so-called pitch angle variable, respectively defined by

(2.4a,b)\begin{equation} j= \frac{1}{2Lv}J,\quad \varLambda = \frac{\mu B_0}{{\mathcal{E}}}. \end{equation}

For $t=0$, where the flux tube is uniform with $B=B_0$, we find that $j=v_{\parallel }/v$ and $\varLambda ={v_\perp }^2/v^2$, such that for this particular time $v^2= v^2(j^2+\varLambda )$. The distribution function during the evolution of the magnetic perturbation can then simply be expressed as

(2.5)\begin{equation} f(x,v_{{\parallel}},{v_\perp},t)=f_0(v^2(j^2+\varLambda)). \end{equation}

This result follows directly from Jean's theorem (Jeans Reference Jeans1915) stating that $f=g(\mu,J)$ is a solution to the kinetic equation for an arbitrary function $g$, where we note that $v^2(j^2+\varLambda )= (J/(2L))^2+2\mu B_0/m_e$ is a simple function of these invariants, and that (2.5) reproduces the initial condition of $f=f_0(v^2)$ for $t=0$.

For a given point in the configuration space $(x,v_{\parallel },{v_\perp },t)$ the values of $\varLambda$ and $j$ are readily evaluated. For this, we notice that

(2.6)\begin{equation} j(\varLambda,t) = \frac{2}{L} \int_0^{L_b} \frac{v_{{\parallel}}}{v} {\textrm{d} x}' = \frac{2}{L} \int_0^{L_b} \sqrt{1 - \varLambda {\tilde{B}}(x',t)} \, {\textrm{d} x}' . \end{equation}

Here, for any time $t$, the bounce points of trapped electrons are obtained from the condition $\varLambda {\tilde {B}}(L_b,t)=1$, such that $j$ is independent of $v$ and a function of only $(\varLambda,t)$. With $\varLambda$ determined by the local values of $v_{\parallel }$, ${v_\perp }$ and $B(x,t)$, numerical integration of (2.6) then provides all the information required for evaluation of $f(x,v_{\parallel },{v_\perp },t)$ in (2.5).

The function $j(\varLambda,t)$ in (2.6) is also fundamental to the pumping model to be developed in the following sections, and figure 2(a) shows $j(\varLambda )$ computed for ${\tilde {B}}_m= 1$, $1.2$ and $1.5$. A value of $\varLambda =0$ is characteristic of passing orbits with $|v_{\parallel }/v|=1$, and given the normalization in (2.6) we have $j=1$. Because $B_{\rm {min}} /B_0= 2- {\tilde {B}}_m(t)$ is a function of time, the maximum value of $\varLambda$, which we denote $\varLambda _m=B_0/B_{\rm {min}}$, also changes in time. For example, for ${\tilde {B}}_m=1$, we find $j=\sqrt {1-\varLambda }$, and we have $\varLambda _m=1$. Meanwhile, for the green curve in figure 2(a) for ${\tilde {B}}_m=1.5$, we observe $\varLambda _m=2$. For any geometry, the most deeply trapped electrons are those with $\varLambda =\varLambda _m$, corresponding to $v_{\parallel }/v=0$, such that $j(\varLambda _m,t)=0$ for all values of $t$.

Figure 2. (a) The function $j(\varLambda )$ in (2.6) plotted as a function of $\varLambda$ and evaluated for three separate values of ${\tilde {B}}_m$. (b) The relative particle energization $v^2/v_0^2=1/(j^2+\varLambda )$. (c,d) Examples of local electron distributions predicted by (2.5) both computed for $B_m/B_0=1.5$, and a Maxwellian $f_0$ with thermal speed $v_t=1$. In (c) the local position is at the point of minimum field strength, $B/B_0=0.5$, whereas $B/B_0=1.25$ for the distribution in (d).

The increase in particle energy during the evolution of the magnetic perturbation can also be deduced from the expression in (2.5). By Liouville mapping, the argument of $f_0$ must be $v_0^2$, such that $v_0^2=v^2(j^2+\varLambda )$. It follows that the relative energization is ${\mathcal {E}}/{\mathcal {E}}_0= v^2/v_0^2 = 1/(j^2+\varLambda )$, and is in figure 2(b) plotted as a function of $\varLambda$ for the same three values of ${\tilde {B}}_m$ considered in figure 2(a). For ${\tilde {B}}_m=1$ there is naturally no heating and $v^2/v_0^2=1$. For ${\tilde {B}}_m>1$, depending on the value of $\varLambda$ both heating and cooling are observed, and as discussed next, energization is observed for electrons near the trapped–passing boundary, whereas the cooling is most significant for deeply trapped electrons.

Physically, the form of $f$ in (2.5) is the result of electron heating proportional to $\mu \left \langle \partial B / \partial t\right \rangle _{\rm {orbit}}$, where $\mu$ is the electron magnetic moment. We can use this insight to account for the anisotropic features observed in the theoretical distributions displayed in figure 2(c,d). For the black passing trajectory in figure 1 we note that the orbit-averaged value $\left \langle \partial B / \partial t\right \rangle _{\rm {orbit}}$ is negligible because this trajectory evenly observes positive and negative values of $\partial B / \partial t$, yielding negligible energy changes during the pump perturbations. Meanwhile, the green orbit is deeply trapped and as the perturbation grows this orbit will observe a negative value of $\left \langle \partial B / \partial t\right \rangle _{\rm {orbit}}$, resulting in a net cooling of deeply trapped electrons. The limiting case between trapped and passing electrons (not shown) is also interesting. Such marginally trapped electrons will with relatively low values of $v_{\parallel }/v$ reach the location where $B$ is the strongest and $\partial B / \partial t$ is positive. As a result, the marginally trapped electrons spend a relatively long time in these regions resulting in positive values of $\left \langle \partial B / \partial t\right \rangle _{\rm {orbit}}$ and are therefore subject to energization as the standing-wave perturbation grows.

The distributions of the type predicted by (2.5) and displayed in figure 2(c,d) are consistent with electron distributions observed by MMS during bow-shock crossings (Lichko & Egedal Reference Lichko and Egedal2020). The described energization of the marginally trapped electrons results in enhanced values of their phase-space density, and this is consistent with the distribution in figure 2(d) computed for the magnetic geometry in (2.2) with ${\tilde {B}}_m=B_m/B_0=1.5$. For a position $B/B_0=1.25$ near the peak magnetic field, the particles with $v_{\parallel }/v\simeq 0$ just inside the trapped–passing boundary are all marginally trapped with enhanced values of $f$. In figure 2(c) the distribution is evaluated at the minimum in $B$ and here the region of $v_{\parallel }/v\simeq 0$ corresponds to deeply trapped (cooled) electrons with reduced values of $f$. Still, in the vicinity of the inferred trapped–passing boundaries, the distributions at reduced $B$ include the marginally trapped orbits which in each case account for the enhanced values of $f$ observed at pitch angles near the trapped–passing boundaries.

3. Simple picture of magnetic pumping with trapped electrons

The description above for the evolution of $f$ does not include any net energizations because $f$ in (2.5) will simply return to its original form after a full fluctuation period is completed. Energization by pumping therefore requires a pitch angle mixing process that changes the magnetic moment of the electrons in time. Considering figure 2(b), for ${\tilde {B}}_m=1.5$ an electron with $\varLambda \simeq 0.7$ will be energized to a level of ${\mathcal {E}}\simeq 1.1{\mathcal {E}}_0$ as the magnetic field increases. For the most efficient scenario, such an electron may then pitch-angle scatter to a value of $\varLambda \simeq \varLambda _m=2$, and as the perturbation reverts back to the uniform configuration with $B=B_0$, this electron will acquire an additional boost to reach ${\mathcal {E}} \simeq (1/0.5) \times 1.1{\mathcal {E}}_0= 2.2{\mathcal {E}}_0$. Meanwhile, another electron may initially start out with $\varLambda =\varLambda _m$, and then scatter to $\varLambda \simeq 0.7$ as the magnetic field decreases, thus ending up with a lower energy of ${\mathcal {E}}\simeq {\mathcal {E}}_0/2.2$. As is characteristic of a second-order Fermi process, an important observation is here that the energy gain/loss is proportional to the particle's initial energy and that the step size (up or down) in energy is random.

The described diffusion process is characterized by a diffusion coefficient $D=(\Delta v)^2/\tau$, where $\Delta v$ is the characteristic diffusive step in particle speed, while $\tau$ is the characteristic time associated with this step. Given $\Delta v\propto v$, we then have $D= \omega v^2 {\mathcal {G}}$, where ${\mathcal {G}}({\tilde {B}}_m,\nu /\omega )$ is a dimensionless function of the magnetic perturbation size and the scattering rate $\nu$ normalized by the magnetic pump angular frequency $\omega$. With Fick's law, the diffusion will cause a flux of particles in velocity space given by $\varGamma _v=-D \boldsymbol {\nabla } f(v)$. From particle conservation it then follows that $\partial f / \partial t +\boldsymbol {\nabla }\boldsymbol {\cdot } \varGamma _v=0$, which in spherical velocity space becomes

(3.1)\begin{equation} \frac{\displaystyle \partial f}{\displaystyle \partial t}= \frac{1}{v^2}\frac{\displaystyle \partial }{\displaystyle \partial v} v^2 D \frac{\displaystyle \partial f}{\displaystyle \partial v}- \frac{1}{3n} \frac{\displaystyle \partial n}{\displaystyle \partial t} v\frac{\displaystyle \partial f}{\displaystyle \partial v} . \end{equation}

In addition to diffusion by pumping, we here also include the term proportional to $\partial n / \partial t$ corresponding to changes in $f$ due to compression of the plasma, a well-known effect that is also present in the Parker equation (Parker Reference Parker1965; Montag et al. Reference Montag, Egedal, Lichko and Wetherton2017). In the following sections we use the drift-kinetic framework to rigorously derive (3.1) and obtain an analytical expression for the function ${\mathcal {G}}$ describing the efficiency of the Fermi process.

Physically, the energy for particle energization is provided by mechanical work on the flux tube as captured by the term $p_{\perp }\boldsymbol {\nabla }_{\perp }\boldsymbol {\cdot } {\boldsymbol {v}}_{\perp }$ in the standard two-fluid energy equation. In particular, for the regions where the magnetic field is decreasing the flux tube is expanding and the mechanical work during this phase will take energy out of the electrons. Meanwhile, as the field is expanding distributions of the type in figure 2(c) will form locally with reduced values of $p_{\perp }$. Pitch angle diffusion will naturally cause some level of isotropization, leading to enhanced values of $p_{\perp }$. Therefore more mechanical work is required to recompress the flux tube. As a consequence, when averaged over a full cycle a positive amount of mechanical work is delivered from the fluctuation to the electrons.

Based on the physical reasoning that leads to (3.1) and the similarity to magnetic pumping in an infinite, uniform flux-tube geometry (Lichko et al. Reference Lichko, Egedal, Daughton and Kasper2017), it is possible to obtain quick approximations for the function ${\mathcal {G}}$ characterizing the rate of speed diffusion. However, the aim of the present paper is to provide a mathematically rigorous derivation starting from the drift-kinetic evolution equation for $f$, and then arrive at an equation similar in form to (3.1). This derivation will then in separate papers be generalized and combined with related heating mechanisms (Egedal et al. Reference Egedal, Schroeder and Lichko2021) to account for magnetic pumping in more general magnetic perturbations including compressional waves of multiple amplitudes and phases.

4. Kinetic equation including pitch angle mixing

Following the analysis in Montag et al. (Reference Montag, Egedal, Lichko and Wetherton2017) and supplemented with additional discussion in Appendix A, for the assumed periodic geometry the drift-kinetic equation can be written as

(4.1)\begin{equation} \frac{\displaystyle \textrm{d}f}{\displaystyle \textrm{d}t}= \left.\frac{\displaystyle \partial f}{\displaystyle \partial t}\right|_{{{\varLambda}},{\mathcal{E}}} + \frac{\displaystyle \textrm{d}{{\varLambda}}}{\displaystyle \textrm{d}t}\frac{\displaystyle \partial f}{\displaystyle \partial {{\varLambda}}} + \frac{\displaystyle \textrm{d}{\mathcal{E}}}{\displaystyle \textrm{d}t} \frac{\displaystyle \partial f}{\displaystyle \partial {\mathcal{E}}} = \nu\left\langle{\mathcal{L}}\right\rangle_x f. \end{equation}

For $\nu =0$, (4.1) provides a representation of the collisionless fast transit-time limit, whereas for finite scattering frequencies $\nu \neq 0$ pitch angle mixing is included through the orbit average operator:

(4.2)\begin{equation} \left\langle{\mathcal{L}}\right\rangle_x= \frac{1}{\tau_b} \oint \frac{{\mathcal{L}}}{v_{{\parallel}}} \textrm{d}l. \end{equation}

Here the local scattering can be characterized by, for example, the standard Lorentz operator ${\mathcal {L}}=\partial / \partial \xi (1-\xi ^2)\partial / \partial \xi$, with $\xi =v_{\parallel }/v$. Furthermore, we have also introduced the orbit bounce time

(4.3)\begin{equation} \tau_b= \oint \frac{ \textrm{d}l}{v_{{\parallel}}} = 4 \int_0^{L_b} \frac{1}{v_{{\parallel}}}{\textrm{d} x}' , \end{equation}

which together with its normalized form,

(4.4)\begin{equation} {\tilde{\tau}_b}=\frac{v\tau_b}{2L}= \frac{2v}{L} \int_0^{L_b} \frac{1}{v_{{\parallel}}}{\textrm{d} x}' , \end{equation}

will become important in the following calculations.

We now introduce the dimensionless drive function $g(\varLambda,t)$ through the definition

(4.5)\begin{equation} g \equiv \frac{\displaystyle 1}{\displaystyle {{\varLambda}}}\frac{\displaystyle \textrm{d}{{\varLambda}}}{\displaystyle \textrm{d}t}={-}\frac{\displaystyle 1}{\displaystyle {\mathcal{E}}} \frac{\displaystyle \textrm{d}{\mathcal{E}}}{\displaystyle \textrm{d}t} . \end{equation}

Here the expression in terms of ${\mathcal {E}}$ follows from the form of ${{\varLambda }}$ in (2.4a,b), and we may then write the kinetic equation in (4.1) in the form

(4.6)\begin{equation} \frac{\displaystyle \partial f}{\displaystyle \partial t}= g \,{\mathcal{E}} \frac{\displaystyle \partial f}{\displaystyle \partial {\mathcal{E}}} -g {{\varLambda}} \frac{\displaystyle \partial f}{\displaystyle \partial {{\varLambda}}} + \nu{\left\langle{\mathcal{L}}\right\rangle_x}f, \end{equation}

which is equivalent to

(4.7)\begin{equation} \frac{\displaystyle \partial f}{\displaystyle \partial t}= g \,\frac{\displaystyle v}{\displaystyle 2} \frac{\displaystyle \partial f}{\displaystyle \partial v} -g{{\varLambda}} \frac{\displaystyle \partial f}{\displaystyle \partial {{\varLambda}}} + \nu{\left\langle{\mathcal{L}}\right\rangle_x}f . \end{equation}

In order for (4.7) to be useful we need to obtain an expression for $g$ related to the time evolution of the magnetic perturbation. For this we explore that the ratio $J^2/\mu$ is a constant of motion which has the special property that it is independent of ${\mathcal {E}}$ (Egedal Reference Egedal2002), and using $j$ in (2.6) it follows that $0=\textrm {d}/\textrm {d}t (J^2/\mu ) = \textrm {d}/\textrm {d}t(j^2B_0/{{\varLambda }})$, such that

(4.8)\begin{equation} \frac{\displaystyle 2}{\displaystyle j}\frac{\displaystyle \textrm{d}j}{\displaystyle \textrm{d}t}- \frac{\displaystyle 1}{\displaystyle {{\varLambda}}}\frac{\displaystyle \textrm{d}{{\varLambda}}}{\displaystyle \textrm{d}t}=0. \end{equation}

To evaluate the total differential of $j$ we further notice that $j=j({{\varLambda }},t)$, while ${{\varLambda }}={{\varLambda }}(t)$ such that

(4.9)\begin{equation} \frac{\displaystyle \textrm{d}j}{\displaystyle \textrm{d}t} = \frac{\displaystyle \partial j}{\displaystyle \partial t} + \frac{\displaystyle \textrm{d}{{\varLambda}}}{\displaystyle \textrm{d}t} \frac{\displaystyle \partial j}{\displaystyle \partial {{\varLambda}}} . \end{equation}

Combining (4.8) and (4.9) we find

(4.10)\begin{equation} g({{\varLambda}},t)= \frac{\displaystyle 1}{\displaystyle {{\varLambda}}}\frac{\displaystyle \textrm{d}{{\varLambda}}}{\displaystyle \textrm{d}t} =\frac{\displaystyle 2\partial j / \partial t}{\displaystyle j- 2{{\varLambda}}\, \partial j / \partial {{\varLambda}}} = \frac{\displaystyle 2}{\displaystyle {\tilde{\tau}_b}}\frac{\displaystyle \partial j}{\displaystyle \partial t}. \end{equation}

Here we used

(4.11)\begin{equation} {\tilde{\tau}_b}=j-2{{\varLambda}}\frac{\displaystyle \partial j}{\displaystyle \partial {{\varLambda}}}, \end{equation}

which is readily shown using

(4.12a,b)\begin{equation} \frac{\displaystyle \tau_b}{\displaystyle m}=\left.\frac{\displaystyle \partial J}{\displaystyle \partial {\mathcal{E}}}\right|_{\mu},\quad\left.\frac{\displaystyle \partial }{\displaystyle \partial {\mathcal{E}}}\right|_{\mu} = \left.\frac{\displaystyle \partial }{\displaystyle \partial {\mathcal{E}}}\right|_{{{\varLambda}}} + 2{{\varLambda}}\frac{\displaystyle \partial }{\displaystyle \partial {{\varLambda}}}. \end{equation}

Based on the definitions of $j$ and ${\tilde {\tau }_b}$ in (2.6) and (4.4), and with the drive term $g$ given in (4.10), the right-hand side of (4.7) is now readily evaluated numerically for any prescribed perturbation of the magnetic field, $B(x,t)$.

5. Particle conservation and a Krook scattering model

In our model, no electrons are permitted to ‘leak’ out of the considered flux tube, and particle conservation becomes a concept important to the derivations of the present framework. Without loss of generality, we assume that the magnetic field lines lie in the $(x,z)$ plane, conveniently characterized by the flux function $\varPsi$. It then also follows that the flux tube has constant width $\Delta y_0$ in the $y$ direction. The total number of particles $N$ within our flux tube of length $L$ must be conserved, and using results of Montag et al. (Reference Montag, Egedal, Lichko and Wetherton2017) it follows that

(5.1)\begin{equation} N = \int f\, \textrm{d}^3v\, \textrm{d}^3x = \frac{2 {\rm \pi}\Delta y_0}{ m^2} \int \textrm{d} {\mathcal{E}} \int \textrm{d} \mu \int \textrm{d} \varPsi \tau_b f . \end{equation}

The width of the tube in the $z$ direction varies as $1/B$, but for the point $x = L/4$ the magnetic field is constant in time $B = B_0$, and we have $\Delta \varPsi = B_0 \Delta z_0$, where $\Delta z_0$ is the width of the flux tube at this location (and also the width of the initially uniform configuration). Furthermore, because $\tilde {\tau }_b = v \tau _b/(2 L)$ and $\textrm {d} {\mathcal {E}} \,\textrm {d} \mu = m^2 v^3 \,\textrm {d}v \,\textrm {d} \varLambda / B_0$ we get

(5.2)\begin{equation} N = 4 {\rm \pi}\Delta y_0 \Delta z_0 L \int v^2 \,\textrm{d}v \int \textrm{d} \varLambda \tilde{\tau}_b f . \end{equation}

Guided by (5.2) we introduce the averaging operator

(5.3)\begin{equation} \left\langle (\dots) \right\rangle_{\varLambda} =\int \textrm{d}{{\varLambda}} \,{\tilde{\tau}_b} (\dots), \end{equation}

where the range of integration over the full domain of $\varLambda$ changes in time, $0\leq {{\varLambda }}\leq {{\varLambda }}_m$ with ${{\varLambda }}_m=1/{B_{\rm {min}} }(t)$. We define

(5.4)\begin{equation} F(v)= \left\langle f \right\rangle_{\varLambda}, \end{equation}

such that

(5.5)\begin{equation} N = \Delta y_0 \Delta z_0 L \int 4{\rm \pi} F(v) v^2 \,\textrm{d}v . \end{equation}

From (5.5) it is clear that $F(v)v^2\,\textrm {d}v$ is proportional to the number of particles of the entire flux tube within a differential velocity interval $\textrm {d}v$. It also follows that the volume $V$ of the flux tube is given by

(5.6)\begin{equation} V= \Delta y_0 \Delta z_0L \left\langle 1 \right\rangle_{\varLambda} = \Delta y_0 \Delta z_0 L \int \textrm{d}{{\varLambda}} \,{\tilde{\tau}_b}, \end{equation}

and $\left \langle 1 \right \rangle _{\varLambda }$ therefore represents the relative change of the volume during the evolution of the magnetic perturbation.

If we assume a static magnetic configuration while integrating (4.1) for a long time compared to $1/\nu$, any initial distribution $f({\mathcal {E}},\varLambda )$ will pitch-angle mix into a form $f_{\rm {mixed}} ({\mathcal {E}})$ independent of $\varLambda$. Particle conservation at each energy then requires that $\left \langle f({\mathcal {E}},\varLambda ) \right \rangle _{\varLambda }= \left \langle f_{\rm {mixed}} ({\mathcal {E}}) \right \rangle _{\varLambda }= f_{\rm {mixed}} ({\mathcal {E}}) \left \langle 1 \right \rangle _{\varLambda }$. This observation motivates the construction of a Krook scattering operator:

(5.7)\begin{equation} {\mathcal{L}}_Kf ={-}C_K \left(f - \frac{\displaystyle \left\langle f \right\rangle_{\varLambda}}{\displaystyle \left\langle 1 \right\rangle_{\varLambda}}\right), \end{equation}

which compared with the orbit-averaged Lorentz operator, $\left \langle {\mathcal {L}}\right \rangle _x$, is much better suited for analytical calculations. The factor $C_K$ was introduced by Lichko & Egedal (Reference Lichko and Egedal2020) to calibrate ${\mathcal {L}}_K$ to have scattering efficiency similar to $\left \langle {\mathcal {L}}\right \rangle _x$, as a function of the velocity space scale length of the anisotropic features in $f({\mathcal {E}},\varLambda )$. For example, for situations where the anisotropy of the local distribution is dominated by the $P_l(v_{\parallel }/v)$ Legendre polynomial of order $l$, the appropriate calibration factor becomes $C_K= l(l+1)$.

6. Kinetic model of magnetic pumping

In the present section we derive an evolution equation for a slowly evolving background distribution energized by the pumping process. Consistent with the Krook operator in (5.7) the isotropic part of the distribution is given by

(6.1)\begin{equation} f_0= \frac{\displaystyle \left\langle f \right\rangle_{\varLambda}}{\displaystyle \left\langle 1 \right\rangle_{\varLambda}} = \frac{\displaystyle 1}{\displaystyle \left\langle 1 \right\rangle_{\varLambda}} F(v) . \end{equation}

During each pump cycle the volume of the flux tube may change as characterized by $\left \langle 1 \right \rangle _{\varLambda }$, and $f_0$ therefore does not necessarily represent a slowly evolving distribution. For this reason it turns out that the pumping process is better characterized in terms of the $\varLambda$-averaged distribution, $F(v)=\left \langle f \right \rangle _{\varLambda }=f_0\left \langle 1 \right \rangle _{\varLambda }$, introduced above in (5.4).

6.1. Limit of strong scattering, $\nu \gg \omega$

Our aim is first to derive an equation for $\partial F / \partial t$, which describes the evolution of $F(v,t)$ due to magnetic pumping in the limit of strong scattering, $\nu \gg \omega$, and small magnetic perturbation amplitudes. By direct differentiation of (5.4) we obtain

(6.2)\begin{equation} \frac{\displaystyle \partial F}{\displaystyle \partial t} =\left\langle \frac{\displaystyle \partial f}{\displaystyle \partial t} \right\rangle_{\varLambda} + \int \textrm{d}{{\varLambda}} \, \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} \, f + \left. {\tilde{\tau}_b} f\right|_{{{\varLambda}}_m}\dot{{{\varLambda}}}_m , \end{equation}

where $\dot {{{\varLambda }}}_m =\partial {{\varLambda }}_m / \partial t$. To evaluate the various terms in (6.2) we introduce an expansion for $f$ for which the first term is the isotropic component:

(6.3)\begin{equation} f=f_0(t,v) + f_1(t,v,{{\varLambda}})+ {\cdots}. \end{equation}

We here assume a small-amplitude ordering $\Delta B/B\ll 1$, such that $|g| \simeq |f_1| \ll |f_0|$, and the kinetic equation in (4.7) with the Krook approximation in (5.7) then yields an approximate equation for $f_1$:

(6.4)\begin{equation} \frac{\displaystyle \partial f_1}{\displaystyle \partial t} = h\frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial f_0}{\displaystyle \partial v} + \nu {\mathcal{L}}_K f_1, \end{equation}

where

(6.5)\begin{equation} h\equiv g-(\left\langle g \right\rangle_{\varLambda}/\left\langle 1 \right\rangle_{\varLambda}) . \end{equation}

Furthermore, by introducing the ordering $\nu /\omega \gg 1$, we can ignore $\partial f_1 / \partial t$ on the left-hand side of (6.4), such that

(6.6)\begin{equation} f_1= \frac{\displaystyle h}{\displaystyle \nu C_K}\frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial f_0}{\displaystyle \partial v}. \end{equation}

Considering again (6.2) we may now evaluate the two last terms on the right-hand side. First we observe that

(6.7)\begin{equation} \int \textrm{d}{{\varLambda}} f \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} = f_0\int \textrm{d}{{\varLambda}} \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} +\frac{\displaystyle 1}{\displaystyle \nu C_K} \int h \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} \textrm{d}{{\varLambda}} \,\frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial }{\displaystyle \partial v}f_0, \end{equation}

while

(6.8)\begin{equation} \left.{\tilde{\tau}_b} f\right|_{{{\varLambda}}_m}\dot{{{\varLambda}}}_m = \left.{\tilde{\tau}_b} \dot{{{\varLambda}}}_m \right|_{{{\varLambda}}_m}f_0 + \left. {\tilde{\tau}_b} h \dot{{{\varLambda}}}_m \right|_{{{\varLambda}}_m} \frac{\displaystyle 1}{\displaystyle \nu C_K}\frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial f_0}{\displaystyle \partial v}, \end{equation}

such that (6.2) becomes

(6.9)\begin{equation} \frac{\displaystyle \partial F}{\displaystyle \partial t} = \left\langle \frac{\displaystyle \partial f}{\displaystyle \partial t} \right\rangle_{\varLambda} + f_0\int \textrm{d}{{\varLambda}} \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} + T \,\frac{\displaystyle 1}{\displaystyle \nu C_K}\frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial f_0}{\displaystyle \partial v} , \end{equation}

where

(6.10)\begin{equation} T= \int h \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} \textrm{d}{{\varLambda}} + \left.{\tilde{\tau}_b} h\right|_{{{\varLambda}}_m} \dot{{{\varLambda}}}_m . \end{equation}

To evaluate $\left \langle \partial f / \partial t \right \rangle _{\varLambda }$ we insert $f=f_0+f_1$ back into (4.7):

(6.11)\begin{equation} \frac{\displaystyle \partial f}{\displaystyle \partial t} = g \,\frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial f_0+f_1}{\displaystyle \partial v} - g{{\varLambda}} \frac{\displaystyle \partial f_1}{\displaystyle \partial {{\varLambda}}}- \nu C_K f_1. \end{equation}

Then, using $\left \langle f_1 \right \rangle _{\varLambda }=0$ and $\left \langle kf_0 \right \rangle _{\varLambda }=f_0\left \langle k \right \rangle _{\varLambda }$ we find

(6.12)\begin{align} \left\langle \frac{\displaystyle \partial f}{\displaystyle \partial t} \right\rangle_{\varLambda} & = \left\langle g \right\rangle_{\varLambda} \frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial }{\displaystyle \partial v}f_0+ \frac{\displaystyle 1}{\displaystyle \nu C_K}\left\langle gh \right\rangle_{\varLambda} \frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial }{\displaystyle \partial v}\frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial }{\displaystyle \partial v}f_0 -\frac{\displaystyle 1}{\displaystyle \nu C_K}\left\langle g {{\varLambda}} \frac{\displaystyle \partial h}{\displaystyle \partial {{\varLambda}}} \right\rangle_{\varLambda}\frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial }{\displaystyle \partial v}f_0\nonumber\\ & = \left\langle g \right\rangle_{\varLambda}\frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial }{\displaystyle \partial v}f_0+ \frac{\displaystyle 1}{\displaystyle \nu C_K}\left\langle gh \right\rangle_{\varLambda}\frac{\displaystyle 1}{\displaystyle 4v^2}\frac{\displaystyle \partial }{\displaystyle \partial v}v^4\frac{\displaystyle \partial }{\displaystyle \partial v}f_0 -\frac{\displaystyle 1}{\displaystyle \nu C_K} T \frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial }{\displaystyle \partial v}f_0. \end{align}

Here the last form of (6.12) assumes that

(6.13)\begin{equation} \frac{\displaystyle 3}{\displaystyle 2}\left\langle gh \right\rangle_{\varLambda} ={-}\left\langle g {{\varLambda}} \frac{\displaystyle \partial h}{\displaystyle \partial {{\varLambda}}} \right\rangle_{\varLambda} + T, \end{equation}

which is proved below. In addition we also verify a similar expression needed in the next few steps:

(6.14)\begin{equation} \frac{\displaystyle 3}{\displaystyle 2}\left\langle g \right\rangle_{\varLambda} = \int \textrm{d}{{\varLambda}} \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} + \left. {\tilde{\tau}_b} \right|_{{{\varLambda}}_m}\dot{{{\varLambda}}}_m=\frac{\displaystyle \partial }{\displaystyle \partial t} \left\langle 1 \right\rangle_{\varLambda} . \end{equation}

To prove (6.13) and (6.14) we recall that ${\tilde {\tau }_b} g=2\partial j / \partial t$, $j= {\tilde {\tau }_b} + 2 {{\varLambda }} \partial j / \partial {{\varLambda }}$ and $j|_{{{\varLambda }}_m}=0$. As an intermediate result, for any $G=G(t,{{\varLambda }})$ it then follows that

(6.15)

We now consider $\left \langle gh \right \rangle _{{{\varLambda }}} = 2 \int \textrm {d}{{\varLambda }} h\partial j / \partial t$, and use (

6.15

) with $G=h$ to simplify the expressions:

(6.16)\begin{align} \left\langle gh \right\rangle_{{{\varLambda}}}& = 2 \int \textrm{d}{{\varLambda}} h\frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} + 4 \int d{{\varLambda}} \,{{\varLambda}} h\frac{\displaystyle \partial }{\displaystyle \partial {{\varLambda}}} \frac{\displaystyle \partial j}{\displaystyle \partial t} \nonumber\\ & =2 \int \textrm{d}{{\varLambda}} h \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} + 4 \int \textrm{d}{{\varLambda}} \frac{\displaystyle \partial }{\displaystyle \partial {{\varLambda}}} h{{\varLambda}} \,\frac{\displaystyle \partial j}{\displaystyle \partial t} - 4 \int \textrm{d}{{\varLambda}} h \frac{\displaystyle \partial j}{\displaystyle \partial t} - 4 \int \textrm{d}{{\varLambda}} \frac{\displaystyle \partial j}{\displaystyle \partial t} \varLambda \frac{\displaystyle \partial }{\displaystyle \partial \varLambda} h \nonumber\\ & =2 \int \textrm{d}{{\varLambda}} h \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} + 2 h \dot{{{\varLambda}}}_m \left. {\tilde{\tau}_b} \right|_{{{\varLambda}}_m} - 2\left\langle gh \right\rangle_{{{\varLambda}}} - 2 \left\langle g\varLambda \frac{\displaystyle \partial h}{\displaystyle \partial \varLambda} \right\rangle_{{{\varLambda}}} . \end{align}

From here (6.13) is readily obtained.

Similarly, this time with $G=1$ in (6.15), we get

(6.17)\begin{align} \left\langle g \right\rangle_{{{\varLambda}}}& =2 \int \textrm{d}{{\varLambda}} \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} + 4 \int \textrm{d}{{\varLambda}} \,{{\varLambda}} \frac{\displaystyle \partial }{\displaystyle \partial {{\varLambda}}} \frac{\displaystyle \partial j}{\displaystyle \partial t} \nonumber\\ & =2 \int \textrm{d}{{\varLambda}} \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} + 4 \int \textrm{d}{{\varLambda}} \frac{\displaystyle \partial }{\displaystyle \partial {{\varLambda}}} {{\varLambda}} \frac{\displaystyle \partial j}{\displaystyle \partial t} - 4 \int \textrm{d}{{\varLambda}} \frac{\displaystyle \partial j}{\displaystyle \partial t} \nonumber\\ & =2 \int \textrm{d}{{\varLambda}} \frac{\displaystyle \partial {\tilde{\tau}_b}}{\displaystyle \partial t} + 2 \dot{{{\varLambda}}}_m \left. {\tilde{\tau}_b} \right|_{{{\varLambda}}_m} - 2 \left\langle g \right\rangle_{{{\varLambda}}}, \end{align}

and from here (6.14) follows directly.

Inserting (6.12) into (6.9) and using (6.14) together with $f_0=F/\left \langle 1 \right \rangle _{\varLambda }$ yields the desired evolution equation for $F(v,t)$:

(6.18)\begin{equation} \frac{\displaystyle \partial F}{\displaystyle \partial t} = \frac{\displaystyle 1}{\displaystyle \nu C_K}\frac{\displaystyle \left\langle gh \right\rangle_{\varLambda}}{\displaystyle \left\langle 1 \right\rangle_{\varLambda}}\,\frac{\displaystyle 1}{\displaystyle 4v^2}\frac{\displaystyle \partial }{\displaystyle \partial v}v^4\frac{\displaystyle \partial }{\displaystyle \partial v}F + \, \frac{\displaystyle \left\langle g \right\rangle_{\varLambda}}{\displaystyle 2\left\langle 1 \right\rangle_{\varLambda}}\,\frac{\displaystyle 1}{\displaystyle v^2}\frac{\displaystyle \partial }{\displaystyle \partial v}v^3F. \end{equation}

Again, the result in (6.18) is obtained in the limit $\nu \gg \omega$ and describes the evolution of a near-isotropic electron plasma. Because $F|_{v=\infty }=0$, through integration by parts it is readily seen that $\int \partial F / \partial t\, v^2 \,\textrm {d}v=0$ and (6.18) is therefore consistent with particle conservation. The anisotropic features in $f$ are given by $f_1$ proportional to $h/(\nu C_k)$ with $h$ given in (6.5). The term proportional to $\left \langle gh \right \rangle _{\varLambda }/(\nu C_k)$ then accounts for the bi-linear interaction between $h/(\nu C_k)$ and the drive $g$, yielding velocity diffusion akin to Fermi acceleration.

To understand the term in (6.18) proportional to $\left \langle g \right \rangle _{\varLambda }/(2\left \langle 1 \right \rangle _{\varLambda })$ we recall (6.14), and note again that $\left \langle 1 \right \rangle _{\varLambda }$ is proportional to the volume of the flux tube, such that $3/2 \, \left \langle g \right \rangle _{\varLambda }/\left \langle 1 \right \rangle _{\varLambda } = \partial \log \left \langle 1 \right \rangle _{\varLambda } / \partial t$ is the rate of relative compression of the flux tube. In terms of the average plasma density $n\equiv N/(\Delta y_0\Delta z_0 L \left \langle 1 \right \rangle _{\varLambda })$, it follows that $\left \langle g \right \rangle _{\varLambda }/(2\left \langle 1 \right \rangle _{\varLambda }) = - (1/3) \partial \log (n) / \partial t$, and it is then clear that this last term of (6.18) accounts for compressional heating of an isotropic plasma, as is familiar from the similar term in the Parker equation (Parker Reference Parker1965) and also included in (3.1) above.

6.2. General case, including arbitrary scattering levels

Above, we solved (6.4) in the limit $\nu /\omega \gg 1$ to get an approximate solution for $f_1$. Still for $\Delta B/B_0\ll 1$, to address the more general case with no restrictions on $\nu /\omega$, we apply Fourier expansions of $g$, $h$ and $f_1$:

(6.19)\begin{equation} \begin{gathered}g(\eta,t)=g^1(\eta) \textrm{e}^{\textrm{i}\omega t} + g^2(\eta) \textrm{e}^{2\textrm{i}\omega t} + {\cdots },\\ h(\eta,t)=h^1(\eta) \textrm{e}^{\textrm{i}\omega t} + h^2(\eta) \textrm{e}^{2\textrm{i}\omega t} + {\cdots},\\ f_1(\eta,t)=f_1^1(\eta) \textrm{e}^{\textrm{i}\omega t} + f_1^2(\eta) \textrm{e}^{2\textrm{i}\omega t} + {\cdots }. \end{gathered}\end{equation}

For this to be meaningful we have introduced the new variable $\eta$, for which the range is constant in time $0\leq \eta \leq 1$. Furthermore, because

(6.20)\begin{equation} \int h {\tilde{\tau}_b} \,\textrm{d}{{\varLambda}} =0, \end{equation}

this constraint must be fulfilled by each of the terms in the sum

(6.21)\begin{equation} \int h^n {\tilde{\tau}_b} \frac{\displaystyle \partial {{\varLambda}}}{\displaystyle \partial \eta} \textrm{d}\eta =0. \end{equation}

For the general case where ${\tilde {\tau }_b}$ depends on $(\varLambda,t)$, the constraint implies that ${\tilde {\tau }_b} \partial {{\varLambda }} / \partial \eta =K(t)$, where $K(t)$ is an arbitrary function independent of $\varLambda$. Imposing further that $0\leq \eta \leq 1$ and $\eta =1$ for $\varLambda =0$, there is then only one definition of $\eta$ available:

(6.22a,b)\begin{equation} \frac{\displaystyle \partial \eta}{\displaystyle \partial {{\varLambda}}} ={-}\frac{\displaystyle {\tilde{\tau}_b}}{\displaystyle \left\langle 1 \right\rangle_{\varLambda}} ,\quad \eta=1-\frac{\displaystyle 1}{\displaystyle \left\langle 1 \right\rangle_{\varLambda}}\int_0^{{{\varLambda}}} {\tilde{\tau}_b}\, \textrm{d}{{\varLambda}}'. \end{equation}

To emphasize some properties of $\eta$, the averaging operator in (5.3) may be expressed as $\left \langle (\dots ) \right \rangle _{\varLambda }= \left \langle 1 \right \rangle _{\varLambda }\int (\dots ) \,\textrm {d}\eta$, and for the uniform flux tube we observe that $\eta =v_{\parallel }/v$.

Returning to (6.4) we recall that this equation for $f_1$ was derived with $(v,{{\varLambda }})$ as the independent velocity space variables. In principle, with $f_1=f_1(v,\eta,t)$ a new term will then occur as we note that

(6.23)\begin{equation} \left.\frac{\displaystyle \partial f_1}{\displaystyle \partial t}\right|_{{{\varLambda}},{\mathcal{E}}} = \left.\frac{\displaystyle \partial f_1}{\displaystyle \partial t}\right|_{\eta,{\mathcal{E}}} + \left. \frac{\displaystyle \partial \eta}{\displaystyle \partial t} \right|_{{{\varLambda}}} \frac{\displaystyle \partial f_1}{\displaystyle \partial \eta}. \end{equation}

However, because $\partial \eta / \partial t |_{{{\varLambda }}}=0$ for $\Delta B=0$, the term $\partial \eta / \partial t |_{{{\varLambda }}}\, \partial f_1 / \partial \eta$ is of second order in $\Delta B/B_0$ and can be ignored. This permits us to write $\partial f_1 / \partial t=\textrm {i}\omega f_1$, and similar to (6.4) for the separate frequencies the equation for $f_1^n$ then reads

(6.24)\begin{equation} \textrm{i}n\omega f_1^n =h^n\frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial f_0}{\displaystyle \partial v} - C_K \nu f_1^n , \end{equation}

such that $f_1=\sum _n f_1^n$, with

(6.25a,b)\begin{equation} f_1^n= K_n \frac{\displaystyle v}{\displaystyle 2}\frac{\displaystyle \partial f_0}{\displaystyle \partial v}, \quad K_n= \frac{\displaystyle h^n(-\textrm{i}n\omega+C_K\nu)}{\displaystyle n^2\omega^2+C_K^2\nu^2}. \end{equation}

Using (6.25a,b) in place of (6.6), we may now follow the same steps that led from (6.6) to (6.18). Analogous to before when pumping was driven by $\left \langle \left \langle gf_1 \right \rangle _{\varLambda }\right \rangle _t$, we find that parts of $f_1^n$ are in phase with $g^n$ causing the time average $\left \langle \left \langle g^nf_1^n \right \rangle _{\varLambda }\right \rangle _t$ to become finite. An equation for the slowly varying $F(v,t)$ is then obtained as

(6.26)\begin{equation} \left.\begin{gathered} \frac{\displaystyle \partial F}{\displaystyle \partial t}= \frac{\displaystyle 1}{\displaystyle v^2}\frac{\displaystyle \partial }{\displaystyle \partial v} \left(v^2D\frac{\displaystyle \partial }{\displaystyle \partial v}F\right),\quad D=\omega v^2{\mathcal{G}},\\ {\mathcal{G}}\left(\frac{\nu C_K}{\omega}\right)\equiv \frac{\displaystyle 1}{\displaystyle 4}\sum_n\frac{\displaystyle \nu C_K/\omega}{\displaystyle n^2+(\nu C_k/\omega)^2} \left\langle\left\langle g^nh^n \right\rangle_{\varLambda}\right\rangle_t ,\\ \left\langle\left\langle g^nh^n \right\rangle_{\varLambda}\right\rangle_t\equiv\left\langle\frac{\displaystyle \left\langle \textrm{Re}(g^n\textrm{e}^{\textrm{i}n\omega t})\textrm{Re}(h^n\textrm{e}^{\textrm{i}n\omega t}) \right\rangle_{\varLambda}}{\displaystyle \omega^2\left\langle 1 \right\rangle_{\varLambda}}\right\rangle_t. \end{gathered}\right\} \end{equation}

The efficiency of magnetic pumping with trapped electrons is characterized by (6.26) and is now readily evaluated numerically for any periodic magnetic perturbation, oscillating as a standing wave as a function of time. We may also generalize (6.26) to account for slow flux-tube compression simply by including the term proportional to $\left \langle g \right \rangle _{\varLambda }$ in (6.18).

6.3. Evaluation of the pumping efficiency

The model in (6.26) is verified by evaluating its predictions as a function of $\nu C_K/\omega$ for a number of discrete perturbation magnitudes, parameterized by $\Delta B/B_0$ in (2.1). These predictions are then compared with results obtained by numerically integrating (4.7) directly, and using equation (7) of Lichko & Egedal (Reference Lichko and Egedal2020) to estimate ${\mathcal {G}}$ from the numerical evolution of $f(v,t)$. This integration can be carried out with $\left \langle {\mathcal {L}}\right \rangle _x$ based on the Lorentz form or the Krook form in (5.7), but for direct comparison with the model in (6.26) we here only consider results obtained using the Krook representation.

The results of the numerical integration of (4.7) are represented by the red crosses in figure 3. In comparison, the lines are obtained from (6.26) based on the following steps. For a given value of $\Delta B/B_0$ the spatial and temporal variation of $B$ is given through (2.1) and (2.2). This then permits direct evaluation of $j$ and ${\tilde {B}}$ as a function of $(\varLambda,t)$ using (2.6) and (4.4). Evaluation of $g$ and $h$ in (4.10) and (6.5) then follows, from which their Fourier expansions in (6.19) are determined. All quantities needed for evaluating the required $\varLambda$ and $t$ averages of (6.26) thereby become available. Given that the time dependency of the configuration is set through the $|\sin (\omega t)|$ term in (2.1), only even terms of $n$ contribute in (6.26). The lines in figure 3 represent the result of the first one, two and five finite terms in the expansion of frequencies. The red lines obtained with $n\in \{2,4\}$ are mostly hidden behind the black lines obtained with $n\in \{2,4,6,8,10\}$ demonstrating the rapid conversion of the expansion.

Figure 3. Solid lines show predictions for ${\mathcal {G}}$ calculated with $\Delta B/B_0 \in \{0.3, 0.5, 0.7, 0.9\}$, and including the first one, two and five finite terms of the sum in (6.26). For comparison, the red crosses are obtained by numerical integration of (4.7), using equation (7) of Lichko & Egedal (Reference Lichko and Egedal2020) to estimate ${\mathcal {G}}$.

For all values of $\Delta B/B_0$ and $\nu C_K/\omega$ the predictions of (6.26) are observed to be in good agreement with the numerical integration results. Lichko & Egedal (Reference Lichko and Egedal2020) further shows that the model based on the Krook scattering operator also is in agreement with numerical results obtained with the Lorentz form of $\left \langle {\mathcal {L}}\right \rangle _x$ (defined as the orbit average of the local ${\mathcal {L}}=\partial / \partial \xi (1-\xi ^2)\partial / \partial \xi$). This agreement, however, requires that the value of $C_K$ be adjusted to become large at small values of $\Delta B/B_0$, where the pitch angle scale size of the features in $f_1$ are the smallest and most efficiently scattered by $\left \langle {\mathcal {L}}\right \rangle _x$. Curve fitting yields the following easy-to-evaluate form accurate for the range of $0\leq \Delta B/B \lesssim 0.9$:

(6.27)\begin{equation} {\mathcal{G}}\left(\frac{\nu C_K}{\omega}\right)\simeq \frac{\displaystyle 0.02\,\nu C_K/\omega}{\displaystyle 4+(\nu C_k/\omega)^2} \left[\left(\frac{\Delta B}{B}\right)^{2.6} +3\left(\frac{\Delta B}{B}\right)^{5.6}\right] . \end{equation}

This form based on numerical results with the Krook scattering operator is somewhat different from the similar expression in Lichko & Egedal (Reference Lichko and Egedal2020) for results based on the Lorentz form of the scattering operator. These differences are mainly due to a normalization mistake, as the expression and curves in Lichko & Egedal (Reference Lichko and Egedal2020) were multiplied up by an erroneous factor of ${\rm \pi}$ (still, the model was correctly applied in the analysis of the MMS observations).

Given the model was derived through an expansion $f=f_0+f_1+\cdots$ with $f_1\ll f_0$, the good agreement between the numerical integration results and those of (6.26) is perhaps surprising. Strictly, agreement should only be expected for $\nu C_K/\omega \gg 1$ and/or $\Delta B/B_0\ll 1$, which ensures that $f_1\ll f_0$. For example, for $\nu C_K/\omega \ll 1$ and $\Delta B/B_0=0.5$ distributions similar to those in figure 2 will occur where the anisotropic part described by $f_1$ is comparable in magnitude to $f_0$. For this limit where $\nu C_K/\omega \lesssim 1$ a model can likely be obtained where again $f=f_0+f_1+\cdots$, but with $f_0=f_0(J,\mu )$ similar to (2.5). We expect that such a framework will provide results similar to (6.26), but the analysis is beyond the scope of the present paper.

7. Summary and conclusions

Lichko & Egedal (Reference Lichko and Egedal2020) developed a new magnetic pumping model which successfully accounts for detailed in situ MMS observations of electron energization at the Earth's bow shock. In the present paper we provide a detailed mathematical derivation of the model obtained by exploring the fast transit-time limit of a drift-kinetic plasma description, including the effects of electrons becoming trapped in the magnetic perturbations. While the model is derived with an assumption of small magnetic perturbations, it also accurately accounts for the heating rates at large magnetic perturbation amplitudes. This is evident with a comparison of the result of the model with those obtained by direct numerical integration of the drift-kinetic equation. Despite the highly idealized form of the magnetic geometry considered, the model provides an accurate description of the anisotropic form and energization of electron distributions observed at the Earth's bow shock. The mathematical framework may also serve as a starting point for additional studies, generalizing the model to more complicated and realistic scenarios where the magnetic perturbations are comprised of multiple waves with variable amplitudes and propagation directions.

We emphasize the role of magnetic trapping for magnetic pumping to become effective for superthermal electrons. Without trapping the rapid thermal streaming of electrons along magnetic field lines renders pumping ineffective for electron energization. Meanwhile, even small-amplitude magnetic perturbations can trap a significant fraction of the electrons, and given that the magnetic mirror force is proportional to ${\mathcal {E}}_{\perp }$, this trapping applies to all energies for which the electrons are magnetized. Particularly important for space and astrophysical applications, the present model represents a second-order Fermi process consistent with the generation of electron distributions with energetic power-law tails.