Self—Consistent Description of nitrogen dc Glow Discharge

A self-consistent hybrid Monte Carlo Fluid model is presented to describe the nitrogen dc glow discharge.The movement of fast electrons is simulated by the Monte Carlo method while the dynamics of slow electrons and ions is by fluid equations.The spatial features of the charged species and the corresponding electric field throughout the discharge have been calculated,which include the creation rates of ions and slow electrons,densities of the charged species,the electric field and the potential distribution.These closely related results can give a selfconsistent explanation of the discharge characteristics throughout the space of nitrogen dc glow discharge.The calculated ion density is also compared with the corresponding experimental result.     

ISSDE: A Monte Carlo implicit simulation code based on Stratonovich SDE approach of Coulomb collision

A Monte Carlo implicit simulation program, Implicit Stratonovich Stochastic Differential Equations(ISSDE), is developed for solving stochastic differential equations(SDEs) that describe plasmas with Coulomb collision. The basic idea of the program is the stochastic equivalence between the Fokker–Planck equation and the Stratonovich SDEs. The splitting method is used to increase the numerical stability of the algorithm for dynamics of charged particles with Coulomb collision. The cases of Lorentzian plasma, Maxwellian plasma and arbitrary distribution function of background plasma have been considered. The adoption of the implicit midpoint method guarantees exactly the energy conservation for the diffusion term and thus improves the numerical stability compared with conventional Runge–Kutta methods. ISSDE is built with C++ and has standard interfaces and extensible modules. The slowing down processes of electron beams in unmagnetized plasma and relaxation process in magnetized plasma are studied using the ISSDE, which shows its correctness and reliability.     

Particle-in-Cell/Monte Carlo Collision simulation of planar DC magnetron sputtering

In this paper a numerical simulation of a planar DC magnetron discharge is performed with the Particle-in Cell/Monte Carlo Collision （PIC/MCC） method. The magnetic field used in the simulation is calculated with finite element method according to experimental configuration. The simulation is carried out under the condition of gas pressure of 0.665 Pa and voltage magnitude of 400V. Typical results such as the potential distribution, charged particle densities, the discharge current density and ion flux onto the target are calculated. The erosion profile from the simulation is compared with the experimental data. The maximum erosion position corresponds to the place where the magnetic field lines are parallel to the target surface.     

Determination of event start time at BESⅢ

The principle of the method for the BESⅢ event start time determination and the code construction are described. The investigation of influence of the noise, the method of rejecting noise and the performance checking by the Monte Carlo physics data sample are also presented. The preliminary results got from the Monte Carlo simulation are presented, the reconstruction efficiency of J/ψ→ anything events at noise level 0~60% can achieve above 99%, and the error rate is below 1%.     

Microwave propagation with the gas breakdown

By combining the microwave propagation theory and the gas breakdown theory, the microwave propagation with the gas breakdown is analyzed theoretically. Particle-in-cell/Monte Carlo collision （PIC/MCC） simulations are carried out to verify the theoretical results. Based on this theoretical method, the breakdown phenomenon of the pulse microwave is analyzed. The results show that the product values of the initial electron density and the propagation length are the criterion to distinguish the pulse peak decline breakdown and the pulse width reduction breakdown. Furthermore, the energy transmission is also studied, which shows that the total output energy is approximately independent of the input electric field if the electric field is not extremely large.     

Charged pion condensation in anti-parallel electromagnetic fields and nonzero isospin density

The formation of charged pion condensate in anti-parallel electromagnetic fields and in the presence of the isospin chemical potential is studied in the two-flavor Nambu-Jona-Lasinio model.The method of Schwinger proper time is extended to explore the quantities in the off-diagonal flavor space,i.e.the charged pion.In this framework,π^± are treated as bound states of quarks and not as point-like charged particles.The isospin chemical potential plays the role of a trigger for charged pion condensation.We obtain the associated effective potential as a function of the strength of the electromagnetic fields and find that it contains a sextic term which possibly induces a weak first order phase transition.The dependence of pion condensation on model parameters is investigated.     

Design studies on the 4π γ-ray calorimeter for the ETF experiment at HIRFL-CSR

A high detection efficiency calorimeter which is used to detect γ-rays with energies from 1 MeV up to 10 MeV as well as light charged particles has been proposed. Design of the geometry, results of the crystal tests and Monte Carlo simulations are presented in this paper. The simulation results confirm that the calorimeter can obtain high detection efficiency and good energy resolution with the current designed geometry. And the calorimeter is competent for the future External Target Facility (ETF) experiments.     

MONTE CARLO STUDIES OF RADIATION DAMAGE INDUCED BY FUSION DYNAMIC SIMULATION OF SPUTTERING

In this paper, the sputtering caused by fusion of ⁴He particles in the first wall materials is investigated by using the dynamic Monte Carlo simulation, which is based on the binary collision approximation. The dependences of sputtering yields on the incident energy and angle, as well as the comparisons of results calculated using the Monte Carlo methed with results from experiments, are discussed. Energy spectrum and angular distributions of sputtered species, the depth of origin of the sputtered particles as well as range distributions of incident ions are given. From a comparison between the related experimental data and the calculated results with static and dynamic simulation it has been found that a better agreement is obtained for dynamic simulation.     

Computer simulation of the collision frequency of two particles in optical tweezers

Optical tweezers have been successfully used in the study of colloid science. In most applications people are concerned with the behaviour of a single particle held in the optical tweezers. Recently, the ability of the optical tweezers to simultaneously hold two particles has been used to determine the stability ratio of colloidal dispersion. This new development stimulates the efforts to explore the characteristics of a two-particle system in the optical tweezers.An infinite spherical potential well has been used to estimate the collision frequency for two particles in the optical trap based on a Monte Carlo simulation. In this article, a more reasonable harmonic potential, commonly accepted for the optical tweezers, is adopted in a Monte Carlo simulation of the collision frequency. The effect of hydrodynamic interaction of particles in the trap is also considered. The simulation results based on this improved model show quantitatively that the collision frequency drops down sharply at first and then decreases slowly as the distance between the two particles increases. The simulation also shows how the collision frequency is related to the stiffness of the optical tweezers.     

Numerical study on the characteristics of nitrogen discharge at high pressure with induced plasma

Based on the fluid theory of plasma, a model is built to study the characteristics of nitrogen discharge at high pressure with induced argon plasma. In the model, species such as electrons, N2+, N4+, Ar+, and two metastable states (N 2(A3∑u+), N2 (a1 ∑u-)) are taken into account. The model includes the particle continuity equation, the electron energy balance equation, and Poisson抯equation. The model is solved with a finite difference method. The numerical results are obtained and used to investigate the effect of time taken to add nitrogen gas and initially-induced argon plasma pressure. It is found that lower speeds of adding the nitrogen gas and varying the gas pressure can induce higher plasma density, and inversely lower electron temperature. At high-pressure discharge, the electron density increases when the proportion of nitrogen component is below 40%, while the electron density will keep constant as the nitrogen component further increases. It is also shown that with the increase of initially-induced argon plasma pressure, the density of charged particles increases, and the electron temperature as well as the electric field decreases.     

电子回旋共振放电的电离特性PIC/MCC模拟(Ⅰ)——物理模型与理论方法

采用粒子模拟与蒙特卡罗相结合(PIC/MCC)的方法对电子回旋共振(ECR)放电中的电离过程进行了模拟，其中带电粒子与微波的相互作用由PIC方法的电磁模型描述，粒子间的碰撞过程由MCC方法描述.考虑的碰撞类型有电子与中性粒子的弹性、激发、电离碰撞，离子与中性粒子的弹性、电荷交换碰撞，碰撞截面均依赖于能量而变化.阐述了理论分析的过程，为数值模拟ECR放电奠定了基础. 关键词： 电子回旋共振放电 粒子模拟 蒙特卡罗 电离     

电子回旋共振放电的电离特性PIC/MCC模拟(Ⅱ)——数值模拟与结果讨论

采用粒子模拟与蒙特卡罗相结合(PIC/MCC)的方法，应用电磁模型，编写了准三维的电子回旋共振(ECR)放电电离过程的模拟程序，得到了ECR放电过程中电子与离子的相空间分布、电磁场分布.通过对这些分布随时间演化的分析，得出ECR加热发生在ω≈ω_c0且垂直于轴向的区域；ECR区域，微波能量几乎全部耦合给电子，获得能量的电子通过与中性粒子的电离碰撞产生了大量的带电粒子；随着放电的进行，大量带电粒子通过频繁的碰撞，分布由各向异性逐渐趋于各向同性. 关键词： 电子回旋共振放电 粒子模拟 蒙特卡罗 电离     

电子回旋共振放电中电子与微波互作用特性的粒子模拟和蒙特卡罗碰撞模拟

对电子回旋共振（ECR）放电电离过程中的电子与微波互作用特性进行了理论分析与数值模拟.采用粒子模拟（PIC）方法描述带电粒子与微波的互作用,采用蒙特卡罗碰撞（MCC）方法描述粒子间碰撞过程及带电粒子与边界的相互作用.编写了准三维的PIC/MCC数值模拟程序,并对放电过程中电子能量与微波场随时间、空间的演化进行了数值诊断. 关键词： 电子回旋共振放电 粒子模拟 蒙特卡罗方法 电离     

基于负氢离子源的全三维PIC/MCC模拟算法研究

在分析负氢离子源中等离子体物理机理基础下, 研究并优化粒子模拟算法, 设计高效的粒子存储方法. 研究并运用粒子碰撞蒙特卡罗方法, 考虑等离子体势以及带电粒子间库仑碰撞, 研制了全三维粒子模拟/蒙特卡罗算法(PIC/MCC). 采用磁荷模型, 运用时域有限差分方法计算多峰磁场, 并结合国外负氢离子源JT-60U, 考虑负氢离子源中主要反应, 对全三维PIC/MCC模拟算法模拟验证.     

ECR微波放电氩离子输运过程的研究

建立了电子回旋共振（ＥＣＲ）微波放电等离子体中离子输运过程的蒙特卡罗模型，考虑了离子与中性原子的电荷交换碰撞和弹性碰撞，以及精确依赖于离子能量的电荷交换和动量转移截面，模拟了源于氩气ＥＣＲ微波放电的氩离子向衬底输运的过程，得到与实验报道相符的模拟结果。     

Influence of an acetylene impurity on the properties of a radio-frequency gas discharge in argon

The influence of an acetylene fraction in a mixture with argon on the parameters of a radio-frequency gas discharge has been investigated using the kinetic simulation by the particle-in-cell Monte Carlo collision method. Special attention has been focused on analyzing the formation of ions and radicals that are precursors for growth of nanoparticles in the gas discharge volume. It has been shown that the concentrations of electrons, positive ions, and negative ions are nonmonotonic functions of the acetylene amount in the mixture. The presence of negative ions with a concentration comparable to the electron concentration changes the parameters of the radio-frequency gas discharge only insignificantly.     

Particle-in-cell charged-particle simulations, plus Monte Carlocollisions with neutral atoms, PIC-MCC

Many-particle charged-particle plasma simulations using spatial meshes for the electromagnetic field solutions, particle-in-cell (PIC) merged with Monte Carlo collision (MCC) calculations, are coming into wide use for application to partially ionized gases. The author emphasizes the development of PIC computer experiments since the 1950s starting with one-dimensional (1-D) charged-sheet models, the addition of the mesh, and fast direct Poisson equation solvers for 2-D and 3-D. Details are provided for adding the collisions between the charged particles and neutral atoms. The result is many-particle simulations with many of the features met in low-temperature collision plasmas; for example, with applications to plasma-assisted materials processing, but also related to warmer plasmas at the edges of magnetized fusion plasmas     

�����Թ�ŵ����������̵�PIC/MCCģ��

采用粒子模拟与蒙特卡罗模型相结合(PIC/MCC)的方法,建立了氮分子气体辉光放电自洽的混合模型,其中带电粒子在电场中的运动及其产生的自洽场由PIC方法的静电模型描述,粒子间的碰撞过程由MC方法描述。模拟了放电过程及带电粒子(e、N2+、N+)在整个放电空间的行为及与之对应的自洽电场和电势的分布,通过计算带电粒子在负辉区的几率分布函数,讨论了带电粒子(e、N2+、N+)在负辉区的行为和分布特征。计算的分子离子密度与实验结果一致。     

圆柱形阳极层霍尔等离子体加速器的质点网格方法模拟

对圆柱形阳极层霍尔加速器内的放电等离子体运用二维质点网格方法（particle in cell）进行数值模拟,用蒙特卡罗碰撞方法处理带电粒子与中性粒子之间的碰撞. 得到了放电通道内离子与电子的分布以及离子流的运动,并且对出口外侧的能量分布进行了统计. 结果发现圆柱形阳极层等离子体加速器的磁场对电子有明显的约束作用,电子集中于阳极附近很小的区域内. 由于电磁场的特殊分布,离子流呈现出双峰式的分布. 离子能量范围从放电电压的20%到接近放电电压,平均能量在放电电压的40%—50%之间. 关键词： 质点网格方法 蒙特卡罗碰撞 数值模拟 阳极层霍尔等离子体加速器     

Self-consistent charge-up simulation for the microscopic feature of SiO2 layer in rf capacitive discharge

The charge-up simulation of the microscopic feature with SiO₂ layer was investigated in various conditions of rf capacitive discharge by using the three-dimensional (3-D) particle-in-cell (PIC) charge-up simulation coupled with the one-dimensional (1-D) particle-in-cell Monte Carlo collision (PIC-MCC) simulation of rf capacitive argon discharge. The result showed that the charge-up effect on the micro-trench was greatly influenced by the conditions of the gas pressure and the discharge voltage in rf capacitive discharge. Based on the analysis of the distributions of electrons and ions arriving at the substrate in various plasma conditions, the charge-up effect and its reduction mechanisms on the micro-trench of capacitive discharge were discussed. This article is expected to provide qualitative and quantitative insight for the understanding of charging and its reduction mechanism on many plasma processes performed by the rf capacitive discharge.