Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

TMP 91 E

Download as pdf or txt
Download as pdf or txt
You are on page 1of 15

Ann. Geophys.

, 32, 975989, 2014


www.ann-geophys.net/32/975/2014/
doi:10.5194/angeo-32-975-2014
Author(s) 2014. CC Attribution 3.0 License.
Plasma wave mediated attractive potentials: a prerequisite for
electron compound formation
R. A. Treumann
1,2,*
and W. Baumjohann
3
1
Department of Geophysics and Environmental Sciences, Munich University, Munich, Germany
2
Department of Physics and Astronomy, Dartmouth College, Hanover, NH 03755, USA
3
Space Research Institute, Austrian Academy of Sciences, Graz, Austria
*
currently at: International Space Science Institute, Bern, Switzerland
Correspondence to: R. A. Treumann (rudolf.treumann@geophysik.uni-muenchen.de)
Received: 26 February 2014 Revised: 16 July 2014 Accepted: 17 July 2014 Published: 22 August 2014
Abstract. Coagulation of electrons to form macro-electrons
or compounds in high temperature plasma is not generally
expected to occur. Here we investigate, based on earlier
work, the possibility for such electron compound formation
(non-quantum pairing) mediated in the presence of vari-
ous kinds of plasma waves via the generation of attractive
electrostatic potentials, the necessary condition for coagula-
tion. We conrm the possibility of production of attractive
potential forces in ion- and electron-acoustic waves, pointing
out the importance of the former and expected consequences.
While electron-acoustic waves presumably do not play any
role, ion-acoustic waves may potentially contribute to for-
mation of heavy electron compounds. Lower-hybrid waves
also mediate compound formation but under different con-
ditions. Buneman modes which evolve from strong currents
may also potentially cause non-quantum pairing among
cavity-/hole-trapped electrons constituting a heavy electron
component that populates electron holes. The number densi-
ties are, however, expected to be very small and thus not vi-
able for justication of macro-particles. All these processes
are found to potentially generate cold compound popula-
tions. If such electron compounds are produced by the at-
tractive forces, the forces provide a mechanism of cooling
a small group of resonant electrons, loosely spoken, corre-
sponding to classical condensation.
Keywords. Space plasma physics (magnetic reconnection;
waveparticle interactions; waves and instabilities)
1 Introduction
Plasmas consist of equal numbers of electrons and ions
forming quasi-neutral uid-like matter at temperatures suf-
ciently high for maintaining ionisation and with very large
numbers of electrons populating the Debye sphere. At such
high temperatures electrons and ions are mutually well sepa-
rated located at instantaneous distances L
N
N
1/3
, where
N N
e
=N
i
is the average plasma density in a singly
charged plasma. The electrostatic Coulomb elds of the
naked electric charges are conned to Debye spheres of ra-
dius
D
by the collective effect of the many particles of oppo-
site charge passing around at their average tangential speeds
within radial distances r
D
. The geometric shape of the
Debye spheres is very close to a sphere, deviating from it
only in very strong magnetic elds and for very high plasma
ow speeds. Outside the Debye sphere the residual particle
eld decays exponentially while contributing to a thermal
uctuation background eld. From a particle point of view
each of the plasma particles is a charged Fermion. In classi-
cal plasmas at the high plasma temperatures the spin has no
importance, and the fermionic property of the particles plays
no role. In quantum plasmas, which for obeying quantum
properties must be dense, this property is rather important.
For, when two electrons form pairs, the spins add up and the
pair becomes a Boson of either zero or integer spin. Many
pairs can occupy the same energy level and, altogether, tend
to condensate in the lowest energy level permitted by suf-
ciently low temperatures (cf., however, the Appendix). This
property is very well known from solid state physics (cf. e.g.
Fetter and Walecka, 1971; Huang, 1987).
Published by Copernicus Publications on behalf of the European Geosciences Union.
976 R. A. Treumann and W. Baumjohann: Electron coagulation
Pairing is not expected under any normal plasma con-
ditions. However, when a plasma wave passes across the
non-quantum plasma, the dielectric properties of the plasma
change. Electrons assuming a relative velocity with respect
to the phase velocity of the plasma wave nding themselves
exposed to the dielectric polarisation which adds to the De-
bye screening that compensates for the naked particle charge.
Such electrons may evolve an attractive electrostatic inter-
particle force acting on its neighbour electrons, an effect
different from classical wave trapping in the wave poten-
tial trough that causes wave saturation and other nonlinear
effects like solitons and holes. The attractive forces are di-
rect current (dc) forces. Experienced by two electrons of ap-
proximately same velocity, they bind these together to form
binary compounds, that is, classical pairs. Experienced
by many electrons of same velocity, they may form large
electron compounds the nature of which is that of massive
macro-particles of same charge-to-mass ratio e / m being
sufciently strongly correlated to behave dynamically like
one single particle. Below we demonstrate that it is not the
electrons themselves which coagulate but their Debye clouds
such that the repulsive forces between the coagulating elec-
trons do grow only very weakly with the number n
com
of
coagulations.
Depending on the number n
com
of particles in the newly
formed conglomerate being odd or even, macro-electrons
in principle behave like Fermions or Bosons. At very low
temperatures Bosons are known to condensate to form a
dense population. At the high plasma temperatures, triv-
ially implying thermal wavelengths
T
=
_
2

h
2
/m
e
T
10
9

1 eV/T m of the order of atomic scales, such effects


are inhibited (for another strong counter argument see the
Appendix). Classical coagulation of particles caused by at-
tractive inter-particle forces is, however, permitted.
In the following we show that, under certain conditions
in high temperature plasmas, attractive potentials between
neighbouring electrons are indeed produced as is schemat-
ically shown in Fig. 1. They are a necessary condition for
subsequent coagulation of electrons. The intention is to nd
out whether it provides a natural mechanism justifying the
still unavoidable assumption of macro-particles in numerical
particle-in-cell simulations.
The nding is that, though electrons may possibly coag-
ulate, their numbers will in all cases remain very small be-
ing just of the order of a fraction of the resonant particles
which are a small fraction of the plasma population only, too
small to be taken as justication of the assumption of macro-
particles.
The conditions for generation of attractive electrostatic po-
tentials in the presence of plasma waves are derived as a nec-
essary condition for coagulation. The numbers of particles in
each of the electron compounds, the compound densities and
temperatures, require studying the sufcient conditions for
coagulation for each kind of waves, which requires knowl-
edge of the wave spectrum.
1
2 Generation of wave-mediated attractive potentials
The method of calculating the potential around a test charge
in plasma was explicated sixty years ago (Neufeld and
Ritchie, 1955). Thirty years later it was revived (Nambu and
Akama, 1985) to include the effect of plasma waves and was
used in this form to suggest the coagulation of particles in
the presence of dust in plasmas in order to explain the for-
mation of dust structure, which for some while became an
industry (cf. e.g. Shukla and Melands, 1997; Shukla et al.,
2001; Nambu and Nitta, 2001, and references therein).
Following Neufeld and Ritchie (1955), the general ex-
pression for the electrostatic potential (x, t ) is obtained
from Poissons law
2
(x, t )(x, t ) =qN(x, t )/4
0
.
The density of a test particle traversing the plasma with ve-
locity v and charge q =q
t
is N
t
=(2)
3
(r). In Fourier
space this yields for the potential (cf. Neufeld and Ritchie,
1955; Sitenko, 1967; Krall and Trivelpiece, 1973; the latter
two for textbook presentations)
(x, t ) =
q
t
8
2

0
_
dkd
(k v)
k
2
(k, )
e
ikr
(1)
Here r =x vt is the distance between the location x

=vt
of the particle and the reference point of measurement of the
potential disturbance. (k) is the frequency of a spectrum of
plasma wave eigenmodes, presumably present in the plasma
as background noise or wave excitations, with wave number
k. The function (k, ) is the complete dielectric plasma re-
sponse function
2
experienced by the particle including the
disturbance caused by the moving charge of the test particle.
(One may note that the function in the numerator can be
used to replace the wave-number component parallel to the
particle velocity in the Fourier exponential exp(ik r) reduc-
ing it to an integration.)
The problem as seen from the test charge is spherically
symmetrical. Thus it makes sense to formulate it in spherical
coordinates r, k,
r
,
k
both in wave number and real space,
with conventional notation for the angular volume elements.
Choosing an expansion into spherical harmonics (as done in
1
We thank the referees for reminding us that generation of at-
tractive potentials is merely the necessary condition for coagula-
tion. Real coagulation will take place only if a number of sufcient
conditions is also met, a point that will be further discussed in the
Discussion and Conclusions section.
2
It may be important to remark that in Poissons law no assump-
tion is made about the linearity or nonlinearity of (x, t ) which,
hence, in principle is the general response function including all
electrostatic nonlinearities, which contribute to the dielectric prop-
erties of the plasma. Thus in full generality one should rather write
[k, , (k)]. In discussing the Buneman mode later, we will in
passing make use of this more general form.
Ann. Geophys., 32, 975989, 2014 www.ann-geophys.net/32/975/2014/
R. A. Treumann and W. Baumjohann: Electron coagulation 977
Neufeld and Ritchie, 1955) one has for the exponential factor
e
ikr
=4

m
i
l
_

2kr
J
l+
1
2
(kr)Y
m
l
(
k
)Y
m
l
(
r
). (2)
The notation for the spherical harmonics Y
m
l
(
k
), Y
m
l
(
r
)
is again conventional in wave number and real space, and the
indicates the conjugate complex version of the azimuthal
exponentials exp(im
k
). The k, r dependence is taken care
of by the half integer Bessel functions J
l+1/2
(kr).
The response function [k, (k)] is a function of fre-
quency and wave number and is taken in the electrostatic
limit. In the above representation it is a scalar function. When
electromagnetic contributions or an external magnetic eld
would have to be taken into account, it becomes a tensor.
Only its longitudinal part
L
=k k/k
2
enters the expres-
sion for the potential, however. In addition, one would have
to consider a variation of the vector potential caused by the
test particle, if transverse waves are included.
In the absence of the latter, it is well known that the po-
tential of the test particle, in our case an electron of el-
ementary charge q =e, consists of its Coulomb poten-
tial
C
(r) =q/4
0
r, which will be Debye screened by
the plasma particles becoming
D
(r) =
C
(r)exp(r/
D
),
and a disturbance caused by the reaction of the plasma eigen-
modes to the presence of the charge the eigenmodes that
are either present or are amplied by the moving test charge
for which the plasma appears as a dielectric to whose normal
modes the charge couples. These eigenmodes contribute via
adding each of them to the vacuum dielectric constant its par-
ticular susceptibilities
s
(, k), where s is the index identi-
fying the particle species which responds to the eigenmodes.
Hence, with s =e, i for electrons and ions, respectively:
(k, ) =1 +(k
D
)
2
+
e
(k, ) +
i
(k, ) (3)
Independent of the wave modes, the Debye term on the right
is included here in order to account for the presence of the
point like test charge. In a non-magnetised plasma the sus-
ceptibilities assume the form

s
(k, ) =(k
Ds
)
2
[1 +
s
Z(
s
)] , (4)

s
=(k u
s
)/kv
s
.

Ds
is the Debye length of species s, Z(
s
) the plasma dis-
persion function, v
s
the thermal speed of species s, and u
s
is a possible bulk streaming velocity of species s which,
for our application, will for simplicity be put to zero but
has to be retained both for streaming and electric currents
J =

s
q
s
N
s
u
s
.
One should note that Poissons law is quite general hold-
ing for both linear and nonlinear interactions. Restriction
to linear response functions only implies small disturbances
caused. For large nonlinear disturbances the linear response
function would have to be replaced by its nonlinear counter-
parts.
C
D
Distance
P
o
t
e
n
t
i
a
l

De

Figure 1. Schematic comparison between the dc Coulomb (


C
),
Debye (
D
) and plasma wave modied (
a
) potentials for the case
when the modied potential evolves an attractive domain. In all
cases considered the potential is repulsive inside the Debye sphere
at distances
De
. In the presence of (electrostatic) plasma waves
investigated in this paper the dc potential
a
outside but close to a
distance
De
develops a region where it becomes negative. At large
distances the effect of the plasma wave decreases rapidly.
These expressions have been partially analyzed in the
available literature with focus on the effect of dust in plas-
mas (adding the susceptibility of dust particles in Shukla and
Melands, 1997; Shukla et al., 2001, and others). In the fol-
lowing we follow some of the lines in these papers in view
of application to space plasma conditions and with the inten-
tion of checking the chances for classical electron pairing
and possible related effects like formation of macro-electron
compounds. Being aware that there is no quantum pairing,
we nevertheless keep the expression pairing in the follow-
ing for the reason that we consider attractive forces between
two electrons. This is trivially generalised later to account for
formation of coagulations of electrons.
3
3
An apparently strong argument against the effectiveness of at-
tractive potentials is based on the comparison of the strengths of the
attractive and the wave potentials. Attractive potentials are weak
and thus apparently negligible, in particular when the wave spec-
trum is unstably excited. This argument fails for two reasons. First,
at the location of the attractive potential the wave potential and the
resulting electric force act both on all electrons in the same direc-
tion; in contrast the attractive potential acts between electrons, com-
pletely independent of the direction of the wave electric eld. It can-
not be compensated by the wave eld. Second, the scale of the wave
electric eld is given by k
De
1, that is, it is large against the
Debye length (Fig. 1). In contrast, the attractive potentials between
electrons, as shown below, act on scales
De
, only slightly larger
than the Debye length, a scale on which the wave potential is practi-
cally constant for all waves under consideration, and no force exists
which could compensate for the attractive force even then, when the
wave is strongly excited, reaches very high amplitudes and evolves
nonlinearly. Its nonlinear evolution being taken care already in the
general response function [, k, (k)].
www.ann-geophys.net/32/975/2014/ Ann. Geophys., 32, 975989, 2014
978 R. A. Treumann and W. Baumjohann: Electron coagulation
2.1 Ion-acoustic pairing potential
Our rst example is the response of the test charge poten-
tial to the presence of a spectrum of ion-acoustic waves in
a plasma, a problem originally treated cursorily by Nambu
and Akama (1985). In this case the linear response function
is well known. Neglecting damping, its real part is given by

ia
(, k) =1 +
1
k
2

2
De

2
i

2
_
1 +3k
2

2
Di
_
, (5)
where in the round brackets we iterated the frequency by ap-
proximating it with the ion plasma frequency
i
, which pro-
duced the ion Debye length
Di
. Putting
ia
=0 yields the
ion-acoustic dispersion relation

2
ia
(k) =

2
i
1 +1/k
2

2
De
_
1 +
3T
i
T
e
_
1 +k
2

2
De
_
(6)
+
k
2

2
De
1 +k
2

2
De
N
N
_
The last term in the bracket on the right results from a possi-
ble nonlinear density modulation N. It vanishes in the long-
wavelength regime k
2

2
De
1. In order to proceed, we need
a treatable form of
1
ia
. It is not difcult to show that this can
conveniently be written as
1

ia
(, k)
=
k
2

2
De
1 +k
2

2
De
_
1 +

2
ia
(k)

2
ia
(k)
_
, (7)
which separates it into two parts. The rst term is inde-
pendent of the presence of ion-acoustic waves. It is thus
completely spherically symmetric resulting in the known
Debye-screening potential eld of the point charge (treated
in Neufeld and Ritchie, 1955). Its contribution to the poten-
tial at distances r
De
, large with respect to the Debye ra-
dius, is exponentially small. The nonlinear term in the wave
dispersion relation is of higher order and can be neglected
meaning that the nonlinear modulation of the wave spectrum
is of too large scale for causing a rst-order effect in the po-
tential disturbance.
The dominant effect of the test electron interaction with
the ion-acoustic wave spectrum is contained in the wave-
mediated part. Its main contribution comes from the resonant
denominator in frequency space the contribution of which
dominates over that of the exponentially decreasing screened
repulsing Coulomb potential
D
outside the Debye sphere.
Equation (7) inserted into the general expression for the
test particle potential Eq. (3) yields for the wave-induced
contribution

ia
(x, t ) =
e
2
De
16
2

0
_
dkd
ia
(k)e
ikr
1 +k
2

2
De
(8)

_
(k v)

ia
(k)

(k v)
+
P
ia(k)
_
.
The argument of the function depends on the direction of
electron velocity v =vz, which we arbitrarily chose in z di-
rection. It is then appropriate to treat the integral in cylin-
drical rather than spherical coordinates with wave number
k
z
parallel to the electron velocity v and k

perpendicu-
lar to it. With being the radius in the plane perpendicu-
lar to v, the argument of the exponential becomes ik r =
ik
z
(zvt ) +ik

sin. Referring to the denition of Bessel


functions, the integration with respect to the azimuthal angle
results in the Bessel function of zero order, and the expres-
sion for the potential reads

ia
(z, , t ) =
e
2
De
16
0
_
dk
z
dk
2

J
0
(k

)
ia
(k
z
, k

)
1 +k
2
z

2
De
+k
2

2
De
(9)

_
(k
z
v)

ia
(k)

(k
z
v)
+
ia
(k)
_
e
ik
z
(zvt )
d.

ia
offers a possible change in sign which opens up the pos-
sibility for the potential of becoming attractive for another
electron, in which case two electrons may form pairs.
One rst makes use of the functions to replace k
z
=/v
in the exponential and elsewhere by performing the k
z
inte-
gration. The two singularities at =
ia
require perform-
ing the integration in the complex =()+i() plane
via the principal values of the two integrals and the two
residua with integration contour now closed in the lower-half
plane, that is, for z vt < 0 and damped ion-acoustic waves
() < 0. One readily shows that the principal value van-
ishes, for each integral contributes lim
0
(lnln)+i =
i which cancel when subtracted. The residua yield the res-
onant result

ia
(z, , t ) =
e
2
De
8v
0
_
dk
2

J
0
(k

)
ia
(v, k

)
1 +
2
ia

2
De
/v
2
+k
2

2
De
(10)
sin
_

ia
(v, k

)
_
z
v
t
__
for the waveparticle interaction part of the electrostatic
potential. In this expression the ion-acoustic frequency is
implicitly dened through the ion-acoustic dispersion rela-
tion. Replacing k
z
=
ia
/v, the latter can be iterated, yield-
ing to lowest order in the long-wavelength regime k

De

1 that
2
ia
(v, k

)
2
i

2
De
k
2

[1 +
2
ia
/k
2

v
2
)] c
2
ia
k
2

/[1
(m
e
/m
i
)T
e
/K
test
]. Here K
test
=m
e
v
2
/2 is the test particle
kinetic energy. This becomes simply
2
ia
(v, k

) k
2

c
2
ia
/(1
c
2
ia
/v
2
) with c
ia

De
the ion sound velocity.
The wave-number integral must be truncated at the Debye
radius k

De
1 because inside the Debye sphere the point
charge potential dominates. This accounts for long wave-
lengths only. Then the integral becomes

ia
(z, , t ) C
1
_
0
d
2
J
0
( )
1 +
2
[1 +1/(v
2
/c
2
ia
1)]
sin() (11)
Ann. Geophys., 32, 975989, 2014 www.ann-geophys.net/32/975/2014/
R. A. Treumann and W. Baumjohann: Electron coagulation 979
with =k

/
De
, =/
De
, and =(c
ia
/v)(1
c
2
ia
/v
2
)
1/2
, =(z vt )/
De
. The constant is
C =(e/4
0

De
)(c
ia
/v)(1 c
2
ia
/v
2
)
1/2
. Strictly speaking,
this integral with respect to is the sum of its principal value
and the contribution of the poles at

=i(1 c
2
ia
/v
2
)
1/2
.
At sufciently large particle speeds the pole contribution is
negligible, and only the principal value counts. This is seen
as follows. For resonant particles v c
ia
the poles are purely
imaginary. Extending the singular integral over the entire do-
main implies that only the positive pole contributes, which is
obvious already from Eq. (8) since the exponential vanishes
at large r for positive imaginary part of k

De
= only. In
performing the path integration the pole is surrounded in
negative direction. Taking the residuum yields a term
2i|
+
|
2
J
0
_
i|
+
|
_
sin
_
i|
+
|
_
= (12)
2|
+
|
2
I
0
_
|
+
|
_
sinh
_
|
+
|
_
.
I
0
(x) is the zero-order modied Bessel function. It is obvious
that this entire term for particles close to resonance with v
c
ia
, |
+
| O(v
2
c
2
ia
) is very small, conrming that it can
safely be neglected.
When calculating the principal part of the integral, we con-
sider the case k

< 1, that is, radial distances perpendicu-


lar to the particle velocity less than the ion-acoustic wave-
length but large with respect to the Debye length. Shortest
distances are thus =1, yielding J
0
( ) =J
0
() a func-
tion of the integration variable only, varying in the interval
0.77 < J
0
() < 1.0. Its average is J
0
0.85, which we ex-
tract from the integral

ia
(z, 1, t ) C

1
_
0

2
d sin() (13)
=
C

3
__
1
1
2

2
_
cos + sin 1
_
with C

=CJ
0
. The integral is of the same form as in
Nambu and Akama (1985). The requirement =z vt < 0
implies < 0. The potential becomes negative whenever the
expression in the brackets is positive. The interesting case
is when the test particle moves at velocity v c
ia
exceeding
the wave velocity only slightly. Then || mod 2 > 1, and
the dominant term is
2
cos thus conrming Nambu and
Akama (1985) and yielding

ia
(z, 1, t )
_

_
cos. (14)
The potential is attractive in all regions cos < 0 (i.e. >
/2). In the moving particle frame (z vt )/
De
the po-
tential is attractive behind the particle in its wake in regions
< 0 for v c
ia
. Here the value of the potential is

ia
()
eJ
0

4
0

De
||

cos
_
c
ia

(v
2
c
2
ia
)
1/2
_

, (15)

2

c
ia

_
v
2
c
2
ia

3
2
. (16)
The effective distance |z vt |
De
over which the poten-
tial is attractive is thus given by
att
/
De
(v
2
/c
2
ia
1)
1/2
>
2/, mod 2, that is, the attraction is strongest just outside
the Debye length which implies that two electrons one De-
bye length apart attract each other. In other words, two Debye
spheres mutually overlapping by one Debye radius attract
each other. For the resonant particle velocity this condition
yields v
res
> 1.1c
ia
.
In order to attract another electron, it is clear that the two
electrons must move close to each other within a distance
> 1 in the region of negative
ia
, both being in resonance
with the wave at velocities v c
ia
. In this case they can form
pairs effectively becoming Bosons of either zero or integer
spin. We may note that in a magnetic eld B with the elec-
trons moving along the eld, the ion-sound wave depends
on the propagation angle cos =k B/kB. In this case we
have for the sound speed c
ia
c
ia
cos, and the potential be-
comes a sensitive function of , maximising along B. More-
over, we can set =0 and J
0
(0) =1 as only the distance z
along B comes into play.
This attractive potential has to be compared to the wave
potential
w
the particles are in resonance with. With the as-
sumption k
De
1, we are in the long-wavelength regime
with the potential assumed being nearly constant over the
range of variability of the attractive potential. Thus the attrac-
tive force of the trapping wave potential is small. In negative
wave phases it adds to that of the particle by conning low-
energy electrons in the potential well. These electrons oscil-
late at the high trapping frequency with their average speeds
in resonance with the wave. Wave trapping, though being dif-
ferent in the average, helps attracting as in the attracting po-
tential only the average trapped speed v c
ia
counts. The
high jitter speed at trapping frequency of the electrons aver-
ages out.
Wave-trapped electrons are the best candidates for form-
ing pairs. Moreover, since a pair of charge 2e that has been
formed in the negative wave potential may well by the same
mechanism produce a negative pair potential
pair
=
2
3

ia
over the distance of 3
De
, it may attract other electrons
or pairs to form larger macro-particles of large mass and
charge but constant mass-to-charge ratio. In the extreme
(though possibly unrealistic) case, the maximum number of
coagulated electrons could about equal the number of elec-
trons trapped in the wave potential well, since all the neg-
ative potentials of the particles involved in producing at-
tractive potentials add to the wave potential. In effect this
www.ann-geophys.net/32/975/2014/ Ann. Geophys., 32, 975989, 2014
980 R. A. Treumann and W. Baumjohann: Electron coagulation
mechanism may produce macro-electrons of large mass and
charge which behave like a single particle and have such
properties as exploited in small mass-ratio numerical PIC
simulations.
2.2 Electron-acoustic pairing potential
Another wave of similar dispersion is the electron-acoustic
wave. It is excited wherever the plasma contains two elec-
tron populations of different temperatures and densities. Its
response function resembles that of ion-acoustic waves with
the only difference that two populations of electrons are in-
volved, and ions are assumed forming a xed charge neutral-
ising background such that for the densities N
i
=N
c
+N
h
where indices c, h refer to the cold and hot electron com-
ponents. Electron-acoustic waves are high-frequency waves
in the sense that kv
h
, kv
c
|k u)
c
|, where v
c
, v
h
are
the thermal speeds of the different electron components. The
electron-acoustic dielectric response function in its simplest
form reads

ea
(k, ) =1 +
1
k
2

2
Dh


2
c
(k u
c
)
2
. (17)
The Debye radius for sufciently large temperature differ-
ences T
h
> T
c
is completely determined by the hot compo-
nent, and for xed ions there is no need to include the in
term. u
c
is the bulk streaming velocity of cold electrons. The
inverse of the dielectric function can again been brought into
the same form as for ion-acoustic waves
1

ea
(k, )

k
2

2
h
1 +k
2

2
h
_
1 +

2
ea
(k u
c
)
2

2
ea
_
. (18)
This is exactly the same form as for ion-acoustic waves,
however, now with the electron-acoustic dispersion relation

2
ea
=k
2
c
2
ea
/(1 +k
2

2
h
) and c
2
ea
=v
2
h
(N
c
/N
h
). For this rea-
son, the analysis is the same as for the ion-acoustic wave.
The result has already been given by Shukla and Melands
(1997) and is listed here for completeness only:

ea

_
e/|z ut |
_
cos
_
|z vt |/
h
(1 c
2
ea
/v
2
)
1/2
_
. (19)
The bulk speed of the electrons has been suppressed here.
As before, there are some ranges in which the wave po-
tential at the test charge can be negative and thus attract
other electrons. This will, however, only happen in a plasma
where two widely separated in temperature electron popula-
tions exist of which the colder one is streaming. Interestingly
this might be the case in conditions when BernsteinGreen
Kruskal (BGK) electron hole modes are excited. In this case
the hole generates a dilute hot electron component which is
traversed by a rather cold component of beam electrons. Pos-
sibly in this case mutually attracting electrons become pos-
sible. Unfortunately, electron-acoustic waves have not been
detected in these cases, however, in numerical simulations
of electron hole formation. As electron-acoustic waves re-
quire strong forcing in order to overcome damping, electron-
acoustic waves are not a primary candidate for generating
attractive wave potentials.
2.3 Lower-hybrid pairing potential
A most important medium frequency wave is the lower-
hybrid mode (Huba et al., 1977; Yoon et al., 2002). It prop-
agates in a plasma under almost all conditions on scales be-
low the ion cyclotron radius and frequency. Hence the ions
behave non-magnetically while the electrons are completely
magnetised being tied to the magnetic eld and drifting in the
electric eld of the wave mode. Lower-hybrid waves can be
excited by density gradients, diamagnetic drifts and all kinds
of transverse currents J

s
q
s
N
s
u
s
in a plasma, where
u
s
is the perpendicular drift velocity of species s. They are
primarily electrostatic, propagating at oblique angle with re-
spect to the magnetic eld though being strongly inclined
with k

< k

. Their response function including the test par-


ticle Coulomb potential term reads

lh
(k, ) =1 +
1
k
2

2
De

2
lh

2
_
1 +
m
i
m
e
k
2

k
2

(20)
+
3k
2
2k
2

_
1 +

2
e

2
ce
_
k
2

2
Di
_
,
k

/k


_
m
e
/m
i
.
The lower-hybrid frequency is dened as
2
lh
=
2
i
(1 +

2
e
/
2
ce
)
1
. The term in brackets results from the large ar-
gument expansion of the derivative of the plasma dispersion
function Z

(
i
) =2[1 +
i
Z(
i
)] with =/(
i
k
Di
) the
argument for the immobile ions. This response function is
formally of the same structure as the ion-acoustic response
function Eq. (5). Thus dening

2
lh
(k) =

2
lh
1 +1/k
2

2
De
_
1 +
m
i
m
e
k
2

k
2

(21)
+
3k
2
2k
2

_
1 +

2
e

2
ce
_
k
2

2
Di
_
the whole formalism developed for ion-acoustic waves can
be applied to lower-hybrid waves. We write for the inverse
response function
1

lh
(, k)
=
k
2

2
De
1 +k
2

2
De
_
1 +

2
lh
(k)

2
lh
(k)
_
. (22)
Again, the Debye-screening term outside the brackets is of
no interest at distances r >
De
. The contribution to the wake
potential comes from the integral in Eq. (9) with
ia
(k) re-
placed by
lh
(k) and k
z
k

=/v

, where k

,
Ann. Geophys., 32, 975989, 2014 www.ann-geophys.net/32/975/2014/
R. A. Treumann and W. Baumjohann: Electron coagulation 981
=m
e
/m
i
, for nearly perpendicular wave propagation. Per-
forming the integration reproduces a form similar to
Eq. (10)

lh
(z, , t ) = (23)
e
2
De
8v
0
_
dk
2

J
0
(k

)
lh
(v, k

)
1 +
2
De

2
lh
(v, k

)/v
2
+k
2

2
De
sin
_

De

lh
(v, k

v
(z vt )

De
_
but now including the more complicated lower-hybrid fre-
quency Eq. (21). We simplify the lower-hybrid frequency
by observing k
2

/k
2

1. In dense plasma the last term in


the brackets becomes k
2

2
Di
(
2
e
/
2
ce
) k
2

r
2
ce
which is of
the order of the electron gyroradius-to-wavelength squared,
being small for completely magnetised electrons. Hence,

2
lh

2
De
2V
2
A
(v
e
/c)
2
c
2
lh
will be used in the factor in front
of the sine function. The lower-hybrid wave in this case prop-
agates at the Alfvn speed V
A
corrected by the ratio of elec-
tron thermal to light velocity. In this approximation and with
=k

De
, we have for the lower-hybrid dispersion relation

2
De

2
lh
(v, k)
c
2
lh

2
(1 +
2

2
De
/v
2

2
)
1 +k
2

2
De
_
1 +
k
2

k
2

_
(24)

2c
2
lh

2
1 c
2
lh
/v
2
which is to be used in the above integral in the long-
wavelength approximation k

De
< 1 and v c
lh,
=
c
lh
/

> c
lh
. The last version on the right results from it-
erating the frequency =
lh
(v, k). Within these approxi-
mations and restricting to the interval 1 for long wave-
lengths, the potential becomes

lh
(z, , t ) C
lh
_

2
d J
0
( )
1 +
2
(1 +2c
2
lh,
/v
2
)
sin
_

lh

_
(25)
C

lh
1/

3
_
0

2
d sin(
lh
)
C
lh
=
C

lh
J
0

=
e
2
0

De
c
lh
(v
2
c
2
lh
)
1/2
, (26)

lh

2c
lh,

(v
2
c
2
lh,
)
1/2
< 0,
where (z vt )/
De
. One may note that in the only in-
teresting long-wavelength regime the factor multiplying
2
in the denominator is at most 3. In order to neglect the en-
tire term
2
(1 +2c
2
lh
/v
2
) 1 and being able to analytically
solve the integral one thus requires that the upper limit of the
integral is taken as < 1/

3. Averaging the Bessel func-


tion over this interval again produces the numerical factor
J
0
. In the argument of the sine function the larger paral-
lel wave velocity c
lh,
=c
lh
/

appears. It is due to the


higher phase velocities of the lower-hybrid waves parallel
rather than perpendicular to the magnetic eld, while the test
electron moves along the magnetic eld at velocity v c
lh,
being in resonance with the wave.
With these assumptions the integration of the sine func-
tion with respect to k

can be performed as before and an


attractive wake potential is obtained under similar conditions
as for the ion-acoustic wave Eq. (15):

lh
()
C

lh
3|
lh
|
cos

lh

3
(27)

eJ
0

12
0

De

_
v
2
c
2
lh,
v
2
c
2
lh
_
1
2
cos
_ 2c
lh,

_
3
_
v
2
c
2
lh,
_
_
.
This potential becomes negative for
1
2
<|
lh
|/

3 <
3
2
mod 2, in which case it attracts a neighbouring par-
allel electron. An attractive potential requires v c
lh
/


43 c
lh
in an electronproton plasma. As a consequence the
fraction under the square root does not shorten out but be-
comes small of the order of o
_
1 c
2
lh,
/v
2
_
O
_

_
. Un-
der the condition on the argument of the cos function the
amplitude of the potential is of the order of
C

lh
3|
lh
|

eJ
0

3
0

De
c
lh
v
(28)
which is small of the order of the ratio c
lh
/ v

. Never-
theless, lower-hybrid waves may attract some resonant elec-
trons in parallel motion along the magnetic eld. In the trans-
verse direction any electrons gyrate and thus are insensitive
to attraction. Any potential generated will just cause a cross-
eld electron drift weakly contributing to local current uc-
tuations.
2.4 Buneman mode-mediated inter-electron potential
A most important plasma wave is the current driven non-
magnetic Buneman mode (Buneman, 1958, 1959). It oc-
curs under conditions of collisionless shocks, in collision-
less guide eld reconnection (Drake et al., 2003; Cattell et
al., 2005), and in auroral physics, in all cases producing
highly dynamical localised electron structures of the type of
BGK modes which trap electrons and cause violent effects
in plasma dynamics (Newman et al., 2001). Again account-
ing for the presence of test electrons, the dielectric response
function of the Buneman mode is
(, k) =1 +
1
k
2

2
De

2
i

2


2
e
(ku)
2
, (29)
www.ann-geophys.net/32/975/2014/ Ann. Geophys., 32, 975989, 2014
982 R. A. Treumann and W. Baumjohann: Electron coagulation
with u the current drift velocity of the electrons, and k the
one-dimensional wave number. For the Buneman mode one
has k
e
/u and
i

e
. Under these conditions the
(nonlinear) version of the Buneman response function be-
comes

B
(, k) =1 +
1
k
2

2
De
+

2
_

_
3
_
1 +
3
2
N
N
_
. (30)
The Buneman dispersion relation is obtained as

3
B
(k) =
3
e

2
k
2

2
De
1 +k
2

2
De
_
1 +
3
2
N
N
_
, (31)
where we retained the nonlinear modulation term propor-
tional to the density variation N. In equilibrium it becomes
N/N =(
0
/4m
i
c
2
ia
N)|E
B
|
2
which is proportional to the
Buneman electric eld intensity causing hole formation. In
the following this term will be neglected. We note that the so-
lution
B
(k) =(
B
) +i(
B
) has a non-negligible imag-
inary part which must be taken into account. Inverting the
response function yields,
1

B
(, k)
=
k
2

2
De
1 +k
2

2
De
_
1 +

3
B
(k)

3
B
(k)
_
. (32)
2.4.1 Attractive potential in linear theory
The structure of this function is more complicated than in
the ion-acoustic case which is due to the higher power in fre-
quency and its imaginary part. This function is to be used
in Eq. (1). Again, the rst term just reproduces the Debye
screening and can thus be dropped. In order to treat the in-
tegral of the second term, we again assume that the electron
moves in z direction at velocity v. Rewriting the integral in
cylindrical coordinates and replacing k

=/v as required
by the delta function, we nd

B
(z, , t ) =
e
2
De
16
0
_
ddk
2

J
0
(k

)
1 +
2

2
De
/v
2
+k
2

2
De
(33)

3
B
(k

, v) exp[i(z vt )/v]
_

3
B
(k

, v)
_ .
Treating the integral is complicated by the third power
of the frequency. It requires expansion of the last term
into a Laurent series. Since we know that
B
is a solu-
tion of the dispersion relation, the denominator can be ex-
panded around =
B
yielding in the denominator 3
2
B
(

B
)
_
1 +(
B
)/
B
+
1
2
(
B
)
2
/
2
B
_
. The bracket can
then be further expanded. Ultimately applying the residuum
theorem, only the rst term survives producing

B
(z, , t ) =
ie
2
De
4
0
_
k

dk

1 +
2
B

2
De
/v
2
+k
2

2
De
(34)
J
0
_
k

B
_
k

, v
_
exp
_
i
B
(z vt )/v
_
and we must, for () > 0, require that z vt > 0 and in-
tegrate over the positive frequency half-space. Indeed, solv-
ing the dispersion relation still, for completeness, keeping the
nonlinear term, we obtain the usual Buneman frequency and
growth rate
(
B
)

e
(1 +1/k
2

2
De
)
1/3
_

16
_1
3
_
1 +
1
2
N
N
_
, (35)
(
B
) =

3 (
B
).
Hence, electrons in resonance with the wave lag slightly be-
hind the wave. The integral may be written as a derivative
with respect to =(z vt )v
e
/v
De
. Further simplifying the
denominator and dening =
e
(/16)
1/3
(1 +N/2N)
0.03
e
(1 +N/2N) the integral becomes

B
(z, , t )
e
4
0

De

1
_
0
dJ
0
( ) (36)
exp
_

2
3
_

3 i
_

_
.
Changing variables and solving for the integral and restrict-
ing to the dominant term, we nd that

B
(, =0)
3
4
e

De

exp
_

3
_
(37)

_
cos
_
+

6
_
+i sin
_
+

6
_
_
holding for > 0. Only the real part of the potential has
physical relevance, the imaginary part causing a spatial undu-
lation along of wavelength 6 /11. We thus nd that the
potential can indeed become attractive when the cos func-
tion is negative, that is, in the interval
1
3

4
3
and
for resonant electrons lagging slightly behind the wave. This
last condition can also be written

3
0.03
v
e
v
|z vt |

De

4
3
. (38)
Such electrons are presumably trapped in the wave potential
well which connes themto the interior of holes generated by
the Buneman mode. For the distance on which the potential
is attractive the last expression yields

(z vt )

att
10 (v/v
e
)
De
. (39)
For the Buneman mode one requires that u > v
e
. Electron
holes arising from Buneman modes extend up to several

De
(Newman et al., 2001). They are thus well capable of
allowing trapped slow electrons of velocity in the narrow in-
terval v
e
< v < u to experience attracting inter-electron po-
tentials and, in principle, form classical pairs or larger
compounds.
Ann. Geophys., 32, 975989, 2014 www.ann-geophys.net/32/975/2014/
R. A. Treumann and W. Baumjohann: Electron coagulation 983
This attractive potential caused by the Buneman mode
is weak. This is obvious from the exponential factor
exp(

3). Inserting for and using the condition


Eq. (38) with v
e
/v < 1, it is found that this factor is of the
order of 0.007. Moreover, though the imaginary part of
the potential plays no role in the sign of the potential, it im-
plies a periodic modulation of the electric eld along z which
is obtained when taking the derivative E
z

B
/ . This
modulation is, however, spatially damped away by the expo-
nential factor. Hence the potential becomes indeed weakly
attractive only in the near zone given by Eq. (38).
2.4.2 Undamped contribution of the singularity in
For completeness we check for the resonant contribution of
the integral. This is complicated by the wave-number de-
pendence of the Buneman dispersion relation induced by the
presence of the test charge. Iteratively, the remaining k de-
pendence of the Buneman dispersion relation is reduced to
k

only. With this in mind, the denominator in Eq. (34) put


to zero becomes

2
+
B

4
3
+1 =0 (40)
where
B
=(/16)
2/3
(1+i

3). Dening

=
2/3
, this be-
comes a third-order equation

3
+
B

2
+1 =0 the solution
of which is complicated by the complexity of the coefcient

B
. In general it has one real and two complex solutions. The
real solution is of no interest as it only contributes to a weak
deformation of the Debye sphere. In order to obtain the com-
plex solutions, we may refer to the smallness of 1 in
the long-wavelength regime and neglect the third-order term.
Solving for

yields four solutions

1,...,4

3
B
_
1
2
=
_
1
i
_
a
1
e
i/4
(41)
with a (/16)
1/2
. Checking with these solutions for the
exponential in Eq. (34) it can be shown that of the solutions
in the upper row only the solution
1
with the + sign con-
verges. Its pole lies in the lower-half plane. The pole of the
converging lower-row solution
3
lies in the upper half plane
and corresponds also to the + sign. The denominator of the
integral can thus be written (
2

2
1,2
)(
2

2
3,4
) where only
the solutions
1
,
3
contribute. The integrand splits into the
two resonant terms

2
1,2

2
3,4

4/3
+1
=
_
1
(
1
)(
2
)

1
(
3
)(
4
)
_
. (42)
The rst term on the right contributes a factor 2i, the
second a factor 2i which, when including the minus sign
in the bracket, yields a common factor 2i. Moreover,

2
1,2

2
3,4
=2i/a
2
. Hence a factor 2ia
2
/2i =a
2
re-
sults. Since
2
=
1
,
4
=
3
, a further factor a/2 appears
which makes a nal common factor a
3
/2. In addition, the
two singular terms are multiplied: the rst by e
i/4
and the
second by e
i/4
.
Solving for the residues at small 0, one again obtains
a complex potential
4
which, after some simple but lengthy
algebra, yields for the real and imaginary parts of the singular
integral contribution to the potential

sg
B
( ) (43)
Acos

2
_
sinh

3
cos

12
cosh

3
sin

12
_
sin

2
_
sinh

3
cos

12
+cosh

3
sin

12
_
,

sg
B
( )
Acos

2
_
sinh

3
cos

12
+cosh

3
sin

12
_
sin

2
_
sinh

3
cos

12
cosh

3
sin

12
_
,
where A ea
2
/4
0

2
De
. For positive and small, 0
< 1, that is, in the domain of largest interest, the dominant
term of the real part becomes

sg
B
( ) A sin(/12)cos
_
1
2

_
cosh
_
/

3
_
(44)
This contribution to the electrostatic potential is both at-
tractive and not exponentially damped. It thus represents an
important, in fact the dominant, contribution to the attrac-
tive electric force exerted by Buneman modes. In contrast to
ion-acoustic wave mediated potentials, the singularity of the
integral in presence of the Buneman mode therefore adds
substantially to the attractive pairing potential in the near
zone 0 which acts on the slow electron component and
causes electron coagulation possibly leading to the formation
of electron compounds or macro-electrons in Buneman tur-
bulence. As before, the imaginary part of the potential con-
tribution merely causes a spatial undulation of the potential.
2.4.3 Weakly nonlinear Buneman mode
The Buneman mode is a strong wave in the sense that it
grows very fast, actually close to explosive growth. This has
a profound effect on the plasma which appears as hole forma-
tion, with N =0 reacting on the wave. In a simplied the-
ory this reaction is most easily described by taking the varia-
tion of the Buneman frequency =(
B
) with respect to
both density and wave number (Treumann and Baumjohann,
1997). The latter is varied with respect to k
B
=
e
/u, yield-
ing
(
B
)
_
1
3
_
u
v
e
_
2
k
2

2
e
+
1
2
N
N
_
,
k < k
B
,
v
e
< u
(45)
4
The case =0 produces a series of Bessel functions of com-
plex argument which just provides another severe mathematical
complication without adding to any further physical insight.
www.ann-geophys.net/32/975/2014/ Ann. Geophys., 32, 975989, 2014
984 R. A. Treumann and W. Baumjohann: Electron coagulation
It is customarily interpreted as an operator equation acting
on the Buneman mode electric eld envelope E(z, t ). This
procedure results in a nonlinear Schrdinger equation
5
_
i

+
1
2

2
z
+

E( z, )

2
_
E( z, ) =0 (46)
where =(
B
)t, z =

6
e
z/u. The coefcient =

0
/8m
i
c
2
ia
N of the nonlinear term results from the density re-
sponse of the plasma to the presence of the nite amplitude
Buneman wave.
The stationary solution in the comoving frame of the
Buneman wave is, in this approximation, a caviton of am-
plitude E( z) =E
m
/cosh( z/L) of width L =1/
_
E
m

1/2
_
and maximum dip amplitude E
m
. In this comoving frame
=
0
/8m
i
(c
ia
u)
2
N
0
/8m
i
u
2
N for u c
ia
. Electrons
trapped in the cavitons have velocities
v < (
0
/m
e
N)
1
2
E
m
. (47)
Oscillating back and forth in the caviton, electrons in their
backward traveling phase of motion or near their turning
points at the boundaries of the cavitons are sensitive to attrac-
tion. Hole-passing electrons in either direction, on the other
hand, are not in resonance and thus do not experience any
attraction.
2.4.4 Strong nonlinearity: electron hole effect
These arguments hold for weakly modulated Buneman
modes. As noted above, the Buneman mode is, however, a
strong wave which during its evolution causes electron holes
to evolve from BernsteinGreenKruskal (BGK) modes
which cannot be described by the above approximate weakly
nonlinear theory. In this case the variation of the density
|N/N| 1 (48)
becomes itself of the order of the density.
Under this condition one may assume that in the Buneman
dispersion relation Eq. (35)
N/N |E( z, )|
2
(49)
in which case the effective plasma frequency
=
e
_
1

N
N

_1
2
_

16
_1
3

e
_

16
_1
3
(50)
5
Strictly speaking, for the strongly growing Buneman mode one
should also account for the variation of the imaginary frequency
(growth rate). This results in a complex nonlinear Schrdinger
equation
_
i

+
1
2

2
z
+
_
1 +i

3
_

E( z, )

2
_
E( z, ) =0.
Equations of this kind are known from LandauGinzburg theory in
many-particle quantum statistics but have not yet been considered
in plasma physics.
becomes very small, yielding that 0 in the exponential
damping factor in electron holes vanishes this is a very im-
portant fact.
The attractive potential under the condition of elec-
tron hole generation becomes undamped, and the condition
Eq. (38) assumes full validity. This is the case when the
Buneman mode evolves into BGK-mode electron holes as
observed in several places in space, the aurora and strong
collisionless reconnection. It then becomes capable of con-
tributing to the proposed classical pairing or coagulation
of electrons inside an electron hole affecting the low-velocity
trapped-electron component. As before, passing electrons are
immune to any attractive potentials and coagulation.
2.5 Summary
In this paper we examined four types of plasma waves for
their capability of causing attraction between two electrons
in close distance. All four wave families can, under certain
conditions, contribute. Attraction is a purely classical ef-
fect which just resembles real quantum pairing of electrons
in electronphonon interaction at low temperatures in solid
state physics. Nevertheless the mechanisms are similar in the
sense that they imply electron-wave interaction. This lets one
ask whether the multiple classical pairing (coagulation)
6
may have observable effects. In the last section we present a
few speculative hypotheses in this direction.
3 Discussion and conclusions: possible effects
Of all the plasma waves checked, the most promising candi-
dates for pairing are ion-acoustic waves. These had been
proposed already by Nambu and Akama (1985) in view of
application in non-magnetised dusty plasma. Such waves
populate the solar wind and magnetosheath where they might
produce attractive potentials and generate a minor compo-
nent of heavy cold coagulated electrons. Electron-acoustic
waves, because of their very strong damping, are no really
good promising candidate. Lower-hybrid waves propagating
into a nearly perpendicular direction have weak parallel po-
tentials only, though we have given arguments for attractive
potentials generated by them as well. Large amplitude linear
Buneman modes, a particularly important wave mode, suffer
from exponential damping.
However, Buneman modes when evolving into electron
holes from BGK modes, the density modulation becomes
large and as argued above the exponential damping factor
6
Another term in place of multiple classical pairing or coag-
ulation would be bunching. However, bunching has the connota-
tion of particles being bunched into a common dynamical phase,
for instance in their gyrational motion in an external magnetic eld
as used in free-electron laser and electron-cyclotron maser theory.
Since coagulation meant in this paper is a different process, we pre-
fer avoiding use of this term.
Ann. Geophys., 32, 975989, 2014 www.ann-geophys.net/32/975/2014/
R. A. Treumann and W. Baumjohann: Electron coagulation 985
is strongly reduced. In view of applications, this is the most
interesting case involving Buneman modes for it causes sus-
ceptible attractive potentials evolving in the interior of an
electron hole. Our calculations do, however, apply only to
single wave modes trapped inside a hole. In order to account
for the effect of the modulated wave spectrum it would be
necessary to integrate over the hole-trapped wave spectrum
under the restrictive condition imposed by the resonance con-
dition limited to the necessary condition for producing at-
tractive potentials. The latter two separate out just a small
group of particular resonant electrons from the trapped elec-
tron component for each of the wave numbers k in the spec-
trum of hole-trapped waves. Electrons at the bottom of the
hole potential are clearly not involved in the resonance and
attraction; they are at rest. This all implies that the number of
resonant electrons ready for attraction will be very small. It
consists of the fraction, say
res
1 of resonant particles sat-
isfying the attractive condition cut out of the trapped electron
distribution located in a shell of (negative) attractive poten-
tial just outside the Debye sphere being of spatial extension
r
D
with < 1. For a trapped electron density N the
fraction of electrons per Debye sphere in this narrow shell
is N. Of these just a fraction
res
is in resonance. This
yields per Debye sphere a fraction of
res
N N avail-
able for compound formation. Clearly this fraction is very
small.
In principle, one could also think of electromagnetic
plasma waves causing attractive potentials. The candidates
would, however, only be electromagnetic waves possessing
sufciently large magnetic eld aligned electric elds. Natu-
rally, low-frequency electromagnetic waves have relativisti-
cally small electric components. Hence, the only candidates
could be highly oblique whistlers, which generally resemble
lower-hybrid modes and need not be discussed further, ki-
netic Alfvn waves which are known to possess large-scale
and comparably strong electric elds, in particular in the au-
roral region, but also on the ion-inertial scale near reconnec-
tion sites, and the extraordinary electromagnetic mode. Of
these, only kinetic Alfvn waves are worth being checked.
This will be reserved for a separate investigation. It requires
an electromagnetic treatment involving the magnetic vector
potential.
3.1 Mass and charge of prospective coagulations
In solid state physics, electron pair formation is related to
super-uid and super-conducting behaviour of matter (Fet-
ter and Walecka, 1971; Huang, 1987; Ketterson and Song,
1999) in metals and semi-conductors which are based on the
fact that pairing electrons become Bosons with either vanish-
ing or integer spin. At low temperatures they are capable of
releasing their kinetic energy until condensing in their lowest
energy level which, in a magnetic eld, is the lowest Landau
level
1
2

h
ce
(Landau, 1930).
Classical pairing produces compounds of electrons
which attract each other. Each electron may become sur-
rounded by other weakly bound electrons. This happens on
the scale of the Debye length (Fig. 1). Such compounds have
large masses and charges
m

=n
com
m
e
, q

=n
com
e (51)
with n
com
the number of electrons in the compound, but con-
stant charge-to-mass ratio e / m
e
. The mass increase affects
thermal speed, momentum and kinetic energy. The charge
will be compensated by the unchanged number of ions.
It remains an open question whether or not classical pair-
ing or coagulation will actually take place. As noted, the
presence of an attractive potential which is responsible for
the attractive force between neighbouring electrons, is just
the necessary condition for subsequent coagulation of elec-
trons to form classical pairs or larger electron compounds.
Real compound production requires, in addition, the observa-
tion of the sufcient conditions. These are more complicated
to investigate than the mere though already quite involved
generation of attractive potentials given in this paper.
The necessary (attractive potential) condition consists of
two parts. In brief, for an electron experiencing the attractive
potential force these are the resonance condition imposed on
the electron and the requirement that the electron is localised
at the right location in space where the potential is attractive.
The former depends on the wave mode. The latter, as has
been noted, says that for becoming attracted a resonant elec-
tron must be located at a radial distance from the attracting
electron outside but very close to the latters Debye sphere.
Inside the Debye sphere the potential is repulsive. At dis-
tance larger than the Debye sphere the attractive force rapidly
decays with distance. Attraction is, hence, limited to a thin
shell of some thickness d located at radial distance
De
from
the attracting electron. Solving the sufcient conditions not
only requires determining the attractive force (taking the ra-
dial gradient of the attractive potential) but also integrating
in momentum space over the resonant particle distribution in
presence of a given wave spectrum, and integrating spatially
over the attractive shell.
Such a calculation can only be done numerically and re-
mains to be a formidable task. Still, it does not yet provide
information about the (average) number n
com
of particles in a
single compound. This number depends on how many Debye
spheres become correlated in the attraction process, a num-
ber which is not known a priori.
3.2 Electron cooling
Since only a small number of electrons participate in attrac-
tion, their distribution function is just a narrow cut out of
the distribution of all electrons available in the volume. Fig-
ure 2 sketches the situation for the case of ion-acoustic waves
which may originally have been unstably excited in a ther-
mally imbalanced ionelectron plasma T
e
> T
i
as shown by
www.ann-geophys.net/32/975/2014/ Ann. Geophys., 32, 975989, 2014
986 R. A. Treumann and W. Baumjohann: Electron coagulation
F(v)
F (v)
i
e
F (v)
com
v, /k 0 c
s
v
d
ia
|E(/k)|
2
Figure 2. Phase space of ion-acoustic waves excited by the ion-
acoustic instability. Shown are the one-dimensional background ion
F
i
(v) and electron F
e
(v) distributions. Ion-acoustic wave with spec-
trum |E(/k)|
2
ia
evolve at phase velocities above the minimum
of c
ia
in the range c
s
< /k < v
d
. The electron pair distribution
function F
com
(v) produced in the high-phase speed range is shown
schematically in blue. One may note the very low velocity spread
of the pair distribution indicating the much lower pair than original
electron background temperature T
pair
T
e
.
the two distributions F
i
(v), F
e
(v) in one-dimensional phase
space. This is the canonical case of ion-acoustic wave ex-
citation. The ion-acoustic wave spectrum exists in a narrow
phase velocity range c
s
<
ia
/k < v
d
as shown in red. c
s
is
the minimum of the ion-acoustic wave phase velocity. At-
tractive potentials can be generated only at nite wave ampli-
tudes and for electron velocities v c
ia
. The resulting low-
density pair distribution is shown in blue.
One may note the very narrow velocity spread of the at-
tracted distribution F
com
(v) which is at most as wide as the
ion-acoustic wave spectrum corresponding to a rather low-
temperature T
com
T
e
of the electrons participating in at-
traction and available for compound formation. Their maxi-
mum speed is sufciently far below v
d
. Compound distribu-
tions are cold.
Figure 3 is for the Buneman case which holds for u v
e
.
The excited spectrum in this case is as well extremely nar-
row with phase velocity spread of the same order as the lin-
early excited Buneman waves, that is, v uv
e
. Buneman
modes are excited for u v
e
just above the electron thermal
speed. One may note that the reactively growing wave read-
ily reduces any speed u v
e
to values marginally exceeding
v
e
. Consequently, the compound distribution which is at most
as wide as the Buneman spectrum, also has low-temperature
T
com
m
e
|u v
e
|
2
T
e
.
Buneman modes are known to evolve into electron holes.
In this case the hole-trapped electrons become heated in the
trapped wave spectrum. Clearly, the prospectively attracted
electrons or compounds formed will, in the long term, par-
ticipate in this heating. However, formation of attractive po-
tentials and attraction are almost immediate processes in the
interaction of resonant electrons with one of the propagat-
ing waves trapped in the hole. This process is much faster
F(v)
F (v)
i
e
F (v)
com
v, /k 0 v
e
v
d
B
|E(/k)|
2
caviton/hole
formation
Figure 3. Phase space of Buneman modes excited with spectrum
|E(/k)|
2
b
and wave number k =
e
/v
d
evolving at phase veloci-
ties above v
e
for u v
d
> v
e
. The spectrum is very narrow in phase
velocity. The electron pair distribution function F
com
(v) produced
(blue) has similar width as the spectrum and is thus much colder
than the original electron distribution. In caviton formation the spec-
trum extends to much larger phase velocities which, however, has
no remarkable effect on the pair distribution.
than any heating. It selects out a small number of resonant
electrons from the trapped electron distribution to form com-
pounds of at least two electrons resulting in a cold elec-
tron compound component. In the long term, when collision-
less heating sets in (due to phase mixing in the hole-trapped
Buneman wave spectrum), the compounds should also par-
ticipate in the heating becoming destroyed (due to internal
oscillations excited by the higher external temperature) when
the compound temperature T > 2e exceeds the potential of
attraction forcing the compound electrons to join back into
the trapped population. From this point of view compound
formation of electrons in holes will occur preferably in the
initial state of the hole before the heating phase sets in. It is
thus questionable whether the cold trapped component will
survive at all. On the other hand, attractive potentials could
as well be generated at later times if only wave modes remain
trapped and survive after phase mixing. Hence, the case re-
mains unclear.
In all these scenarios the possibly generated compound
plasmas turn out to be of low temperature. Classical pair-
ing in collisionless plasma is a non-radiative cooling mech-
anism acting on a small number of resonant plasma electrons
being sensitive to the attractive potential.
3.3 Secondary electron-acoustic wave excitation
The rst side-effect of cooling is that the plasma after pair-
ing consists of a two-temperature electron plasma of con-
stant charge-to-mass ratio and cold particle density less than
plasma density. Such a two-electron temperature plasma ex-
cites high frequency/high velocity electron-acoustic waves
which are radiated away from the coagulation region. In prin-
ciple the electron-acoustic waves could be observed if exci-
tation is strong enough to overcome the strong damping of
the electron-acoustic waves.
Ann. Geophys., 32, 975989, 2014 www.ann-geophys.net/32/975/2014/
R. A. Treumann and W. Baumjohann: Electron coagulation 987
3.4 Electron re-magnetisation in Buneman waves
In the case of the Buneman instability, coagulations form
from the slower electron component (Fig. 3) with attrac-
tive elds being exponentially damped on large scales.
Hence attraction will preferentially be relevant inside BGK-
mode electron holes affecting the trapped electron compo-
nent in the cavity/hole which forms when the Buneman
mode evolves nonlinearly. The coagulations constitute a low-
density electron population of temperature substantially be-
low T
e
which remains hole trapped, unable to escape.
It is interesting to speculate on the importance of this kind
of Buneman-induced compound formation in reconnection.
Guide eld simulations and observations strongly suggest
that the Buneman mode causes generation of electron holes
during reconnection (Drake et al., 2003; Cattell et al., 2005).
In the geomagnetic tail reconnection region, electron temper-
atures are lower than ion temperatures inhibiting ion-acoustic
wave excitation. Electrons in this case are nonmagnetic in-
side the electron diffusion reconnection site (electron ex-
haust) being accelerated in the cross-tail eld. In presence of
a guide eld this acceleration causes high guide eld aligned
velocities exceeding the thermal electron speed, a situation
favouring the excitation of Buneman modes and generation
of chains of electron holes along the guide eld.
Production of a surviving cold dilute compound-electron
plasma in the Buneman excited electron holes in the ion-
diffusion region and near the reconnection site implies re-
magnetisation of the hole-trapped nonmagnetic electrons un-
til their gyroradius r
com
e
drops below the inertial scale of
the plasma. This is easily conrmed by forming the ratio of
the compound gyroradius to the bulk electron inertial scale

e
=c/
e
. Accounting for the constancy of the compound
charge-to-mass ratio, this ratio can be written as
r
com
e

e
T
com
T
e
_1
2
, (52)
where
e
=2
0
NT
e
/B
2
is the bulk plasma electron-. Since
T
com
/T
e
is substantially less than one, the compound elec-
trons regain magnetisation in the reconnection electron ex-
haust where the bulk electrons remain to be nonmagnetic, an
effect which necessarily affects the evolution of reconnection
in several ways. One effect is that magnetised electrons trans-
port magnetic ux into the bulk-electron diffusion region
thereby enhancing reconnection. Their stronger magnetisa-
tion also modies reconnection. Moreover, electron holes
forming chains along the guide eld naturally contribute to
amplication and deformation of the guide eld on the spa-
tial scale of the holes, a process which self-consistently gen-
erates localised non-zero B
z
components in the current sheet
centre. However, because of the expected very low number
of compounds formed, the effect will be rather small if not
completely negligible.
In summary, though attractive potentials will certainly
arise in various wave particle interactions in plasma, the
number of electrons which may under favourable circum-
stances coagulate and cool down to low temperatures will in
all cases be very small and therefore ineffective for plasma
processes. Unfortunately, attraction though a natural process
does not provide any natural mechanism of large macro-
particle number generation. It would be interesting to inves-
tigate whether particle bunching in low-frequency electro-
magnetic waves (whistlers, kinetic Alfvn waves, etc.) might
be another option of imposing a common dynamic behaviour
on large numbers of electrons to perform correlated dynam-
ics and appear as macro-particles. Observation of very dilute
cold electrons in the presence of high levels of plasma wave
activity would, however, indicate ongoing attraction and co-
agulation.
www.ann-geophys.net/32/975/2014/ Ann. Geophys., 32, 975989, 2014
988 R. A. Treumann and W. Baumjohann: Electron coagulation
Appendix A: No macro-quantum condensation effects
Here we demonstrate that the classical condensate will not
undergo any quantum condensation (i.e. it will not become a
BoseEinstein macro-condensate) inside a caviton. The low-
est reachable energy level is at the bottomof the caviton. This
can be taken as zero-energy level for the composed electrons.
Furthermore, composition temperatures are low, the order of
a fraction of an eV in the classical condensate. Assuming
that the composed electrons obey a Bose distribution in the
caviton potential we thus write for their density
dn
com
d
=
1
4
2
_
4m
e

h
2
_3
2

1
2
e
(2e)
1
. (A1)
The chemical potential 0 is compensated by the cavi-
ton potential at inverse temperature
0
which is calculated
fromthe total density of the trapped particles. The upper limit
of the integral can be assumed at innity. Hence (Fetter and
Walecka, 1971)
T
0

1
0
1.6

h
2
m
e
n
2
3
com
10
19
n
2
3
com
eV. (A2)
Since any densities are very low in space plasmas this limit
temperature on any BoseEinstein condensation is practi-
cally zero in comparison to the estimated lowest compound
temperatures T
com
0.1 eV. This precludes that in ordinary
space plasmas any composed electrons produced would form
macro-quantum BoseEinstein condensates, indeed an intu-
itive reasoning. Under the extreme conditions in the plasma
of neutron star crusts with their high nuclear densities such
condensates could possibly occur if Buneman modes would
evolve along the neutron star magnetic eld.
Ann. Geophys., 32, 975989, 2014 www.ann-geophys.net/32/975/2014/
R. A. Treumann and W. Baumjohann: Electron coagulation 989
Acknowledgements. This research was part of a Visiting Scientist
Program at ISSI, Bern, executed by RT in 2006/2007. Hospitality of
the ISSI staff is thankfully acknowledged. The intriguing remarks of
the two anonymous referees concerning pairing, quantum effects,
compound formation, and role of unstable waves were greatly ap-
preciated.
Topical Editor G. Balasis thanks two anonymous referees for
their help in evaluating this paper.
References
Buneman, O.: Instability, turbulence, and conductivity
in current-carrying plasma, Phys. Rev. Lett., 1, 89,
doi:10.1103/PhysRevLett.1.8, 1958.
Buneman, O.: Dissipation of currents in ionized media, Phys. Rev.,
115, 503517, doi:10.1103/PhysRev.115.503, 1959.
Cattell, C., Dombeck, J., Wygant, J., Drake, J. F., Swisdak,
M., Goldstein, M. L., Keith, W., Fazakerley, A., Andr, M.,
Lucek, E., and Balogh, A.: Cluster observations of elec-
tron holes in association with magnetotail reconnection and
comparison to simulations, J. Geophys. Res., 110, A01211,
doi:10.1029/2004JA010519, 2005.
Drake, J. F., Swisdak, M., Cattell, C., Shay, M. A., Rogers, B. N.,
and Zeiler, A.: Formation of electron holes and particle ener-
gization during magnetic reconnection, Science, 299, 873877,
doi:10.1126/science.1080333, 2003.
Fetter, A. L. and Walecka, J. D.: Quantum theory of many-particle
systems McGraw-Hill Publ. Comp., New York, USA, 1971.
Huang, K.: Statistical Mechanics, 2nd Edn., John Wiley & Sons,
New York, USA, 1987.
Huba, J. D., Gladd, N. T., and Papadopoulos, K.: The lower-hybrid-
drift instability as a source of anomalous resistivity for mag-
netic eld line reconnection, Geophys. Res. Lett., 4, 125128,
doi:10.1029/GL004i003p00125, 1977.
Ketterson, J. B. and Song, S. N.: Superconductivity, Cambridge
University Press, Cambridge UK, 1999.
Krall, N. A. and Trivelpiece, A. W.: Principles of Plasma Physics,
McGraw-Hill, New York, USA, 1973.
Landau, L. D.: Diamagnetismus der Metalle, Z. Phys., 64, 629637,
1930.
Nambu, M. and Akama, H.: Attractive potential between resonant
electrons, Phys. Fluids, 28, 23002301, doi:10.1063/1.865284,
1985.
Nambu, M. and Nitta, H.: On the Shukla-Nambu-Salimullah
potential in magnetized plasma, Phys. Lett. A, 300, 8285,
doi:10.1016/S0375-9601(02)00741-7, 2001.
Neufeld, J. and Ritchie, R. H.: Passage of charged par-
ticles through plasma, Phys. Rev., 98, 16321642,
doi:10.1103/PhysRev.98.1632, 1955.
Newman, D. L., Goldman, M. V., Ergun, R. E., and Mangeney,
A.: Formation of Double Layers and Electron Holes in a
Current-Driven Space Plasma, Phys. Rev. Lett., 87, 255001,
doi:10.1103/PhysRevLett.87.255001, 2001.
Shukla, P. K. and Melands, F.: Test charge potential in the presence
of electron-acoustic waves in plasmas, Phys. Scripta, 55, 239
240, doi:10.1088/0031-8949/55/2/013, 1997.
Shukla, P. K., Nambu, M., and Salimullah, M.: Effect of ion po-
larization drift on shielding and dynamical potentials in magne-
tized plasmas, Phys. Lett. A, 291, 413416, doi:10.1016/S0375-
9601(01)00762-9, 2001.
Sitenko, A. G.: Electromagnetic uctuations in a plasma, Academic
Press, New York N.J., 1967.
Treumann, R. A. and Baumjohann, W.: Advanced Space Plasma
Physics, Imperial College Press, London, ch. 11.4, 1997.
Yoon, P. H., Lui, A. T. Y., and Sitnov, M. I.: Generalized lower-
hybrid drift instabilities in current-sheet equilibrium, Phys. Plas-
mas, 9, 15261538, doi:10.1063/1.1466822, 2002.
www.ann-geophys.net/32/975/2014/ Ann. Geophys., 32, 975989, 2014

You might also like