Slow test charge response in a dusty plasma with Kappa distributed electrons and ions

The electrostatic potential around a slowly moving test charge is studied in a dusty plasma where the ions and electrons follow a powerlaw Kappa distribution in velocity space. A test charge moving with a speed much smaller than the dust thermal speed gives rise to a short-scale Debye–Hückel potential as well as a long-range far-field potential decreasing as inverse cube of the distance to the test charge along the propagation direction. The potentials are significantly modified in the presence of high-energy tails, modeled by lower spectral indices in the ion and electron Kappa distribution functions. Plasma parameters relevant to laboratory dusty plasmas are discussed.


Introduction
Fluid and kinetic models are two common approaches to study waves and instabilities in plasma [1,2]. Test charge techniques [3][4][5][6] are often utilized for studying shielding phenomena such as electrostatic potential and energy loss caused by a moving test charge in plasma. In this technique, a test charge with a given charge and velocity is externally inserted into the plasma and its charge is coupled with the plasma charge density through the space charge effects, leading to the shielding of the charge by a cloud of opposite charges in the plasma. Problems involving shielded electrostatic potentials around charged particles and energy loss of particles moving in plasma occur in many different areas of physics, such as particle acceleration [7,8], the formation of new materials in low-temperature dusty plasma [9][10][11], heating and ignition in inertial confinement fusion [12][13][14][15][16], heavy ion energy deposition in dense plasma [17], and the interaction of ion clusters with condensed matter [18][19][20][21].
Taking into account moving and stationary test charges, a far-field (FF) potential distribution depending as the inverse third power of the distance to the test charge was found around a test charge in a Maxwellian plasma of dynamical electrons with static ions [4]. Later, the effects of ionic motions [22,23], collisions [24][25][26] and turbulence [27,28] were taken into account. A kinetic description was used to study the electrostatic potential around a test charge with low and high velocities [29,30], revealing the excitation of longrange wakefields in a Maxwellian plasma. The effects of magnetic field, collisions and plasma density gradients have also been studied for slowly and rapidly moving test charges [31][32][33]. The Debye shielding has been studied in a dusty plasma with Boltzmann distributed electrons and ions in the presence of static negatively charged dust grains [34], taking into account both small and finite amplitude potential distributions. Collisions, Landau damping, and dust-charge perturbation effects on the dipole-like FF potential have been taken into account for a slow test charge in a Maxwellian dusty plasma [35,36]. Oscillatory wake-field, Debye-Hückel (DH) and FF potentials caused by a test charge in a colloidal Maxwellian plasma in the presence of streaming ions and dust grains have also been studied theoretically [37] and numerically [38]. Recently, the shielding of two slow comoving teset charges in collisional plasma has also been studied theoretically [39].
There is a growing interest in studying non-Maxwellian power-law distributions involving superthermal particles, accelerated by wave-particle interactions and turbulence. It has opened up new research avenues for researchers to model plasma collective interactions and instabilities in astrophysical environments [40]. Suprathermal electrons often occur in plasmas where modulational instabilities and Langmuir turbulence take place [41], such as beam-plasma interactions [42], in the solar wind where type III solar radio emissions occurs [43], intense microwave-plasma interactions [44] and ionospheric heating experiments [45,46]. Superthermal energetic tails in the velocity distribution function of particle species s (s=e for electrons and i for ions) can be modeled by a three-dimensional Kappa distribution function [47] k q = + Ts s s s s 1 2 is the equilibrium density and effective thermal speed, respectively. The Kappa distribution is valid for k > 3 2 s . The temperature is denoted by T s and the mass by m s . The Kappa distribution function tends to a Maxwellian distribution function, Vasyliunas [47] first pointed out the relevance of the Kappa distribution function by fitting empirical data from the solar wind to show its importance for low spectral index κ∼2-4. After that, the investigation of several dispersive modes and instabilities have been carried out for power-law Kappa-distributed plasmas [48][49][50][51]. Recently, the potential around a slowly moving test charge in a Kappadistributed plasma was studied [52] using a Vlasov model, containing supra-thermal hot and cold electrons with immobile ions. It was shown that the short-range DH potential decreases as an exponential decay of the distance and longrange FF potential as the inverse cube of the distance from the test charge.
In this paper, we present a theory for the potential distributions around a slowly moving test charge in an unmagnetized dusty plasma having supra-thermal electrons and ions, as well as negatively charged dynamical dust grains.

Vlasov-Poisson model
The electrostatic potential due to a moving test charge in a collisionless dusty plasma is in the small-amplitude limit governed by the linearized Vlasov-Poisson system where Q T is the charge of the test particle and δ is Dirac's delta function representing the test charge located at r. We have used the equilibrium charge-neutrality condition 0 0 0 , where N j0 is the particle number density of the jth species ( j=e for electrons, i for positive ions, and d for the negatively charged dust grains), and Z d0 is the equilibrium dust charge state. We assume that all dust grains are spherical in shape with radius r d and mass m d . The linearized electron and ion density fluctuations f for Kappa-distributed inertialess electrons and ions with temperatures T e and T i . Furthermore, the perturbation of the dust distribution function is assumed to be much smaller compared to the equilibrium Maxwellian distribution, is the dust thermal speed and T d the dust temperature. Note that the effects of suprathermal tails are accounted for through the parameters The shielded test charge Q T gives rise to the screened potential f 1 , which depends on its speed in comparison with the thermal speeds.

Dielectric response to a test charge
Applying the space-time Fourier transform to equations (1) and (2), we obtain Here ω and k are the angular frequency and wave vector, f 1 is the Fourier transformed potential, the plasma dispersion function [53], and Taking the inverse Fourier transform of equation (3) gives w = · k v T , and the potential is obtained as [2,54] The dielectric constant can in the limit In equation (6), we have made an expansion of ( ) W C d function in the limit T d , to account for a slowly moving test charge, ignoring the higher order terms. Note that the plasma response function, and in turn the test charge distribution, are modified by the dust Landau damping term [4] proportional to C d in right-hand side of equation (6). For = v 0 T , the dust Landau damping vanishes, and the shielded test charge gives rise to the standard DH potential [52,[55][56][57].
Taking the inverse of equation (6) in the limit  [34] in the limit of static dust grains. Figure 1 displays the impact of the dust temperature and charge state on the normalized effective Debye length given by equation (8). Since all plasma species participate in the shielding process when the test charge is moving slowly, the effective shielding length is also influenced by the dust temperature and charge state. We note that the effective shielding length is reduced in a dusty plasma compared to an electron-ion plasma due to the cold dust population shielding the test charge. . Consistent with previous results [60], the effective Debye length is smaller for smaller values of the spectral indices.

Profiles of the DH and FF potentials
Choosing a coordinate system in which the test charge is moving in the z-direction, Equation (9) describes the short-range DH potential, accounting for small distances between the test charge and the observer, and equation (10) the long-range FF potential derived in the limit l k  r 1 Def , which to leading order decays as the inverse cube of the distance to the test charge.
Here ), who also considered the effects of collisions and dust charge fluctuations on the FF potential.
As a numerical illustration, we have normalized equations  give a relatively small impact on the DH potentials (see figure 3(b)). This is expected, since the shielding is primarily due to the colder species (ions and dust), while the hot electrons do not significantly participate in the DH shielding. The shielding FF potential is more localized for smaller ion spectral indices k = 1.6 i and 2 compared to the near-Maxwellian case k = 100 i (see figure 4(a)), and the shielding is also more efficient for lower values of the dust temperatures ( figure 4(d)). The values of the electron spectral index and the dust concentration have only minor effects on the FF shielding potential in figures 4(b) and (c).

Summary
To summarize, we have presented slow response of a moving test charge in a non-Maxwellian dusty plasma, whose constituents are supra-thermal and inertialess electrons and ions in the presence of Maxwell-distributed inertial dust grains. For this purpose, linearized Vlasov equation for the dust,   coupled with the electron and ion density responses and Poisson's equation for the potential are solved with a Fourier transform technique to obtain the modified electrostatic potential modified by the supra-thermal tails of the ions and electrons. When the test charge moves with a constant speed, slowly in comparison with the dust thermal speed, it gives rise to a shielded potential that can be decomposed into a shortrange DH and long-range FF potential. The variation of FF potential is to leading order proportional to the inverse cube of the distance to the test charge. Both the DH and FF shielding potentials are more localized (i.e. more efficient shielding) for smaller spectral indices. For large spectral indices, the DH and FF potentials approach the Maxwellian case.