A Strategy for Extracting Full Material Coefficients of AlN Thin Film Based on Resonance Method

AlN thin film is widely used in piezoelectric MEMS devices, and the accurate characterizations of its material coefficients are critical for the optimization of the AlN thin film process and the design of AlN thin-film-based devices. However, it is difficult to extract the material coefficients of AlN in the form of thin film. This paper reports a strategy for systematically extracting full elastic coefficients, piezoelectric coefficients and dielectric constants of c-axis-oriented AlN thin film based on the resonance method outlined in IEEE Standard on Piezoelectricity Std 176-1987. In this strategy, four self-suspended resonators with length thickness extension (LTE), thickness extension (TE), radial extension (RE), lateral electric field excited thickness shear (LEF-TS) modes together with a lamb wave resonator (LWR) are specifically adopted, and the material coefficients of AlN thin film are extracted by measuring the impedance spectra of these resonators. In addition, the effects of the pad and electrodes on the resonators were systematically studied, and the corresponding procedures to eliminate their influences on the extraction accuracy of material coefficients were proposed. Finally, a complete extraction process based on the above strategy was established. The simulation results show that the strategy can achieve high accuracy for AlN thin film with different thicknesses and electrode configurations, and it can also be applied to other materials belonging to the 6 mm piezoelectric crystal class such as ZnO, ScAlN, etc.


Introduction
AlN thin film has attracted extensive attention because of its superiorly material properties, such as fairly good piezoelectric properties, high sound velocity and stability, low deposition temperature (<400 • C) through RF or pulsed DC magnetron sputtering process and compatibility with CMOS process. Benefiting from these properties, the AlN thin film has become a promising material for piezoelectric MEMS devices such as sensors [1,2], actuators [3,4], resonators [5,6], filters [7,8], energy harvester [9] and many other fields. For the design and performance prediction of these devices, accurate material coefficients of AlN thin film such as elastic coefficients, piezoelectric coefficients and dielectric constants are highly required. Moreover, an accurate and complete characterization method for these material coefficients is necessary for the development and optimization of new thin film processes [10,11] since the characterization of material properties will contribute to the assessment of the thin film process.
Until now, many methods have been reported to predict the material coefficients of AlN. First-principles method based on density functional theory (DFT) can accurately calculate the material constants of AlN single crystals [12]. While the actual material coefficients of the AlN thin film vary greatly with its thickness, substrate and deposition process because the AlN thin film will first start growing with a precision of alignment as defined by the underlying substrate, which results in a disordered layer at the bottom of the thin film [13]. Hence, accurate material coefficients of the AlN thin film need to be extracted experimentally. Tonisch et al. measured the piezoelectric coefficients of the AlN thin film using the piezoresponse force microscopy and interferometric technique [14], which is called the direct measurement method. While this method can only be used to measure the effective piezoelectric coefficients d 33, f and d 31, f due to the clamping of the substrate. Tsubouch et al. determined ten independent material coefficients of the AlN thin film at once by fitting the theoretical Rayleigh wave velocity v-f curve to the measured one [15], but this method has a certain degree of uncertainty because it is susceptible to model errors, and sensitive to the selection of fitted regions and constraints on fitting results during curve fitting. Ohashi et al. calculated all material constants of AlN single crystal from the basic acoustic coefficients evaluated using ultrasonic microspectroscopy technology [16]. However, this method requires a certain thickness of the samples for measurement; therefore, it cannot be applied to thin film materials. Konno et al. determined full material constants of 40% ScAlN thin film by measuring resonance frequencies and electromechanical coupling coefficients of four MEMS resonators and fitting theoretical v-f curves of a lamb wave resonator to the measured ones [17], while this method also has a certain degree of uncertainty because half of the material coefficients are obtained by curve fitting. Hence, there is currently a lack of systematic and accurate characterization method for the material coefficients of the AlN thin film.
The IEEE Standard on Piezoelectricity Std 176-1987 outlines a resonant method [18], which can systematically calculate the material coefficients of the piezoelectric materials by the resonant frequencies f r , anti-resonance frequencies f a and electromechanical coupling coefficients k 2 measured from the impedance spectra of different resonators. This method provides fairly high accuracy and is easy to implement as long as samples have specific shapes according to this standard, which makes it widely used for coefficient extraction of piezoelectric materials [17,19,20]. However, when it is applied to the AlN thin film, many factors will seriously affect extraction accuracy, mainly including: (1) the clamping of substrate will severely limit the motion of magnetic domain walls inside crystals, and influences the measured material coefficients [21,22]; (2) the actual measurement results contain the parasitic parameters of the pad, which will cause serious distortion of the measured impedance spectra [23]; and (3) the mass loading effect of the electrodes will shift the frequencies and electromechanical coupling coefficients of the resonators, and the extent of which depends on the density and acoustic impedance of the electrode material [24,25].
Thus, this paper proposes a strategy for systematically extracting full elastic coefficients, piezoelectric coefficients and dielectric constants of the AlN thin film, where the influence of the above factors on extraction accuracy of material coefficients are studied and eliminated by the corresponding procedures. Furthermore, a complete extraction process is established in which four bulk acoustic wave (BAW) resonators with length thickness extension (LTE), thickness extension (TE), radial extension (RE), lateral electric field excited thickness shear (LEF-TS) modes together with a lamb wave resonator (LWR) are adopted with the light of [17]. The strategy has high accuracy for the AlN thin film with different thickness and electrode configuration, and it could be applied to other materials of 6 mm crystal class.

The Material Coefficients of AlN
Besides the density ρ, AlN with a wurtzite structure has ten independent material coefficients, including five elastic coefficients c E 11 , c E 12 , c E 13 , c E 33 , c E 44 , three piezoelectric coefficients e 31 , e 33 , e 15 and two dielectric constants ε T 11 , ε T 33 , which can be expressed as the following matrix: Piezoelectric coefficient matrix: Dielectric constant matrix: where the superscript E and T represent the measurement conditions of a constant electric field and a constant stress, respectively.

Extraction Method
The density ρ can be determined by the standard water displacement method described in [19]. The weight of the sample is first measured by a precision balance. Then, the sample is suspended in the water by a thin wire, and the weight of the water at this time is the sum of its actual weight and the buoyancy. Thus, the weight of the water before and after adding the suspending sample can be measured by the precision balance, and the difference in the two cases give the buoyancy of the water, which is divided by the density of the water to obtain the volume of the sample according to Archimedes' principle. Based on this method, the weight and volume of the AlN thin film can be obtained from the difference in ones of the substrate sample before and after thin film deposition, and the ratio of weight to volume gives the density.
Based on the resonance method outlined in the IEEE Standard on Piezoelectricity Std 176-1987 [18], while considering the little thickness of the thin film and the fabrication of the resonators, five c-axis-oriented AlN thin-film-based resonators with different modes were proposed to extract a complete set of material coefficients of the AlN thin film: four BAW resonators with LTE, TE, RE and LEF-TS mode and a LWR as shown in Figure 1. The geometries of these resonators are set to specific shapes to reduce spurious modes, and the ten material coefficients described in the above matrices are obtained by the following method, where the effects of the substrate, pad and electrodes are not considered temporarily, and they will be studied and eliminated in Section 2.3.
The dielectric constant ε T 33 can be simply converted from the plate capacitance value C T p , which is best measured at a frequency well below the lowest resonant frequency (typically 1 kHz), of any BAW resonator by Equation (4), where the RE mode resonator is recommended because its large plate area A can reduce the measurement error to some extent. The dielectric constants ε T 11 can be calculated from the off-resonance capacitance C T LWR of LWR by Equation (5) [17].
where, ε 0 , A, t, N, W and λ = 2(w 1 + g 1 ) are the vacuum dielectric constant, the plate area and thickness of the AlN thin film, the number, aperture and wavelength of the IDT pairs, respectively. It should be noted that A includes the surface area of an anchor because the anchor on one side also covers the electrodes for the practical reason: the bottom electrode is patterned by the anisotropic etching, so the electrode at the bottom of the anchor is difficult to be etched due to the shadowing of AlN.
Micromachines 2022, 13, 4 of 13 where, , , , , and = 2 + are the vacuum dielectric constant, the plate area and thickness of the AlN thin film, the number, aperture and wavelength of the IDT pairs, respectively. It should be noted that includes the surface area of an anchor because the anchor on one side also covers the electrodes for the practical reason: the bottom electrode is patterned by the anisotropic etching, so the electrode at the bottom of the anchor is difficult to be etched due to the shadowing of AlN.
The eight elastic and piezoelectric coefficients ( , , , , , , , ) are calculated by the measured resonant frequencies and anti-resonant frequencies extracted from the impedance spectra of the five resonators. In detail, by using the LTE mode resonator, the elastic coefficient is calculated from by Equation (6), and the piezoelectric coefficient is extracted from the calculated electromechanical coupling coefficient by Equation (7).
Similarly, the TE mode resonator is used to calculate the elastic coefficient from and by Equation (8), and the piezoelectric coefficient can be extracted by Equation (9), while the dielectric constant is temporarily unknown, so it will be obtained below by combining it with the equations of LWR. It is worth mentioning that this resonator is designed to be elliptical to reduce the spurious mode.
(9) The eight elastic and piezoelectric coefficients (c E 11 , c E 12 , c E 13 , c E 33 , c E 44 , e 31 , e 33 , e 15 ) are calculated by the measured resonant frequencies f r and anti-resonant frequencies f a extracted from the impedance spectra of the five resonators.
In detail, by using the LTE mode resonator, the elastic coefficient s E 11 is calculated from f r by Equation (6), and the piezoelectric coefficient d 31 is extracted from the calculated electromechanical coupling coefficient k 2 31 by Equation (7).
Similarly, the TE mode resonator is used to calculate the elastic coefficient c E 33 from f a and k 2 t by Equation (8), and the piezoelectric coefficient e 33 can be extracted by Equation (9), while the dielectric constant ε S 33 is temporarily unknown, so it will be obtained below by combining it with the equations of LWR. It is worth mentioning that this resonator is designed to be elliptical to reduce the spurious mode.
As mentioned in [26], Poisson's ratio σ depends on the ratio r σ = f I I r / f I r of the first resonant frequency f I r and second resonant frequency f I I r of the RE mode. Thus, the elastic coefficient s E 12 can be extracted by Equation (10), where σ is obtained by looking up the table in [26] from the measured r σ . In addition, Equation (11) is the calculation equation of f I r and f I I r , which will be used in Section 2.3.3.
where η I r is a coefficient related to Poisson's ratio σ, which can be obtained by looking up the table in [26].
Next, c E 44 and e 15 are extracted using the LEF-TS mode resonator instead of the common TS mode resonator since the latter requires electrodes to be added on both sides of the thin film, which is difficult to fabricate. While in LEF-TS mode resonator, both electrodes are on the top surface of the AlN thin film, and they will generate a lateral electric field to excite the TS mode [27]. In the same way, the elastic coefficient c E 44 is calculated from f a and k 2 15 by Equation (12), and the piezoelectric coefficient e 15 is extracted by Equation (13).
where the above equations are all derived from the conversion equations of the material coefficients including c E ij = s E ji −1 , d mj = e mi s E ij and ε T mn − ε S mn = d mj e nj in [18]. In other word, s E 33 is the last coefficient that needs to be extracted, and it is originally obtained through a length extension (LE) mode resonator [20]. However, such resonator is not available because its thickness (generally a few µm) should be more than 6 times other sizes to correctly excite the mode, which results in other sizes about 1 µm or less. It will pose great challenges to the fabrication process. Instead, the LWR is used in this method, and s E 33 is extracted by the measured velocity v S 0 ,m (= f r λ) of S 0 mode lamb wave. The detailed numerical extraction method is as follows: Firstly, the s E 33 value in [15] or other reported values is adopted as the initial one. Then, the material coefficients c E 11 , c E 12 , c E 13 , e 33 and e 31 are converted by Equations (14)- (18), and the dielectric constants ε S 11 and ε S 33 can be obtained by Equations (19) and (20). In this way, all material coefficients of the AlN thin film have been calculated, and they are used to calculate the theoretical velocity v S 0 ,t according to the similar process in [28], where the model shown in Figure 2 is adopted, and eight boundary conditions including four electrical boundary conditions at and at x 3 = −t, ε 0 E I I I 3 = D I I 3 , E I I I 1 = E I I 1 , and four elastic boundary conditions at both x 3 = 0 and x 3 = −t, T I I 3 = 0 and T I I 5 = 0 are used. The calculated v S 0 ,t will be compared with the measured v S 0 ,m , and finally a suitable s E 33 can be found by adjusting its value so that v S 0 ,t is equal to v S 0 ,m . Using this s E 33 , the material coefficients c E 11 , c E 12 , c E 13 , e 33 and e 31 can be extracted by Equations (14)- (18).
Micromachines 2022, 13, 6 of 13 AlN thin film have been calculated, and they are used to calculate the theoretical velocity , according to the similar process in [28], where the model shown in Figure 2 is adopted, and eight boundary conditions including four electrical boundary conditions at = 0, = , = , and at = − , = , = , and four elastic boundary conditions at both = 0 and = − , = 0 and = 0 are used. The calculated , will be compared with the measured , , and finally a suitable can be found by adjusting its value so that , is equal to , . Using this , the material coefficients , , , and can be extracted by Equations (14)- (18). So far, the 11 independent material coefficients of the AlN thin film have been completely extracted.

Self-Suspended Resonators
The small sizes of the resonators based on the AlN thin film determine that it must be supported by a substrate, while for ones based on the AlN single crystal, only electrodes are added to the bulk piezoelectric material [20]. The clamping of the substrate will greatly affect the extraction accuracy of material coefficients [21,22]. To solve this problem, as shown in Figure 1, the resonators are suspended in the air, and the anchors on both sides are used for supporting and introducing electrodes. Since the connection points of the anchors and the resonators are set at the minimum displacement of the resonance, the anchors have little effect on the frequencies, which is confirmed in FEM simulation. It is worth mentioning that the sizes of these resonators are best to be designed carefully to reduce the spurious modes, which will help to improve the extraction accuracy of the material coefficients.

De-Embedding of the Pad
In actual measurements, a measurement result from the probe includes parasitic parameters associated with the pad [23]. In order to measure the accurate responses of the resonators, these parasitic effects must be removed, i.e., de-embedding the pad. The following method is used in this paper to complete the de-embedding of the pad: A simple pad model composed of an equivalent resistor and an equivalent capacitor was established. According to transmission line theory, this resistor and capacitor can be obtained from the measured reflection coefficient Γp of a pad-only device (without resonator) through Equation (21), then the S parameters including , , and of So far, the 11 independent material coefficients of the AlN thin film have been completely extracted.

Self-Suspended Resonators
The small sizes of the resonators based on the AlN thin film determine that it must be supported by a substrate, while for ones based on the AlN single crystal, only electrodes are added to the bulk piezoelectric material [20]. The clamping of the substrate will greatly affect the extraction accuracy of material coefficients [21,22]. To solve this problem, as shown in Figure 1, the resonators are suspended in the air, and the anchors on both sides are used for supporting and introducing electrodes. Since the connection points of the anchors and the resonators are set at the minimum displacement of the resonance, the anchors have little effect on the frequencies, which is confirmed in FEM simulation. It is worth mentioning that the sizes of these resonators are best to be designed carefully to reduce the spurious modes, which will help to improve the extraction accuracy of the material coefficients.

De-Embedding of the Pad
In actual measurements, a measurement result from the probe includes parasitic parameters associated with the pad [23]. In order to measure the accurate responses of the resonators, these parasitic effects must be removed, i.e., de-embedding the pad. The following method is used in this paper to complete the de-embedding of the pad: A simple pad model composed of an equivalent resistor R and an equivalent capacitor C was established. According to transmission line theory, this resistor R and capacitor C can be obtained from the measured reflection coefficient Γ p of a pad-only device (without resonator) through Equation (21), then the S parameters including S 11 , S 12 , S 21 and S 22 of the pad can be calculated according to [29]. Thus, the reflection coefficient Γ r of the resonators without the pad can be obtained from the measured reflection coefficient Γ m including the effect of the pad by Equation (22). Finally, the accurate impedance spectra of the resonators can be converted from Γ r , and the de-embedding of the pad is completed.
where Z and Z 0 are the impedance of the pad and characteristic impedance, respectively.

Mass Loading Effect of Electrodes
The mass of the electrodes added to the piezoelectric material will seriously shift the frequency of the resonator [24,25], which is called the mass loading effect. This effect cannot be ignored for the AlN thin-film-based resonators because the thickness of the electrodes is comparable to the piezoelectric layer. Based on the principle of piezoelectric mass sensor [30], the frequency without electrodes ( f 0 ) and the frequency with electrodes ( f ) have the relationship as shown in Equation (23).
where m p and m e are the mass of the piezoelectric layer and electrodes, respectively. n is a constant related to the resonance mode and material coefficients of AlN because f / f 0 is basically linear with the mass of the electrodes [24]. It can be seen that there is a constant relationship γ between f 0 and f for a resonator with a certain size and mode. And the γ can be obtained by FEM simulation if the material coefficients are known. Thus, we can perform FEM simulation of these resonators with assumed material coefficients in COMSOL Multiphysics to obtain these constant scale factors γ r and γ a for each mode by Equation (24).
where f r0 and f a0 are the frequencies calculated from the assumed material coefficients by Equations (6)-(13), while f r0 = v S 0 /λ in LWR is calculated through the similar process in [28]. f r and f a are the frequencies obtained from the FEM simulation with the same material coefficients, which correspond to the actual measured frequencies. In the actual measurement, the mass loading effect of the electrodes can be eliminated by multiplying γ r and γ a with the measured frequencies. However, the assumed material coefficients used in simulation is not necessarily consistent with the actual values, which will lead to inaccurate γ and thus inaccurate extracted material coefficients. To solve this problem, an iterative approach is used, which will be explained in detail in Section 2.4.

Extraction Process
Based on the above analysis, five self-suspended resonators with different modes and a pad-only device will be adopted, and the 11 material coefficients of the AlN thin film can be extracted by the following steps as shown in Figure 3: (a) Firstly, the density ρ is determined by measuring the thin film sample based on the standard water displacement method; (b) The reflection coefficient Γ p of the pad-only device is then measured to calculted its S parameters by Equation (21); (c) Thus, the accurate impedance spectra of the five resonators can be measured after de-embedding the pad based on Equation (22); (d) The dielectric constants ε T 11 and ε T 33 can be calculated from the off-resonance capacitance of the LWR and the plate capacitance of BAW resonator by Equations (4) and (5); (e) Afterward, eight elastic and piezoelectric coefficients will be extracted by the following iterative approach. The FEM simulation of the five resonators are first performed with a set of initial material coefficients, which could be reported coefficients or other reasonable values, and the scale factors γ r , γ a of different modes can be obtained; (f) Using these scale factors, the frequencies without electrodes f r0 = γ r f r and f a0 = γ a f a can be calculated; (g) Next, the material coefficients s E 11 , s E 12 , c E 33 , c E 44 , d 31 and e 15 can be calculated directly from these frequencies by Equations (6)-(13); (h) Additionally, the material coefficient s E 33 can be obtained from the measured v S 0 by the numerical method; (i) Finally, the remaining material coefficients c E 11 , c E 12 , c E 13 , e 33 and e 31 can be calculated by Equations (14)- (18); (j) Since the initial material coefficients used may be different from the actual values, there will be a certain error in the extracted material coefficients. In order to improve the extraction accuracy, we can reestablish the FEM simulation with the extracted material coefficients and repeat steps (e-i) for iteration until the change in the extracted material coefficients between two iterations is less than a certain threshold, for example 1%, which depends on actual demand. of initial material coefficients, which could be reported coefficients or other reasonable values, and the scale factors , of different modes can be obtained; (f) Using these scale factors, the frequencies without electrodes = and = can be calculated; (g) Next, the material coefficients , , , , and can be calculated directly from these frequencies by Equations (6)-(13); (h) Additionally, the material coefficient can be obtained from the measured by the numerical method; (i) Finally, the remaining material coefficients , , , and can be calculated by Equations (14)-(18); (j) Since the initial material coefficients used may be different from the actual values, there will be a certain error in the extracted material coefficients. In order to improve the extraction accuracy, we can re-establish the FEM simulation with the extracted material coefficients and repeat steps (e-i) for iteration until the change in the extracted material coefficients between two iterations is less than a certain threshold, for example 1%, which depends on actual demand.

Results and Discussion
The above extraction strategy will be verified by the FEM simulation below. The sizes of the five resonators shown in Table 1 were designed with a typical 1, 2 μm-thick AlN thin film. Using the material coefficients of the AlN thin film in [15], the FEM simulations of these resonators were performed in COMSOL Multiphysics with 50 and 100 nm-thick Mo as the electrode material, where a rectangular pad of 100 × 100 μm and a routing of 20 × 30 μm based on 10 μm-thick Mo were used in all resonators. In the simulation, the "piezoelectric effect" physics interface was selected, and the 3D model of each resonator is established according to the sizes in Table 1, where the "fixed constraint" was added on the anchors to simulate the effect of the substrate on the resonator. The mesh is built by "swept", and its element size was set to the software predefined "semiconductor-normal", then the impedance spectrum was obtained by the analysis of the "frequency domain".

Results and Discussion
The above extraction strategy will be verified by the FEM simulation below. The sizes of the five resonators shown in Table 1 were designed with a typical 1, 2 µm-thick AlN thin film. Using the material coefficients of the AlN thin film in [15], the FEM simulations of these resonators were performed in COMSOL Multiphysics with 50 and 100 nm-thick Mo as the electrode material, where a rectangular pad of 100 × 100 µm 2 and a routing of 20 × 30 µm 2 based on 10 µm-thick Mo were used in all resonators. In the simulation, the "piezoelectric effect" physics interface was selected, and the 3D model of each resonator is established according to the sizes in Table 1, where the "fixed constraint" was added on the anchors to simulate the effect of the substrate on the resonator. The mesh is built by "swept", and its element size was set to the software predefined "semiconductor-normal", then the impedance spectrum was obtained by the analysis of the "frequency domain". Table 1. The sizes of the resonators used in the simulation, where the symbols correspond to those in Figure 1.

Mode Resonator (µm) Anchor (µm)
According to the size of the pad, the equivalent resistance R of 1400 kΩ and the equivalent capacitance C of 445.8 f F are calculated by R = ρ Mo l/s and C = ε r ε 0 s/d, where ρ Mo = 5.6 Ω·m, ε r = 9.5 and d = 2 µm are used. Based on the above data, the FEM simulations of the five resonators with the pad are performed. The measured impedance spectra are the blue curves in Figure 4, while the ones after de-embedding the pad are the red curves. As we can see, it is difficult to observe the resonant frequencies f r and anti-resonant frequencies f a in the blue curves. On the contrary, although the red curves still have some spurious modes, the main resonant modes can be clearly observed, and the two frequencies can be obviously distinguished, which proves that the method of Equation (22) has a good de-embedding effect. It should be stated that these spurious modes can be reduced by a better design, which is beyond the scope of this paper.   According to the size of the pad, the equivalent resistance of 1400 kΩ and the equivalent capacitance of 445.8 are calculated by = / and = / , where = 5.6 Ω • m, = 9.5 and = 2 μm are used. Based on the above data, the FEM simulations of the five resonators with the pad are performed. The measured impedance spectra are the blue curves in Figure 4, while the ones after de-embedding the pad are the red curves. As we can see, it is difficult to observe the resonant frequencies and anti-resonant frequencies in the blue curves. On the contrary, although the red curves still have some spurious modes, the main resonant modes can be clearly observed, and the two frequencies can be obviously distinguished, which proves that the method of Equation (22) has a good de-embedding effect. It should be stated that these spurious modes can be reduced by a better design, which is beyond the scope of this paper. Figure 4. The impedance spectra of resonators before (blue curves) and after (red curves) de-embedding the pad: (a) length thickness extension, (b) thickness extension, (c) first-order radial extension, (d) second-order radial extension, (e) lateral electric field excited thickness shear mode resonators and (f) lamb wave resonator. Figure 5 shows the variation of the scale factor γ in padless resonators with Mo electrodes thickness from 20 to 200 nm. Since the density of the electrode material is a constant, this figure also represents the relationship between γ and the mass of the electrodes. As shown, γ is linearly proportional to the thickness of electrodes at first, while with further increased thickness, it gradually deviates from the original straight line. It means that when the thickness of the electrodes is small ( / < 0.05), is indeed a constant, and the introduction of γ can effectively eliminate the influence of the electrodes. However, will deviate from the original constant value when the thickness  Figure 5 shows the variation of the scale factor γ in padless resonators with Mo electrodes thickness from 20 to 200 nm. Since the density of the electrode material is a constant, this figure also represents the relationship between γ and the mass of the electrodes. As shown, γ is linearly proportional to the thickness of electrodes at first, while with further increased thickness, it gradually deviates from the original straight line. It means that when the thickness of the electrodes is small (t Mo /t AlN < 0.05), n is indeed a constant, and the introduction of γ can effectively eliminate the influence of the electrodes. However, n will deviate from the original constant value when the thickness gradually increases, which may cause a certain degree of extraction error. Therefore, it is better to design relatively small thicknesses of the electrodes when using this strategy. gradually increases, which may cause a certain degree of extraction error. Therefore, it is better to design relatively small thicknesses of the electrodes when using this strategy. In order to verify the effect of iteration times on eight elastic and piezoelectric coefficients, the simulated iteration processes are performed using three sets of resonator samples with different AlN and Mo thicknesses (sample 1: 2 μm AlN, 50 nm Mo; sample 2: 2 μm AlN, 100 nm Mo; and sample 3: 1 μm AlN, 100 nm Mo) and two different material coefficients of AlN in [20,31] as initial coefficient 1 and initial coefficient 2, where the material coefficients of the AlN thin film in [15] are assumed to be the actual data that needs to be extracted. Moreover, the FEM simulations of these resonators are established with the material coefficients of the AlN thin film [15], and the obtained impedance spectra are considered as actual measured ones. The material coefficients of the AlN thin film are extracted according to the steps shown in Figure 3, where the dielectric constants = 8.2, = 10.8 and density = 3260 kg/m [15] are treated as known values since they are unchanged during iterations. The six sets of extraction results under different iteration times are shown in Figure 6, and the final extraction coefficients are compared in Table 2. It can be seen that each coefficient gradually approaches the actual value with the increase in iteration times, and basically coincides with it at the fourth iteration except for , whose extraction error is within 2% at the 8th, 13th and 14th iterations for sample 1, sample 2 and sample 3, respectively. It means that the iteration can indeed play a role in improving the extraction accuracy, and this strategy can achieve extremely high accuracy without considering actual error including measurement errors, manufacturing errors, etc. Furthermore, accurate and similar material coefficients are obtained in all six cases, which shows that this strategy can be applied to the AlN thin film with different thickness and electrode configuration. It is worth mentioning that the speed at which approaches the actual value is negatively correlated with the ratio / of Mo to AlN thickness, i.e., a relatively small electrode thickness is beneficial to the reduction in the In order to verify the effect of iteration times on eight elastic and piezoelectric coefficients, the simulated iteration processes are performed using three sets of resonator samples with different AlN and Mo thicknesses (sample 1: 2 µm AlN, 50 nm Mo; sample 2: 2 µm AlN, 100 nm Mo; and sample 3: 1 µm AlN, 100 nm Mo) and two different material coefficients of AlN in [20,31] as initial coefficient 1 and initial coefficient 2, where the material coefficients of the AlN thin film in [15] are assumed to be the actual data that needs to be extracted. Moreover, the FEM simulations of these resonators are established with the material coefficients of the AlN thin film [15], and the obtained impedance spectra are considered as actual measured ones. The material coefficients of the AlN thin film are extracted according to the steps shown in Figure 3, where the dielectric constants ε T 11 = 8.2, ε T 33 = 10.8 and density ρ = 3260 kg/m 3 [15] are treated as known values since they are unchanged during iterations. The six sets of extraction results under different iteration times are shown in Figure 6, and the final extraction coefficients are compared in Table 2. It can be seen that each coefficient gradually approaches the actual value with the increase in iteration times, and basically coincides with it at the fourth iteration except for e 15 , whose extraction error is within 2% at the 8th, 13th and 14th iterations for sample 1, sample 2 and sample 3, respectively. It means that the iteration can indeed play a role in improving the extraction accuracy, and this strategy can achieve extremely high accuracy without considering actual error including measurement errors, manufacturing errors, etc. Furthermore, accurate and similar material coefficients are obtained in all six cases, which shows that this strategy can be applied to the AlN thin film with different thickness and electrode configuration. It is worth mentioning that the speed at which e 15 approaches the actual value is negatively correlated with the ratio t Mo /t AlN of Mo to AlN thickness, i.e., a relatively small electrode thickness is beneficial to the reduction in the number of iterations required for e 15 . In addition, Table 2 also lists the results of direct extraction without eliminating the mass loading effect of the electrodes when the thicknesses of AlN and Mo are 2 µm and 100 nm, it obviously has a great error, c E 13 is even calculated as 0, which proves the importance of eliminating the influence of electrodes once again. machines 2022, 13,1 number of iterations required for . In addition, Table 2 also lists the results of extraction without eliminating the mass loading effect of the electrodes when the nesses of AlN and Mo are 2 μm and 100 nm, it obviously has a great error, is calculated as 0, which proves the importance of eliminating the influence of elect once again.  In summary, the above strategy can achieve a high extraction accuracy, withou of generality, and can even be used for other materials of 6 mm crystal class due  In summary, the above strategy can achieve a high extraction accuracy, without loss of generality, and can even be used for other materials of 6 mm crystal class due to its consistency with the crystal structure of AlN. The current study is carried out under ideal conditions, but in fact, the imperfections in the production process such as the inhomogeneity of the AlN thin film thickness may have a potential impact on the extraction accuracy, which will be further studied in subsequent work to evaluate the strategy more fully.

Conclusions
In this paper, a strategy for systematically extracting full elastic coefficients, piezoelectric coefficients, dielectric constants and density of the AlN thin film is proposed based on the resonance method, where five resonators with different modes are used, and TS and LE mode resonators in the conventional method are replaced by LEF-TS mode resonator and LWR, respectively. Furtherly, a complete extraction process is established in which the structure of the self-suspended resonator and the pad-only device are adopted to eliminate the effects of the substrate and pad. In addition, the scale factor γ (frequency without electrodes/frequency with electrodes) is introduced and iterated to eliminate the mass loading effect of electrodes. The simulation results show that the above methods can effectively decrease the extraction error and the material coefficients extracted by the strategy have high accuracy. This strategy can be used in c-axis-oriented AlN thin film of different thicknesses, and it can also be applied to other materials belonging to the 6 mm crystal class.  Data Availability Statement: All data and models generated or used during the study appear in the submitted article.

Conflicts of Interest:
The authors declare no conflict of interest.