1. Introduction
Pb(Zr,Ti)O3 (PZT) films are used for commercialized piezoelectric micro electro mechanical systems (MEMS), such as ink-jet printer heads and vibratory gyroscopes, due to their well-balanced performance. However, PZT includes lead (Pb) at high concentration, which may pollute the environment and eventually harm the human body. At present, there is no good replacement of PZT from a practical point of view. Therefore, lead-free piezoelectric films are important targets of research.
(K,Na)NbO
3 (KNN) is a replacement candidate for PZT [
1,
2,
3]. Shibata et al. investigated KNN films deposited by RF magnetron sputtering [
4,
5,
6,
7,
8,
9]. The sputtered KNN film on Pt/Ti/SiO
2/Si exhibited excellent transverse piezoelectric coefficient
e31* from −10.0 to −14.4 C/m
2 compared to other lead-free piezoelectric materials [
8]. Besides pure KNN, different kinds of KNN-based materials have also been widely studied [
10,
11]. Ikeuchi et al. reported that a CaTiO
3-doped KNN (KNN-CT) film deposited by RF sputtering showed excellent piezoelectric characteristics (
e31* = −11.7 C/m
2) [
12]. This level of performance can open the application of the KNN-CT film to piezoelectric MEMS. However, material constants (elastic, piezoelectric and dielectric) of the KNN-CT film, which are required for simulation to design MEMS, have not been revealed. For ceramic and single crystal KNN, full material constants have been reported [
13,
14]. The full material constants of a piezoelectric film, however, are difficult to determine, and no one has reported them for a KNN-CT film.
The methods to determine all the material constants of a piezoelectric films have been proposed. Konno et al. reported full material constants of a ScAlN film [
15], which were determined using MEMS-based bulk acoustic resonators and a line-focus-beam ultrasonic material characterization (LFB-UMC) system [
16,
17]. This method is partly based on that for bulk ceramics [
18].
Figure 1 shows each resonance mode relative to the polarization direction. Some of the material constants are measured by a length expander mode resonator (
Figure 1a), a radial expander mode resonator (
Figure 1b) and a thickness longitudinal mode resonator (
Figure 1c). However, all the material constants of a film cannot be determined as those of ceramics, because it is impossible to fabricate a cylinder longitudinal mode resonator (
Figure 1d) or a thickness shear mode resonator (
Figure 1e) using the c-axis oriented film. Therefore, the phase velocity of the leaky Lamb wave on a self-suspended film was measured using the LFB-UMC system, and parameter fitting using a theoretical model was then done to determine the material constants that were not directly measured [
15].
Parsapour et al. determined
c33 and
c44 of ScAlN, which were derived from the cylinder longitudinal mode and thickness shear mode, by using a dual mode bulk acoustic wave resonator with a tilted c-axis oriented film [
19]. The tilted c-axis oriented ScAlN film was grown on a non-electrode, i.e., insulating layer. Both methods by Konno et al. and Parsapour et al., however, have a limitation in deposition; the bottom electrode layer must be patterned before the deposition of a piezoelectric film. These methods cannot be used for a blanket piezoelectric film on a fully metal-coated bottom layer, which is often better in terms of film quality.
In this study, we investigated the material constants of a blanket KNN-CT film on a 4-inch wafer using acoustic resonators and the LFB-UFC system, inspired by [
15]. In a previous study, ScAlN was deposited on a Si substrate, where Pt bottom electrodes were patterned in advance. In this study, however, a blanket KNN-CT film deposited on a fully Pt-coated Si substrate without a pattern was characterized. This structure of sample is more general than the previous ones.
In addition, genetic algorithm was applied for the fitting of the phase velocity of the leaky Lamb wave. The fitting is an optimization technique to minimize errors between a theoretical model and measured data using several fitting parameters. For such an optimization problem, initial values should be close to global solutions to avoid the risk of local convergence. Compared to ScAlN, few material constants of KNN-CT were reported at present. Genetic approaches such as crossover, mutation and selection were introduced to escape local convergence and reach optimized fitting parameters from relatively far initial values.
2. Measurement Methods and Theoretical Equations
2.1. Piezoelectric Resonators
Piezoelectric resonators are used to measure the material constants for C
6v and C
∞v crystal structures, including c-axis oriented KNN. All material constants of bulk ceramics can be determined from series and parallel resonance frequencies (
fs and
fp) of the following modes of resonators: length expander, radial expander, cylinder longitudinal, thickness longitudinal and thickness shear [
18]. For the length expander mode,
fs and electromechanical coupling coefficient (coupling factor)
k31 are expressed as:
where
l and
ρ are the length of the resonator and density, respectively.
s11E,
d31, and
ε33T are elastic compliance at electric field
E = 0, piezoelectric
d constant and dielectric constant in stress-free condition (
T = 0), respectively.
For the radial expander mode,
fs and coupling factor
kr are expressed as:
where
r and
σE are the radius of the resonator and Poisson’s ratio, respectively.
s12E is elastic compliance at
E = 0.
is a constant determined by
σE. Elastic compliance
s66E is calculated from
s11E and
s12E as
s66E = 2(
s11E −
s12E).
For the thickness longitudinal mode,
fp and coupling factor
kt are expressed as:
where
h,
c33D and
c33E are the thickness of the resonator and elastic stiffness constants at electrical displacement
D = 0 and
E = 0, respectively.
By measuring free capacitance
Cf of the radial expander mode resonators, the relative dielectric constant at
T = 0,
ε33T⁄
ε0, is calculated as:
where
r and
ε0 are the radius of the resonator and the dielectric constant of vacuum, respectively.
Using the above resonators,
s11E,
s12E,
s66E,
c33E,
k31,
kr,
kt,
d31,
σE and
ε33T are determined. Elastic compliance
s33E, coupling factor
k33 and piezoelectric d constant
d33 are obtained by the cylinder longitudinal resonator, and elastic stiffness
c44E, coupling factor
k15, elastic compliance
s44E, dielectric constant
ε11T and piezoelectric d constant
d15 are obtained by the thickness shear resonator. All other constants (e.g., elastic compliance
s13E, elastic stiffness
c11E,
c12E,
c13E and
c66E, piezoelectric e constants
e31,
e33 and
e15, dielectric constants in strain-free condition (
S = 0)
ε11S and
ε33S) are obtained from the above constants. The detailed flow is described in [
18]. However, both cylinder longitudinal and thickness shear resonators cannot be fabricated using a c-axis oriented KNN-CT film. Instead, the phase velocity of leaky Lamb wave in a suspended film is used, as described in
Section 3.3.
2.2. Phase Velocities of Leaky Lamb Wave
The phase velocity of leaky Lamb wave is affected by most material constants. The relationship between the phase velocity and the material constants is not explicit; however, the material constants not available from the resonators can be determined by fitting the calculated phase velocity to the measured one in terms of different modes of leaky Lamb wave [
15].
The phase velocity of leaky Lamb wave is measured by the LFB-UMC system.
Figure 2a shows the schematics of the LFB-UMC system. Longitudinal bulk wave is generated by applying RF pulse to the transducer, focused on by the cylindrical lens at the end of the sapphire rod. The focused beam is incident on the specimen, i.e., the self-suspended KNN-CT film, via the water coupler (Path #1), and then converted into leaky Lamb waves. The leaky Lamb waves propagate in the suspended film, emitting acoustic energy into water. The emitted bulk waves return to the transducer via Path #2. On the other hand, the normal incident wave reflects vertically at the specimen surface and returns via Path #0.
By changing the relative distance
z between the sapphire rod lens and the specimen, the transducer output
V(
z) curve is obtained, as illustrated in
Figure 2b. After the removal of background and lens characteristics and the extraction of an interference component from
V(
z) curve, the oscillation interval
z is obtained by fast Fourier transform (FFT). Finally, the phase velocity of LSAW
VL is given as:
where
VW and
f are the longitudinal wave velocity of water and measurement frequency, respectively.
2.3. Theoretical Phase Velocity of Lamb Wave
The theoretical phase velocity of Lamb waves is calculated by Campbell’s method [
20]. In the LFB-UMC system, a water coupler exists on one side of the membrane, as shown in
Figure 2a, and thus Campbell’s method is extended to calculate the phase velocities of leaky Lamb wave.
Figure 3 shows an analysis model used in this study. The top of the KNN-CT membrane is covered with Al, which is exposed to the water coupler.
H and
hAl are thickness of the KNN-CT membrane and Al, respectively. Acoustic wave propagation in this model is governed by Hook’s law (Equation (11)) for non-piezoelectric materials (air, Al and water) or piezoelectric fundamental equations (Equations (12) and (13)) for KNN-CT, Maxwell’s equation (Equation (14)) and Newton’s equation of motion (Equation (15)).
T,
S,
E,
D,
u and
t are stress, strain, electric field, electric displacement, particle displacement and time, respectively.
x is the propagation direction of acoustic wave. The suffixes (
i,
j), (
n,
m) and (
h,
k) take 1–6, 1–3 and 1–3, respectively.
S and
E are expressed as:
where
φ is an electric potential. Note that
Sj is the abbreviated notation of
Shk. The particle displacements
u1 and
u3, and
φ of laterally propagating acoustic waves in each layer are assumed as follows:
Aα,
Cα and
Dα are amplitudes, Ω
α and
λα are decay constants, and
k and
ω are wave number and angular frequency, respectively. The superscript
α indicates each layer (I–IV) in
Figure 3.
For non-piezoelectric film layers (I, III and IV), Equations (18) and (19) are substituted for the displacements in Newton’s equation of motion (Equation (15)) using Equations (11) and (16). The eigenvalues are obtained from the characteristic equation of Ωα as ΩlI (l = 1–2), ΩmIII (m = 1–4) and ΩIV. Note that only negative ΩlI and positive ΩIV are adopted to prevent divergence at positive and negative infinity in x3 direction, respectively. By substituting Equation (20) for the electric potential in Maxwell’s equation (Equation (14)) using Equation (17), the eigenvalues are obtained from the characteristic equation of λα as λI, λlIII and λIV. Note that only negative λI and positive λIV are adopted for the same reason. For the piezoelectric film layer (II), Equations (18)–(20) are substituted for the displacements and electric potential in Equations (14) and (15) using Equations (12), (13), (16) and (17). In a similar way, eigenvalues of ΩnII (n = 1–6) and λlII are obtained.
The general solutions of Equations (14) and (15) are given by Equations (18)–(20) with the eigenvalues for each layer. The boundary conditions summarized in
Table 1 are applied to the general solutions, and a simultaneous equation in terms of
AlI,
AnII,
AmIII and
AIV are obtained. This equation has non-zero, i.e., nontrivial, solutions at particular
v =
ω/
k, where
v is a complex phase velocity with an imaginary part of leaky component.
v is numerically searched assuming the material constants of the KNN-CT, Al and water. By repeating this calculation at different
ω, the frequency characteristic of the phase velocity is obtained.
2.4. Parameter Fitting Using Genetic Algorithm
All of the material constants are determined by solving the optimization problem to minimize the error between the measured phase velocities of leaky Lamb wave and the theoretical model. The parameters to be fitted are related to the cylinder longitudinal mode and the thickness shear mode, which are not measured by the piezoelectric resonators. Unfortunately, proper initial values are not available in this study because a few material constants of KNN-CT have been reported. Least squares methods (Newton method, Gauss-Newton method or Levenberg-Marquardt method) are not options to prevent local convergence.
In this study, genetic algorithm was employed for parameter fitting. Genetic algorithm explores the best solutions of such optimization problems based on the mechanisms of natural selection, genetic combining, and mutation [
21,
22]. The role of mutation is to escape local convergence. Here, let fitting parameters be specific material constants,
k33,
k15,
s13E,
c44E and
d15. A solution vector of integer number and that of binary number to operate genetic algorithm are called
individual and
chromosome, respectively. The chromosomes are made of discrete units, called
genes.
Population is used as a collection of the individuals or the chromosomes.
Figure 4 shows the flow chart of the developed algorithm using genetic algorithm. The algorithm starts with setting the initial values of the fitting parameters (step 1). Initial population with size
M is then defined for each fitting parameter (step 2). The initial value of each fitting parameter is transformed into 3-digits integer part and index part. The integer part is used for individuals in the population, whereas the index part is used later to transform individuals into the fitting parameters. Random integer matrix (
M × 1), which is 100 times multiplied Gauss-Laplace random number (mean 0, variance 1), is defined, and then each row of the matrix is added to the integer part to form the initial population. The initial population of each fitting parameter is formed independently in this way. Next, the individuals in the initial populations are transformed into binary digits to express as chromosomes (chromosome encoding, step 3). The populations of binary chromosomes are defined in this way.
Next, two of the chromosomes are selected randomly for crossover among the respective populations (step 4). The selected chromosomes become parents to form new chromosomes (offspring) with specific crossover probability Pc. The selection is repeated until the number of pair of parents becomes M. As a crossover method, uniform crossover is employed, which forms the chromosomes of offspring by exchanging each gene in a chromosome of parents at specific probability Pu (50%). The size of the population for offspring becomes 2M × 1 since two offspring are formed from single pair of parents. The crossover is conducted on each population independently. After that, the genes of chromosomes among the offspring are randomly changed at specific probability Pm as mutation, which makes it possible to escape local convergence (step 5).
In step 6, natural selection is done by selecting chromosomes with the number
M among the offsprings toward a parameter, which is called
fitness. Fitness is defined with theoretical phase velocities, thus the fitting parameters of the offspring are necessary. To obtain fitting parameters, the individuals of the offspring are calculated by transforming the binary chromosomes into a decimal integer. The fitting parameters of the offsprings are then calculated by multiplying the individuals with the index part, which was removed at the definition of initial populations. The populations of fitting parameters of offspring are defined in this way. After that, the combinations of fitting parameters (e.g.,
k33,
k15,
s13E,
c44E and
d15), which are the candidates of the solutions in this generation, are formed in the same row of the populations. The theoretical phase velocities are then calculated for each combination, according to the calculation described in
Section 2.3.
Fitness is defined as an inverse function of the squared error between the calculated frequency characteristic and measured one, as Equation (21).
where the suffixes “meas”, “cal”, A
0 and S
0 represent measured value, calculated value, antisymmetric (A
0) mode, and symmetric (S
0) mode, respectively. Tournament selection was employed to select an offspring with the number
M, in which combinations of fitting parameters with a specific number (tournament size) are randomly selected; a combination with the best fitness was selected as a winner [
23,
24]. The selection is repeated until the winning number, which is also called the next generation, becomes
M. Finally, maximum fitness
Fitnessmax in the next generation is compared with threshold
Fitnessthres (step 7). If the maximum fitness is less than the threshold, the populations of the next generation return to the chromosome encoding for the next calculation. Else, the fitting parameters of the best fitness are output as the best solution (step 8) and then substituted to calculate all material constants (step 9) [
18]. The loop is also broken if the number of generation
Ngen become a specific number
Nthres.
5. Discussion
We determined all elastic constants, piezoelectric constants, and dielectric constants of the KNN-CT film. In this section, the accuracy of the determined constants is discussed.
First, a clear resonance response is needed to achieve reasonable accuracy in the measurement of resonance frequency. If the quality of a piezoelectric film is poor, the exact resonance frequencies are often difficult to determine, and the material constants become inaccurate. This is the case with the KNN-CT film prepared in this study. The thickness longitudinal mode, especially, showed a very small resonance response. As a result, the measured
c33E and coupling factors are significantly smaller than those of hot-pressed pure KNN, single crystal and KNN-based ceramics [
13,
14,
27].
The second major source of inaccuracy is the measurement error of phase velocities by LFB-UFC. The measured
V(
z) curve of leaky Lamb wave, especially A
0 mode, (
Figure 12a), had smaller interference peaks than that of bulk materials [
16,
17]. This is probably because deliquescence occurs on the KNN-CT membrane during measurement by the water coupler, which fills a gap on the cylindrical lens. It is well known that alkaline niobate-based ceramics, such as KNN, shows deliquescence and a reaction with moisture in air [
28,
29,
30].
To confirm the deliquescence, two of the KNN-CT-coated chips were prepared. One of them was immersed in DI water for 36 h. Surface morphology and crystal orientation were then observed by Atomic Force Microscopy (AFM) and XRD, respectively.
Figure 15a,b show the surface morphology of the as-coated and water-dipped films, respectively. The tent-like crystalline structures in
Figure 15b were thinner than that in
Figure 15a. In
Figure 16, the FWHM of the water-immersed chip and another one is 2.699° and 2.365°, respectively, which clearly shows the deterioration of crystal orientation by DI water. The KNN-CT film, therefore, has deliquescence, and thus easily reacts with moisture in the air. Although the KNN-CT membrane was covered with the Al film, the membrane was exposed to the water coupler for a long time during the LFB-UMC measurement. The crystalline deterioration may have occurred by moisture during the measurement, which decreased the interference peak. Our proposed method, therefore, may not be suitable for materials that react with water.
The third major source of inaccuracy is in the fabrication of the self-suspended Al-metalized KNN-CT membrane for the LFB-UMC measurement. The ion milling process to remove SiO
2/Ti/Pt layers does not have a selectivity to the KNN-CT film. Therefore, the thickness of the KNN-CT film could be thinner than the original one (1.75 μm), which affects the calculation results. Moreover, the self-suspended KNN-CT membrane is easy to break during the measurement. We also prepared the self-suspended KNN-CT membrane without an Al metallization layer, but it broke during the measuring process, probably because of the deliquescence and stress caused by the water coupler. The KNN-CT membrane with the Al metallization layer survived without cracks, but we needed to consider the Al metallization layer in the calculation described in
Section 2.3. The error of the material constants of Al also affects the material constants. One solution is to measure leaky Rayleigh wave instead of leaky Lamb wave, because it is not necessary to fabricate the self-suspended membrane. However, numerical calculation is more complicated, because it requires several layers (seed, adhesion, isolation and substrate). The inaccuracy of the material constants of such layers is also a source of error.
As shown in
Figure 14a–e, increasing fitness is an index of the convergence of the fitting parameters. Fitness is theoretically infinite when the calculated phase velocities fit the measured sample completely. After the calculation of 500 generations, fitness was actually limited to the 10
−4 range, where
is less than 20 m/s. This may be because the variance of
k33 and
d15 is relatively large compared to that of
s13E,
k15 and
c44E, which are almost convergent over the fitness range of 10
−4. Further convergence of
k33 and
d15 may improves the fitness.
6. Conclusions
In this study, we investigated the full set of material constants of a blanket KNN-CT film for the first time. First, some of the material constants were measured from the resonance responses of piezoelectric resonators, which were fabricated by the MEMS process. Next, the phase velocities of leaky Lamb wave (A0 and S0 modes) were measured on a self-suspended Al-metalized KNN-CT membrane using a LFB-UMC system. The other constants were determined by fitting the measured phase velocities with theoretical ones, which were calculated from the material constants including multiple fitting parameters, i.e., unknown constants. After the fitting using genetic algorithm, the theoretical phase velocities were well fitted with error of 16.0 m/s. The KNN-CT film measured in this study was not of high quality, and the resonance response were small and unclear, especially for the thickness longitudinal mode. In addition, the KNN-CT film has deliquescence, which affected LFB-UMC measurement with a water coupler. The sample of the LFB-UMC measurement, the suspended Al-metalized KNN-CT membrane, also involves unconfident properties such as thickness and material constant of electrodes. Some of the determined material constants, especially d15, may be affected by such problems. Nevertheless, all material constants were successfully determined, and the viability of the proposed method was demonstrated.