Subscribe Now Subscribe Today
Fulltext PDF
Research Article

Effect of Nonthermal Electron on Dust-acoustic Shock Waves in Dusty Plasma

Louis E. Akpabio and Akpan N. Ikot

A dusty plasma system containing non-thermal electron distributions, Boltzmann distributed ions and mobile charge fluctuating positive dust has been considered. The nonlinear propagation of the dust-acoustic (DA) waves in such dusty plasma has been investigated by employing the reductive perturbation method. The effect of non-thermal electrons on the height and thickness of DA shock waves are also studied. It has been found that the thickness of the Dust acoustic shock composition decreases as the non-thermal parameter increases, while the amplitude of the shock composition thickness varies with the charge fluctuating dust.

Related Articles in ASCI
Similar Articles in this Journal
Search in Google Scholar
View Citation
Report Citation

  How to cite this article:

Louis E. Akpabio and Akpan N. Ikot, 2012. Effect of Nonthermal Electron on Dust-acoustic Shock Waves in Dusty Plasma. Research Journal of Physics, 6: 41-49.

DOI: 10.3923/rjp.2012.41.49

Received: February 09, 2012; Accepted: June 13, 2012; Published: July 19, 2012


Dust and plasmas exist together in the universe and they make dusty plasmas. Dusty plasmas are found in cometary tails, asteroid zones, planetary ring, interstellar media, lower part of Earth’s ionosphere and magnetosphere, etc. (Goertz, 1989; Mende’s and Rosenberg, 1994; Shukla, 2001; Horanyi and Mendis, 1986; Verheest, 2001; Horanyi, 1996). Noticeable applications of dusty plasmas are also found in laboratory devices (Barkan et al., 1995; Merlino et al., 1998; Homann et al., 1997). There has been a rapidly growing interest in the field of dusty plasmas because of its great variety of new phenomena associated with wave and instabilities (Verheest, 1992; Pieper and Goree, 1996; Bliokh and Yaroshenko, 1985). The existing plasma wave spectra are not only modified by the presence of charged dust grains in a plasma (De-Angels et al., 1988; Shukla and Stenflo, 1992), it also brings about new novel eigen modes such as Dust Acoustic (DA) wave (Rao et al., 1990; Barkan et al., 1996), Dust Ion Acoustic (DIA) waves (Shukla and Silin, 1992; Barkan et al., 1996) etc.

From the first theoretical study on the ultra low frequency DA waves in dusty plasma by Rao et al. (1990) and motivated by the experimental observations of these waves (Barkan et al., 1996; Pieper and Goree, 1996), numerous investigations have been carried out to study the different aspects of the physics of dusty plasmas during the past few years. However, most of the investigations were mainly on dusty plasmas with negatively charged dust grains (Amin et al., 1998; Popel and Yu, 1995; Ma and Liu, 1997), in this regard, nonlinear solutions and double layers in dusty plasmas have been investigated by several authors (Bharuthram and Shukla, 1992; Mamun et al., 1996). However, in the space plasmas environments, some plasma systems are found with positively charged dust grains (Mendis and Horanyi, 1991; Chew et al., 1993; Haunes et al., 1996; Horanyi et al., 1993). Such dust grains with net positive charge are due to processes such as irradiation by Ultraviolet (UV) light; thermionic emission produced by radiative heating as well as secondary emission of electrons from the surface of the dust grains (Verheest, 1992; Shukla and Mamun, 2002).

Recently, Paul et al. (2009), investigated the nonlinear propagation of DA waves accounting for the charge fluctuating positive dust and Boltzmann-distributed electrons and ions. For this purpose, they derived the Burgers equation, by employing the reductive perturbation method (Washimi and Taniuti, 1996). They showed that, the dust charge fluctuation is a source of dissipation and is responsible for the formation of collisionless DA shock waves in such dusty plasma. Since in a real dusty plasmas; the electron behaviour can be powerfully modified by the nonlinear potential of the localized DA composition by generating a population of fast vigorous electrons, the present paper is mainly to determine how the electron non-thermality effect can be expected to modify the result of Paul et al. (2009). This simplification involves a little increase in algebraic intricacy of the pertinent formulas. This notwithstanding, the basic principles do not change.


We consider unmagnetized collisionless dusty plasma consisting of non-thermal electrons, Boltzmann-distributed ions and charge fluctuating positively charged mobile dust. We assume for simplicity that all the grains have the same charge, equal to qd = +zde, with zd representing the charge state of the dust component. Hence, charge neutrality at equilibrium is given by ne0 = ni0+zd0 nd0, where ne0(ni0) is the equilibrium electron (ion) number density, nd0 is the dust density at equilibrium, zd0 represent equilibrium charge state of the dust component. All the dust grain is assumed to be spheres of radius rd. The basic equations for one-dimensional DA waves for such a dusty plasma is given as:




where, φ is the electrostatic potential, nd, ne, ni are respectively, the number density for the plasma species for dust, electrons and ions, ud is the dust fluid speed. The non-thermal electron distribution is given as (Carins et al., 1995):




and the Boltzmann distributed ion as:


where, α1 is a parameter determining the number of non-thermal electrons present in our plasma model, kB is the Boltzmann constant and Te (Ti) is the electron (ion) temperature. Neglecting all other charging processes, we assume that the dust is charged by photoemission current (I+p), thermionic emission current (I+t) and electron absorption current (I¯e) only. The charge state zd of the dust component is not constant but varies according to the following equation (Paul et al., 2009; Shukla and Mamun, 2002):






where, h is the Planck’s constant, Tph is the photon temperature, we is the work function, J is the UV photon flux, Y is the yield of photons. The typical values of we, J and Y are given respectively as 8.2 eV, 5.0x104 photons /cm2/s and 0.1. For convenience, we express the set of Eq. 1 to 7 in normalized form by introducing the following normalized variables: Nd = nd/ndo, ud-ud/Cd, φ = eφ/kBTe, Zd = zd/zdo, X = x/λDd, T = twpd, λDd = (kBTe/4πz2do)1/2, Cd ≡ (zdokBTe/md)1/2 and wpd = (4πz2dondoe2/md)1/2 to obtain the following equations:





where, σ is Te/Ti; μi is nio/zdondo; μi is nio/zdondo; μ is π r2d/zdowpd; P is JY; Q is 2e-we/kBTe(2πmekBTe)/h2 3/2(8kBTe/πme)1/2, R is ne0(8kBTe/πme)1/2, α is zdoe2/rdkBTph; β is zdoe2/rdkBTe.


To study the dynamics of nonlinear dust acoustic shock waves in the presence of non-thermal electrons, Boltzmann distributed ions and charge fluctuating positive dust grains; we employ the reductive perturbation technique (Washimi and Taniuti, 1996). We introduce the stretched coordinates (Das et al., 1997) ξ = ε(X-V0T) and τ = ε2T, where ∈ is a small parameter and V0 is the DA shock waves velocity normalized by Cd. The variables Nd, Ud, Zd and φ are expanded as:


Now, substituting these expansions in to Eq. 11-14 and collecting the terms of different powers of∈, in the lowest order, we obtain:







The next order in ε,O(ε2) yields a system of equations that leads to Burgers equation as follows:





Making using of Eq. 16-23 we eliminate Nd2, Ud2, Zd2 and φ2 to obtain the following equation:


where, the nonlinear coefficient A and the dissipation coefficient B are given by:




The Burgers equation which describes the nonlinear propagation of the DA shock waves in the dusty plasma under consideration is given as Eq. 24. It can be observed that, the right hand side of Eq. 24 which represent the dissipative term is due to the presence of non thermal parameter (βn), the ratio of electron and ion temperature (σ) and the charge fluctuating positive dust (μi).


Our expression for Zd1 as Eq. 18 agrees with what is obtained by Paul et al. (2009) when nonthermal parameter βn is set to zero and σ set to 1. We strongly feel that the last term in the denominator for Eq. 18 should be [+μR(1+β)] as against what is obtained by Paul et al. (2009) as [-μR (1+β)]. Likewise, the last term for f as (-μRβ) reported by Paul et al. (2009) should be (+μRβ) as in our report for Γ. Equation 19 which gives the linear desperation relation for DA waves is greatly altered by the presence of the electron nonthermal parameter, ratio of electron and ion temperature, as well as the positive dust charge fluctuation. For stationary shock wave solution of Eq. 24, we set ζ = ξ-U0τ’ and τ’ = τ to obtain the equation:


The latter equation can be integrated, using the conduction that φ is bounded as wζ→ or by the application of Tanh method (Malfliet, 1992, 2004; Malfliet and Hereman, 1996a, b) to yield:




Equation 29 represents a monotonic shock-like solution with the shock speed, the shock height and the shock thickness given by U0, φ0 and Δsh, respectively. It is obvious from Eq. 29 that, the presence of electron non- thermal parameter significantly modifies the shock wave amplitude and its width.

To see the influence of the non-thermal parameter on the DA shock waves, we chose σ as 1.5 and vary βn. The following parameters; U0 = 0.1 m sec-1, P = 5.00x1017 m-2 sec-1, Q = 1.07x1031 m-2 sec-1, R = 2.48x1013 m-2 sec-1, V0 = 0.8, β = 1.2x10-9 c2 kg-1 m-3 sec-2, μ = 2.5x10-12 m2 sec-1 corresponding to the mesosphere event has been chosen from Paul et al. (2009). Figure 1 shows that as the βn increases; the positive shock width decreases, while the amplitude of the positive shock width varies as μi increases.

Fig. 1: Variation of positive shock thickness (Δsh) potential profile with non-thermal electron parameter (βn) for different values of μi

Fig. 2: Variation of amplitude φ0 of the positive shock wave with (βn) for different values of μi

Fig. 3: Variation of amplitude φ0 of the negative shock wave with βn for different values of μi

The variation of amplitude (φ0) of the positive shock height waves are presented in Fig. 2. From this Fig. 2, it can be observed that, the amplitude of the positive shock height decrease with increase in βn up to 0.5 and after this point it increases. Also, we can observe that, when μi the amplitude of the shock height wave, is smaller than what is obtained for μi = 0.8 due to the fact that the dust charge fluctuation is a sources of dissipation DA waves.

Figure 3 demonstrate the variation of negative amplitude of the shock height with βn. It shows that the shock height increases with βn up to 1.75 after which, there is no variation of the shock height with βn for μi = 0.8. When μi = 1.2, the amplitude of the negative shock height suddenly increases with βn from 1.78 and then drops to exactly the same amplitude on the negative axis. After which the negative amplitude of the shock height increase with βn to infinity. Since space plasma are more realistically modeled by making use of non-thermal velocity distributions as discussed by Maharaj et al. (2006), we have seen that the non-thermal electron distribution function significantly modifies the result obtained by Paul et al. (2009).


We have extended the recent work of Paul et al. (2009) to see under what conditions the electron non-thermality effect can be expected to modify the results of their analysis. We have shown here that, the basic feature; of the non-linear DA waves are modified by both presence of the non-thermal electron and the charge fluctuating dust in dusty plasmas. Present results are summarized as follows:

The width of the DA shock structures decreases as the non- thermal parameter increases, while the amplitude of the shock structures width varies as μi increases. This is due to the fact that dust charge fluctuation is a source of dissipation and lead to the development of DA shock waves in the dust plasma
It is also shown that, the positive amplitude of the shock height decreases with increase in βn up to a point and then increases. While, when the charge fluctuation (μi) is increased to 1.2, the negative shock height exhibits the occurrence of kink as in Fig. 3. The findings in this paper are important in understanding nonlinear DA wave phenomena in space plasmas
Amin, M.R., G.E. Morfill and P.K. Shukla, 1998. Modulational Instability of dust acoustic and dust-ion-acoustic waves. Phys. Rev. E, Stat. Phys. Plasmas Fluids Relat. Interdis Top, 58: 6517-6523.

Barkan, A., N. D'Angelo and R.L. Merlino, 1996. Experiments on ion-acoustic waves in dusty plasmas. Planet Space Sci., 44: 239-242.

Barkan, A., R.L. Merlino and N. D'Angelo, 1995. Laboratory observation of the dusty acoustic wave mode. Phys. Plasma, 2: 3563-3565.
CrossRef  |  Direct Link  |  

Bharuthram, R. and P.K. Shukla, 1992. Large Amplitude ion-acoustic solutions in a dusty plasma. Planet Space Sci., 40: 973-977.

Bliokh, P.V. and V.V. Yaroshenko, 1985. Electrostatic waves in Saturn's rins. Sov. Astron., 29: 330-336.
Direct Link  |  

Carins, R.A., A.A. Mamun, R. Bingham, R.O. Bostrom, C.M.C. Nairn and P.K. Shukla, 1995. Electrostatic solitary structures in non-thermal plasmas. Geophys. Res. Lett., 22: 2709-2712.
CrossRef  |  Direct Link  |  

Chew, V.W., D.A. Mendis and M. Rosenberg, 1993. Role of grain size and particle velocity distribution in secondary electron emission in space plasma. J. Geophys. Res., 98: 19065-19076.
CrossRef  |  Direct Link  |  

Das, G.C., C.B. Dwivedi, M. Talukdar and J. Sharma, 1997. Anew mathematical approach for shock-wave solution in a dusty plasma. Phys. Plasmas, 4: 4236-4239.
CrossRef  |  Direct Link  |  

De-Angels, U., V. Formisano and M. Giordano, 1988. Ion plasma waves in dusting plasmas. J. Plasma Phys., 40: 399-406.

Goertz, C.K., 1989. Dusty plasmas in the solar system. Rev. Geophys., 27: 271-292.
CrossRef  |  Direct Link  |  

Haunes, O., J. Troim, T. Blix, W. Moretensen, L.I. Naeshein, E. Thrane and T. Tonnesen, 1996. First detection of charged dust particles in the Earth's mesosphere. J. Geophy. Res., 10: 10,839-10,847.
CrossRef  |  Direct Link  |  

Homann, A., A. Melzer, S. Peters ans A. Piel, 1997. Determination of the dust screening length by loser-excited lattice waves. Phys. Rev. E. Staf. Phys. Plasmas Fluids Relate Intercrossing Top, 56: 7138-7141.

Horanyi, M. and D.A. Mendis, 1986. The dynamics of charged dust in the tail of comet giacobini-zinner. J. Geophys. Res., 91: 355-361.
CrossRef  |  Direct Link  |  

Horanyi, M., 1996. Charged dust dynamics in the solar system. Annu. Rev. Astron. Astrophys, 34: 383-418.
CrossRef  |  Direct Link  |  

Horanyi, M., G. Morfill and E. Griin, 1993. Mechanism for the acceleration and ejection of dust grain from Jupiter's Magnetosphere. Nature, 363: 144-146.
CrossRef  |  Direct Link  |  

Ma, J.X. and J. Liu, 1997. Dust-acoustic soliton in dusty plasma. Phys. Plasmas, 4: 253-255.
CrossRef  |  Direct Link  |  

Maharaj, S.K., S.R. Pillay, R. Bharuthram, R.V. Reddy, S.V. Singh and G.S. Lakhina, 2006. Arbitrary amplitude dust acoustic double layers in a non-thermal plasma. J. Plasma Phys., 72: 43-58.
CrossRef  |  Direct Link  |  

Malfliet, W. and W. Hereman, 1996. The tanh method: II. Perturbation technique for conservative systems. Phys. Ser., 54: 569-575.
CrossRef  |  Direct Link  |  

Malfliet, W. and W. Hereman, 1996. The Tanh method: I. Exact solutions of nonlinear evaluation and wave equations. Phys. Scr., 54: 563-568.

Malfliet, W., 1992. Solitary wave solution wave solutions of nonlinear wave equations. Am. J. Phys., 60: 650-654.
Direct Link  |  

Malfliet, W., 2004. The Tanh method: A tool for solving certain classes of nonlinear evolution and wave equations. J. Comp. Applied Math. 164-165: 529-541.
CrossRef  |  Direct Link  |  

Mamun, A.A., R.A. Cairns and P.K. Shukla, 1996. Solitary potentials in dusty plasmas. Phys. Plasmas, 3: 702-704.
CrossRef  |  Direct Link  |  

Mende's, D.A. and M. Rosenberg, 1994. Cosmic dusty plasm. Annu. Rev. Astron. Astrophys., 32: 419-463.
CrossRef  |  Direct Link  |  

Mendis, D.A. and M. Horanyi, 1991. Dust-Plasmas Interaction in the cometary environment. Geophys. Monogr. Ser., 61: 17-25.
Direct Link  |  

Merlino, R.L., A. Barkan, C. Thompson and N. D'Angela, 1998. Laboratory Studies of wave and instabilities industry plasmas. Phys. Plasmas, 5: 1067-1614.
CrossRef  |  Direct Link  |  

Paul, K.S., G. Mandal, A.A. Mamun and M.R. Amin, 2009. Dust-Acoustic Shock waves in a dusty plasma with charge fluctuating positive dust. IEEE Trans. Plasma Sci., 37: 627-631.
CrossRef  |  

Pieper, J.B. and J. Goree, 1996. Dispersion of plasma dust acoustic waves in the strong-coupling regime. Phys. Rev. Lett., 77: 3137-3140.
CrossRef  |  Direct Link  |  

Popel, S.I. and M.Y. Yu, 1995. Ion acoustic solution in impurity-containing plasmas. Contrib. Plasma Phys., 35: 103-108.

Rao, N.N., P.K. Shukla and M.Y. Yu, 1990. Dust-Acoustic waves in dusty plasmas. Planet. Space. Sci., 38: 543-546.
CrossRef  |  Direct Link  |  

Shukla, P.K. and A.A. Mamun, 2002. Introduction to Dusty Plasma Physics. lop publishing, Bristol, UK.

Shukla, P.K. and L. Stenflo, 1992. Stimulated scattering of electromagnetic waves in dusty plasma. Austrophys. Space Sci., 190: 23-32.
Direct Link  |  

Shukla, P.K. and V.P. Silin, 1992. Dust ion acoustic wave. Phys. Scr., Vol. 45. 10.1088/0031-8949/45/5/015

Shukla, P.K., 2001. A survey of dusty plasma physics. Phys. Plasmas, 8: 1791-1803.
CrossRef  |  Direct Link  |  

Verheest, F., 1992. Nonlinear dust-acoustic waves in multispecies dusty plasmas. Planet. Space Sci., 40: 1-6.
CrossRef  |  Direct Link  |  

Verheest, F., 2001. Waves in Dusty Space Plasmas. Kluwer Academic, Dordrecht, Netherlands..

Washimi, H. and T. Taniuti, 1996. Propagation of ion-acoustic solitary wave of small amplitude. Phys. Rev. Lett., 17: 996-998.
CrossRef  |  Direct Link  |  

©  2019 Science Alert. All Rights Reserved
Fulltext PDF References Abstract