Next Article in Journal
Ultra-Precision Cutting and Characterization of Reflective Convex Spherical Blazed Grating Elements
Next Article in Special Issue
Frequency Scanning Dual-Mode Asymmetric Dual-OAM-Wave Generation Base on Broadband PB Metasurface
Previous Article in Journal
Ergonomics Design and Assistance Strategy of A-Suit
Previous Article in Special Issue
Design and Measurement of a Novel Overmoded TE01 Mode Converter for a Rectangular Gyro-TWT
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A General FEM Model for Analysis of Third-Order Nonlinearity in RF Surface Acoustic Wave Devices Based on Perturbation Theory

1
College of Information, Mechanical and Electrical Engineering, Shanghai Normal University, Shanghai 200234, China
2
Key Laboratory of Optoelectronic Material and Device, Department of Physics, Shanghai Normal University, Shanghai 200234, China
3
The School of Materials Science and Engineering, Tsinghua University, Beijing 100084, China
4
Shoulder Electronics Limited, Wuxi 214124, China
*
Authors to whom correspondence should be addressed.
Micromachines 2022, 13(7), 1116; https://doi.org/10.3390/mi13071116
Submission received: 16 June 2022 / Revised: 8 July 2022 / Accepted: 13 July 2022 / Published: 15 July 2022
(This article belongs to the Special Issue Microwave Passive Components)

Abstract

:
This article presents a general-purpose model that enables efficient and accurate calculation of third-order nonlinear signals in surface acoustic wave (SAW) devices. This model is based on piezoelectric constitutive equations combined with perturbation theory, which can be analyzed by full finite element method (FEM). For validation, third-order harmonic (H3) responses and intermodulation distortions (IMD3) in SAW resonators are simulated, and their calculation results fit well to experimental data in the literature. Then, the generation mechanisms of the third-order nonlinearity in SAW resonators are discussed. The dominant generation mechanisms for different nonlinear signals and the relation between electrode materials and H3 peak magnitude are revealed, which provides an important guideline for further nonlinear suppression.

1. Introduction

Surface acoustic wave and bulk acoustic wave (SAW/BAW) devices used as radio frequency (RF) front-end filters and duplexers are mass-produced and applied in telecommunication systems [1,2]. Recent advances in mobile communication technologies have led to multiple signal handling, downsizing, and higher power level operations [3]. In this circumstance, the nonlinearity in SAW/BAW has become an important issue, as it generates harmonic and intermodulation distortion (IMD) products, which result in noticeable signal distortions [4,5,6].
Numerous groups have made efforts to investigate the nonlinearity in SAW/BAW devices and its suppression. For BAW devices, the one-dimensional (1D) Mason equivalent circuit model and modified Butterworth Van Dyke (MBVD) model are commonly employed [3,7,8]. Shim and Feld [9] proposed a 1D nonlinear Mason model using the harmonic balance (HB) technique [10,11]. The proposed Mason model is applicable to arbitrary piezoelectric nonlinear sources, and it simulates well for nonlinear signals generated in RF BAW. Considering the nonlinearity in RF BAW devices is extremely weak, Hashimoto et al. [12,13] pointed out that the perturbation method should be more applicable than the HB method. They demonstrated that the nonlinear behaviors of the thickness extensional resonators for both 1D and 2D cases are accurately simulated by the MBVD model with the first-order perturbation method. Chen et al. [14,15] also applied the above method to calculate third-order nonlinear distortion of SAW duplexers. The validity and accuracy were established by comparing simulation and measurement results.
As for SAW devices, the classical coupling-of-modes (COM) and P-matrix model methods are widely used to characterize their nonlinear behaviors. Nakagawa et al. [16] derived a COM model by introducing the nonlinear stress and electric displacement to identify the contributions from different mechanisms to nonlinearity in SAW devices. Furthermore, Chauhan et al. [17,18] extended it to P-matrix formalism to analyze that of the temperature compensated SAW (TC-SAW) devices on 128° YX LiNbO3 substrate. In the latest publications, a finite element method (FEM) model was applied by Mayer et al. [19], Guan et al. [20] and Pang et al. [21] to study the intermodulation and harmonic generation in LiNbO3 based SAW resonators. However, the above nonlinear FEM models are limited to SAW resonators on piezoelectric substrate such as Quartz and LiNbO3, whose higher-order material constants are known. This means they have difficulty in analyzing the cases when higher-order material constants are unknown, for example, the commercial commonly used LiTaO3.
With the coming of the 5G era, SAW/BAW and hybrid SAW-BAW devices with complicated structures are emerging one after another. Although the above publications have illustrated the generation mechanisms of nonlinearity and proposed several nonlinearity suppression methods, a general-purpose and efficient simulation tool for both SAW and BAW devices is still absent. In our previous work [22], we proposed a full FEM model for analyzing nonlinearity in BAW resonators and verified its effectiveness. Therefore, in this paper, we extend this model to calculate nonlinearity in SAW resonators. Theoretical derivations of this method are based on piezoelectric constitutive equations and perturbation theory, by which multiple piezoelectric nonlinearity can be considered simultaneously and combined arbitrarily. Compared to the previous FEM model [19,20,21], the advantage of the proposed model is that it not only inherently has the universality of a finite element, but also remains effective even if higher-order material constants of piezoelectric substrate are unknown. Validations were made by comparing the third-order nonlinear responses obtained by the proposed method with experimental results in the literature qualitatively. In addition, possible nonlinearity generation mechanisms and their suppression are discussed.

2. Theoretical Background and Analysis Procedures

2.1. Linear Equations of Piezoelectricity

We start by linear cases for piezoelectric resonators, the linear constitutive equations with charge-stress form (namely e-form) describe piezoelectric coupling as follows:
T = c E S e E
And
D = e S + ε S E
where T, S, E and D are the stress, strain, electric field and electric displacement, respectively, and c E , e and ε S are the elastic stiffness constant under constant E, piezoelectric coefficient and dielectric constant under constant S, respectively. Note that the subscripts of variables and material constants were omitted to simplify the following derivation.
The motion of the piezoelectric body is governed by Newton’s second law:
  ·   T = ρ 2 u t 2
where ρ is the mass density and u is displacement, and the charge equation of electrostatics is given by:
  ·   D = 0
Therefore, the electric displacement D is spatially uniform.

2.2. Nonlinear Equations of Piezoelectricity

On the basis of the perturbation theory [23], we modify the linear piezoelectric constitutive equations by introducing the nonlinear stress T N and the nonlinear electric displacement D N perturbations as follows:
T = c E S e E + T N
And
D = e S + ε S E + D N
where T N and D N can be represented as a series consisting of integer powers of their variables S and E .
According to the thermodynamics of solids [24], we derived the expression of perturbation terms T N and D N by expanding Helmholz free energy A until the fourth-order terms of strain S and electric field E :
T N = 1 2 ! ( 3 A S 3 ) S 2 + ( 3 A S 2 E ) S E + 1 2 ! ( 3 A S E 2 ) E 2 + 1 3 ! ( 4 A S 4 ) S 3 + 1 2 ( 4 A S 3 E ) S 2 E + 1 2 ( 4 A S 2 E 2 ) S E 2 + 1 3 ! ( 4 A S E 3 ) E 3
And
D N = 1 2 ! ( 3 A S 2 E ) S 2 ( 3 A S E 2 ) S E 1 2 ! ( 3 A E 3 ) E 2 1 3 ! ( 4 A S 3 E ) S 3 1 2 ( 4 A S 2 E 2 ) S 2 E 1 2 ( 4 A S E 3 ) S E 2 1 3 ! ( 4 A E 4 ) E 3
Equations (7) and (8) are then rewritten as Equations (9) and (10), respectively, as follows:
T N = 1 2 χ 20 T S 2 + χ 11 T S E + 1 2 χ 02 T E 2 + 1 6 χ 30 T S 3 + 1 2 χ 21 T S 2 E + 1 2 χ 12 T S E 2 + 1 6 χ 03 T E 3
And
D N = 1 2 χ 11 T S 2 χ 02 T S E 1 2 χ 02 D E 2 1 6 χ 21 T S 3 1 2 χ 12 T S 2 E 1 2 χ 03 T S E 2 1 6 χ 03 D E 3
where χ i j T and χ i j D are nonlinear coefficients, and the superscripts “T” and “D” indicate contributions of T N and D N . Note that the dominated terms in T N and D N expressions are related to the order of nonlinearity under consideration.
Substituting Equations (5) and (6) into Equations (3) and (4) derives:
· ( c E S e E ) ρ 2 u t 2 = · T N
And
· ( e S + ε S E ) = · D N
Nonlinear responses in the piezoelectric resonator can be determined by linear partial differential equations of Equations (11) and (12) where the terms on the right-hand side are taken as perturbations.

2.3. Derivation of Third-Order Nonlinear Responses

For simplicity, it is assumed that the excitation with two different frequency components f 1 and f 2 are applied to a piezoelectric resonator. Therefore, the linear strain and electric field in piezoelectric material are expressed as:
S = S f 1 + S f 2
And
E = E f 1 + E f 2
The perturbations caused by nonlinear products including intermodulation distortions and harmonics responses can be estimated by substituting Equations (13) and (14) into Equations (9) and (10). The nonlinear stress T N and nonlinear electric displacement D N for third-order harmonics (H3) are derived as:
T N 3 f 1 = 1 6 χ 30 T S f 1 3 + 1 2 χ 21 T S f 1 2 E f 1 + 1 2 χ 12 T S f 1 E f 1 2 + 1 6 χ 03 T E f 1 3
And
D N 3 f 1 = 1 6 χ 21 T S f 1 3 1 2 χ 12 T S f 1 2 E f 1 1 2 χ 03 T S f 1 E f 1 2 1 6 χ 03 D E f 1 3
In a similar way, the expressions of T N and D N for third-order intermodulation distortions (IMD3) cases are derived as:
T N 2 f 1 ± f 2 = 1 2 χ 30 T S f 1 2 S f 2 + 1 2 χ 21 T ( 2 S f 1 E f 1 S f 2 + S f 1 2 E f 2 ) + 1 2 χ 12 T ( 2 S f 1 E f 1 E f 2 + S f 2 E f 1 2 ) + 1 2 χ 03 T E f 1 2 E f 2
And
D N 2 f 1 ± f 2 = 1 2 χ 21 T S f 1 2 S f 2 1 2 χ 12 T ( 2 S f 1 E f 1 S f 2 + S f 1 2 E f 2 ) 1 2 χ 03 T ( 2 S f 1 E f 1 E f 2 + S f 2 E f 1 2 ) 1 2 χ 03 D E f 1 2 E f 2
It is noted that perturbation terms T N and D N in Equations (5) and (6) are composed by the combination of the electrostatic field, the strain field, and the mixing of the strain field and the electrostatic field [25], which can be selectively estimated by choosing specific expressions from Equations (15)–(18) for different nonlinear responses.

2.4. Analysis Procedures of Nonlinear Signals

Figure 1a shows a schematic of a SAW resonator under nonlinear test with 50 Ω matching impedance. The SAW resonator with an infinitely long interdigital transducers (IDTs) structure, comprised of periodic metal electrodes, is considered. For modeling of this device, a quasi-3D periodic FEM model (shown in Figure 1b) is built and used for numerical simulations by using a built-in piezoelectric module of FEM software COMSOL Multiphysics 5.6, as two-dimensional (2D) models are insufficient in analyzing the horizontal shear field component along the aperture y -direction. As shown, one period with periodic metal electrodes is considered, the continuity periodic boundary condition is applied to field variables at the left surface ( Γ L ) and the right surface ( Γ R ). For modeling accuracy, an air layer overlay is added to the top of the IDTs to consider the parallel capacitance effect. The perfectly matched layer (PML) is applied to the bottom to reduce the model size and suppress the unwanted boundary reflection. The periodic metal electrodes are applied with a terminal of one voltage and ground, respectively, for linear analysis, and then the terminal is changed to zero voltage for nonlinear analysis. As for the calculation of nonlinear responses, co-simulation of the quasi-3D FEM model of the resonator with its peripheral circuit is performed by using LiveLink of COMSOL with MATLAB.
The SAW resonator is characterized by the harmonic analysis, and the harmonic admittance Y per IDT period is estimated by Y = 2 π f j Q / U , where f is the driving frequency, U is the applied electric potential and Q is the total charge induced on the electrode. Figure 2 shows a flow chart illustrating the procedures for analysis of nonlinear responses. As shown, the linear input admittance Y 11 of the resonator at both driving frequency and output frequency is first evaluated by a quasi-3D FEM model using COMSOL. Meanwhile, the linear strain S and electric field E at same frequency spectrum are also obtained. Then, the nonlinear terms T N and D N , combination terms of linear strain S and electric field E with nonlinear coefficients χ i j , are estimated and added into the linear model as perturbations. To be specific, T N and D N can be added as mechanical and electric loading into “solid mechanics” and “electrostatics” interfaces in COMSOL, respectively. Next, the effects of a peripheral circuit, such as die, package and matching impedance effects, are taken into account using LiveLink with a MATLAB interface. Finally, different nonlinear responses can be obtained by solving the nonlinear piezoelectric constitutive equations at an aiming output frequency range, provided that a good priori of nonlinear coefficients χ i j are given. It is noted that priori values of nonlinear coefficients χ i j are obtained by fitting simulation results with experimental results. In the fitting procedure, nonlinear coefficients χ i j were tried one by one and optimized by minimizing the absolute difference between the simulated and measured results. Additionally, multiple nonlinear coefficients can be considered simultaneously to find the best agreement with the measurement.

3. Simulation Results and Validations Examples

To confirm validity for SAW devices of the proposed model, a SAW resonator on 42° YX LiTaO3 (42-LT) substrate is taken as an example. Figure 3 shows the reflection coefficient S 11 of the 42-LT SAW resonator. As shown, the S 11 curve calculated by the quasi-3D FEM model fits fairly well with the measured results in [16]. A resonant frequency f r exists at 837.1 MHz and an anti-resonant frequency f a at 866.8 MHz.
For H3 simulation, a continuous wave (CW) signal was used and its incident power level is 15 dBm. The driving frequency is swept from 800 to 900 MHz, and an H3 signal appears at triple the driving frequency, namely from 2.4 to 2.7 GHz. As shown in Figure 4a, the calculated H3 frequency dependence also compares well with that of the experiment in [16], particularly, the peak and notch shape. A peak with the maximum H3 level at 2.56 GHz can be seen and a steep notch occurs at about 2.6 GHz. In this case, nonlinear parameters χ 21 T , χ 03 T and χ 03 D in Equation (16) are used for calculation.
As for simulations of IMD3, two excitation CW signals with frequencies f 1 and f 2 are applied to the quasi-3D FEM model of the 42-LT SAW resonator. The incident power levels of the two input tones are 15 dBm as well. In this case, the driving frequency f 1 is swept from 824 to 849 MHz, and f 2 is equal to f 1 45 MHz. Thus, IMD3 2 f 1 f 2 response can be found at 869 to 894 MHz and IMD3 2 f 1 + f 2 response will appear at 2427 to 2502 MHz. The simulated results of IMD3 signals fitted to each measured one are shown in Figure 4b,c, respectively. As shown, simulations of both the IMD3 signal with 2 f 1 f 2 and the IMD3 signal with 2 f 1 + f 2 exhibit decent agreement with the measured results in [16]. In these simulations, χ 30 T in Equation (17) is used to predict the IMD3 2 f 1 f 2 response, and χ 21 T , χ 03 T in Equation (18) are applied for the calculation of IMD3 2 f 1 + f 2 response.

4. Generation Mechanisms and Suppression of Nonlinearity

The proposed model can study the contributions of different generation mechanisms to nonlinearity in SAW resonators by setting the corresponding nonlinear parameters χ i j in Equations (9) and (10) to zero and non-zero. Thus, discussions on generation mechanisms of nonlinearity can be given based on simulation results of the above-mentioned SAW resonator on 42-LT substrate.

4.1. Generation Mechanisms of Third-Order Nonlinearity

Nonlinear effects in elasticity, dielectric, and electro-mechanical coupling such as piezoelectricity and electrostriction will contribute to nonlinear signals generated in SAW devices. To find the dominant contributions to third-order nonlinearity, nonlinear terms composed of different nonlinear coefficients χ i j combined with strain S and electric field E were investigated, respectively. Linear strain and electric field components of S y z and E z were used as they are expected to be the predominant components for SH-type SAW. Similarly, T x y and D z were selected and added as perturbations.
Figure 5a–c illustrate the separate contributions of different nonlinear terms to the simulation results shown in Figure 4a–c, respectively. Perturbations in D N are considered for the H3 and IMD3 2 f 1 + f 2 cases as their output frequency f f r , and that in T N are used for the IMD3 2 f 1 f 2 case as the output frequency is adjacent to the resonant frequency f r . For the H3 simulation shown in Figure 5a, nonlinear terms χ 21 T S y z 3 , χ 03 T S y z E z 2 and χ 03 D E z 3 in D z are of significant contributions, whereas the effect of χ 12 T S y z 2 E z is negligible. χ 21 T S y z 3 term contributes to a simple peak dependency of the H3 curve, and the notch at about 2.6 GHz in Figure 5a is mainly caused by the coupling effect of nonlinear terms χ 21 T S y z 3 and χ 03 D E z 3 . Thus, nonlinear piezoelectricity and dielectric represented by χ 21 T S y z 3 and χ 03 D E z 3 , respectively, are predominant for H3 generation. For the simulation of IMD3 2 f 1 f 2 response shown in Figure 5b, only χ 30 T S y z 3 in T x y matters and the contributions from the other three nonlinear terms are insignificant. It means nonlinear elasticity induced by χ 30 T S y z 3 is responsible for the IMD3 2 f 1 f 2 response in the SAW resonator. As for the simulation of the IMD3 2 f 1 + f 2 response shown in Figure 5c, χ 21 T S y z 3 , χ 03 T S y z E z 2 in D z are dominant in this case, and effects from the other two nonlinear terms can be neglected. Namely, nonlinear piezoelectricity caused by χ 21 T S y z 3 are also the major source of IMD3 2 f 1 + f 2 generation. It is concluded that the nonlinear effect of acoustic strain, i.e., S 3 terms in Equations (9) and (10), contributes to both H3 and IMD3 responses considerably. Therefore, nonlinear elasticity and piezoelectricity generated by acoustic strain are the dominant sources for third-order nonlinear responses of SAW devices.

4.2. Nonlinearity Suppression

In the last section, we point out that the third-order nonlinear signals of SAW devices are generated by nonlinear elasticity dominantly. Nakagawa et al. [26] also investigated the effect of different Ti layer thicknesses on H3 generation in SAW resonators with Al/Ti layered electrodes, which consist of metals with different elastic constants. Furthermore, they demonstrated that the H3 response of the SAW resonator decreases obviously with an increase in the thickness of the adhesive layer, Ti.
Similarly, we calculated the H3 signals of SAW resonators on 42-LT with three different Cu/Ti electrode structures, and the designed layer thicknesses are given in Table 1. As shown in Figure 6, the variation in H3 magnitude with Ti layer thickness is consistent with [26].
Furthermore, H3 calculation is performed for the SAW resonators with electrode materials of different Young’s modulus. Table 2 gives the parameters related to the elasticity of different electrode materials, where E is Young’s modulus, v is Poisson’s ratio, c 12 and c 44 are the two independent elastic constants for isotropic materials.
In simulations, thicknesses of different electrode materials are adjusted to keep the same resonant frequency f r   850 MHz. Figure 7 shows the simulated H3 results of the SAW resonators on 42-LT substrate with different electrode materials listed in Table 2. It is seen that the H3 curves in Figure 7 exhibit similar frequency dependence with those curves shown in the above Figure 4a. Then, Figure 8a compares the H3 peak values corresponding to different electrode materials in this calculation. These values are arranged in ascending order of Young’s modulus, and the values of c 12 and c 44 are organized in same order as well in Figure 8b. As shown, the H3 peak curve has the same trend as the c 12 curve, except for the values of c 44 being closer to or greater than c 12 . This means the larger one of the two independent elastic constants c 12 and c 44 possesses a more dominant contribution for H3 generation. It should be noted that the H3 peak can be suppressed by about 25 dBm by choosing proper electrode materials, which provides a vital insight into improving the linearity of SAW devices.

5. Conclusions

In this paper, we proposed a general FEM model for analyzing third-order nonlinear signals in RF SAW resonators based on perturbation theory. For validation, simulations of the H3 response and IMD3 response for a SAW resonator on 42-LT substrate are performed. The comparison of simulation results with measured results in the literature demonstrated the accuracy of the proposed model. The generation mechanisms of third-order nonlinearity in SAW resonators are discussed in detail and the dominant mechanisms are distinguished to give a guideline for further nonlinear suppression. In particular, the relation between elastic constants of electrode materials and H3 peak value is concluded. Additionally, due to the generality of the proposed nonlinear FEM model, this model could be extended to analyze the nonlinearity of SAW/BAW devices based on arbitrary structure configurations and materials for further study.

Author Contributions

Conceptualization, Q.Z.; methodology, Q.Z. and L.Q.; validation, B.L., S.Z. and Q.Z.; writing—original draft preparation, B.L.; writing—review and editing, B.L., Q.Z. and S.F.; supervision, Q.Z., S.F. and W.W.; project administration, X.Z.; funding acquisition, Q.Z. and X.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the National Natural Science Youth Foundation of China [Grant No. 11904233] and the National Natural Science Foundation of China [Grant No. 51772192].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Zhibin Xu, Shaoxu Dou and Qingshan Niu of Shoulder Electronics Limited. for their generous help and technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hagelauer, A.; Fattinger, G.; Ruppel, C.C.; Ueda, M.; Hashimoto, K.-Y.; Tag, A. Microwave acoustic wave devices: Recent advances on architectures, modeling, materials, and packaging. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2018, 66, 4548–4562. [Google Scholar] [CrossRef]
  2. Warder, P.; Link, A. Golden age for filter design: Innovative and proven approaches for acoustic filter, duplexer, and multiplexer design. IEEE Microw. Mag. 2015, 16, 60–72. [Google Scholar] [CrossRef]
  3. Liu, Y.; Cai, Y.; Zhang, Y.; Tovstopyat, A.; Liu, S.; Sun, C. Materials, Design, and Characteristics of Bulk Acoustic Wave Resonator: A Review. Micromachines 2020, 11, 630. [Google Scholar] [CrossRef] [PubMed]
  4. Hashimoto, K.-Y.; Omori, T.; Maruta, K.; Ahn, C.-J. Nonlinearity in RF front-end as a bottleneck in high speed mobile communications. In Proceedings of the 2017 International Symposium on Nonlinear Theory and Its Applications (NOLTA2017), Cancún, Mexico, 4–7 December 2017; pp. 193–196. [Google Scholar]
  5. Hashimoto, K.-Y. Theoretical considerations on influence of circuit impedance to IMD2 measurement of radio-frequency bulk acoustic wave duplexers. In Proceedings of the 2010 IEEE International Frequency Control Symposium, Newport Beach, CA, USA, 1–4 June 2010; pp. 146–150. [Google Scholar]
  6. Collado, C.; Rocas, E.; Mateu, J.; Padilla, A.; O’Callaghan, J.M. Nonlinear distributed model for bulk acoustic wave resonators. IEEE Trans. Microw. Theory Techn. 2009, 57, 3019–3029. [Google Scholar] [CrossRef] [Green Version]
  7. Rocas, E.; Collado, C.; Mateu, J.; Orloff, N.D.; Booth, J.C.; Aigner, R. Electro-thermo-mechanical model for bulk acoustic wave resonators. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2013, 60, 2389–2403. [Google Scholar] [CrossRef] [PubMed]
  8. Larson, J.D.; Bradley, P.D.; Wartenberg, S.; Ruby, R.C. Modified Butterworth-Van Dyke circuit for FBAR resonators and automated measurement system. In Proceedings of the 2000 IEEE Ultrasonics Symposium, San Juan, PR, USA, 22–25 October 2000; pp. 863–868. [Google Scholar]
  9. Shim, D.S.; Feld, D.A. A general nonlinear Mason model and its application to piezoelectric resonators. Int. J. RF Microw. Comput-Aid. Eng. 2011, 21, 486–495. [Google Scholar] [CrossRef]
  10. Brachtendorf, H.G.; Welsch, G.; Laur, R. Fast simulation of the steady-state of circuits by the harmonic balance technique. In Proceedings of the ISCAS’95-International Symposium on Circuits and Systems, Seattle, WA, USA, 30 April–3 May 1995; pp. 1388–1391. [Google Scholar]
  11. Gourary, M.; Ulyanov, S.; Zharov, M.; Rusakov, S.; Gullapalli, K.; Mulvaney, B. A robust and efficient oscillator analysis technique using harmonic balance. Comput. Methods Appl. Mech. Engrg. 2000, 181, 451–466. [Google Scholar] [CrossRef]
  12. Hashimoto, K.-Y.; Li, X.; Bao, J. Perturbation analysis of nonlinear signal generation in radio frequency bulk acoustic wave resonators. In Proceedings of the 2017 Joint Conference of the European Frequency and Time Forum and IEEE International Frequency Control Symposium (EFTF/IFCS), Besancon, France, 9–13 July 2017; pp. 538–541. [Google Scholar]
  13. Hashimoto, K.-Y.; Li, X.; Bao, J.; Qiu, L.; Omori, T. Perturbation Analysis of Nonlinearity in Radio Frequency Bulk Acoustic Wave Resonators Using Mass-Spring Model. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 2020, 67, 1479–1484. [Google Scholar] [CrossRef] [PubMed]
  14. Chen, L.; Solal, M.; Briot, J.; Hester, S.; Malocha, D.; Wahid, P. A nonlinear mason model for 3 rd order harmonic and intermodulation simulations of SAW duplexers. In Proceedings of the 2012 IEEE International Ultrasonics Symposium, Dresden, Germany, 7–10 October 2012; pp. 56–60. [Google Scholar]
  15. Chen, L.; Briot, J.; Girard, P.; Ledesma, C.; Solal, M.; Cheema, K.; Malocha, D.; Wahid, P. Third order nonlinear distortion of SAW duplexers in UMTS system. In Proceedings of the 2010 IEEE International Ultrasonics Symposium, San Diego, CA, USA, 11–14 October 2010; pp. 283–286. [Google Scholar]
  16. Nakagawa, R.; Suzuki, T.; Shimizu, H.; Kyoya, H.; Nako, K.; Hashimoto, K.-Y. Discussion about generation mechanisms of third-order nonlinear signals in surface acoustic wave resonators based on simulation. Jpn. J. Appl. Phys. 2016, 55, 07KD02. [Google Scholar] [CrossRef]
  17. Chauhan, V.; Mayer, M.; Ruile, W.; Ebner, T.; Bleyl, I.; Wagner, K.; Weigel, R.; Hagelauer, A. A P-Matrix Model for Third Order Electric and Acoustic Nonlinearities in TC-SAW Devices. In Proceedings of the 2018 IEEE International Ultrasonics Symposium (IUS), Kobe, Japan, 22–25 October 2018. [Google Scholar]
  18. Forster, T.; Mayer, M.; Chauhan, V.; Ebner, T.; Wagnery, K.-C.; Hagelauer, A. A general P-matrix model to calculate second-order nonlinearity in TC-SAW devices. In Proceedings of the 2020 IEEE International Ultrasonics Symposium (IUS), Las Vegas, NV, USA, 7–11 September 2020; pp. 1–4. [Google Scholar]
  19. Mayer, A.; Mayer, E.; Mayer, M.; Jäger, P.; Ruile, W.; Bleyl, I.; Wagner, K. Full 2D-FEM calculations of third-order intermodulations in SAW devices. In Proceedings of the 2016 IEEE International Ultrasonics Symposium (IUS), Tours, France, 18–21 September 2016; pp. 1–4. [Google Scholar]
  20. Guan, P.; Shi, R.; Yang, Y.; Qin, P.; Han, T. Mechanisms of Third-order Harmonic in TC-SAW Resonators Using a Nonlinear FEM Model. In Proceedings of the 2021 IEEE International Ultrasonics Symposium (IUS), Xi’an, China, 11–16 September 2021; pp. 1–4. [Google Scholar]
  21. Pang, X.; Yong, Y.K. Simulation of Nonlinear Resonance, Amplitude-Frequency, and Harmonic Generation Effects in SAW and BAW Devices. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2020, 67, 422–430. [Google Scholar] [CrossRef] [PubMed]
  22. Li, B.; Zhang, Q.; Liu, H.; Zhao, X.; Fu, S.; Wang, W. A General Fem model for Analysis of Nonlinearity in RF BAW Devices. In Proceedings of the 2020 15th Symposium on Piezoelectrcity, Acoustic Waves and Device Applications (SPAWDA), Zhengzhou, China, 16–19 April 2021; pp. 54–57. [Google Scholar]
  23. Bajenitchev, A. A numerical procedure for a non-linear elastic problem for incompressible material based on a perturbation method. Comput. Methods Appl. Mech. Engrg. 1996, 131, 31–39. [Google Scholar] [CrossRef]
  24. Auld, B.A. Thermodynamics of solids. In Acoustic Fields and Waves in Solids; Wiley: New York, NY, USA, 1973; Volume 1, pp. 275–278. [Google Scholar]
  25. Nakagawa, R.; Kyoya, H.; Shimizu, H.; Kihara, T.; Hashimoto, K.-Y. Study on generation mechanisms of second-order nonlinear signals in surface acoustic wave devices and their suppression. Jpn. J. Appl. Phys. 2015, 54, 07HD12. [Google Scholar] [CrossRef]
  26. Nakagawa, R.; Suzuki, T.; Shimizu, H.; Kyoya, H.; Hashimoto, K.-Y. Influence of electrode structure on generation of third-order nonlinearity in surface acoustic wave devices. Jpn. J. Appl. Phys. 2015, 54, 07HD11. [Google Scholar] [CrossRef]
Figure 1. (a) Schematic of a SAW resonator with a peripheral circuit (b) quasi-3D model for the SAW resonator used in simulation. (not to scale).
Figure 1. (a) Schematic of a SAW resonator with a peripheral circuit (b) quasi-3D model for the SAW resonator used in simulation. (not to scale).
Micromachines 13 01116 g001
Figure 2. Analysis procedures for simulation of nonlinear signals.
Figure 2. Analysis procedures for simulation of nonlinear signals.
Micromachines 13 01116 g002
Figure 3. Measured and simulated reflection coefficient S 11 of the SAW resonator.
Figure 3. Measured and simulated reflection coefficient S 11 of the SAW resonator.
Micromachines 13 01116 g003
Figure 4. Measured [16] and simulated nonlinear responses of the SAW resonator: (a) H3 response, (b) IMD3 2 f 1 f 2 response, (c) IMD3 2 f 1 + f 2 response.
Figure 4. Measured [16] and simulated nonlinear responses of the SAW resonator: (a) H3 response, (b) IMD3 2 f 1 f 2 response, (c) IMD3 2 f 1 + f 2 response.
Micromachines 13 01116 g004
Figure 5. Contributions of employed nonlinear terms to different nonlinear responses: (a) H3 response, (b) IMD3 2 f 1 f 2 response, (c) IMD3 2 f 1 + f 2 response.
Figure 5. Contributions of employed nonlinear terms to different nonlinear responses: (a) H3 response, (b) IMD3 2 f 1 f 2 response, (c) IMD3 2 f 1 + f 2 response.
Micromachines 13 01116 g005
Figure 6. Simulated H3 responses of the SAW resonators in Table 1.
Figure 6. Simulated H3 responses of the SAW resonators in Table 1.
Micromachines 13 01116 g006
Figure 7. Simulated H3 responses of the SAW resonators with different electrode materials.
Figure 7. Simulated H3 responses of the SAW resonators with different electrode materials.
Micromachines 13 01116 g007
Figure 8. (a) Simulated H3 peak magnitude of the SAW resonators with different electrode materials. (b) Elastic constants of different electrode materials arranged in ascending order of Young’s modulus.
Figure 8. (a) Simulated H3 peak magnitude of the SAW resonators with different electrode materials. (b) Elastic constants of different electrode materials arranged in ascending order of Young’s modulus.
Micromachines 13 01116 g008
Table 1. Layer Thicknesses of Different Electrodes.
Table 1. Layer Thicknesses of Different Electrodes.
Electrode TypeCu Thickness (nm)Ti Thickness (nm)
A136.515
B10087.6
C80127.3
Table 2. Parameters Related to the Elasticity of Different Electrode Materials.
Table 2. Parameters Related to the Elasticity of Different Electrode Materials.
Metal E [GPa] v c 12 [GPa] c 44 [GPa]
Al70.000.3351.0826.32
Ag83.000.3786.2230.29
Cu110.000.3595.0640.74
Ti115.700.3278.5343.79
Fe152.000.2770.2559.84
Pt168.000.38192.7560.87
Ni219.000.31136.3883.59
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, B.; Zhang, Q.; Zhao, X.; Zhi, S.; Qiu, L.; Fu, S.; Wang, W. A General FEM Model for Analysis of Third-Order Nonlinearity in RF Surface Acoustic Wave Devices Based on Perturbation Theory. Micromachines 2022, 13, 1116. https://doi.org/10.3390/mi13071116

AMA Style

Li B, Zhang Q, Zhao X, Zhi S, Qiu L, Fu S, Wang W. A General FEM Model for Analysis of Third-Order Nonlinearity in RF Surface Acoustic Wave Devices Based on Perturbation Theory. Micromachines. 2022; 13(7):1116. https://doi.org/10.3390/mi13071116

Chicago/Turabian Style

Li, Baichuan, Qiaozhen Zhang, Xiangyong Zhao, Shaotao Zhi, Luyan Qiu, Sulei Fu, and Weibiao Wang. 2022. "A General FEM Model for Analysis of Third-Order Nonlinearity in RF Surface Acoustic Wave Devices Based on Perturbation Theory" Micromachines 13, no. 7: 1116. https://doi.org/10.3390/mi13071116

APA Style

Li, B., Zhang, Q., Zhao, X., Zhi, S., Qiu, L., Fu, S., & Wang, W. (2022). A General FEM Model for Analysis of Third-Order Nonlinearity in RF Surface Acoustic Wave Devices Based on Perturbation Theory. Micromachines, 13(7), 1116. https://doi.org/10.3390/mi13071116

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