Conductivity,effects,during,the,transition,from,collisionless,to,collisional,regimes,in,cylindrical,inductively,coupled,plasmas
Wei YANG (杨唯), Fei GAO (高飞) and Younian WANG (王友年)
1 College of Science, Donghua University, Shanghai 201620, People’s Republic of China
2 School of Physics, Dalian University of Technology, Dalian 116024, People’s Republic of China
Abstract A numerical model is developed to study the conductivity effects during the transition from collisionless to collisional regimes in cylindrical inductively coupled argon plasmas at pressures of 0.1-20 Pa.The model consists of electron kinetics module, electromagnetics module, and global model module.It allows for self-consistent description of non-local electron kinetics and collisionless electron heating in terms of the conductivity of homogeneous hot plasma.Simulation results for non-local conductivity case are compared with predictions for the assumption of local conductivity case.Electron densities and effective electron temperatures under non-local and local conductivities show obvious differences at relatively low pressures.As increasing pressure, the results under the two cases of conductivities tend to converge, which indicates the transition from collisionless to collisional regimes.At relatively low pressures the local negative power absorption is predicted by non-local conductivity case but not captured by local conductivity case.The two-dimensional (2D) profiles of electron current density and electric field are coincident for local conductivity case in the pressure range of interest, but it roughly holds true for non-local conductivity case at very high pressure.In addition,an effective conductivity with consideration of non-collisional stochastic heating effect is introduced.The effective conductivity almost reproduces the electron density and effective electron temperature for the non-local conductivity case,but does not capture the non-local relation between electron current and electric field as well as the local negative power absorption that is observed for nonlocal conductivity case at low pressures.
Keywords: inductively coupled plasmas, conductivity effects, electron kinetics, plasma parameters, electromagnetic wave characteristics, electrodynamics
Low pressure inductively coupled plasmas (ICPs) are widely utilized for ion sources and material processing.Their popularity for both scientific and industrial applications depends on the capability to achieve high density even for relatively low power and high plasma uniformity, as well as on the possibility of independent control of the ion energy at the plasma boundary.
In the past few decades, substantial experimental and theoretical efforts on the electron kinetics and the electron heating mechanisms have been made in order to understand low pressure ICPs and optimize the discharge characteristics in terms of the individual requirements [1, 2].The common feature of low pressure ICPs is their non-equilibrium nature.In particular,electrons are not only far from equilibrium with ions and neutrals, but even not in equilibrium within their own ensemble [3].Non-local electrodynamics and non-local electron kinetics as well as anomalous skin effect and local negative power absorption are distinctive features of the lowpressure regime [4-8].The differences between non-local electrodynamics and non-local electron kinetics have been introduced in detail by Kolobov et al[8].The electron energy probability function (EEPF) can deviate from a Maxwellian distribution in these discharges (non-local electron kinetics)[9, 10].Due to the smaller collision frequency of electrons with neutrals than excitation frequency, the local relation between electron current and electric field (Ohm’s law) does not hold true (non-local electrodynamics).In this case, the electron current relies on the distribution of the electric field in the whole plasma and size of the discharge device[11-13].It is due to the spatial dispersion of conductivity caused by electron thermal motion.Therefore, a non-local conductivity is required to determine the electric field penetration into the plasma.
The accurate description of the discharge properties under low pressures relies on kinetic simulations, where the EEPF,non-local conductivity and plasma density need to be self-consistently calculated.The non-local conductivities for homogeneous plasmas were derived in planar-type [14, 15] and solenoidal-type [16] ICP reactors by Yoon et al.Due to assuming a constant electron density and a Maxwellian EEPF,only the averaged electron Boltzmann equation under linear approximation and the 2D Maxwell equations were solved in their models.Therefore, the interaction between electron properties and electric field cannot be self-consistently considered.Regarding the complexity of theoretical studies in cylindrical geometry,a self-consistent kinetic model with the calculations of nonlocal conductivity, EEPF, and non-uniform profiles of electron density in a one-dimensional (1D) slab geometry was developed by Kaganovich et al [12].Their model contains the averaged electron Boltzmann equation for EEPF, the Maxwell equation for the electric field, the quasi-neutrality condition for electrostatic potential, and the fluid conservation equations for ion density profile.Subsequently,the kinetic model was used to study the influences of non-local conductivity and non-Maxwellian EEPF on plasma power absorption and plasma density profiles in a 1D planar inductively coupled argon discharge under the non-local regime (pressure ≤10 mTorr) [17, 18].Recently, Kushner et al reported the influence of non-local conductivity on plasma absorption power in a 2D planar-type argon ICP at 10 mTorr,where electron Monte Carlo simulation instead of electron Boltzmann equation was used to calculate the electron current density [19].
The aforementioned theoretical modellings in the collisionless regime have been well established, while the consideration of inhomogeneous plasmas in a 2D cylindrical geometry is very cumbersome.The fluid equations need to be solved and meanwhile the EEPF as a function of not only the kinetic energy but the potential energy needs to be taken into account.It might limit the studies to some fixed discharge conditions regarding the computational expense.To the best of our knowledge,the effects of plasma conductivity have not been fully reported over a wide pressure range that involves the transition from collisionless to collisional regimes.As mentioned earlier, plasma conductivity couples plasma kinetics to the electromagnetic fields.In this study,a numerical model based on previous work[20]that has been validated is developed to study the conductivity effects during the transition from collisionless to collisional regimes in cylindrical inductively coupled Ar plasmas.In the model, the electron kinetics are addressed by the electron Boltzmann equation, electromagnetic wave properties are obtained by the Maxwell’s equations,and plasma density is calculated from global model.The non-local conductivity,local conductivity, and effective conductivity are respectively considered in the model and their effects on plasma parameters and electromagnetic wave characteristics are compared under a wide pressure range of 0.1-20 Pa.The non-local conductivity is computed by the electron Boltzmann equation, while the local conductivity is obtained based on the Ohm’s law.The effective conductivity is introduced based on the Ohm’s law but with consideration of non-collisional stochastic heating effect by defining a stochastic collision frequency.In principle, the nonlocal conductivity model should be valid no matter which heating mechanisms (i.e.collisional, collisionless or both of them) dominate the discharge in a given situation.
This manuscript is structured as follows.Section 2 shows a brief description of the model.Verification and validation of the model are shown in section 3.Simulation results and discussions are presented in section 4, where the non-local conductivity effects are compared with local conductivity case and with effective conductivity case.Section 5 presents the conclusions.
In this work, we consider a cylindrical ICP source with a solenoidal-type coil, and the structure of the chamber is the same as that of the driver chamber of a two-chamber radiofrequency (RF) ICP in the experiments by Li et al [21].The radius R is 6 cm and the length H is 14 cm.The model is based on the assumption of homogeneous plasma, since the discharge is assumed to be under inductive H mode to produce a thin sheath near the walls.
The model mainly includes three modules of the electron kinetics module (EKM), electromagnetics module (EMM),and global model module(GMM).Harmonic electromagnetic field is generated by the alternating current in the coil by solving the frequency domain wave equation (i.e.EMM),where all physical quantities are assumed to be axisymmetric.For local model,electron current density is obtained based on Ohm’s law where the conductivity is local.For non-local model, the determination of electron current density needs a non-local conductivity that is deduced from the EKM (linearized electron Boltzmann equation for perturbed distribution function) with EMM.Electric field from the EMM is used in the EKM (electron Boltzmann equation under quasilinear approximation for isotropic distribution function) to produce normalized EEPF.The normalized EEPF is used to calculate electron impact rate coefficients and other related parameters for the use in the GMM.Particle balance equations for densities of particle species, and power balance equation for density of electrons are solved in the GMM.It is worth noting that only the parts involved in this work as well as variations and improvements in terms of our previous model are shown here and more details can be found in our previous work[20].
2.1.Electromagnetics module
The solution of the inductive electric field can be obtained by the mode analysis method using the Fourier series in a rectangular coordinate system(x,z) [20]
whereknandqmcan refer to the previous work [20], and the expansion coefficientEnmof the electric field is given as[20]
whereμ0,ω,I0,N, andzkare respectively the vacuum permeability,excitation frequency,current amplitude in the coil,number of turns of the coil,and location of thekth turn in the axial direction.k0iswith speed of light c.The conductivityσwill be discussed in the EKM.
Based on the solution of electric field, the electron current densityJpcan be expressed in terms ofEyvia a conductivityσ, and it is given as [20]
2.2.Electron kinetics module
In this part, the non-local conductivity, local conductivity,and effective conductivity, as well as normalized EEPF will be presented.
2.2.1.Non-local conductivity.In the non-local limit, where the electron-neutral collision frequencyνenis of the same order as (or even smaller than) the excitation frequencyω,collisionless heating mechanism might compete with or even dominate over collisional heating in the low-pressure regime.The Ohm’s law for electron current densityJpis therefore considered to be not valid anymore.To account for collisionless electron heating,Jpis calculated based on the perturbed distribution functionf1,and it is given by
whereeis the elementary charge,neis the electron density,andvyis the velocity component in the y direction.
The non-local conductivity can be self-consistently computed by equations (3) and (4).σnis thenth coefficient of non-local conductivity in the Fourier series,and it is given as [20]
whereεis the electron kinetic energy,meis the electron mass,is the electron velocity, andf0is isotropic distribution function(i.e.the normalized EEPF).Ψ andΦ can refer to the previous work [20].
2.2.2.Local conductivity.When the electron-neutral collision frequencyνenis much larger than the excitation frequencyωat higher pressures, the conductivity is defined by the Ohm’s law.It indicates that the collisional heating is the only manner that the electromagnetic field transfers random thermal energy to the plasma.In cases where the local relation between electron current and electric field is valid,the local conductivity for an arbitrary EEPF is given by[17, 22, 23]
The expansion coefficientEnm(see equation (2)) of electric fieldEyand electron current densityJp(see equation(3))mentioned in the EMM are therefore determined based on the three cases of conductivities.
2.2.4.Normalized EEPF.The normalized EEPFf0can be solved by
where the energy diffusion coefficientDeis deduced in appendix in detail and is expressed as
Table 1.Reactions considered in the model.
where the first term on the RHS has the same form as the term on the LHS of equation (8).DeeandVeein equation (10) are respectively given as [25]
whereveeis the frequency of electron-electron collision [25].
The normalization condition foris required.Equation (8) forf0can be written in the form of an ordinary nonlinear differential equation[25],and solved using the tridiagonal speed-up method.
2.3.Global model module
Global model is also called zero-dimensional model, since it ignores the spatial variation of plasma properties.As known,the global model usually includes particle balance equations for each species except for the background gas, the quasineutrality condition, the equation of state for total particle number density, and the power balance equation.Therefore,the global model can be used for the calculation of electron temperature and densities of all particle species.In this study,the normalized EEPFf0(effective electron temperatureTeff)is calculated using the electron Boltzmann equation under quasilinear approximation, so the global model is only used for the calculation of the densities of particle species.We discard one of the particle balance equations due to a reduced unknown,i.e.electron temperature in the global model.In our previous work [20], only the power balance equation that is used to determine the electron density is considered.In this work, the particle balance is also included to calculate densities of other particle species, because the Ar metastable is also considered.Even if including Ar metastable has very small effects on the plasma parameters and electromagnetic field properties in this work, the consideration of particle balance equation will help to extend the model to more complex gas discharges for future work.The particle balance equation other than background gas is given as [1, 26]
whereniis the density of the ith reactant.For the electronimpact volume reactions,kis calculated based on the normalized EEPF [1]
For the surface reaction for neutrals and positive ions,the rate coefficients can be found in [26].In this work, the reaction processes of Ar gas are shown in table 1.The cross sections of electron-impact elastic scattering, excitation, de-excitation and ionization are in [27-29].
The density of background gas is determined by the equation of state.Meanwhile, the discharge is assumed to satisfy the quasi-neutrality condition.These two criteria are respectively given by
whereTgasis the gas temperature andkBis the Boltzmann constant.In this work, equation (17) is simplified asne=
The power balance equation in our previous work[20]is based on the assumption of particle balance by equating the total surface particle loss to the total volume ionization[1].In fact, it is not quite appropriate, since the EEPF is not strictly Maxwellian and the electron temperature is not determined by the particle balance equation but by the electron Boltzmann equation under quasilinear approximation.It leads to the underestimation of electron density predicted by simulations[20]as increasing pressure when comparing with experiments[21], due to the overestimation of electron loss rate at relatively high pressures.In this work, the power balance equation is utilized in a more general form.The absorption power equals to the bulk power losses caused by electronneutral collisions, and surface power losses due to charged species flowing to the walls [26]:
where the absorption powerPabsis obtained from the EMM.In this work, we only have one kind of ions, i.e.Ar+, so the ion densityniequals to electron densityne.PVis the bulk power loss of an electron per unit volume andPSis the surface power loss of an ion on the walls per unit volume, respectively.They are given as [26]
3.1.Verification of the model
In order to verify the model,it is compared with the non-local approach to the solution of the Boltzmann equation developed by Kolobov et al [34].In their work, electron kinetics in a cylindrical ICP sustained by a long solenoidal coil was studied for the near-collisionless regime when the electron mean free path was large compared to the radius of the chamber.For simplicity, they prescribed 1D profiles of the fields and skin depth rather than self-consistently solving the Maxwell equations with electron kinetics and the plasma density.The amplitude of magnetic field (B0= 2 G), skin depth(δ= 1 cm) related to plasma density, radius of the chamber(R = 5 cm), and RF frequency (f= 13.56 MHz) are the input parameters in their model.The former two quantities are self-consistently calculated in our model.In order to present more reasonable comparison, we try to find the discharge condition capable of making the results of field and plasma density closer to the input ones used in the model of Kolobov et al and the electron mean free path larger than the radius of the chamber [34].In our simulation, the same radius and RF frequency are used.The length of the chamber is 15 cm, the pressure is 0.5 Pa and the absorption power is 80 W.Under the discharge condition,the calculated skin depthδis around 0.95 cm and the amplitude of electric fieldE0is 4 V cm−1that corresponds to the amplitude of magnetic fieldB0of around 1.9 G in our model.
Figure 1 shows that the normalized EEPF predicted by our model is compared with that by the model of Kolobov et al[34].In their model,the electron-electron collisions were excluded.The normalized EEPFs with and without consideration of electron-electron collisions are predicted by our model.Note that two models achieve reasonable agreement with each other.The model of Kolobov et al predicts a more remarkable depletion at the electron energy higher than 18 eV and more cold electrons with the energy lower than 2.5 eV.The differences between the two models can be due to the different collision processes considered in the two models.The collision processes included in the model of Kolobov et al were not mentioned [34].If more inelastic collision processes were included in the model of Kolobov et al, it would lead to more depletion in the energetic electrons and enhancement in low-energy ones.As expected, our model without consideration of electron-electron collisions agrees better with the model of Kolobov et al.The inclusion of the electron-electron collisions decreases cold electrons with the energy lower than 2.5 eV, and slightly increases energetic electrons with the energy higher than 18 eV.In general,electron-electron collisions tend to Maxwellize the EEPF.
Figure 1.Normalized EEPF versus electron energy.Predictions by the model in this work with and without electron-electron collisions and by the model of [34] without electron-electron collisions.Reprinted with permission from [34].Copyright (1997) by the American Physical Society.
In addition,the effects of chamber lengths(i.e.5,10,15,and 20 cm)have also been studied,even if they are not shown here.The simulated EEPF shows that the energetic electrons in inelastic energy range slightly increase with increasing chamber length, because the increased electron energy relaxation length decreases electron energy loss via the inelastic collisions.It leads to slightly higher effective electron temperature at larger chamber length.Shorter chamber length, e.g.5 cm, makes its scale smaller than the energy relaxation length, but the increased electron-electron collisions due to the increased plasma density still tend to Maxwellize the EEPF.This is because in the model the smaller chamber length produces larger power deposition density and thus higher plasma density for fixed total absorption power.
3.2.Validation of the model
In the previous study [20], the EEPF (electron kinetics)as well as the electron temperature and electron density (basic plasma properties) predicted by the model, was validated with the experiments[21]under the same discharge conditions.Here,we validate the model against experimental measurements in a cylindrical ICP discharge, where the data on electromagnetic wave characteristics under different discharge conditions are available[7].The experiment is composed of a Pyrex glass tube(Φ 26.3 cm (ID)×45 cm×0.5 cm) and a lower stainless-steel chamber (Φ 45 cm×25 cm×0.5 cm).The ICP source was sustained by a one-turn RF coil made of a copper strip(Φ 32 cm×6 cm×0.05 cm),and the vertical distance between the lower side of the copper strip and the upper side of the stainless-steel chamber was 14.9 cm.The plasma source was shielded with an aluminium cage to suppress RF electromagnetic interference.It is worth noting that the plasma can be very non-uniform especially along the axial direction for this kind of coil position and chamber size.The radial distribution of the RF magnetic field(axial component)was measured using the RF magnetic probe at the center of strip coil along the axial direction.The azimuthal electric field was obtained based on the measured magnetic field.
Figures 2(a) and (b) show the 1D radial distributions of electric fields versus power obtained from the experiments and simulations, respectively.Both experiments and simulations show that the electric field rapidly decays from the boundary near the coil to the center of chamber.The electric field decreases faster in the experiment than that in the simulation.This is because the measurements have been performed at the center of strip coil along the axial direction where the plasma density is almost the maximum.Nevertheless, the uniform plasma density assumed in the simulations can possibly make the volume-averaged density smaller than the one in the experiments.Therefore, the lower plasma density in the simulations leads to the deeper penetration and slower decay of the electric field into the plasma.As expected,the lower power facilitates the penetration of electric field into the plasma, which is captured by both simulations and experiments.Note that the higher power produces the nonmonotonic field decay,which is typically the anomalous skin effect caused by the collisionless electron thermal motion.The phenomenon has also been observed and discussed in the experiments of an ICP source driven by a planar coil [35].The non-monotonic evolution of the electric field appears at higher power in the simulation, due to higher power required in the simulation to produce the comparable plasma density to that in the experiment.
Figure 2.Radial distribution of normalized electric field under different absorption powers at 0.1 Pa.Experimental data (a)reprinted with the permission from [7], copyright (2015), AIP Publishing LLC and predictions by the model in this work (b).
4.1.Non-local conductivity effects
In this section, the effects of non-local and local conductivities on plasma parameters and electromagnetic wave characteristics are compared.The non-local and local conductivities are respectively obtained by solving the electron Boltzmann equation and by using the local approximation.In the simulations, the absorption power is 50 W, the RF excitation frequency is 13.56 MHz, and the pressure varies from 0.1 to 20 Pa.The conductivity is expressed asσin the following figures.
The characteristic lengths play a key role in the definition of non-local and local electrodynamics [8].If the characteristic size of plasma is much larger than electron mean free pathλbut is comparable to electron energy relaxation lengthλε, in the collisional (but nonlocal) regime the isotropic part of the electron energy distributionf0at a given point relies not only on the electric field at this point but also on plasma properties in the vicinity of the point of the size, while the anisotropic partf1is a local function of the field (i.e.local electrodynamics).Ifλbecomes comparable or even larger than the characteristic size of the plasma.In this nearly collisionless regime, thef1at a point is determined not only by the local value of the electric field at this point,but also by the profile of the electric field in the vicinity of the point of the size along the electron trajectory (i.e.nonlocal electrodynamics) [8].
Figure 3.Electron energy relaxation length and mean free path versus pressure.
The electron density and effective electron temperature versus pressure under non-local and local conductivities are shown in figure 4.In figure 4(a), electron density monotonously increases with increasing pressure, which achieves qualitatively agreement with the previous experiments [21].The electron density of non-local conductivity is higher than that of local conductivity at relatively low pressures.It also shows one of the important reasons why the fluid simulation is not valid anymore at very low pressures.The fluid simulations are based on the Ohm’s law where only the collisional heating considered is not enough to sustain discharge in lowpressure regime.For example, the electron density of nonlocal conductivity is far higher than that of local conductivity at 0.1 Pa in our simulation.As increasing pressure, the electron densities under two cases of conductivities tend to converge, which indicates the transition from non-local to local electrodynamics regime.As expected, the effective electron temperature of non-local conductivity is lower than that of local conductivity at relatively low pressures,as shown in figure 4(b).Similar to the variation of electron density,the effective electron temperatures under two cases of conductivities also tend to converge as increasing pressure.The simulation results indicate that the consideration of the nonlocal conductivity in a low-pressure regime is very important for the development of a numerical model.
Figure 4.Electron density(a)and effective electron temperature(b)versus pressure under non-local and local conductivities.The excitation frequency is 13.56 MHz and the absorption power is 50 W.
The effective electron temperature is obtained based on the integral of normalized EEPF.Figure 5 shows the normalized EEPF versus electron energy at different pressures under non-local and local conductivities.Two cases of conductivities result in a big difference in the normalized EEPF,especially at 0.1 Pa, where the local conductivity produces more electrons in the inelastic energy range and fewer electrons in the elastic energy range.The shape of EEPF is determined by the expansion coefficients in the Fourier series of the electric field, according to equation (9), where the energy diffusion coefficient is proportional to the square of the expansion coefficients of electric field.The effects of collisionless heating under non-local conductivity case allow the electromagnetic wave to couple energy to plasma more efficiently.As a result, the electric field required to deposit a given power into the plasma is smaller than that needed when only considering the collisional heating under local conductivity case.The simulation results of electric fields will be shown later.The stronger electric field under local conductivity case leads to a larger energy diffusion coefficientDein the simulations and therefore more energetic electrons penetrate into the RF field region.Consequently, local conductivity produces more electrons in the inelastic energy range.
Figure 5.Normalized EEPF versus electron energy at different pressures under non-local and local conductivities.The discharge conditions refer to figure 4.
In the simulations, the energy diffusion coefficientDeis smaller than the energy diffusion coefficient related to electron-electron collisionsDee(see equation (11)) at low electron energy.Therefore, the electron-electron interactions are the only mechanism of energy exchange for the low energy electrons that cannot reach the region of strong electric field where the heating occurs.Deeunder non-local conductivity case is larger than that under local conductivity case due to higher electron density (see figure 4(a)).It leads to more low energy electrons and meanwhile tends to Maxwellize the EEPF in elastic energy range at low pressure of 0.1 Pa under non-local conductivity case.As increasing pressure, the collisionless heating is less efficient and the thermalization of EEPF occurs.The normalized EEPFs evolve into a Maxwellian distribution and are almost identical at 20 Pa for two cases of conductivities.This is because the model-predictedDeandDeefor two cases of conductivities are almost the same at 20 Pa.
Figure 6 shows the 2D profile of the absorption power density versus pressure under non-local conductivity case.The absorption power density is taken the absolute values and plotted on a log scale to emphasize the negative power absorption.The local negative power absorption in the plasma is predicted by the model with consideration of non-local conductivity at low pressures, for example, 0.1 and 0.5 Pa,and it still exists at 1 Pa even if it is not shown in figure 6.It is well known that the negative power absorption is caused by the thermal motion of electrons.For the case of non-local conductivity, the phase difference between electric field and electron current is not constant.The energetic electrons originating from the discharge edge can be out of phase with an electric field outside the skin layer (close to the center), and lose energy to electric field, thus resulting in negative power absorption.For example, negative power absorption appears close to the discharge center outside the skin layer at 0.1 and 0.5 Pa.The simulation results confirm that there is a spatial region in the discharge where the electrons, instead of taking energy from the field,give back energy to the electromagnetic wave.The size and the number of negative and positive power absorption region depend on the power, which has been experimentally confirmed by Cunge et al[36]and Ding et al [7].At low power, there is one area of positive and one area of negative power absorption in the radial direction.As increasing power,the number of negative and positive power absorption region increases.Due to a low absorption power of 50 W used in our simulations, there is only one area of positive and one area of negative power absorption.
At higher pressures,for example,3 and 20 Pa,the region of negative power absorption vanishes.It implies that the dominant power transfer mechanism is collisional heating at higher pressures.Even if homogeneous plasma is assumed in the model, the important physical property, i.e.the local negative power absorption, is still captured by the model.It has been demonstrated in the work of Godyak and Kolobov,when comparing the 1D profile of the absorption power density obtained in the experiments with that in the calculation assuming homogeneous plasma [6].The position and amplitude of positive and negative power densities predicted by the model can be very possibly different from the reality in the experiments where the plasma is inhomogeneous.The uniformity of plasma is better at lower pressure and power.It is well known that the skin depth inversely relies on the electron density [1].The magnitude and profile of plasma density can affect the magnitude and profile of the electromagnetic fields, electron current, and absorption power density.The absorption power density can be underestimated at the center of the discharge chamber and overestimated at the boundary of the discharge chamber.In addition, the power density profile shrinks and its value becomes larger towards the boundary near the coil, due to the increase in electron density (see figure 4(a)) as increasing pressure.
The 2D profile of the absorption power density versus pressure under local conductivity case is shown in figure 7.Very different from the non-local conductivity case,the local negative power absorption is not captured by the local conductivity case at low pressures.For the case of local conductivity, the phase difference between electron current and electric field is constant.Thus,the profile of absorption power relies only on the magnitudes of electric field and local conductivity.At higher pressures, for example, 20 Pa, the 2D profile of the absorption power density under local conductivity case is almost identical with that under non-local conductivity case.It indicates that the collisional heating dominates the discharge at higher pressures.Similar to the non-local conductivity case, the amplitude of absorption power density increases with increasing pressure under local conductivity case due to the increased electron density.As expected,the difference of amplitude in the absorption power density under two cases of conductivities significantly decreases with increasing pressure.For example, for nonlocal and local conductivities the amplitudes of absorption power density are 0.62 W cm−3and 0.28 W cm−3at 0.1 Pa and 1.6 W cm−3and 1.8 W cm−3at 20 Pa.
Figure 6.2D profiles of absorption power density versus pressure under non-local conductivity.The discharge conditions refer to figure 4.Power density is taken the absolute values and plotted on a log scale to emphasize the negative power absorption.
Figure 7.2D profiles of absorption power density versus pressure under local conductivity.The discharge conditions refer to figure 4.
The 2D profiles of electric field and electron current density versus pressure under non-local conductivity are shown in figure 8.As increasing pressure, the 2D profiles of electric field and electron current density shrink towards the boundary near the coil,and the amplitude of the electric field decreases that is consistent with variation of the effective electron temperature with pressure, but the amplitude of the electron current density increases due to the increase in electron density.The non-local relation between electron current density and electric field is more pronounced at lower pressures represented by very different 2D profiles of them.The skin depthsδpredicted by the model are respectively 1.5,1.0, 0.7 and 0.4 cm at 0.1, 0.5, 3 and 20 Pa.At lower pressures,more electrons contribute to the electron current outside the skin layer and near the plasma bulk, which is a non-local behavior.As increasing pressure,the electron current tends to be determined by the local electric field.Due to homogeneous plasma assumed in the model, the 2D profiles of electron current density and electric field can be slightly different from the reality where the electron density usually has maximum in the discharge center and decreases at the boundary,especially for higher pressures.
The 2D profiles of electric field and electron current density versus pressure for local conductivity case are shown in figure 9.It is clearly shown that the electron current is determined by the local electric field in the investigated pressure range due to the assumption of Ohm’s law between electron current density and electric field.The 2D profiles of electron current density and electric field are identical to each other.The skin depthsδpredicted by the model are respectively 2.3,1.2,0.7 and 0.4 cm at 0.1, 0.5, 3 and 20 Pa.Note thatδunder local conductivity case with lower electron density (see figure 4(a)) is larger than that under non-local conductivity case (mentioned in figure 8) at lower pressures.As increasing pressure to 20 Pa,the skin depths are the same(i.e.0.4 cm)under two cases of conductivities due to the same electron density.
Figure 8.2D profiles of induced electric field (first row) and plasma current density (second row) versus pressure under non-local conductivity.The discharge conditions refer to figure 4.
The amplitudes of electric field under local and non-local(see figure 8) conductivities are very closer with each other at higher pressures and very different at lower pressures.The amplitude of the electric field under local conductivity case is far higher than that under non-local conductivity case at lower pressures.As mentioned in figure 5, for the case of non-local conductivity the effects of warm-electron and long mean free path cause the electromagnetic wave to couple energy to the plasma more efficiently.The lower electric fields are therefore required to deposit the same total power.It leads to the higher electron density and lower effective electron temperature under non-local conductivity case (see figure 4).For example, at 0.1 Pa,the amplitudes of electric field under local and non-local conductivities are respectively 12.4 V cm−1(here) and 5.02 V cm−1(see figure 8), therefore the effective electron temperatures are respectively 8.6 eV and 5.4 eV (see figure 4(b)).The consideration of non-local conductivity by solving electron Boltzmann equation in the model is very important for the prediction of non-local relation between electron current and electric field as well as the local negative power absorption in low-pressure regime.
In order to estimate the variation of heating mechanisms with pressure,the characteristic frequencies of plasmas versus pressure are presented in figure 10.These frequencies include electron-neutral collision frequencyνen,stochastic collision frequencyνstoc, and effective collision frequencyνeff=νen+νstoc.In addition, the value of the excitation frequencyωis also shown.Bothνenandνstocare obtained based on the results of the model with consideration of non-local conductivity.The collision frequencyνenalmost linearly increases with increasing pressure because of its proportionality to the density of neutral particles.νenis larger than excitation frequencyωat higher pressures (≥6 Pa), thus the Ohmic heating dominates the discharge.The stochastic collision frequencyνstocshows a slight increase with pressure.Comparingνenwithνstoc, the Ohmic heating is expected to be dominant at higher pressures asνenexceedsνstocby almost one order of magnitude.Due to almost negligible effect of stochastic heating at higher pressures,νefftends to be equivalent toνen.
Figure 9.2D profiles of induced electric field (first row) and plasma current density (second row) versus pressure under local conductivity.The discharge conditions refer to figure 4.
Figure 10.Characteristic frequencies of electrons and excitation frequency (ω )versus pressure.The characteristic frequencies include the electron-neutral collision frequency (ν en ),the stochastic heating frequency (νs toc ),and the effective collision frequency (ν eff ).The discharge conditions refer to figure 4.
Towards the lower pressures (≤3 Pa), the influence of stochastic heating is important.As revealed by the simulation results of absorption power density (see figure 6), the local negative power absorption vanishes at 3 Pa.At 0.8 Pa,νenandνstocare the same magnitude.At the lowest pressure limit of 0.1 Pa, because ofνstocfar higher thanνen,νeffsignificantly deviates fromνenand tends to be equivalent toνstoc.Therefore,the evaluation of the characteristic frequencies shows a transition of the heating mechanism from non-collisional stochastic heating to collision heating as increasing pressure from 0.1 to 20 Pa.It is necessary to consider stochastic heating in low-pressure regime.
4.2.Effective conductivity effects
Due to the importance of collisionless heating caused by the thermal motion of electrons at low pressures, the effective conductivity (equation (7)) is introduced based on the local conductivity (equation (6)) but accounts for the stochastic heating similar to collisional heating by defniing a stochastic collision frequencyνstoc.In order to evaluate the feasibility of effective conductivity at low pressures, its effects on plasma parameters and electromagnetic wave characteristics are compared with non-local conductivity case.
Figure 11.Electron density (a) and effective electron temperature (b) versus pressure under non-local and effective conductivities.The discharge conditions refer to figure 4.
Figure 12.Normalized EEPF versus electron energy at different pressures under non-local and effective conductivities.The discharge conditions refer to figure 4.
The electron density and effective electron temperature versus pressure under effective and non-local conductivities are shown in figure 11.Different from the comparison between local and non-local conductivities,i.e.figure 4,where very pronounced discrepancies of effective electron temperature and electron density at lower pressures are presented, the effective conductivity almost reproduces the effective electron temperature and electron density under non-local conductivity case in the investigated pressure range.The addition of collisionless heating to collisional heating by introducing the stochastic collision frequency effectively increases electron density compared with the case of only collisional heating taken into account in low-pressure regime.Figure 12 shows the normalized EEPF versus electron energy at different pressures under effective and non-local conductivities.As expected, the normalized EEPF under non-local conductivity case is also almost reproduced by the effective conductivity.It confirms that the effective collision frequency instead of electronneutral collision frequency taken into account in the solver of local conductivity achieves a more accurate description of electron density and effective electron temperature in low-pressure regime.
The 2D profile of absorption power density versus pressure under effective conductivity case is shown in figure 13.The effective conductivity almost reproduces the skin depth and the 2D profile of electric field under non-local conductivity case at low pressures of 0.1 and 0.5 Pa, and therefore the effective electron temperature and electron density (see figure 11).However,it does not capture the local negative power absorption that is revealed by the non-local conductivity at low pressures.Meanwhile, the non-local relation between electron current and electric field at low pressures is also not captured by the effective conductivity,as shown in figure 14.This is because despite the effective collision frequencyνeffinstead of the collisional frequencyνenused in the calculation of conductivity, the determination of electron current density still depends on the Ohm’s law.The nature of the local relation between electron current and electric field is not changed for an effective conductivity case.
Figure 13.2D profiles of absorption power density versus pressure under effective conductivity.The discharge conditions refer to figure 4.
Figure 14.2D profiles of induced electric field (first row) and plasma current density (second row) versus pressure under effective conductivity.The discharge conditions refer to figure 4.
Plasma conductivity is a fundamental parameter defining the electrical characteristics of a plasma device.It couples plasma kinetics to electromagnetic fields.In this study, a numerical model for homogeneous plasma is developed to study the conductivity effects in cylindrical inductively coupled Ar discharges at pressures of 0.1-20 Pa.The discharges experience a transition from collisionless to collisional regimes as the pressure increases from 0.1 to 20 Pa.The model consists of an electron kinetics module, electromagnetics module, and global model module.The non-local conductivity, local conductivity, and effective conductivity are respectively considered in the model and their effects on plasma parameters and electromagnetic wave characteristics are compared.The non-local conductivity is addressed by the electron Boltzmann equation.The local conductivity is defined by relating electron current to the electric field using Ohm’s law.The effective conductivity is introduced based on the Ohm’s law but with consideration of non-collisional stochastic heating effect by defining a stochastic collision frequency.
Simulation results under non-local conductivity cases are first compared with predictions under the assumption of local conductivity case (Ohm’s law).The lower effective electron temperature, as well as higher electron density, is obtained for non-local conductivity case than that for local conductivity case at relatively low pressures.This is because the collisionless heating effects under non-local conductivity case allow the electromagnetic wave to couple energy to the plasma more efficiently.Consequently,the electric field required to deposit a given power into the plasma is smaller than that needed when only considering the collisional heating under local conductivity case.The weaker electric field under non-local conductivity case leads to lower effective electron temperature and therefore higher electron density at lower pressures.As increasing pressure, the effective electron temperature and electron density,respectively,tend to converge under two cases of conductivities.It displays a transition of electron kinetics from non-local to local regime.The local conductivity predicts that there is no negative power absorption and the absorption power density monotonically decays from the boundary near the coil towards the discharge center in the investigated pressure range.However,the non-local conductivity predicts a local negative power absorption region and two peaks respectively for positive and negative power absorption regions at relatively low pressures.Almost the same 2D profiles of the absorption power density for non-local conductivity as that for local conductivity are obtained at relatively high pressures.This difference between two cases of conductivities at relatively low pressures can be understood by the prediction of the 2D profiles of electron current density and electric field.More spread-out 2D profile of electron current density than that of electric field under non-local conductivity case indicates a large fraction of electrons coming from the skin layer contributing to the electron current outside the skin layer.The 2D profiles of electron current density and electric field are coincident for local conductivity case in the pressure range of interest,but it only roughly holds true for non-local conductivity case at very high pressure.It verifies that the local relation between electron current and electric field is not valid anymore in the low-pressure regime.The non-local conductivity model is valid no matter which heating mechanisms(i.e.collisional,collisionless or both of them) dominate the discharge in a given situation.
In addition, an effective conductivity with consideration of non-collisional stochastic heating effect is introduced.The effective conductivity almost reproduces the electron density,effective electron temperature,and normalized EEPF under nonlocal conductivity case in the investigated pressure range.However, the 2D profiles of electric field and electron current density are coincident and there is no negative power absorption in the investigated pressure range.Therefore, the non-local relation between electron current and electric field is not captured by the effective conductivity.This is because the introduction of effective conductivity is still based on the Ohm’s law and therefore the nature of the local relation between electron current and electric field is not changed.Even so, for the calculations of effective electron temperature and electron density,the effective conductivity can be alternative to non-local conductivity under the non-local regime in complex modellings such as multi-dimensional fluid simulations for simplicity.
In this work,the effect of the ponderomotive force is not investigated, because it has been proven to be negligibly small at the industrial frequency 13.56 MHz [37].In future work, the model will be developed to describe the molecular gas discharges, such as hydrogen gas discharge that are extensively used in negative hydrogen ion sources operated at lower RF frequency (1 MHz).In that case, the magnetic(ponderomotive)force can be rather large,because it is almost inversely proportional to the RF frequency [34, 38].
Acknowledgments
This work was sponsored by National Natural Science Foundation of China (Nos.12105041, 11935005 and 12035003), Fundamental Research Funds for the Central Universities (No.2232020D-40) and Shanghai Sailing Program (No.20YF1401300).
Appendix.Derivation of the energy diffusion coefficient in the kinetic equation
In order to clarify the derivation of energy diffusion coefficientDεin equation (8) of the main text, the deduction is shown in detail here.A quasilinear equation for isotropic distribution functionf0is given as [12, 20]
whereEandBare respectively the electric field and magnetic field.For arbitrary physical quantities pertaining to simple harmonic oscillation with timeA(t)=Ae-iωtandB(t)=Be-iωt,by averaging their product over the oscillation periodT=we have
Therefore, the kinetic equation forf0averaged over the RF period is give as
Then,using the orthogonality conditions[20],the kinetic equation forf0averaged over the space is given as
Substituting equation (A6) into (A5), we have
where the flux in the velocity spaceΓ is given as
Using the spherical coordinates in the velocity space
we have
Averaging the∇v·Γover azimuthal angle in the velocity space, we obtain
where
Substituting equation (A15) into (A14), we have
where the functionΨ can be found in [20].
Substituting equation (A16) into (A12) and combining equation (A8), we can obtain
where
For convenience,the RHS of equation(A17)is expressed as a function of energyε.It is given as
where the energy diffusion coefficientDein equation (8) of the main text is obtained as follows
Finally, the effect of magnetic field is cancelled out, because the Lorentz force term (v×B) is averaged over azimuthal angle in the velocity space (see equation (A12)).
ORCID iDs
相关热词搜索:collisionless collisional transition
- 范文大全
- 说说大全
- 学习资料
- 语录
- 生肖
- 解梦
- 十二星座
-
2022年4月主题党日活动记录范文15篇
2022年4月主题党日活动记录范文15篇2022年4月主题党日活动记录范文篇1一个崇尚阅读的民族,必然精神饱满、意气风发、活力四射。习近平总书记强调:“学习
【活动总结】 日期:2022-08-01
-
少先队的光荣历史故事 队前教育-光辉历程
2017-2018学年队前教育1光辉历程一、劳动童子团1924——1927二、三十年代年的中国是一个
【法律文书】 日期:2020-06-23
-
家乡赋|最美的家乡赋
家乡赋 孙传志 今安康市,白河双丰镇,吾之家乡也。三环沃土,山水环抱。其北依山,山系五岭,山
【调研报告】 日期:2020-04-01
-
党支部1-12月全年主题党日活动计划表
2022年党支部主题党日活动计划表序号活动时间活动方式活动内容12022年1月专题学习研讨集中观看2022年新年贺词,积极开展学习研讨交流。组织生活会组织党员认真对照党章...
【活动总结】 日期:2022-10-14
-
【人教版1-6年级数学上册知识点精编】1-6年级数学人教版教材
人教版二年级数学上册知识点汇总第一单元长度单位一、米和厘米1、测量物体的长度时,要用统一的标准去测量
【调研报告】 日期:2020-11-08
-
2022年2月份主题党日活动记录5篇
2022年2月份主题党日活动记录5篇2022年2月份主题党日活动记录篇1尊敬的党组织:在今年的开学初,本人积极参加教研室组织的教研活动,在学校教研员的指
【活动总结】 日期:2022-08-12
-
2023年平安校园建设方案13篇
平安校园建设方案“平安校园”创建工作,我们幼儿园全体教职员工一直把它当作头等大事来抓。领导高度重视,以“平安校园”创建活动为抓手,建立和规范校园安全工作机制
【规章制度】 日期:2023-11-02
-
医院最佳主题党日活动11篇
医院最佳主题党日活动11篇医院最佳主题党日活动篇1 医院最佳主题党日活动篇2为隆重纪念中国共产党成立100周年,进一步巩固党的群众路线教育实践活动成果,切实
【活动总结】 日期:2022-10-29
-
主题党日活动记录202210篇
主题党日活动记录202210篇主题党日活动记录2022篇12021年是中国共产党成立100周年,为广泛开展爱国主义宣传教育,铭记党的历史,讴歌党的光辉历程,
【活动总结】 日期:2022-08-02
-
南京大屠杀国家公祭日悼念文案句子11篇
南京大屠杀国家公祭日悼念文案精选句子1、惟有民魂是值得宝贵的,惟有他发扬起来,中国才有真进步。——鲁迅2、我爱我的祖国,爱我的人民,离开了它,离开了他们,我
【企划文案】 日期:2023-10-20
-
正式的晚宴邀请函 公司晚宴邀请函
尊敬的先生 女士: 我公司谨定于xxxx年xx月xx日xx:xx在xxxx店隆重举行xx市xx届xxxx晚宴(宴会地址:xx区xx路xxxx) 敬请届时光临!xxxxxx集团股份有限公司xxxx有限公司敬邀xxxx年xx月xx日
【简历资料】 日期:2019-08-03
-
一年级新学期目标简短_一年级学生新学期打算
新学期到了,我是一年级下册的小学生了。 上课的时候,我要认真学习,不做小动作,认真听讲。我要认真学习,天天向上,努力学习,耳朵要听老师讲课,眼睛要瞪得大大的看老...
【简历资料】 日期:2019-10-26
-
入少先队员改正的缺点有哪些_少先队申请书
敬爱的少先队组织:我们是共产主义接班人,继承革命先辈的光荣传统,爱祖国,爱人民,鲜艳的红领巾飘扬在前胸 我叫xx,是一年级(x)班的小学生。每当听到这首...
【简历资料】 日期:2019-07-28
-
[信访复查复核制度作用探讨]信访复查复核有用吗
作为我国特有的一项制度,信访制度的出现并长期存在不是偶然的,虽然一些法学专家认为信访制度具有“人治”
【职场指南】 日期:2020-02-16
-
[党员干部2019年主题教育个人问题检视清单及整改措施2篇] 党员干部
2019年主题教育问题检视清单及整改措施根据主题教育领导小组办公室《关于认真做好主题教育检视问题整改
【求职简历】 日期:2019-11-08
-
网络维护工作内容_(精华)国家开放大学电大专科《网络系统管理与维护》形考任务1答案
国家开放大学电大专科《网络系统管理与维护》形考任务1答案形考任务1理解上网行为管理软件的功能【实训目
【职场指南】 日期:2020-07-17
-
党委会与局长办公会的区别_局长办公会制度
为进一步加强xxx局工作的规范化、制度化建设,提高行政效能,规范议事程序,特制定本制度。一、会议形式1、局长办公会议由局长、副局长参加。由局长召集和主持。根据工作需要...
【求职简历】 日期:2019-07-30
-
学生会组织部部长竞选稿5篇
学生会组织部部长竞选稿以“三制”为统领推进农村党的建设中共**市委组织部近年来,**市认真落实中央、省和徐州市委的部署,积极适应发展要求,从加强领导体制、运
【求职简历】 日期:2023-11-06
-
如何凝心聚力谋发展【坚定信心谋发展凝心聚力促跨越】
当前,清河正处于在苏北实现赶超跨越基础上全面腾飞的战略机遇期,处于在全市率先实现全面小康基础上率先实
【简历资料】 日期:2020-03-17
-
《铁拳砸碎“黑警伞”》警示教育片观后感
影片深刻剖析了广西北海市公安局海西派出所原所长张枭杰蜕变堕落的轨迹。观看警示教育片后,做为一名党员教
【简历资料】 日期:2020-08-17
-
员工的培训与奖励对于心理授权及工作满意度的影响
席雯,温红丽,张耀,吕浩(太原师范学院管理系,山西晋中030619)在体验经济时代和知识经济时代到来
【其他范文】 日期:2023-01-28
-
2023年某镇年度落实全面从严治党主体责任工作情况报告【优秀范文】
2022年以来,**镇党委坚持以习近平新时代中国特色社会主义思想为指导,深入贯彻学习党的二十大思想和习近平总书记来*考察重要讲话重要指示精神,认真执行省委、市委和区委各项...
【其他范文】 日期:2023-09-07
-
高密度互通立交出口匝道驾驶人视觉搜索行为特征
徐进孙子秋王思棋杨子邈(1重庆交通大学交通运输学院,重庆400074)(2重庆交通大学山区复杂道路环
【其他范文】 日期:2023-02-19
-
2023年在市政府第六次全体(扩大)会议上讲话
XX市“两会”昨天下午刚闭幕,今天上午,XX市政府就第一时间召开全体(扩大)会议,主要是全面学习贯彻XX市“两会”精神,对《政府工作报告》确定的目标任务进行分解,齐心协...
【其他范文】 日期:2023-08-28
-
当前农村留守孩子教育问题的几点思考:留守农村妇女人
当前随着农村劳动力的不断外流,特别是上世纪六七十年代以后出生的人已成为当前外出务工的主流。而他们留守
【其他范文】 日期:2020-02-06
-
互动式教学法应用于大专化学教学中的实践价值
道寿樱(和田地区教育学院,新疆和田848000)针对大专化学教学,互动式教学是一种必然运用趋势。按照
【其他范文】 日期:2023-01-26
-
【关于学习心得】
通过学习与讨论,我想谈谈减负与教改的关系。《纲要》指出:要切实减轻中小学生课业负担。要减轻课业负担,不需进行教学改革。对于教学改革,教师、学生包括家长都反映强烈,...
【其他范文】 日期:2019-07-05
-
农家乐项目申请报告
农家乐项目申请报告第一章农家乐项目申请报告总论一、项目背景项目名称、项目承办单位、承担编制申请报告的咨询意念及其资质、申请报告工作依据及其工作概括。二、项目单位申...
【其他范文】 日期:2023-10-01
-
2022【关于公司尾牙晚会策划方案】关于公司尾牙晚会策划案
当前位置:>>>2021-10-0320XX年尾牙晚会方案紫气满堂万象新、桃花盛开报春来。20XX年公司的发展走过了许多坎坷,崭新的20XX年已开始。在此辞旧迎新之际,为让员工欢聚一堂,共庆佳节,感受
【其他范文】 日期:2022-09-19
-
2022“之”字型排比句40例(1月8日)(完整)
1 开拓视野、解放思想之行,自我加压、负重奋进之行,坚定信心、鼓舞斗志之行,凝聚力量、加快发展之行。2 加强党性修养、筑牢信仰之基,加强政德修养、打牢从政之基,加强...
【其他范文】 日期:2022-10-09
-
军转座谈会交流发言4篇
军转座谈会交流发言4篇军转座谈会交流发言篇1大家好,我叫贺丽,2015届选调生,来自康定市委组织部,现在省委编办跟班学习。今天,非常荣幸向大家汇报我的学习收
【发言稿】 日期:2022-10-27
-
12岁生日小寿星发言4篇
12岁生日小寿星发言4篇12岁生日小寿星发言篇1各位来宾、各位朋友:大家好!今天,我们欢聚在这里,共同庆祝**十二周岁生日。首先,我代表**的父母以
【发言稿】 日期:2022-07-31
-
廉政大会总结发言稿7篇
廉政大会总结发言稿7篇廉政大会总结发言稿篇1各位领导,同志们:根据会议安排,我就党风廉政建设工作做表态发言,不妥之处,请批评指正。一、提高认识,切实
【发言稿】 日期:2022-10-30
-
【企业疫情风险控制方案】 2020企业复工疫情方案
企业疫情风险控制方案2020新冠病毒肺炎疫情防控工作总结汇报3篇 关于新型冠状病毒感染的肺炎疫
【演讲稿】 日期:2020-02-27
-
纪委书记工作表态发言4篇
纪委书记工作表态发言4篇纪委书记工作表态发言篇1在镇党委政府正确领导下,在全村干部和群众的共同努力下,紧紧围绕建设社会主义新农村工作为重点,尽职尽责,与时俱
【发言稿】 日期:2022-09-30
-
党员教育培训总结交流发言12篇
党员教育培训总结交流发言12篇党员教育培训总结交流发言篇1根据市委组织部《关于开展我市〈20XX
【发言稿】 日期:2022-12-19
-
我最敬佩的人开头_我敬佩的一个人作文20篇2020年
我敬佩的一个人作文20篇 我敬佩的一个人作文一): 我身边有很多值得我们敬佩的人,但我最敬佩的一
【发言稿】 日期:2020-11-10
-
[钻井队队长(副队长、指导员)岗位HSE应知应会试题(1863)]
钻井队队长(副队长、指导员)岗位HSE应知应会试题(判断题:771;单选题:626;多选题:466)
【贺词】 日期:2020-09-23
-
话剧《家》剧本 话剧剧本:爱的空间
找文章到更多原创-(http: www damishu cn)人物介绍:刘伟,男,32岁,某购物广
【演讲稿】 日期:2020-01-21
-
五言绝句大全500首古诗_五言绝句144首
五言绝句144首 五言绝句(一): 1《春夜喜雨》唐朝·杜甫 好雨知时节,当春乃发生。随风潜入
【祝福语】 日期:2021-03-13
-
2023年中国行政区划调整方案(设想优秀3篇
中国行政区划调整方案(设想优秀民政部第二次行政区划研讨会会议内容一、缩省的意义与原则1.意义1)利于减少中间层次中国行政区划层级之多为世界之最,既使管理成本
【周公解梦】 日期:2024-02-20
-
2023年和儿媳妇在一起幸福的句子3篇
和儿媳妇在一起幸福的句子1、假如人生不曾相遇,我还是那个我,偶尔做做梦,然后,开始日复一日的奔波,淹没在这喧嚣的城市里。我不会了解,这个世界还有这样的一个你
【格言】 日期:2023-11-10
-
XX老干局推进党建与业务深度融合发展工作情况调研报告:党建调研报告
XX老干局推进党建与业务深度融合 发展工作情况的调研报告 党建工作与业务工作融合发展始终是一个充满生
【成语大全】 日期:2020-08-28
-
中国共产党第三代中央领导集体的卓越贡献
中国共产党第三代中央领导集体的卓越贡献 --------------继往开来铸就辉煌 【摘要】改
【成语大全】 日期:2020-03-20
-
信息技术2.0能力点 [全国中小学教师信息技术应用能力提升工程试题题库及参考答案「精编」]
全国中小学教师信息技术应用能力提升工程试题题库及答案(复习资料)一、判断题题库(A为正确,B为错误)
【格言】 日期:2020-11-17
-
党建工作运行机制内容有哪些_构建基层党建工作运行机制探讨
党的基层组织是党在社会基层组织中的战斗堡垒,是党的全部工作和战斗力的基础。加强和改进县级以下各类党的
【经典阅读】 日期:2020-01-22
-
电大现代教育原理_最新国家开放大学电大《现代教育原理》形考任务2试题及答案
最新国家开放大学电大《现代教育原理》形考任务2试题及答案形考任务二一、多项选择题(共17道试题,共3
【成语大全】 日期:2020-07-20
-
集合推理_七,推理与集合
七推理与集合1 期中考试数学成绩出来了,三个好朋友分别考了88分,92分,95分。他们分别考了多少分
【名人名言】 日期:2020-12-18
-
基层党务工作基本内容_党建基本工作有哪些
党建基本工作有哪些(一) 基层党建工作包括哪些内容 选择了大学生村官这条路,你就与农村基层党
【名人名言】 日期:2020-08-06
-
【2020-2021学年高一英语外研版(2019)选择性必修第一册Unit3Faster,higher,strongerSectionⅠ导学讲义】
Unit3 Faster,higher,stronger背景导学MichaelJordan—Head
【歇后语】 日期:2021-04-19
-
关于三农工作重要论述心得体会3篇
关于三农工作重要论述心得体会3篇关于三农工作重要论述心得体会篇1习近平总书记指出:“建设现代化国家离不开农业农村现代化,要继续巩固脱贫攻坚成果,扎实推进乡村
【学习心得体会】 日期:2022-10-29
-
【福生庄隧道坍塌处理方案】 福生庄隧道在哪里
(呼和浩特铁路局大包电气化改造工程指挥部,内蒙古呼和浩特010050)摘要:文章介绍了福生庄隧道
【学习心得体会】 日期:2020-03-05
-
五个一百工程阅读心得体会13篇
五个一百工程阅读心得体会13篇五个一百工程阅读心得体会篇1凡益之道,与时偕行。在全国网络安全和信
【学习心得体会】 日期:2022-12-07
-
城管系统警示教育心得体会9篇
城管系统警示教育心得体会9篇城管系统警示教育心得体会篇1各党支部要召开多种形式的庆七一座谈会,组织广大党员进行座谈,回顾党的光辉历程,畅谈党的丰功伟绩,
【学习心得体会】 日期:2022-10-09
-
发展对象培训主要内容10篇
发展对象培训主要内容10篇发展对象培训主要内容篇1怀着无比激动的心情,我有幸参加了__新区区委党校20__年第四期(区级机关)党员发展对象培训班。这次的学习
【培训心得体会】 日期:2022-09-24
-
凝聚三种力量发展全过程人民民主心得体会12篇
凝聚三种力量发展全过程人民民主心得体会12篇凝聚三种力量发展全过程人民民主心得体会篇1新民主主义革命是指在帝国主义和无产阶级革命时代,殖民地半殖民地国家中的
【学习心得体会】 日期:2022-08-31
-
2022年全国检察长会议心得7篇
2022年全国检察长会议心得7篇2022年全国检察长会议心得篇1眼睛是心灵上的窗户,我们通过眼睛才能看到世间万物,才能看到眼前这美好的一切。拥有一双明亮的眼
【学习心得体会】 日期:2022-10-31
-
在街道深化作风建设推动高质量发展走在前列动员会上讲话
在2023年街道深化作风建设推动高质量发展走在前列动员会上的讲话同志们:今天我们召开“街道深化作风建设推动高质量发展走在前列动员会”,这次会议是街道三季度召开的第一场...
【军训心得体会】 日期:2024-03-17
-
全面从严治党的心得体会800字7篇
全面从严治党的心得体会800字7篇全面从严治党的心得体会800字篇1中国特色社会主义是我们党领导
【学习心得体会】 日期:2022-12-14
-
2月教师党员个人思想汇报5篇
2月教师党员个人思想汇报敬爱的党组织:最近这一个月的时间对于我来说是极不平凡的,在这段时间里我认真学习了文化部网上党校的相关内容,经过长达40小时的
【教师心得体会】 日期:2023-10-15
-
2024年主题教育民主生活会批评与自我批评意见(38条)(范文推荐)
2023年主题教育民主生活会六个方面个人检视、相互批评意见:1 理论学习系统性不强。学习习近平新时代中国特色社会主义思想不深不透,泛泛而学的时候多,深学细照的时候少,特...
【邓小平理论】 日期:2024-03-19
-
2024年交流发言:强化思想理论武装,增强奋进力量(完整)
习近平总书记指出:“一个民族要走在时代前列,就一刻不能没有理论思维,一刻不能没有思想指引。”党的十八大以来,伴随着新时代中国特色社会主义思想在实践中形成发展的历程...
【三个代表】 日期:2024-03-19
-
2024年度镇年度县乡人大代表述职评议活动总结
xx镇20xx年县乡人大代表述职评议活动总结为响应县级人大常委会关于开展县乡两级人大代表述职评议活动,进一步激发代表履职活力,加强代表与人民群众的联系,提高依法履职水平...
【马克思主义】 日期:2024-03-19
-
“千万工程”经验学习体会(研讨材料)
“千万工程”是总书记在浙江工作时亲自谋划、亲自部署、亲自推动的一项重大决策,也是习近平新时代中国特色社会主义思想在之江大地的生动实践。20年来,“千万工程”先后经历...
【三个代表】 日期:2024-03-19
-
2024年在市政协机关工作总结会议上讲话
同志们:刚才,XX同志对市政协机关20XX年工作进行了很好的总结,很精炼,很到位,可以感受到去年机关工作确实可圈可点。XX同志宣读了表彰决定,机关优秀人员代表、先进集体代...
【邓小平理论】 日期:2024-03-18
-
在全区防汛防涝动员暨河长制工作推进会上讲话提纲【完整版】
区长,各位领导,同志们:汛期已经来临,我区城区防涝工作面临强大考验,形势不容乐观。年初,区城区防涝排渍指挥部已经召开专题调度会,修订完善应急预案,建立网格化管理机...
【马克思主义】 日期:2024-03-18
-
2024年镇作风整治工作实施方案(完整文档)
XX镇作风整治工作实施方案为深入贯彻落实党的二十大精神及省市区委深化作风建设的最新要求,突出重点推进干部效能提升,坚持不懈推动作风整治工作纵深发展,根据《关于印发《2...
【毛泽东思想】 日期:2024-03-18
-
2024市优化法治化营商环境规范涉企行政执法实施方案【优秀范文】
xx市优化法治化营商环境规范涉企行政执法实施方案为持续优化法治化营商环境,激发市场主体活力和社会创造力,规范行政执法行为,创新行政执法方式,提升行政执法质效,着力解...
【毛泽东思想】 日期:2024-03-18
-
2024年度关于开展新一轮思想状况摸底排查工作通知(完整)
关于开展新一轮思想状况摸底排查工作的通知为深入贯彻落实关于各地开展干部职工思想状况大摸底大排查情况上的批示要求和改革教育第二次调度会议精神,有针对性做好队伍教育管...
【三个代表】 日期:2024-03-18
-
2024年公路养护中心主任典型事迹材料(完整文档)
“中心的工作就是心中的事业”——公路养护中心主任典型事迹材料**,男,1976年6月出生,1993年参加工作,2000年4月调入**区交通运输局工作,大学本科学历,中共党员,现任**...
【马克思主义】 日期:2024-03-17