Published: 31 March 2018

Thickness shear vibration of a ZnO film structure covered with magnetic fluid

Jing Liu1
Xiangyang Li2
Weipeng Zhang3
1Piezoelectric Device Laboratory, School of Mechanical Engineering and Mechanics, Ningbo University, Ningbo, Zhejiang, 315211, China
2Magnetoelectronic Laboratory, Nanjing Normal University, Nanjing, 210023, China
3, 2School of Mechanical and Electrical Engineering, Ningbo Dahongying University, Ningbo, Zhejiang, 315175, China
Corresponding Author:
Jing Liu
Views 141
Reads 54
Downloads 1452

Abstract

This paper aims to reveal the effect of magnetic liquid on the resonant frequency shift and the output impedance. For this purpose, a thin film sensor was designed with a finite-thickness layer of magnetic fluid, and the thickness shear vibration of the film was analysed in details. The main findings are as follows. First, the relative resonance frequency is positively correlated with the intensity of magnetic field, and the viscosity coefficient, density and volume fraction of magnetic fluid. The effect of magnetic field on the shift of the relative resonant frequency depends on various parameters of magnetic liquid. With the increase of the volume fraction, density and viscosity coefficient, the resonant frequency plunges, but the output impedance peak grows in size. Therefore, the relative resonant frequency shift and the output impedance are sensitive to various parameters in the design. The research findings shed new light on the application of thin piezoelectric film in sensors.

1. Introduction

The last couple of decades has seen the successful deposition of an ultra-thin piezoelectric film of Aluminium-Nitride (AlN) or Zinc oxide (ZnO) on a silicon layer, laying the basis for thin-film bulk acoustic resonator (FBAR) working in the GHz frequency range [1-3]. The c-axis of the FBAR orients itself along the normal/in-plane direction of the film, deviates from the direction by a certain angle, or zigzags in a multilayer film [4-8]. It is possible to excite a pure shear vibration mode at a specific c-axis angle. The shear vibration can be detected with high sensitivity in the liquid environment, for the vibration is not damped due to the lack of compression. Hence, thickness shear mode sensors have been developed to detect compounds coated with surface bound receptive layers in gases or liquids. In particular, the FBAR has been utilized to make acoustic wave sensors for mass and fluid sensing [9-12].

For sensor applications, the resonator is usually capped with a thin chemical/physical sensitive film, e.g. the thin multilayer piezoelectric film. Over the years, much research has been done on magnetoelectric structures containing magnetic and ferroelectric/piezoelectric layers [13], because of their interesting properties for data storage [14]. As a magnetoelectric structure, the magnetic fluid hosts widely dispersed magnetic particles. In a certain sense, the fluid is “smart” because its viscoelasticity can be controlled by the intensity of the magnetic field. Thanks to the unique feature, the fluid has been increasingly applied in the sensor field. For instance, the magnetic fluid has been successfully implemented in vehicle acceleration sensor [15], angle sensor [16] and volumetric sensor [17]. Nevertheless, reports on the application of multilayer magnetic fluid structures few and scattered at present [18].

To make up for the gap, this paper proposes a thin film sensor coated with magnetic fluid layer (Figure 1). In general, a resonator consists of a thin layer of piezoelectric material and two electrodes on both sides. When an AC electric field is applied across the electrodes, acoustic waves come into being and start to propagate across the device. The propagation amplitude and phase depend on the interfacial reactions and the thin film parameters, including thickness, density, and complex shear modulus. Under the action of magnetic field, the magnetorheological properties change, and act on the piezoelectric thin film by sound wave. The variation in amplitude and phase, in turn, changes the electrical properties of the resonator. Based on piezoelectric elastic theory, viscoelastic mechanics and liquid mechanics, specifically, the author analysed the thickness shear vibration of ZnO thin film covered with a finite-thickness layer of magnetic liquid. The author investigated the propagation properties of bulk acoustic waves in the multilayer piezoelectric thin film with a tilted c-axis, seeking to derive the analytical method for bulk acoustic waves and aiming to disclose the effect of magnetic liquid on the shift of resonant frequency.

2. Governing equation

As shown in Fig. 1, the research object is the thickness shear vibration of a tilted c-axis ZnO film structure covered with magnetic fluid. The top and bottom surfaces of the film are positive and negative electrodes, respectively. The two layers were treated separately before applying the continuity and boundary conditions. The edge effect was neglected in view of the large length/thickness and width/thickness ratios. Then, a basic 1D theoretical model was constructed depending on the plate thickness and time. The 1D model depicts the basic behaviours of FBARs, and predicts the pure thickness resonant frequencies and modes. The model has nothing to do with the in-plane coordinates of the film, and varies only along the film thickness with one spatial variable along the normal direction of the layers. The model was constructed in the following steps.

Fig. 1A ZnO film structure with inclined c-axe covered with magnetic liquids. 2h is the thickness of the ZnO film. θ is the c-axis tilt angle of the ZnO film

A ZnO film structure with inclined c-axe covered with magnetic liquids.  2h is the thickness of the ZnO film. θ is the c-axis tilt angle of the ZnO film

Consider the following displacement fields that describe the thickness shear motion of interest, the displacement and electric potential fields are assumed to take the following forms [19]:

1
u1= u1(x3)exp(iωt) , u2=u3=0.

The differential equations of thickness shear motion and ZnO electrostatics [20] with a compact matrix notion[21] can be written as:

2
T5,3= ρu¨1,
3
D3,3=0,

where T5, u1 and D3 are thickness shear stress T13, horizontal displacement, and vertical electric displacement, respectively; ρ is the mass density of ZnO. In the following analysis, the mass and stiffness of the electrodes are neglected, the subscript comma denotes a partial derivative with respect to the coordinates, and a superimposed dot represents the derivative with respect to the time. Then, the strain-displacement and electric field-potential relations [20, 21] are:

4a
S5=u1,3,
4b
E3=-ϕ,3.

The nontrivial components of the stress and electric displacement [20, 21] are expressed as:

5
T5=c55S5-e35E3,
6
D3=e35u1,3+ε33E3.

Considering time-harmonic motions, the commonly used complex notation is adopted for this research. All the fields sharing the same time dependence with a common factor exp(iω t) are presented below. Substituting Eqs. (4a)and (5) into Eq. (2), and substituting Eqs. (4b) and (6) into Eq. (3), we have:

7
c-55u1,33+ρω2u1=0,
8
e35u1,33-ε33ϕ,33=0,

where c-55=c55+e352/ε33.

The general solution to Eqs. (7), (8) can be written as:

9
u1x3=Acosξx3+Bsinξx3,
10
ϕx3=e35ε33u1+Cx3+D,

where ξ=ωρ/c-55 is the wave number for the thickness shear mode; A, B, C and D are the unknown amplitudes determined by the boundary conditions.

As mentioned above, the magnetic fluid hosts widely dispersed magnetic particles, and its viscoelasticity can be controlled by the intensity of the magnetic field. The pressure gradient is negligible because only the shear deformation occurs. Hence, the governing equations for the liquid are simplified as:

11
μL2νL=ρLν˙L,
12
T5L=μLνL,3,

where μL, νL, ρL and TL are the dynamic viscosity coefficient, velocity, density and stress of the magnetic liquid, respectively.

The magnetic field was considered to be uniform, and the external magnetic field vector was regarded as perpendicular to the vorticity. Then, the dynamic viscosity coefficient can be expressed as [18]:

13
μL=η01+32ϕψLψ2+ψLψ,

where η0 is the magnetic fluid viscosity coefficient without external magnetic field; ϕ is the fraction of solid; ϕ and Lψ is the Langevin function. In addition, η0=ηc[1-52ϕ~+1.55ϕ~2]-1 [22], where ϕ~ is the volume fraction of the suspended material and ηc is the kinetic viscosity coefficient of the current carrier. The parameter ψ is calculated by ψ=μ0mH/k0T, where T is the ambient temperature, m is the magnetic moment of single solid particle, k0is the Boltzmann constant, μ0 is the vacuum permeability, and H is the external magnetic field. Based on the value of ψ, it is easy to obtain the numerical value of Langevin function Lψ=cosψ-1/ψ.

The general solutions to Eqs. (10), (11) can be written as:

14
νL=ALexp-λ(x3-h)+iωt+BLexpλ(x3-h)+iωt,
15
T5L=-μLλALexp-λx3-h+iωt,

where:

16
λ=iωρLμ=1+i2ρωμ.

Since x3+, νL is finite, the second term of Eq. (14) must be removed if BL= 0.

3. Boundary conditions and solutions

The continuity conditions at the interface between the thin film and the magnetic fluid are:

17
T5h=T5Lh, u˙1(h)=νL(h).

The stress boundary conditions at the bottom of the film is:

18
T5(-h)=0.

When the driving electrode is shorted to the ground electrode, the electric potential boundary conditions are:

19
ϕ(-h)=0,
20
ϕ(h)=0.

Substituting Eqs. (9), (10) and (14), (15) into Eqs. (17)-(20), it is possible to derive the algebraic equations of the integration constants:

21
c-55ξsin(ξh)A+c-55ξcos(ξh)B+e35C=0,
22
e35ε33cos(ξh)A-e35ε33sin(ξh)B-hC+D=0,
23
e35ε33cos(ξh)A+e35ε33sin(ξh)B+hC+D=0,
24
-c-55ξsinξhA+c-55ξcosξhB+e35C=-μLλAL,
25
iωcosξhA+iωsinξhB=AL.

To find the nontrivial solutions to the integration constants, the coefficient matrix of these linear algebraic equations must be equal to zero. This yields the frequency equation below, which determines the resonant frequencies of the plate carrying the magnetic fluid:

26
iωμλh-2e352ε33tan2(ξh)+iωμλe352c-55ξε33+2hξc-55tan(ξh)-iωμλh=0.

For forced vibration, the boundary conditions are:

27
e35ε33cosξhA-e35ε33sinξhB-hC+D=-ϕ0,
28
e35ε33cosξhA+e35ε33sinξhB+hC+D=ϕ0.

The current of the FBAR is expressed by I=dQ/dt=-iωε33CS. The output impendence can be expressed as:

29
Z=2ϕ0I=2ϕ0iωε33CS
=1iωε33S2iωμLλe352tanξh-4e352c-55ξtan2ξh2c-552ξ2ε33tanξh+iωμLλε33c-55ξtan2ξh-1+2h,

where S is the electrode area; i is the imaginary unit. The numerical solutions of Eqs. (26) and (29) are obtained with Matlab by selecting appropriate parameters.

4. Numerical results and discussion

In our numerical example, the thickness of the ZnO film is 2h= 3 μm and the electrode area is S= 9×10-8m2. These are the typical values for thin film AlN and ZnO resonators. Whereas S/2h=100, the length/thickness ratio is sufficiently large to make pure thickness shear vibration with negligible edge effect. The c-axis inclination is θ= 43° [23], whose e35 is large and e33 is virtually zero [5]. Hence, thickness shear vibration can be excited easily without stretching the thickness.

The constants in the x,y,z system of ZnO are derived from a widely accepted reference for material data [24]. In the x,y,z system, the material constants are as follows:

30
ε0=8.75×10-12 F/m, ρ=5700 kg/m3,
cpq=210121105000121210105000105105211000000430000004300000044.5×109 N/m2,
eip=0000-0.480000-0.4800-0.57-0.571.32000 C/m2,
[εij]=ε07.610007.610008.85F/m.

The material constants of ZnO with c-axis tilted 43° were obtained by tensor transformation. The constants in our calculation are listed in Table 1. Some of the geometric parameters of magnetic fluid fixed for all calculations are presented in Table 2. The other geometric parameters are specified separately in individual figures. Based on Eqs. (26), (29) and material constants, the resonator spectrum with magnetic fluids was simulated.

Table 1Material coefficients of the piezoelectric layer Zinc oxide

c55 (109 Pa)
e35 (C/m2)
ε33 (10-11C2/Nm2)
ρ (103 kg/m3)
Zinc oxide
52.703
0.7123
6.9949
5.7

Table 2Material coefficients of the magnetic fluid

T (K)
m (A/m)
k0
μ0 (N/A-2)
Magnetic fluid
300
9.273×10-24
1.380×10-23
π×4×10-7

For free vibration, the frequency shift curves are shown in Figs. 2-5. Here into, Fig. 2 (ρL= 800 kg/m3) and Fig. 3 (ϕ~= 0.45) depict the shift of the resonant frequency versus viscosity coefficient, when the external magnetic field is fixed at H= 1×104 A/m. There are three different curves of volume fraction ϕ~ in Fig. 2 and three different curves of density ρL in Fig. 3. The relative resonant frequency increases with the viscosity coefficient. In other words, the frequent shift is more prominent at greater values of ϕ~ and ρL. Moreover, it can be seen that the three volume fraction curves are more separated than the three density curves.

Fig. 4 (ρL= 800 kg/m3) and Fig. 5 (ϕ~= 0.45) illustrate the shift of the relative resonant frequency versus uniform magnetic field at the fixed viscosity coefficient of ηc= 10 N·s/m2. The three density curves (ϕ~= 0.2; 0.45; 0.6) in Fig. 4 are almost parallel to the x-axis, and are virtually unchanged with the increase in the uniform magnetic field. This means the magnetic field has a weak effect on the acoustic waves of the ZnO piezoelectric film. In Fig. 5, the relative resonant frequency increases with the uniform magnetic field, that is, the frequency shift is positively correlated with density. In comparison, the effect of magnetic field is more pronounced in Fig. 5.

Fig. 2Relative resonant frequency versus viscous coefficient with different volume fraction

Relative resonant frequency versus viscous coefficient with different volume fraction

Fig. 3Relative resonant frequency versus viscous coefficient with different density

Relative resonant frequency versus viscous coefficient with different density

Fig. 4Relative resonant frequency versus uniform magnetic field with different volume fraction

Relative resonant frequency versus uniform magnetic field with different volume fraction

Fig. 5Relative resonant frequency versus uniform magnetic field with different density

Relative resonant frequency versus uniform magnetic field with different density

For forced vibration, Fig. 6 shows the output impedance versus the driving frequency at ρL= 800 kg/m3, ηc= 10 N·s/m2 and H= 1×104 A/m, and the three values of ϕ~; Fig. 7 presents the output impedance versus the driving frequency at ηc= 10 N·s/m2, ϕ~= 0.45 and H= 1×104 A/m and the three values of ρL; Fig. 8 illustrates the output impedance versus the driving frequency at ρL= 800 kg/m3, ϕ~= 0.45 and H= 1×104 A/m and the three values of ηc.

In Fig. 6, there are three resonant peaks of the curves, and the curve (ϕ~= 0.6) differs from the other two. It is special that the curve (ϕ~= 0.6) stays at frequency about 0.8 GHz. The larger the volume fraction, the stronger the interface interaction and the more resonant peaks. In the given frequency range, there is a resonant peak much higher than the other peaks, and the output impedance reaches the maximum value in Figs. 6-10. The resonant frequencies of ZnO FBARs decrease, because the resonator is being loaded with magnetic fluid. With the increase of the volume fraction ϕ~, the density ρLand the ηc, the resonant frequency plunges, but the resonant peaks grow in size. Due to the small damping of the material, the resonance peaks are sharp, and the output impedance is significant near resonance. In Fig. 10, it seems that there are more than three resonant peaks of curve (ϕ~= 0.6). As shown inFig. 6 and Fig. 10, the resonant peakschange alike. The volume fraction of magnetic particles apparently has a stronger impact on the acoustic waves of the ZnO piezoelectric film than the other three factors.

Fig. 6Output impedance spectrum with different volume fraction

Output impedance spectrum  with different volume fraction

Fig. 7Output impedance spectrum with different density

Output impedance spectrum  with different density

Fig. 8Output impedance spectrum with different viscous coefficient

Output impedance spectrum  with different viscous coefficient

Fig. 9Output impedance spectrum with different magnetic field intensity coupling

Output impedance spectrum with different magnetic field intensity coupling

Fig. 10Output impedance spectrum with different coupling of the volume fraction, density, viscous coefficient, magnetic field intensity coupling

Output impedance spectrum with different coupling of the volume fraction,  density, viscous coefficient, magnetic field intensity coupling

5. Conclusions

This paper analyses the thickness shear vibration of ZnO film covered with a finite-thickness layer of magnetic liquid. The analysis reveals that relative resonant frequency and the output impedance are sensitive to various design parameters. Specifically, the relative resonant frequency is positively correlated with the intensity of magnetic field, and the viscosity coefficient, density and volume fraction of magnetic fluid. Under the same viscosity coefficient and density, the magnetic field has a weak effect on the shift of the relative resonant frequency. By contrast, the effect of magnetic field is obvious when the density varies and the other two parameters are constant. For forced vibration, a resonant peak is much higher than other peaks, and the output impedance reaches the maximum in the given frequency range. Moreover, with the increase of the volume fraction ϕ~, the density ρL and the kinetic viscosity coefficient ηc, the resonant frequency plunges, but the resonant peaks grow in size.

References

  • Lakin K. M., Wang J. S. Acoustic wave bulk composite resonator. Applied Physics Letters, Vol. 38, Issue 3, 1981, p. 125-127.
  • Nakamura K., Sasaki H., Shimizu H. ZnO/SiO2 diaphragm composite resonator on a silicon wafer. Electronics Letters, Vol. 17, Issue 14, 1981, p. 507-509.
  • Ruby R., Bradley P., Larson J. D., Oshmyansky Y. PCS 1900 MHz duplexer using thin film bulk acoustic resonators (FBARs). Electronics Letters, Vol. 35, Issue 17, 1999, p. 794-795.
  • Yanagitani T., Mishima N., Matsukawa M., Yatanabe Y. Electromechanical coupling coefficient k15 of polycrystalline ZnO films with the c-axes lie in the substrate plane. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, Vol. 54, Issue 4, 2007, p. 701-704.
  • Du J. K., Xian K., Wang J., Yang J. S. Thickness vibration of piezoelectric plates of 6 mm crystals with tilted six-fold axis and two-layered thick electrodes. Ultrasonics, Vol. 49, Issue 2, 2009, p. 149-152.
  • Yanagitani T., Morisato N., Takayanagi S., Matsukawa M., Watanabe Y. C-axis zig-zag ZnO film ultrasonic transducers for designing longitudinal and shear wave resonant frequencies and modes. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, Vol. 58, Issue 5, 2011, p. 1062-1068.
  • Zhang H. F. Analysis of thickness vibrations of C-axis inclined zig-zag multi-layered zinc oxide thin film resonators. Integrated Ferroelectrics, Vol. 113, Issue 1, 2010, p. 95-108.
  • Zhang H., Bao Y. Sensitivity analysis of multi-layered C-axis inclined zigzag zinc oxide thin-film resonators as viscosity sensors. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, Vol. 61, Issue 3, 2014, p. 525-534.
  • Satoh Y., Nishihara T., Yokoyama T., Ueda M., Miyashita T. Development of piezoelectric thin film resonator and its impact on future wireless communication systems. Japanese Journal of Applied Physics, Vol. 44, Issue 5, 2005, p. 2883-2894.
  • Link M., Schreiter M., Weber J., Primig R., Pitzer D., Gabl R. Solidly mounted ZnO shear mode film bulk acoustic wave resonators for sensing applications in liquids. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, Vol. 53, Issue 2, 2006, p. 492-496.
  • Fu Y. Q., Luo J. K., Du X. Y., Flewitt A. F., Li Y., Markx G. H., Walton A. J., Milne W. I. Recent developments on ZnO films for acoustic wave based bio-sensing and microfluidic applications: a review. Sensors and Actuators B Chemical, Vol. 143, Issue 2, 2010, p. 606-619.
  • Liu J., Du J. K., Wang J., Yang J. S. Long thickness-extensional waves in thin film bulk acoustic wave filters affected by interdigital electrodes. Ultrasonics, Vol. 75, 2017, p. 226-232.
  • Li X. Y., Liu J., Zhang N. Magnetoelectric coupling by acoustic wave guide. Journal of Applied Physics, Vol. 119, Issue 13, 2016, p. 441-282.
  • Nan W., Bichurin M. I., Dong S., Viehland D., Srinivasan D. Multiferroic magnetoelectric composites: Historical perspective, status, and future directions. Journal of Applied Physics, Vol. 103, Issue 3, 2008, p. 31101.
  • Kato T., Hatakeyama M. Acceleration Sensor. US, 5248861, 1993.
  • Donatini F., Monin J., Noyel G. Investigation of an inclination sensor using magnetic fluid for angle measurement up to 90 degrees. Measurement Science and Technology, Vol. 6, Issue 1, 1995, p. 1-3.
  • Stanci A., Iusan V., Buioca C. D. Magnetofluidic sensor for volume measurement. Sensors and Actuators A Physical, Vol. 84, Issue 3, 2000, p. 246-249.
  • Ye J. J., Liu G. X., Liu H. Analysis of SH-SAW propagation characteristics in piezoelectric layered structure based on viscous loading effect of magnetic fluid. Science Technology and Engineering, Vol. 11, Issue 2, 2011, p. 4614-4617.
  • Tiersten H. F. Thickness vibrations of piezoelectric plates. Journal of the Acoustical Society of America, Vol. 35, Issue 1, 1963, p. 53-2.
  • Yang J. S. A review of analyses related to vibrations of rotating piezoelectric bodies and gyroscopes. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, Vol. 52, Issue 5, 2005, p. 698-706.
  • Yang J. Analysis of Piezoelectric Devices. World Scientific, 2006.
  • Popplewell J., Rosensweig R. E., Siller J. K. Magnetorheology of ferrofluid composites. Journal of Magnetism and Magnetic Materials, Vol. 149, Issue 149, 1995, p. 53-56.
  • Qin L., Chen Q., Cheng H., Wang Q. M. Analytical study of dual-mode thin film bulk acoustic resonators (FBARs) based on ZnO and AlN films with tilted c-axis orientation. IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, Vol. 57, Issue 8, 2010, p. 1840-1853.
  • Auld B. A. Acoustic Fields and Waves in Solids. John Willey and Sons, New York, 1973, p. 357-382.

Cited by

Frequency shifts in a laterally finite piezoelectric resonator loaded with a viscous liquid layer of finite size
Xueli Zhao | Feng Zhu | Peng Li | Zhenghua Qian | Iren Kuznetsova | Andrey Smirnov
(2023)

About this article

Received
26 October 2017
Accepted
12 January 2018
Published
31 March 2018
SUBJECTS
Mechanical vibrations and applications
Erratum
The title of affiliation was wrong in the paper finally approved (after the acceptance) by the authors. For more information read Editor's Note.
Keywords
piezoelectric film
resonator
vibration
magnetic liquid
Acknowledgements

This research was supported by the Ningbo Natural Science Foundation in Zhejiang province, China under Grant No. 2017A610105, and the Public Projects of Zhejiang Province (Nos. 2015C31150 and LGN18F020003).