Collective Mode Interactions in Lorentzian Space Plasma

Plasmas exhibit a vast variety of waves and oscillations in which moving charged particle produce fields which ultimately give rise to particle motion. These wave-particle effects are used in the acceleration heating methods of plasma particles, and in wave generation as well. Plasmas are often manipulated with EM waves, e.g., Alfvén waves are long-wavelength modes (drift-waves) where fluid theory is most reliable, while for short wavelength modes (e.g., Kinetic Alfvén waves), collisionless effects becomes important. In this chapter, the properties of kinetic Alfvén waves are aimed to study by employing two potential theory by taking particle streaming and Weibel instability with temperature anisotropy in a Lorentzian plasma.


Introduction
This chapter addresses one of the intriguing topics of Astrophysics-the existence of kinetic Alfvén wave (KAW) and the important consequences for astrophysical and space science to explore and investigate the new avenues. Due to the fact that KAWs have non-zero electric field E ∥ which is parallel to background magnetic field and possess anisotropic polarized and spatial structures which contribute to particle energization. It is an interesting mechanism that KAWs can accelerate the field-aligned charged particles and has been applied in the dissipation of solar wind turbulence, the acceleration, and heating of charged particles in both the filed-aligned and perpendicular directions and is anticipated to play a vital role in the particle energization in laboratory, space and astrophysical plasma. The progress reported here would have immense impact and hence a small step in particular direction.
The solar wind plasma is hot and weakly collisional, existing in a state far from thermal equilibrium [1] as observed in situ in the solar wind through its nonthermal characteristics of velocity distribution function (VDF). The electron VDFs measured at 1 AU have been used as boundary condition to determine the VDFs at different altitudes. It has been confirmed that for several solar radii, the suprathermal population of particles is present in the corona [2]. For low collision rates in such plasmas, the particles can develop temperature anisotropy and the VDFs become slanted and build up high energy tails and heat fluxes along the magnetic field direction especially in fast winds and energetic interplanetary shocks. Various processes in a collisionless solar wind plasma lead to the development of particle temperature anisotropy to generate plasma instabilities which are often kinetic in nature. The free energy sources associated with the deviation from the thermodynamical equilibrium distribution function could also excite plasma waves [3][4][5][6][7][8][9][10][11].
In general, the study of plasma waves and micro-instabilities in the solar wind shows that proton VDFs are prone to anisotropic instability and originate to be stable or marginally stable. Marsch [12] has discussed four significant electrostatic and electromagnetic wave modes and free energy sources to make them unstable. For example, the electrostatic ion acoustic wave may be destabilized by the ion beams and electrons and electron heat flux, [13] the electromagnetic ion Alfvén-cyclotron wave needs proton beam and temperature anisotropy, magnetosonic wave requires proton beam and ion differential streaming and whistler-mode and lower-hybrid wave [14] unstable solutions. Among several electromagnetic instabilities, the kinetic Alfvén wave instability is the most important one.
The satellite missions in space and astrophysical plasmas have confirmed the presence of non-Maxwellian high energy and velocity tails in the particle distribution function and found in the magnetosphere of Saturn, Mercury, Uranus and Earth [2,[15][16][17]. The non-Maxwellian distribution of charged particles has been observed to give a better fit to the thermal and superthermal part by employing kappa distribution, since it fits both thermal and suprathermal parts in the energy velocity spectra.
The subject area of this chapter involves the basic research of space plasma physics and in particular, focuses the investigations of electrostatic and electromagnetic waves in a multicomponent dusty (complex) Maxwellian and non-Maxwellian plasmas. In the last few years, various power-law distribution functions (in velocity space), i.e., kappa and r; q ðÞ have been used to investigate collective phenomena and associated instabilities, such as dust-acoustic waves, kinetic Alfvén waves, Weibel instabilities, dust charging processes (in linear and nonlinear regimes) in space and astrophysical situations for better fitting the observational data in comparison to Maxwell distribution. These distributions have relevance to space plasmas containing solar wind, interstellar medium, ionosphere, magnetosphere, auroral zones, mesosphere, lower thermosphere, etc.
When the intense radiations interact with plasmas, it ends up with many applications like instabilities, inertial confinement fusion [18], and pulsar emissions [19]. These instabilities further generate turbulent electromagnetic fields in plasma regimes. We can characterize instabilities as electrostatic as well as electromagnetic according to the conditions provided by nature [20]. In this chapter, we shall also discuss electromagnetic instability called Weibel instability in a Lorentzian plasma. The free energy source available for Weibel instability is temperature anisotropy and can be developed in magnetically confined and magnetic free plasma environment as well. First time Weibel [21] came up with the calculations of impulsive growing transverse waves with anisotropic velocity distribution function in 1958. This instability developed when the electrons in the fluctuating magnetic field generates momentum flux, this flux sequentially effects velocity < v > (and ultimately current density < J >)astoincreasethe fluctuating field [22]. The property of Weibel instability is that it is different from normal resonant wave-particle instabilities because it depends on effects in bulk plasma without any resonant particle contributions [23]. The particle distribution functions in kinetic model adequately describe a physical phenomenon in terms of time and phase space configurations providing more information to investigate plasma waves, instabilities, plasma equilibrium, Landau damping phenomena, etc. In this chapter, we shall review the kinetic/Inertial Alfvén waves and instabilities, the effects of dust grain charging as well as field aligned/cross field currents, streaming velocity and the non-Maxwellian power-law distribution and its effect on various electromagnetic modes. We intend to show that the presence of dust grains introduces a new cutoff frequency Ω dlh which is associated with the motion of mobile charged particles. Moreover, an interesting feature is to show that the employed model inhibits the temperature anisotropy and supports the velocity anisotropy. Further, we shall calculate the linear dispersion relation for Weibel instability in Lorentzian plasma B 0 ¼ 0; B 0 6 ¼ 0 ðÞ by using linearized, nonrelativistic Vlasov equation. We shall solve Z * κ α ðÞ by assuming α < 1orα > 1forκ ¼ 3, 5, 7.

Model and methodology
In long-wavelength modes the fluid theory is most reliable, while for short wavelength modes (like KAWs), collisionless effects are important, for example, Landau damping due to finite ion Larmor radius explains observed damping rate and in dusty plasmas and charge fluctuations. Kinetic Alfvén waves (KAWs) are small scale dispersive Alfvén waves (AWs) which plays a significant role in particle acceleration and plasma heating. A coupling mechanism between small-scale KAWs and large-scale AWs in the presence of superthermal particles has been discussed which in turns giving rise to the excitation of KAWS in a solar/stellar wind plasma have been studied in the past. In this chapter, we intend to show the relationship between the growth rates of excited anisotropic KAWs and perpendicular wavelength by taking charge fluctuation and Landau damping variations into account. Moreover, when the perpendicular component of the wavelength, when comparable to the ion gyroradius, a magnetic field aligned electric field plays a significant role in the plasma acceleration/heating. Utilizing a two potential theory along with kinetic description, the properties of kinetic Alfvén waves are aimed to investigate different modes in low beta plasmas by incorporating the streaming effects. We present overview of electromagnetic KAW streaming instability in a collisionless dusty magnetoplasma, whose constituents are the electrons, ions and negatively charged dust particles. The interaction between monochromatic electron/ion beam with plasma is also discussed under various conditions. Further, to calculate the linear dispersion relation for Weibel instability in unmagnetized Lorentzian plasma, we shall employ linearized, nonrelativistic Vlasov equation.

Two potential theory
In a low beta plasma, β < 1, the electric field can be described by two potential theory or fields expressing the electromagnetic perturbations with shear perturbations only in the magnetic field. We may neglect the electromagnetic wave compression along the direction of magnetic field (B 1z ¼ 0), which leads to the coupling of Alfvén-acoustic mode. Thus, we adopt a two potential theory which represents both the transverse and parallel components of the electric field as E ⊥ ¼À∇ ⊥ φ and E ∥ ¼À∂ψ=∂z, with φ 6 ¼ ψ [23]. We shall also consider the charge on the dust grain which may fluctuate according to the plasma conditions. At equilibrium, the charge neutrality imposes the condition n e0 À n i0 þ Z d0 n d0 ¼ 0, where n e0 n i0 ðÞ is the electron (ion) number density and Z d0 is the equilibrium dust charging state.
The linearized Poisson and Maxwell equations in terms of parallel and perpendicular operators can be expressed as and where E 0 μ 0 ÀÁ is the permittivity (permeability) of the free space and J j1z represents the field aligned current density for jth species j ¼ e ð for electrons, i for ions and d for dust grainsÞ.I n obtaining Eq. (2), we have ignored the factor ∇ ⊥ ∇ z E z À ∇ 2 E ⊥ : The main idea is to decouple the compressional Alfvén mode under the assumption ∇ Â E ðÞ z ¼À ∂B=∂t ðÞ z ¼ 0, i.e., to highlight bending of line of force and minimize any change in field strength due to wave compression. Moreover, ∇ Â E ðÞ z engross only E ⊥ and the least restrictive assumption for ∇ Â E ðÞ z to vanish is E ⊥ ¼À∇ ⊥ φ in which the perpendicular electric field E ⊥ is electrostatic, leaving an incompressible mode. When φ ¼ ψ, the twist of the magnetic field lines vanishes, therefore, the incompressible shear modes have ∇:u ⊥ ¼ 0 ¼ B 1z , and E ⊥ ¼À∇ ⊥ φ: where f j0 is the equilibrium distribution function. The dynamics of cold and magnetized dust is governed by set of fluid equations, i.e., and

Number density and current density perturbations
Here, we may define the number density and current density as Kinetic Theory

Dispersion and damping of kinetic Alfvén waves (KAW)
Kinetic Alfvén waves (KAWs) are small scale dispersive Alfvén waves (AWs) which plays a significant role in particle acceleration, turbulence, wave particle interaction and plasma heating. Kinetic processes prevail in the regimes where plasma is dilute, multi-component, and nonuniform. A coupling mechanism between small-scale KAWs and large-scale AWs with superthermal plasma species which in turns gives rise to the excitation of KAWS in a solar/stellar wind plasma has proved dispersive Alfvén waves responsible for the solar wind turbulence especially when the turbulence cascade of these electromagnetic waves transfer from larger to smaller scale as compared to proton gyro radius. Moreover, from spacecraft observations in ionospheric plasma, it is evident that Alfvénic Poynting flux is responsible to transfer the energy for particle acceleration. All the energized auroral particles accelerate in ionosphere, initiate Joule heating phenomenon and stream out into the magnetosphere [25][26][27][28].
There are number of studies to show the relationship between the growth rates of excited anisotropic KAWs and perpendicular wavelength by taking charge fluctuation and Landau damping variations into account. Moreover, the perpendicular component of wavelength, when comparable to ion gyroradius, a magnetic field aligned electric field plays a significant role in plasma acceleration/heating.
One of the important features in astrophysical plasma is the transportation of electromagnetic energy through the wave interaction with thermal plasma ions [29][30][31]. The KAW plays a vital role to transfer the wave energy through Landau damping (when thermal electrons travel along the magnetic field lines), which is regarded as collisionless damping of low-frequency waves and during this process the particles gain kinetic energy from the wave. This process can only happen when the distribution function has a negative slope which results in the heating of plasmas or acceleration of electrons along the magnetic field direction [24]. Recent studies also suggest the impact of non-Maxwellian distribution functions on the dynamics of solar wind and auroral plasma [32]. This study shows that the plateau formation in the parallel electron distribution functions minimize the Landau damping rate significantly.
In this chapter, the properties of kinetic Alfvén waves would be discussed by employing two potential theory, Maxwell equations and Vlasov model to study different plasma modes and by taking streaming of charged particles along and across the field direction in a Maxwellian and Lorentzian plasma.

Kinetic Alfvén waves in Maxwellian plasma
The propagation of kinetic Alfvén waves in a dusty plasma with finite Larmor radius effects will be discussed using a fluid-kinetic formulation by taking charge variations of dust particles. The coupling of Alfvén-acoustic mode results in the formations of kinetic Alfvén wave which would be discussed in forth coming subsections. In a magnetized plasma, we shall consider the electrons are thermal and strongly magnetized obeying an equilibrium Maxwellian distribution, while ions are hot and magnetized so that finite Larmor radius can be taken into account. For ions, we may employ Vlasov equation by utilizing guiding center approach to obtain the perturbed distribution function for an electromagnetic wave when the electric field and the wave vector k lie in the xz plane and, where J n is the Bessel function of first kind, having order n and f i0 is the equilibrium distribution function. On using Eq. (7), we obtain the modified number and current densities for hot and magnetized ions and thermal electrons, i.e., and where I n is the modified Bessel function with argument ϑ i, e ¼ k ⊥ v 2 ti, e =2Ω ci, e and Z ξ in ðÞ is the usual dispersion function for a Maxwellian plasma with ξ in ¼ ω À nΩ ci ðÞ =k z v ti and Z 0 is the derivative of Z with respect to its argument.
The dust component is considered to be cold and unmagnetized such that ω << ω cj, k z v te >> ω and k z v ti << ω, therefore, we use hydrodynamical model with momentum balance equation and continuity equation For cold and unmagnetized dust and thus we obtain and To find the relation between φ and ψ, the expressions of n e1 ,n i1 and n d1 are used into Eq. (1) and J i1z ,J e1z , and J d1z into Eq. (2) to obtain the following coupled equations: The coefficients in Eq. (15) are given by The solution of homogeneous Eq. (15) in the form of a biquadratic equation, i.e., where V Ai ¼ cω ci =ω pi is the Alfvén velocity of ions. The solution of biquadratic equation in the form of kinetic Alfvén wave is as follows, where, T 0 ¼ T e T i : and Λ i ¼ n i0 =n e0 This shows the dispersion relation of kinetic Alfvén waves in the presence of mobile dust that are the extension of shear Alfvén waves in the range of small perpendicular wavelength. The first term on the R. H. S appears due to dust dynamics, i.e., a new cut off frequency due to the hybrid dynamics of cold dust and magnetized ions which provides a limit to the propagation of electromagnetic wave. In a dustless plasma, i.e., ω pd ¼ 0, we obtain usual dispersion relation in electron-ion plasma. Expressing ω in terms of real and imaginary part, ω ¼ ω r þ iγ, with ω r >> γ, we either obtain growth or damping of KAW satisfying the condition, ω=k z ¼ v A ≤ v z through wave particle interactions [33,34]. In a dusty plasma with dust charge fluctuation effect, the main mechanism of wave damping is associated with dust charge fluctuation effects as compared to Landau damping [34]. It is a wellknown fact that if the particle thermal velocity exceeds the Alfvén velocity, then the particles interact with Alfvén wave as the result of wave particle interaction/resonance, the linear Landau damping prevails. In a dusty plasma, the massive dust grains move slowly as compared to Alfvén velocity, therefore they may interact with Alfvén wave through linear Landau damping (which is negligible in case of dust species) or charge fluctuation effects.

Lorentzian distribution function
A number of processes in a space based plasma lead to the development of particle anisotropy through streaming or temperature and are responsible for plasma instabilities in collision-free plasma which are frequently kinetic in nature and their persistent features have been confirmed by many spacecraft measurements, e.g., the electron energy spectra and the near-earth environment observations have witnessed the presence of superthermal populations. It is a well-known fact that the equilibrium Maxwell-Boltzmann distributions are associated with the Boltzmann collision term, but on the large scale Fokker-Plank model is not appropriate due to strong interaction and correlation in a collisionless plasma. The kinetic foundations of generalized Lorentzian statistical mechanics has been remarkably established by [35] with the generalization of Boltzmann collision term that is not based on binary collisions. The long range correlation between particles vindicates that power law distributions posses a particular thermodynamical equilibrium state. The mathematical form on isotropic Lorentzian distribution function is given by where A κ ¼ n i0 Due to the stated fact, the deviation from the Maxwellian equilibrium distribution function could also excite plasma waves by using free energy sources. Such distributions are frequently observed in solar and terrestrial environments and can be represented by anisotropy in temperature and velocity, i.e., [36] where is the thermal speed of jth plasma component, the number densities are represented by n and anisotropic temperatures components are represented as moment of second order In the limit κ ! ∞, the bi-Lorentzian function is reduced to bi-Maxwellian, f κ v ðÞ!f M v ðÞ :

Lorentzian current and number density perturbations
Many space and astrophysical plasmas have been found to have generalized Lorentzian particle distribution functions. It is of some interest to observe the impact of the high energy tail on the current and number densities of plasma species. By using Eqs. (4), (7) and (20), we get the modified expressions of number and current densities based on kappa distribution function, i.e., and where Z κ ξ j0 ÀÁ ðÞ κþ1 , is the plasma dispersion function and κ 0 ¼ 2κ À 1 ðÞ =2κ:.

KAW and instability in Lorentzian plasma
In a low β plasma, the kinetic Alfvén wave instability driven by field aligned currents has dependence on plasma β and streaming velocity of current carrying species which can be responsible for particle energization. In this subsection, we extend the above scenario of electromagnetic kinetic Alfvén wave by introducing the streaming of Lorentzian ions along an external magnetic field B 0ẑ ðÞ with constant ion drift velocity V 0 ∥B 0 ðÞ , strongly magnetized and hot electrons to be Maxwellian and cold unmagnetized dust. The plasma beta β e is assumed to be very small. The electric field and the wave vector k lie in the xz plane, i.e., B 0 ¼ 0; 0; B 0ẑ ðÞ , V 0 ¼ 0; 0; V 0ẑ ðÞ , k ¼ k ⊥x ; 0; k zẑ ðÞ : We again solve the Vlasov equation For hot and magnetized electrons [33] to get the number and current density of electrons as obtained in the previous section. Making use of Eqs. (4), (7) and (21), we get the perturbed number density of Lorentzian type streaming ions, The longitudinal components of current density perturbation [7,19,37] is given by where η ¼ ω À k z V 0 ðÞ =k z v tiκ : By incorporating the values of n i1 and J i1 in Eqs. (1), (2) and using (15), the dispersion relation of KAW streaming instability in a Lorentzian dusty plasma is obtained as where K ¼ k 2 ⊥ =k 2 z and A visible modification can be noticed by the effect of superthermality via the kappa-modified plasma dispersion function and the appearance of dust lower hybrid frequency due to dust effects on the dispersion characteristics. Numerous standard wave modes can originate from the above dispersion equation by applying particular limits, i.e., (i) k∥B 0 ; V 0 ðÞ : For n ¼ 0, ϑ e ≪ 1, a dispersion relation two stream instability (TSI) in unmagnetized plasma is obtained [37], i.e., In the limit κ ! ∞, our results approach to its classical Maxwellian counterpart in a dustless plasma environment [38].
(ii) k∥B 0 ðÞ ,V 0 ¼ 0, Ω ci ≪ ω ≪ Ω ce ðÞ : In a dustless plasma, we get whistler-like mode whose frequency is below the electron cyclotron frequency, i.e., where v ph ¼ c 2 k z Ω ce =ω 2 pe is the phase velocity of whistler waves which is obviously not susceptible to the Lorentzian index κ. Again, in the limiting case ω ≪ Ω ce , ϑ e ≪ 1, and expanding plasma dispersion function Eq. (15) depicts the coupling of electromagnetic and electrostatic mode, i.e., shear Alfvén-acoustic mode due to thermal kinetic effects due to which shear Alfvén wave builds a longitudinal component, e.g., where and c κ ¼ 2κ À 1 ðÞ = 2κ À 3 ðÞ : In the limit ω 2 =k 2 z V 2 A ! 0,V 0 ¼ 0,k 2 ⊥ r 2 e ≪ 1 we get the dispersion relation of Lorentzian dust-acoustic waves, 1=2 is the electron Alfvén speed with electron mass density. For a low beta plasma, the coupling between dustacoustic and shear Alfvén wave becomes weak and two modes would decouple. In the limit κ ! ∞, we approach to a Maxwellian DAW [40]. It is worthy to mention here that due to the contribution of Lorentzian particles, the KAW instability suppresses. As a matter of fact, the coupling mechanism enhances the unstable regions as the wave exchanges the energy, and we can deduce that in case of generalized Lorentzian plasma, the coupling between two modes becomes weak to some extent. Moreover for non-zero streaming velocity of ions, the unstable regions tend to grow. After simplifying Eq. (30), we get the mixed shear Alfvén-acoustic mode, i.e., where r 2 e ¼ T i =m e Ω 2 ce , and V Ae ¼ cΩ ce =ω pe : In the limit β i ≪ 1 and for k 2 ⊥ r 2 e ≪ 1, the two modes decouple and we get, where E z ¼ ω 2 pd þ k 2 z c 2 and Ω 2 dlh ðÞ e ¼ ω 2 pd =Ϝ e is the dust lower hybrid frequency which arises due to the hybrid motion of magnetized electrons and unmagnetized dust grains and is referred as a cutoff frequency which gives rise to a limit for the propagation of electromagnetic waves in the presence of dust grains. For graphical representation, we have chosen parameters typical to space dusty environment, for example, we consider n i0 ¼ 10 À 10 4 cm À3 ,n d0 ¼ 10 À 10 À2 cm À3 , Z d0 ¼ 10 À 10 4 ,m d ¼ 10 5 À 10 8 m i . For computational convenience, we introduce dimensionless parameters which are as follows: It has been observed that the growth rates of KAW instability are significantly affected by the presence of superthermal population, i.e., instability suppresses due to energetic particles possessed by kappa distribution when compared to its Maxwellian counterpart as shown in Figure 1. Similarly, the effect of streaming velocity, dust number density and charge on the growth rates is depicted in Figures 2-4 respectively. The free energy is associated with the drift motion of ions along the field direction which is responsible for the excitation of KAW. In a streaming plasma the velocity of ions is directly coupled to dust-acoustic waves and through this coupling the maximum growth rate is obtained when the wave exchanges energy through the streaming of ions. Moreover, the presence of dust particles has a noticeable effect on the wave dynamics through dust charge Z d and number density n d , i.e., it modifies the wave propagation and excitation. We can observe that Z d and n d enhances the growth rates of KAW due to the reason that when dust concentration in plasma is introduced they attach the plasma electrons toward them and the electron loss rate increases which in particular enhances the drift velocity to facilitate the unstable wave structure.

Dust kinetic Alfvén waves (DKAWs)
DKAWs arise when the dispersion relation of ordinary Alfvén waves is modified by the finite Larmor radius effect of dust. This process is dominated by the collective dynamics of magnetized dust particles. We have investigated shear Alfvén waves and their coupling with dust-acoustic wave by considering magnetized dust and Lorentzian electrons and ions.
The perturbed current and number densities of cold and magnetized dust are obtained by using Eqs. (5) and (6) The parallel component of perturbed dust current density turns out to be from Eq. (14) J d1z ¼ ω 2 pd =ω E 0 k 2 z ψ and the dispersion relation of kinetic Alfvén wave in the presence of magnetized dust is given by  Kinetic Theory In the limit κ ! ∞, we obtain classical results in a Maxwellian plasma.

Lorentzian-type charging currents
The charging equation containing Lorentzian electron and ion currents is where the electron and ion currents are calculated using a surface integral through the dust grain surface of radius r d having potential φ d are given as and where Te, i 4πn e, i0 e 2 1 2 is the Debye wavelength in superthermal plasma which is much smaller than found for a Maxwellian plasma and has been shown by [39,41] and a d is the radius of dust grain.
The Lorentzian charging currents are derived by using Vlasov-kinetic model whose fluid version by Rubab and Murtaza [41] and in the limit κ ! ∞, our results matched with Das et al., [32]. Now, by putting the value of perturbed dust grain charge, Q d1 ¼AE i ω Ωψ, in Eq. (1), the dispersion equation of DKAW becomes, which clearly shows that charge fluctuation effects are insensitive to the form of the distribution function.

Modified dust-acoustic wave
In the limit ω 2 =k 2 z V 2 A ! 0, the Eq. (16) after simplification turns out to be where ðÞ =n i0 n eo : Eq. (39) is the dispersion relation of dust-acoustic wave in a magnetized plasma whose Maxwellian version without dust charge fluctuation effects is given by Mahmood and Saleem [42]. It could be seen that the component of dust velocity in the direction of magnetic field V dz ðÞ , which finally turns out to be dust gyroradius, is responsible for the coupling of Lorentzian type DKAW with DAW. When the dust-acoustic wave frequency is very large compared to the dust gyroradius, then the dust is considered to be unmagnetized. In an unmagnetized plasma B 0 ¼ 0 ðÞ with T d ¼ 0, we get the dispersion relation of Lorentzian dust-acoustic wave (without dust charge fluctuation effects) which is exactly equal as discussed by [40]. The effect of Lorentzian index when growth rates are plotted as function of parallel and perpendicular wave number are depicted through graphical representation in Figure 5 and Figure 6 and shows that Maxwellian distribution functions are supportive to enhance the wave frequency.

DKAW: Perpendicular streaming
We consider an electromagnetic dust kinetic Alfvén wave streaming instability in a collisionless electron-ion dusty magnetoplasma. The motion of DKAW is followed by considering thermal and magnetized Lorentzian electrons to be Maxwellian and Lorentzian ions drifting across the external magnetic field B 0 ∥ẑ ðÞ with a constant drift velocity V 0x , i.e., V 0 ⊥B 0 ðÞ : The dust is considered to be cold and magnetized ω ≪ ω cd ðÞ and the charge on the dust grain surface is taken to be constant. The wave vector associated with the electromagnetic wave lies within xz plane where The distribution function of Lorentzian ions where ions are streaming perpendicular to the field direction is given as, where η 0 ¼ ω À k ⊥ V 0 ðÞ k ⊥ v tiκ : As there are no ions along the field direction due to perpendicular streaming, therefore we may neglect the ions current density J i1z ¼ 0: In the limit κ ! ∞, our results reduce to Maxwellian distribution.
The dispersion relation with the aid of Eq. (15) is obtained by using Eqs. (23), (24), (33) and (41) in Eqs. (1) and (2), i.e., which is the general dispersion relation of kinetic Alfvén waves in the presence of perpendicular streaming ions and cold and magnetized dust. In the above equation, Ϝ D ¼ ω 2 pd =ω 2 cd is responsible for the magnetized dust part.
For parallel propagation and in the limit ω 2 pd =c 2 k 2 ⊥ ≪ 1, Ϝ D ≪ 1, we get dispersion relation of two stream instability (TSI) in an unmagnetized dusty plasma. In a dust free plasma where Ω 2 dlh ðÞ d ¼ ω pi =Ϝ D , Λ d ¼ n d0 z d0 =n e0 and r d ¼ C D =ω cd : In the limit k 2 ⊥ r 2 d ≪ 1, we obtain modified shear Alfvén wave associated with the hybrid dynamics of the ions and magnetized dust through Ω 2 dlh ðÞ d which provides a cut-off for the EM wave propagation, i.e., where λ i ¼ c=ω pi : The dispersion relation for the DKAW instability is found to be dependent on the spectral index κ which means Lorentzian plasma is able to support a number of unstable branches. Lorentzian index is found to be more effective in large wave length limit as compared to small wavelength where the tail of unstable region remains independent of κ. When a large number of dust grains are introduced, it will enhance the loss rate of electrons by attachment on a dust grain surface which reduces the wave activity. At the same time the electron loss rate increases the drift velocity which in turns helps to excite the DKAWs and a further increase will help to stabilize the system. Due to particular choice of equations which involves parallel current density, the ions electromagnetic response cant not take part which limits the existence of ions Weibel instability.
By using the same parameters as above, we have plotted the growth rates as the function of propagation vector for different values of kappa. We have seen that the cross-field streaming of superthermal ions inhibit the growth rate of instability as shown in Figure 7. Similarly, β d is found to support the unstable structure and the instability increases with the value of β d as shown in Figure 8.

Weibel instability in a Lorentzian plasma
The Weibel plasma instability has so many applications in astrophysical [43], and in laboratory plasmas as well [44]. The generation of magnetic field can be explained in the domain of gamma ray burst, galactic cosmic rays and supernovae [45,48]. For the case of unmagnetized plasma, the Weibel instability [20] has been widely discussed in relativistic and nonrelativistic regimes. In 1989, Yoon [46,47] generalized his work by using relativistic bi-Maxwellian plasma. Later, Schaerfer [48] have discussed this instability in relativistic regimes of plasma with arbitrary distributions and presented comparison with his previous works which was based on bi-Gaussian distribution functions. The Weibel instability was investigated by Califano [49,50] with temperature anisotropy, produced by two counterstreaming electron populations. Davidson probed the multi species Weibel instability for the charged beam and intense ions in plasma [51].
In our work, we have derived the analytical expressions and compared the results numerically for the real and imaginary parts of the dielectric constant with the Maxwellian and kappa κ ðÞ distributions under two conditions i.e., α ¼ ωΘ kz ≫ 1 and ≪ 1. By using kinetic model, the linear dispersion relation for Weibel instability in unmagnetized plasma has been derived after solving the linearized, nonrelativistic Vlasov equation as below [52], where f 0 is the distribution function and here we will discuss the different velocity distributions, i.e., Maxwellian distribution and κÀ distribution functions.

Kinetic Theory 42
To calculate Weibel instability in a Lorentzian plasma, we use Eq. (21), for zero streaming velocity of particle, i.e., V 0 ¼ 0 and using ∂f 0 ∂v z in Eq. (45), and performing perpendicular integration, we are left with parallel integral which is called modified plasma dispersion function for kappa distribution.
where x ¼ Θ À1 z v z , α ¼ ωΘ z kz and Θ z, ⊥ ¼ 2κÀ3 ðÞ T z, ⊥ κm : Applying same procedure as above and again using Plemelj's formula, the integration of principal part yields The dispersion relation will be solved under two following conditions For α > 1 ðÞ : which shows the real part of Weibel instability is insensitive to the value of Lorentzian index and the imaginary part iπ 1 þ α 2 κ Àκ ! 0: For α < 1 ðÞ : The dispersion relation takes the form Now, we define a new plasma dispersion function, i.e., Z * κ α ðÞ ¼ 1 ffiffiffi ffi π p Γ κ ðÞ the corresponding dispersion relation can be expressed as is shown to play a vital role on the wave dynamics. Moreover, the dust grain charging yield some additional plasma currents, which depends on the streaming velocity, Lorentzian index and plasma beta. The Lorentzian index is found to either enhance or quench the electromagnetic instabilities. The dust component is found to play an essential role in wave dynamics, i.e., introducing dust lower hybrid frequency when mobile dust particles are included in the plasma. We have seen that the temperature anisotropy in the distribution function has no effect on the wave characteristics, i.e., the employed model inhibits the temperature anisotropy, but supports the velocity anisotropy. Moreover, a brief analysis on Weibel instabilities in a non-Maxwellian plasma in is also presented.
Kinetic Alfvén turbulence are always present in the streaming solar wind near 1 AU and in situ measurements have confirmed the presence of non-Maxwellian proton distribution function. The present investigations show that the Lorentzian charged particle distributions in space lead to a essentially new physical situation as compared to the plasma with equilibrium distribution functions. Our results of the present analysis opens a new window of investigation to study various streaming and anisotropic modes in different plasma scenarios when Lorentzian distribution function is employed.