Surface Waves in a Relativistic Plasma Stream Propagating in a Duct: Kinetic Theory ()
1. Introduction
We investigate surface waves generated by a relativistic plasma beam travelling in a duct interfaced with vacuum by using relativistic Vlasov equation. Surface waves propagate along the interface between two different media while being attenuated in the perpendicular direction. Surface waves are the normal modes that are given rise to in bounded plasmas by satisfying the kinetic and the electrodynamic boundary conditions on the interface between two media. The electrodynamic boundary conditions are the connection formulas matching the fields of the two media, which can be mathematically worked out from the governing equations themselves if the density gradient across the plasma and the other side is very steep. “A sharp interface” is synonymous to a theoretically infinite density gradient across the boundary. In this case, the connection formula can be obtained by “infinitesimal integration” across the interface, which is the operation performed on a certain relevant equation in the manner
, where
is a positive infinitesimal. If the quantity
is a perfect differential, this operation yields a non-vanishing surface term that contributes to the connection formula. Usually the surface term is surface charge or surface current, In this way, the well-known electromagnetic and dynamic boundary conditions on the boundary can be derived [1]. We might say that the surface wave in the plasma and the vacuum side wave are two different manifestations of “the same wave” created in an extreme inhomogeneous plasma. The electromagnetic boundary conditions are the same in both non-relativistic and relativistic plasmas.
Integrating the Maxwell equation
over an infinitesimal segment
across the interface
, we obtain
where
signifies the jump across the interface,
is the surface current. In a cold plasma, the normal component of the electric field is discontinuous by the amount of surface charge
:
. On the other hand, we have the relation
due to the charge conservation law, where u is the drift velocity of the plasma in the z direction, and c is the speed of light. Thus we have the characteristic boundary condition in a cold streaming plasma [1] [2]:
(1)
or equivalently,
(2)
where
is the normal to the interface component of
, the electric displacement,
is the z component of the wave vector, and
is the wave frequency. The casual use of
or
in a drifting plasma leads to erroneous results as discussed in earlier works [1] [3]. The physical origin of the boundary relation in Equation (1) is due to the surface current formed in a cold streaming plasma as is evident in the above derivation. Equation (1) is valid in non-relativistic as well as in relativistic plasmas [2].
In bounded Vlasov plasmas, the kinematic boundary condition that is usually referred to as the specular reflection condition is assumed to be satisfied on a sharp boundary [4], regardless of whether the plasma is described by non-relativistic distribution function
or relativistic distribution function
(see Equation (9) below). Most of the works on bounded kinetic plasmas were dealt with non-relativistically. Non-relativistic kinetic theory of surface waves in semi-infinite plasmas is well-known [5] [6]. Non-relativistic kinetic dispersion relation of surface wave in a slab plasma was worked out earlier [7].
In this work, we investigate surface waves of a moving relativistic Vlasov plasma in a duct. We consider an infinite duct formed by intersections of four planes:
and
, with
. We shall discuss the specular reflection boundary condition in terms of
, rather than in terms of
, since the reflection of
is entirely equivalent to the reflection of
in the kinetic equation (see Equation (12) below). Thus the specular reflection conditions require for the distribution function
to satisfy
on
plane
on
plane
on
plane
on
plane.
In our duct-bounded plasma, the kinematic conditions on the four planes are satisfied by introducing extended electric field in the fashion:
(3)
(4)
This scheme is workable if
, the zero order distribution function, is invariant with respect to the reflections
and
, and clearly this reflectional property is satisfied by the moving Maxwellian to be introduced later. Equations (3) and (4) are the important conclusion of the above discussion that is valid in relativistic as well as in non-relativistic kinetic equations.
The function
as defined in Equation (3) is a periodic function of a piecewise continuous function of period “a” extending over the range
with discontinuity at
with a jump of
(say) and with discontinuity at
with a jump of
(say), where n is integer. The profile of the piecewise function
is plotted in the book by Lee [8]. The algebra involved in carrying out the Fourier transform of the piecewise discontinuous functions with aforementioned discontinuous jumps is quite laborious [7]. However, it turns out that, after all the algebraic hard work, the discontinuities that are present in the extended field components
and
in Equations (3) and (4) at the locations
and
and
and
are mathematically (as well as physically) tantamount to placing fictitious surface charges at the corresponding jump locations in the form
(5)
This is the crucial part of improvement in this work as compared with the earlier work [7]. The surface charges are associated with the surface currents by satisfying the charge conservation equation
(6)
Therefore, we can assume the presence of the fictitious surface currents
(7)
The surface charges in Equation (5) and the surface currents in Equation (7) should be included in the Maxwell equations for our duct plasma wave analysis.
In this work, we consider a streaming plasma which moves along the axial direction (z-direction) with a relativistic speed u. In Section 2, we introduce the relativistic Vlasov equation. In Section 3, the boundary value problem satisfying the kinetic and electromagnetic boundary conditions are solved to find the dispersion relation of the duct surface wave. Section 4 furnishes discussions, and compares the dispersion relation with the recent work of Lee and Cho on non-relativistic duct flow [9].
2. Basic Equations
We begin with the relativistic equation of motion for electrons
(8)
where
, the relativistic momentum, is
(9)
with m being the rest mass and
being the particle velocity, and
(10)
An easy way to write down the relativistic Vlasov equation is to note that the particle orbit as given by Equation (8) is the characteristics of the Vlasov equation. Indeed, the characteristic equation of the following equation is Equation (8):
(11)
Linearizing (11) gives
(12)
where
is the zero order equilibrium distribution function which will be specified later. We also have the Maxwell equations for electrons. Ions are assumed to be stationary and only form the neutralizing background.
(13)
(14)
where the plasma current
is
(15)
with N being the zero order equilibrium number density. The fictitious surface current
corresponding to the surface charge S in Equation (5) takes the form
(16)
(17)
(18)
We Fourier transform all the dependent variables in the basic equations, including f,
, and
, and assume the form
. Then Equation (12) gives
(19)
where
reads as the function of
as given in Equation (9). Using Equation (19) in the Fourier transform of Equation (15) yields
(20)
where
is the Cartesian index and
is the Cartesian tensor called Levi-Civita symbol, and repeated indexes are summed over. We evaluate the current in Equation (20) for a cold streaming plasma whose zero order distribution function is
(21)
where
(22)
which corresponds to the streaming velocity
. Integrating by parts Equation (20) gives
(23)
The following relations are useful for evaluating the above integral:
(24)
(25)
Using Equation (21) in Equation (23), we obtain after a long algebra,
(26)
where
is Doppler-shifted frequency.
Using
and writing in components, we have
(27)
(28)
(29)
The array formed by the coefficients of the electric field components is the conductivity tensor, as defined by
.
Combining Equations (13) and (14) with the aid of Equations (27)-(29) yields a vector equation for
in the form
(30)
where
(31)
is the dielectric tensor, and
is the plasma frequency, and
Putting Equation (31) into Equation (30) gives a 3 by 3 matrix equation for
:
(32)
where
is the refractive index, and
is the Fourier transform of the fictitious surface currents in Equation (16):
(33)
where A’s and B’s may be functions of
and
, and the double signs are summed over, and the notations
and
are the summations in Equation (16).
Equation (32) can be inverted after spending a considerable amount of time to solve for
:
(34)
(35)
(36)
where
(37)
After further algebra,
becomes
(38)
Also we obtain
(39)
(40)
In addition, we put
since we investigate the transverse magnetic mode.
3. Boundary Equations
In order to apply the boundary conditions, the electric and magnetic field components in the Fourier
space should be inverted to the fields in the ordinary
space by performing
. The integrals involve infinite series through the surface charge
, but the infinite series are nicely summed at the particular positions corresponding to
and
. Thus, we apply the boundary conditions along the two infinite lines:
and
with
. The two lines correspond to the two seams of the duct which are diagonally opposite.
When the inversion integrals are performed, the following formulas are useful, which can be verified by a simple change of variable, as is shown in earlier work [7]. We have integrals of the type in the inversion integrals
(41)
where
is an even function of
. Then, we have
(42)
(43)
We also have integrals of the type
(44)
Then, we have
(45)
(46)
where
(47)
(48)
Formulas in Equations (42), (43), (45), and (46) are useful for evaluating the integrals. Let us evaluate:
(49)
where we used Equation (42). In the last (also in the later) integral,
hidden in
is
. We have
(50)
where we used Equation (43). Analogous integrations yield
(51)
(52)
In the above (also in the later)
integral,
hidden in
is
.
(53)
(54)
(55)
(56)
We encounter with different type of integral in
which becomes
which we write in the form
(57)
where we used Equation (45), and
(58)
Analogously, we obtain
(59)
where we used Equation (46).
· Vacuum solution.
Vacuum solutions should solve
(60)
and
(61)
Equation (60) is solved by
(62)
with constraint
and
. Furthermore, we assume
since we consider the TM (transverse magnetic) mode.
The vacuum regions corresponding to (or exterior to) the lines
and
, which we designate as (i) and (ii), respectively, are:
Vacuum region (i)
, where we have
(63)
(64)
(65)
(66)
(67)
(68)
Vacuum region (ii)
, where
(69)
(70)
(71)
(72)
(73)
(74)
Putting
or
in the above equations gives the vacuum side values of the relevant quantities.
4. Dispersion Relation
We enforce the following boundary conditions to connect the plasma and the vacuum fields:
,
,
Along line
gives, per Equations (57) and (66),
(75)
where
(76)
gives
(77)
where
(78)
gives
(79)
Along line
gives
(80)
gives
(81)
gives
(82)
In addition, we have, per
and
,
(83)
(84)
Thus, we have 8 equations for 8 unknowns;
.
Eliminating
gives
(85)
(86)
(87)
(88)
where
(89)
Eliminating
and
gives
(90)
(91)
The above two equations for
and
yield the dispersion relation in the form
(92)
where we used
(93)
In regard to the Fourier variables
and
outside the integrals, we imposed the constraint
. Therefore it is convenient to transform
(94)
[
,
inside the integrals are dummy and let them be there as they are.] The transform in Equation (94) satisfies the constraint and the relation
. In fact it can be derived from the latter and the constraint. Then the dispersion relation takes the form
(95)
If either a or
, we recover the slab dispersion relation in the non-relativistic limit [7].
(96)
where
and
are given by Equation (89) with
.
It is recalled that
hidden in
is:
in
-integral and
in
-integral. Thus, let us change the integration variables, both
and
, in Equation (95) to
:
(97)
where
. In regard to the double signs in Equation (97), the upper (lower) signs correspond to the symmetric (anti-symmetric) mode which also occurs in a slab plasma. For a square duct (
), Equation (97) reduces to the form identical with the slab dispersion equation Equation (96), except for factor
. This reduction is due to the x-y symmetry. To recover slab dispersion relation from Equation (92), we take
,
, and put
. We can take
since the y-direction has a translational invariance in a slab.
The duct dispersion relation in Equation (97) can be contour-integrated for a cold plasma, giving
(98)
where
. For the anti-symmetric mode, tanh-function above is replaced by coth-function.
5. Discussion
In a bounded plasma, one way of solving Vlasov equation by satisfying the specular reflection condition is to extend the plasma electric field in the manner of Equation (3). The job of Fourier transforming such as piecewise continuous periodic function, extending to infinity, is laborious. In this work, we present an alternative way of avoiding the hard algebra by placing sheets of fictitious surface charges at the location of discontinuities of the electric field. The magnitudes of the surface charges are undetermined constants, but they can be determined through the connection formula with the vacuum side field—resulting in the dispersion relation of the surface wave. This method enables one to deal with semi-infinite, slab, and duct plasmas in a common work-frame. Taking
in Equation (98) gives
(99)
which is the slab (
) dispersion relation. Taking
in Equation (99) gives the semi-infinite plasma dispersion relation
(100)
Equation (100) agrees with the semi-infinite dispersion relation obtained by Lee [3]. If
, Equation (99) agrees with the slab dispersion relation obtained from the fluid theory worked out by Gradov and Stenflo [10].
For a square duct, putting
in Equation (98) yields
(101)
which is similar to the slab dispersion relation. This is because the complete symmetry between x and y coordinates makes the three-dimensional problem a two-dimensional problem practically.
The Doppler-shifted frequency
appearing in Equation (98) represents the streaming effect in the lowest order of
. The higher order effect enters through the relativistic factor
and
in Equation (31). The relativistic effect manifests itself through the attenuation constant
per
and in the dispersion relation per
. Recently Lee and Cho [9] investigated surface waves in a non-relativistically streaming Vlasov plasma in a duct. Their result is identical with Equation (98) upon putting
.
A visual understanding of the extended electric field may be grasped by plotting the extended function [8]. A useful reference for relativistic Vlasov equation is, among others, Momtgomery and Tidman [11], in which the velocity-version of Vlasov equation (in contrast to the momentum-version employed in this work) is presented in detail. This work may find applications in laboratory or astrophysical situation where electromagnetic waves propagate through certain channels. This work might be useful for analysis of a proton beam travelling in a duct.
Acknowledgements
This research was supported by National Cancer Center Grant (NCC-2110370), Korea.
NOTES
*Former professor of physics at Hanyang University.