Numerical analysis of a mixture of Ar/NH3 microwave plasma chemical vapor deposition reactor

A two-dimensional fluid model has been used to investigate the properties of plasma in Ar/NH3 microwave electron cyclotron resonance discharge at low pressure. The electromagnetic field model solved by the three-dimensional Simpson method is coupled to a fluid plasma model. The finite difference method was employed to discrete the governing equations. 40 species (neutrals, radicals, ions, and electrons) are consisted in the model. In total, 75 electron-neutral, 43 electron-ion, 167 neutral-neutral, 129 ion-neutral, 28 ion-ion, and 90 3-body reactions are used in the model. According to the simulation, the distribution of the densities of the considered plasma species has been showed and the mechanisms of their variations have been discussed. It is found that the main neutrals (Ar*, Ar**, NH3*, NH, H2, NH2, H, and N2) are present at high densities in Ar/NH3 microwave electron cyclotron resonance discharge when the mixing ratio of Ar/NH3 is 1:1 at 20 Pa. The density of NH is more than that of NH2 atom. And NH3+ are the most important ammonia ions. But the uniformity of the space distribution of NH3+ is lower than the other ammonia ions.