O2中α型射频放电的理论研究

根据射频放电中电流连续性方程、稳态时电子数密度的连续性方程以及电子的能量平衡方程,建立了气体放电氧碘激光器中α型射频放电等离子体的理论模型,通过数值求解得到了射频放电等离子体中电场、电子数密度的空间分布,分析了放电参数对放电特性的影响.结果表明采用频率高的射频放电会使放电空间电场降低,电子数密度增加,从而有利于单重氧的生成,为提高气体放电中单重氧的产率提供了理论依据.     

射频容性耦合Ar/O_2等离子体的轴向诊断

实验通过朗缪尔探针辅助激光诱导解离技术对27.12 MHz驱动的不同含氧量条件下容性耦合Ar等离子体进行了诊断研究.通过质量流量计改变通入Ar与O_2的流量以得到不同含氧量的等离子体.结果表明,由于氧气的分解吸附反应需要消耗电子,致使朗缪尔探针测得的电子能量概率函数(EEPF)的中能部分随着含氧量的上升而下降.射频输入功率增加时电子密度的上升引起了电子-电子碰撞热化,从而使EEPF由Druyvesteyn分布向麦克斯韦分布转变.在功率电极附近,由于鞘层边界附近的电子氧气分子碰撞时的分解吸附反应使得鞘层区附近的负离子密度较高.另外,负离子密度沿着轴向呈现勺型分布的特征.这主要是由于负离子在双极电场作用下向等离子体放电中心扩散的过程中所存在的负离子产生与损失的反应过程导致.     

双频容性耦合CF4等离子体放电的二维混合模拟

研究了等离子体反应装置内的等离子体密度、粒子能量与角度分布等参量在装置径向与轴向上的分布特性。在研究过程中应用二维混合模型对CF4气体放电进行模拟。计算结果显示:在电极表面与侧壁附近的鞘层区特性有明显的区别。由于装置侧壁处受电源产生的射频电场的影响较小,侧壁处的鞘层主要由双极扩散机制形成,其产生的径向电场强度较弱,鞘层厚度也较薄。而在电极附近,由于受到射频电场的影响,鞘层的厚度显著增加,指向电极方向的轴向电场强度也远大于指向侧壁方向的径向电场强度。在电极区域内,离子通量分布均匀;在电极边缘与侧壁的间隙内,因电场强度减小,离子通量则发生迅速衰减。在射频电极覆盖的范围内离子能量分布大体上保持不变,电极与侧壁的交界处,由于受到侧壁处径向电场的影响,离子能量分布略有不同。在放电装置的中心区域,入射到电极上的离子角度分布基本保持一致,而在电极边界与装置侧壁的交界处,由于径向电场的影响,离子的垂直入射角增加,以大角度轰击电极的离子数量也显著增加。     

双频容性耦合等离子体相分辨发射光谱诊断

采用相分辨发射光谱法, 对双频容性耦合纯Ar和不同含O₂量的Ar-O₂混合气体放电等离子体的鞘层激发模式进行了探究. 在射频耦合电源上极板的鞘层区域处观察到两种电子激发模式: 鞘层扩张引起的电子碰撞激发模式和二次电子引起的电子碰撞激发模式; 并发现这两种激发模式均受到低频射频电源周期的调制. 在纯Ar放电等离子体中, 两种激发模式的激发轮廓相似; 而在Ar-O₂混合气放电等离子体中, 随着含O₂量的增加, 二次电子的激发轮廓变弱. 此外, 利用相分辨发射光谱法对不同含O₂量的Ar-O₂混合气放电下Ar的 750.4 nm谱线的平均低频电源周期轴向分布进行了研究, 得到了距耦合电源上极板约3.8 mm处为双频容性耦合射频等离子体的鞘层边界. 关键词： 双频容性耦合等离子体 等离子体鞘层 发射光谱     

ICP等离子体鞘层附近区域发光光谱特性分析

为了独立控制鞘层附近区域离子密度和离子能最分布,采用光发射谱(OES)测量技术,对不同射频功率、放电气压和基底偏压下感应耦合等离子体鞘层附近区域辉光特性进行了研究.原子谱线和离子谱线特性分析表明,在鞘层附近区域感应耦合等离子体具有较高的离子密度和较低的电子温度.改变放电气压和射频功率,对得到的光谱特性分析表明,鞘层附近区域离子密度随射频功率的增大而线性增大,在低压下随气压的升高而增大.低激发电位原子谱线强度增加迅速,高激发电位原子谱线强度增加缓慢,而离子谱线强度增加很不明显.改变基底直流偏压,对得到的发射光谱强度变化分析表明,谱线强度随基底正偏压的增加而增大.随着基底负偏压的加入,谱线强度先减小而后增大;直流偏压为-30 V时,光谱强度最弱.快速离子和电子是引起Ar激发和电离过程的主要能量来源.     

微腔等离子体放电过程的数值研究

采用2维自洽完全流体模型,数值研究了阳极为通孔的高气压微腔放电结构中等离子体参数的变化过程。模拟结果获得了当氩气压强为13.3 kPa时,放电中的电势分布、等离子体密度分布、径向电场分布和电子温度分布等重要参数的演化过程。模拟结果表明在放电过程中,阴极附近的电场由轴向电场逐步转变为径向电场,等离子体密度最大值位于放电腔中间处,并随时间推移由阳极附近向阴极附近移动,电子温度的最大值出现在阴极环形鞘层区域。     

��������Ƶ���Ƶ������Ϸŵ�������������ıȽ�

开发了氢气甚高频(60MHz)容性耦合放电的PIC/MC模型.在模型中考虑了带电粒子(e,H+,H2+,H3+)与H2的21种碰撞反应过程,模拟了氢气甚高频放电射频电场和电势分布以及电子和氢离子(H+,H2+,H3+)粒子密度和平均能量分布,并与频率为13.56MHz的放电结果进行了比较.结果表明,相对于频率13.56MHz的放电,氢气甚高频放电等离子体电势增高,导致两电极附近的电场增强;另外,两鞘层厚度变窄并且电子和H3+离子平均能量减小,其总密度却增加.H3+离子为氢气甚高频放电空间的主要离子,H2+离子密度比H3+离子小约2～3个数量级.     

LEAD螺旋波射频等离子体电子能量分布函数的测量

通过电压扫描探针测量了LEAD装置螺旋波射频等离子体伏安特性曲线和电子能量分布函数,并对等离子体基本参数剖面如电子温度、密度、等离子体电位和鞘电位降系数等进行了计算.实验发现在LEAD装置半径5cm以内的位置电子能量分布函数都遵循麦克斯韦分布.而在半径7cm附近,即与射频源线圈处于相同径向位置,电子能量呈现出典型的双麦...     

直流空心阴极放电中鞘层区电子的输运过程

采用蒙特－卡洛方法对空心阴极放电中鞘层区电子的输运过程进行了研究，电子在鞘层区被非均匀电场加速，两次碰撞之间的步长是由电子与中性粒子的碰撞频率确定。模型中三种碰撞截面积是由实验和理论数据拟合而来。研究了电子平均能量、电子密度和电离系数在径向的分布。电子能量的空间分布结果与实验很好地符合。 关键词：     

碰撞效应对入射到射频偏压电极上离子能量分布和角度分布的影响

考虑了离子与中性粒子的弹性碰撞和电荷交换碰撞效应,建立了一套描述射频等离子体鞘层动力学特性的自洽模型,并利用MonteCarlo模拟方法研究了入射到电极上的离子的能量分布和角度分布.数值结果表明:随着放电气压增加,入射到电极上离子的能量分布逐渐地由双峰分布变成单峰分布,而且低能离子的数目也逐渐地增加.入射到电极上的离子呈小角分布,而且放电气压等参数对角度分布的影响不是太明显. 关键词： 射频放电 等离子体 离子 鞘层     

Tonks-Langmuir problem for a bi-Maxwellian plasma

An analytical solution of the Tonks-Langmuir (TL) problem with a bi-Maxwellian electron energy distribution function (EEDF) is obtained for a plasma slab. The solution shows that the ambipolar potential, the plasma density distribution, and the ion flux to the wall are mainly governed by the cold electrons, while the ionization rate and voltage drop across the wall sheath are governed by the hot electrons. The ionization rate by direct electron impact is found to be spatially rather uniform, contrary to the T-L solution where it is proportional to the plasma density distribution. The temperature of hot electrons defined by the ionization balance is found to be close to that of the T-L solution for a mono-Maxwellian EEDF, and is in reasonable agreement with experiments carried out in a low pressure capacitance RF discharge. The energy balance for cold electrons in this discharge shows that their heating by hot electrons via Coulomb interaction is equalized by the cold electrons' escape to the RF electrodes during collapse of the RF sheath     

Modeling of plasma-etch processes using well stirred reactorapproximations and including complex gas-phase and surfacereactions

A 0-D or well stirred reactor model determines spatially and time-averaged species composition in plasma-etch reactors, through solution of species, mass, and electron-energy balance equations. The use of well stirred reactor approximations reduces the computational expense of detailed kinetics calculations and allows investigation of the dependence of plasma chemistry on etch-process parameters. The reactor is characterized by a chamber volume, surface area, net mass flow or residence time, pressure, energy loss to surroundings, power deposition, and inlet-gas composition. The electron-energy equation includes a detailed power balance with losses to ions and electrons through the sheath, as well as inelastic and elastic collision losses. The model employs reaction-rate coefficients for electron-impact reactions, which require an assumption of the electron energy distribution function (EEDF). We compare model results using Maxwellian EEDF's, as well as reaction-rate coefficients determined as a function of average electron energy through solution of the Boltzmann equation, for chlorine chemistry. The Boltzmann rates are determined by time-lagging the equilibration of electrons with applied electric fields. The Maxwellian reaction rates give higher ionization fractions than the Boltzmann rates, affecting the predicted electronegativity and positive ion composition for chlorine plasmas. The model also shows a strong sensitivity of the plasma composition to the assumed surface-recombination probability of atomic chlorine     

Modeling of the electrostatic corona discharge reactor

A model is presented for the electrostatic corona discharge reactor (ECDR) in a pin-plate configuration. The main objective is to describe the fundamental chemistry and physics governing the discharge behavior and to predict the ECDR performance under various operating conditions. The electric field strength is estimated assuming a space-charge-free field. A two-term spherical harmonic expansion is used to solve the Boltzmann equation for the electron energy distribution function (EEDF) and calculate the electron-molecule reaction rates using collision cross-section data. Species continuity equations are solved for the dry and wet air systems to predict ozone and NO_x at various feed flow rates (1630, 4890, 14, 670 cm/s) and an applied voltage of 10 kV. Among the various results reported, it is noted that the calculations indicate the Maxwell EEDF cannot be used because it overpredicts the electron-molecule rate coefficients by several orders of magnitude     

射频低压等离子体电子能量分布函数的探针诊断

本文对于无内电极放电的射频低压等离子体,采用探针电流调制法,设计并建立了三探针诊断电子能量分布函数的测量系统,测量了压力在10^-3—10^-1Torr下氮气等离子体的电子能量分布函数。从理论和实验上研究了探针鞘层上射频干扰电压对测量电子能量分布函数的影响,给出了确定射频干扰值以及对射频干扰的影响进行修正的方法,并采用该方法对实验数据进行了修正,得到了正确的电子能量分布函数。根据所得到的电子能量分布函数计算的电子平均能量,与由探针伏安特性计算得到的电子平均能量相差不超过5％,该结果间接证明了本文实验测量系统的可靠性以及对射频干扰影响进行修正的正确性。     

Dust particles in collisionless plasma sheath with arbitrary electron energy distribution function

Dust particles often appear in industrial plasmas as undesirable product of the plasma-wall interactions. Large particles of several micrometers in diameter are concentrated in a thin layer (the sheath) above the lower electrode of the rf driven parallel plate device, where the electric force is strong enough to compensate particle’s gravity. Experimental and theoretical uncertainties are significantly increased in the plasma sheath. Common models of dust charging in the plasma sheath suppose the Maxwellian electron energy distribution function (EEDF) in conjunction with a flux of cold ions satisfying classical Bohm criterion at the sheath edge. In this paper we generalize this model to arbitrary EEDF with adapted Bohm criterion. We limit our considerations to collisionless or slightly collisional plasma, where the EEDF inside the sheath is expressed through the EEDF in the plasma bulk. Derived theoretical formulas are incorporated into numerical model, describing collisionless radio frequency (rf) plasma sheath together with the electrical charge, various kinds of forces, balancing radius and oscillation frequency of particles.     

A hybrid Monte Carlo/fluid model of RF plasmas in aSiH4/H2 mixture

A fluid model for simulating the capacitively coupled steady RF plasmas in a monosilane-hydrogen mixture under plasma-processing conditions is developed. This model includes Monte Carlo treatment for electron transport properties in the mixture. The Monte Carlo simulation shows that the electron transport phenomena in the nonuniform RF field of 10 MHz are not in equilibrium with the local electric field. The nonequilibrium transport properties are incorporated in the fluid model (hybrid model) by modifying the equilibrium values of the swarm parameters using data from the Monte Carlo simulation. Using this model, the spatiotemporal variations of the charged species and the electric field in the sheath region of the RF plasma are calculated. For obtaining the steady RF plasma structure, ion-induced slow processes such as recombination and diffusion of ions are calculated by combining the hybrid model with rate equations for ions. Using the calculated steady RF plasma structure, a preliminary calculation of the silyl (SiH ₃) and hydrogen (H) radical distributions caused by the generation, diffusion, and reaction of the radicals are carried out. The effect of sticking on the profile of the radical distribution is presented     

Numerical simulation of super-short pulsed discharge in Helium with particle-in-cell Monte--Carlo collisions technique

This paper reports that a simulation of glow discharge in pure helium gas at the pressure of 1.333×10³ Pa under a high-voltage nanosecond pulse is performed by using a one-dimensional particle-in-cell Monte Carlo collisions (PIC--MCC) model. Numerical modelling results show that the cathode sheath is much thicker than that of anode during the pulse discharge, and that there exists the phenomenon of field reversal at relative high pressures near the end of the pulse, which results from the cumulative positive charges due to their finite mobility during the cathode sheath expansion. Moreover, electron energy distribution function (EEDF) and ion energy distribution function (IEDF) have been also observed. In the early stage of the pulse, a large amount of electrons can be accelerated above the ionization threshold energy. However, in the second half of the pulse, as the field in bulk plasma decreases and thereafter the reverse field forms due to the excessive charges in cathode sheath, although the plasma density grows, the high energy part of EEDF decreases. It concludes that the large volume non-equilibrium plasmas can be obtained with high-voltage nanosecond pulse discharges.     

Electron temperature in a decaying krypton plasma in the presence of a low electric field

The influence of low electric fields on the average electron energy in an afterglow krypton plasma is studied by means of probe diagnostics and theoretical analysis. It is shown that, when the average electron energy is lower than the energy corresponding to the minimum scattering transport cross section, the degree of plasma ionization substantially affects the shape of the electron energy distribution function (EEDF). The nonequlibrium character of the EEDF results in the density dependence of the coefficient of ambipolar diffusion, which leads to a change in the radial profile of the charged particle density, an increase in the drop in the ambipolar potential across the plasma, and an increase in the rate of diffusive plasma decay. These effects substantially enhance the diffusive cooling of electrons, which is probably a decisive factor influencing the electron energy balance in high-Z noble gases.     

弱电离大气等离子体电子能量分布函数的理论研究

使用球谐展开的方法求解玻尔兹曼方程,得到了弱电离大气等离子体(79％氮气和21％的氧气)的电子能量分布函数(EEDF).发现当约化电场较小时(E/N＜100 Td),EEDF在2-3 eV急剧下降,在此情况下,高能尾部比麦氏分布要小;当约化电场增加,E/N＞ 400 Td,分布函数趋近于麦氏分布;当约化电场进一步增加,E/N＞ 2000 Td,EEDF的高能尾部(超过200 eV)相对于麦氏分布增加,在高频场作用下,EEDF更倾向于麦氏分布.当ω》vm时,有效电子温度只依赖于E/ω,而与碰撞频率无关;当ω《vm时,有效电子温度只依赖于E/N,与微波频率无关.与一些单原子分子等离子体中电子-电子碰撞在电离度大于10-6时就会影响EEDF不同,空气等离子体中,只有当电离度大于0.1％时,电子-电子碰撞才会对EEDF有明显影响.     

Langmuir Probe in Non‐Maxwellian Plasma: Correlation Between the Electron Energy Distribution Function,the Ion Energy at the Sheath Edge and the Floating Voltage

This work is devoted to the study of the Langmuir probe in non‐Maxwellian plasma, assuming mono‐energetic singly charged ions and a collisionless sheath. Using a general analytical equation for the Electron Energy Distribution Function (EEDF), we study the effect of the EEDF profile on: a) The ion energy at the sheath edge of a negatively biased collector, b) the I‐V probe characteristic and c) the floating voltage (V_p‐Vf). Different methods are used and compared to determine these parameters or characteristics. A correlation is given between the floating voltage, the ion energy at the sheath edge and the EEDF profile. The study is also extended to distribution functions with several components. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)