Two-dimensional Plasma Research Articles

Identification of Active Magnetic Reconnection Using Magnetic Flux Transport in Plasma Turbulence

Abstract Magnetic reconnection has been suggested to play an important role in the dynamics and energetics of plasma turbulence by spacecraft observations, simulations, and theory over the past two decades, and recently, by magnetosheath observations of MMS. A new method based on magnetic flux transport (MFT) has been developed to identify reconnection activity in turbulent plasmas. This method is applied to a gyrokinetic simulation of two-dimensional (2D) plasma turbulence. Results on the identification of three active reconnection X-points are reported. The first two X-points have developed bidirectional electron outflow jets. Beyond the category of electron-only reconnection, the third X-point does not have bidirectional electron outflow jets because the flow is modified by turbulence. In all cases, this method successfully identifies active reconnection through clear inward and outward flux transport around the X-points. This transport pattern defines reconnection and produces a new quadrupolar structure in the divergence of MFT. This method is expected to be applicable to spacecraft missions such as MMS, Parker Solar Probe, and Solar Orbiter.

Self-organized pattern formation in radio frequency capacitively coupled discharges

Self-organized structures that are spread perpendicular to the radio frequency (RF) current direction have been observed in low temperature RF capacitively coupled plasmas. A fluid plasma model that includes thermoelectric electron energy transport is used to understand how these structures form. The electron thermoelectric transport coefficient is calculated using Bolsig+ for different chemistries and is found to be large for Ar plasma. Thermoelectric electron energy transport, which is driven by particle diffusion, opposes electron thermal conduction and can localize the plasma, leading to periodic structures. To examine these structures in radio frequency (RF) capacitive plasmas, two-dimensional Ar plasma at 13.5 MHz is first simulated without and then with thermoelectric electron energy transport. The charged species densities are perturbed in the simulations, and the growth or decay of different modes with time is observed. The periodicity of the structure is found to be determined by the relative strength of thermoelectric electron energy transport compared to energy conduction and losses. The effect of operating variables such as chemistry and pressure and design variables such as inter-electrode gap and steps in the electrode have been studied. For Ar plasma as pressure is decreased, the plasma peaks become stronger since thermoelectric electron energy transport is enhanced. Within limits, steps in the electrodes can be used to control the location of the periodic structures. For N2 plasma, the periodic structure does not appear as thermoelectric electron energy transport is weak. The spacing between plasma peaks is found to be dependent on pressure, chemistry, and inter-electrode gap.

Relaxation of Plasma Excitations in Two-Dimensional Electron Systems

An investigation is performed of the relaxation of plasma excitations in two-dimensional electron systems (2D-ESes) in GaAs/AlGaAs heterostructures. The effect radiational and incoherent collisional relaxation mechanisms have on the linewidth of cyclotron magnetoplasma resonance (CR) is studied. CR arises as a pure resonance that does not hybridize with dimensional magnetoplasma excitations. It is shown how magnetoplasma resonances form a CR’s fine structure due to the interaction between coherent radiative and incoherent collisional mechanisms of two-dimensional plasma relaxation. A comparative analysis is performed of the dependences of cyclotron and transport periods of relaxation on both the temperature and electron density of 2D-ESes. It is demonstrated that the period of cyclotron relaxation could exceed that of transport at low electron densities.

Determination of viscosity in shear-induced melting two-dimensional dusty plasmas using Green-Kubo relation

Langevin dynamical simulations of shear-induced melting two-dimensional (2D) dusty plasmas are performed to study the determination of the shear viscosity of this system. It is found that the viscosity calculated from the Green-Kubo relation, after removing the drift motion, well agrees with the viscosity definition, i.e., the ratio of the shear stress to the shear rate in the sheared region, even the shear rate is magnified ten times higher than that in experiments. The behaviors of shear stress and its autocorrelation function of shear-induced melting 2D dusty plasmas are compared with those of uniform liquids at the same temperatures, leading to the conclusion that the Green-Kubo relation is still applicable to determine the viscosity for shear-induced melting dusty plasmas.

Theoretical analysis of injection driven thermal light emitters based on graphene encapsulated by hexagonal boron nitride

We develop the device model for the proposed injection (electrically) driven thermal light emitters (IDLEs) based on the vertical hexagonal boron nitride layer/graphene layer/ hexagonal boron nitride layer (hBNL/GL/hBNL) heterostructures and analyze their dynamic response. The operation of the IDLEs is associated with the light emission of the hot two-dimensional electron-hole plasma (2DEHP) generated in the GL by both the lateral injection from the side contacts and the vertical injection through the hBNL (combined injection) heating the 2DEHP. The temporal variation of the injection current results in the variation of the carrier effective temperature and their density in the GL leading to the modulation of the output light. We determine the mechanisms limiting the IDLE efficiency and the maximum light modulation frequency. A large difference between the carrier and lattice temperatures of the IDLEs with an effective heat removal enables a fairly large modulation depth at the modulation frequencies about dozen of GHz in contrast to the standard incandescent lamps. We compare the IDLEs with the combined injection under consideration and IDLEs using the carrier Joule heating by lateral current. The obtained results can be used for the IDLE optimization.

Laser-induced melting of two-dimensional dusty plasma system in RF discharge

We present a detailed analysis of experimental study, which shows clear evidence of a two-stage melting process of a quasi-two-dimensional dusty plasma system in a high-frequency gas discharge. We accurately calculated global parameters of the orientational and translational order, as well as their susceptibilities to determine two critical points, related to “solid-to-hexatic” and “hexatic-to-liquid” phase transitions. The nature of the emerging defects and changes in their mutual concentration, in addition to the estimate of core energy of free dislocations also counts in favor of the formation of an intermediate hexatic phase. These results are fully consistent with the Berezinsky–Kosterlitz–Thouless theory.

Head-on collision of compressional shocks in two-dimensional Yukawa systems.

The head-on collision of compressional shocks in two-dimensional dusty plasmas is investigated using both molecular dynamical and Langevin simulations. Two compressional shocks are generated from the inward compressional boundaries in simulations. It is found that, during the collision of shocks, there is a generally existing time delay of shocks τ, which diminishes monotonically with the increasing compressional speed of boundaries, corresponding to the time resolution of the studied system. Dispersive shock waves (DSWs) are generated around the shock front for some conditions. It is also found that the period of the DSW decreases monotonically with the inward compressional speed of boundaries, more substantially than the time delay of shocks τ. When the inward compressional speed of boundaries increases further, the DSWs gradually vanish. We speculate that, for these high compressional speeds of boundaries, the period of the DSW might be reduced to a comparable timescale of the time delay of shocks τ, i.e., the time resolution of our studied system, or even shorter, thus the DSW reasonably vanishes.

Dispersion relation of square lattice waves in a two-dimensional binary complex plasma

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.

Thermal conduction in two-dimensional complex plasma layers

A simple vibrational model of heat transfer in two-dimensional (2D) fluids relates the heat conductivity coefficient to the longitudinal and transverse sound velocities, specific heat, and mean interatomic separation. This model is demonstrated not to contradict the available experimental and numerical data on heat transfer in 2D complex plasma layers. Additionally, the heat conductivity coefficient of a 2D one-component plasma with a logarithmic interaction is evaluated.

The distribution of OH radicals in a pulsed DC He + H2O plasma jet with the presence of a dielectric material

The distribution and density of OH radicals in a pulsed DC He + H2O plasma jet are studied by a two-dimensional (2D) plasma jet model. This article shows the effect of adding H2O to the pure He plasma jet on the OH density distribution and the changes of the OH density under different H2O contents. The process of streamer propagation toward the dielectric surface and development of the streamer on the dielectric surface are discussed in this article. The total OH radicals are generated from two parts in a He + H2O atmospheric pressure plasma jet (APPJ). The first part is produced with the water in ambient air outside the tube, and the other part is generated with the water in the helium flow. At the time when the streamer approaches the dielectric surface, it is found that electron-impact dissociation of H2O and electron neutralization of H2O+ are the main reactions to generate OH species and three-body combination of OH and OH with He is the main reaction to consume OH species. The OH density shows a donut shape in a pure He APPJ, while the density distribution of OH radicals under the nozzle gradually becomes uniform with the increase in the H2O content added to the helium flow, but the maximum value is at the both ends of the streamer outside the tube. The electron density initially increases with the increase in the proportion of H2O added to the He plasma jet and then gradually decreases with the H2O ratio greater than 0.5%. The enhancement of Penning ionization with the increasing H2O content is the main reason for the increase in electron density, while the increasing consumption of electrons by the excitation energy loss causes the decrease in the electron density with the H2O ratio more than 0.5%. The increase in the H2O content will bring both positive and negative effects. On the one hand, the increase in the H2O content added to the He gas flow will promote the generation of OH species. On the other hand, the decrease in the electron density caused by the increasing excitation energy loss slows the rapid growth of the reaction rates of the two main reactions related to the production of OH species. Eventually, OH density reaches saturation with the H2O ratio at 1.5% under the combined effect mentioned above.

Propagation of positive discharges in an air bubble having an embedded water droplet

Low temperature plasmas generated inside gas bubbles immersed in water is an effective method of rapidly transferring plasma generated reactive species to the water for applications in biomedicine, agriculture and environment. Reactive species are generally produced in the gas phase plasma and then solvate into the liquid. The large surface-to-volume ratio (SVR) of the bubble accelerates this process. In generating bubbles in water, aerosols and droplets are also contained within the bubble. These droplets also have a large SVR and so can be rapidly plasma activated. However, the presence of the droplets can also impact the propagation of the plasma in the bubble. In this paper, results are discussed from computational and experimental investigations of the formation and evolution of discharges in an air bubble immersed in water with an embedded water droplet. The computations were performed with a two-dimensional plasma hydrodynamics model. Experiments were performed with a quasi-2D bubble apparatus. In bubbles having a droplet, a plasma filament typically bridges from the powered electrode to the droplet, and then from the droplet to the bubble surface. A surface-hugging streamer also occurs on the inner bubble surface and on the surface of the droplet. Both surface streamers result in part from surface charge accumulation and can dominate the formation of reactive species that transport into the droplet. Increasing droplet conductivity suppresses propagation of the surface discharge and leads to a lower density of aqueous reactive species. Increasing conductivity of the surrounding water does not change the overall structure of the discharge but does slightly elevate the discharge intensity. The size and shape of the embedded droplet can significantly affect the formation and propagation of the streamer.

A tunable narrow-band plasma photonic crystal filter based on bound state

Photonic crystals are widely used in a class of narrow-band frequency selective filter due to their excellent ability to control electromagnetic waves, in which the working frequency depends on the structural parameters of the point defect resonant cavity of the photonic crystal, and the introduction of some dispersive media into the cavity makes the filter adjustable. In general, this kind of cavity-filter is very sensitive to the parameter disturbance of the cavity, and the quality factor of the filter can be reduced significantly by material loss. On the other hand, some studies have shown that there may be bound states at the interface between two different photonic crystals, and the bound state is often accompanied by narrow band and high transmittance, which implies that a narrow-band filter based on bound states is feasible. Importantly, filters based on bound states may be resistant to material loss to some degree. In this paper, a bound state related tunable narrow-band filter composed of a one-dimensional photonic crystal and a two-dimensional plasma photonic crystal is proposed, and the working frequency of the filter is located in the common band gap of the two photonic crystals. The COMSOL Multiphysics finite element simulation software is used to study the influences of geometric parameters of the one-dimensional photonic crystal and plasma parameters on the performance of the filter. It is found that the closer to each other the center frequencies and depths of the two different forbidden bands are, the greater the peak transmittance of the filter, in which the center frequency dominates, will be. On the other hand, the working frequency of the filter is directly proportional to plasma density and inversely proportional to collision frequency. The quality factor of the filter first increases and then decreases with the increase of plasma density, and decreases with the increase of collision frequency. The peak transmittance of the filter first increases and then decreases with the increase of plasma density, and decreases with the increase of plasma collision frequency. Finally, with the increase of collision frequency, both the peak transmittance and the quality factor decrease slightly, which indicates that the filter has a certain resistance to plasma loss. We believe that this work is helpful in investigating some new plasmonic photonic crystal filters.

Structure optimization of carbon nanotube ionization sensor based on fluid model

Compared with the traditional ionization sensor, triple electrode ionization gas sensor based on carbon nanotubes features excellent performances of small size and low operation voltage, and plays an important role in developing smart grids and the ubiquitous Internet of Things. However, they exhibit the disadvantages of small collecting current and low sensitivity, and need to be optimized in structure. Based on the Townsend discharge mechanism and three governing equations of particle mass conservation, electron energy conservation and Poisson equations, a two-dimensional plasma discharge fluid simulation model of the sensor is established by using the COMSOL Multiphysics software. According to gas composition, component concentration, structure of the sensor and voltage applied to the electrodes, we determine the chemical reaction, reaction rate coefficient, field control equation, initial values and boundary conditions. In order to calculate the interior distribution of charged particles inside the boundary, the inter-electrode space region is meshed. Initial values of mesh are set, such as a dense mesh of nano-meter length around the electrode having micro-nano structure, and a rough mesh of micro-meter length in the other region of inter-electrode space. The initial values of time-step are set for each mesh grid, and the discharge model of each mesh grid is numerically solved by the finite volume method. The three linkage governing equations are discretized and calculated iteratively in time and space, and the mesh values and time-steps are adjusted to make the results converged. The electrostatic field distribution and the average collecting current density of the eight kinds of sensors are obtained in the background gas of nitrogen. The optimal sensor structure is obtained by comparing the electric field intensities and current densities under different structure parameters. The eight kinds of the sensors are prepared for experimentally verifying the optimization method. The optimal structure of the 7# sensor has the highest collecting current density in the eight kind structures, which is consistent with the simulation results, and proves the feasibility of the structure optimization method proposed in the paper. Based on the optimal structure, the sensors, respectively with the electrode spacings of 100 and 120 μm, are fabricated, and the response characteristics of NO/SO<sub>2</sub> mixture gas are obtained. Compared with other technologies, the sensor with optimal structure exhibits 1－2 orders of sensitivity higher than the others. The optimized triple electrode ionization sensor based on carbon nanotubes exhibits many potential applications.

Continuous and discontinuous transitions in the depinning of two-dimensional dusty plasmas on a one-dimensional periodic substrate.

Langevin dynamical simulations are performed to study the depinning dynamics of two-dimensional dusty plasmas on a one-dimensional periodic substrate. From the diagnostics of the sixfold coordinated particles P_{6} and the collective drift velocity V_{x}, three different states appear, which are the pinning, disordered plastic flow, and moving ordered states. It is found that the depth of the substrate is able to modulate the properties of the depinning phase transition, based on the results of P_{6} and V_{x}, as well as the observation of hysteresis of V_{x} while increasing and decreasing the driving force monotonically. When the depth of the substrate is shallow, there are two continuous phase transitions. When the potential well depth slightly increases, the phase transition from the pinned to the disordered plastic flow states is continuous; however, the phase transition from the disordered plastic flow to the moving ordered states is discontinuous. When the substrate is even deeper, the phase transition from the pinned to the disordered plastic flow states changes to discontinuous. When the depth of the substrate further increases, as the driving force increases, the pinned state changes to the moving ordered state directly, so that the disordered plastic flow state disappears completely.

The boundary density profile of a Coulomb droplet. Freezing at the edge

We revisit the problem of computing the boundary density profile of a droplet of two-dimensional one-component plasma (2D OCP) with logarithmic interaction between particles in a confining harmonic potential. At a sufficiently low temperature, but still in the liquid phase, the density exhibits oscillations as a function of the distance to the boundary of the droplet. We obtain the density profile numerically using Monte-Carlo simulations of the 2D OCP. We argue that the decay and period of those oscillations can be explained within a picture of the Wigner crystallization near the boundary, where the crystal is gradually melted with the increasing distance to the boundary.

Dynamics of two-dimensional binary dusty plasma in liquid-like state

Using a driven-dissipative Langevin dynamics simulation, the dynamics of two-dimensional binary dusty plasma in the liquid-like state has been investigated. The system consists of charge-asymmetric microparticles with an identical charge-mass ratio. To include the effect of particle size on the dynamics of the system, a modified Yukawa potential with consideration of the particle radius is employed to account for the electrostatic interaction between the charged particles. It is found that the mean squared displacement (MSD) of the system changes little when the ratio of particle radius to interparticle spacing is up to the order of magnitude of one tenth. By analyzing MSD at different time scales, coexistence of the distinct self-diffusion regimes for the two species is observed. The compositions of the two species together with the effective coupling strength of the system are responsible for the motion separation of the two species.

A forced Korteweg–de Vries model for nonlinear mixing of oscillations in a dusty plasma

Nonlinear mixing of oscillations in a dusty plasma due to the harmonic time varying modulation of a nonlinear compressional oscillation is analyzed using a simple mathematical model consisting of a forced Korteweg–de Vries equation. An exact analytical solution of this equation is found to exhibit nonlinear mixing in the system. The model solution can be usefully employed to predict the existence of nonlinear mixing of oscillations in a two-dimensional dusty plasma system of a particular experimental configuration.

Streamer propagation along a profiled dielectric surface

We investigate the propagation of positive streamers along a profiled dielectric surface in air at atmospheric pressure. Results from experiments and two-dimensional planar low-temperature plasma fluid simulations are presented and analysed. The test object consists of a disk-shaped high voltage electrode and a dielectric slab with 0.5 mm deep corrugations. The corrugated surface has a 47% larger surface area than the smooth reference surface. The experiments and simulations are performed at voltage levels that lead to either gap-bridging or arrested streamers. In both experiments and simulations, the streamers take a longer time to reach the ground electrode when propagating along the profiled surface than along the smooth reference surface. Also, arrested streamers stop closer to the high voltage electrode when a profiled surface is used. Streamers propagate closely along the surface profile in the simulations, which suggests that the observed surface profile effect is mainly a result of elongated streamer channels. Compared to the streamers propagating along the smooth surface, the elongated streamers on the profiled surface have less residual voltage at the streamer front to drive the streamer advancement.

Buckling of two-dimensional plasma crystals with nonreciprocal interactions.

Laboratory realizations of two-dimensional (2D) plasma crystals typically involve monodisperse microparticles confined into horizontal monolayers in radio-frequency (rf) plasma sheaths. This gives rise to the so-called plasma wakes beneath the microparticles. The presence of wakes renders the interactions in such systems nonreciprocal, a fact that can lead to a quite different behavior from the one expected for their reciprocal counterparts. Here we examine the buckling of a hexagonal 2D plasma crystal, occurring as the confinement strength is decreased, taking explicitly into account the nonreciprocity of the system via a well-established point-wake model. We observe that for a finite wake charge, the monolayer hexagonal crystal undergoes a transition first to a bilayer hexagonal structure, unrealizable in harmonically confined reciprocal Yukawa systems, and subsequently to a bilayer square structure. Our theoretical results are confirmed by molecular dynamics simulations for experimentally relevant parameters, indicating the potential of their observation in state-of-the-art experiments with 2D complex plasmas.

Thermal conductivity analysis of two-dimensional complex plasma liquids and crystals

A novel homogeneous perturbed non-equilibrium molecular dynamics (HPMD) scheme, proposed by Evan-Gillan, has been employed to calculate the thermal conductivity of two-dimensional (2D) complex plasma liquids and crystals (CPLCs). The thermal conductivity has been reported using an improved HPMD method under the influence of constant external perturbation with different system sizes (N) and combinations of plasma parameters (Γ, κ). The current HPMD scheme provides precise outcomes with fast convergence for small-to-large N effects over a complete range of (Γ, κ). Temperature scaling law is tested for 2D thermal conductivity with appropriate Einstein frequency and found excellent behaviors. New simulations show that the thermal conductivity of CPLCs depends on (Γ, κ) and N and a slightly decreasing behavior is noted for thermal conductivity with increasing Γ and N, but, overall, thermal conductivity becomes constant at intermediate-to-large Γ. The reported thermal conductivity obtained from present HPMD method, in the limit of low equilibrium perturbation, has established a reasonable agreement with that obtained from earlier known 2D numerical and experimental data. It is demonstrated that the present HPMD method is an alternative efficient tool to compute the thermal conductivity of 2D CPLCs and can be a suitable method for complete trends of complex plasmas.