The plasma emission process, which is the mechanism for solar type II and type III radio bursts phenomena, is studied by means of particle-in-cell and weak turbulence simulation methods. By plasma emission, it is meant as a loose description of a series of processes, starting from the solar flare associated electron beam exciting Langmuir and ion-acoustic turbulence, and subsequent partial conversion of beam energy into the radiation energy by nonlinear processes. Particle-in-cell (PIC) simulation is rigorous but the method is computationally intense, and it is difficult to diagnose the results. Numerical solution of equations of weak turbulence (WT) theory, termed WT simulation, on the other hand, is efficient and naturally lends itself to diagnostics since various terms in the equation can be turned on or off. Nevertheless, WT theory is based upon a number of assumptions. It is, therefore, desirable to compare the two methods, which is carried out for the first time in the present paper with numerical solutions of the complete set of equations of the WT theory and with two-dimensional electromagnetic PIC simulation. Upon making quantitative comparisons it is found that WT theory is largely valid, although some discrepancies are also found. The present study also indicates that it requires large computational resources in order to accurately simulate the radiation emission processes, especially for low electron beam speeds. Findings from the present paper thus imply that both methods may be useful for the study of solar radio emissions as they are complementary.

We study the dispersion relation for parallel propagating ion-cyclotron (IC) waves in a dusty plasma, considering situations where the velocity dispersion along perpendicular direction is greater than along the parallel direction, and considering the use of product-bi-kappa (PBK) velocity distributions for the plasma particles. The results obtained by numerical solution of the dispersion relation, in a case with isotropic Maxwellian distributions for electrons and PBK distribution for ions, show the occurrence of the electromagnetic ion-cyclotron instability (EMIC), and show that the decrease in the kappa indexes of the PBK ion distribution leads to significant increase in the magnitude of the growth rates and in the range of wavenumber for which the instability occurs. On the other hand, for anisotropic Maxwellian distribution for ions and PBK distribution for electrons, the decrease of the kappa index in the PBK electron distribution contributes to reduce the growth rate of the EMIC instability, but the reduction effect is less pronounced than the increase obtained with ion PBK distribution with the same $\kappa$ index. The results obtained also show that, as a general rule, the presence of a dust population contributes to reduce the instability in magnitude of the growth rates and range, but that in the case of PBK ion distribution with small kappa indexes the instability may continue to occur for dust populations which would eliminate completely the instability in the case of bi-Maxwellian ion distributions. It has also been seen that the anisotropy due to the kappa indexes in the ion PBK distribution is not so efficient in producing the EMIC instability as the ratio of perpendicular and parallel ion temperatures, for equivalent value of the effective temperature.

In this paper we derive the dielectric tensor for a plasma containing particles described by an anisotropic superthermal (bi-kappa) velocity distribution function. The tensor components are written in terms of the two-variables kappa plasma special functions, recently defined by Gaelzer and Ziebell [Phys. Plasmas 23, 022110 (2016)]. We also obtain various new mathematical properties for these functions, which are useful for the analytical treatment, numerical implementation and evaluation of the functions and, consequently, of the dielectric tensor. The formalism developed here and in the previous paper provides a mathematical framework for the study of electromagnetic waves propagating at arbitrary angles and polarizations in a superthermal plasma.