Next Article in Journal
“Green” Extraction and On-Site Rapid Detection of Aflatoxin B1, Zearalenone and Deoxynivalenol in Corn, Rice and Peanut
Previous Article in Journal
Metabolite Analysis of Alternaria Mycotoxins by LC-MS/MS and Multiple Tools
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Isotropic and Anisotropic Monolayer Structures in RF Discharge Plasma

by
Anastasiya A. Alekseevskaya
1,
Elena V. Vasilieva
1,
Anatoly V. Filippov
1,2,
Mikhail M. Vasiliev
1,* and
Oleg F. Petrov
1
1
Joint Institute for High Temperatures, Russian Academy of Sciences, Izhorskaya St. 13 Bld. 2, Moscow 125412, Russia
2
Troitsk Institute for Innovation and Fusion Research, Pushkovykh St., vl. 12, Troitsk 108840, Russia
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(7), 3259; https://doi.org/10.3390/molecules28073259
Submission received: 1 March 2023 / Revised: 31 March 2023 / Accepted: 3 April 2023 / Published: 6 April 2023
(This article belongs to the Section Physical Chemistry)

Abstract

:
We present the results of an experimental and analytical study of the structural and dynamic properties of a monolayer consisting of dust grains in an electrostatic trap in an RF discharge plasma. The possibility of forming a monolayer with an isotropic distribution for interparticle distance and kinetic energy of particles in the structure has been experimentally shown. Isotropy has crucial importance for the study of various processes in such systems, including the kinetics of phase transitions, the formation of directed flows, wave propagation, and others.

1. Introduction

Recently, there has been a surge of interest in active matter studies, including research of systems consisting of active Brownian particles in colloidal plasma. Active Brownian particles are particles capable of converting the energy of an external source into the energy of their own (not thermal) motion [1,2,3]. Colloidal plasma is an open dissipative system in which charged granules of micron and submicron sizes levitate in plasma and, due to strong electrostatic interactions, form ordered structures of various types depending on conditions: chains, vortices, self-oscillatory systems, crystalline-like structures, etc. Such systems can be studied at the kinetic level by observing and analyzing the motion of individual particles of the system. The fundamental property of such systems is their ability to self-organize [4], namely, the formation of collective motions (vortices [5], solitons [6], etc.), or evolution of structures, passing through transformations during non-equilibrium phase transitions [7,8,9]. To keep like-charged granules from scattering, a metal ring laying on the electrode is usually used, which creates an electrostatic trap (confinement). This confinement can significantly affect the properties of systems. For the first time, density inhomogeneity for systems of particles, whose interaction is described by a screened Coulomb potential was discussed in [10,11]. In [12], the influence of a parabolic trap on the radial distribution of the interparticle distance in the system was discussed. In a number of numerical and theoretical papers [12,13,14,15,16,17,18], the authors consider that radial inhomogeneity of a dust monolayer (namely, the interparticle distance at the periphery of the structure is greater than at its center) is a fundamental feature. Consequently, the term “dusty plasma crystal” is questioned as applied to ordered systems of dust particles in the plasma of gas discharges.
In this paper, we present the experimental results confirming the fact that the interparticle distance in a dusty monolayer can widely vary, depending on the parameters of the gas discharge. It can be heterogenic, both with a compact center and more rarefied periphery, as well as with a loose packing of particles in the center and a closely packed periphery. Moreover, it can be isotropic as well over the entire diameter of the structure under certain conditions. This fact is of crucial importance and should be taken into account when designing experimental studies of various processes in such systems, including the kinetics of phase transitions, the formation of directed flows, the development of an oscillatory process, etc. In this case, one should first set the discharge parameters to obtain a monolayer with isotropic properties, and then, without changing the discharge parameters, affect the properties of the structure by applying an external perturbation, for example, laser radiation in the case of kinetic heating of active Brownian macroparticles [19].

2. Results and Discussion

2.1. Experimental Results

Figure 1 shows a video frame of a monolayer structure of charged particles in RF discharge plasma. The evolution of the structure was experimentally observed while changing the discharge parameters. Based on the particle coordinates reconstructed from the video recording, we obtained the radial distribution of interparticle distances in the structure.
It was experimentally observed that a change of discharge power led to a change of interparticle distance under a constant value of plasma-forming gas pressure P = 5 Pa (see Figure 2).
At a power of W = 4.4 W, the interparticle distance in the central part of the structure appeared to be about 1.5 times less than the interparticle distance at the periphery of the structure. However, with increasing power, the dust monolayer changed and at W = 6.9 W, the interparticle distances in the structure leveled off. A further increase in the discharge power to W = 16.7 W led to the formation of a sparse area in the central part of the structure, where the interparticle distances were significantly greater than at the periphery of the structure (Figure 3). A further increase in power led to the formation of a void in the central part of the structure. The appearance of voids in dusty plasma structures was previously observed in [20,21,22,23].
A similar trend was observed with a fixed discharge power W = 13.1 W while varying argon pressure in the discharge chamber (Figure 4). At P = 3 Pa, the interparticle distance in the central part of the structure was about 2 times less than the interparticle distance at the periphery of the structure. With an increase in pressure to P = 4 Pa, the interparticle distances in the structure were equalized. A further increase in pressure to P = 5 Pa led to the formation of a sparse area in the central part of the structure (Figure 5).
In work [12], the effect of a parabolic trap on the confined dusty plasma structure of particles interacting by a screening Coulomb potential was theoretically investigated. The authors came to a conclusion, based on the simplified model of dusty plasma, that the structures should be always heterogeneous. However, one can observe various structure configurations, including homogeneous, in the laboratory experiment (see Figure 2 and Figure 4). In the next section, we present theoretical models which can explain the obtained experimental results.

2.2. Basic Equations of the Hydrodynamic Model for Dusty Plasma in a Capacitive RF Discharge

The electron component in the dust plasma of an RF discharge can be described by the equations of the balance of the particles’ number and electron energy in terms of a hydrodynamic approximation with the values, such as electron concentration, average energy, and drift velocity, averaged over the period of the RF field (the derivation of these equations is presented, for example, in [24,25,26]):
n e t + div j e = k i o n n e N J d n d
( n e ε ) t + div h e = e E p grad ( n e D e ) + e μ e n e E p 2 + 1 2 e μ e , Re n e E 0 2 n e N W e l k e x n e N E e x k i o n n e N E i o n 1 3 ε J d n d ,
where ne is the electron number density, je is the electron flux density in the drift-diffusion approximation defined by
j e = grad ( n e D e ) μ e n e E p ,
De, μe are the electron diffusion coefficient and the mobility, respectively, Ep is the strength of the polarization constant field of the electrode sheaths, kion is the total gas ionization coefficient, Jd is the total sink of electrons and ions on a dust particle, nd is the number density of dust particles, the N is the number density of atoms, 〈ε〉 is the mean electron energy, he is the electron energy flux density defined as
h e = grad ( n e D ε ) μ ε n e E p
μe,Re is the real part of mobility, in phase with the RF field, με is the energy mobility, Dε is the energy diffusion coefficient, Wel is the total net energy loss rate due to elastic electron-neutral collisions, Eion is the ionization potential, kex is the atom excitation coefficient, Eex is the atom excitation potential.
We added to electron number balance Equation (1) a term that takes into account the absorption of electrons by dust particles, where an electron flux Je and an equal ion flux Ji (since the characteristic charging times are much less than the characteristic times of the particle motion and the formation of dusty plasma structures) can be determined from orbital-motion-limited (OML) theory [27]:
J e = π a 2 n e ( 16 ε 3 π m e ) 1 / 2 exp ( 3 2 e ϕ d ε ) , J i = π a 2 n i ( 8 T i π m i ) 1 / 2 ( 1 e z i ϕ d T i ) ,
where Je, Ji—the number of plasma electrons and ions absorbed by a dust particle of radius a per unit time, ϕd—the potential of the particle surface, which can be found from the equality of the fluxes of electrons and ions:
J e ( ϕ d ) = J i ( ϕ d ) J d .
Other processes of loss (for example, recombination) and creation of electrons (for example, stepwise ionization with the participation of metastable excited atoms) can be neglected under conditions of low pressures and a low degree of ionization.
The electron transfer coefficients and constants for ionization and excitation rates were calculated using the BOLSIG+ program [28,29] with cross sections from the LXCat set [30]. Figure 6 and Figure 7 show the calculated dependences of the average electron energy on the specific field and the electronic coefficients as functions of the average electron energy, respectively. Note that in Figure 6  E = E 0 / 2 , where E0 is the RF field amplitude. We also note that at the pressures indicated in the captions to the figures, the values of the specific circular frequency of the RF field ω/N form a uniform grid.
Figure 6 shows that the average electron energy increases with pressure growth and the ratio of the transport frequency νm to the circular frequency of the RF field ω decreases (see Figure 7a). Note that the electron energy distribution function (EEDF) in a radio-frequency field when the condition ωνu is met (νu is the characteristic frequency of electron energy relaxation, ν u 2 m e ν m / M , me is the mass of the electron, M is the mass of the atoms of the working gas) is determined by the effective value of the field [26]:
E e f f = ν m ω 2 + ν m 2 E 0 2 .
Since the transport frequency depends on the energy of electrons, the effective field also depends on the energy of the electrons. The behavior of the average electron energy and the ratio νm/ω noted above is a consequence of the complex behavior of the cross section of the collision of electrons with argon atoms, which passes through the deep Ramsauer minimum [30]. It can be seen in Figure 7 that the diffusion coefficient De is a weak function of argon pressure (the values με, Dε, and Wel also had a similar behavior), and the electron mobility μe and the ionization rate constant kion noticeably change with pressure changes (the dependences μe,Re, kexc, and Pe were similar).
The polarization field of the near-electrode layers (so-called ambipolar field) is determined by the Poisson equation:
div E p = 4 π e ( z i n i n e + z d n d ) ,
where e is the proton charge, zi is the charge number of ions (hereafter is assumed that zi = 1), ni is the ion number density, zd is the dust particle charge in elementary charges which is connected with the dust surface potential by relation
z d = ϕ d a ( 1 + k D a ) ,
where kD is the inverse Debye radius.
The ion component is described by the equations of balance for ion number and flux density:
n i t + div j i = k i o n n e N J i n d ,
where ji is the ion flux density defined without allowance for the diffusion:
j i = μ i n i E p ,
μi is the ion mobility.
The experimental mobilities of the argon ion (see Figure 8) via reduced field under normal conditions at a temperature of 300 K [30,31,32,33,34,35,36] are well approximated by the dependence:
μ i = 0.286 + 0.669 exp ( E / N 179.5 ) + 0 . 679 exp ( E / N 1305 ) .
Here μi is valued in cm2/(V·s), and E/N—in Td.
Figure 8 shows only part of the literature data (more complete data are given in [32,33,34,35]). Note that in a number of papers (see, for example, [35]), to determine the ion drift velocity at high values of E/N, its relationship with the measured average ion energy Ei in the form [36,37] was used:
E i = π m i v i , d r 2 2 ,
where mi is the ion mass, vi,dr is the ion drift velocity. In our paper, we use the Wannier approximation [38,39],
E i = 3 2 T g + ( m i + m g ) v i , d r 2 2 ,
where mg is the gas atom mass and Tg is the gas temperature.
The dust component is described by the equations of the balance for the number of particles and the velocity of motion (see, for example, [40,41,42]):
n d t + div ( n d v d ) = 0 ; v d t = g + 1 m d ( e z d E p + F i d F n f )
where g is a free fall acceleration, Fid is an ion drag or the force, induced by ion flux focusing, Fnf is the gas resistance force, which is determined at low pressures under the assumption of complete accommodation of the momentum of the gas atoms by the Epstein formula [43]:
F n f = 4 3 π a 2 δ M N v g v d ,
vg—thermal velocity of gas atoms: v g = 8 T / π M , , δ ≈ 1 is a coefficient depending on the nature of the reflection of gas atoms from the surface of a dust particle after a collision.
In [27] (see also [44,45] and references therein), the expression for the ion drag force is the following:
F i d , φ = n i v s m i v i π ( b c 2 + 4 b L 2 Λ ) ,
where vTi, vs, Λ, bL, bc, RDe is thermal and total ion velocities, taking into account the drift motion, the Coulomb logarithm, the Landau radius, the ion capture radius, and the electron Debye radius, respectively, defined by the expressions:
v T i 2 = 8 T i π m i ,   v s 2 = v i 2 + v T i 2 ,   Λ = 1 2 ln ( R d 2 + b L 2 b c 2 + b L 2 ) , b L = e 2 | z d | m i v s 2 ,   b c = a ( 1 2 e ϕ d m i v s 2 ) ,   R D 2 = ε 6 π e 2 n e .
In this work, to study the dynamics of dust particles, we conducted simulations in a one-dimensional approximation along the z axis, directed upward towards gravity. As a first approximation, we assumed that the electron transfer coefficients and the gas ionization rate constant were not changing. This was justified, since in the regime under consideration, the field in the positive column was weak, and the electron density in the near-cathode layers was small.
In the steady state, the gas velocity can be set equal to zero, since the gas heating rate in the considered low-current RF discharge glowing mode was small. Additionally, the characteristic time for the establishment of the velocity of dust particles significantly exceeded the characteristic time for the establishment of their concentration; therefore, the inertia of dust particles can be neglected and their velocity can be determined from the stationary equation of motion of dust particles.
Let us show that the ion drag force can be comparable in magnitude to the force of gravity. In this case, taking into account the ion drag force is of fundamental importance, since it is responsible for changing the distribution of dust particles along the radius of the RF discharge.

2.3. Analytical Solution of RF Discharge

There are many works in the scientific literature devoted to the analytical approximate solution of an RF discharge (see, for example, [46,47,48,49,50]), but all of them are not free from shortcomings and have a limited area of applicability, which either lies in the region of rare collisions, or in the region of frequent collisions. In our experiments, we observe the case when the parameter αc, which determines the collisional or collisionless nature of ion transport in the near-electrode layers, defined by the relation [48]
α c = π R D e l i ,
is valued as the order of unity (see Tables). Here, RDe is the electron Debye radius, li is the mean free path of ions. Therefore, we propose the following method for determining the parameters of an RF discharge under our conditions.
Let us consider the stationary solution of Equations (1), (6), and (8) for a symmetric RF discharge near the center of the discharge gap, where the polarization field is still sufficiently small and, therefore, dust particles cannot levitate. Let us assume that the transfer coefficients and the frequency of gas ionization in this region are constant. Equations (1), (6), and (8) for a one-dimensional problem take the form:
D e 2 n e z 2 μ e ( n e E p ) z = ν i o n n e , μ i ( n i E p ) z = ν i o n n e , E p z = 4 π e ( z i n i n e ) .
From Equation (18), after simple algebra we find
E p 2 z 2 = 8 π e [ ( 1 μ e + 1 μ i ) ν i o n n e + D e μ e 2 n e z 2 ] .
We choose the origin of the z axis at the center of the discharge gap. Now, we seek a solution for the electron density in the form (see, for example, [51]):
n e = n e 0 ( z ) cos ( α z ) ,
where α is an unknown constant. We are interested in a solution near the center of the discharge gap, so we assume that in this region ne0 is a weak function of the z coordinate and neglect its derivatives. Now from (19), using the discharge symmetry condition and the zero value of the field in the center, we find:
E p 2 = 8 π e [ 1 α 2 ( 1 μ e + 1 μ i ) ν i o n D e μ e ] n e 0 [ 1 cos ( α z ) ] .
From Equation (14), using (16), we can find:
n e E p = ( D e α μ e ν i o n α μ e ) n e 0 sin ( α z ) , n i E p = ν i o n α μ i n e 0 sin ( α z ) .
From here, for the ratio of the concentrations of electrons and ions near the center of the discharge, we obtain:
n e n i = α μ i ν i o n ( D e α μ e ν i o n α μ e ) = μ i μ e ( D e α 2 ν i o n 1 ) .
Since there must be ne ≈ ni, near the center of the RF discharge for the formation of a quasi-neutral plasma, then from (23) we find that for this equality the following expression must be satisfied:
D e α 2 ν i o n = 1 + μ e μ i .
In general, the value of α is not defined here, but it can be expected that it will be of the order of 2π/L, where L is the inter-electrode distance. In this case, the commonly used zero boundary conditions for the electron concentration at the electrodes will be satisfied. Then the condition for the formation of a quasi-neutral plasma at the center of the RF discharge will have the form:
ζ μ i μ e ( 4 π 2 D e ν i o n L 2 1 ) = 1 .
Note that although we have four equations for determining ne, ni, and Ep at the center of the discharge, but due to their approximate nature, they do not allow us to determine the value of the parameter α. Expression (25) allows one to determine the value of the reduced field E/N for which this condition is met. Further, the value of this field will be referred to as ERF.
Figure 9 shows the values of the parameter ζ as a function of the reduced field E/N at different values of the argon pressure. It can be seen from Figure 9 that the field required to ignite the RF discharge decreases with increasing pressure, which is associated with a decrease in the frequency of diffusive electron escape.
The power of the energy deposition into the RF discharge is the sum of the electronic and ionic parts, which can be written as [48]:
W d = W e + W i = 1 2 ( V p I + R s h I 2 ) ,
where Rsh is an ohmic resistance of the near-electrode layers, for which we use the relation RshSHF ≈ 104 Ω cm2 [48], SHF is an electrode area, Vp is the ohmic part of the RF voltage across the RF discharge plasma, which also takes into account the so-called stochastic heating of electrons due to their reflection in the near-electrode layers. At this stage, we neglect stochastic heating, since there are debates in the literature about the mechanism of this heating, the analytical model of which is developed in [46,47,48] (see, for example, [51,52,53,54,55,56]).
Since the thickness of the near-electrode layers is several electron screening radii, which is negligibly small compared to the length of the discharge gap L, we take for estimates
V p = E 0 L ( 1 + ω 2 ν m 2 ) 1 / 2 ,
Here E0 is amplitude of RF field. Then, knowing a value of the field E, found from the condition (25), and the experimental input power of the RF field Wd, one can determine the RF discharge current and obtain the following expression for the ratio of ion and electron heating:
η W i W e = 1 2 ( 1 + 4 R s h W d E 2 L 2 ( ν m 2 + ω 2 ) ν m 2 1 ) .
As a result, we derive the following equation to determine the electron number density:
n e = W d P e ( 1 + η ) ,
where Pe is the energy per unit time absorbed by the electrons from the electric field: Pe/N = (μe,Re N)(E/N)2. Table 1, Table 2, Table 3, Table 4 and Table 5 show the parameters of the plasma in the RF discharge and the dust component in the levitation region, where the condition is met:
e z d E p = m d g .
This field will be further denoted as Elev. The coefficients related to the electron component in Table 1, Table 2, Table 3 and Table 4 are determined by the value of the reduced E/N field found from (25), and the electron number density is determined from (29). It is assumed that the number density of ions is equal to the number density of electrons. The ion energy in the levitation region of dust particles El,ion is determined from the expression (12), while the drift velocity is determined using a field value equal to E 2 + E lev 2 . This value of ion energy is used to determine the potential of dust particles ϕd from the transcendental Equation (4) and their charge zd from the ratio (7).
Table 3 and Table 4 also show the values of the charge zd0 and the floating potential of dust particles ϕ0, determined with the effective ion temperature: T i = 2 3 E i , where Ei is the ion energy determined by Formula (12) taking into account only the drift motion in the RF field. It can be seen that this leads to a noticeable decrease in the charge of dust particles (in absolute magnitude). This is a consequence of the fact that the ion flow to a dust particle at high values of grain charge depends on the ion temperature as T i 1 / 2 , since the last multiplier in Equation (3) becomes of order of unit. As a result, when taking into account the ion drift only in the RF field, the effective temperature turns out to be lower, and therefore, the ion flux increases. In order to maintain its balance with the electron flux, the potential of the dust particle must decrease in absolute magnitude. As can be seen from Table 5, the electric polarization field in the levitation region of the dust particles significantly exceeds the intensity of the RF field.
Table 5 shows the specific values of the power of energy deposition and stochastic heating, the ion drag force, the transport frequency, and the ratios of the ion and electron plasma frequencies to the circular frequency of the high-frequency field in the experiments described above. It can be seen that the ion drag force is only slightly lower than the gravity force and should be taken into account when determining the levitation height of dust particles. These assumptions will be considered in our future work, together with other RF discharge parameters determined from the numerical solution of the system of Equations (1), (2), (6), (8), and (13). In addition, one can see from Table 5 that the electron plasma frequency turns out to be much higher than the circular frequency of the RF field in our experiments, while the ion frequency is much lower. These conditions are considered in the analytical theory, developed in [48].
The power of stochastic heating of electrons in Table 5 was determined according to [51]:
W s t = 3 8 n e v T t h , e e ( u ˜ e v t h , e ) 2 H G k ( H ) ,
where u ˜ e is the electron drift velocity amplitude, vth,e is the electron thermal velocity: v t h , e = 8 T e / π m e , s ˜ e is the amplitude of electron oscillations in the bulk plasma, H is the parameter defined by the relation H = s ˜ e 2 / π R D 2 , Gk(H) is an integral that is well approximated by [52]
G k ( H ) 24 H + 55 .
The amplitudes of the drift velocity and oscillations in the RF field were determined from the ratios:
u ˜ e = E ω 2 μ e . Re 2 + μ e . Im 2 , s ˜ e u ˜ e ω .
With the values of the screening radius found, the pair correlation function (PCF) of dust particles was calculated for different charges or effective coupling parameters [57]:
Γ * = e z 2 d e a d T d ( 1 + k D a d + 1 2 k D 2 a d 2 ) e k D a d ,
where ad is the mean interparticle distance of the dust particles: a d = n d 1 / 3 . Figure 10a shows the dependence of the PCF on the distance found by the numerical solution of the Ornstein–Zernike (OZ) iteration method [58,59] in the hypernetted chain (HNC) approximation. The figure caption indicates the charge value at Γ* = 240. The charge of the dust particles shown in Table 5 is an order of magnitude higher than this value. As can be seen from Figure 10, the height of the first peak increases with the growth of the charge, so a slight difference in the height of the experimental first peak may be due to the higher value of particle charge. Unfortunately, the iterative method of solving the OZ equation in the HNC approximation ceases to converge at higher values of the coupling parameter.
The determination of the fraction of the RF discharge power associated with the ionic component is approximate, so the pair correlation function is also calculated at zero ionic power η = 0 (see Figure 10b). In this case, the electron number density increases more than 6 times, while the Debye radius decreases and the value of the structural parameter λ = kDad grows from 1.7 to 4.4. Since in numerous experiments on the study of dusty plasma, the value of the structural parameter does not exceed two, we can say that taking into account the ionic component of the power gives a more consistent result when compared with the experimental data.
The dependences of the PCF on the distance shown In Figure 10 are similar to those determined experimentally, which allows us to conclude that the analytical approach developed in this work is applicable for estimating the parameters of dusty plasma in an RF discharge. As can be seen from Table 5, according to analytical estimates, the force of ion drag in the experimental conditions turns out to be comparable to gravity; therefore, it must be taken into account when finding conditions for the levitation of dust particles.

3. Experimental Setup

For the experimental study of quasi-two-dimensional structures of charged macroparticles formed in the RF discharge plasma, an experimental setup described in [19] was used (see Figure 11). The main element of the setup is a vacuum chamber with optical windows in the horizontal plane and in the upper part of the chamber. The atmosphere is pumped out from the chamber by a consequently connected fore-vacuum and turbomolecular pumps to a residual pressure of 10−4 Pa. Then, the chamber is filled with plasma-forming gas—argon, up to a working pressure of 3–5 Pa. In the center of the chamber, there are two flat horizontal electrodes, the distance between which is 5 cm. A voltage of 350 V with a frequency of 13.56 MHz is applied to them, as a result of which low-pressure plasma ignites in the chamber.
After that, melamine–formaldehyde (MF) grains with a diameter of d = 10.6 μm and a density of ρ = 1.51 g/cm3 are injected into the discharge. To create an electrostatic trap, a ring 8 cm in diameter and 0.2 cm in height is installed on the lower electrode, so that it keeps the charged dust grains from scattering in the radial direction. Being illuminated by a flat horizontal beam of an argon laser, the grains can be easily seen with naked eye. We use the CCD—camera placed above the chamber to record top view of the structure in dynamics. The obtained video data is processed then using a specialized program code, as a result of which various characteristics of dust particles can be obtained, such as their coordinates and motion trajectories, pair correlation functions, the coupling parameter, etc.
Two series of experiments with a monolayer structure were carried out. In the first one, we fixed the working gas pressure at 5 Pa and varied the discharge power from 4.4 to 16.7 W. In the second series, we fixed the discharge power at the level of 13.1 W and varied the argon pressure from 3 to 5 Pa. Then, we processed the corresponding video recordings in order to obtain the above-mentioned parameters of the dusty plasma structure, and therefore, to calculate the radial density distribution of dust particles depending on pressure and discharge power.

4. Conclusions

Our experiments unambiguously show that in the plasma of an RF discharge, it is possible to observe structures of charged grains not only rarefied at the periphery and dense in the center, but also with an isotropic interparticle distance, as well as more rarefied in the central part. This fact is of crucial importance and should be taken into account when designing experimental studies of various processes in such systems, including the kinetics of phase transitions, the formation of directed flows, the development of an oscillatory process, etc.
It is shown that by changing the discharge power or the pressure of the plasma-forming gas, one can influence the radial distribution of the interparticle distance in the dust monolayer and the uniformity of the structure.
Previously, based on the hydrodynamic model of dusty plasma in a capacitive RF discharge without taking into account the ion drag force, it was assumed that being placed in a parabolic trap, particles form a structure with a rarefied periphery. However, in this work, it is shown that the ion drag force can be comparable in magnitude to the force of gravity. The ion drag force is believed to be taken into account in the model, since it is responsible for changing the distribution of dust particles along the radius of the RF discharge. It is known that this force leads to the formation of voids in RF discharges, and in our experiments, the formation of voids also took place at certain plasma parameters.
Verification of the hypothesis about the influence of the ion drag force on the inhomogeneity in the distribution of dust particles requires more accurate calculations, which is planned to be done by self-consistent numerical simulation of dusty plasma in an RF discharge within the framework of the hydrodynamic model developed in this work.

Author Contributions

Conceptualization and methodology M.M.V. and O.F.P.; formal analysis and validation, A.V.F. and E.V.V.; investigation, A.A.A.; data curation and resources, M.M.V.; writing—original draft preparation, A.V.F., M.M.V. and A.A.A.; writing—review and editing, M.M.V. and E.V.V.; visualization, E.V.V., M.M.V. and A.V.F.; supervision, O.F.P.; project administration, M.M.V.; funding acquisition, O.F.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was carried out as part of JIHT RAS state assignment with funds from the Ministry of Science and Higher Education of the Russian Federation.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

Samples of the compounds are not available from the authors.

References

  1. Bechinger, C.; Di Leonardo, R.; Löwen, H.; Reichhardt, C.; Volpe, G.; Volpe, G. Active particles in complex and crowded environments. Rev. Mod. Phys. 2016, 88, 045006. [Google Scholar] [CrossRef]
  2. Elgeti, J.; Winkler, R.G.; Gompper, G. Physics of microswimmers single particle motion and collective behavior: A review. Rep. Prog. Phys. 2015, 78, 056601. [Google Scholar] [CrossRef] [Green Version]
  3. Sriram Ramaswamy. Active matter. J. Stat. Mech. 2017, 2017, 054002. [Google Scholar] [CrossRef]
  4. Hagan, M.F.; Grason, G.M. Equilibrium mechanisms of self-limiting assembly. Rev. Mod. Phys. 2021, 93, 025008. [Google Scholar] [CrossRef] [PubMed]
  5. Vasiliev, M.M.; Antipov, S.N.; Petrov, O.F. Large-scale vortices in dc glow discharge dusty plasmas. J. Phys. A Math. Gen. 2006, 39, 4539–4544. [Google Scholar] [CrossRef]
  6. Trukhachev, F.M.; Vasiliev, M.M.; Petrov, O.F.; Vasilieva, E.V. Microdynamic and thermodynamic properties of dissipative dust-acoustic solitons. J. Phys. A Math. Theor. 2021, 54, 095702. [Google Scholar] [CrossRef]
  7. Vasilieva, E.V.; Petrov, O.F.; Vasiliev, M.M. Laser-induced melting of two-dimensional dusty plasma system in RF discharge. Sci. Rep. 2021, 11, 523. [Google Scholar] [CrossRef]
  8. Vaulina, O.S.; Vasilieva, E.V. Orientational order and topological defects in two-dimensional Yukawa systems. Phys. Lett. A 2014, 378, 719–722. [Google Scholar] [CrossRef]
  9. Petrov, O.F.; Vasiliev, M.M.; Tun, Y.; Statsenko, K.B.; Vaulina, O.S.; Vasilieva, E.V.; Fortov, V.E. Two-dimensional phase transition in a strongly nonideal dusty plasma. J. Exp. Theor. Phys. 2015, 120, 327–332. [Google Scholar] [CrossRef]
  10. Totsuji, H. Structure and melting of two-dimensional dust crystals. Phys. Plasmas 2001, 8, 1856–1862. [Google Scholar] [CrossRef]
  11. Henning, C.; Baumgartner, H.; Piel, A.; Ludwig, P.; Golubnichiy, V.; Bonitz, M.; Block, D. Ground state of a confined Yukawa plasma. Phys. Rev. E 2006, 74, 056403. [Google Scholar] [CrossRef] [Green Version]
  12. Timofeev, A.V.; Nikolaev, V.S.; Semenov, V.P. Inhomogeneity of structural and dynamical characteristics of dusty plasma in a gas discharge. J. Exp. Theor. Phys. 2020, 130, 153–160. [Google Scholar] [CrossRef]
  13. Nikolaev, V.S.; Timofeev, A.V. Inhomogeneity of a harmonically confined Yukawa system. Phys. Plasmas 2019, 26, 073701. [Google Scholar] [CrossRef]
  14. Hariprasad, M.; Bandyopadhyay, P.; Agora, G.; Sen, A. Experimental observation of a dusty plasma crystal in the cathode sheath of a DC glow discharge plasma. Phys. Plasmas 2018, 25, 123704. [Google Scholar] [CrossRef] [Green Version]
  15. Hariprasad, M.G.; Bandyopadhyay, P.; Nikolaev, V.S.; Kolotinskii, D.A.; Arumugam, S.; Arora, G.; Singh, S.; Sen, A.; Timofeev, A.V. Self-sustained non-equilibrium co-existence of fluid and solid states in a strongly coupled complex plasma system. Sci. Rep. 2022, 12, 13882. [Google Scholar] [CrossRef]
  16. Kolotinskii, D.A.; Nikolaev, V.S.; Timofeev, A.V. Effect of Structural Inhomogeneity and Nonreciprocal Effects in the Interaction of Macroparticles on the Dynamic Properties of a Dusty Plasma Monolayer. JETP Lett. 2021, 113, 510–517. [Google Scholar] [CrossRef]
  17. Nikolaev, V.S.; Timofeev, A.V. Nonhomogeneity of phase state in a dusty plasma monolayer with nonreciprocal particle interactions. Phys. Plasmas 2021, 28, 033704. [Google Scholar] [CrossRef]
  18. Klumov, B.A. On the Effect of Confinement on the Structure of a Complex (Dusty) Plasma. JETP Lett. 2019, 110, 715–721. [Google Scholar] [CrossRef]
  19. Petrov, O.F.; Statsenko, K.B.; Vasiliev, M.M. Active Brownian motion of strongly coupled charged grains driven by laser radiation in plasma. Sci. Rep. 2022, 12, 8618. [Google Scholar] [CrossRef]
  20. Vrancken, R.; Paeva, G.V.; Kroesen, G.M.W.; Stoffels, W.W. Dust void formation above rectangular and circular potential traps in an rf plasma. Plasma Sources Sci. Technol. 2005, 14, 509. [Google Scholar] [CrossRef]
  21. Tsytovich, V.N.; Vladimirov, S.V.; Morfill, G.E. Size of dust voids as a function of the power input in dusty plasma. J. Exp. Theor. Phys. 2006, 102, 334. [Google Scholar] [CrossRef]
  22. Hu, Z.; Chen, Y.; Zheng, X.; Huang, F.; Shi, G.-F.; Yu, M.Y. Theory of void formation in dusty plasmas. Phys. Plasmas 2009, 16, 063707. [Google Scholar] [CrossRef]
  23. Wang, Y.; Zhu, X.; Chen, Z.; Huang, F. Numerical study on the influence of electron temperature on void formation in dusty plasma. Phys. Scr. 2020, 95, 055606. [Google Scholar] [CrossRef]
  24. Holstein, T. Energy distribution of electrons in high frequency gas discharges. Phys. Rev. 1946, 70, 367–384. [Google Scholar] [CrossRef]
  25. Huxley, L.G.H.; Crompton, R.W. The Diffusion and Drift of Electrons in Gases; Wile-Interscience Pub.; John and Wiley and Sons: Hoboken, NJ, USA, 1974. [Google Scholar]
  26. Raizer, Y.P.; Shneider, M.N.; Yatsenko, N.A. Radio-Frequency Capacitive Discharges; CRC Press: Boca Raton, FL, USA, 2017. [Google Scholar] [CrossRef]
  27. Barnes, M.S.; Keller, J.H.; Forster, J.C.; O’Neill, J.A.; Coultas, D.K. Transport of dust particles in glow-discharge plasmas. Phys. Rev. Lett. 1992, 68, 313–316. [Google Scholar] [CrossRef]
  28. Hagelaar, G.J.M. Coulomb collisions in the Boltzmann equation for electrons in low-temperature gas discharge plasmas. Plasma Sources Sci. Technol. 2016, 25, 015015. [Google Scholar] [CrossRef] [Green Version]
  29. Hagelaar, G.J.M.; Pitchford, L.C. Solving the Boltzmann equation to obtain electron transport coefficients and rate coefficients for fluid models. Plasma Sources Sci. Technol. 2005, 14, 722–733. [Google Scholar] [CrossRef]
  30. Pitchford, L.C.; Alves, L.L.; Bartschat, K.; Biagi, S.F.; Bordage, M.-C.; Bray, I.; Brion, C.E.; Brunger, M.J.; Campbell, L.; Chachereau, A.; et al. LXCat: An Open-Access, Web-Based Platform for Data Needed for Modeling Low Temperature Plasmas. Plasma Process. Polym. 2017, 14, 1600098. [Google Scholar] [CrossRef]
  31. Dashdorj, J.; Pfalzgraff, W.C.; Trout, A.; Fingerlow, D.; Cordier, M.; Viehland, L.A. Determination of mobility and diffusion coefficients of Ar+ and Ar2+ ions in argon gas. Int. J. Ion Mobil. Spec. 2020, 23, 143–151. [Google Scholar] [CrossRef]
  32. Basurto, E.; de Urquijo, J.; Alvarez, I.; Cisneros, C. Mobility of He+, Ne+, Ar+, N2+, O2+ and CO2+ in their parent gas. Phys. Rev. E 2000, 61, 3053–3057. [Google Scholar] [CrossRef]
  33. Helm, H.; Elford, M.T. The mobility of Ar+ ions in argon and the effect of spin-orbit coupling. J. Phys. B At. Mol. Phys. 1977, 10, 3849. [Google Scholar] [CrossRef]
  34. Hegerberg, R.; Elford, M.T.; Skullerud, H.R. The cross section for symmetric charge exchange of Ne+ in Ne and Ar+ in Ar at low energies. J. Phys. B At. Mol. Phys. 1982, 15, 797–811. [Google Scholar] [CrossRef]
  35. Chicheportiche, A.; Lepetit, B.; Gadéa, F.X.; Benhenni, M.; Yousfi, M.; Kalus, R. Ab initio transport coefficients of Ar+ ions in Ar for cold plasma jet modeling. Phys. Rev. E 2014, 89, 063102. [Google Scholar] [CrossRef]
  36. Radovanov, S.B.; Van Brunt, R.J.; Olthoff, J.K.; Jelenkovic, B.M. Ion kinetics and symmetric charge-transfer collisions in low-current, diffuse (Townsend) discharges in argon and nitrogen. Phys. Rev. E 1995, 51, 4036–4046. [Google Scholar] [CrossRef] [Green Version]
  37. Rao, M.V.V.S.; Van Brunt, R.J.; Olthoff, J.K. Resonant charge exchange and the transport of ions at high electric-field to gas-density ratios E/N in argon, neon, and helium. Phys. Rev. E 1996, 54, 5642–5657. [Google Scholar] [CrossRef] [Green Version]
  38. Ellis, H.W.; Pai, R.Y.; McDaniel, E.W.; Mason, E.A.; Viehland, L.A. Transport properties of gaseous ions over a wide energy range. At. Data Nucl. Data Tables 1976, 17, 177. [Google Scholar] [CrossRef]
  39. Wannier, G.H. On the Motion of Gaseous Ions in a Strong Electric Field. I. Phys. Rev. 1951, 83, 281. [Google Scholar] [CrossRef]
  40. Filippov, A.V.; Babichev, V.N.; Dyatko, N.A.; Pal’, A.F.; Starostin, A.N.; Taran, M.D.; Fortov, V.E. Formation of plasma dust structures at atmospheric pressure. J. Exp. Theor. Phys. 2006, 102, 342–354. [Google Scholar] [CrossRef]
  41. Filippov, A.V.; Babichev, V.N.; Pal’, A.F.; Starostin, A.N.; Cherkovets, V.E.; Rerikh, V.K.; Taran, M.D. Dust Trap Formation in a Non-Self-Sustained Discharge with External Gas Ionization. Plasma Phys. Rep. 2015, 41, 895–904. [Google Scholar] [CrossRef]
  42. Filippov, A.V.; Babichev, V.N.; Pal’, A.F.; Starostin, A.N.; Cherkovets, V.E. A Dusty Plasma in a Non-Self-Sustained Gas Discharge at Atmospheric Pressure. Contr. Plasma Phys. 2016, 56, 286–295. [Google Scholar] [CrossRef]
  43. Epstein, P.S. On the Resistance Experienced by Spheres in their Motion. Phys. Rev. 1924, 23, 710–733. [Google Scholar] [CrossRef] [Green Version]
  44. Vladirnirov, S.V.; Ostrikov, K.; Sarnarian, A.A. Physics and Applications of Complex Plasmas; Imperial College Press: London, UK, 2005; p. 36. [Google Scholar]
  45. Morfill, G.E.; Ivlev, A.V. Complex plasmas: An interdisciplinary research field. Rev. Mod. Phys. 2009, 81, 1353–1404. [Google Scholar] [CrossRef]
  46. Lieberman, M.A. Analytical solution for capacitive RF sheath. IEEE Trans. Plasma Sci. 1988, 16, 638–644. [Google Scholar] [CrossRef]
  47. Lieberman, M.A. Dynamics of a collisional, capacitive RF sheath. IEEE Trans. Plasma Sci. 1989, 17, 338–341. [Google Scholar] [CrossRef]
  48. Godyak, V.A.; Sternberg, N. Dynamic model of the electrode sheaths in symmetrically driven rf discharges. Phys. Rev. A 1990, 42, 2299. [Google Scholar] [CrossRef]
  49. Czarnetzki, U. Analytical model for the radio-frequency sheath. Phys. Rev. E 2013, 88, 063101. [Google Scholar] [CrossRef]
  50. Godyak, V.A.; Piejak, R.B.; Alexandrovich, B.M. Electrical characteristics of parallel-plate RF discharges in argon. IEEE Trans. Plasma Sci. 1991, 19, 660–676. [Google Scholar] [CrossRef]
  51. Turner, M.M. Collisionless heating in radio-frequency discharges: A review. J. Phys. D Appl. Phys. 2009, 42, 194008. [Google Scholar] [CrossRef]
  52. Lafleur, T.; Chabert, P. Is collisionless heating in capacitively coupled plasmas really collisionless? Plasma Sources Sci. Technol. 2015, 24, 044002. [Google Scholar] [CrossRef]
  53. Schulze, J.; Donkó, Z.; Lafleur, T.; Wilczek, S.; Brinkmann, R.P. Spatio-temporal analysis of the electron power absorption in electropositive capacitive RF plasmas based on moments of the Boltzmann equation. Plasma Sources Sci. Technol. 2018, 27, 055010. [Google Scholar] [CrossRef] [Green Version]
  54. Vass, M.; Wilczek, S.; Lafleur, T.; Brinkmann, R.P.; Donkó, Z.; Schulze, J. Observation of dominant Ohmic electron power absorption in capacitively coupled radio frequency argon discharges at low pressure. Plasma Sources Sci. Technol. 2020, 29, 085014. [Google Scholar] [CrossRef]
  55. Wilczek, S.; Schulze, J.; Brinkmann, R.P.; Donkó, Z.; Trieschmann, J.; Mussenbrock, T. Electron dynamics in low pressure capacitively coupled radio frequency discharges. J. Appl. Phys. 2020, 127, 181101. [Google Scholar] [CrossRef]
  56. Kaganovich, I.D.; Polmarov, O.V.; Theodosiou, C.E. Revisiting the anomalous RF field penetration into a warm plasma. IEEE Trans. Plasma Sci. 2006, 34, 696–717. [Google Scholar] [CrossRef] [Green Version]
  57. Vaulina, O.S.; Khrapak, S.A. Scaling Law for the Fluid-Solid Phase Transition in Yukawa Systems (Dusty Plasmas). J. Exp. Theor. Phys. 2000, 90, 287–289. [Google Scholar] [CrossRef]
  58. Filippov, A.V.; Reshetnyak, V.V.; Starostin, A.N.; Tkachenko, I.M.; Fortov, V.E. Investigation of Dusty Plasma Based on the Ornstein—Zernike Integral Equation for a Multicomponent Fluid. JETP Lett. 2019, 110, 659–666. [Google Scholar] [CrossRef]
  59. Filippov, A.V.; Fortov, V.E.; Reshetniak, V.V.; Starostin, A.N.; Tkachenko, I.M. Electrostatic interactions and stability of dusty plasmas and the multicomponent Ornstein–Zernike equation. AIP Adv. 2020, 10, 045232. [Google Scholar] [CrossRef]
Figure 1. Video frame of isotropic dust monolayer structure in plasma of RF discharge at W = 6.9 W, P = 5 Pa.
Figure 1. Video frame of isotropic dust monolayer structure in plasma of RF discharge at W = 6.9 W, P = 5 Pa.
Molecules 28 03259 g001
Figure 2. Quasi-two-dimensional dust structures in argon plasma of RF discharge at P = 5 Pa and various values of discharge power W: (a) 4.4 W, (b) 6.9 W, (c) 16.7 W. Color bar shows radial distribution of interparticle distances (upper row), Voronoi diagrams (lower row) represent different degrees of structure inhomogeneity.
Figure 2. Quasi-two-dimensional dust structures in argon plasma of RF discharge at P = 5 Pa and various values of discharge power W: (a) 4.4 W, (b) 6.9 W, (c) 16.7 W. Color bar shows radial distribution of interparticle distances (upper row), Voronoi diagrams (lower row) represent different degrees of structure inhomogeneity.
Molecules 28 03259 g002
Figure 3. A dependence of the radial distribution of the average interparticle distance in the quasi-two-dimensional structure in an argon plasma of RF discharge at P = 5 Pa and various values of discharge power W: blue curve is for W = 4.4 W; red is for 6.9 W; green is for 16.7 W.
Figure 3. A dependence of the radial distribution of the average interparticle distance in the quasi-two-dimensional structure in an argon plasma of RF discharge at P = 5 Pa and various values of discharge power W: blue curve is for W = 4.4 W; red is for 6.9 W; green is for 16.7 W.
Molecules 28 03259 g003
Figure 4. Quasi-two-dimensional dust structures in argon plasma of RF discharge at W = 13.1 W and various values of argon pressure P: (a) 3 Pa, (b) 4 Pa, (c) 5 Pa. Color bar shows radial distribution of interparticle distances (upper row), Voronoi diagrams (lower row) represent different degrees of structure inhomogeneity.
Figure 4. Quasi-two-dimensional dust structures in argon plasma of RF discharge at W = 13.1 W and various values of argon pressure P: (a) 3 Pa, (b) 4 Pa, (c) 5 Pa. Color bar shows radial distribution of interparticle distances (upper row), Voronoi diagrams (lower row) represent different degrees of structure inhomogeneity.
Molecules 28 03259 g004
Figure 5. A dependence of the radial distribution of the average interparticle distance in the quasi-two-dimensional structure in an argon plasma of RF discharge at W = 13.1 W and various values of argon pressure P: (a) blue curve is for 3 Pa, (b) red is for 4 Pa, (c) green is for 5 Pa.
Figure 5. A dependence of the radial distribution of the average interparticle distance in the quasi-two-dimensional structure in an argon plasma of RF discharge at W = 13.1 W and various values of argon pressure P: (a) blue curve is for 3 Pa, (b) red is for 4 Pa, (c) green is for 5 Pa.
Molecules 28 03259 g005
Figure 6. The dependence of the average electron energy on the specific intensity of the RF electric field with a frequency of 13.56 MHz in argon at pressure 1—3 Pa, 2—3.6 Pa, 3—4.5 Pa, 4—6 Pa.
Figure 6. The dependence of the average electron energy on the specific intensity of the RF electric field with a frequency of 13.56 MHz in argon at pressure 1—3 Pa, 2—3.6 Pa, 3—4.5 Pa, 4—6 Pa.
Molecules 28 03259 g006
Figure 7. Dependences of the ratio of the averaged electron transport frequency to the circular frequency of the RF field (a), the mobility (b) and the electron diffusion coefficient (c), the ionization rate constant of atoms (d) on the average electron energy in the RF electric field with a frequency of 13.56 MHz in argon at pressure 1—3 Pa, 2—3.6 Pa, 3—4.5 Pa, 4—6 Pa.
Figure 7. Dependences of the ratio of the averaged electron transport frequency to the circular frequency of the RF field (a), the mobility (b) and the electron diffusion coefficient (c), the ionization rate constant of atoms (d) on the average electron energy in the RF electric field with a frequency of 13.56 MHz in argon at pressure 1—3 Pa, 2—3.6 Pa, 3—4.5 Pa, 4—6 Pa.
Molecules 28 03259 g007
Figure 8. The mobility of Ar+ ions reduced to normal conditions in argon at a temperature of 300 K as a function of the reduced field E/N: 1—experimental data [29], 2—approximation from [30], 3—work [33], 4—[34], 5—approximation Equation (10).
Figure 8. The mobility of Ar+ ions reduced to normal conditions in argon at a temperature of 300 K as a function of the reduced field E/N: 1—experimental data [29], 2—approximation from [30], 3—work [33], 4—[34], 5—approximation Equation (10).
Molecules 28 03259 g008
Figure 9. A dependence of parameter ζ on a reduced field E/N in RF discharge with a frequency of 13.56 MHz and inter-electrode distance L = 5 cm at various argon pressures: curve 1—3 Pa, curve 2—4 Pa, curve 3—5 Pa, curve 4—6 Pa.
Figure 9. A dependence of parameter ζ on a reduced field E/N in RF discharge with a frequency of 13.56 MHz and inter-electrode distance L = 5 cm at various argon pressures: curve 1—3 Pa, curve 2—4 Pa, curve 3—5 Pa, curve 4—6 Pa.
Molecules 28 03259 g009
Figure 10. Pair correlation functions g(r), (a,b) determined by numeric solution of the Ornstein–Zernike equation in HNC approximation under various values of effective coupling parameters Γ* and nd = 8.0 × 103 cm−3, and (c,d) obtained in experiments under the following parameters: (a) η = 5.53, RD = 1/kD = 293.3 μm, λ = 1.7, zd = −1687.7; (b) η = 0 (Wi = 0), RD = 113.4 μm, λ = 4.4, zd = −3420.6; (c) Discharge power W = 13.1 W, pressure P = 3–5 Pa; (d) Discharge power W = 4.4–16.7 W, pressure P = 5 Pa.
Figure 10. Pair correlation functions g(r), (a,b) determined by numeric solution of the Ornstein–Zernike equation in HNC approximation under various values of effective coupling parameters Γ* and nd = 8.0 × 103 cm−3, and (c,d) obtained in experiments under the following parameters: (a) η = 5.53, RD = 1/kD = 293.3 μm, λ = 1.7, zd = −1687.7; (b) η = 0 (Wi = 0), RD = 113.4 μm, λ = 4.4, zd = −3420.6; (c) Discharge power W = 13.1 W, pressure P = 3–5 Pa; (d) Discharge power W = 4.4–16.7 W, pressure P = 5 Pa.
Molecules 28 03259 g010
Figure 11. A scheme of an experimental setup.
Figure 11. A scheme of an experimental setup.
Molecules 28 03259 g011
Table 1. Parameters of the RF discharge in the first series of experiments at the fixed discharge power (Win = 15 W, Wout = 4.7 W) and at different pressures.
Table 1. Parameters of the RF discharge in the first series of experiments at the fixed discharge power (Win = 15 W, Wout = 4.7 W) and at different pressures.
P (Pa)E/N (Td)ne (cm−3)Te (eV)RDi (μm)RDe (μm)le (cm)li (cm)
3142.463.39 × 1084.457379.4851.92.3040.078
4120.452.95 × 1084.472354.1915.01.6630.061
5107.892.61 × 1084.478336.7974.11.2970.050
699.662.34 × 1084.476324.11028.41.0630.042
Table 2. Parameters of the RF discharge in the second series of experiments at the fixed pressure P = 5 Pa (E/N = 107.9 Td, Te = 4.478 eV, le = 1.297 cm, li = 0.05 cm).
Table 2. Parameters of the RF discharge in the second series of experiments at the fixed pressure P = 5 Pa (E/N = 107.9 Td, Te = 4.478 eV, le = 1.297 cm, li = 0.05 cm).
Win (W)Wout (W)we (eV cm−3 s−1)ne (cm−3)RDi (μm)RDe (μm)η = Wi/We
304.604.03 × 10154.21 × 108264.5767.010.14
203.303.24 × 10153.37 × 108295.7856.48.13
81.102.01 × 10152.10 × 108375.51085.55.06
50.601.57 × 10151.64 × 108425.41228.63.95
Table 3. Parameters of the RF discharge and dust particles (with radius r0 = 5.3 μm and density ρ = 1.5 g/cm3) in the levitation area during the first series of experiments at the fixed discharge power.
Table 3. Parameters of the RF discharge and dust particles (with radius r0 = 5.3 μm and density ρ = 1.5 g/cm3) in the levitation area during the first series of experiments at the fixed discharge power.
P (Pa)ϕ0 (V), Ti = Tgaszd0 (e) Ti = Tgasϕd (V)zd (e)Elev (V/cm)E (V/cm)El,ion (eV)
3−9.906−3.94 × 104−15.53−5.81 × 1049.871.0681.326
4−9.934−3.94 × 104−15.14−5.66 × 10410.121.2051.005
5−9.945−3.93 × 104−14.80−5.54 × 10410.351.3490.803
6−9.941−3.91 × 104−14.49−5.43 × 10410.561.4950.667
Table 4. Parameters of the RF discharge and dust particles (with radius r0 = 5.3 μm and density ρ = 1.5 g/cm3) in the levitation area during the second series of experiments at the fixed pressure.
Table 4. Parameters of the RF discharge and dust particles (with radius r0 = 5.3 μm and density ρ = 1.5 g/cm3) in the levitation area during the second series of experiments at the fixed pressure.
we (eV/cm3 s)ϕ0 (V) Ti = Tgaszd0 (e) Ti = Tgasϕd (V)zd (e)Elev (V/cm)E (V/cm)El,ion (eV)
4.03 × 1015−9.945−3.99 × 104−14.79−5.56 × 10410.311.3490.799
3.24 × 1015−9.945−3.96 × 104−14.80−5.55 × 10410.331.3490.801
2.01 × 1015−9.945−3.90 × 104−14.80−5.53 × 10410.361.3490.804
1.57 × 1015−9.945−3.87 × 104−14.80−5.52 × 10410.381.3490.805
Table 5. Ion drag force, power of stochastic heating of electrons, transport frequency (νm), and ratio of ion and electron plasma frequencies to the circular frequency of the RF discharge ω = 2π × 13.56 MHz = 8.52 × 107 s−1.
Table 5. Ion drag force, power of stochastic heating of electrons, transport frequency (νm), and ratio of ion and electron plasma frequencies to the circular frequency of the RF discharge ω = 2π × 13.56 MHz = 8.52 × 107 s−1.
P (Pa)Win (W)we
(eV cm−3 s−1)
Wst
(eV cm−3 s−1)
Fid (dyn) *νm (s−1) ( ω p i ω ) 2 ( ω p e ω ) 2
3152.79 × 10151.23 × 10168.36 × 10−73.46 × 1072.04 × 10−3149
4152.62 × 10157.77 × 10157.51 × 10−72.80 × 1071.78 × 10−3129
5152.50 × 10155.32 × 10156.60 × 10−72.34 × 1071.57 × 10−3114
6152.41 × 10153.86 × 10155.74 × 10−72.00 × 1071.41 × 10−3103
5304.03 × 10151.02 × 10158.57 × 10−72.34 × 1072.53 × 10−3184
5203.24 × 10157.61 × 10157.64 × 10−72.34 × 1072.03 × 10−3148
582.01 × 10153.89 × 10155.80 × 10−72.34 × 1071.26 × 10−392.0
551.57 × 10152.68 × 10154.96 × 10−72.34 × 1079.86 × 10−371.9
* mdg = ezd Elev = 9.17 × 10−7 dyn.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Alekseevskaya, A.A.; Vasilieva, E.V.; Filippov, A.V.; Vasiliev, M.M.; Petrov, O.F. Isotropic and Anisotropic Monolayer Structures in RF Discharge Plasma. Molecules 2023, 28, 3259. https://doi.org/10.3390/molecules28073259

AMA Style

Alekseevskaya AA, Vasilieva EV, Filippov AV, Vasiliev MM, Petrov OF. Isotropic and Anisotropic Monolayer Structures in RF Discharge Plasma. Molecules. 2023; 28(7):3259. https://doi.org/10.3390/molecules28073259

Chicago/Turabian Style

Alekseevskaya, Anastasiya A., Elena V. Vasilieva, Anatoly V. Filippov, Mikhail M. Vasiliev, and Oleg F. Petrov. 2023. "Isotropic and Anisotropic Monolayer Structures in RF Discharge Plasma" Molecules 28, no. 7: 3259. https://doi.org/10.3390/molecules28073259

APA Style

Alekseevskaya, A. A., Vasilieva, E. V., Filippov, A. V., Vasiliev, M. M., & Petrov, O. F. (2023). Isotropic and Anisotropic Monolayer Structures in RF Discharge Plasma. Molecules, 28(7), 3259. https://doi.org/10.3390/molecules28073259

Article Metrics

Back to TopTop