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 [203]. 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
(81)

(82)

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
(83)

where is the number of oxidant molecules incorporated into a unit of oxide and is the Si SiO expansion coefficient. An up to third order Total Variation Diminishing (TVD) Runge-Katta scheme is used for discretization [66]. - 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

(86)

where 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
(87)

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