Next Article in Journal
A Temperature-Compensated Single-Crystal Silicon-on-Insulator (SOI) MEMS Oscillator with a CMOS Amplifier Chip
Previous Article in Journal
Microhotplates for Metal Oxide Semiconductor Gas Sensor Applications—Towards the CMOS-MEMS Monolithic Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation of Material Constants of CaTiO3 Doped (K,Na)NbO3 Film by MEMS-Based Test Elements

1
Graduate School of Engineering, Tohoku University, 6-6 Aramaki aza, Aoba, Aoba-ku, Sendai, Miyagi 980-8579, Japan
2
New Industry Creation Hatchery Center, Tohoku University, 6-6 Aramaki aza, Aoba, Aoba-ku, Sendai, Miyagi 980-8579, Japan
3
Graduate School of Electrical Engineering, Tohoku University, 6-6 Aramaki aza, Aoba, Aoba-ku, Sendai, Miyagi 980-8579, Japan
4
Murata Manufacturing Co., Ltd., 1-10-1 Higashikotari, Nagaokakyo, Kyoto 617-8555, Japan
*
Author to whom correspondence should be addressed.
Micromachines 2018, 9(11), 558; https://doi.org/10.3390/mi9110558
Submission received: 28 September 2018 / Revised: 18 October 2018 / Accepted: 22 October 2018 / Published: 29 October 2018

Abstract

:
A CaTiO3-doped (K,Na)NbO3 (KNN-CT) film is a lead-free piezoelectric film that is expected to substitute Pb(Zr,Ti)O3 (PZT) film in piezoelectric micro electro mechanical systems (MEMS). However, the full set of the material constants (elastic constants, piezoelectric constants and dielectric constants) of the KNN-CT film have not been reported yet. In this study, all the material constants of a sputter-deposited blanket KNN-CT film were investigated by the resonance responses of MEMS-based piezoelectric resonators and the phase velocities of leaky Lamb waves on a self-suspended membrane. The phase velocities measured by a line-focus-beam ultrasonic material characterization (LFB-UMC) system at different frequencies were fitted with theoretical ones, which were calculated from the material constants, including fitting parameters. A genetic algorithm was used to find the best-fitting parameters. All the material constants were then calculated. Although some problems arising from the film quality and the nature of deliquescence are observed, all the material constants were obtained exhibiting accuracy within 16 m/s in the phase velocity of leaky Lamb wave.

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)NbO3 (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/SiO2/Si exhibited excellent transverse piezoelectric coefficient e31* from −10.0 to −14.4 C/m2 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 CaTiO3-doped KNN (KNN-CT) film deposited by RF sputtering showed excellent piezoelectric characteristics (e31* = −11.7 C/m2) [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 C6v 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:
f s = 1 2 l ρ s 11 E ,
k 31 2 1 k 31 2 = π 2 f p f s cot ( π 2 f p f s )   and
k 31 = d 31 s 11 E ε 33 T ,
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:
f s = ϕ 1 2 π r 1 ρ s 11 E ( 1 σ E 2 )   and
k r = d 31 2 ( s 11 E + s 12 E ) ε 33 T = k 31 2 1 σ E ,
where r and σE are the radius of the resonator and Poisson’s ratio, respectively. s12E is elastic compliance at E = 0. ϕ 1 is a constant determined by σE. Elastic compliance s66E is calculated from s11E and s12E as s66E = 2(s11Es12E).
For the thickness longitudinal mode, fp and coupling factor kt are expressed as:
f p = 1 2 h c 33 D ρ ,
k t 2 = π 2 f s f p cot ( π 2 f s f p )   and
c 33 E = c 33 D 1 k t 2 ,
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:
ε 33 T ε 0 = C f h ε 0 π r 2 ,
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:
V L = V w 1 ( 1 V w 2 f Δ z ) 2 ,
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 i = c i j E S j
T i = c i j E S j e m i E m
D n = e n j S j + ε n m T E m
div   D = 0
ρ 2 u h t 2 = T h k x k
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:
S h k = 1 2 ( u h x k + u k x h )   and
E m = φ x m ,
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:
u 1 α = A α exp ( Ω α k x 3 ) exp { j ( k x 1 ω t ) }
u 3 α = C α exp ( Ω α k x 3 ) exp { j ( k x 1 ω t ) }
φ α = D α exp ( λ α k x 3 ) exp { j ( k x 1 ω t ) }
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).
Fitness = 1 ( v A 0 ,   meas v A 0 ,   cal ) 2 + ( v S 0 ,   meas v S 0 ,   cal ) 2 ,
where the suffixes “meas”, “cal”, A0 and S0 represent measured value, calculated value, antisymmetric (A0) mode, and symmetric (S0) 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.

3. KNN-CT Film Deposition and Device Fabrication

3.1. KNN-CT Film Deposition

KNN-CT was deposited by RF sputtering on a 4 inch Si wafer of 500 µm thickness. As a seed layer, 150 nm thick SiO2, 10 nm thick Ti and 100 nm thick Pt were deposited in advance. The sputtering was done at a substrate temperature of 550 °C and a chamber pressure of 0.3 Pa using Ar and O2 with an O2/(Ar + O2) ratio of 0.01. The substrate temperature was optimized with film stress and crystallinity. Post annealing followed at 700 °C in O2 atmosphere for 10 min. The as-deposited film was self-polarized, as is often the case with sputtered piezoelectric films.
The crystal orientation of the KNN-CT film was analyzed by X-ray diffraction (XRD). Figure 5a,b show the XRD pattern and rocking curve, respectively. Although the main peak is (001), (110) and (210) peaks are also found. The full-width of half-maximum (FWHM) of (001) peak is 2.2°. This result suggests that the KNN-CT film is partially c-axis oriented. In addition, scanning nonlinear dielectric microscopy (SNDM) was conducted to analyze the distribution of polarity. The detail of SNDM is described in [25]. Figure 6 shows the image of SNDM, where the blight and dark areas show c-domain and other domains, respectively. c-domain occupies 50–60% of the area.
The density of the KNN-CT film was calculated from the difference of the mass and volume of the 4-inch wafer, before and after KNN-CT deposition, measured by the Archimedes method [26]. As a result, the density and thickness of the KNN-CT film were obtained as 4.46 × 103 kg/m3 and 1.75 µm, respectively.

3.2. Piezoelectric Resonators

The resonators of length expander mode, radial expander mode and thickness longitudinal mode were fabricated using MEMS process. The radial expander mode and thickness longitudinal mode were excited on the same resonator at different frequencies. The Pt seed layer was utilized as the bottom electrode of the resonators.
Figure 7 summarizes the fabrication process of the piezoelectric resonators. First, the KNN-CT-coated wafer was diced into 20 × 20 mm2 chips (1 of Figure 7). The KNN-CT film was etched by ion milling using a photoresist mask (2). Al top electrodes and the contacts of the bottom electrode of 200 nm thickness were formed by EB evaporation and lift-off process (3 and 4). Next, Pt/Ti/SiO2 layers were removed by ion milling to open etching windows for the following step (5). XeF2 etching was used to release the resonators at the final step (6). Figure 8a,b show the piezoelectric resonators of the length expander mode and the radial expander/thickness longitudinal mode, respectively. The distortion of the resonators happened after the release step, due to the film stress of the KNN-CT film.

3.3. Leaky Lamb Wave Device

A self-suspended Al-metalized KNN-CT membrane with a diameter of 2 mm was fabricated as a specimen for the LFB-UMC measurement. Figure 9 shows the fabrication process. First, the topside of the Si substrate, i.e., the KNN-CT film, was covered with photoresist. After photolithography on the backside, the Si substrate was etched by deep reactive ion etching (RIE) (1 and 2 of Figure 9). The etching was stopped before it reached the SiO2 layer to avoid damaging KNN-CT. Si left on the bottom of the etched hole was removed by XeF2 etching, and then the SiO2/Ti/Pt layers were removed by ion milling (3) and (4). The photoresist on both sides was removed by O2 ashing (5). Figure 10a shows microscope images of the KNN-CT membrane after O2 ashing. Finally, Al was deposited on the topside by EB evaporation (6), as shown in Figure 10b. The Al layer works as protection of the KNN-CT film from water absorption.

4. Results

4.1. Measurement of Piezoelectric Resonators

fs and fp of the piezoelectric resonators were measured by a network analyzer (E5071B, Agilent, Santa Clara, CA, USA). Figure 11a shows the relationship between the resonator length (l in Figure 8a) and fs for the length expander mode. As predicted by theory, fs is in inverse proportion to the length. From the best fit of Equation (1), s11E = 9.27 pm2/N is obtained, and k312 = 0.06% is obtained from Equation (2).
Figure 11b shows the relationship between the resonator radius (r in Figure 8b) and fs for the radial expander mode. fs is in inverse proportion to the radius. From the best fit of Equation (4), σE = 0.33, s12E = −σEs11E = −3.06 pm2/N and s66E = 24.7 pm2/N are obtained, and kr2 = 0.72% is obtained from Equation (5).
Figure 11c shows the relationship between the resonator radius (r in Figure 8b) and fp for the thickness longitudinal mode. The resonator response is very small, probably due to poor orientation of the KNN-CT film. fp is independent of the radius, showing that the measured response is for the thickness longitudinal mode. c33D = 89.8 GPa and kt2 = 2.08% are obtained from Equations (6) and (7), respectively, and c33E = 89.8 GPa is obtained from Equation (8).
Figure 11d shows the free capacitance Cf of the same type of resonators in terms of ε0πr2/h. As found in Equation (9), the relationship should be linear, and ε33T/ε0 = 1.27 × 103 is obtained by the least-squares method. Accordingly, d31 = −78.9 pC/N is calculated from Equation (3).

4.2. Phase Velocity Measurement by the LFB-UMC System

Figure 12a shows a typical V(z) curve measured for the KNN-CT film at 225 MHz. By applying FFT to the V(z) curve, Δz of 9.80 μm and 79.71 μm are obtained for A0 and S0 modes, respectively. From Equation (10), the phase velocities of 1988.91 m/s and 5225.88 m/s are calculated for A0 and S0 modes, respectively. Figure 12b shows the repeatability of the measured phase velocity at a single point. The frequency of the transducer is 225 MHz. The average values of phase velocities are 1988.46 m/s with ±2σ of ±0.71 m/s (±0.036%) and 5226.64 m/s with ±2σ of ±1.92 m/s (±0.037%) for A0 and S0 modes, respectively. σ is the standard deviation. Phase velocity was measured 30 times and was found to be almost stable. A similar process was repeated at different frequencies from 180 MHz to 270 MHz, where clear V(z) curve was obtained. Figure 13 shows the measured phase velocities of A0 and S0 modes at this frequency range.

4.3. Determination of Material Constants

k33, k15, s13E, c44E and d15 were determined by the algorithm described in Section 2.4 using the measured phase velocities shown in Figure 13. The initial values were set as k33 = 0.61, k15 = 0.64, s13E = −3.43 pm2/N, c44E = 37.0 GPa and d15 = 306 pC/N. These values are for hot-pressed pure KNN materials reported in [13]. M was set to 10, and the tournament size was set to 5. Pc and Pm were set to 90% and 5%, respectively. The calculation was repeated three times at different numbers of generation, and all of them were terminated after calculating all generations. The total number of generations were 500. The relationship between the fitness and the distribution of the fitting parameters in all generations were summarized in Figure 14a–e. There are many possible combinations of fitting parameters in the range of fitness less than 2 × 10−5. As shown in Figure 14a–e, possible combinations decreased as the fitness increased. The fitting parameters seem to converge over the fitness range of 10−4. The best fitness was obtained as 1.99 × 10−4 when k33 = 0.516, k15 = 0.776, s13E = −4.15 pm2/N, c44E = 10.5 GPa and d15 = 985 pC/N. Table 2 summarizes all the material constants of the KNN-CT film calculated by the fitting parameters, compared to [13]. By doping CT, piezoelectric and dielectric properties are enhanced, compared to the hot-pressed pure KNN. These properties are also higher than the reported ScAlN film [15]. The calculated phase velocities with the calculated material constants are in good agreement with the measured ones, as shown in Figure 13. The average of the root mean squared error of all the phase velocities over the measured frequencies Δ V 2 was 16.0 m/s at the best fitness.

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 A0 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 SiO2/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 Δ V 2 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 Δ V 2 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.

Author Contributions

R.K. performed all experiments except deposition of KNN-CT film, and made programs of genetic algorithm. M.K. provided scientific and technical advice. Y.O. provided technical support during LFB-UMC and density measurement. J.-i.K. designed the LFB-UMC system and provided technical support during measurement. S.I. deposited the KNN-CT film and provided technical advice for the film. S.T. supervised the entire work.

Funding

This study was supported by the Program for Leading Graduate Schools, ‘Inter-Graduate School Doctoral Degree Program on Global Safety’ of MEXT, Japan.

Acknowledgments

We would like to appreciate the help of Professor Yasuo Cho (Tohoku University) and his student in measuring the distribution of polarity of the KNN-CT film with SNDM.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Maeder, M.D.; Damjanovic, D.; Setter, N. Lead Free Piezoelectric Materials. J. Electroceram. 2004, 13, 385–392. [Google Scholar] [CrossRef]
  2. Shrout, T.R.; Zhang, S.J. Lead-free piezoelectric ceramics: Alternatives for PZT? J. Electroceram. 2007, 19, 111–124. [Google Scholar] [CrossRef]
  3. Panda, P.K. Review: Environmental friendly lead-free piezoelectric materials. J. Mater. Sci. 2009, 44, 5049–5062. [Google Scholar] [CrossRef]
  4. Shibata, K.; Oka, F.; Ohishi, A.; Mishima, T.; Kanno, I. Piezoelectric Properties of (K,Na)NbO3 Films Deposited by RF Magnetron Sputtering. Appl. Phys. Express 2008, 1, 011501. [Google Scholar] [CrossRef]
  5. Shibata, K.; Oka, F.; Nomoto, A.; Mishima, T.; Kanno, I. Crystalline Structure of Highly Piezoelectric (K,Na)NbO3 Films Deposited by RF Magnetron Sputtering. Jpn. J. Appl. Phys. 2008, 47, 8909–8913. [Google Scholar] [CrossRef]
  6. Shibata, K.; Suenaga, K.; Nomoto, A.; Mishima, T. Curie Temperature, Biaxial Elastic Modulus, and Thermal Expansion Coefficient of (K,Na)NbO3 Piezoelectric Thin Films. Jpn. J. Appl. Phys. 2009, 48, 121408. [Google Scholar] [CrossRef]
  7. Suenaga, K.; Shibata, K.; Watanabe, K.; Nomoto, A.; Horikiri, F.; Mishima, T. Effect of Lattice Strain and Improvement of the Piezoelectric Properties of (K,Na)NbO3 Lead-Free Film. Jpn. J. Appl. Phys. 2010, 49, 09MA05. [Google Scholar] [CrossRef]
  8. Shibata, K.; Suenaga, K.; Watanabe, K.; Horikiri, F.; Nomoto, A.; Mishima, T. Improvement of Piezoelectric Properties of (K,Na)NbO3 Films Deposited by Sputtering. Jpn. J. Appl. Phys. 2011, 50, 041503. [Google Scholar] [CrossRef]
  9. Shibata, K.; Suenaga, K.; Watanabe, K.; Horikiri, F.; Mishima, T.; Shiratani, M. Evaluation of Crystal Orientation for (K,Na)NbO3 Films Using X-ray Diffraction Reciprocal Space Map and Relationship between Crystal Orientation and Piezoelectric Coefficient. Jpn. J. Appl. Phys. 2012, 51, 075502. [Google Scholar] [CrossRef]
  10. Safari, A.; Abazari, M. Lead-Free Piezoelectric Ceramics and Thin Films. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 2010, 57, 2165–2176. [Google Scholar] [CrossRef] [PubMed]
  11. Li, P.; Chen, X.; Wang, F.; Shen, B.; Zhai, J.; Zhang, S.; Zhou, Z. Microscopic Insight into Electric Fatigue Resistance and Thermally Stable Piezoelectric Properties of (K,Na)NbO3-Based Ceramics. ACS Appl. Mater. Interfaces 2018, 10, 28772–28779. [Google Scholar] [CrossRef] [PubMed]
  12. Ikeuchi, S.; Yoneda, T.; Matsuki, Y.; Endo, N.; Takeshima, Y.; Horiuchi, H.; Kishimoto, Y.; Yamamoto, K.; Fujimoto, K. Preparation of (K,Na)NbO3-CaTiO3 Film by RF Magnetron Sputtering. In Proceedings of the IEEE International Urbana-Champaign, Chicago, IL, USA, 3–6 September 2014. [Google Scholar]
  13. Berlincourt, D. Piezoelectric Crystals and Ceramics. In Ultrasonic Transducer Materials, 1st ed.; Mattiat, O.E., Ed.; Plenum Press: London, UK, 1971; pp. 100–110. [Google Scholar]
  14. Tian, H.; Hu, C.; Meng, X.; Zhou, Z.; Shi, G. Dielectric, piezoelectric, and elastic properties of K0.8Na0.2NbO3 single crystals. J. Mater. Chem. C 2015, 3, 9609–9614. [Google Scholar] [CrossRef]
  15. Konno, A.; Kadota, M.; Kushibiki, J.; Ohashi, Y.; Esashi, M.; Yamamoto, Y.; Tanaka, S. Determination of Full Material Constants of ScAlN Thin Film from Bulk and Leaky Lamb Waves in MEMS-based Samples. In Proceedings of the IEEE International Urbana-Champaign, Chicago, IL, USA, 3–6 September 2014. [Google Scholar]
  16. Kushibiki, J.; Chubachi, N. Material Characterization by Line-Focus-Beam Acoustic Microscope. IEEE Trans. Sonics. Ultrason. 1985, SU-32, 189–211. [Google Scholar] [CrossRef]
  17. Kushibiki, J.; Ono, Y.; Ohashi, Y.; Arakawa, M. Development of the line-focus-beam ultrasonic material characterization system. IEEE Trans. Ultrason. Ferroelect. Freq. Contr. 2002, 49, 99–113. [Google Scholar] [CrossRef]
  18. Kadota, M.; Saito, Y.; Yoshida, T.; Ozeki, H.; Yamaguchi, T.; Hasegawa, M.; Sato, H. Electrical Test Methods for Piezoelectric Ceramic Vibrators; Japan Electronics and Information Technology Industories Association Chiyoda: Tokyo, Japan, 1993; pp. 39–55. (In Japanese) [Google Scholar]
  19. Parsapour, F.; Pashchenko, V.; Nicolay, P.; Muralt, P. Material constants extraction for AlScN thin films using a dual mode baw resonator. In Proceeding of the IEEE MEMS, Belfast, UK, 21–25 January 2018. [Google Scholar]
  20. Campbell, J.J.; Jones, W.R. Propagation of Surface Waves at the Boundary between a Piezoelectric Crystal and a Fluid Medium. IEEE Trans. Sonics. Ultrason. 1970, SU-17, 71–76. [Google Scholar] [CrossRef]
  21. Goldberg, D.E. Genetic Algorithms in Search, Optimization, and Machine Learning, 1st ed.; Addison-Wesley Publishing Company Inc.: Boston, MA, USA, 1989; pp. 1–25. [Google Scholar]
  22. Konak, A.; Coit, D.W.; Smith, A.E. Multi-objective optimization using genetic algorithms: A tutorial. Reliab. Eng. Syst. Saf. 2006, 91, 992–1007. [Google Scholar] [CrossRef]
  23. Goldberg, D.E.; Deb, K. A Comparative Analysis of Selection Schemes Used in Genetic Algorithms. In Foundations of Genetic Algorithms, 1st ed.; Rawlins, G.J.E., Ed.; Morgan Kaufmann Publishers: Burlington, MA, USA, 1991; Volume 1, pp. 69–93. [Google Scholar]
  24. Miller, B.L.; Goldberg, D.E. Genetic Algorithms, Tournament Selection, and the Effects of Noise. Complex Syst. 1995, 9, 193–212. Available online: http://wpmedia.wolfram.com/uploads/sites/13/2018/02/09-3-2.pdf (accessed on 27 September 2018).
  25. Cho, Y.; Kirihara, A.; Saeki, T. Scanning nonlinear dielectric microscope. Rev. Sci. Instrum. 1996, 67, 2297–2303. [Google Scholar] [CrossRef]
  26. Bowman, H.A.; Schoonover, R.M.; Carroll, C.L. A Density Scale Based on Solid Obiects. J. Res. Natl. Bur. Stand. A Phys. Chem. 1974, 78A, 13–40. Available online: https://nvlpubs.nist.gov/nistpubs/jres/78A/jresv78An1p13_A1b.pdf (accessed on 10 September 2018). [CrossRef]
  27. Zheng, L.; Huo, X.; Wang, R.; Wang, J.; Jiang, W.; Cao, W. Large size lead-free (Na,K)(Nb,Ta)O3 piezoelectric single crystal: Growth and full tensor properties. CrystEngComm 2013, 15, 7718–7722. [Google Scholar] [CrossRef]
  28. Kakimoto, K.; Masuda, I.; Ohsato, H. Lead-free KNbO3 piezoceramics synthesized by pressure-less sintering. J. Eur. Ceram. Soc. 2005, 25, 2719–2722. [Google Scholar] [CrossRef]
  29. Tanaka, D.; Tsukada, T.; Furukawa, M.; Wada, S.; Kuroiwa, Y. Thermal Reliability of Alkaline Niobate-Based Lead-Free Piezoelectric Ceramics. Jpn. J. Appl. Phys. 2009, 48, 09KD08. [Google Scholar] [CrossRef]
  30. Marcos, F.R.; Marchet, P.; Mejean, T.M.; Fernandez, J.F. Role of sintering time, crystalline phases and symmetry in the piezoelectric properties of lead-free KNN-modified ceramics. Mater. Chem. Phys. 2010, 123, 91–97. [Google Scholar] [CrossRef]
Figure 1. Piezoelectric resonance modes relative to polarization direction and displacement. (a) Length expander mode. (b) Radial expander mode. (c) Thickness longitudinal mode. (d) Cylinder longitudinal mode. (e) Thickness shear mode.
Figure 1. Piezoelectric resonance modes relative to polarization direction and displacement. (a) Length expander mode. (b) Radial expander mode. (c) Thickness longitudinal mode. (d) Cylinder longitudinal mode. (e) Thickness shear mode.
Micromachines 09 00558 g001
Figure 2. Schematic representation of the LFB-UMC system. (a) Cross-sectional geometry of the LFB acoustic lens and specimen to show V(z) curve measurement. (b) Typical V(z) curve as a transducer output.
Figure 2. Schematic representation of the LFB-UMC system. (a) Cross-sectional geometry of the LFB acoustic lens and specimen to show V(z) curve measurement. (b) Typical V(z) curve as a transducer output.
Micromachines 09 00558 g002
Figure 3. Analysis model of leaky Lamb wave propagation on self-suspended Al-metallized KNN-CT membrane. Lamb wave propagates to x1 direction.
Figure 3. Analysis model of leaky Lamb wave propagation on self-suspended Al-metallized KNN-CT membrane. Lamb wave propagates to x1 direction.
Micromachines 09 00558 g003
Figure 4. Flow chart of genetic algorithm for parameter fitting.
Figure 4. Flow chart of genetic algorithm for parameter fitting.
Micromachines 09 00558 g004
Figure 5. XRD results of KNN-CT film. (a) XRD pattern. (b) Rocking curve.
Figure 5. XRD results of KNN-CT film. (a) XRD pattern. (b) Rocking curve.
Micromachines 09 00558 g005
Figure 6. Polarity distribution of KNN-CT film. Bright and dark areas show c-axis orientation and other orientation, respectively.
Figure 6. Polarity distribution of KNN-CT film. Bright and dark areas show c-axis orientation and other orientation, respectively.
Micromachines 09 00558 g006
Figure 7. Process flow of piezoelectric resonators.
Figure 7. Process flow of piezoelectric resonators.
Micromachines 09 00558 g007
Figure 8. Piezoelectric resonators. (a) Length expander mode, l = 200–600 µm. (b) Radial expander mode and thickness longitudinal mode, r = 15–100 µm.
Figure 8. Piezoelectric resonators. (a) Length expander mode, l = 200–600 µm. (b) Radial expander mode and thickness longitudinal mode, r = 15–100 µm.
Micromachines 09 00558 g008
Figure 9. Process flow of self-suspended Al-metallized KNN-CT membrane for leaky Lamb wave.
Figure 9. Process flow of self-suspended Al-metallized KNN-CT membrane for leaky Lamb wave.
Micromachines 09 00558 g009
Figure 10. Topside microscopy images of KNN-CT membrane. The diameter is designed as 2 mm. (a) After O2 ashing. (b) After Al deposition.
Figure 10. Topside microscopy images of KNN-CT membrane. The diameter is designed as 2 mm. (a) After O2 ashing. (b) After Al deposition.
Micromachines 09 00558 g010
Figure 11. Measurement results of piezoelectric resonators. (a) Relationship between fs and length of length expander mode resonator. (b) Relationship between fs and radius of radial expander mode resonator. (c) Relationship between fp and radius of thickness longitudinal mode resonator. (d) Relationship between free capacity Cf and radial function of radial expander mode resonator. ε33T/ε0 is obtained from the tilt of the fitting curve.
Figure 11. Measurement results of piezoelectric resonators. (a) Relationship between fs and length of length expander mode resonator. (b) Relationship between fs and radius of radial expander mode resonator. (c) Relationship between fp and radius of thickness longitudinal mode resonator. (d) Relationship between free capacity Cf and radial function of radial expander mode resonator. ε33T/ε0 is obtained from the tilt of the fitting curve.
Micromachines 09 00558 g011
Figure 12. Measurement results by the LFB-UMC system. (a) Typical V(z) curve measured for the KNN-CT film at 225 MHz. (b) Transition of phase velocity of A0 and S0 modes, measured 30 times.
Figure 12. Measurement results by the LFB-UMC system. (a) Typical V(z) curve measured for the KNN-CT film at 225 MHz. (b) Transition of phase velocity of A0 and S0 modes, measured 30 times.
Micromachines 09 00558 g012
Figure 13. Phase velocities of A0 and S0 mode relative to measurement frequency f and the KNN-CT thickness H. Black broken lines are calculated phase velocities with best fitted parameters, as shown in Table 2.
Figure 13. Phase velocities of A0 and S0 mode relative to measurement frequency f and the KNN-CT thickness H. Black broken lines are calculated phase velocities with best fitted parameters, as shown in Table 2.
Micromachines 09 00558 g013
Figure 14. Plot of all generated fitting parameters related to fitness. (a) s13E. (b) k33. (c) k15. (d) c44E. (e) d15.
Figure 14. Plot of all generated fitting parameters related to fitness. (a) s13E. (b) k33. (c) k15. (d) c44E. (e) d15.
Micromachines 09 00558 g014
Figure 15. Different surface morphology of two test samples. (a) as-coated KNN-CT film. (b) after dipping into DI water for 36 h.
Figure 15. Different surface morphology of two test samples. (a) as-coated KNN-CT film. (b) after dipping into DI water for 36 h.
Micromachines 09 00558 g015
Figure 16. Deterioration of crystal orientation of KNN-CT film by dipping DI water. Deliquescence occurred at KNN-CT film.
Figure 16. Deterioration of crystal orientation of KNN-CT film by dipping DI water. Deliquescence occurred at KNN-CT film.
Micromachines 09 00558 g016
Table 1. Boundary conditions.
Table 1. Boundary conditions.
LayersBoundary ConditionsComment
Layer I and II
(x3 = H)
T31I = T31II = 0Air
T33I = T33II = 0Air
φI = φII
D3I = D3II
Layer II and III
(x3 = 0)
u1II = u1III
u3II = u3III
T31II = T31III
T33II = T33III
φII = φIII =0Metalize
D3II = D3III
Layer III and IV
(x3 = −hAl)
u3III = u3IV
T31III = T31IV
T33III = T33IV
Table 2. Determined material constants of KNN-CT film, compared to a hot-pressed pure KNN [13].
Table 2. Determined material constants of KNN-CT film, compared to a hot-pressed pure KNN [13].
This StudyRef. [13]Unit
Elastic compliance constants: sijE
s11E9.258.20pN/m2
s12E−3.05−2.54pN/m2
s13E−4.15−3.43pN/m2
s33E16.710.10pN/m2
s44E95.227.00pN/m2
s66E24.621.48pN/m2
Elastic stiffness constants: cijE
c11E162196.7GPa
c12E80.3103.6GPa
c13E60.1102.0GPa
c33E89.8168.0GPa
c44E10.537.0GPa
c66E40.646.6GPa
Piezoelectric constants: eij, dij
e1510.311.30C/m2
e31−5.67−2.36C/m2
e3310.69.80C/m2
d15985306pC/N
d31−79.0−51pC/N
d33223127pC/N
Relative dielectric constants: ɛijS0), ɛijT0)
ɛ11S760545-
ɛ33S901306-
ɛ11T1911496-
ɛ33T1270938-

Share and Cite

MDPI and ACS Style

Kaneko, R.; Kadota, M.; Ohashi, Y.; Kushibiki, J.-i.; Ikeuchi, S.; Tanaka, S. Investigation of Material Constants of CaTiO3 Doped (K,Na)NbO3 Film by MEMS-Based Test Elements. Micromachines 2018, 9, 558. https://doi.org/10.3390/mi9110558

AMA Style

Kaneko R, Kadota M, Ohashi Y, Kushibiki J-i, Ikeuchi S, Tanaka S. Investigation of Material Constants of CaTiO3 Doped (K,Na)NbO3 Film by MEMS-Based Test Elements. Micromachines. 2018; 9(11):558. https://doi.org/10.3390/mi9110558

Chicago/Turabian Style

Kaneko, Ryosuke, Michio Kadota, Yuji Ohashi, Jun-ichi Kushibiki, Shinsuke Ikeuchi, and Shuji Tanaka. 2018. "Investigation of Material Constants of CaTiO3 Doped (K,Na)NbO3 Film by MEMS-Based Test Elements" Micromachines 9, no. 11: 558. https://doi.org/10.3390/mi9110558

APA Style

Kaneko, R., Kadota, M., Ohashi, Y., Kushibiki, J. -i., Ikeuchi, S., & Tanaka, S. (2018). Investigation of Material Constants of CaTiO3 Doped (K,Na)NbO3 Film by MEMS-Based Test Elements. Micromachines, 9(11), 558. https://doi.org/10.3390/mi9110558

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop