3.1.3 Simulating Oxide Growth using Volume Expansion
Silvaco, Inc., as mentioned in Section 3.1.1, offers a three-dimensional process simulation tool which includes
empirical and mechanical oxidation simulations. Their approach to simulating the mechanics of silicon oxidation includes the use
of the LS method . Instead of using unstructured meshes, their approach makes use of a fixed Cartesian mesh in
a LS environment, solving the problem of moving boundaries which arise when unstructured meshes are used.
The model includes four major steps which work together to generate a moving oxide interface:
- Oxidant diffusion through the oxide is modeled using the well-known diffusion equation
where is the diffusion coefficient, is the reaction rate, is the gas-phase mass-transfer coefficient, is the
equilibrium bulk concentration in the oxide, and is the normal to the corresponding interface. The diffusion equation is
obviously identical to the one presented in Section 2.3.1.
- Propagation of the oxide-silicon interface is solved after the diffusion equation for a time step .
To find the new position of the SiO-Si interface, the LS is solved for the distance function
where is the number of oxidant molecules incorporated into a unit of oxide and is the Si
coefficient. An up to third order Total Variation Diminishing (TVD) Runge-Katta scheme is used for discretization .
- The volume expansion resulting from the chemical reaction is calculated using the creep equation
where is a Cauchy stress tensor. Assuming a Maxwell visco-elastic fluid, the Cauchy stress tensor becomes
where is a velocity component.
The system of equations (3.13) and (3.14) make up the Stoke's equations with the boundary conditions
denotes the jump across a liquid-liquid interface, is the surface tension coefficient, and
is the surface curvature.
- Propagation of the interfaces in the deformation velocity field is the step during which all interfaces above the Si-SiO interface
are updated. This includes the SiO-Nitride interface, SiO-ambient interface, and Nitride-ambient interface. The interfaces
are updated using the LS representation
The main advantage of the LS method when compared to the FEM for solving visco-elastic problems is that all potential errors which arise
due to unstructured mesh irregularities are removed. The LS also makes it very easy to follow shape topologies which change
with time. The LS method natively handles complex surfaces, which can split or recombine during a simulation process.
L. Filipovic: Topography Simulation of Novel Processing Techniques