Potential Surface Waves at The Vacuum-Radiative Collisional Plasma Interface
Potential Surface Waves at The Vacuum-Radiative Collisional Plasma Interface
Potential Surface Waves at The Vacuum-Radiative Collisional Plasma Interface
vðxÞ exp iðky y xtÞ, where vðxÞ is a constant, ky is the Bp ¼ Cp exp ðbp xÞ;
component of the wave vector along y-axis, and x is the per-
BVac ¼ CVac expðbVac xÞ; (18)
turbation frequency. After linearizing and taking the diver-
gence of Eq. (2) further using Eqs. (1) and (8), we obtain where
1 @ 2 dn 4pe2 dn 8 ar 4 5 n0 kB T0e r2 dn 1=2
¼ T þ x2
n0 @t2 m 9 q 0e 3 q n0 bVac ¼ ky2 2 : (19)
c
fcol @dn
: (9)
n0 @t The indices p and Vac are used for plasma medium and vac-
uum, respectively, while Cp and CVac are constants. We pro-
For our system d is for perturbed quantity and
ceed similarly for the electric field equations in plasma
" # medium and vacuum from the Maxwell equations
2 @2 @2 @ 2 dn
r dn ¼ þ 2 n0ðxÞ expðiky y ixtÞ ¼ dnky2 :
@x 2 @y @x2 @ 2 dE 2 4pn0 e @dV 1 @ 2 dE
þ k y dE ¼ 2 : (20)
(10) @x2 c2 @t c @t2
After linearizing and temporal Fourier transformation of
Substituting Eq. (10) into Eq. (9), the wave equation for the
Eq. (2), we get
density perturbation of electron becomes as
2 3
8 ar 4 5 n0 kB T0e
@ 2 dn 6 T þ 7
g2p dn ¼ 0; (11) i 4e 9 q 0e 3 q 5
@x2 dV ¼ dE þ rdn :
0 11=2 ðx þ ifcol Þ m n0
x2p x2 ixfcol (21)
gp ¼ Bky2 þ C ; (12)
@ 8 ar 4 5 n0 kB T0e A
T0e þ Substitution of Eq. (21) into Eq. (20) develops the elec-
9q 3 q
tric field equation of semi-bounded collisional radiative
where dn is the perturbed number density of electrons; it can plasma as
be noted that while deriving Eq. (11), very slow variations
@2 2 4peA expðgp xÞ gp x^ þ iky y^
@x2 ky are neglected. Eq. (11) has usual second order dif- dE ¼ Mp exp bp x
ferential equation’s solution x2p x2 ixfcol
8 ar 4 5 n0 kB T0e
dn ¼ C exp ðgp xÞ; (13) T þ ; (22)
9 q 0e 3 q
where C is a constant. Taking Curl of Eq. (2), using Maxwell
where rdn¼ @n @x þinky ¼Aexpð gp xÞ gp þiky . Similarly
Eqs. (3)–(6) and applying the sinusoidal characteristic of
for vacuum coupling of Eqs. (3) and (4) gives the second
B ¼ B0 ðxÞ exp½iky y ixt ; leads to
order differential equation and its solution
@2 x2 Bc r2 EVac b2Vac EVac ¼ 0; (23)
2
ky2 þ 2 ¼ r dV: (14)
@x c 4pn0 e
EVac ¼ MVac expðbVac xÞ; (24)
The coupling of Eqs. (2) and (14) gives
where MVac is a constant. At the interface, i.e., at x ¼ 0; we
@2 x2 Bðx þ ifcol Þ xB sustain only the part of solution that decay away from the
2
ky2 þ 2 ¼ 2 ; (15) interface in both regions. The parallel and perpendicular
@x c x2p c
components of the electric field and electric displacement
2
where x2p ¼ 4pnm0e e . After some straightforward mathematics, vector are continuous for both plasma and vacuum regions,
the wave equation for the magnetic field of the collisional- i.e., Epy ¼ EVacy and vpx ¼ vVacx ; where vpx and vVacx are
radiative plasma is the electric displacement vectors. Applying the boundary
condition upon Eqs. (22) and (24) yields the component
@2B separation
b2p B ¼ 0; (16)
@x2
MVacx expð bVac xÞ ¼ Mpx exp bp x
where 4pegp A expð gp xÞ
0 11=2 þ f2 ; (25)
x x2p x2 ixfcol
@ x2p 2
x A
x þ ifcol
bp ¼ ky2 þ : (17) MVacy expð bVac xÞ ¼ Mpy exp bp x
c2
4peAiky expð gp xÞ
The solution of the second order differential equation in the f2 ; (26)
plasma region and vacuum are x2p x2 ixfcol
112111-4 Khan, Rozina, and Jamil Phys. Plasmas 23, 112111 (2016)
here f2 ¼ 8 ar 4
T0e þ 53 n0 kqB T0e whereas on interface (x ¼ 0) plasmas and vacuum is first time derived. Applying the elec-
9q
trostatic limit c ! 1, and supercritical
dense
plasma condi-
above equations reduce to 8 ar 2 4
tion 9 q ky T0e þ 53 n0 kqB T0e x2p x2 ixfcol , dispersion
4pegp A equation reduces to
MVacx ¼ Mpx þ n2 ; (27)
x2p x2 ixfcol 2
2
x2p ky
4peAiky x þ ixfcol ¼ 4 1 þ qffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
MVacy ¼ Mpy n2 : (28) 2 x2 x2 ixf
p col
x2p x2 ixfcol
sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi#
Again from the Maxwell equations 8 ar 4 5 n0 kB T0e
T þ : (35)
9 q 0e 3 q
@Ey
iky ¼ a2Vac Ex : (29) It may be noted here that while deriving the above dis-
@x
persion equation, we have taken into account that bp ¼ ky
By substitution of Ex and Ey from Eq. (22) into Eq. (29) and bVac ¼ ky for electrostatic field case, whereas gp , bp,
and using Eq. (28) for the plasma region, we obtain and bVac are defined in Eqs. (12), (17), and (19), respec-
" # tively. Eq. (35) is the general equation that provides the dis-
4peAiky 2 persion relation. If the radiation pressure is ignored and
ap iky exp bp x MVacy þ 2 n
xp x2 ixfcol thermal pressure replaced by the Fermi pressure, the results
4peAky2 gp expð gp xÞ of Khorashadizadeh et al.19 may be retrieved. The results of
n2 Shahmansouri33 are retrieved if the radiation pressure and
x2p x2 ixfcol
the collisional effects are omitted after replacing thermal
4pegp Ab2Vac expð gp xÞ plasmas by the degenerate plasmas, while the results of
¼ b2P Mpx exp bp x n2 : (30) Lazar et al.34 may be obtained on ignoring collision fre-
x2p x2 ixfcol
quency in quantum plasmas. It may be noted that Kaw and
Similarly for the vacuum case McBride18 have studied the effect of density variation on the
! surface waves in half-space plasmas using full set of
ib 4pegp A Maxwell equations and have shown that the inclusion of den-
MVacy ¼ Vac Mpx þ 2 2
n : (31)
ky xp x2 ixfcol sity variations leads to the propagation of backward surface
waves in plasmas, i.e., surface waves with opposite group
Henceforth, Eq. (31) at the interface, i.e., x ¼ 0, becomes and phase velocities.
In order to obtain dispersion relation of surface waves
1 4peA on a radiative half space plasma, Eq. (35) is solved under the
Mpx ¼
bVac bp þ x 2 x2 ixf
b2Vac
p col condition of weakly coupled plasma, i.e., fxcol 1; hence, the
h i dispersion relation reduces to
ky2 bp þ bp bVac gp ky2 gp þ gp b2Vac n2 : (32) 2 pffiffiffi 3
i 2fcol sffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi
6 1þ ky 7
Applying the continuity condition for electromagnetic fcol xp 6 x 8 ar 4 5 n0 kB T0e 7
x ¼ i 6 pffiffiffi 41 þ pffiffiffi p T0e þ 5:
field components along with the boundary condition Vx 2 2 2xp 9q 3 q
(x-component of V) ¼ 0 at interface plane located at x ¼ 0
4peAg (36)
for plasma particles, Eq. (2) leads to Epx ¼ x2 p n2 , required 3
p ar T0e
to find x-component of Eq. (22), which in turn gives after Introducing the normalization coefficients as R ¼ ; n0 kB
some simple mathematics ky n0 kB T0e x fcol
K ¼ xp q ; W ¼ xp , and H ¼ xp ; the normalized distribu-
" # tion function reduces to
2 1 1
Mpx ¼ 4peAgp n 2 : (33) " pffiffiffi rffiffiffiffiffiffiffiffiffiffiffiffiffiffiffi#
x2p xp x2 ixfcol H 1 1 þ i 2H 5 8
W ¼ i 6 pffiffiffi 1 þ pffiffiffi K þ R : (37)
2 2 2 3 9
To the best of our knowledge, this dispersion relation The above equation shows that electromagnetic (thermal)
for surface waves at the interface of semi-bounded radiative radiations along with the total collisional rate of electron
112111-5 Khan, Rozina, and Jamil Phys. Plasmas 23, 112111 (2016)
FIG. 1. Imaginary frequency with change in temperature in space (1(A)) and laboratory plasma (1(B)).
causes radiative plasma surface unstable and leads to an Physically, rise in temperature of the larger electron number
instability if we have Im(W) > 0 for some particular values density contributes less to enhance the oscillation of electron
of wave vector k, or surface may become unstable if plasma frequency which in turn supports the comparatively
Im(W) < 0 for some other values of wave vector k. lesser instability growth rate of an electrostatic plasma wave
in the laboratory plasma, because the increasing temperature
IV. NUMERICAL RESULTS AND CONCLUSION in the bounded systems may enhance the collisional effects
For numerical analysis, Eqs. (38) and (39) are plotted to while space plasmas are extremely large dimensional natural
examine the effects of the radiation and collision frequency on system that supports higher temperature effects over the
the surface waves at the interface of collisional-radiative half growth rate of surface waves. Furthermore, Fig. 2 shows the
space plasma. The dispersion relation of electrostatic radiative growing surface wave behavior with the normalized collision
surface waves is plotted with typical parameters of space frequency H. These figures show that the growth rates of the
plasma35 n0 ¼ 105 cm3, ar ¼ 7:57 1015 K4 erg=cm3 ; T surface wave are increased by increasing the collision
¼ 103 K, whereas the parameters n0 ¼ 5:9 108 cm3, T ¼ 6 frequency. For the space plasmas (Fig. 2(A)) and laboratory
104 K, are for laboratory plasma36 in c.g.s system of (Fig. 2(B)), the normalized frequency probably affects the
measurements. growth rates in almost similar fashion. However, the labora-
In Fig. 1, we have numerically displayed the variations tory plasma facilitates the higher growth rate of interface sur-
of normalized dispersion frequencies with normalized wave face wave at higher collision frequency, while space plasmas
number at different temperature variations expressed through facilitate comparatively less growth of the plasmon signal.
the radiation parameter R, where R is the ratio of radiation to Fig. 3 emphasizes the behaviour of phase speed of surface
the thermal pressure. Further, it is evident from Fig. 1 that wave with varying temperature through the radiation parame-
with the increase in R, the growth rate of collisional radiative ter R for space plasma (Fig. 3(A)) and laboratory plasma(Fig.
surface waves increases. The left side, Fig. 1(A), describes 3(B)). For both the plasma systems, it is observed that the
the growth rate of the space vacuum interface, whereas the phase speed increases with increasing the thermal tempera-
right side, Fig. 1(B), shows the unstable behaviour of the sur- ture; however, the increasing rate of space plasmas is higher
face in laboratory plasmas. The surface waves in both the than that of the laboratory plasmas. It is clear from Figs. 1
plasmas show growing trends; however, in space plasmas the and 3 that the slope of the line increases with the radiation
growth rates are much higher than the laboratory plasmas. parameter R, both for laboratory and space plasmas. In other
FIG. 2. The normalized imaginary frequency instability against the normalized wave number with different collisions in space (2(A)) and laboratory (2(B)) plasmas.
112111-6 Khan, Rozina, and Jamil Phys. Plasmas 23, 112111 (2016)
FIG. 3. The normalized real frequency against the normalized wave vector with temperature variations in space (3(A)) and laboratory (3(B)) plasmas.
words, the propagation velocity increases by increasing the To conclude, we have studied the propagation of surface
radiation pressure. It is also clear from these figures that the waves in semi bounded half space radiative collisional
radiation pressure has a significant role for high wave number plasma. Following the fluid model, which includes the full
or for shorter wavelength and is less effective for the long set of Maxwell equations and applying particular boundary
wavelength case. Conclusively, it can be said that surface conditions, we have developed the dispersion relation of
waves show a similar trend in both laboratory and space plas- electrostatic surface waves in collisional radiative electron-
mas, i.e., both phase velocity and the growth rate increase as ion plasma and vacuum. The specific case of space plasma is
the parameter R increases. Fig. 4 shows the surface plot of considered to show that the electrostatic surface waves are
growth rate as a function of temperature. The temperature significantly affected by electromagnetic radiation pressure
range in space plasmas is 103 K–3 103 K (Fig. 4(A)). The and total collisional effects. The dispersion relation is plotted
least value of temperature contributes to the smallest range of numerically to show that the electromagnetic thermal radia-
growth rate. The range of growth rate increases on increasing tion along with the collisional frequency leads to an instabil-
the temperature. However, comparatively small increasing ity for some particular value of the normalized wave vector
trends are observed up to the medium range of the wave vec- k and normalized growth rate increases with the increase of
tor spectrum. The surface plot for laboratory plasmas is radiation parameter R collisional frequency fcol ¼ Hxp : On
shown in Fig. 4(B). The growth rate attains the highest value the other hand, the collisionless plasma radiation always sta-
at the maximum temperature effect. For medium spectrum bilizes the surface waves and helps in its propagations. Our
range of normalized wavevector (0–0.5), the growth rate result are of practical interest in order to understand the
observes small change on increasing temperature; however, dispersion properties of electron plasma oscillation at the
for the higher wavevector spectrum, the growth rate as well interface of radiative collisional plasma and vacuum for both
as the change in growth increases. At the least value of ther- laboratory and ionosphere plasmas. From the comparative
mal temperature, the growth rate attains smallest range while analysis of surface wave propagation in lab and ionosphere
at highest temperature the growth rate attains highest range of plasmas, it may be suggested that the bounded systems sup-
growth rate frequency. Figures 4(A) and 4(B) describe the port less instability growth of signals with increasing the
same competition of space and laboratory parameters. temperature through the surface wave mechanism.
FIG. 4. Surface plot as a function of temperature. (4(A)) Space parameters and (4(B)) laboratory plasma.
112111-7 Khan, Rozina, and Jamil Phys. Plasmas 23, 112111 (2016)
16
ACKNOWLEDGMENTS N. L. Tsintadze, G. Murtaza, and Z. Ehsan, Phys. Plasmas 13, 022103
(2006).
17
We would like to thank the anonymous referee(s) for R. H. Ritchie, Phys. Rev. 106, 874 (1957).
18
constructive suggestions to improve the manuscript. One of P. K. Kaw and J. B. McBride, Phys. Fluids 13, 1784 (1970).
19
S. M. Khorashadizadeh, S. Taheri Boroujeni, E. Rastbood, and A. R.
the authors (A. A. Khan) is thankful to HEC Pakistan for the Niknam, Phys. Plasmas 19, 032109 (2012).
partial financial support against the letter No. 2323/NRPU. 20
Ya. Zeldovich and Yu. Raizer, Physics of Shock Waves and High-
Temperature Hydrodynamic Phenomena, edited by W. D. Hayes and R. F.
1 Protstein (Academic Press, New York and London, 1966), Vol. 2.
S. A. Magnitskii, V. T. Platonenko, and A. V. Tarasishin, Laser Phys. 21
L. D. Landau and E. M. Lifshitz, Statistical Physics, Part 1 (Butterworths-
9(2), 522 (1999). Heinemann, Oxford, 1998).
2
J. Hubert, S. Bordeleau, K. C. Tran, S. Michaud, B. Milette, R. Sing, J. 22
S. S. Chandrasekhar, An Introduction to the Study of Stellar Structure
Jalbert, D. Boudreau, R. Moisan, and J. Margot, Fresenius, J. Anal. Chem. (Dover Publications, Inc., 1957).
355, 494 (1996). 23
J. M. Pitarke, V. M. Silkin, E. V. Chulkov, and P. M. Echenique, Rep.
3
A. Shivarova, T. Stoychev, and S. Russeva, J. Phys. D: Appl. Phys. 8, 383 Prog. Phys. 70, 1 (2007).
(1975). 24
W.-P. Hong and Y.-D. Jung, Phys. Lett. A 380(11), 1193 (2016).
4
Yu. M. Aliev and G. Brodin, Phys. Rev. A 42, 2374 (1990). 25
C.-G. Kim and Y.-D. Jung, J. Phys. D: Appl. Phys. 39, 119 (2006).
5
B. Buti, Advances in Space Plasma Physics (World Scientific, Singapore, 26
M. Ayub, H. A. Shah, M. N. S. Qureshi, and M. Salimullah, Commun.
1985), p. 167. Theor. Phys. 60(5), 623 (2013).
6
M. Y. Yu and I. Zhelyazkov, J. Plasma Phys. 20(2), 183 (1978). 27
H. J. Lee and M.-J. Lee, J. Mod. Phys. 7, 1400 (2016).
7
A. W. Trivelpiece and R. W. Gould, J. Appl. Phys. 30, 1784 (1959). 28
C. D. Wilson, R. K. Agarwal, and F. G. Tcheremissine, AIP Conf. Proc.
8
N. A. Azarenkov, A. N. Kondratenko, and K. N. Ostrikov, Radiophys. 1333, 146 (2011).
Quantum Electron. 36(5), 213 (1993). 29
M.-J. Lee and Y.-D. Jung, Phys. Plasmas 23, 104501 (2016).
9 30
A. Shivarova and I. Zhelyazkov, Phys. Plasmas 20, 1049 (1978). L. N. Tsintsadze, Phys. Plasmas 2, 4462 (1995).
10 31
M. Opher, L. O. Silva, D. E. Dauger, V. K. Decyk, and J. M. Dawson, L. N. Tsintsadze, D. K. Callebaut, and N. L. Tsintsadze, J. Plasma Phys.
Phys. Plasmas 8, 2454 (2001). 55, 407 (1996).
11 32
S. H. Glenzer, O. L. Landen, P. Neumayer, R. W. Lee, K. Widmann, S. W. H. M. Barlow and J. Brown, Radio Surface Waves (Clarendon Press,
Pollaine, R. J. Wallace, G. Gregori, A. Holl, T. Bornath, R. Thiele, V. Oxford, 1962).
33
Schwarz, W.-D. Kraeft, and R. Redmer, Phys. Rev. Lett. 98, 065002 (2007). M. Shahmansouri, Phys. Plasmas 22, 092106 (2015).
12 34
G. S. Selwyn, Jpn. J. Appl. Phys., Part 1 32, 3068 (1993). M. Lazar, P. K. Shukla, and A. Smolyakov, Phys. Plasmas 14, 124501
13
D. Mendis, Plasma Sources Sci. Technol. 11, A219 (2002). (2007).
14 35
V. N. Tsytovich, G. E. Morfill, and H. Thomas, Phys. Plasma Rep. 28, 675 W. Baumjohann and R. A. Treumann, Basic Space Plasma Physics
(2002). (Imperial College Press, London, 1996).
15 36
G. E. Morfill, V. N. Tsytovich, and H. Thomas, Phys. Plasma Rep. 29, M. A. Hellberg, R. L. Mace, R. J. Armstrong, and G. Karlstad, J. Plasma
2901 (2003). Phys. 64, 433 (2000).