Access the full text.

Sign up today, get DeepDyve free for 14 days.

Astrophysics
, Volume 2019 (1909) – Sep 23, 2019

/lp/arxiv-cornell-university/electron-pairing-in-mirror-modes-surpassing-the-quasilinear-limit-geK1HptcpS

- ISSN
- 1432-0576
- eISSN
- ARCH-3330
- DOI
- 10.5194/angeo-37-971-2019
- Publisher site
- See Article on Publisher Site

1,3 2 Rudolf A. Treumann and Wolfgang Baumjohann International Space Science Institute, Bern, Switzerland Space Research Institute, Austrian Academy of Sciences, Graz, Austria Geophysics Department, Ludwig-Maximilians-University Munich, Germany, Correspondence to: Wolfgang.Baumjohann@oeaw.ac.at Abstract. The mirror mode evolving in collisionless magnetised high-temperature thermally anisotropic plasmas is shown to develop an interesting macro-state. Starting as a classical zero frequency ion ﬂuid instability it saturates quasi-linearly at very low magnetic level, while forming elongated magnetic bubbles which trap the electron component to perform an adiabatic bounce motion along the magnetic ﬁeld. Further evolution of the mirror mode towards a stationary state is determined by the bouncing trapped electrons which interact with the thermal level of ion sound waves, generate attractive wake potentials which give rise to formation of electron pairs in the lowest-energy singlet state of two combined electrons. Pairing takes preferentially place near the bounce-mirror points where the pairs become spatially locked with all their energy in the gyration. The resulting large anisotropy of pairs enters the mirror growth rate in the quasi-linearly stable mirror mode. It breaks the quasilinear stability and causes further growth. Pressure balance is either restored by dissipation of the pairs and their anisotropy or inﬂow of plasma from the environment. In the ﬁrst case new pairs will continuously form until equilibrium is reached. In the ﬁnal state the fraction of pairs can be estimated. This process is open to experimental veriﬁcation. To our knowledge it is the only process where in high temperature plasma pairing may occur and has an important observable macroscopic effect: breaking the quasilinear limit and generation of localised diamagnetism. Keywords. Mirror modes, magnetosheath, solar wind, turbulence, plasma diamagnetism 1 Introduction There seems to be nothing particularly interesting left about a very low frequency effect in high temperature magnetised plasma known as the mirror mode (see, e.g., Tsurutani et al., 2011, for a more recent observational review). It was formally discovered some sixty years ago (Chandrasekhar, 1961; Hasegawa, 1975; Gary, 1993) as a theoretical complement to the zero-frequency hose instability, two purely growing linear instabilities in the presence of pressure anisotropies. The hose instability excites propagating Alfvén waves when the magnetically parallel temperature T > T exceeds the perpendicular temperature, the k ⊥ mirror mode grows under the opposite condition T < T that the perpendicular temperature is higher than the parallel by a k ⊥ certain amount, passing a threshold. The mirror mode generates magnetically elongated magnetic bottles in pressure balance thereby providing the plasma a local texture. In the presence of weak plasma gradients the mirror mode formally assumes a small but ﬁnite real frequency (Hasegawa, 1969). Various properties of the instability were added under different plasma conditions and in different wavenumber ranges like ﬁnite gyroradius effects, dependencies on electron temperature and electron arXiv:1909.10163v1 [physics.space-ph] 23 Sep 2019 anisotropies, as well as on plasma convection. The instability saturates quasilinearly at a rather low level by exhausting the bulk thermal anisotropy (cf., e.g., Treumann & Baumjohann, 1997; Noreen et al., 2017, for the quasilinear numerics). Finally, the trapped particle components give rise to the excitation of ion cyclotron and electron whistler waves if only a thermal anisotropy of resonant particles evolves. Identiﬁcation in real plasmas became possible when measuring the pressure balance between external magnetic ﬁeld and increased internal plasma pressure. In observations the two pressures in ion-mirror modes are anti-correlated, a condition which generally serves as the key identiﬁer of mirror modes. The disturbing barely understood remaining point is how in an ideally conducting plasma at high temperature the magnetic ﬁeld can become expelled to a high degree from the interior of the magnetic mirror bottle, an effect resembling the Meissner effect in low-temperature superconductivity which, however, is forbidden in classical physics as it requires the presence of quantum correlations known to be restricted to very dense and low-temperature conditions only. In superconductivity the Meissner effect arises from the electron-phonon interaction in the crystal lattice when spin-compensated electron-Cooper pairs (Cooper, 1956) form, occupy the same quantum state and, though remaining fermions, together with the interacting phonons assume bosonic properties. These quasi-particles condensate in the lowest energy level above Fermi energy and evolve into a Landau-Fermi ﬂuid (Landau, 1941; Ginzburg & Landau, 1950) culminating in BCS theory (Bardeen et al., 1957). The condensate current becomes capable of generating up to ∼ 100% bulk diamagnetism rather different from pressure balance which in BCS is warranted by the stiffness of the crystal lattice. In mirror modes the quasi-diamagnetic effect amounts to roughly 50%, many times more than the ∼ 1% magnetic amplitudes quasilinear theory predicts (cf., e.g., Noreen et al., 2017). For its explanation no weak kinetic plasma turbulence theory (cf., e.g., Yoon, 2007) is in sight. This discrepancy suggests that in the evolution of mirror modes some fundamental effect is still missing. Recently we exhumed kind of a parallelism between superconductivity and the growth of the ion mirror mode (Treumann & Baumjohann, 2018a). Here we demonstrate that the mirror mode can be understood as a combination of the classical plasma ion effect which generates magnetic bottles at the low quasilinear saturation level, while the main large mirror effect may possibly be caused by the trapped bouncing particle component in a similar though classical way as in the BCS theory. In this case it are the mirror- ing particles (electrons and possibly also ions), which interact with the available general thermal ion-sound wave population in the plasma at either thermal or non-thermal level to produce trapped pair singlets which are dynamically distributed over the volume of a mirror bottle. Together with the ion sound ﬂuctuations they form kind of quasiparticles which become locked near their centre-of-mass mirror points, condensate at perpendicular energy close to the electron temperature. This leads to electron boundary currents causing some weak diamagnetism which is probably incapable of explaining the further growth of mirror modes. However, evolving from the quasilinear stable state, the pairs generate a large perpendicular electron anisotropy that enters the mirror growth and breaks the quasilinear stability. Further mode growth is either pressure compensated by dissolu- tion of pairs and plasma heating or by quasi-neutral plasma inﬂow from the environment along the magnetic ﬁeld. This kind of pairing is proper to the mirror instability and may develop for both ions and electrons. For being speciﬁc we concentrate on electrons only in the following. The physical mechanism behind this effect is the interaction of the trapped electrons with the thermal background noise of the ion-acoustic wave spectrum excited in the mirror unstable plasma at low but sufﬁciently large amplitude. This noise is always 2 present at measurable intensity, for instance in the magnetosheath (cf., e.g., Rodriguez & Gurnett, 1975; Treumann & Baumjohann, 2018b, for additional arguments concerning mirror mode observations). 2 Electron trapping Once the ion-mirror mode starts growing at the well known (cf., e.g., Noreen et al., 2017) ion-mirror growth rate γ (k) k λ β T k k i k m e⊥ ≈ A + A − (1) ω 1 + A π T k β ci i⊥ ⊥ T T i⊥ e⊥ A ≡ − 1 & 0, A ≡ − 1 T T ik ek with approximately vanishing real frequency ω ≈ 0, neglecting the effect of density gradients which would cause a ﬁnite real part on the frequency, λ = c/ω ion inertial length, a magnetic bottle evolves in slightly oblique direction k ≪ k with i i k ⊥ magnetic disturbances |δB | ≫ |δB |. This bottle is elongated along the ambient magnetic ﬁeld B and has a narrow opening angle θ given by tanθ = k /k ≪ 1. Instability corresponds to a second order phase transition in plasma which happens whence the magnetic ﬁeld locally drops below a critical threshold value B < B ≈ 2µ NT A |sinθ| (2) c 0 i⊥ where we neglected the electron contribution. Though substantial, this growth rate is just a fraction of the ion cyclotron fre- quency ω = eB/m for k λ ≪ 1 and B respectively A near threshold, the usual case (Treumann & Baumjohann, 2018a). ci i i As noted above, the instability readily stabilises quasi-linearly (Noreen et al., 2017; Treumann & Baumjohann, 1997) at very 2 2 low level |δB| ≪ B via depleting the anisotropy A. One may note that a marginally stable case is obtained for 2 2 Aβ = k /k (3) More interesting is that for initially negligible electron contribution A ≪ A at quasilinear stability the electron anisotropy may become important to surpass the stable quasilinear state in which case the ion mirror instability may further grow. Thus it is of particular interest whether a substantial electron anisotropy A can self-consistently be generated in the mirror mode as this could cause growth of the ion mirror mode to large amplitudes. In the following we propose such a mechanism. We may, however note in passing, that this same mechanism can also itself drive an electron mirror mode after quasilinear stability of the ion mode has been reached. This is a short-scale magnetic ﬂuctuation on the ion mode background with growth rate γ k e ek ≈ k λ A − (4) e e ω π k β ce e⊥ a fraction of the electron cyclotron frequency ω . The marginal stability condition requires vanishing of the bracket. The ce electron mode also stabilises quasi-linearly causing short wavelength depletions of the magnetic ﬁeld on the course of the above ion mirror mode as have been observed previously (cf., e.g., Treumann & Baumjohann, 2018b). In the following we do not further investigate this electron mirror mode even though it will also be affected by the proposed electron pairing as it might be responsible for its surprisingly large amplitudes. Our focus here is on the main ion mode. 3 2.1 Electron dynamics, energy limit, trapped density fraction The conventional ion mirror mode provides a quasi-stationary magnetic bottle (see, e.g., Constantinescu, 2002, for a geo- metric analytical model) structure, which necessarily traps electrons of sufﬁciently small magnetic moment µ . Because the mirror mode frequency practically vanishes and the mirror mode grows slowly compared to the electron dynamics, electrons react adiabatically to the presence of the mirror instability. They conserve their magnetic moment µ = E /B = const when e e⊥ moving along the magnetic ﬁeld B(s). Trapping occurs between the two mirror magnetic ﬁelds B = B(±s ), with s co- ±m m ordinate along the magnetic ﬁeld. The trapped-electron perpendicular kinetic energy E (s) = µ B(s) ≡ V (s) plays the role e⊥ e of a retarding potential E (s) = E − V (s) (5) ek e At the mirror points the parallel energy of trapped electrons vanishes, and E (±s ) = E . Thus trapping occurs for all e⊥ m e µ ≤ µ ≡ E /B , a well-known fact. Though it does not bunch them, mirroring keeps these electrons together by conﬁning e m e ±m them to the volume of the bottle, inside which they perform the oscillatory bounce motion between mirror points ±s . The parallel electron equation of motion is dv (s,t) µ ∂ k e = − ∇ B(s), ∇ = (6) k k dt m ∂s For symmetric bottles and motion around and not too far away from the minimum B(s ) = min{B(s)} ≡ B of the magnetic 0 0 ﬁeld we have ∂ B 1 ′′ 2 ′′ B(s) ≈ B + B (s− s ) , B = (7) 0 0 0 0 2 2 ∂s s which immediately gives the bounce frequency ′′ ω = B ≪ ω (8) b ce a frequency much less than the electron cyclotron frequency ω = eB/m . We shall show below that this kind of trapping, in ce e the case of the mirror mode, becomes advantageous for electron pairing, an effect otherwise observed only under solid state conditions in superconducting metals. In order to get an idea on the trapping energy condition we consider the mirror point s = ±s . Here all the energy is in the perpendicular direction, i.e. the local gyro motion of electrons. Hence de-trapping of electrons occurs once their gyroradius exceeds the opening radius of the bottle neck r & R = L tanθ, where L ∼ 2π/k is the half length of the ion mirror ce,s s k k k bottle. This yields immediately that electrons remain trapped as long as their energy satisﬁes the condition 1 T E . ≡ E , r = v /ω (9) e trap ce e ce,s 2 2 2 4π k r ce with r is the electron gyro-radius and k the perpendicular wave number of the ion-mirror mode. All electrons with this ce ⊥ energy remain trapped in the magnetic mirror bottle. Larger energy electrons escape from the bottle along the magnetic ﬁeld. (We do not discuss the subtle problem that quasi-neutrality requires them to be replaced by low-perpendicular-speed electron inﬂow along the magnetic ﬁeld.) 4 2.2 Fractional trapping conditions The fractional number density of maxwellian mirror trapped electrons is trap N ǫ trap e = C dǫ ǫ exp − e e N T 0 e trap −1 3 3 = Γ γ , (10) 2 2 r r 1 E E E trap trap trap = erf − exp − 4 T T T e e e where C is a normalisation constant, and γ(a,b) is the incomplete Gamma function. Constancy of the trapped electron magnetic moment (cf., e.g., Baumjohann & Treumann, 1996, for a textbook presentation), yields the parallel energy B(s) E (s) = E 1− (11) k e B(s ) which deﬁnes the angle between velocity and magnetic ﬁeld for the trapped electrons B(s) −1 θ(s) = cos 1− (12) B(s ) These or the equivalent expressions we will need below. 3 Single electron wake potential In preparing for the investigation of ion-mirror-mode trapped electrons, we consider the interaction of an electron with the bath of ion sound waves. This is most easily done in the naked test particle picture, assuming that we grab one of the electrons and ask for its reaction to the presence of the dielectric in which it moves. The electron is a point charge −e with velocity v that is located at its instantaneous position x = x− vt in the observers frame (x,t). This is represented by the point charge density function N(x,t) = −eδ(x− vt). We assume that the electron is non-relativistic which for trapped electrons under the conditions in the magnetosheath (Lucek et al., 2005) or the solar wind is good enough. The relative dielectric constant of the plasma it experiences is ǫ(ω,k) where ω,k are frequency and wavenumber of the plasma wave which changes the dielectric properties. In general, we have a whole spectrum of waves which is taken care of below by integrating over the entire spectrum. The naked charged electron polarises the plasma. The total electric potential the moving non-relativistic charge at location x causes at location x is obtained from Poisson’s equation with above charge density and has the form e δ(ω − k · v) ik·(x−vt) Φ(x,t) = − dωdk e (13) 3 2 (2π) ǫ k ǫ(k,ω) This can easily be shown (originally given by Neufeld & Ritchie, 1955; Krall & Trivelpiece, 1973, for a textbook description) by Fourier-transformation. One may note that this expressions also holds for ions. In this representation the action of the 5 δ-function on the exponential has been taken care of. Integration is over wave numbers and frequencies, the wave spectrum responsible for the dielectric properties experienced by the test electron. Integration with respect to frequency ω implies ω → k · v also in the dielectric response function ǫ(ω,k), which we shift until having discussed the latter. In solid state physics it is assumed that the oscillations of the ion lattice generate a thermal spectrum of phonons. In plasmas these waves are not restricted to the Brillouin zones but are freely propagating waves either forming a thermal background noise or, for T ≫ T providing a broad spectrum spectrum of unstably excited ion sound for which the plasma response function e i accounts. In high temperature thermal equilibrium it is the jitter motion of electrons which leads to spontaneous emission of sound and modiﬁes the dielectric properties of the plasma. The general electrostatic response function reads ǫ(ω,k) = 1 + + χ (ω,k) + χ (ω,k) (14) e i 2 2 k λ with χ (ω,k) the electron and ion susceptibilities. Under nonlinear conditions the susceptibilities depend in addition on the e,i wave amplitude. This electrostatic dispersion relation contains both the effects of electrons and ions. One may wonder why for wavelengths usually much longer than the Debye length λ ≪ λ the second term in this expression is not neglected. The reason that it must be retained here, is that the uncompensated charge of the test particle when immersed into the plasma excites short wavelengths waves on the Debye scale in order to screen the charge. Therefore, independent of the wavelength of plasma waves, the test particle dielectric response must include the Debye term. The dielectric response function of the thermal spectrum of ion-sound waves at frequencies far below the electron plasma frequency ω ≪ ω is 1 ω ǫ(ω,k) = 1 + − (15) 2 2 k λ ω where ω is the ion plasma frequency, λ ≈ v /ω the Debye screening distance, and the frequency of ion sound waves ω is i e e e k obtained putting the real part of this expression to zero, which as usually yields 2 2 2 2 ω k λ c k e 2 s 2 = or ω = < ω (16) k i 2 2 2 2 ω 1 + k λ 1 + k λ e e 2 2 2 2 Here c = ω λ ≈ (m /m )v ≈ 2T /m is the ion-sound speed square. It is simple matter to show that the inverse response e i e i s i e e function becomes 2 2 2 1 k λ ω = 1 + (17) 2 2 2 ǫ(ω,k) 1 + k λ ω − ω Actually, this is also the general inverse form of any dielectric response, if only ω is understood as the solution of the general 2 2 2 response function ǫ(ω,k) = 0 for electrostatic waves, and ǫ(ω,k)−k c /ω = 0 for very-low frequency electromagnetic waves like magneto-sonic or (kinetic) Alfvén waves. In the latter case one has (cf., e.g., Treumann & Baumjohann, 1997), including kinetic effects, −1 1 c 2 3 T ǫ (ω,k) = 1 + + 1 + k · r + (18) A ci 2 2 2 k λ V 4 T e A 6 2 2 2 with r = v /ω the vectorial ion gyro-radius. The relevant wave frequency is ω ≈ k V for the ordinary Alfvén wave, ci i⊥ ci kA A with V ≪ c the Alfvén speed (if wanted including the bracketed modiﬁcation factor). The kinetic Alfvén wave propagates fast along the magnetic ﬁeld almost at Alfvén speed, and roughly ten times slower perpendicular to it. The weak wave electric potential resulting from its kinetic nature is along the magnetic ﬁeld and is being believed to be responsible for electron acceleration. Any attractive pairing effect in resonance with those fast parallel electrons will be in this direction as well, a very interesting fact in itself which we do not investigate here, leaving it for a separate investigation. It may be applicable to particles in the auroral magnetosphere with its magnetic trapping conﬁguration. Inserting Eq. (17) into the above electrostatic potential of the test electron eλ dω k dk dk ⊥ ⊥ Φ(x,t) = − × 2 2 2 (2π) ǫ 1 + k λ k ik·(x−vt) × 1 + δ ω − k · v e (19) 2 2 ω − ω shows that Φ consists quite generally of two contributions, the screened Coulomb potential of the test electron, and another wave induced term which multiplies the screened potential by the frequency dependent term in the last expression. (We note again that a similar form trivially holds as well for ions.) This form demonstrates the well known self-screening Debye effect of the naked point charge, which leads to the ﬁrst term in the above expression and causes the Debye-Yukawa potential to exponentially compensate for the electron charge ﬁeld in a spherical region of radius λ . We are not interested here in the deformation of the Debye sphere introduced by the electron motion as this is a higher order effect. The zero order effect of the test electron contained in the wave-independent term, the proper self-screening is, in the wave- 2 2 2 2 dependent term, multiplied by the wave-induced factor. For frequencies ω = k · v > ω higher than ion sound, this factor 2 2 2 2 is positive adding to the screening but changes sign for frequencies ω = k · v < ω , thereby indicating the possibility of over-screening at wavelengths larger than the Debye radius λ (cf. Treumann & Baumjohann, 2014, their Fig. 1). Under certain conditions it may come into play outside the Debye radius where the charge-electric ﬁeld is practically already compensated, and the long range wave electric ﬁeld adds up over some distance, may dominate and cause a spatially restricted deﬁciency of repulsion. In this case the potential may even turn negative, eliminates the repulsive nature of the electron locally and becomes attractive for electrons. This was ﬁrst shown (Neufeld & Ritchie, 1955) for high frequency Langmuir waves even before the discovery of Cooper pairs in superconductivity and solid state physics. In a bath of Langmuir waves this attraction turned out to be unimportant however, while in an isotropic non-magnetic plasma it survives for low-frequency ion sound, ﬁrst suggested by Nambu & Akama (1985). With θ the angle between electron speed and wavenumber, it happens at resonant electron speeds 2 2 2 2 v cos θ . ω /k (20) requiring the parallel electron speed to be less than the wave phase velocity. The above expression depends on angle θ between velocity v and wavenumber k, which in our case will turn out to be of crucial importance. For completeness we note that in magnetised plasma the ion acoustic wave is azimuthally symmetric with respect to the magnetic ﬁeld B. However, its frequency depends itself on the angle of propagation between k = (k ,k ) and B according ⊥ k 7 to (Baumjohann & Treumann, 1996) 2 2 c Λ (η )k 0 i s k ω = (21) 2 2 Λ (η ) + k λ 0 e 1 2 2 with Λ (η ) = I (η )exp(−η ), η = k r , and the index j = e,i on the gyroradius is for electrons and ions. I (η ) is the 0 j 0 j j j 0 j cj 2 ⊥ Bessel function of imaginary argument. r = v /ω is the gyroradius, and ω is the cyclotron frequency. One has that, cj ⊥,j cj cj moreover, k λ ≪ k r ≪ 1 and k /k < 1. Long-wavelength ion sound in magnetised plasma thus propagates essentially ⊥ e ⊥ ci k ⊥ along the magnetic ﬁeld, a well known fact which in observations, for instance in the magnetosheath (Rodriguez & Gurnett, 1975), manifests itself as a complete drop out of the electrostatic low frequency thermal ion noise spectrum when the antenna points strictly perpendicular to the ambient magnetic ﬁeld (cf., e.g., Treumann & Baumjohann, 2018b, for an example and discussion). The interaction between electrons and ion sound waves thus opens up the option that electrons in a Debye-screened potential may, under certain conditions, experience an attractive potential which compensates and overcomes the Coulomb repulsion between two negatively charged electrons, resembling a the famous effect of Cooper pairing in solid state physics though here in the realm of classical physics. The paired electrons and the propagating ion sound wave form a quasiparticle in both these cases. It is important to insist that this attraction is not due to trapping of the electron by a large amplitude wave in the wave potential trough; on the contrary, it is an electron-induced change in the dielectric properties of the wave-carrying plasma causing the electron to evolve an attractive electrostatic wake potential which it carries along when moving across the plasma. We have previously shown (Treumann & Baumjohann, 2014) that this can happen also with other waves than ion-sound. Below we demonstrate that it becomes crucial in the evolution of mirror modes to which plasma wave trapping does not contribute in no sense. Since the waves are propagating along the magnetic ﬁeld and the bounce motion of the electrons is as well along the magnetic ﬁeld, the coordinate s of interest is parallel to the magnetic ﬁeld sˆkB, and the gyration of the electrons decouples from the interaction. In this case we have for the wave number k = (k ,k ) and velocity k ⊥ k ≡ k · sˆ, v (s) ≡ v · sˆ = v cosθ(s) k k parallel to the local magnetic ﬁeld. The problem then consists in solving Eq. (19) under the conditions of a bouncing test electron. This task resembles the solution under non-magnetised conditions which had been given in our previous paper (Treumann & Baumjohann, 2014). In the known form it cannot be applied here but has to be substantially modiﬁed in or- der to become adapted to the conditions of electron trapping in mirror modes. 8 3.1 Conditions for an attractive potential In the light of the previous discussion we rewrite Eq. (19) in the magnetic ﬁeld as 2 ik·(x−vt) eλ ω dω dk e Φ(x,t) = − × 2 2 2 2 2 2(2π) ǫ (1 + k λ + k λ ) ⊥ e e δ(ω − k v ) δ(ω − k v ) k k k k × − (22) (ω − ω ) (ω + ω ) k k Here we left the Debye-potential term out as it is of no interest, and resolved the denominator. We also refer to the parallel particle velocity v = v cosθ which in our case of mirror trapped test particles is along the magnetic ﬁeld. It selects the parallel wavenumber of the wave in the Dirac δ-function to replace the frequency ω. In the same spirit the argument of the exponential becomes ik · (x− vt) = ik ρsinφ + ik (s− tv cosθ) with ρ the independent perpendicular spatial coordinate. It is assumed ⊥ k that the magnitude v of the velocity remains constant in this kind of interaction, which holds for the adiabatic motions along the magnetic ﬁeld where no further external force acts on the electron except for the stationary restoring magnetic force. (Note also that the wave frequency ω depends on k ,k but not anymore on angle φ because it has been determined independently k k ⊥ from kinetic wave theory not using the test particle picture.) These assumptions reduce the integral to integrations over the perpendicular wavenumber k ,φ, and frequency ω. Moreover, since the problem has become cylindrically symmetric with respect to B, integration over φ can easily be performed by using the representation of the exponential as a series of Bessel functions (Gradshteyn & Ryzhik, 1965) which reduces to the zero- order Bessel function J (k ρ). The formal result before ﬁnal integration is 0 ⊥ ω dω k dk dk eλ k ⊥ ⊥ Φ(s,ρ,t) = − × 2 2 2 2 2 2(2π) ǫ (1 + k λ + k λ ) ⊥ e e δ(ω − k v ) δ(ω − k v ) k k k k × − × (23) ω − ω ω + ω k k ik (s−v t) k k × J (k ρ)e 0 ⊥ where one understands v = v cosθ(s). We note again that this form is still valid also for electromagnetic waves if only the frequency is understood as the solution of the electromagnetic dispersion relation. In view of later application to the mirror mode we now restrict to purely electrostatic waves, in our case ion sound which, in contrast to other electrostatic waves like Bernstein modes, have the right property of propagating along the magnetic ﬁeld with parallel electric ﬁeld. The ion sound wave frequency is 2 2 2 2 Λ (η )c k Λ (η )k c 0 i 0 i s s k k ω ≈ ≈ (24) 2 2 2 2 Λ (η ) + k λ 1 + k λ 0 e e e with the right-hand side holding since the electron term in the denominator is Λ (η ) ≈ 1. In the low frequency approximation 0 e applicable here, the frequency is proportional to the parallel wave number. In the following we simplify this dispersion relation setting Λ (η ) ≈ 1, which is its maximum value, and in the resonant denominators neglecting the inverse dependence of ω on 0 i k kλ , only keeping it in the nominator of the integral. Then one may perform the integration with respect to k which gives, e ⊥ 9 ′ with ξ = λ k , ρ = ρ/λ , ζ = k λ , e ⊥ e e ∞ p ′ ′ exp − ρ 1 + ζ ξdξJ (ξρ ) I(ρ ,ζ) ≡ = (25) 3/2 2 2 1 + ζ 1 + ζ + ξ In order to perform the integral, its singular properties have to be elucidated. The dominant contribution will come from the resonant denominators in the bracketed terms. Any possible resonances in the Coulomb factor do not play any role here. The Dirac δ-functions prescribe replacing the frequency everywhere with k v . It is, however, convenient to delay this action until k k integrating out the singularities in the complex ω plane. To see their effect, one temporarily replaces k in the argument of the exponential with ω as the δ-function prescribes as an inverse action. Then we have for ik (s− v t) = iω(s/v − t). Since k k k the waves are damped, the imaginary part of the frequency is required to be negative. This forces demanding s/v − t < 0, consequently taking the ω-integration over the lower complex ω-half plane, which in surrounding the poles in the positive sense adds a factor 2πi to the integral and includes the sum of residua ω = ±ω in the integral in this order. The result is Z ′ −ρ 1+ζ +iζ(s−v t)/λ k e ie c ζ dζ e Φ(s,ρ,t) = − p 4πǫ λ v 0 e k 1 + ζ ω λ ω λ k e k e × δ ζ − − δ ζ + (26) v v k k Performing the substitution prescribed by the delta functions in the exponential only yields the sum of two exponentials which turns into a sine function. One then obtains for the potential of the particle in the presence of ion sound waves Z ′ h i −ρ 1+ζ e c ζ dζ e Φ(s,ρ,t) = p sin ζ(s− v t)/λ (27) k e 2πǫ λ v 0 e k 1 + ζ What remains is the ζ integration with ζ = k λ < 1 limited. To simplify, we can either neglect ζ or replace it by unity in the k e arguments of the roots. To be conservative and decide for the weakest case, we chose the latter, what yields the integral h i e c − 2ρ 2 Φ(s,ρ,t) = − √ e dζ sin ζ s− v t /λ (28) k e 4 2πǫ λ 0 e The argument of the sine function is negative. So we have taken its sign out and use its absolute value. Integration gives − 2ρ/λ e c e Φ(s,ρ,t) = − √ sin|σ|−|σ|cos|σ| (29) v |σ| 2 2πǫ 0 k where −1 σ = (v t− s)λ > 0 The condition for an attractive potential follows immediately as tan|σ| > |σ| or 0 < σ < mod (2π) (30) 10 Depending on the parallel velocity v > 0 there is an entire range of distances s < v t < πλ /2 in which the conditions for k k an attractive potential are satisﬁed. We may note that for negative velocities v < 0 there is no range where the potential can become attractive as the braced expression is always positive. It is the scalar product k · v between the wave number of the ion-sound and the test particle velocity which selects those speeds which are parallel to the sound velocity, not anti-parallel. One should keep in mind that this attraction has nothing in common with wave trapping, however! It is the over-screening effect of the particle, which is moving on the background of the wave noise and experiences the modiﬁed dielectric properties of the plasma. It should also be noted that in this condition the time explicitly appears because the test electron is seen from the stationary observers frame in which the electron moves. Instead, σ is measured in the moving electron frame. This distinction is important to make as it will be picked up again below. 2 2 The restriction on the velocity is obtained from that ω < ω when referring to the replacement ω = k v = k v cosθ pre- k k k scribed by the δ-functions. Rescaling ω ∼ c k , it follows that the parallel particle speed is limited as k s k |v | . c or |cosθ| . (31) This is in fact a condition on the angle θ. For small speeds v < c the condition is trivial. The largest effect is caused when the particle speed is parallel, below and close to the phase velocity c of the ion-sound wave. For large velocities v > c the s s angle between the phase speed and velocity must be close to π/2, in agreement with the above requirement on the potential becoming attractive. This is an important point in application to a plasma. In thermal plasmas we have generally c ≈ m /m v which is far s e i e below the thermal speed. Hence there are only few electrons in the distribution sufﬁciently far below thermal speed which would satisfy the resonance condition v < c . Higher speed electrons can be in resonance and thus contribute to attraction only at strongly oblique wave and electron speeds. Consequently under normal conditions in a plasma the generation of attractive potentials becomes obsolete, a point which had been missed in previous work (Neufeld & Ritchie, 1955; Nambu & Akama, 1985). In the particular case of mirror modes it becomes the crucial ingredient, as will be demonstrated below. 3.2 Correlation length In all cases the attraction exceeds the repulsion outside the Debye sphere of the electron in its wake and, therefore and most important, can be felt by other electrons. From here it is clear that two electrons must move at distance somewhat larger than λ and at nearly same speed in the same direction in order to be held together by their attractions and form a pair. This is the important point when applying our model to the mirror mode below. Having obtained the conditions under that the wake potential behind the moving test electron becomes attractive, we would like to know the distance over that the negative potential extends. This distance is measured in the instantaneous frame of the electron and is, hence, given by the above absolute normalised value of |σ| < π/2 which repeats itself periodically. It is, however, clear that it is only the zeroth period which counts as the effect of the dielectric polarisation on the electron diminishes 11 ′ with increasing distance s = σλ . In absolute numbers this distance becomes λ = |s− v t| < λ ≈ 1.57λ (32) corr e e which can be understood as an electron “correlation length” between neighboured electrons. Any electrons within such a distance will behave about coherently, an important conclusion which, however, has to be extended below to many electrons. This correlation length is to be compared with the particle spacing in the plasma. Plasmas are deﬁned for particle densities Nλ ≫ 1, which implies that the distance between the particles is ≪ λ . Consequently the extension of the attractive potential in the electron wake is much larger than the spatial distance between two electrons. It thus affects many electrons, an effect which cannot be neglected when speaking about attraction. As for an example, in the magnetosheath which is the preferred domain where the mirror mode is permanently excited, the 7 −3 average density is, say, N ≈ 3×10 m at temperature T ≈ 50−100 eV. For the Debye length we have λ ≈ 10 m, while the e e inter-particle distance is a mere of ≈ 0.005 m. Roughly ≈ 10 electrons should experience the presence of the attraction behind the test particle, which thus becomes a many-electron effect. Because pair formation depends on the quite severe condition on the particle velocity, not all those electrons of course will form pairs though however, in reality, the attractive potential involves a substantial fraction of electrons which necessarily will cause modiﬁcations of the plasma conditions. Normally such modiﬁcations will only cause minor effects in the wave spectrum and will be negligible. Below we show that in the evolution of the mirror mode they become important. 3.3 Ensemble averaged potential If we understand the plasma as a compound of a large number of electrons, we can ask for the ensemble averaged potential hΦi of the single electron averaging over the particle energy distribution. In an isotropic plasma this is the Boltzmann distribution. Writing for the parallel velocity v = v cosθ the average potential becomes e Cc 2 2 − 2ρ/λ −v /v Φ = √ e vdv e × 0 2λ (33) (s+πλe/2)/tv dcosθ × sinσ − σ cosσ σ cosθ s/tv which immediately tells that the mean potential taken over the full Boltzmann distribution in repulsive. This is clear, however, because it accounts for all electrons in the Debye sphere. To calculate the cos-integral we expand the trigonometric functions to obtain π π π π vt σdcosθ ≈ vt− slog 1 + (34) 6 6λ 12 2 s 12 We now exclude the Debye sphere by restricting the integration with respect to v over a shell between the thermal and trapped speeds. This gives trap r r π t 2E πt 2E −E/T dEe − slog 1 + ≈ 12 m 2s m e e (35) π πt 2T 1 e −x − 1− x e dx 6 2 m For a mean attractive potential the last integral should be positive. Doing it yields (Gradshteyn & Ryzhik, 1965) 1 2 −x 2 −y −1 3 2 x e dx = y e − e ≈ −y + 2y − 1 (36) which is positive only if y = 1 + Δ and Δ = E − T /T < 1 in which case there is a narrow energy range (or energy trap e e “gap”) for trapped electrons where the mean potential Φ < 0 becomes attractive for the electrons when averaging over their energy distribution and warranting that they behave coherently. The latter we will show can under certain condition be the case. 4 Two-electron potential We saw that, under a certain condition, an electron moving in the plasma in resonant interaction with an ion-sound background may give rise to an attractive potential in its wake where another electron can be captured and thus be forced to accompany the ﬁrst electron. First of all, in plasma all electrons are in permanent motion. Hence, if an electron satisﬁes the resonance condition with an ion sound ﬂuctuation, it acts attracting on another one moving nearly at same speed. We have seen that this attractive potential in the presence of a large number of thermally distributed electrons becomes depleted. This holds when just one electron contributes to the potential. We now extend this to the combined effect of two electrons in the interaction, in which case we can immediately use the above solution when, however, accounting for the slightly different velocities v ,v k1 k2 and initial locations s ,s of the electrons along the magnetic ﬁeld. In view of the later application to mirror modes, we again 1 2 consider only motion along the magnetic ﬁeld not yet specifying to the pecularities introduced by bouncing in the mirror ﬁeld. Then the two-electron potential becomes − 2ρ/λ X e e c e Φ(s,ρ,t) = − √ × v |σ | 2 2πǫ j 0 kj (37) n o × sin|σ |−|σ |cos|σ | j j j with j = 1,2 counting the electrons. Here −1 σ = (v t− s )λ > 0 j kj j 13 As before, the requirement σ > 0 results from the condition that the waves in resonance with the electrons must be damped. In order to obtain the combined effect of the two electrons, we transform to their centre-of-mass frame 2Z = s + s , 2z = s − s 1 2 1 2 (38) 2U = v + v , 2u = v − v k1 k2 k1 k2 From the previous we saw that the large correlation length implies that many electrons are affected. Any attractive potential couples two ore more particles together. The most probable state to be formed is the two-particle (singlet) state. These will be distributed over the plasma, resembling the Cooper states in solid state superconductivity while not being a quantum effect here. Rather it is the polarisation effect moving particles produce in the high temperature collisionless plasma which causes singlet states of pairs. In order to be realistic, we now derive the condition for singlet states to evolve. To simplify the algebra, let us deﬁne 1 −1 Σ =: σ + σ ≡ Ut− Z λ > 0 1 2 (39) ′ 1 −1 σ =: σ − σ λ ≡ ut− z λ 1 2 e The restriction on Σ > 0 maps the ω-resonance onto the new variables. At the contrary, σ can be positive or negative. With these expressions and after some rather tedious though simple calculations, Eq. (37) can be brought into the form − 2ρ/λ 2e c t e Φ(s,ρ,t) ≈ − × π ǫ (λ Σ + Z) |Σ| 0 e (40) n o × sin|Σ|−|Σ|cos|Σ| cosσ where we made use of the above representations and replaced 1 ′ v t = λ Σ± σ + Z ± z (41) k1,2 This expression holds under the reasonable assumptions σ ≪ Σ and z ≪ Z that the difference between the two electrons in location is small enough to be found within the correlation length. Only under this condition one expects that the electrons will be correlated. Interestingly, the form of the potential remains the same as that for the one-particle case with the only difference ′ ′ that the potential is multiplied by cosσ . Hence the condition for attraction depends on the value of σ . Closely spaced electrons of similar and, as required, resonant speeds not differing too much from the phase speed of the ion-sound give indeed rise to attraction between the two electrons if the following conditions are satisﬁed: tan Σ > Σ if cosσ > 0 tan Σ < Σ if cosσ < 0 which yields 0 < Σ < if cosσ > 0 (42) − < Σ < 0 if cosσ < 0 14 ′ These conditions are essentially the same as in the one-electron case. There modiﬁcation is due to cosσ being positive or negative and that they apply to the centre of mass coordinate system Z and mean particle speed U which both are contained in the variable Σ. We remark that these conditions are very general. They substantially generalise the conditions found earlier by Nambu & Akama (1985) to the much more important interaction between two electrons, the lowest order singlet state and thus most realised state in a plasma. Higher order states like interaction of three electrons leading to triplets and so on are in principle also possible but will not play any important role because the interaction decays with distance even though they may be located within the correlation length and form “quasi-particles”. In the singlet state the electrons behave like one particle of double charge and double mass for the time of their interaction, the time they remain inside one correlation length. This length for the singlet is the same as given above that produced by one electron with the difference that it now applies to the centre of mass of the two electrons. Measured from the centre of mass it extends to its both sides over a length of roughly λ ≈ 1.5λ . corr e In physical units the ﬁrst singlet state, for instance, is realised for π π 0 < Ut− Z < λ , |ut− z| ≪ λ (43) e e 2 2 In these cases resonant electrons in the presence of an ion-sound wave background will arrange into loosely bound electron pairs. In high temperature plasma a substantial number of such pairs will exist. However, they will mostly not play any role in the dynamics. In order to do so the plasma must offer additional ways for the bound singlet pair states to cause any susceptible effect in the plasma. Such conditions are provided by the quasilinearly stable mirror mode. 5 Mirror bottle and pairs Being in the possession of the conditions under which electron pairs can form in a high temperature plasma in interaction with a thermal background of ion sound waves, we now intend to apply them to the case of mirror modes. We saw that the correlation length between electrons provided by one single test electrons is of the order of λ ∼ 1.5λ . This value is only corr e slightly increased by the interaction of two electrons, such that we can roughly take λ ∼ 2λ for the singlet. A mirror bottle corr e is a preferred place for pair formation. This is in contrast to a spatially extended plasma. Firstly, the bottle conﬁnes trapped electrons which cannot easily escape. Secondly, the parallel velocity of a bouncing electrons varies along the mirror magnetic ﬁeld and at some place may get in resonance with the thermal ion-sound spectrum present in the entire plasma volume. If this happens at some location along the mirror magnetic ﬁeld, electrons might form pairs and remain correlated for some time, bunch, perform like bound orbits and thus represent resonant correlated states which due to the correlation become coherent states. 5.1 Centre of mass pair bounce motion The application of these ﬁnding to mirror modes is not an easy task. The electrons perform a complicated bounce motion along the inhomogeneous magnetic ﬁeld with periodically changing bounce velocity and bounce frequency depending on the value of 15 their constant magnetic moment. Under these conditions we need to know the variation of their bounce velocity as function of the location along the magnetic ﬁeld between the mirror points. We moreover need to satisfy the common resonance condition of the pairs with respect to the phase velocity of the ion sound. Since the electron velocity is generally much larger than the latter this immediately suggests that the best conditions for attraction will be found near the mirror points s . There the parallel velocity of the electrons drops to zero, and there will be a certain range Δs at distance s . s where the resonance condition is satisﬁed most easily. Here, near s one expects that attraction will become important. In order to understand this process we thus need to transform to the moving frame of the pairing electrons. For this purpose we use the electron bounce motion to deﬁne the new pair-electron quantities 1 1 M =: µ + µ , µ =: µ − µ (44) 1 2 1 2 2 2 2 ′′ 2 ′′ Ω =: MB /m , ̟ =: µB /m (45) e e 0 0 2 2 2 1 U =: E − Ω Z , E =: E +E − 2MB (46) 1 2 0 2 2 2 1 u =: ǫ˜− ̟ z , ǫ˜=: E −E − 2µB (47) 1 2 0 The mean bounce velocity U of the pair becomes a function of the location Z of the centre of mass along the magnetic ﬁeld. This requires knowledge of its displacement as a function of the bounce phase which again requires solution of the two dynamics of the two electrons. Note the adiabatic constants E,M,µ, Ω,̟. The only variables are the mean and difference velocities U(Z),u(z). In the magnetic mirror symmetry U(t) is the bounce velocity of the trapped electron pair, and Z(t) is its location along the magnetic ﬁeld at time t. The difference speed u(z) is measured in the centre of mass frame relative to Z and U. The mean speed U along the magnetic ﬁeld must be expressed either as function of time t or distance Z. For this to accomplish one needs to solve the parallel equation of motion: dU dZ 2 2 2 = −Ω Z − ̟ z ≈ −Ω Z, U = (48) dt dt which is given in the reasonable approximation of small ̟ z. Obviously the mean speed along the ﬁeld obeys the mean bounce 2 2 equation, an oscillation at frequency Ω. Integrating the bounce equation of motion with U(Z) = 2E/m− Ω Z yields πt 2E Z(t) = Z sin , Z = Ω (49) m m 2t m m e Z is the distance of the centre of mass mirror point reached by the pair at mirror time Ωt = π along the magnetic ﬁeld. m m Symmetric mirror bottles have been assumed, implying time symmetry ±t . For the lag in distance z as measured relative to the centre of mass Z we obtain similarly πt 2ǫ˜ z(t) = z sin , z = ̟ (50) ǫ ǫ 2t m ǫ e with ̟t = π the lag in time in the electron pair to reach the mirror point at relative location z to the mirror location of the ǫ m centre of mass Z . 16 These expressions give the centre of mass and jitter velocities as functions of time π Z πt U(t) = cos (51) 2 t 2t m m π z πt u(t) = cos (52) 2 t 2t ǫ ǫ It is important to remark that at this point we have only transformed to the centre of mass motion. We have not yet determined the paring conditions. The ﬁrst of the former expressions thus gives the two-electron bounce velocity in the centre of mass frame, the second the jitter velocity around this centre. In order to apply to the quasilinear stable mirror mode one must transform to the bounce motion. This can easily be done inside the mirror bottle by using the parallel equation of motion and the magnetic trapping conditions. It implies substituting the parallel bounce solution for t(s). We shall do it later below. Instead, here, it is ﬁrst more important to obtain the condition of pair formation in the two-electron potential. 5.2 Condition for pair formation Electron pair formation proceeds if, in addition to the conditions for attraction which have been given above, the pair electron are in resonance with the ion sound. This condition is non-trivial. We mentioned already that electrons participating in attraction move at speed comparable to the thermal speed v which exceeds c substantially. Under non-mirror conditions pair formation e s will thus barely take place. However, magnetic mirrors as provided by the mirror instability are a rare exception. The resonance condition is in fact not a condition on U(Z) but on the angle between the pair velocity and the direction of the magnetic ﬁeld, as the latter is the direction of the propagation of the ion sound. During the bounce motion the particle velocities are adiabatically conserved. It is only the angle θ(s) that changes along the magnetic ﬁeld. Thus writing U(s) = v cosθ (s) + cosθ (s) , 1 2 assuming that v ≈ v ,µ ≈ µ , we have 1 2 1 2 U + u U − u cosθ = , cosθ = (53) 1 2 v v 1 1 Introducing the mean angles Θ = (θ + θ ),ϑ = (θ − θ ) ≪ Θ we obtain 1 2 1 2 2 2 hUi Θ ϑ Θ = cos cos ≈ cos v 2 2 2 (54) hui Θ ϑ ϑ Θ = −sin sin ≈ − sin v 2 2 2 2 Note that u can be negative as it is measured in the centre of mass frame. The condition of resonance U(Z) . c then reduces to Θ c hUi . c −→ cos . ≪ 1 2 v (55) ϑ hui ∼ ≪ 1 2 v This condition shows that our assumption of about equal magnitudes v ≈ v is not crucial because of the smallness of this 1 2 ratio. It shows moreover that the resonance condition is nicely satisﬁed near the mirror points s where the average angle 17 p Θ/2 ≈ . As for an example, c /v ≈ m /m ≈ 0.023 which shows that the average cosine is very small, and the effective s e e i ◦ ◦ angle Θ/2 ≈ 88.7 is close to 90 . Allowing for a deviation in the ratio of ∼ 0.002 the angular variation would amount to ϑ ≈ 0.2 as obtained from the average jitter velocity hui, as is suggested by the second above condition with sinΘ/2 ≈ 1, which gives the angular spread in case of attraction. This analysis shows that once a mirror bottle evolves there is a narrow spatial range near the mirror points ±s along the magnetic ﬁeld for the trapped electrons to generate attractive potentials in their wake during their bounce motion inside the magnetic mirror trap. This is the range of resonance where the parallel speed matches the slow speed of ion sound. This attractive potential extends over approximately one to two Debye lengths along the magnetic ﬁeld outside the Debye sphere of the acting electrons (roughly some ten meters in the magnetosheath!) whose charge ﬁelds are compensated by the bulk of the surrounding electrons populating the Debye sphere. This length is however much larger than the mutual particle distance. It thus affects a substantial number of electrons which, in case their velocities do not differ much, form pairs within a correlation length which the attractive potential attributes to them. As a consequence, there is a substantial number of paired electrons inside the mirror bottle along the magnetic ﬁeld around all the many mirror points of trapped electrons of different initial angle and velocity. The distribution of those mirror points depends on the (equatorial) pitch angle distribution of the electrons trapped in the ﬁeld minimum B = B at the centre of the mirror bottle. One thus expects that over a certain length along the mirror magnetic ﬁeld an almost homogeneous distribution of electron pairs will evolve. The potential acts to combine two electrons into a pair at location close to Z = s , where the centre of mass velocity of the electrons goes into resonance with U ≈ c ≪ v . This effect attaches the pair to the particular group of ion sound waves s e in resonance with the pair and thus locks the pair to them. In addition there is a small jitter velocity u around the resonance left over which is required to be small in order not to destroy the resonance. Investigation of the stability of pairs in this case is a separate problem which we do not attack in this communication, even though it is important for the further considerations. However, we may note that even if pairs will not be stable for long time not becoming locked completely to the group of ion sound waves with resonant wave numbers, they will stay for some ﬁnite time in resonance and will be replaced by other newly formed pairs when dispersing. Since the electrons and pairs are all indistinguishable, a ﬂuctuating population of pairs will thus always be present if only the conditions for pairing exist. In this sense a mirror mode is the ideal and possibly the only place in high temperature plasma physics where pairing can and probably will occur. We may also point on a similar mechanism for ions. Our basic expression for the potential for ions would be similar to that for the electrons. It moreover holds for any kind of waves while pairing in addition requires that the particle candidates for pairing must be capable of resonating with the waves. This latter condition is rather severe and therefore the necessary condition to have pairing. 5.3 Dynamics of pair population So far we just derived the conditions for pairing in the centre of mass system of the bouncing electrons. It is most interesting what happens after pairing. Naively thinking nothing. Pairs form and may dissociate. Those formed and do not disintegrate may perform a combined bounce motion as prescribed by the centre of mass bounce equation given above. In this case they remain pairs but bounce between their centre of mass mirror points. Thus for some pairs both disintegration and common 18 bounce may be the case. In addition there will be a number of pairs who are stable at least for some longer time. They neither disintegrate nor bounce because they remain stable at resonance with the group of phonons of resonant wave numbers k near Z ≈ s and centre of mass speed U ≈ c . As long as this interaction holds they will not be able to return into bouncing simply m s because the phonons which are independent, not being subject to bounce or any mirror force, do not allow the paired electrons to return into bouncing. These phonons continue their slow motion along the magnetic ﬁeld and in this way lock the pair. Vice versa, they by locking the pair become themselves locked near the pair’s mirror point. Thus such pairs drop out of bouncing and for their life time become locked near s . The length of this time is a question of their stability in which not only the two electrons but also the ion sound phonons are involved, while these are independent of the bounce. Stable pairs still posses a small jitter speed u around s which is insufﬁcient of re-injecting them into the bouncing. They also move together with the phonons further up the ﬁeld at reduced parallel speed. The condition on the jitter energy m|u| for stability is that it should not ∗ 2 exceed the trapping potential m u . e|Φ|. Pairs of smaller jitter energy should thus be stable. A more precise selection rule requires the solution of the stability problem. We consider a mirror bottle. The pitch angle distribution of trapped electrons in a mirror bottle is not known a priori. The equatorial pitch angle θ is given as sin θ = B /B(s ) (56) 0 0 m the ratio of minimum magnetic ﬁeld to the mirror ﬁeld of trapped electrons. Electrons with large equatorial pitch angle mirror very close to the minimum magnetic ﬁeld. Electrons with small equatorial pitch angle mirror near the end of the bottle. It is thus clear that there is practically a continuous distribution of mirror points along the mirror magnetic ﬁeld in the bottle depending on the given initial distribution of equatorial pitch angles. Moreover this applies to all magnetic ﬁeld lines inside the bottle, not only the central one. The dependence of the mirror points s on electron velocity v < E /m , location s, and pitch angle m tr e θ is obtained from the bounce frequency ω , the location along the ﬁeld 0 b s(t) = s sin ω t (57) m b and the time to reach from s = 0 to s within the bounce motion ds/cosθ t(s) = p (58) ′′ 2 1− s (B /2B )tan θ 0 0 Solving the latter integral and resolving for the mirror point yields the wanted expression −1 ′′ −1 0 s (v,s,θ ) = s sin v sin s tan θ (59) m 0 0 2B This is a complicated though continuous dependence on s and θ . The distribution of mirror points on the mirroring electron velocity respectively its energy can, in principle, be integrated over the range of velocities 0 < v−v < v −v contributing e trap e to pairs. Assuming an equatorial isotropic gaussian-velocity distribution ∼ f(θ )exp(−ǫ/T ) yields 0 e π/2 ΔE s dθ sinθ f(θ ) pair 0 0 0 s (s) ∝ p (60) −1 ′′ sin s B /B tanθ 0 0 19 with ΔE = E − T > 0. For any given equatorial pitch angle distribution f(θ ) there is a corresponding distribution of pair trap e 0 mirror points inside the mirror bottle along s, i.e. along the magnetic ﬁeld B(s). As a consequence, the entire narrow volume of the quasi-linearly stable mirror bottle will be subject to the presence of pairs each of which is located at and along the magnetic ﬁeld line centred around its common pair-mirror point. (One may note that the point s = 0 does not contribute to the integral, neither mathematically nor physically because it is not a real mirror point. Also θ = 0 does not contribute because those particles do not participate in mirroring.) Under these conditions the pairs become an important part of the population of a mirror bottle and might contribute to its ∗ ∗ dynamics. The locked pairs have twice the electron mass m = 2m and twice the electron charge q = −2e. Energetically ∗ 2 their parallel motion condensates in the lowest energy bounce level ǫ˜∼ (m /2)hu i ≪ E , the energy in their average jitter trap motion around the mirror point, negligibly small with respect to their gyration energy which near the mirror point has absorbed almost all kinetic energy into the gyration. At same magnetic ﬁeld strength they thus have equal gyroradii concentrating in a current shell which represents a surface current J = −e N v Δr on the time-scale of the bounce which averages over ⊥pair pair e all particle phases. This shell has width Δr ∼ (v /ω )ΔB/B. The surface current it carries might give rise to an integrated e ce orbital diamagnetism. In the last subsection below we explore its effect. The promising macroscopic effect is the direct contribution of pairs to both, the evolution of an electron and the ion mirror modes. A gyrating locked pair population has large perpendicular energy ∼ 2T and very small parallel energy ∼ m |u| . It e e hence possesses a large anisotropy 2T 2T e e A = − 1 ≈ ≫ 1 (61) pair 2 2 m |u| m |u| e e which contributes to the evolution of mirror modes through the growth rate (1) of ions as this growth rate contains the anisotropy of all electron components, including pairs. It is thus clear that a large electron pair anisotropy A will directly contribute pair to the growth of the ion mirror mode once quasilinear stabilisation of the ion mode has eaten up the ion anisotropy and the conditions for trapping and pair formation have emerged. (In addition pair anisotropy will also destabilise the electron mirror mode in the stable quasilinear state of the ion mirror mode.) In both stable quasilinear cases the remaining ion (electron) anisotropy is negligibly small, and the quasilinear growth rate of the ion mirror mode is zero by deﬁnition. Thus, starting from the quasilinear level one must use the pair anisotropy in Eq. (1) ki e⊥ γ (k) ≈ k V A (62) m,pair k A,ql e π T i⊥ including the expression for the electron anisotropy A . This applies to the quasilinear level where the depletion of the magnetic ﬁeld is remarkable though not large, and the Alfvén speed V is based on the weak quasilinear magnetic ﬁeld, while β A,ql ki is on the quasi-linearly heated plasma level. Since, quasi-linearly, T ≈ T has decreased to approach T , the temperature i⊥ ik ik ratio T /T ≈ T /T increases. However, one has to determine A as function of the pairs. Since A is a pressure ratio, e⊥ i⊥ e⊥ ik e e one must take into account also the non-paired isotropic electron component. Let the fraction of pairs be α, then one has approximately (1− α)T + 2αT 2α e e A ≈ − 1 ≈ (63) (1− α)T + αm |u| 1− α e e 20 which gives for the electron-pair generated ion-mirror growth rate 2α 1 + α β e,ql γ (k) ≈ k V (64) m,pair A,ql 1− α π where β = 2µ NT /B is the electron-β based on B at quasilinear stability, and the root 1 + α arises from the com- e,ql 0 e ql ql bination of electron and pair pressures. Since α < 1, the pair induced electron anisotropy is not large but will be sufﬁcient to break the quasilinear state and cause further growth of the ion mirror mode. This further growth violates pressure equilibrium and, if not immediately restored by breaking the pairs off and heating the plasma to restore isotropy, pressure equilibrium must be restored otherwise. Under closed boundary conditions stabilisation proceeds again quasi-linearly through heating the plasma on the expense of the pair anisotropy A and, in addition, radiation of resonant whistlers. On the other hand, if open boundary condition allow inﬂow of cold charge-neutralised plasma from the environment, pressure balance will achieved in this way selecting particles of low perpendicular energy. Mirror modes will in both cases quickly restore pressure balance while growing to substantially larger than quasilinear amplitudes because at any quasilinear level new electron pairs will be newly produced. The magnetic energy density W = B /2µ will then about exponentially decrease n o W (t− t ) = W 1− exp 2γ t− t (65) ql ql m,pair ql from its quasilinear level W (t ) until reaching its ﬁnal minimum value W (t ) at ﬁnal time t which is in equilibrium ql ql min ﬁn ﬁn with the environmental pressure P = W + β according to ext ext ⊥,ext W 1 + β min ⊥ext = (66) W 1 + β ext ⊥ﬁn Combining the last two equations gives an estimate for the time t when the mirror mode reaches its ﬁnal equilibrium state ﬁn W 1 + β min ⊥ext 2γ (t − t ) = log 1− (67) m,pair ﬁn ql W 1 + β ql ⊥ﬁn With t ≫ t and W ≪ W this expression simpliﬁes to become ﬁn ql min ql W 1 + β min ⊥ext γ t ≈ (68) m,pair ﬁn 2W 1 + β ql ⊥ﬁn Moreover, β = (1 + α)β W /W ≫ 1. One then obtains the following estimate for the ﬁnal saturation time ⊥ﬁn ql ql min √ √ 2 2 t π 1− α 1 + β W π W ﬁn ext min min ≈ ≈ (69) 3/2 3/2 τ 2α W 2αβ W (1 + α) A ql ql ql ql −1 where τ = k V is the parallel Alfvén time based on the quasilinear saturation ﬁeld B . In a mirror bubble of a few A k Aql ql 3 −3 10 km parallel length, density a few 10 cm and some 10 nT internal ﬁeld the Alfvén time is of the order of τ ∼ 1 min, sufﬁciently long. Hence, if the ﬁnal saturation time is about one Alfvén time we have for an estimate of the fraction of pairs in the mirror bubble π W min −4 α ≈ ∼ O(10 ) (70) 2β W ql ql 21 contributing both to the growth beyond quasilinear saturation limit and maintenance of pressure balance with the external plasma and ﬁeld, with measurable right-hand side. This fraction is sufﬁciently small. Not more than roughly ≈ 10 pairs per m sufﬁce for producing the desired effect for instance in the magnetosheath, a rather small, unfortunately probably barely detectable fraction. Again similar to superconductivity, it is only their effect which indicates their presence. Here we have not included any highly probable inﬂow of plasma along the magnetic ﬁeld to help restore the required pressure balance. Such an inﬂow is expected to occur because the different mirror bubbles in the magnetosheath evolve on different time scales. Being magnetically connected they compete and tend to exchange low energy plasma along the magnetic ﬁeld on their mutual expenses. An analogous mechanism will work for the electron mirror mode which, on the quasilinear level, will grow fast and attract cold electrons from the environment inside the ion mirror mode. In addition, the perpendicular pressure anisotropy of the pairs may in both cases, the ion as well as the electron mirror mode, also excite resonant whistlers as well as Bernstein modes, the latter propagating in the perpendicular direction and having a characteristic banded structure following the electron cyclotron harmonics. Observation of these wave spectra, in particular Bernstein modes, should provide direct information on the presence of electron pairs. It is clear that ion pairs, if formed, would behave similarly in the ion mirror mode. We have, however, not checked this possibility here. 5.4 Magnetic susceptibility Above we developed a dynamical physical evolution model for the mirror mode which allows it to grow beyond the quasilinear limit. The ﬁnal state is an equilibrium which should be treated in the framework of thermodynamics or statistical mechanics of open systems. This does not mean open boundary conditions. It means that the system is embedded into a larger system with that it exchanges information and possibly energy and for that the distortion it causes is negligibly small, say a mirror mode train in the magnetosheath. In this ultimately achieved and observed or measured equilibrium state the mirror mode should be described by thermodynamic quantities, and the localised expulsion of the magnetic ﬁeld should be accounted for by a diamagnetic susceptibility χ. This susceptibility ∂M pair χ = µ (71) pair 0 ∂B is deﬁned as the derivative of the total pair-magnetic moment with respect to the magnetic ﬁeld, with the former not known and very difﬁcult to determine. It requires knowledge of the grand partition function Q of which it is the logarithmic derivative ∂ logQ M = T N (72) pair e pair ∂B N,T,µ with µ chemical potential. Thus the susceptibility is the second derivative of the logarithm of the grand partition function with respect to the magnetic ﬁeld. This requires knowledge of all energy states of the pairs in the volume and in addition their spatial distribution. The former can be restricted to only one state, the perpendicular ﬁnal temperature T , while the spatial distribution requires model assumptions on the geometry of the ultimate mirror bubble and the diamagnetic surface current. Both barely 22 known. Thus we are basically unable to solve 0this basic physics problem whose solution would provide insight into the most important thermodynamical connections. We can, however, referring to a heuristic approach, determine the sign of the susceptibility by assuming that to some extent the total magnetic moment of all pairs that are localised in the mirror bubble is proportional to the sum of all single moments of the pairs, while the dependence on the magnetic ﬁeld is reasonably unchanged. This may not be completely true though, because in analogy to superconductivity and superﬂuidity, the mirror mode diamagnetism is a second order phase transition, and therefore the dependence on B must be modiﬁed to some power B , with δ 6= 1 which cannot be an integer number. Neglecting this complication, we can write hE i pair M = αN (73) pair where essentially hE i ∼ T is the electron temperature, and as above α accounts for the fraction of pairs in the volume. pair e Then we obtain that hE i pair χ = −µ αN (74) pair 0 The important result is that this magnetic susceptibility is negative which is required for observing a local diamagnetic effect in excess of the general weak about unremarkable bulk diamagnetism of any plasmas which is caused by the orbital gyration of the charged particles in the magnetic ﬁeld and is of the order of Landau diamagnetism (Huang, 1987). Though the precise dependence on the magnetic ﬁeld is not known, this argument suggests that the process leading to pair formation in mirror modes causes the diamagnetism required for blowing up the mirror bubble, expelling the magnetic ﬁeld and causing chains of mirror modes. Of course, this estimate suffers from the involved uncertainties and the impossibility of constructing the condi- tions of phase transition in the high temperature mirror plasma which also accounts for pressure balance with the surrounding plasma and determines the ﬁnal magnetic amplitude B in saturated thermodynamical equilibrium from the general relation ﬁn B = B 1 + χ . This is an implicit third order expression for the ﬁnal magnetic ﬁeld ﬁn ql pair 3 β ql x ≈ x− α (75) where x = W /W and β = NT /W . Neglecting the left-hand side yields approximately ﬁn ql ql e ql W β ﬁn ql ≈ α (76) W 2 ql for a very rough estimate. From the previous subsection we conclude that the left-hand side is of same order as α, say. Thus β /2 ∼ O(1) which, in spite of the rough assumptions made, is not an unreasonable value for the quasilinear equilibrium ql where pair formation sets on. However, the main important result is the negative sign of the susceptibility which. independent on the real numerical value of χ, suggests that the mirror mode is not a simple plasma instability as such. Rather it is a particular plasma state occurring in high temperature plasmas. This kind of diamagnetism resembles a phase transition whose precise physics has still to be developed. Our discovery of the possibility of the involvement of pairing in this case is an interesting though but a ﬁrst step into that direction. 23 6 Conclusions Mirror modes seem to be an exception in high-temperature collisionless plasmas. They start from a simple magnetohydrody- namic instability in an anisotropic pressure conﬁguration far from thermal equilibrium that has been produced for instance in the magnetosheath by the forced ﬂow across the bow shock (for a comprehensive review cf., e.g., Tsurutani et al., 2011; Balogh & Treumann, 2013, for the relations around the bow shock) and may be a general property of shocked plasma ﬂows (Lucek et al., 2005). Linear theory shows that this instability produces magnetic ﬁeld-elongated magnetic bottles which sta- bilise by quasilinear interaction between the anisotropic ions and the magnetic ﬁeld in which course the thermal anisotropy is depleted and settles at a low stable rudimentary value. The amplitude of the magnetic depletion, as numerical simulations with periodic boundary conditions demonstrate, is very low. It is in fact so low that the quasilinear mirror mode would in ob- servations not be noticed but added to the ordinary thermal ﬂuctuations of the magnetic ﬁeld and thermal pressure. It does not explain the notorious though not persistent observation of very large amplitude chains of mirror modes of up to 50% magnetic depletion. Open boundary simulations (Shoji et al., 2012) show the evolution of large amplitude mirror modes with electron dynamics reduced to a neutralising ﬂuid. Dynamics is purely ionic then, and pressure balance is provided by forced inﬂow of plasma from the external surrounding. This suggests that ions are capable of generating large amplitude structures, and that pressure balance is achieved by external plasma inﬂow. However, the assumption of ﬂuid electrons is strong in this case and may not apply to real natural plasmas where electrons are far from being a ﬂuid. For this to happen ion anisotropies must be large enough for pressure deﬁcit to enforce inﬂow prior to quasilinear stabilisation. As noted above ion pairing would be a reason to generate such a deﬁcit. In the present communication we have not considered ion pairing except for pointing on its possibility. Instead we demon- strated that the physics of large amplitude mirror modes can be affected by a pairing mechanism which is unique in application to mirror modes as these provide the rare conditions where it may happen. We demonstrated this on the example of electrons but the theory can indeed be extended to include ions as well. A possible resolution of the large mirror amplitude problem is thus found when accounting for the dynamics of the electron (or possibly ion) component trapped in the magnetic mirror bottle. Electrons perform their bounce motion and can, in the vicinity of their mirror points, where the parallel speed drops to near zero, get into resonance with the always present thermal ion-acoustic noise spectrum. Experiencing a modiﬁed dielectric constant they generate an attractive potential difference in their wake outside the charge-compensating Debye sphere that may affect another close electron, attract it and form an electron singlet pair which consists of the primary and attracted electrons in interaction with the resonant spectrum of ion sound background waves. Formation of triplet pairs are of higher order and thus substantially less probable. The conditions for this to happen have been obtained. The electron pairs are not spin compensated because at the high temperatures in classical plasmas the Pauli principle plays no role. The remaining electron jitter energies are still high above spin energies but may be in many cases too low for letting the pair, which near its common mirror point Z(s ) is at average velocity U(z ) = c ≪ v , return into bounce m m s e 24 motion as this would require that the jitter energy exceeds the attractive pair potential. We have not solved the stability problem of this process as this would become a separate investigation. However, a group of pairs fulﬁls an important function in mirror modes. Trapped in attraction, they drop out of the bounce motion, become locked near mirror points along the quasi-linearly stable mirror bottle, and spend all their kinetic energy into their gyration. Thus the pair distribution in mirror modes becomes highly and narrowly peaked just above and near the perpendicular thermal velocity, an effect which is very interesting to investigate in all its further consequences. These may be manifold, since the pairs contribute a highly anisotropic population which, as noted, may become unstable to electrostatic and electromagnetic plasma modes. Among those are both whistlers, and electrostatic Bernstein modes. Whistlers are actually known to exist in observations and may be related to electron pairing. In the case of ion paring one would expect similar effects, in particular kinetic Alfvén waves. Closer investigation has been given to two effect, the direct production of diamagnetism via the magnetic susceptibility (see the last subsection), and the contribution of the temperature anisotropy of the pairs to the envisaged further evolution of ion-mirror modes. It has turned out that this effect may be realistic. It is based on two observations, the dependence of the ion mirror growth rate on the electron anisotropy which under normal conditions would not be of any interest as it provides just a small negligible electron contribution. However, in the case of quasilinear stabilisation the ion anisotropy is depleted and the growth rate vanishes. It is just this case when electron (and possibly ion pairs) are produced and come into play, as we have demonstrate above. (Considering the ion effect which we did not do leads to the result that ion pairs would restore an ion anisotropy and in this way break the quasilinear state and cause further growth of the mirror mode. This may have happened in the simulations of Shoji et al. (2012) where it neither has been discussed nor taken into account and where just large mirror modes have been demonstrated to develop in open boundary simulations.) We suggested that electron pairs even in ion mirrors produce such an effect. This seems reasonable because of the high vulnerability of the mobile electrons to interact strongly with a background ﬁeld. The production of even a small fraction of electron pairs breaks the quasilinear stability condition and causes further growth of the ion mirror instability. Pressure balance is restored as we have shown even for a small fraction of electron pairs per unit volume being created. It can also be warranted by forced inﬂow of cold plasma from the environment along the magnetic ﬁeld, which probably happened in the cited simulations, and effect we are aware of but have not included in any detail here. Altogether, the present communication discovered an interesting new effect in high-temperature plasma which might have other consequences as well. It brings the theory of mirror modes to an intermediate physical conclusion by contributing to the so far badly understood generation of large magnetic amplitudes of mirror bubbles, the deep diamagnetic holes in the magnetic ﬁeld. It also provides an important unexpected, at least interesting application of the apparently superﬁcial attractive electron potentials in a plasma. Here we have demonstrated its possible importance in the evolution of mirror modes when electron-pair singlets form in close analogy to superconductivity. It would be of interest investigating which effects the process elucidated here might have in turbulence theory and as well in astrophysical applications, in particular in view of our above ﬁnding that kinetic Alfvén waves would also be capable of 25 generating attractive potentials to form electron pairs. Since they naturally have high phase speeds, satisfaction of the resonance condition seems to be natural for them with bouncing electrons for instance in the auroral magnetosphere. An example where kinetic Alfvén waves could become involved is auroral-plasma sheet coupling where electrons are naturally trapped in the geomagnetic ﬁeld and perform large-scale bounce motions. The relaxed kinetic Alfvén pairing condition may generate a fraction of trapped pairs along the auroral magnetic ﬁeld in this case to produce observable effects for instance in the aurora, generation of radiation, and in reconnection. In mirror-modes kinetic Alfvén waves are of little interest, as there is no obvious reason for them to be generated. They have, moreover, never been identiﬁed in relation to mirror observations while ion sound waves are generally present within and outside them. In ion-inertial range turbulence kinetic Alfvén waves seem to play some role as various observations indicate and theory also supports for the reason that the scale of the ion-inertial range coincides with the perpendicular wave lengths of kinetic Alfvén waves. Our pair-singlet mechanism should work in those cases as well and might have consequences for turbulence, entropy generation, and turbulent dissipation. Acknowledgement. This work was part of a brief Visiting Scientist Programme at the International Space Science Institute Bern. We ac- knowledge the interest of the ISSI directorate as well as the generous hospitality of the ISSI staff, in particular the assistance of the librarians Andrea Fischer and Irmela Schweitzer, and the Systems Administrator Saliba F. Saliba. 26 References Ahmadi N, Wilder FD, Ergun RE, Argall M, Usanova ME, Breuillard H, Malaspina D, Paulson K, Germaschewski K, Eriksson S, Goodrich K, Torbert R, Le Contel O, Strangeway RJ, Russell CT, Burch J, and Giles B: Generation of electron whistler waves at the mirror mode magnetic holes: MMS observations and pic simulations, accessible at arXiv:1803.06399v1 [physics.space-ph], J. Geophys. Res., 123, 6383-6393, https://doi.org/10.1029/2018JA025452, 2018. Balogh A and Treumann RA: Physics of Collisionless Shocks (Springer: New York) pp. 500, https://doi.org/10.1007/978-1-4614-6099-2, Bardeen J, Cooper LN, and Schrieffer JR: Theory of superconductivity, Phys. Rev. 108, 1175-1204, https://doi.org/10.1103/PhysRev.108.1175, 1957. Baumjohann W and Treumann RA: Basic Space Plasma Physics, London 1996, Revised and enlarged edition, Imperial College Press, London, https://doi.org/10.1142/P850, 2012. Binney JJ, Dowrick NJ, Fisher AJ, and Newman MEJ: The Theory of Critical Phenomena (Oxford University Press: Oxford, 1999. Breuillard H, Le Contel O, Chust T, Berthomier M, Retino A, Turner DL, Nakamura R, Baumjohannn W, Cozzani G, Catapano F, Alexan- drova A, Mirioni L, Graham DB, Argall MR, Fischer D, Wilder FD, Gershman DJ, Varsani A, Lindqvist P-A, Khotyaintsev YuV, Marklund G, Ergun RE, Goodrich KA, Ahmadi N, Burch JL, Torbert RB, Needell G, Chutter M, Rau D, Dors I, Russell CT, Magnes W, Strangeway RJ, Bromund KR, Wei H, Plaschke F, Anderson BJ, Le G, Moore TE, Giles BL, Paterson WR, Pollock CJ, Dorelli JC, Avanov LA, Saito A, Lavraud B, Fuselier SA, Mauk BH, Cohen IJ and Fennell JF: The properties of Lion roars and electron dynamics in mirror mode waves observed by the Magnetospheric MultiScale Mission, J. Geophys. Res. 123, 93-103, https://doi.org/10.1002/2017JA024551, 2018. Chandrasekhar S: Hydrodynamic and hydromagnetic stability, Clarendon Press: Oxford UK, 1961. Constantinescu OD: Self-consistent model of mirror structures, J. Atmosph. Solar-Terr. Phys. 64, 645-649, https://doi.org/10.1016/S1364- 6826(02)00024-X, 2002. Constantinescu OD, Glassmeier K-H, Treumann RA, and Fornaçon K-H: Magnetic mirror structures observed by Cluster in the magne- tosheath, Geophys. Res. Lett. 30, ID 1802, https://doi.org/10.1029/2003GL017313, 2003. Cooper LN: Bound electron pairs in a degenerate Fermi gas, Phys. Rev. 104, 1189-1190, https://doi.org/10.1103/PhysRev.104.1189, 1956. Gary SP: Theory of space plasma microinstabilities, Cambridge Univ. Press, Cambridge UK, 1993. Ginzburg VL and Landau LD: On the theory of superconductivity, J. Exp. Theoret. Phys. 20, 1064-, 1950. Ginzburg VL: On the theory of superconductivity, Nuovo Cim. 2, 1234-1250, https://doi.org/10.1007/BF02731579, 1955. Gradshteyn IS and Ryzhik IM: Tables of Integrals, Series, and Products (Academic Press: New York) pp. 1089, 1965. Hasegawa A: Drift mirror instability in the magnetosphere, Phys. Fluids 12, 2642-2650, https://doi.org/10.1063/1.1692407, 1969. Hasegawa A: Plasma instabilities and nonlinear effects (Springer-Verlag: Berlin-New York, 1975. Huang K: Statistical Mechanics. Second Edition (John Wiley & Sons: New York, 1987. Krall NA and Trivelpiece AW: Principles of Plasma Physics (McGraw-Hill: New York, 1973. Landau L: Theory of the superﬂuidity on Helium II, Phys. Rev. 60, 356-358, https://doi.org/10.1103/PhysRev.60.356, 1941. Lucek EA, Constantinescu D, Goldstein ML, Pickett J, Pinçon JL, Sahraoui F, Treumann RA, and Walker SN: The magnetosheath, Space Sci. Rev. 118, 95-152, https://doi.org/10.1007/s11214-005-3825-2, 2005. 27 Lucek EA, Dunlop MW, Horbury TS, Balogh A, Brown P, Cargill P, Carr C, Fornaçon KH, Georgescu E, and Oddy T.: Cluster magnetic ﬁeld observations in the magnetosheath: four-point measurements of mirror structures, Ann. Geophysicae 19, 1421-1428, https://doi.org/10.5194/angeo-19-1421-2001, 2001. Nambu M and Akama H: Attractive potential between resonant electrons, Phys. Fluids 28, 2300-2301, https://doi.org/10.1063/1.865284, Neufeld J and Ritchie RH: Passage of charged particles through plasma, Phys. Rev. 98, 1632-1642, https://doi.org/10.1103/PhysRev.98.1632, Noreen N, Yoon PH, López RA, and Zaheer S: Electron contribution in mirror instabiity in quasi-linear regime, J. Geophys. Res. A (Space Phys.) 122, 6978-6990, https://doi.org/10.1002/2017JA024248, 2017. Rodriguez P and Gurnett DA: Electrostatic and electromagnetic turbulence associated with the earth’s bow shock, J. Geophys. Res. 80, 19-31-386, https://doi.org/10.1029/JA080i001p00019, 1975. Shoji M, Omura Y, and Lee LC: Multidimensional nonlinear mirror-mod structures in the Earth’s magnetosheath, J. Geophys. Res. 117, A08208, https://doi.org/10.1029/2011JA017420, 2012. Treumann RA and Baumjohann W: Advanced Space Plasma Physics, London 1997, Imperial College Press, London, https://doi.org/10.1142/p020, 1997. Treumann RA and Baumjohann W: Plasma wave mediated attractive potentials: a prerequisite for electron compound formation, Ann. Geophys. 32, 975-989, https://doi.org/10.5194/angeo-32-975-2014, 2014. Treumann RA and Baumjohann W: The mirror mode: A "superconducting" space plasma analogue, Ann. Geophys. 36, 1015-1026, https://doi.org/10.5194/angeo-36-1015-2018, 2018. Treumann RA and Baumjohann W: Electron mirror branch: observational evidence from “historical” AMPTE-IRM and Equator-S measure- ments, Ann. Geophys. 36, 1563-1576, https://doi.org/10.5194/angeo-36-1563-2018, 2018. Tsurutani BT, Lakhina GS, Verkhoglyadova OP, Echer E, Guarnieri FL, Narita Y, and Constantinescu OD: Magnetosheath and heliosheath mirror mode structures, interplanetary magnetic decreases, and linear magnetic decreases: Differences and distinguishing features, J. Geophys. Res. Space Phys. 116, ID A02103, https://doi.org/10.1029/2010JA015913, 2011. Yao ST, Shi QQ, Yao ZH, Li JX, Yue C, Tao X, Degeling AW, Zong QG, Wang XG, Tian AM, Russell CT, Zhou XZ, Guo RL, Rae IJ, Fu HS, Zhang H, Li L, Le Contel O, Torbert RB, Ergun RE, Lindqvist PA, Pollock CJ, and Giles BL: Waves in Kinetic-Scale Magnetic Dips: MMS Observations in the Magnetosheath, Geophys. Res. Lett. 46, 523-533, https://doi.org/10.1029/2018GL080696, 2019. Yoon PH: Kinetic theory of hydromagnetic turbulence. I. Formal results for parallel propagation, Phys. Plasmas 14, 102302, doi: 10.1063/1.2780139, 2007.

Astrophysics – arXiv (Cornell University)

**Published: ** Sep 23, 2019

Loading...

You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!

Read and print from thousands of top scholarly journals.

System error. Please try again!

Already have an account? Log in

Bookmark this article. You can see your Bookmarks on your DeepDyve Library.

To save an article, **log in** first, or **sign up** for a DeepDyve account if you don’t already have one.

Copy and paste the desired citation format or use the link below to download a file formatted for EndNote

Access the full text.

Sign up today, get DeepDyve free for 14 days.

All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.