The Berg model is used for modelling the process of reactive sputtering [
9]. The results of the modelling can be used in a variety of ways, from predicting processing behavior to scaling up the deposition and automated control of the manufacturing process [
10]. Hegedüs et al. used Berg modelling for the calculations of the hydrogen incorporation into radio frequency sputtered amorphous silicon films [
11]. In this work, the Berg model was applied in order to fit the deposition temperature and the partial pressure of oxygen during the reactive sputtering of SiO
xN
y layers. In the Berg model, the target material and the gases involved in the sputter deposition are characterized by three quantities, namely the metal sputtering yield Y
m, the compound sputtering yield Y
c, and the sticking coefficient α. The actual values of the material constants in each sample were determined from a fit at the beginning of the variable composition layers, where the partial pressure of oxygen was 8 × 10
−3 Pa and the deposition temperature was 310 K. The point-by-point changes in both oxygen partial pressure and deposition temperature were fitted using the measured compositional variation of the sample, while knowing the material constants. Although, the Berg model can handle two different reactive gases, we present here a different approach. Our model takes into account a single reactive gas, oxygen, and a composite target, silicon nitride, instead of silicon, given the effect of nitrogen gas applied in the chamber. This approach is based on the two facts that the Berg model considers only the erosion of the top atomic layer of the target and that enough nitrogen gas has been introduced to form silicon nitride over the entire surface of the target. Furthermore, the fact that nitrogen accounts for 90% of the total pressure in the deposition chamber and argon for only 10% also confirms that nitrogen and oxygen can be considered a sputtering and reactive gas, respectively.
The section is structured into two subsections; the equations of the Berg model used for the calculations have been briefly summarized and the fitting results for the material constants, the deposition temperature, and the oxygen partial pressure have been presented, respectively.
4.1. Berg’s Model Approach
The reactive gas flow leads to a uniform bombardment of the neutral reactive gas molecules to the inner surface of the deposition chamber. The
F flux of the molecules in molecules/unit area and unit time can be expressed as follows:
where
p,
k,
T, and m refer to the partial pressure of the reactive gas, the Boltzmann constant, the absolute temperature in the processing chamber, and the mass of the reactive gas molecule, respectively.
Since reactions take place between the elemental target atoms and the reactive gas molecules, compound molecules will cover a certain,
fraction of the target surface. This phenomenon is usually called “target poisoning” in the literature [
12].
At the same time, the entire target surface (including the
compound fraction and the
elemental metal fraction) will be continuously sputter eroded. The rate of target poisoning (compound formation) and the rate of compound sputtering are identical in the equilibrium state, which is expressed by the following equation:
where
J,
q,
Yc, and
refer to the ion current density over the target surface, the elementary electronic charge, the sputtering yield of compound molecules, and the probability of the reaction between a reactive gas molecule and a target atom due to collision, respectively.
From Equations (1) and (2),
can be expressed as follows:
As a result of the above, the sputter-eroded material from the target surface can be devided into two groups: (i) compound material sputtered from the
fraction, and non-reacted target atoms sputtered from the
fraction of the target.
Fc and
Fm refer to the total number of eroded molecules or atoms from the
and the
fraction of the target, respectively, and can be expressed by the following equations:
where
and
refer to the total area of the target surface and the metallic sputtering yield, respectively.
In the Berg model’s approach, the term “collecting surface” is used to refer to the Ac total surface area (including the substrate) that receives sputtered material during the deposition. All eroded materials from the target are assumed to be uniformly deposited on the entire collecting surface.
The material deposition on the collecting surface is considered to be the result of five different material fluxes, as is presented in
Figure 6.
The Fc compound material flux eroded from the part of the target leads to compound formation on the compound part and metallic part of the collecting surface. Similarly, the Fm flux of eroded metallic material promotes metallic deposition on the and parts of the collecting surface.
In addition to the target erosion, compound formation by reactions between the reactive gas molecules, and the non-reacted atoms on the collecting surface will also contribute to the layer growth on the collecting surface. These reactions will consume
Qc amount of reactive gas molecules, which can be expressed by as follows based on
Figure 6:
where
refers to the compound fraction of the collecting surface.
From the material fluxes presented in
Figure 6,
Fc does not change the stoichiometry of the layer, and thus the value of
, since this material flux refers to compound material erosion from the target to the
part of the collecting surface which is already in a compound state. Similarly,
Fm will also not change the value of
due to the metallic material erosion to the
metallic part of the collecting surface. On the other hand,
Fm will decrease, while
Fc and Q
c will increase the value of
Assuming steady state at the collecting surface, contributions to the increase and decrease of
should be identical, which can be mathematically expressed by:
In Equation (8), Qc is multiplied by two based on the fact that one reactive gas molecule contains two reactive gas atoms; therefore, one reactive gas molecule will contribute with two compound molecules to the increase of
From Equations (1) and (2) and Equations (5)–(8),
can be solved as a function of the
p partial pressure of the reactive gas and the
T deposition temperature:
From Equations (3) and (9)
can be expressed in the following form:
The
J ion current density over the target was calculated based on the following equation [
13]:
where
P and
U refer to the measured power and effective voltage, respectively.
In the following subsections, we present the results of two fitting processes, aiming to fit (i) the
Ym,
Yc, and
material constants and (ii) the O
2 partial pressure and the
T temperature based on the Berg model approach. In the case of all constants fitting, the
experimental value of
and the
modelled value of
were compared at certain points of the sample.
was determined from EDS measurements based on the following equation:
where [
O] and [
N] refer to the concentration of oxygen and nitrogen, respectively.
was calculated from Equation (10), where At and Ac were considered to be 44 cm2 and 80 cm2, respectively.
During the fit, the mean squared error between
and
was minimized by varying the free parameters of
. Since two different fitting processes were performed, we describe the used measurement points and the free parameters for both fittings separately in
Section 4.2 and
Section 4.3.
4.2. Fitting of the Material Constants
At the beginning of the four SiO
xN
y combinatorial layer deposition, the values of
T (310 K) and p (8 × 10
−3 Pa) are known. Therefore, the first measured points of the four samples were used to fit the
Ym,
Yc, and
material constants. The value of F was calculated from Equation (1), knowing T and p. The J ion current density was calculated based on Equation (11) from the measured values of the RF power and the DC potential presented in
Section 2.1. Knowing
p,
T,
F,
J,
At, and
Ac, the only unknown parameters in Equation (10) is
Ym,
Yc, and
, which were considered to be free parameters. During the fit, these free parameters were varied in order to minimize the mean squared error between the four
data and the corresponding
values.
Figure 7 presents
data and the corresponding
values obtained from the fit results.
The fitting results of the
Ym,
Yc, and
material constants are reported in
Table 1.
The fitted value of
Ym (0.7) is very close to the value reported in the literature [
14] for nitrogen (~0.5–0.6). The resulted value of α is also in the range of its typical reported value (0.5–1) [
10]. However, we are not aware of any reported data in the literature for
Yc for nitrogen in case of SiO
xN
y compound material, and the fact that the resulting
Yc value is lower than the
Ym is in good alignment with the known correlation that
Ym is always higher than
Yc [
15].
The good alignment between the
and
values as well as the agreement between the obtained material constants and their values reported in the literature prove that the Berg model equations of
Section 4.1 are suitable for describing the reactive sputtering process of oxygen-rich SiO
xN
y layers.
4.3. Fitting of Temperature and Oxygen Partial Pressure
The variations of p and T during the deposition of SiO
xN
y layers sputtered with 1.62 and 1.95 kV were determined by the second fitting process. However it is known that the values of
Ym, Yc, and
could vary slightly under different gas processing conditions; in the application of the Berg model approach they are usually considered to be constant [
10], characterizing the affected materials and sputtering equipment independently from the actual process parameters. Based on this assumption, the values of
Ym, Yc, and
were considered to be the values reported in
Table 1 throughout the entire deposition process of all samples. Similar to the first fitting process, the value of J at each measurement point of the samples was calculated based on Equation (11) from the measured RF power and the DC potential presented in
Section 2.1. Knowing
Ym, Yc, and
and
At and
Ac, as well as
J at a certain point of the sample,
, referring to the stoichiometry of the layer at the affected point, can be calculated as a function of
p and
T based on Equations (1) and (10). Varying
p and
T as free parameters, the difference between
and
at the affected point can be minimized. Based on this method,
p and
T were fitted individually in each measured points of the samples sputtered with 1.62 and 1.95 kV.
Figure 8 represents a very good agreement between the fitted
and
data in case of all points of the samples, proving that the introduced Berg model equations are able to correctly follow the lateral compositional variation of the samples.
Figure 9 shows the variation of the fitted oxygen partial pressure and the deposition temperature as a function of sputter time. Essentially, the shape of the oxygen partial pressure functions (
Figure 9a) follows the characteristic behavior of the peristaltic pump [
6] applied to the oxygen injection during the deposition. For both samples, the variation of the fitted O
2 partial pressure is in a good agreement with the lateral compositional variation measured by EDS (
Figure 2).
Figure 9b represents the fitted T values of the deposition process carried out at 1.62 and 1.95 kV DC potential. At the beginning of O
2 injection, T increases regardless of the applied DC potential and reaches its maximum after about 45 min of sputtering. Such an increase in temperature could be attributed to the exotherm heat of oxidation [
16,
17]. It is also supported by the fact that after ~45 min of sputtering,
T starts to decrease in agreement with the reduced O
2 partial pressure shown by
Figure 9a. The maximum of
T for the process carried out with 1.62 kV was ~345 K, while increased DC potential, and thus increased RF power, resulted in a higher maximum of the temperature (~350 K), which suggests that the effect of the oxidation on the temperature defined in Equation (1) was more pronounced when a higher DC potential was applied.