Elasticity Constants of a Two-Phase Tungsten Thin Film

The IET was used to determine the macroscopic elasticity constants of the multiphase coating. In order to determine the macroscopic elasticity constants of the film firstly, a critical assessment of Young’s modulus determination was done by comparing all the models proposed in the literature. The best model was identified and a study was performed to identify and quantify the most influent factors on the global uncertainty. Secondly, an enhanced formulation to determine the shear modulus of coating by IET was developed. The methodology was applied on a tungsten thin film deposited by DC magnetron sputtering.


Introduction
In order to improve the physicochemical and mechanical properties of a material, new surface treatment processes were developed.Among the most common processes we can mention: PVD, CVD and thermal spray.The elastic properties of thin films are different from those of bulk materials and depend on the technique of elaboration and its parameters.Moreover, a coating is generally polycrystalline, textured and composed of different phases.This makes it necessary to control the elasticity of the coated components.
The present paper proposes a methodology to determine the elasticity constants of a multiphase polycrystalline textured thin film.Firstly, the macroscopic elasticity constants of a coating will be determined by using IET.Therefore, a study of the models proposed in the literature to determine the Young's modulus of films will be done with a view to identify the best model.Secondly, an uncertainty study will be performed to identify the errors or uncertainty sources and evaluate their contributions.Then an enhanced analytical formulation to determine the shear modulus of thin films will be developed.It is valid for any ratio of coating to substrate thicknesses.Finally, an example of application on a tungsten thin film deposited by DC magnetron sputtering will be presented.

Impulse Excitation Technique
The principle of the measurement is to perform a two-steps frequency measurement, a first one for the substrate without deposition and a second one for the coated substrate.Then, knowing the elastic modulus of the substrate, the density of both coating and substrate and the resonance frequencies of the coated beam and the substrate, the Young's modulus of the coating can be determined through analytical expressions.

Analytical Models to Determine the Young's Modulus of a Coating
Lopez has developed a model for a rectangular beam, which does not take into account the shift of the neutral axis after deposition (Equation ( 1)) [1].
Pautrot and Mazot, in their model, take into account the shift of the neutral axis (Equation ( 2)) [2]: where: is the Young's modulus and the indices c and s mean respectively the coating and the substrate, Berry's model (Equation ( 6)) [3] is a simplified approach which can be obtained from the first order Taylor series expansion of either Pautrot's or Lopez's model.
The CLBT model is based on the Classical Laminated Beam Theory [4]: For a one layer coating (substrate + coating): A divergence between Berry's and Lopez's models has been noticed by Mazot [1] for a thicknesses ratio in the order of 0.1.Moreover, the difference between the models becomes larger for higher ratios.To analyze quantitatively the differences between the models, a comparison with a finite element model, taken as a reference, will be performed to assess the trueness of the method.Then, an uncertainty analysis will be done to compare the behavior of the different models with regard to their sensitivity to various uncertainty sources.

Identification of the Best Model to Measure Young's Modulus
In order to identify the best model, a parametric comparison of the analytical models with the developed finite element model (more details can be found in [5]) was done as a function of the thickness ratio , the Young's modulus ratio and the density ratio .Four different and ratios were chosen and was varied from 0 to 0.5.Figure 1 represents the evolution of the frequency ratio as a function of , and ratios.A good agreement between the Pautrot's model and the Finite Element Model (FEM) for any , and ratios can be seen.The divergence of the other models can be clearly seen and it may be due to the different assumptions on which is based each model.Based on these results, Pautrot's model is recommended to determine the Young's modulus of coatings.

Uncertainty Analysis
The uncertainty analysis was performed following the guidelines of the Guide to the expression of uncertainty in measurement [6].More explanation about the uncertainty calculation can be found in our previous work [5].The results of an example of application on an aluminum thin film deposited on a steel substrate are listed in Table 1.From Table 1, we can conclude that the main sources of uncertainty are the density and the thickness of the coating and the frequency before and after deposition.

Shear Modulus
The principle of the measurement is the same that used to determine the Young's modulus.The shear modulus can be determined using Gadaud's model (Equation ( 11)) [7]: where: is the first torsional resonance frequency of the whole bilayer beam.G, h, b and are respectively the shear modulus, thickness, width and density.
The indices c and s mean respectively the coating and the substrate.This model is valid only if the thickness of the coating is very small compared to that of the substrate.It does not take into account the shift of the neutral axis after deposition.We developed a new formulation considering the shift of the neutral axis; applicable with any shear moduli, thicknesses and density ratios.The proposed relation was established thanks to the application of the Hamilton's principle based on the minimization of the strain and kinetic energy [8] (Equation ( 12)): where: To check the reliability of the developed formulation, a comparison between Gadaud's model, the new enhanced formulation and the FEM was done (Figure 2).The comparison was performed with a combination of four different and ratios with variation of the ratio between 0 and 0.5.From Figure 2, a good agreement between the FEM and the new formula can be seen.

Determination of the Macroscopic Elasticity Constants
The Pautrot's model and the new enhanced formulation have been applied to determine the macroscopic elasticity constants of a tungsten thin film deposited on glass substrates by DC magnetron sputtering (Table 2).For more details about the deposition parameters, the reader can refer to our previous work [8].The resonance frequencies of the substrate and coated beam have been determined using the IET.The uncertainty analysis was performed following the methodology detailed in our previous work [5].The elasticity constants are listed in Table 2.The uncertainty on the substrate and coating elasticity constants are respectively presented in Tables 3 and 4.

Conclusions
In order to obtain the macroscopic elasticity moduli of a film, we have applied literature models after measurements by the IET.We have identified the best model and have done an uncertainty analysis to determine the sensitivity of the models to the different sources of errors.Then we have developed a new formulation to determine the shear modulus.
Tungsten film has been deposited by DC magnetron sputtering and the macroscopic elasticity constants of the coating have been determined by IET.

Table 1 .
Contribution of each uncertainty source (in GPa and %) on the coating Young's modulus uncertainty.

Table 2 .
Elasticity constants of the glass substrates and the as-deposited tungsten film.

Table 3 .
Uncertainties on the glass substrate elasticity constants obtained by IET.

Table 4 .
Uncertainties on the coating elasticity constants obtained by IET.