Skip to main content

ORIGINAL RESEARCH article

Front. Astron. Space Sci., 11 October 2022
Sec. Low-Temperature Plasma Physics
This article is part of the Research Topic Wakefield Generation and Particle Acceleration in High-Intensity Laser Plasma and Beam-Plasma Interactions View all 5 articles

Test charge driven response of a dusty plasma with polarization force

  • 1National Centre for Physics at Quaid-e-Azam University Campus, CoE Physics, Islamabad, Pakistan
  • 2Department of Mathematics, College of Arts and Sciences, Khalifa University of Science and Technology, Abu Dhabi, United Arab Emirates
  • 3Department of Physics, COMSATS University, Islamabad Campus, Islamabad, Pakistan

By using the framework of kinetic model, the test charge driven response of a dusty plasma is evaluated in the presence of polarization force. The plasma containing electrons, singly ionized positive ions, and highly charged negative dust particulates is known as a dusty plasma, which can be perturbed by a test particle of charge qT moving with velocity vT along the z-axis. The polarization force purely comes from the high charging state of micron-sized dust particles, causing a deformation of shielding length due to density gradient and modifying the collective motion and particulates’ acceleration. The application of Fourier transformation technique to the set of dynamical equations may result in the shielded potential for a test charge with a modified dielectric constant. Several possibilities have been explored to evaluate the shielded potentials by imposing different conditions on the test charge speed in comparison with the thermal speeds of plasma species. It is found that the profiles of wakefield, Debye-Hückel and farfield (FF) potentials are strongly modified by the polarization force coefficient via the dust charging state and dust concentration. The present findings are useful in the study of strongly coupled dusty plasma, where micron-sized negatively charged dust grains are characterized by a high charging state.

1 Introduction

Most of the solid matter in the solar system indicates the presence of dust particles. In the gaseous form, the matter is ionized and dust particles coexist with it to give rise to the formation of dusty plasmas (Goertz, 1989). The dust particles are at least 10 to 12 orders of magnitude more massive than the ions (Goertz, 1989) and therefore, a unique situation arises that leads to the emergence of new spatio-temporal scales and hence new normal modes in the system. Dust particles are originally neutral, while the fragile and slippery electrons (whose electron thermal velocity is comparatively larger than the ion thermal velocity) stick onto the surface of dust to acquire a negative charge. However, there are instances where positively charged dust particles are also reported, for instance, on Jupiter (Horányi et al., 2004). Many new phenomena are associated with waves and instabilities in astrophysical and space environments due to the inclusion of dust components (Rao et al., 1990; Shukla, 2001). The most fundamental normal modes in a dusty plasma are the dust-acoustic (DA) (Rao et al., 1990) and dust-ion-acoustic (DIA) (Shukla and Silin, 1992) modes, which have been observed in laboratory experiments (Barkan et al., 1995; Barkan et al., 1996). Several distinct features of the linear waves have been confirmed later by analytical and experimental analyses (Shukla and Mamun, 2002) and new characteristics for nonlinear waves were explained by using various plasma compositions including the nonthermal character and magnetic field effects in dusty plasmas (Barkan et al., 1995; Barkan et al., 1996; Shukla and Mamun, 2002; Shukla and Mamun, 2003; Masood et al., 2010a; Masood et al., 2010b; Masood and Mirza, 2010; Masood et al., 2012; Sabeen et al., 2017; Nawaz et al., 2022; Shohaib et al., 2022).

The electrostatic potential in the close proximity of a free electron is the single-particle Coulomb potential (Clemmow and Dougherty, 1969; Chen, 1974). In case of plasmas, there are a number of moving charges that contribute to the field at any point and, therefore, the combined potential field at a certain location for a particular electron rapidly changes with time. It is extremely difficult to give an exact description of the fluctuating microfield in the vicinity of any given electron, however, the averaged, say over a time, value of the microfield and a volume adequate to contain a large number of particles, the fluctuating field can become sizeable and of practical interest (Clemmow and Dougherty, 1969). It is well-known that the effect of a single electron and its associated field in the presence of many other electrons is heavily shielded in the plasma. Such shielding is termed Debye shielding and it wards off the effect of an externally applied potential much quicker in the plasma (Neufeld and Ritchie, 1955; Joyce and Montgomery, 1967; Montgomery et al., 1968; Clemmow and Dougherty, 1969; Cooper, 1969; Kan, 1971; Laingal et al., 1971; Stenflo et al., 1973; Yu et al., 1973; Chen, 1974). This effect is intimately related to the effect of charged grains in dusty plasmas (Chen et al., 2001; Stoffels et al., 2001; Denysenko et al., 2005), in satellites and re-entry vehicles in space and ionosphere (Al’Pert et al., 1963; Liu, 1965; Gurevich et al., 1969).

It is to be noted that the expression for the Debye shielding due to a charged particle was obtained by considering it as static charge and infinitesimally small. In such a situation, the ambient plasma often acts as a Maxwellian plasma having the plasma species in thermal equilibrium. In a variety of physical situations, there could be a partial breakdown of Debye shielding for test charges, for instance, when the charge is moving or the plasma is inhomogeneous, anisotropic, collisional, magnetized, etc. (Neufeld and Ritchie, 1955; Joyce and Montgomery, 1967; Montgomery et al., 1968; Cooper, 1969; Kan, 1971; Laingal et al., 1971; Stenflo et al., 1973; Yu et al., 1973). Stenflo and Yu (Stenflo et al., 1973) utilized the Bhatnager-Gross-Krook (BGK) model (Bhatnagar et al., 1954) to investigate the electrostatic shielding of a slowly moving test charge in a collisional plasma and showed a parametric regime in which the Debye shielding could break down and lead to the possibility of completely unshielded charge along the direction of its motion. Such an occurrence could seriously affect the theories of charged-particle collisions that rely on the Debye cutoff in the interaction range. In other cases, the resonant interaction of the test charge with a plasma wave gives rise to an oscillatory wakefield behind the test particle. The wake potential was initially introduced by Nambu and Akama (Nambu and Akama, 1985). Since then, many investigations have been made to explain the wakefield potential distribution in different scenarios (Lemons et al., 2000; Ali, 2009; Ali, 2016). The dielectric response function and modified potential distribution around a test charge have also been studied in strongly coupled unmagnetized and magnetized dusty plasmas (Shahmansouri et al., 2017; Shahmansouri and Khodabakhshi, 2018; Shahmansouri, 2019).

The polarization force for dust particles arises on account of the deformation of the Debye sheath for density gradient. The polarization force is negligibly small for electrons and ions but becomes significant for dust species due to the higher charging state in complex plasmas. The polarization force on the dust particles can be written (Hamaguchi and Farouki, 1994a; Hamaguchi and Farouki, 1994b) as Fp=qd22λDλD2, where qd=zd0e is the charge of the dust particles, λD=λDiλDe/(λDe2+λDi2)1/2 is the effective Debye length due to the electron and ion species with electron and ion Debye radii λDe=Te/4πe2ne1/2, and λDi=Ti/4πe2ni1/2, respectively. The electron and ion temperatures are represented by Te and Ti (in energy units), whereas ne and ni refer to the electron and ion densities. If there are no gradients in temperature and the plasma is assumed to be weakly nonuniform in density, then λD=λD2(ni+TiTene)/(ni+TiTene). The latter can be simplified by considering the electron and ion species as inertialess. Hence, the polarization force becomes Fp=qdRE with its coefficient R=14zd0e2λDTi1TiTe(1zd0η), where η=nd0/ni0 represents the dust concentration and zd0 the dust charging state. Note that this force only accounts for high dust charging states or dust charges.

It is worth mentioning that polarization force acts on the dust particles when a low-frequency wave propagates and plasma background turns into a nonuniform plasma locally. This force always acts opposite to an electric field, regardless of the sign of charged dust particles. The polarization effect was first studied by Hamaguchi and Farouqi (Hamaguchi and Farouki, 1994a; Hamaguchi and Farouki, 1994b) and later, the idea was implemented by Kharapak et al. (Khrapak et al., 2009) to investigate the linear properties of DA waves in dusty plasmas. They showed that the phase speed of the DA waves decreases as long as the dust polarization force coefficient increases. Since then, numerous investigations have been carried out to examine the effect of polarization force on the profiles of linear and nonlinear wave propagation in polarized dusty plasmas (Mamun et al., 2010; Bandyopadhyay et al., 2012; Singh et al., 2018a). Recently, the modification of polarization force has been carried out with non-Maxwellian hybrid distribution (Singh et al., 2018b; Mehdipoor, 2022) and Tsallis statistics (Bentabet et al., 2017) for applications in dusty plasmas. Furthermore, the dispersive properties of dust acoustic waves have also been analyzed with charge-fluctuating dust grains (Shahmansouri and Mamun, 2016). In the case of degenerate plasmas, the magnitude of polarization force is shown to be enhanced and the DA wave frequency is exhibited to be reduced as compared to classical plasmas (Shahmansouri and Misra, 2019).

The layout of this manuscript is as follows: In Section 2, we adopt a kinetic approach to evaluate the electrostatic potential distribution around a test charge, moving with a constant speed in a (polarized) dusty plasma. On account of the dust polarization effect, the dielectric constant of the DA wave becomes significantly modified. Several limiting cases are derived for shielded potentials by imposing certain conditions on the test charge speed in comparison with the thermal speeds of plasma species. The profiles of the Debye–Hückel and wakefield and farfield potentials are analyzed for varying the polarization force coefficient and dust concentration. Section 3 presents numerical findings and discussion while Section 4 summarizes the main results and concludes this manuscript.

2 Kinetic model and shielded potentials

We consider a collisionless, unmagnetized multicomponent plasma containing electrons, singly ionized positive ions, and negatively charged dust particulates. On a dust timescale, an extremely low-frequency DA wave is produced due to dust motion in the presence of inertialess electrons and ions. The dusty plasma is globally a quasineutral gas of charged particles which holds an equilibrium charge-neutrality condition of the form ni0 = ne0 + zd0nd0, where ne0, ni0, nd0 and zd0 denote the equilibrium electron density, equilibrium ion density, equilibrium dust density and dust charging state, respectively. Dust particles attain a negative charge on their surface due to the impingement of ambient electrons and ions, resulting in the charge variation. However, for simplicity, the charge on the dust surface is treated as a constant (fixed) quantity in the limit when the dust charging period is much smaller than the dust plasma period.

The test charge-driven response of a (polarized) dusty plasma can be governed by the following linearized set of coupled Vlasov–Poisson equations (Krall and Trivelpiece, 1973; Nasim, 1999):

fd1t+vfd1r+zd0e1Rmdϕ1rfd0v=0,(1)

and

2ϕ1=4πene1+ni1zd0fd1dv4πqTδrvTt,(2)

where

ne1ne0eϕ1Te,  and   ni1ni0eϕ1Ti,(3)

where md and e, are the customary notations for dust mass and electronic charge, respectively. The perturbed (equilibrium) dust distribution function is denoted by fd1fd0 and perturbed electron (ion) number density by ne1 (ni1). It is to be noted that electrostatic and polarization forces are added up to express the total force as F=zd0e1Rϕ1/r acting on the dust particles with induced potential ϕ1. Dust polarization coefficient R(<1) not only modifies the Vlasov equation but also its associated force acts opposite to the electrostatic force and contributes to the dust susceptibility term for high dust charging state. For R>1, the dust polarization force dominates over the electrostatic force causing a null net restoring force on dust particles and leads to unstable oscillations. In the limit TeTi, the polarization force coefficient may reduce to R=14zd0e2λDTi in cgs units. In the present investigation, we shall focus our attention on the situation R<1. On the other hand, the first term in Eq. 2 indicates the plasma species density and the last term corresponds to the test charge density. A 3D Dirac delta function is denoted by the symbol δ and an observation point by r.

For the ES potential caused by a test charge, we now apply the space − time Fourier transformation to Eqs 13 and obtain the transformed quantities into ωk space, as

f̃d1k,ω=zd0e1Rmdωkvϕ̃1k,ωkfd0v,(4)
k2ϕ̃1k,ω=4πeñi1k,ωñe1k,ωzd0f̃d1k,ωdv+4πqTexpikvTωtdt,(5)
ñe1k,ωne0eϕ̃1k,ωTe,(6)

and

ñi1k,ωni0eϕ̃1k,ωTi,(7)

where ω and k pinpoint the wave frequency and wave number. Solving Eqs 47 together, we perform time integration to eventually come up with a transformed potential in this form

k2Dk,ωϕ̃1k,ω=8π2qTδωkvT.(8)

For instance, if qT = 0 in Eq. 8, then we immediately get the dispersion relation for DA waves from Dk,ω=0, because ϕ̃1k,ω0. See that dielectric constant Dk,ω is now modified by the dust-polarization effect, which appears in the dust-susceptibility term, can be given by:

Dk,ω=1+1k2λD2+1Rωpd2k2kωkufd0uudu,(9)

where fd0(u)=[1/2π1/2vTd]expu2/2vTd2 is the reduced 1D Maxwellian dust distribution function with dust-thermal speed vTd[=Td/md1/2]. Eq. 9 shows the presence of inertialess electrons and ions besides the mobile dust species. Expressing the dust-susceptibility term into the well-known dust dispersion function (Fried and Conte, 1961), one obtains from Eq. 9 as

Dk,ω=1+1k2λD2+1Rk2λDd2WCd,(10)

where λDd=Td/4πe2zd02nd01/2 stands for the dust Debye length and WCd=2π1/2qexpq2/qCddq the dust dispersion function with its argument Cd=ω/kvTd. The expansion of W − function for Cd > 1 exactly recovers the previous result (Bandyopadhyay et al., 2012) from (10).

To proceed further, we simplify Eq. 8 by applying the inverse Fourier transformation and integrating over ω, to eventually arrive at the standard result (Krall and Trivelpiece, 1973)

ϕ1r,t=qT2π2expikrvTtk2Dk,kvTdk.(11)

This is the ES potential of a test charge, which is moving along the z-axis with a constant speed vT in the polarized dusty plasma. The ES potential shall be solved either in spherical polar or cylindrical coordinate systems (depending upon the situation and mathematical tractability). The choice of the coordinate system may rely on the motion of the test charge; if the test charge is at rest or slowly moving, it exhibits spherically symmetric potential distribution around it. However, if the test charge moves with a finite speed in a certain direction, the potential distribution appears behind the test charge and may be appropriately solved by using cylindrical coordinates. Some limiting cases are to be discussed in the upcoming subsections by imposing certain conditions on the test charge speed in regard to the plasma thermal and phase speeds, for example (A) vTω/k, (B) vTvTd (C) vTvTe, (D) vT = 0, and (E) vTdω/kvTe, vTi.

2.1 Potential distribution around the resonating test charge

In this subsection, we shall decompose the total ES potential into two parts to investigate the Debye-Hückel and wakefield potential distributions around the test charge. Thus, we first simplify the plasma dispersion function [given in Eq. 10] in the limit Cd > 1 to immediately obtain the dielectric constant (Bandyopadhyay et al., 2012) in the form D(k,ω)=1+1k2λD21Rωpd2ω2 with its reciprocal form (Ali et al., 2003)

D1=k2λ02A+k2λ021+ωD2ω2ωD2,(12)

and the modified DA frequency

ωD=1Rωpdkλ0k2λ02+A1/2.

Observe that effective Debye length can be expressed in the usual form of the Debye length for an electron-ion plasma, as

λD=λDe2+λDi21/2λ0A,(13)

where λ0=Te/4πe2n01/2, A=ne0n0+TeTi with n0 = ni0ne0 + zd0nd0. It is worth mentioning that dust polarization force coefficient is always taken less than unity (viz., R<1) for micron-sized charged particulates and tends to unity for large-sized charged particulates. The variation of dust polarization force coefficient gives rise to a reduction of the phase speed of DA waves as compared to the conventional DA wave speed. Substituting Eq. 12 into (11), the Debye-Hückel (DH) and wakefield (WF) potentials can be expressed separately, as

ϕDHr,t=qT2π2λ02A+k2λ02expikrvTtdk,(14)

and

ϕWFr,t=qT2π2λ02A+k2λ02ωD2ω2ωD2expikrvTtdk.(15)

Now solving the DH potential first, we need to express the wave, position, and velocity vectors into spherical polar coordinates (Yaqoob et al., 2002; Ali et al., 2003), respectively, as k=ksinθkcosϕk,ksinθksinϕk,kcosθk, r=rsinθrcosϕr,rsinθrsinϕr,rcosθr, and vT=0,0,vT. Hence one has krvTt=kXcosϕk+kYsinϕk+kμZ, where X=ρ1μ21/2cosϕr, Y=ρ1μ21/2sinϕr and Z = ξzvTt represent the separation of the test charge from an observation point r in a moving frame. Here μ = cos θk, sinθk=1μ21/2, ρ = r sin θr, and z = r cos θr. After some more simplifications, Eq. 14 is finally reduced to the following result (Ali et al., 2003)

ϕDH=qTrexprλD,(16)

This is the DH potential for the test charge in a new moving frame having the radial and axial distances ρ=rsinθr andξ=zvTt, respectively, and r[=ρ2+ξ21/2] refers to the distance of the test charge from an observer. Notice that DH potential decreases exponentially with distance r and has no influence on the dust polarization effect because (14) does not rely on the dust dynamics.

To solve Eq. 15, we utilize cylindrical coordinates and adopt the same standard procedures (Nambu et al., 1995; Vladimirov and Nambu, 1995; Salimullah and Nambu, 2000; Nambu et al., 2001), to eventually obtain the main contribution to the WF potential, as

ϕWFρ=0,ξ2qTξA1+1Rωpd2λ02AvT211Rωpd2λ02AvT21×cos1RωpdξvTA.(17)

This is the WF potential of a test charge moving with a uniform speed vT along the z-axis in the polarized dusty plasma. The WF potential also leads to the existence of an attractive oscillatory potential (Nambu et al., 1995; Vladimirov and Nambu, 1995; Salimullah and Nambu, 2000; Nambu et al., 2001) for holding the limits vT>CD(1Rωpdλ0/A) and cos1Rωpdξ/vTA<0, where ξ=rvTt is the axial distance from the test charge. It is now clear that attractive WF potential is significantly influenced by the dust polarization force coefficient R and could have important consequences on the attraction forces between the same polarity charges to forming the ordered crystalline structures in (polarized) dusty plasmas (Fortov et al., 1996). The numerical analysis also confirms that attractive WF potential dominates over the repulsive DH potential as the DH potential fastly decreases beyond the shielding length or cloud.

2.2 Potential distribution around the slow test charge

Here we calculate the potential distributions around the slow test charge qT having speed vT in comparison with the dust thermal speed vTd in a polarized dusty plasma. Thus, assuming that vTvTd, which further implies that vTvTe,Ti as well. Consequently, the test charge is shielded by all the plasma species and there will be no resonance and no wakefields in this case. For slow test charge, the dielectric function takes the following form:

D=1+1k2λDe2+1k2λDi2+1Rk2λDd21+iπ2CdexpCd22.(18)

One observes that if the test charge is static (vT = 0), then the Landau damping term gets vanished and ultimately the DH potential (similar to Eq. 16 is obtained but with a modified shielding length as (Lakshmi et al., 1993)

λDef=λDe2+λDi2+1RλDd21/2

or

λDef=λ0/A+1RTeTdzd02η1/2(19)

Moreover, it was examined (Montgomery et al., 1968) later that Landau damping term cannot be rigorously neglected, especially in the situation when test charge speed is much slower compared to the dust thermal speed. Thus, keeping the Landau damping term, we simplify Eq. 18 to obtain its reciprocal form (Ali and Eliasson, 2017), as

D1=k2λDef2k2λDef2+1iπ21RμvTvTdλDd2k2λDef4k2λDef2+12.(20)

Here, μ is the angle between the wave number (k) and the test charge velocity vT. On substituting Eq. 20 into (11) and decomposing the total ES potential into two parts, as ϕ1 = ϕDH + ϕFF, where ϕDH=qT/rexpr/λDef is the usual short-range DH potential (Debye et al., 1923) with a modified shielding length λDef. Recalling Eqs 13, 19, it is found that there is only a difference in the finite dust temperature and polarization coefficient. By adopting the standard technique (Ali and Eliasson, 2017), the second part of the ES potential can be solved, as

ϕFF=4qT2πr1RvTλDefvTdλDdξλDdλDef3r3.(21)

This is the farfield (FF) potential caused by a slow test charge, derived under the limit rλDef. The FF potential decays as the inverse cube of the axial and radial distances from the test charge is strongly affected by the polarization coefficient R and shielding length λDef.

2.3 Potential distribution around the rapidly moving test charge

Here, we assume that if the test charge is moving very fast as compared to the electron thermal speed (viz., vTvTe) in the polarized dusty plasma. This would further imply that vTvTi and vTvTd. Hence we need to consider all the plasma species (i.e., the electrons, ions, and dust particulates) as dynamical in the dielectric constant (9) for holding the large velocity limit and so it becomes as

D11k2λDe22πvTekkvTuexpu22vTe2du1k2λDi22πvTikkvTuexpu22vTi2du1Rk2λDd22πvTdkkvTuexpu22vTd2du.(22)

In Eq. 22, we have only collected up to first order terms in the expansions from (9). Using the fact uexpβu2du=0 in (22), the dielectric constant eventually reduces to Dk,kvT1 and the corresponding potential yields the Coulomb potential

ϕ1ρ,ξ=ϕcρ,ξ,tqTr,(23)

with r=ρ2+ξ21/2. This confirms that if the test charge is moving very fast as compared to all the surrounding plasma species then it would not be screened anymore by the plasma species and hence no shielding effect in the dusty plasma.

2.4 Potential distribution around the static test charge

For a static test charge, we assume that vT = 0 in the expression of the total ES potential (11), which turns out to be

ϕ1r,0=qT2π2expikrk2Dk,0dk,(24)

with

D1+1k2λDe22πvTeexpu22vTe2du+1k2λDi22πvTiexpu22vTi2du+1Rk2λDd22πvTdexpu22vTd2du.(25)

In obtaining (25), we have assumed that electrons and ions are also dynamical besides dust particles in Eq. 9. After solving the above integrals using uexpβu2du=2π/β, we immediately come up with the result Dk,0=k2λDef2+1/k2λDef2, where the modified effective shielding length λDef is already defined in Eq. 19. Thus Eq. 24 finally turns out to be

ϕ1ρ,z,0=ϕDHQtrexprλDef.(26)

This confirms that for a null test charge speed, the test particle acts as static, and is screened by all the plasma species in the polarized dusty plasma.

2.5 Potential distribution around an intermediately moving test charge

For an intermediate velocity regime, like vTdvTvTe, vTi, the dielectric function (9) can be expressed in this form

D=1+1k2λD21Rk2λDd22πvTdkkvTuexpu22vTd2du.(27)

Solving the above integral, we get the dielectric function, as

D=1+1k2λD21Rk2λDd22πvTdkkvT0=1+1k2λD2,(28)

where λD is already given in (13). The dust particulates act as the inertial term, which contains dust polarization force coefficient R and gets vanished because only collecting the first order term in its expansion. Therefore, the velocity-dependent contribution does not appear in Eq. 28. As a consequence, the test charge potential finally reduces to the DH potential (Debye et al., 1923)

ϕ1ρ,ξ,t=ϕDHqTrexprλD,(29)

where r=ρ2+ξ21/2. Observe that only electrons and ions take part in the shielding process when the test charge is moving slower than the electron and ion species but faster than the mobile dust species.

3 Results and discussion

For numerical analysis, we first normalize the resultant Eqs 17, 19, and 21 by using the scaling parameters such as λ̄Def=λ̄Def/λ0, ρ̄=ρ/λ0, ξ̄=ξ/λ0, ϕ̄WF=ϕWF/(qT/λ0), ϕ̄DH=ϕDH/(qT/λ0), ϕ̄FF=ϕFF/(qT/λ0) and v̄T=vT/CD. We also choose some recent simulation data (Bandyopadhyay et al., 2012) that has been used in the studies of Mach cones in dusty plasmas, namely, the equilibrium ion density ni0 ∼ 108cm−3, the equilibrium electron density ne0 ∼ 1.9 × 107 cm−3, the dust charging state zd0 = 4,007, the dust radius rd ≃ 4.5 μm, the average inter-particle distance d ≃ 230 μm, the ion temperature Ti = 0.1eV, and the electron temperature Te = 1.5eV. These values further help to find other physical quantities, e.g., the dust plasma frequency ωpd ∼ 68Hz, the dust thermal speed vTd ∼ 0.025 cm/s, the DA speed CD = ωpdλ0 ≡ 1.36 cm/s with λ0 ∼ 0.02 cm and dust polarization force coefficient R=0.05. The latter strongly depends on the dust charging state and leads to the reduction of the phase speed of DA waves in a polarized dusty plasma.

The modifications due to dust-charge fluctuations and dust-neutral collisions (Ali et al., 2003) have been ignored in the present model, just to make more simple the analytical calculations and analysis. The dust charging process is actually a complex phenomenon (Popel et al., 2005) that occurs on a very fast scale in laboratory plasmas (Hazelton and Yadlowsky, 1994; Praburam and Goree, 1996) as well as in numerical simulations (Choi and Kushner, 1994). Hence, ignoring the dust charge variation effects, the primary focus is to examine the impact of dust polarization force on the shielding and dynamical potentials by imposing certain conditions on the test charge speed in relation to the thermal speeds of plasma species. However, an effort has already been made to investigate the characteristics of the Debye shielding with dust-charge fluctuation and dust-neutral collisional effects in negatively charged dusty plasmas (Ali et al., 2003). Later, Ali (Ali, 2016) continued his investigation to study a positive-dust electron plasma and examined the influence of secondary electron and photoelectron emissions on the potential profiles.

Figure 1 displays the normalized WF potential (ϕ̄WF) of a test charge against the normalized axial distance (ξ̃) as a function of (A) dust charging state 1,000 ≤ zd0 ≤ 6,000 with η ≃ 10–6 and (B) dust concentration 10–4η ≤ 8 × 10–6 with zd0 = 4,007 at constant speed vT=0.3CD. When the test charge resonates with a DA mode, an asymmetric wakefield is excited behind the test charge in the form of potential distribution in an axial direction. Since the dust polarization coefficient R (zd0, η) is the function of dust charging state and dust concentration, therefore, the variation in these parameters significantly mitigates (increases) the magnitude of the WF potentials, respectively. The wakefield is oscillatory in nature and appears in terms of positive and negative potential regions that damp behind the test charge particle for several Debye lengths [see Figure 1A]. The strength of the WF potential decreases because the dust polarization force coefficient effectively reduces the phase speed of the collective DA mode. It is important to note that Eq. 17 is only valid in the limit when R <1 for micron-sized dust grains. However, the variation of dust concentration which appears through an equilibrium charge-neutrality condition, significantly modifies the amplitudes of wakefield [as shown by Figure 1B] and its potential regions. The increase in dust concentration essentially makes the depletion of electrons in the system and consequently, the phase speed of the DA wave increases. This effect pronounces the potential associated to the wakefield, and also has fundamental importance in attracting the same polarity dust grains and forming the ordered crystalline structures in the (polarized) dusty plasma (Fortov et al., 1996). If a test charge moves with a speed slower than the dust thermal speed (viz., vTvTd), then all the plasma species (viz., the electrons, ions, and negatively charged dust grains) would shield the test particle, which may lead to significant modification of the shielding length.

FIGURE 1
www.frontiersin.org

FIGURE 1. The profiles of WF potential (ϕ̄WF) are displayed against the normalized axial distance (ξ̄) as a function of (A) dust charging state 1000 ≤ zd0 ≤ 6000 with η ≃ 10–6 and (B) dust concentration 10–4η ≤ 8 × 10–6 with zd0 = 4007 at constant speed vT=0.3CD. All other parameters are the same as mentioned in Section 3.

In Figure 2, we show how effective shielding length varies as function of electron-to-dust temperature ratio for different values of dust polarization force coefficient R=0.057 (black dashed curve), 0.058 (blue dashed curve), 0.059 (red solid curve). It is observed that dust polarization coefficient is effectively modified by the dust concentration parameter η(= 10–5, 3 × 10–5, 9 × 10–5) at fixed charging state zd0 = 4,007 and so the shielding effect diminishes in dusty plasmas. Moreover, the variation of dust charging states zd0=2000,3000,4990 not only affects the dust polarization coefficient R but also the effective shielding length, showing a reduction in the shielding length as long as the electron-to-dust temperature ratio varies as a function.

FIGURE 2
www.frontiersin.org

FIGURE 2. The profiles of the effective-shielding length (λ̄Def) are shown as a function of electron-to-dust temperature ratio (Te/Td) for varying (A) the polarization coefficient R=0.057 (black dashed curve), 0.058 (blue dashed curve), 0.059 (red solid curve) with zd0 = 4007, and (B) the dust charging state zd0 = 2000, 3000, 4990 at η = 10–5. All other parameters are the same as given in Section 3.

Figure 3 depicts the effect of the dust polarization force coefficient R on the profiles of normalized DH potential (ϕ̄DH) against the normalized axial distance ξ̄ at fixed radial distance ρ = 0.15 for (A) zd0 = 3500 and (B) zd0 = 5000. It may be noticed that the values of dust polarization coefficient R[ = 0.05053 (black dashed curve), 0.05054 (blue dashed curve), 0.05057 (red solid curve)] are not only dependent on the dust concentration parameter η(= 10–6, 2 × 10–6, 4 × 10–6) but also on the dust charging state and leading to the reduction of DH potential, as can be seen from Figure 3A. However, keeping the same dust concentrations at high dust charging state zd0 (= 5000) would result in the modification of dust polarization coefficient R (= 0.07148, 0.07151, 0.07156) with more decreased DH potential profiles [see Figure 3B]. The DH potentials are mainly caused by the static or slow moving test charges in a polarized dusty plasma. Figure 4 represents how a normalized FF potential (ϕ̄FF) is affected by the normalized axial distance ξ̄ for varying (A) the dust polarization R [ = 0.05766 (black dashed curve), 0.05769 (blue dashed curve), 0.05772 (red solid curve)] at zd0 = 4007 and (B) the dust polarization R [ = 0.06454 (black dashed curve), 0.06459 (blue dashed curve), 0.06463(red solid curve)] at zd0 = 4500. Such type of potentials often arise from very slow motion of test charges as compared to dust thermal speed, viz., vT = 0.02vTd. The magnitudes of the FF potentials become lower (higher) at high (low) dust charging state in both the cases and are significantly altered by the impact of dust polarization coefficient. The FF potentials are only derived in the long-range limit when rλDef.

FIGURE 3
www.frontiersin.org

FIGURE 3. The variation of DH potential (ϕ̄DH) is displayed against the normalized axial distance (ξ̄) for changing dust polarization coefficient (R) at low and high charging states (A) zd0 = 3500 and (B) zd0 = 5000 with ρ̄=0.15. All other parameters are the same as given in Section 3.

FIGURE 4
www.frontiersin.org

FIGURE 4. The profiles of the FF potential (ϕ̄FF) are plotted against the normalized axial distance (ξ̄) for varying the polarization coefficient (R) at low and high charging states (A) zd0 = 4007 and (B) zd0 = 4500 with ρ̄=0.1, v̄T=0.02 and Te/Td = 20. All other parameters are the same as given in Section 3.

4 Summary and conclusion

To summarize, we have studied the test charge-driven response of a polarized dusty plasma, showing distinct features of interaction potentials. In this context, we have solved the coupled linearized set of Vlasov-Poisson equations by using the Fourier transformation technique and derived an expression for the total electrostatic potential with a modified dielectric constant that accounts for dust polarization force. The latter acts as an opposite force to the electrostatic force and plays a vital role in interaction potentials for micron-sized dust particulates. Imposing certain conditions on the test charge speed in comparison with thermal speeds of plasma species, the total electrostatic potential is analyzed for various limiting cases, obtaining oscillatory wakefield (WF), Debye-Hückel (DH), Farfield (FF) and Coulomb potentials. It is found that dust polarization force coefficient is essentially dependent on the dust charging state and dust concentration to modifying the profiles of DH potential, which decay exponentially with the distance, the FF potential which decays as the inverse cube of the distance and WF potential behind a test charge due to resonant interaction. Moreover, dust polarization and dust concentration do not influence the Coulomb potential because the test charge moves very fast relative to plasma species and so no shielding occurs around it in the polarized dusty plasma. Many years ago (Nambu et al., 1995; Vladimirov and Nambu, 1995; Salimullah and Nambu, 2000; Nambu et al., 2001), the concept of wakefield was proposed for making new materials with help of attractive forces between the same polarity dust grains in negative regions of the potential. The negatively charged dust grains are attracted in the similar way as that of Cooper pairing of electrons in superconductors (De Gennes-Pierre, 1966) and in turn, may lead to the possibility to have dust crystallization and dust-coagulation in dusty plasmas.

Test charge technique is one of the effective techniques, which is used as a diagnostic tool for studying interaction potentials in dusty plasmas. The present findings are important to understand the physics of the shielding phenomenon and to explore new features of dynamical potentials caused by the wave-particle interaction in the presence of dust polarization force. Various effects like dust-charge perturbations and dust-neutral collisions (Ali et al., 2003; Ali et al., 2005a; Ali et al., 2005b), external magnetic field and ion-streaming (Vladimirov and Nambu, 1995; Nambu et al., 2001), two-body correlations, etc., have been investigated with significant modifications to potential distribution both analytically and numerically. In many laboratories, the massive charged dust grains are assumed to be strongly coupled due to their high dust charge, low temperature, and small intergrain spacing. Relying on the specific simulation data (Bandyopadhyay et al., 2012), the test charge potential profiles are analyzed in a negatively charged dusty plasma, which corresponds to a strongly coupled regime owing to a large coupling parameter, i.e., Γc ≫ 1. However, dust temperature can play a decisive role in the transition from strongly coupled to weakly coupled dusty plasmas and vice versa. In particular, dust-charge variation and dust-neutral collisions have significantly modified the profiles of the shielding and wakefield potentials elsewhere (Ali et al., 2003), but in the present model, it is beyond the scope of our interest and would be addressed in future studies.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

Author contributions

SA: Introduced the main idea and performed both analytical and numerical calculations. WM: Involved in write-up and analysis. KS: Discussed modeling and findings. RJ: Analyzed calculations and prepared analysis.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

Ali, S., and Eliasson, B. (2017). Slow test charge response in a dusty plasma with Kappa distributed electrons and ions. Phys. Scr. 92, 084003. doi:10.1088/1402-4896/aa7c09

CrossRef Full Text | Google Scholar

Ali, S., Murtaza, G., and Nasim-Maqsood, H. (2005). Energy loss for the assemblies of charged projectiles in a dusty plasma. Phys. Plasmas 12, 072104. doi:10.1063/1.1933739

CrossRef Full Text | Google Scholar

Ali, S., Nasim, M. H., and Murtaza, G. (2003). Effects of dust-charge fluctuations on the potential of an array of projectiles in a partially ionized dusty plasma. Phys. Plasmas 10, 4207–4216. doi:10.1063/1.1619976

CrossRef Full Text | Google Scholar

Ali, S., Nasim-Maqsood, H., and Murtaza, G. (2005). Correlation effects due to an axial propagation of projectiles in a dusty plasma. Phys. Plasmas 12, 033502. doi:10.1063/1.1844496

CrossRef Full Text | Google Scholar

Ali, S. (2016). Potential distribution around a test charge in a positive dust-electron plasma. Front. Phys. (Beijing). 11 (3), 115201–115209. doi:10.1007/s11467-015-0545-2

CrossRef Full Text | Google Scholar

Ali, S. (2009). Test charge potential in the presence of electron acoustic waves in multispecies dusty plasma. Phys. Plasmas 16 (5), 054502. doi:10.1063/1.3122941

CrossRef Full Text | Google Scholar

Al’Pert, Y. L., Gurevich, A. V., and Pitaevskiĭ, L. P. (1963). Effects produced by an artificial satellite rapidly moving in the ionosphere or in an interplanetary medium. Sov. Phys. Usp. 6 (1), 13–46. doi:10.1070/PU1963v006n01ABEH003492

CrossRef Full Text | Google Scholar

Bandyopadhyay, P., Jiang, K., Dey, R., and Morfill, G. E. (2012). Effect of polarization force on the Mach cones in a complex plasma. Phys. Plasmas 19 (12), 123707. doi:10.1063/1.4773194

CrossRef Full Text | Google Scholar

Barkan, A., D’angelo, N., and Merlino, R. L. (1996). Experiments on ion-acoustic waves in dusty plasmas. Planet. Space Sci. 44 (3), 239–242. doi:10.1016/0032-0633(95)00109-3

CrossRef Full Text | Google Scholar

Barkan, A., Merlino, R. L., and D’Angelo, N. (1995). Laboratory observation of the dust-acoustic wave mode. Phys. Plasmas 2, 3563–3565. doi:10.1063/1.871121

CrossRef Full Text | Google Scholar

Bentabet, K., Mayout, S., and Tribeche, M. (2017). Generalized polarization force acting on dust grains in a dusty plasma. Phys. A Stat. Mech. its Appl. 466, 492–501. doi:10.1016/j.physa.2016.09.055

CrossRef Full Text | Google Scholar

Bhatnagar, P. L., Gross, E. P., and Krook, M. (1954). A model for collision processes in gases. I. Small amplitude processes in charged and neutral one-component systems. Phys. Rev. 94 (3), 511–525. doi:10.1103/PhysRev.94.511

CrossRef Full Text | Google Scholar

Chen, F. F. Introduction to plasma physics 1974. New York: Plenum, 249.

Chen, Z., Yu, M. Y., and Luo, H. (2001). Molecular dynamics simulation of two-dimensional dust clusters. Phys. Scr. 64 (5), 476–481. doi:10.1238/Physica.Regular.064a00476

CrossRef Full Text | Google Scholar

Choi, S. J., and Kushner, M. J. (1994). A particle-in-cell simulation of dust charging and shielding in low pressure glow discharges. IEEE Trans. Plasma Sci. IEEE Nucl. Plasma Sci. Soc. 22 (2), 138–150. doi:10.1109/27.279017

CrossRef Full Text | Google Scholar

Clemmow, P. C., and Dougherty, J. P. Electrodynamics of particles and plasmas 1969. London: Addison-Wesley.

Cooper, G. (1969). Shielding of slow test particles in a plasma. Phys. Fluids (1994). 12 (12), 2707–2710. doi:10.1063/1.1692416

CrossRef Full Text | Google Scholar

De Gennes-Pierre, G. (1966). Superconductivity of metals and alloys. New York: Benjamin.

Google Scholar

Debye, P., Hückel, E., and De La Theorie Des Electrolytes, (1923). I. Abaissement du Point de Congelation et Phenomenes Associes. Phys. Z 24, 185.

Google Scholar

Denysenko, I., Yu, M. Y., Stenflo, L., and Azarenkov, N. A. (2005). Ion drag on dust grains in electronegative plasmas. Phys. Plasmas 12 (4), 042102. doi:10.1063/1.1862629

CrossRef Full Text | Google Scholar

Fortov, V. E., Nefedov, A. P., Petrov, O. F., Samarian, A. A., and Chernyschev, A. V. (1996). Particle ordered structures in a strongly coupled classical thermal plasma. Phys. Rev. E 54, R2236–R2239. doi:10.1103/physreve.54.r2236

PubMed Abstract | CrossRef Full Text | Google Scholar

Fried, D. B., and Conte, S. D. (1961). The plasma dispersion function. New York: Academic Press.

Google Scholar

Goertz, C. K. (1989). Dusty plasmas in the solar system. Rev. Geophys. 27 (2), 271–292. doi:10.1029/RG027i002p00271

CrossRef Full Text | Google Scholar

Gurevich, A. V., Pitaevskii, L. P., and Smirnova, V. V. (1969). Ionospheric aerodynamics. Space Sci. Rev. 9, 805–871. doi:10.1007/BF00226263

CrossRef Full Text | Google Scholar

Hamaguchi, S., and Farouki, R. T. (1994). Plasma–particulate interactions in nonuniform plasmas with finite flows. Phys. Plasmas 1 (7), 2110–2118. doi:10.1063/1.870608

CrossRef Full Text | Google Scholar

Hamaguchi, S., and Farouki, R. T. (1994). Polarization force on a charged particulate in a nonuniform plasma. Phys. Rev. E 49 (5), 4430–4441. doi:10.1103/PhysRevE.49.4430

PubMed Abstract | CrossRef Full Text | Google Scholar

Hazelton, R. C., and Yadlowsky, E. J. (1994). Measurement of dust grain charging in a laboratory plasma. IEEE Trans. Plasma Sci. IEEE Nucl. Plasma Sci. Soc. 22 (2), 91–96. doi:10.1109/27.279009

CrossRef Full Text | Google Scholar

Horányi, M., Hartquist, T. W., Havnes, O., Mendis, D. A., and Morfill, G. E. (2004). Dusty plasma effects in Saturn’s magnetosphere. Rev. Geophys. 42 (4), RG4002. doi:10.1029/2004RG000151

CrossRef Full Text | Google Scholar

Joyce, G., and Montgomery, D. (1967). Shielding in anisotropic plasmas. Phys. Fluids (1994). 10 (9), 2017–2020. doi:10.1063/1.1762401

CrossRef Full Text | Google Scholar

Kan, J. R. (1971). A note on the formation of Debye potential well in current-carrying plasmas. J. Plasma Phys. 6 (1), 223–228. doi:10.1017/S0022377800025800

CrossRef Full Text | Google Scholar

Khrapak, S. A., Ivlev, A. V., Yaroshenko, V. V., and Morfill, G. E. (2009). Influence of a polarization force on dust acoustic waves. Phys. Rev. Lett. 102 (24), 245004. doi:10.1103/PhysRevLett.102.245004

PubMed Abstract | CrossRef Full Text | Google Scholar

Krall, N. A., and Trivelpiece, A. W. (1973). Principles of plasma physics. New York: McGraw-Hill.

Google Scholar

Laingal, E. W., Lamont, A., and Fielding, P. J. (1971). Interaction of test particles with a plasma. J. Plasma Phys. 5 (3), 441–454. doi:10.1017/S0022377800005924

CrossRef Full Text | Google Scholar

Lakshmi, S. V., Bharuthram, R., and Shukla, P. K. (1993). Debye shielding in a dusty plasma. Astrophys. Space Sci. 209, 213–219. doi:10.1007/BF00627441

CrossRef Full Text | Google Scholar

Lemons, D. S., Murillo, M. S., Daughton, W., and Winske, D. (2000). Two-dimensional wake potentials in sub-and supersonic dusty plasmas. Phys. Plasmas 7 (6), 2306–2313. doi:10.1063/1.874065

CrossRef Full Text | Google Scholar

Liu, V. C. (1965). A wave mechanical approach to the plasma interaction problems. Nature 208 (5013), 883–884. doi:10.1038/208883a0

CrossRef Full Text | Google Scholar

Mamun, A. A., Ashrafi, K. S., and Shukla, P. K. (2010). Effects of polarization force and effective dust temperature on dust-acoustic solitary and shock waves in a strongly coupled dusty plasma. Phys. Rev. E 82 (2), 026405. doi:10.1103/PhysRevE.82.026405

PubMed Abstract | CrossRef Full Text | Google Scholar

Masood, W., and Mirza, A. M. (2010). Electron thermal effect on linear and nonlinear coupled Shukla–Varma and convective cell modes in dust-contaminated magnetoplasma. Phys. Plasmas 17 (11), 113702. doi:10.1063/1.3505094

CrossRef Full Text | Google Scholar

Masood, W., Mirza, A. M., and Nargis, S. (2010). Revisiting coupled Shukla–Varma and convective cell mode in classical and quantum dusty magnetoplasmas. J. Plasma Phys. 76 (3-4), 547–552. doi:10.1017/S0022377809990791

CrossRef Full Text | Google Scholar

Masood, W., Rizvi, H., Hasnain, H., and Haque, Q. (2012). Rotation induced nonlinear dispersive dust drift waves can be the progenitors of spokes. Phys. Plasmas 19 (3), 032112. doi:10.1063/1.3688077

CrossRef Full Text | Google Scholar

Masood, W., Shah, H. A., Tsintsadze, N. L., and Qureshi, M. N. (2010). Dust Alfven ordinary and cusp solitons and modulational instability in a self-gravitating magneto-radiative plasma. Eur. Phys. J. D. 59 (3), 413–419. doi:10.1140/epjd/e2010-00180-x

CrossRef Full Text | Google Scholar

Mehdipoor, M. (2022). Nonlinear propagation of dust-acoustic modes in a dusty plasma with the Kappa-Cairns polarization force effects. Waves in Random and Complex Media. doi:10.1080/17455030.2022.2066217

CrossRef Full Text | Google Scholar

Montgomery, D., Joyce, G., and Sugihara, R. (1968). Inverse third power law for the shielding of test particles. Plasma Phys. 10 (7), 681–686. doi:10.1088/0032-1028/10/7/304

CrossRef Full Text | Google Scholar

Nambu, M., and Akama, H. (1985). Attractive potential between resonant electrons. Phys. Fluids (1994). 28 (7), 2300–2301. doi:10.1063/1.865284

CrossRef Full Text | Google Scholar

Nambu, M., Saikia, B. J., and Hada, T. (2001). Wake potential around a test dust particulate in a magnetized plasma with streaming ions. J. Phys. Soc. Jpn. 70, 1175–1177. doi:10.1143/JPSJ.70.1175

CrossRef Full Text | Google Scholar

Nambu, M., Vladimirov, S. V., and Shukla, P. K. (1995). Attractive forces between charged particulates in plasmas. Phys. Lett. A 203, 40–42. doi:10.1016/0375-9601(95)00380-l

CrossRef Full Text | Google Scholar

Nasim, M. H. (1999). Energy loss of charged projectiles in a dusty plasma. Islamabad, Pakistan: Ph.D thesis, Quaid-i-Azam University.

Google Scholar

Nawaz, H., Jahangir, R., Masood, W., and Siddiq, M. (2022). Cubic nonlinearity driven dust ion acoustic solitons with superthermal two-temperature electrons in the Saturn’s magnetosphere. Phys. Scr. 97 (7), 075601. doi:10.1088/1402-4896/ac735e

CrossRef Full Text | Google Scholar

Neufeld, J., and Ritchie, R. H. (1955). Passage of charged particles through plasma. Phys. Rev. 98 (6), 1632–1642. doi:10.1103/PhysRev.98.1632

CrossRef Full Text | Google Scholar

Popel, S. I., Losseva, T. V., Golub’, A. P., Merlino, R. L., and Andreev, S. N. (2005). Dust ion-acoustic shocks in a Q-Machine device. Contrib. Plasma Phys. 45 (7), 461–475. doi:10.1002/ctpp.200510052

CrossRef Full Text | Google Scholar

Praburam, G., and Goree, J. (1996). Experimental observation of very low-frequency macroscopic modes in a dusty plasma. Phys. Plasmas 3 (4), 1212–1219. doi:10.1063/1.871745

CrossRef Full Text | Google Scholar

Rao, N. N., Shukla, P. K., and Yu, M. Y. (1990). Dust-acoustic waves in dusty plasmas. Planet. Space Sci. 38 (4), 543–546. doi:10.1016/0032-0633(90)90147-I

CrossRef Full Text | Google Scholar

Sabeen, A., Masood, W., Qureshi, M. N., and Shah, H. A. (2017). Nonlinear coupling of kinetic Alfven waves with acoustic waves in a self-gravitating dusty plasma with adiabatic trapping. Phys. Plasmas 24 (7), 073704. doi:10.1063/1.4990700

CrossRef Full Text | Google Scholar

Salimullah, M., and Nambu, M. (2000). Crystallization in a magnetized and inhomogeneous dusty plasma with streaming ions. J. Phys. Soc. Jpn. 69, 1688–1691. doi:10.1143/JPSJ.69.1688

CrossRef Full Text | Google Scholar

Shahmansouri, M., Farokhi, M. B., and Khodabakhshi, N. (2017). Modified potential around a moving test charge in strongly coupled dusty plasma. Commun. Theor. Phys. 68 (1), 111. doi:10.1088/0253-6102/68/1/111

CrossRef Full Text | Google Scholar

Shahmansouri, M., and Khodabakhshi, N. (2018). On the dielectric response function and dispersion relation in strongly coupled magnetized dusty plasmas. Chin. Phys. B 27 (10), 105206. doi:10.1088/1674-1056/27/10/105206

CrossRef Full Text | Google Scholar

Shahmansouri, M., and Mamun, A. A. (2016). Generalized polarization force acting on charge fluctuating dust grains and its effects on propagation of dust-acoustic waves in a dusty plasma. Eur. Phys. J. Plus 131 (9), 321. doi:10.1140/epjp/i2016-16321-2

CrossRef Full Text | Google Scholar

Shahmansouri, M., and Misra, A. P. (2019). Polarized Debye sheath in degenerate plasmas. Commun. Theor. Phys. 71 (11), 1341–1345. doi:10.1088/0253-6102/71/11/1341

CrossRef Full Text | Google Scholar

Shahmansouri, M. (2019). Potential distribution in a strongly coupled dusty magnetoplasma. IEEE Trans. Plasma Sci. 47 (11), 5108–5112. doi:10.1109/TPS.2019.2946182

CrossRef Full Text | Google Scholar

Shohaib, M., Masood, W., Siddiq, M., Alyousef, H. A., and El-Tantawy, S. A. (2022). Formation of electrostatic solitary and periodic waves in dusty plasmas in the light of Voyager 1 and 2 spacecraft and Freja satellite observations. J. Low Freq. Noise, Vib. Act. Control 12, 896–909. doi:10.1177/14613484221091340

CrossRef Full Text | Google Scholar

Shukla, P. K. (2001). A survey of dusty plasma physics. Phys. Plasmas 8 (5), 1791–1803. doi:10.1063/1.1343087

CrossRef Full Text | Google Scholar

Shukla, P. K., and Mamun, A. A. (2002). Introduction to dusty plasma physics. Bristol: Institute of Physics. doi:10.1201/9781420034103

CrossRef Full Text | Google Scholar

Shukla, P. K., and Mamun, A. A. (2003). Solitons, shocks and vortices in dusty plasmas. New J. Phys. 5 (1), 17. doi:10.1088/1367-2630/5/1/317

CrossRef Full Text | Google Scholar

Shukla, P. K., and Silin, V. P. (1992). Dust ion-acoustic wave. Phys. Scr. 45 (5), 508. doi:10.1088/0031-8949/45/5/015

CrossRef Full Text | Google Scholar

Singh, K., Arab, P., and Saini, N. S. (2018). Effect of polarization force on head-on collision between multi-solitons in the the dusty plasma. Phys. Plasmas 25 (3), 033705. doi:10.1063/1.5020194

CrossRef Full Text | Google Scholar

Singh, K., Ghai, Y., Kaur, N., and Saini, N. S. (2018). Effect of polarization force on dust-acoustic cnoidal waves in dusty plasma. Eur. Phys. J. D. 72 (9), 160–168. doi:10.1140/epjd/e2018-90228-2

CrossRef Full Text | Google Scholar

Stenflo, L., Yu, M. Y., and Shukla, P. K. (1973). Shielding of a slow test charge in a collisional plasma. Phys. Fluids (1994). 16 (3), 450–452. doi:10.1063/1.1694361

CrossRef Full Text | Google Scholar

Stoffels, E., Stoffels, W. W., Kersten, H., Swinkels, G. H., and Kroesen, G. M. (2001). Surface processes of dust particles in low pressure plasmas. Phys. Scr. 2001 (T89), 168. doi:10.1238/Physica.Topical.089a00168

CrossRef Full Text | Google Scholar

Vladimirov, S. V., and Nambu, M. (1995). Attraction of charged particulates in plasmas with finite flows. Phys. Rev. E 52, R2172–R2174. doi:10.1103/physreve.52.r2172

PubMed Abstract | CrossRef Full Text | Google Scholar

Yaqoob, F., Ali, S., Murtaza, G., and Nasim, M. H. (2002). Debye shielding distortion of dust grains in dusty plasmas. Phys. Plasmas 9, 3629–3632. doi:10.1063/1.1491956

CrossRef Full Text | Google Scholar

Yu, M. Y., Tegeback, R., and Stenflo, L. (1973). On test charge potentials in collisional plasmas. Z. Phys. 264, 341–348. doi:10.1007/BF01398860

CrossRef Full Text | Google Scholar

Keywords: dusty plasmas, dust polarization force, kinetic description of plasma particles, test charge technique, interaction potentials

Citation: Ali S, Masood W, Singh K and Jahangir R (2022) Test charge driven response of a dusty plasma with polarization force. Front. Astron. Space Sci. 9:987561. doi: 10.3389/fspas.2022.987561

Received: 06 July 2022; Accepted: 08 September 2022;
Published: 11 October 2022.

Edited by:

Amar Prasad Misra, Visva-Bharati University, India

Reviewed by:

Sergey Popel, Space Research Institute (RAS), Russia
Mehran Shahmansouri, Arak University, Iran

Copyright © 2022 Ali, Masood, Singh and Jahangir. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: W. Masood, waqasmas@gmail.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.