Open access peer-reviewed chapter

Evolutions of Growing Waves in Complex Plasma Medium

By Sukhmander Singh

Submitted: May 10th 2020Reviewed: June 19th 2020Published: November 2nd 2020

DOI: 10.5772/intechopen.93232

Downloaded: 176


The purpose of this chapter to discuss the waves and turbulence (instabilities) supported by dusty plasma. Plasmas support many growing modes and instabilities. Wave phenomena are important in heating plasmas, instabilities, diagnostics, etc. Waves in dusty plasma are governed by the dynamics of electrons, ions and dust particles. Disturbances in solar wind, shocks and magnetospheres are the sources of generation of plasma waves. The strong interest in complex plasma provides us better understanding of physics of dusty universe, solar winds, shocks, magnetospheres, dust control in plasma processing units and surface modifications of materials. The theory of linearization of fluid equation for small oscillation has been introduced. The concept of fine particles in complex plasma and its importance is also explained. The expressions for the growth rate of the instabilities in turbulence plasma have been derived.


  • plasma oscillations
  • dispersion
  • turbulence
  • instabilities
  • dusty plasma
  • fine particles
  • Hall thrusters
  • resistive plasma
  • growth rate

1. Introduction to dusty plasma

The presence of fine particles of mass 10−10 to 10−15 kg and size 1 to 50 micrometer in an electron-ion plasma is called dusty plasma. Dusty plasma also termed complex plasma, plasma crystals, colloidal crystals, fine particle plasma, coulomb crystal or aerosol plasma and has been found in naturally in solar system, planetary rings, interplanetary space, interstellar medium, molecular clouds, circum-stellar clouds, comets, Earth’s environments, etc. Manmade plasmas are ordinary flames, dust in fusion devices, rocket exhaust, thermonuclear fusion, Hall thruster, atmospheric aerosols etc. [1, 2, 3, 4, 5, 6, 7]. The detail of existence of dusty plasma is given in Table 1. Earlier works shows that dust in plasmas has been considered as unwanted constituents and researchers had tried various methods to eliminate dust particles from plasma-processing units. For the moment, Positive aspects of dusty plasmas emerged, dust particles are playing various positive roles in plasma processing devices. The dust particles experience different forces in plasma. The Gravitational force, drag force, electromagnetic forces, polarization force and radiation pressure [1, 2, 3, 4, 5, 6, 7, 8].

Cosmic dusty plasmasDusty plasmas in the solar systemDusty plasmas on the earthMan-made dusty plasmas
Solar nebulaeCometary tails and comaeOrdinary flamesRocket exhaust
Planetary nebulaePlanetary ring Saturn’s ringsAtmospheric aerosolsDust on surfaces of space vehicle
Supernova shellsDust streams ejected from Jupitercharged snowMicroelectronic fabrication
Interplanetary mediumZodiacal lightlightning on volcanoesDust in fusion devices
Molecular cloudsCometary tails and comaeThermonuclear fireballs
Circumsolar ringsDust precipitators used to remove pollution from

Table 1.

Classification of dusty plasmas.


2. Physical processes in dusty plasma

There are many circumstances when astrophysical plasma and dust particles are found to coexist together. The study of dusty plasmas systems has an exciting properties which has attracted researchers over the world. These fine particles acquire some charges from the electrons to get charged. Moreover, in ordinary plasma, the charge considered to be constant on each particle, whereas, the charge on the dust particle varies with time and position [9, 10]. The charge on the dust particle generally depends on the type of dust grain, the surface properties of dust grain, the dust dynamics, the temperature, density of plasma and the wave motion in the medium. The plasma environments around these particles determine the nature of the charge (positive or negative) of these dusty plasmas. Although, most of the cases, these charged dust particles are negatively charged through different charging process. Their electric charge is determined by the size and composition of the grains [9, 10]. The fact that the frequencies associated with dust particles are smaller than those with electrons -ions and presences of fine particles modifies the dynamics of plasma motions and give rise to new types of propagating modes. For Dust acoustic waves, where ions and electrons are supposed to be inertia less pressure as gradient is balanced by the electric force, leading to Boltzmann electron and ion number density perturbations, whereas the mass of the dust play an important role in dust dynamics. In the dust acoustic wave the inertia is provided by the massive dust particles and the electrons and ions provide the restoring force. The effect of dust is to increase the phase velocity of the ion acoustic waves. This can be interpreted formally as an increase in the effective electron temperature which has important consequences for wave excitation. The dusty plasma also has a trend to oscillate at its plasma frequency [9, 10].

3. Parameters of dusty plasma

Dusty plasma and ordinary (electron-ion) plasma are different from each other due to the charge to mass ratio difference.

3.1 Dust plasma frequency

The electron and ion plasma frequency is much greater than the dust plasma frequency and it is defined as


3.2 Gyro frequency in dusty plasma

When charged dust grain/particle executes a spiral motion about the magnetic lines of force, then dust particle moves perpendicular to the magnetic field with Gyro frequency of the dust particle. The centrifugal force is balanced by the Lorentz force. In mathematically,


The radius of gyration is


where, Ωcd=eZBmd, is called the dust cyclotron frequency.

3.3 Macroscopic neutrality

The quasi-neutrality condition is obtained for the negatively charged dusty plasma by ne0=ni0+Znd0, here nd0is equilibrium dust particle density and Zis the electric charge number on the dust particles.

3.4 Strongly vs. weakly coupled dusty plasma (Coulomb correlation parameter)

The property of dust particles in the plasma is expressed by coupling parameter Γ. It is the ratio of the interparticle Coulomb potential energy to the thermal energy of the particles. When the value of Γexceeds unity, the species are termed to be strongly coupled otherwise weakly coupled dusty plasma. if rdis the average interparticle separation between particles, then coupling parameter


The interparticle separation can be found out by the relation rd=4πnd313. For the typical values of Ze=5000e, kBTd=0.05eV, nd=1010m−3, the coupling parameter comes out to be 1500. It is also experienced that, when Γ170, the dust particles are found in arranged fashion and said to be Coulomb crystals.

4. Applications of ordinary plasma

Plasma technology is safe, less costly and playing important roles in every fields of daily life. Some of these applications are discussed in Table 2.

TelecommunicationThe Global Positioning System (GPS) use ionosphere’s plasma layer to reflect the signal transmitted by GPS satellite for further communication usage.
SterilizationTo sterilize the surgical equipments, which are directly connected with patient’s immune system, where cleanliness is difficult
Medical treatmentPlasma treatment is contact-free, painless hardly damage tissue.
In dentistryPlasmas treatment are used inside the root canal to kill the bacteria
Pollution controllingPlasma technology is used to control gaseous and solid pollutions.
Water Purificationfor destroying viruses and bacteria in a water, Ozone (O3) generated by plasma technology is more effective and less costly at large scale than existing chlorination method
Etching and cleaning of materialsTo removes contaminants and thin layers of the substratum by bombarding with the plasma species which break the covalent bonds. It is also used to control the weight of the exposed substrate.
fusion researchPlasma is used to achieve high temperature to run the controlled thermonuclear fusion reactors
nanotechnologyPlasma discharges are helpful in growing the nanoparticles for nano world.

Table 2.

Applications of plasma in different fields.

4.1 Applications of dusty plasma

The presence of dust particles in a system also has positive impacts and has many applications in nanotechnology to synthesize the desired shape and size of the particles by controlling the dynamics of charged dust grains. Surface properties of the exposed materials could be improved by coating with plasma enhanced chemical vapor deposition method. Methane plasma is used to synthesize productive Carbon nanostructures which have like high hardness and chemical inertia. Dusty plasmas are also used for the fabrication of semiconductor chips, solar cells and flat panel displays.

5. Current status of the research

As we know that plasma support electrostatic as well as electromagnetic waves because of the motions of the charged particle. Studies of these waves provide the useful information about the state of the system. The resonance frequencies of plasmas waves can be used as diagnostics tool to characterize the plasma parameters. Plasma waves are generated for acceleration of energetic particles and heating plasmas. The exponential growing waves and modes in plasma removes the free energy from the system and permit the system to become unstable. The study of dusty plasma has gained interest in the last few decades due to its observations [1, 2, 3, 4, 5, 6, 7, 8, 9, 10] and applications in the space and laboratory [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. Many authors studied the linear and nonlinear electrostatic wave in the presence and absence of the external magnetic field [11, 12, 13]. Sharma and Sugawa studied the effect of ion beam in dusty plasma on ion cyclotron wave instability [13]. The presence of the charged dust grains in the plasma modifies the collective behavior of a plasma and excites the new modes [12, 13, 14].

The present charged dust particles introduces dust acoustic and dust ion acoustic waves in the plasma after altering the dynamics of electrostatic and electromagnetic waves of ordinary plasma [11, 12, 13, 14]. The charged dust grain also introduces growing and damped modes. Tribeche and Zerguini studied the dust ion-acoustic waves in collisional dusty plasma [15]. Rao et al. [16] predicted the existence of dust-acoustic wave in an unmagnetized plasma that has inertial dust and Maxwellian distributed electrons and ions. Shukla and Silin [17] showed the existence of dust-ion acoustic wave in a plasma. Barkan experimentally investigated that negatively charged dust grains enhances the growth rate of the electrostatic ion cyclotron instability [18]. Akhtar et al. [19] studied the dust-acoustic solitary waves in the presence of hot and cold dust grains. The existence of dust-acoustic wave and dust ion acoustic wave has been confirmed by many investigators in a laboratory experiments [18, 19, 20]. Ali [21] reported the electrostatic potential due to a test-charge particle in a positive dusty plasma. Bhukhari et al. derived generalized dielectric response function for twisted electrostatic waves in unmagnetized dusty plasmas [22]. Mendonça et al. showed that a modified Jeans instability lead to the formation of photonbubble in a dusty plasma which in turn form two different kinds of dust density perturbations [23]. Pandey and Vranjes predicted that growth rate of the instability is proportional to the whistler frequency in a magnetized dusty plasma [24].


6. Plasma model and basic equations

Phase and group velocity can be calculated by finding the relation between ωand k. This relation ω=ωk, is called the dispersion relation and contains all the physical parameters of the given medium in which wave propagates. If the frequency has an imaginary part, that indicates an instability. Plasma instability involves some growing modes, whose amplitude increases exponentially. In other words instability represents the ability of the plasma to escape from a configuration of fields [8, 9].

6.1 Electrostatic and electromagnetic waves in ordinary (electron-ion) plasma

Charged particles in a plasmas couples to electromagnetic field. Because of this effect various kinds of waves are formed in plasmas. Plasma waves are electrostatic or electromagnetic based on perturbed (oscillated) magnetic field. If there is a perturbed magnetic field (B10), plasma support electromagnetic waves. If the oscillating magnetic field associated with the wave is absent (B1=0), then only electrostatic waves are supported by plasma. In addition, Electrostatic waves may have longitudinal and transverse component depending on the direction of propagation with the perturbed electric field [8, 9].

A thermal unmagnetized plasma support many modes as discussed by Tonks and Langmuir in 1929. One is transverse waves in a plasma have dielectric constant εrω=εωε0=1ωpe2ω2. In case of a lower frequency wave (ω<ωpe), the dielectric constant would be negative. It turn out that if refractive index become imaginary, then waves cannot propagate but are damped (absorbed). Therefore plasma behaves like a waveguide with propagation and cut-off regions depending on the range of frequencies [8, 9].

The electron plasma wave and ion acoustic wave are the examples of electrostatic longitudinal modes, that is particle oscillate parallel to the direction of wave propagation. Ion acoustic waves are electrostatic waves, when both ions and electrons are allowed to oscillate in the wave-field. IA waves are low frequency longitudinal wave and we can use the plasma approximation, ne1ni1n0. The electron plasma wave (Langmuir mode) satisfy the dispersion relation ωk=ωpe2+3k2kBTeme, whereas the ion acoustic mode satisfy the dispersion relation ωk=k1+λDe2k2kBTemi, here λDeis the Debye length of electron [8, 9]. The propagation of electromagnetic waves in the unmagnetized plasma yield the dispersion relation ωk=ωpe2+k2c2.

7. Theoretical formulation for the studies of waves in dusty plasma

We consider a unmagnetized collisionless plasma consisting of electrons, ions and dust particles. Here we use the fluid equations and Maxwell’s equations to derive the dispersion relations in dusty plasma corresponding to ordinary electron –ion plasma. We denote υαand nαare the plasma velocity and density of the different species (α=e,i,d) having mass mα, temperature Tαin electron-volt. We write the equations of continuity and equation of motion of particles to derive the dispersion relation. Then the equations of motion governing the plasma can be written as


If we define thermal velocities V=Tαmα.

In the above equation, the derivative dυαdtis called the convective derivative. dυαdtcan be viewed as the time derivative of υαtaken in a “fluid” frame of reference moving with a velocity of υαrelative to a rest frame. υαtrepresents the rate of change of υαat a fixed point in space and υαυαrepresents the change of υαmeasured by an observer moving in the fluid frame into a region where υαis inhomogeneous.

7.1 Linearization of fluid equations

We consider the perturbed density nα1and velocity υα1indicated by subscript 1 along with their unperturbed density nα0and velocity uα0. The unperturbed electric field (magnetic field) as E0(B0) and the perturbed value of the electric field (magnetic field) is taken as E1(B1). To linearize all the equations, let us write nα=nα0+nα1, υα=υα1+υα0and E=E1+E0. If the amplitude is chosen to be much smaller than the wavelength of the instability, the equations of motion can be linearized. The perturbed quantities fα1are much smaller than their unperturbed values fα0, that is fα1<<fα0. If υα0and nα0are constant, the terms υα0nα0, nα0υα0and nα1υα0are equal to be zero. Further the terms υα1nα1, and nα1υα1are neglected as they are quadratic in perturbation. The linearized form of the fluid equations can be written as


Let us define δ=nd0ni0is the relative dust density, then quasi-neutrality condition follow


Thus, the assumptions of small oscillation give a set of linear equations.

7.2 Dust-acoustic waves (DAW)

The DAW is an electrostatic wave generated in dusty plasma, where inertia is provided by the dust grains. It is same to the ion-acoustic wave in general plasma, where inertia is provided by the ions. The frequency ωpdof dust acoustic wave is very low due to the high dust mass than the ion (electron) plasma frequency (ωpi,ωpe). That why dusty plasma supports low frequencies waves. In mathematically, ωpd=e2Z2nd0εomd<<ωpi,ωpe, where nd0is the equilibrium dust density.

Let us consider a situation, when dust density is get disturbed. This change will alter the charge on the dust particles and results to an enhancing negative space charge due to the process of negative dust charging. This total space charge density of dust ρd=eZnd0is shielded by the surrounding plasma ions and electrons. Therefore an electric field is generated due to the space charge by the fluctuations of dust charge density. This oscillating electric field imparts the force on the dust particle, which further pushes the fluctuations in the direction of the electric field and thus the wave propagates.

We consider that fluctuations are plane wave, which propagating inside the dusty plasma having the form f=foexpik·rωt. Them the time derivative /tcan be replaced by and the gradient by ik. Here f1nα1,υα1, E1, B1. The electrons and ions are assumed to inertia less as compared with mass of the dust grains and should have a Boltzmann distribution, namely


Here, ne0and ni0denote the unperturbed values of the electron and ion density respectively. Let us limit that, all the unperturbed velocities are zero, then equation of continuity and equations of motion follows.

The above three equations can be written as


Eqs. (13) and (14) gives


After substituting into Poisson’s equation, we obtain


The first term reduces to zero under the quasi-neutrality condition (ni0=ne0+Znd0). Let, the relative dust density is defined by δ=nd0/ni0, then we get


After simplification for the nontrivial solution, we readily obtain


the dispersion relation of the DAW shows depends on dust density, temperature of electron, ion and dust. It also shows depends on the inertia of electron and ion.

7.2.1 Limiting cases

In the limit of cold dust (Td=0) and cold ions (Ti<<Te). Then, the dispersion relation simplifies into the dispersion relation of ion-acoustic wave in an ordinary plasma


which is the same as for the ion-acoustic wave in classical plasma.

7.2.2 Behavior at low wave length

For small wave numbers k2λD,i2<<1the wave is acoustic ω=kCDAWwith the dust-acoustic wave speed


7.2.3 Behavior at high wave length

For large wave numbers k2λD,i2>>1, the wave is not propagating and just oscillates at the dust plasma frequency.

7.3 Dust ion acoustic wave (DIAW)

In the previous expression of ion-acoustic wave, the wave speed depends on ion temperature and on the mass of the dust. In the DIAW, the dust particles are supposed to immobile. We write the equation of motion, continuity and Poisson’s equation


The Poisson’s equation contains the perturbed electron and ion densities. The electrons are treated as Boltzmann distributed as follow


The only place, where the dust properties enter is the quasi-neutrality condition


Eqs. (24), (25) and (26) gives the dispersion relation for the DIAW as


Using Eq. (25)


It is clear from Eq. (27), that phase speed of the DIAW is increase as dust charge density increases. But the electron density has opposite effect on the speed of the DIAW.

8. Dissipative turbulence/instabilities in Hall thruster plasma

The section is devoted to the existing instabilities in a Hall thruster plasma. The principle of thrusters is the ionization of a Noble gas (propellant) in a crossed filed discharge channel. The accelerated heavy ions of inert gas are used to generate a thrust by the use of electrostatic forces. Xenon is used as an ion thruster propellant because of its low reactivity with the chamber and high molecular weight [25, 26, 27, 28, 29, 30, 31, 32, 33]. These types of devices support many waves and instabilities because of the turbulence nature of the plasma. These instabilities affect the performance and the efficiency of the device. In order to control these instabilities and further consequences, it has become necessary to study the growth rate of these instabilities. In a Hall thruster, the electrons experiences force along the azimuthal direction because of E×Bdrift. The collision momentum transfer frequency (v) between the electrons and neutral atoms are also taken into account to see the resistive effects in the plasma. Since ions do not feel magnetic field because of their larger larmor radius compared to length of the device. There equation of motion for ions can be written as


Motion of electrons under the electric and magnetic fields


8.1 Linearization of fluid equations

Let us denote the perturbed densities for ions and electrons as ni1and ne1velocities as υi1and υe1respectively. The unperturbed velocities υ0and u0are taken in the x- and y-direction respectively. The amplitude of oscillations of the perturbed densities are taken small enough. The linearized form of Eq. (29) and Eq. (30) are written as


The continuity equations of electrons and ions can be linearized as below


Fourier analysis: We seek the sinusoidal solution of the above equations, therefore the perturbed quantities are taken as f1f0expiωtikr. Them the time derivative /tcan be replaced by and the gradient by ik, here f1ni1,ne1, ϕ1, υi1, υe1, E1together with ωas the frequency of oscillations and kas the propagation vector.

By using Fourier analysis from Eq. (31)(34), the perturbed ion and electron densities are given as follows,


The expression for the electron density ne1is derived under the assumptions that Ω>>ω, kyu0and vin view of the oscillations observed in Hall thrusters.

8.2 Dispersion equation and growth rate of electrostatic oscillations

Finally, we use the expressions for the perturbed ion density ni1and electron density ne1in the Poisson’s equation ε02ϕ=ene1ni1in order to obtain


For the nontrivial solution of the above equation, the perturbed potential ϕ10, we have from Eq. (37)


This is the dispersion relation that governs the electrostatic waves in the Hall thruster’s channel. In the above equations, we introduced parameter ωei=e2n0mMε0, Ω=eB0mand Vth=YeTem.

After simplification of Eq. (38) we obtain


This is the dispersion equation that governs the electrostatic waves in the Hall thruster’s channel. It is clear from the above equation that Hall thruster support different waves and instabilities which satisfies the dispersion relation (39).

8.3 Results and discussion

To estimate the growth rates of the instability, we numerically solve Eq. (39) by giving typical values of all parameters used for the thruster [25, 26, 27, 28, 29, 30, 31, 32]. Therefore, for investigating the growths of the waves, we plot the negative imaginary parts of the complex roots (correspond to the instabilities) in Figures 13.

Figure 1.

Growth rate versus azimuthal wave number, with parameters of Hall plasma thrusters asB0100200G,n05×10171018m3,Te=1015eV,u0106m/s,v106/s andυ02×1045×104m/s.

Figure 2.

Growth rate versus collision frequency.

Figure 3.

Variation of growth rate with electron temperature.

It is found that the growth of the wave is enhanced for the larger values of the wave number of the oscillations as shown in Figure 1. In others words, the oscillations of larger wavelengths are stable. The findings are consistent to the results predicted by Kapulkin et al. [34]. In Figure 2, the variation of growth of the wave with the momentum transfer collision frequency is shown and it has been depicted that instability grow at faster rates in the presence of more electron collisions. This is mainly due to the resistive coupling with electrons’ drift’s in the presence of more collisions. The similar results are also predicted by Fernandez et al. [35] in the simulation studies of dissipative instability which shows that the growth of the instability is directly proportional to the square root of the collision frequency.

In Figure 3 it is noted that the wave grows faster if the electrons carry higher temperature. The present finding matched with the same results observed by other investigators [36, 37] of dissipative instability in an plasma.

8.4 Conclusion

In summary, we can say that dusty plasma physics has vital role in novel material processing and diagnostics tools. The nanoparticles of desired shape can be synthesized by controlling the dynamics of charged particles in the semiconductor industry. For example, the rotation of the dust particles can extract the electron flux in the magnetron sputtering unit. Thus, dusty plasma is a remarkable field in all areas of natural sciences. Though, the elimination of dust particles in the semiconductor industry is still a main alarm. The advanced development tools to learning dusty plasma will lead to additional discoveries about the astrophysical events and new scientific findings. The different dispersion relations are derived to see the behavior of the wave with wave number in complex plasma. The theory of linearization has been used for the smaller amplitude of the oscillations to derive the perturbed quantities. In the last section, the growth rate of dissipative instability has been depicted in a Hall thruster turbulence plasma. The instability grow faster with the collision frequency, azimuthal wave number and the electron temperature.


The University Grants Commission (UGC), New Delhi, India is thankfully acknowledged for providing the startup Grant (No. F. 30-356/2017/BSR).

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Sukhmander Singh (November 2nd 2020). Evolutions of Growing Waves in Complex Plasma Medium, Computational Overview of Fluid Structure Interaction, Khaled Ghaedi, Ahmed Alhusseny, Adel Nasser and Nabeel Al-Zurf, IntechOpen, DOI: 10.5772/intechopen.93232. Available from:

chapter statistics

176total chapter downloads

2Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Oseen’s Flow Past Axially Symmetric Bodies in Magneto Hydrodynamics

By Deepak Kumar Srivastava

Related Book

First chapter

Thermodynamic Properties of Ionic Liquids - Measurements and Predictions -

By Zhi-Cheng Tan, Urs Welz-Biermann, Pei-Fang Yan, Qing-Shan Liu and Da-Wei Fang

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us