6.5 Model Calibration

A look to the model (see Section 3.2.1) shows that there are three available parameters, namely the diffusion coefficient , the maximal possible strength of the spatial sink , and the oxidant concentration at surfaces which have contact to the oxidizing atmosphere. As displayed in (3.9) and (7.1) the diffusion coefficient has a physical background. It is temperature and stress dependent and its real physical value can be determined correctly. Therefore to use for calibration is not appropriate.

The next parameter has more mathematical and modeling origin, but it is also not an optimal paramenter for calibration. At first the thickness of the reaction layer changes with , because it is inversely proportional to (see Section 3.2.4). This can be a problem for small values which lead to thick reaction layers.

For a better understanding of the second trouble with large values the following is worth mentioning: Simulations have shown that with regard to the finite elements for the value of the following choice is reasonable

Here, is the average diameter of the finite elements in the used mesh, and is a constant value independent of the mesh fineness.

Due to the mesh dependence of its variation is limited. In the experiments it was found out that the value of can not be increased arbitrarily. For larger values than suggested in (6.4) the numerical formulation becomes instable. In contrast to this, small values are not a problem. Therefore, is not a suitable parameter for the model calibration of a potentially because of thick reaction layer (small value) and numerical instability (large value).

After excluding two of the three parameters, the last parameter which is the surface oxidant concentration is investigated. On surfaces which have contact with the oxidizing atmosphere the oxidant concentration is used as a Dirichlet boundary condition. The key idea is to modify in order to calibrate the oxide thickness of the simulated oxidation process over time for different oxidation conditions. From the physical aspect a higher surface oxidant concentration means that a larger number of oxidants diffuse to the Si/SiO-interface and react with silicon, which results in a faster oxidation rate.

It was found with experiments that the best results are obtained if consists of a constant part and an -dependent part so that the effective surface concentration can be written as a function of

Here is the standard oxidant concentration in the gas atmosphere as used in the Deal-Grove model. , , and are the calibration parameters. Because the value of is changed during the oxidation process, the value of is also changing with time. For the onward process goes toward 0 on the surface and so the second term for disappears.

As example for the above described calibration concept a (111) oriented and 0.4m height silicon block is wet oxidized and the oxide thickness over time for different temperatures is calibrated. The bottom surface is fixed, the lateral surfaces can only move vertically and on the upper surface a free mechanical boundary condition is applied. Only the upper surface of the body has contact with the oxidizing atmosphere. The oxide thickness is measured between the upper surface and the -level of 0.5.

In the calibration process the values of the parameters , , and are determined with the help of the in-house tool SIESTA (Simulation Environment for Semiconductor Technology Analysis) [108], so that the thickness values of the simulated oxide layers agree with the calculated physical reference values up to approximately 500nm at any time for a temperature range of 900-1100C. The temperature dependent diffusion coefficient is calculated as explained in (3.13). The other two model parameter and are kept constant over the whole temperature range.

It was found that in case of wet oxidation the value of the parameter in (6.5) can be hold constant for the temperature range of T=900-1100C. Furthermore, the experiments show that the parameter can be brought to a linear and the parameter can be brought to a parabolic dependence on temperature, which is described by

The expressions for and were found empirically. Their values over temperature are plotted in Fig. 6.11. and do not have a physical background, they are pure fitting parameters.

In case of wet oxidation the three Figs. 6.12-6.14 show vicegerent for all other temperatures that the formula for with its parameters , , and leads to an excellent agreement between the calculated reference curves and the measured simulation curves. The oxide thickness values for the reference curves are calculated with the Deal-coefficients [51]. The coefficients , , and can be also found without problems for other process conditions (e.g. dry oxidation), and so the calibration with always works well.

Ch. Hollauer: Modeling of Thermal Oxidation and Stress Effects