Physics

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. These results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

Binary complex plasmas consist of microparticles of two different species and can form two-dimensional square lattices under certain conditions. The dispersion relations of the square lattice waves are derived for the longitudinal and transverse in-plane modes, assuming that the out-of-plane mode is suppressed by the strong vertical confinement. The results are compared with the spectra obtained in Langevin dynamics simulations. Furthermore, we investigate the dependence of the dispersion relation on the charge ratio and mass ratio of the two particle species.

In this paper we compare different theoretical approaches to describe the dispersion of collective modes in Yukawa fluids when the inter-particle coupling is relatively weak, so that kinetic and potential contributions to the dispersion relation compete. Thorough comparison with the results from molecular dymamics simulation allows us to conclude that in the regime investigated the best description is provided by the sum of the generalized excess bulk modulus and the Bohm-Gross kinetic term.

The physical foundations of the dissipation of energy and the associated heating in weakly collisional plasmas are poorly understood. Here, we compare and contrast several measures that have been used to characterize energy dissipation and kinetic-scale conversion in plasmas by means of a suite of kinetic numerical simulations describing both magnetic reconnection and decaying plasma turbulence. We adopt three different numerical codes that can also include interparticle collisions: the fully kinetic particle-in-cell VPIC, the fully kinetic continuum Gkeyll, and the Eulerian Hybrid Vlasov-Maxwell (HVM) code. We differentiate between (i) four energy-based parameters, whose definition is related to energy transfer in a fluid description of a plasma, and (ii) four distribution function-based parameters, requiring knowledge of the particle velocity distribution function. There is an overall agreement between the dissipation measures obtained in the PIC and continuum reconnection simulations, with slight differences due to the presence/absence of secondary islands in the two simulations. There are also many qualitative similarities between the signatures in the reconnection simulations and the self-consistent current sheets that form in turbulence, although the latter exhibits significant variations compared to the reconnection results. All the parameters confirm that dissipation occurs close to regions of intense magnetic stresses, thus exhibiting local correlation. The distribution function-based measures show a broader width compared to energy-based proxies, suggesting that energy transfer is co-localized at coherent structures, but can affect the particle distribution function in wider regions. The effect of interparticle collisions on these parameters is finally discussed.

Recently the drift-reduced Landau fluid six-field turbulence model within the BOUT++ framework has been upgraded. In particular, this new model employs a new normalization, adds a volumetric flux-driven source option, the Landau fluid closure for parallel heat flux and a Laplacian inversion solver which is able to capture n=0 axisymmetric mode evolution in realistic tokamak configurations. These improvements substantially extended model's capability to study a wider range of tokamak edge phenomena, and are essential to build a fully self-consistent edge turbulence model capable of both transient (e.g., ELM, disruption) and transport time-scale simulations.

The drift-Alfven instabilities in the magnetic field aligned (parallel) sheared flow of a finite beta plasma with comparable inhomogeneous ion temperature and homogeneous electron temperature are examined. The development of instabilities are quantitatively discussed on the basis of numerical solution of a set of equations for the electrostatic and electromagnetic potentials. It is found that the accounting for the electromagnetic ion kinetic response, which has been ignored usually in existing discussions of the drift-Alfven instabilities of a steady plasma, reveals new drift-Alfven instability driven by the coupled action of the ion temperature gradient, the flow velocity shear, and the ion Landau damping. The excited unstable waves have the phase velocities along the magnetic field comparable with the ion thermal velocity, and the growth rate comparable with the frequency.

Very high frequency (VHF) driven capacitive discharges are now being increasingly adopted for plasma-based materials processing due to their high processing rates and lower substrate damage. Past studies related to complex plasma dynamics and higher harmonics generation in such systems were limited to constant voltage/current conditions, whereas, industrial systems are mostly driven by constant power density sources. In the present study, using particle-in-cell (PIC) simulation, we explore the dynamics of collisionless symmetric capacitive discharges that is operated at constant power densities. Our focus is on the effect of the driving frequency on the discharge parameters like the electron density/temperature, the electron energy distribution function (EEDF), the ion energy distribution function (IEDF), and the generation of higher harmonics in the device. The simulations are performed for a driving frequency from 27.12-100 MHz in argon plasma at a gas pressure of 1 Pa and for two values of the power density, namely, 2 kW/m3 and 20 kW/m3. It is observed that the required discharge voltage for maintaining constant power density decreases and discharge current increases with an increase in the driving frequency. A transition frequency is observed at both power densities. The density decreases (electron temperature increases) before the transition frequency and the trend is reversed after crossing the transition frequency. The EEDF shows an enhancement in the population of the mid-energy range of electrons as the driving frequency increases up to the transition frequency thereby changing the shape of EEDF from bi-Maxwellian to nearly Maxwellian, and then transforms into a nearly bi-Maxwellian at higher driving frequencies.

Analysis of driven dust vortex flow is presented in a weakly magnetized plasma. The 2D hydrodynamic model is applied to the confined dust cloud in a non-uniform magnetic field in order to recover the dust vortex flow driven in a conservative force field setup, in absence of any non-conservative fields or dust charge variation. Although the time independent electric and magnetic fields included in the analysis provide conservative forcing mechanisms, when the a drift based mechanism, recently observed in a dusty plasma experiment by [M. Puttscher and A. Melzer, Physics of Plasmas, 21,123704(2014)] is considered, the dust vortex flow solutions are shown to be recovered. We have examined the case where purely ambipolar electric field, generated by polarization produced by electron E*B drift, drives the dust flow. A sheared E*B drift flow is facilitated by the magnetic field gradient, driving the vortex flow in the absence of ion drag. The analytical stream-function solutions have been analyzed with varying magnetic field strength, its gradient and kinematic viscosity of the dust fluid. The effect of B field gradient is analyzed which contrasts that of E field gradient present in the plasma sheath.

A theoretical investigation of the modulational instability (MI) of dust-acoustic waves (DAWs) by deriving a nonlinear Schrödinger equation in an electron depleted opposite polarity dusty plasma system containing non-extensive positive ions has been presented. The conditions for MI of DAWs and formation of envelope solitons have been investigated. The sub-extensivity and super-extensivity of positive ions are seen to change the stable and unstable parametric regimes of DAWs. The addition of dust grains causes to change the width of both bright and dark envelope solitons. The findings of this study may be helpful to understand the nonlinear features of DAWs in Martian atmosphere, cometary tail, solar system, and in laboratory experiments, etc.

The standard nonlinear Schrödinger equation (NLSE) is one of the elegant equations to find the information about the modulational instability criteria of dust-ion-acoustic (DIA) waves (DIAWs) and associated DIA rogue waves (DIARWs) in a three-component dusty plasma medium having inertialess super-thermal kappa distributed electrons, and inertial warm positive ions and negative dust grains. It can be seen that under the consideration of inertial warm ions along with inertial negatively charged dust grains, the plasma system supports both fast and slow DIA modes. The charge state and number density of the ion and dust grain are responsible to change the instability conditions of the DIAWs and the configuration of DIARWs. These results are to be considered the cornerstone for explaining the real puzzles in space and laboratory dusty plasmas.