TMP 91 E
TMP 91 E
TMP 91 E
h
2
/m
e
T
10
9
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
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
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
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
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
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/
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(
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
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