A mathematical model and simulation results of plasma enhanced chemical vapor deposition of silicon nitride films

We developed a mathematical model of Plasma Enhanced Chemical Vapor Deposition (PECVD) of silicon nitride thin films from SiH4-NH3-N2-Ar mixture, an important application in modern materials science. Our multiphysics model describes gas dynamics, chemical physics, plasma physics and electrodynamics. The PECVD technology is inherently multiscale, from macroscale processes in the chemical reactor to atomic-scale surface chemistry. Our macroscale model is based on Navier-Stokes equations for a transient laminar flow of a compressible chemically reacting gas mixture, together with the mass transfer and energy balance equations, Poisson equation for electric potential, electrons and ions balance equations. The chemical kinetics model includes 24 species and 58 reactions: 37 in the gas phase and 21 on the surface. A deposition model consists of three stages: adsorption to the surface, diffusion along the surface and embedding of products into the substrate. A new model has been validated on experimental results obtained with the "Plasmalab System 100" reactor. We present the mathematical model and simulation results investigating the influence of flow rate and source gas proportion on silicon nitride film growth rate and chemical composition.