1. Introduction
We consider surface waves of wave number $k_z$ and frequency $\omega$, generated by a plasma beam travelling in a duct interfaced with vacuum by using the Vlasov equation. Surface waves are given rise to by satisfying the kinetic and electromagnetic boundary conditions on the interface between the plasma and the vacuum surrounding it. The electric and magnetic fields in the plasma are connected with the vacuum side fields by the appropriate connection formula, and the latter can be deciphered from the basic equations themselves that we adopt. The connection formula can be easily worked out mathematically if the density gradient across the plasma and the other side is very steep. We talk about ‘a sharp interface’ if the density gradient is theoretically infinite. 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 $\int _{-\epsilon }^\epsilon (\cdots )\,{{\rm d} x}$. If the quantity $(\cdots )$ is a perfect differential, this operation yields a non-vanishing surface term that contributes to the connection formula. Usually, the surface term is the surface charge or surface current, In this way, the well-known electromagnetic and dynamic boundary conditions on the boundary can be derived (Lee & Cho Reference Lee and Cho1997). In a gross picture, the surface wave and the vacuum side wave are two different manifestations of ‘the same wave’ given rise to in an extreme inhomogeneous plasma.
A cold drifting plasma has a characteristic boundary condition (Lee & Cho Reference Lee and Cho1997, Reference Lee and Cho1999), which can be expressed in the form
or equivalently,
where $[\cdots ]$ signifies the jump across the interface, $D_x$ is the normal to the interface component of ${\boldsymbol D}$, the electric displacement, $u$ is the drift velocity in the $z$-direction, $B_y$ is the tangential component of the magnetic field and $c$ is the speed of light, and the other symbols have the usual meaning. The casual use of $[D_x]=0$ or $[B_y]=0$ in a drifting plasma leads to erroneous results as, discussed in earlier works (Lee & Cho Reference Lee and Cho1997; Lee Reference Lee2005). The physical origin of the boundary relation in (1.1) is due to the surface current formed in a cold drifting plasma, as is evident in the equation.
In a bounded Vlasov plasma, the kinematic boundary condition that is usually referred to as the specular reflection condition is assumed to be satisfied on a sharp boundary (Landau Reference Landau1946). The kinetic theory of surface waves in semi-infinite plasmas is well known (Barr & Boyd Reference Barr and Boyd1972; Alexandrov, Bogdankevich & Rukhadze Reference Alexandrov, Bogdankevich and Rukhadze1984). The kinetic dispersion relation of a surface wave in a slab plasma was worked out earlier (Lee & Lim Reference Lee and Lim2007). In this work, we investigate surface waves of a moving Vlasov plasma in a duct. We consider an infinite duct formed by the intersections of four planes: $x=0,a$ and $y=0,b$, with $-\infty < z<\infty$. Thus the specular reflection condition on the $x=a$ plane, for example, requires for the distribution function $f({\boldsymbol r},{\boldsymbol v},t)$ to satisfy $f(a,y,z, v_x,v_y,v_z,t)=f(a,y,z, -v_x,v_y,v_z,t)$ or, on the $y=b$ plane, $f(x,b,z, v_x,v_y,v_z,t)=f(x,b,z, v_x,-v_y,v_z,t)$, with similar equations for $x=0$ and $y=0$. In our duct-bounded plasma, the kinematic conditions on the four planes are satisfied by introducing an extended electric field in the fashion
This scheme is workable if $f_0({\boldsymbol v})$, the zero-order distribution function, is invariant with respect to the reflections $v_x\rightarrow -v_x$ and $v_y\rightarrow -v_y$, and clearly this reflectional property is satisfied by the moving Maxwellian, to be introduced later.
The function $E_x(x)$, as defined in (1.3a,b), is a periodic function of a piecewise continuous function of period $'a'$ extending over the range $-\infty < x<\infty$ with discontinuity at $x=\pm 2na$ with a jump of $A_1$ (say) and with discontinuity at $x=\pm (2n-1)a$ with a jump of $A_2$ (say), where $n$ is an integer. The profile of the piecewise function $E_x(x)$ is plotted in Lee (Reference Lee2019). The algebra involved in carrying out the Fourier transform of the piecewise discontinuous functions with the aforementioned discontinuous jumps is quite taxing (Lee & Lim Reference Lee and Lim2007). However, it turns out that, after all the algebraic hard work, the discontinuities that are present in the extended field components $E_x(x,y)$ and $E_y(x,y)$ in (1.3a,b) and (1.4a,b) at the locations $x=\pm 2na$ and $x=\pm (2n-1)a$ and $y=\pm 2nb$ and $y=\pm (2n-1)b$ are mathematically (as well as physically) tantamount to placing fictitious surface charges at the corresponding jump locations in the form
This is the vital point of improvement in this work as compared with the earlier work (Lee & Lim Reference Lee and Lim2007). The surface charges are associated with the surface currents by satisfying the charge conservation equation
Therefore, we can assume the presence of the fictitious surface currents
The surface charges in (1.5) and the surface currents in (1.7) should be included in the Maxwell equations for our duct plasma wave analysis.
The basic equations are the linearized Vlasov equation and the Maxwell equations for electrons. Ions are assumed to be stationary and only form the neutralizing background
where
It should be noted that the last term in (1.8) involving the magnetic field vanishes if $f_0$ is isotropic, but in a beam plasma, that term contributes to complexity. The zero-order distribution function is a moving Maxwellian
We Fourier transform the above equations by performing $\int _{-\infty }^\infty \,{\rm d}^3r \exp ({{\rm i}{\boldsymbol k}\boldsymbol {\cdot }{\boldsymbol r}})(\cdots )$ and by assuming $\partial /\partial t={\rm i}\omega$. Then, the wave has a phasor $\exp ({{\rm i}k_zz-{\rm i}\omega t})$. The $(k_z,\omega )$ dependency in the Fourier amplitudes will be suppressed. Equation (1.8) is Fourier transformed to give
where the repeated indexes are summed over. Here, and in the following, the Fourier-transformed variables are expressed by the argument; for example, $f(t)$ and $f(\omega )$ have different dimensions. If $f_0$ is isotropic, ${\partial f_0}/{\partial v_s}\sim f_0v_s$, and (1.16) reduces to the expression for the isotropic case.
Eliminating ${\boldsymbol B}$ between (1.9) and (1.10) gives
where $\varepsilon _{ij}$ is the dielectric tensor in a moving medium
with $\omega _p$ being the plasma frequency. In (1.17), $J_s({\boldsymbol k},\omega )$, the Fourier transform of the surface currents in (1.11), is
where the $A$ and $B$ values may be functions of $k_z$, the double signs are summed over and the notations $\varSigma _0$ and $\varSigma _1$ are the summations in (1.11).
The foregoing formulation can be applied for simpler boundaries: in a semi-infinite plasma, we can suppose a single surface charge sheet $S(x)=A\delta (x)$, and in a slab geometry $0< x< a$, the surface charges are assumed in the form of (1.5) with $B_1=B_2=0$. The dispersion relation of surface waves in a slab plasma was worked out earlier without introducing the fictitious surface charges by directly Fourier transforming the extended electric field in (1.3a,b) (Lee & Lim Reference Lee and Lim2007).
2. Solutions for plasma fields
First, let us calculate $\varepsilon _{ij}$ in (1.18) by using the moving Maxwellian as given in (1.15). The velocity integral therein can be carried out by transforming ${\boldsymbol w}={\boldsymbol v}-{\boldsymbol u}$, and thus performing the integral over the isotropic Maxwellian distribution $f_0(w)$. Using
the velocity integral becomes
where $\omega '=\omega -k_zu$ is the Doppler-shifted frequency. Using (2.1) again for $v_if_0$, eventually (2.2) is integrated to yield
can be written in terms of the plasma dispersion function $Z(\zeta )$
We have $I_1=-\zeta Z/\omega '$, $I_2=-2\zeta ^2(1+\zeta Z)/\omega '^2$ and $I_3=\zeta ^2[\zeta Z-2\zeta ^2(1+\zeta Z)]/\omega '^3$. Then, (1.18) can be written in the form
Equation (2.6) agrees with the dielectric tensor in a moving medium obtained from Lorentz transform (Alexandrov et al. Reference Alexandrov, Bogdankevich and Rukhadze1984). When the beam velocity ${\boldsymbol u}=0$, the dielectric tensor in (2.6) reduces to the well-known result for an isotropic medium. It is convenient to introduce the longitudinal and transverse dielectric permittivities,
In the above expressions, two frequencies are involved: the wave frequency $\omega$ and the Doppler-shifted frequency $\omega '$ hidden in the variable $\zeta$. In terms of $\varepsilon _{L,T}$, we have
where $U_{ij}$ is the tensor associated with the beam velocity components
whose non-zero components are
and the other elements are zero.
2.1. Solution for ${E}({k},\omega )$
where $n^2=c^2k^2/\omega ^2$ is the refractive index.
We inverted the 3$\times$3 matrix in (2.12) to obtain ${\boldsymbol E}$
where
We use (2.18) since we investigate a transverse magnetic mode.
The above expressions for ${\boldsymbol E}({\boldsymbol k})$ need to be Fourier inverted to ${\boldsymbol E}({\boldsymbol r})$ to apply the boundary conditions. To make further development easier and the equations more transparent, we employ the cold plasma approximation. Then we have the cold plasma dielectric permittivities
In a cold plasma, $\varDelta$ in (2.20) becomes
The above expression reduces to, for non-relativistic beam velocity ($u^2/c^2<<1)$,
The field components are
where $J_s$ is given in (1.19).
The Fourier inversion integrals must be carried out to obtain the plasma fields in ordinary coordinate space before the boundary conditions are applied. The integrals involve infinite series through the surface charge $J_s$, but the infinite series are nicely summed at the particular positions corresponding to $x=0,a$ and $y=0,b$. Thus, we apply the boundary conditions along the two infinite lines $(x,y,z)=(0,0,z)$ and $(a,b,z)$ with $-\infty < z<\infty$. 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 (Lee & Lim Reference Lee and Lim2007). We have integrals of the type in the inversion integrals
where $\varPhi (k)$ is an even function of $k_x$. Then, we have
When we have integrals of the type
we have
where
Integrals (2.30), (2.31), (2.33) and (2.34) are useful for evaluating the integrals. We have
where we have used (2.30). In the last (also in the later) integral, $k^2$ hidden in $\varDelta$ is $k^2=k_z^2+k_x^2$. We have
where we used (2.31). Analogous integrations yield
In the above (also in the later) $\int \,{\rm d}k_y$ integral, $k^2$ hidden in $\varDelta$ is $k^2=k_z^2+k_y^2$
We encounter a different type of integral in
which becomes
which we write in the form
where we have used (2.33), and
Analogously, we obtain
where we have used (2.34).
2.2. Vacuum solution
Vacuum solutions should solve
and
Equation (2.50) is solved by
with constraint $k_x^2+k_y^2=k_z^2-{\omega ^2}/{c^2}\equiv \lambda ^2$ and $\boldsymbol {\nabla }\boldsymbol {\cdot }{\boldsymbol B}=0$. Furthermore, we assume $B_z=0$ since we consider the transverse magnetic mode.
The vacuum regions corresponding to (or exterior to) the lines $(0,0,z)$ and $(a,b,z)$, which we designate as $(i)$ and $(ii)$, respectively, are:
Vacuum region $(i)$ $~x<0,~y<0$, where we have
Vacuum region $(ii)$ $x>a,y>b$, where
Putting $(x,y)=(0,0)$ or $(a,b)$ in the above equations gives the vacuum side values of the relevant quantities.
3. Dispersion relation
We enforce the following boundary conditions to connect the plasma and the vacuum fields: $[E_z]=0,\ [B_y]=({u}/{c})[E_x],\ [B_x]=-({u}/{c})[E_y]$
Along line $(0,0,z)$
$[E_z]=0$ gives, per (2.47) and (2.56),
$[B_y]=({u}/{c})[E_x]$ gives
$[B_x]=-({u}/{c})[E_y]$ gives
Along line $(a,b,z)$
$[E_z]=0$ gives
$[B_y]=({u}/{c})[E_x]$ gives
$[B_y]=-({u}/{c})[E_x]$ gives
In obtaining (3.3), (3.5), (3.7) and (3.8), we neglected $u^2/c^2\ll 1$ as compared with unity. In addition, we have, per $\boldsymbol {\nabla }\boldsymbol {\cdot }{\boldsymbol B}=0$ and $B_z=0$,
Thus, we have 8 equations for 8 unknowns; $A_1,\,A_2,\,B_1,\,B_2,\,H_x,\,H_y,\,G_x,\,G_y$.
Eliminating $H_x$, $H_y$, $G_x$, $G_y$ gives
where
Eliminating $B_1$ and $B_2$ gives
and (3.16) and (3.17) yield the dispersion relation in the form
where we used
In regard to the Fourier variables $k_x$ and $k_y$ outside the integrals, we imposed the constraint $k_x^2+k_y^2=k_z^2-\omega ^2/c^2\equiv \lambda ^2$. Therefore, it is convenient to transform
($k_x$, $k_y$ inside the integrals are dummy variables and are left as they are). The transform in (3.20a,b) satisfies the constraint and the relation $ak_x=bk_y$. In fact, it can be derived from the latter and the constraint. Then, the dispersion relation takes the form
If either $a$ or $b\rightarrow \infty$, we recover the slab dispersion relation (Lee & Lim Reference Lee and Lim2007)
It is recalled that the $k^2$ hidden in $\varDelta$ is $k^2=k_z^2+k_x^2$ in $\int \,{\rm d}k_x$-integral and $k^2=k_z^2+k_y^2$ in $\int \,{\rm d}k_y$-integral. Thus, let us change the integration variables, both $k_x$ and $k_y$, in (3.21) to $\kappa$
where $k^2=k_z^2+\kappa ^2$. In regard to the double signs in (3.23), the upper (lower) signs correspond to the symmetric (anti-symmetric) mode which also occurs in a slab plasma. For a square duct ($a=b)$, (3.23) reduces to the form identical to the slab dispersion equation (3.22), except for the factor $\sqrt 2$. This reduction is due to the $x$–$y$ symmetry. To recover the slab dispersion relation from (3.18), we take $k_y\rightarrow 0$, $b\rightarrow \infty$, and put $k_x=\lambda$. We can take $k_y\rightarrow 0$ since the $y$-direction has a translational invariance in a slab.
The duct dispersion relation in (3.23) can be contour integrated for a cold plasma, giving
where $\gamma =\sqrt {k_z^2-({\omega ^2-\omega _p^2})/{c^2}}$. For the anti-symmetric mode, the $\tanh$-function above is replaced by a $\coth$-function.
4. Discussion
In a bounded plasma, one way of solving the Vlasov equation by satisfying the specular reflection condition is to extend the plasma electric field in the manner of (1.3a,b). The job of Fourier transforming such a piecewise continuous periodic function, extending to infinity, is taxing. In this work, we present an alternative way of avoiding the hard labour 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 duct plasmas in a common work frame. Taking $b\rightarrow \infty$ in (3.24) gives
which is the slab ($0< x< a$) dispersion relation. Taking $a\rightarrow \infty$ in (4.1) gives
Equation (4.2) agrees with the semi-infinite dispersion relation obtained by Lee (Reference Lee2005) without introducing the fictitious surface charge sheet. If $u=0$, or $\omega '=\omega$, (4.2) agrees with the slab dispersion relation obtained from the fluid theory worked out by Gradov & Stenflo (Reference Gradov and Stenflo1983).
For a square duct, putting $a=b$ in (3.24) yields
which is similar to the slab dispersion relation. This is because the complete symmetry between the $x$ and $y$ coordinates makes the three-dimensional problem a two-dimensional problem practically.
The Doppler-shifted frequency $\omega '$ appearing in (3.24) represents the streaming effect. This is the first-order effect of the ratio $u/c$. The second-order effect enters through $U_{zz}$ in (2.11a–c), which produces terms of order $u^2/c^2$. If we chase those terms, we end up with replacing in (3.24) $(1-{\omega _p^2}/{\omega '^2}) \rightarrow 1-{\omega _p^2}/{\omega '^2}(1-{u^2}/{c^2})$. Relativistic treatment may be desirable to include the higher-order effect of $u^2/c^2$ in a systematic way.
This work may find applications in a laboratory or astrophysical situation where electromagnetic waves propagate through certain channels.
Acknowledgements
Editor P. Helander thanks the referees for their advice in evaluating this article.
Funding
This research received no specific grant from any funding agency, commercial or not-for-profit sectors.
Declaration of interests
The authors report no conflict of interest.