SIMULATION OF METAL LAYERS THICKNESS INFLUENCE ON PHASE VELOCITIES OF ACOUSTIC WAVES IN THE PIEZOELECTRIC PLATES

Texto integral

Introduction. At the present time the properties of the elastic waves propagating in the thin plate which thick- ness has the same order as longitudinal and transverse waves lengths in the considered material are studied quite explicitly [1]. There is a large number of pilot and theo- retical studies of characteristics of various Lamb elastic waves and SH waves modes in plates from isotropic, pie- zoelectric and composite materials [2-4]. Metal layer application (metallization) on piezoelec- tric plates surfaces leading to the electric short circuit of a surface, as a rule, decreases the velocity of acoustic The executed simulation results can be used when de- veloping various acousto-electronic devices including those used as components of the rocket and space technology electronic base. Theoretical bases of elastic waves propagation in layered piezoelectric medium. Let us consider the propagation of small amplitude elastic waves in unper- turbed by external influences piezoelectric crystal. The motion equation, the electrostatics equation and the state of piezoelectric medium equations have the following form [14]: waves propagation [5; 6]. Acoustic waves in piezoelectric plates are a perspec- ρ0UɺA = τ AB,B ; DM ,M = 0; tive object for the development of microwave resonators τ = cE η - e E ; (1) AB ABCD CD MAB M and sensors of the next generation. Due to high interest in this issue there appeared many pilot and theoretical stud- ies of acousto-electronic devices on the basis of the pie- zoelectric plates operating at various modes of elastic waves [7-10]. As a rule, for of elastic waves in devices D = ε E + e η , η M MN N MAB AB c ε - η where ρ0 - crystal density in an undeformed state; UA - vector of dynamic elastic displacements; τAB - tensor of thermodynamic stresses; DM - electric induction vector; electrodes (interdigitated transducers, IDT) placed on a plate surface are used. ηCD - tensor of small deformations; E ABCD , eMAB , MN In devices created on the basis of a single-layer or a multilayer piezoelectric structure, the thickness of elec- trodes can be comparable to the acoustic line plate thick- ness. At the same time the precious metals, such as gold or platinum which are often used in sensors for various biochemical applications, have essential values of an acoustic impedance [11]. Therefore, metal layers thickness and mass cannot be neglected. Thus, considering the impact of metal layers on the change of dispersion characteristics of various modes of an elastic wave in the piezoelectric plate is important for the improvement of devices characteristics on similar structures [12]. In this work the computer simulation of load influence in the form of two metal layers on dispersion characteris- tics of Lamb elastic wave modes in zincum oxide plates (Me/ZnO/Me structure) and aluminum nitride (Me/AlN/Me structure) was executed. These piezoelec- tric materials have great values of the electromechanical coupling coefficient and considerable values of phase elastic, piezoelectric and dielectric constants of the second order. For elastic displacements and electric potential in the form of plane monochrome waves of small amplitude the set of equations (1) is recorded as Green-Christoffel equation which has to be solved for each layered struc- ture medium [15]. Let axis X3 be directed along an outward normal to the layer surface in a working orthogonal system of coordi- nates and axis X1 coincide with the wave propagation direction (fig. 1, а). The elastic wave propagation in such system has to meet the corresponding boundary condi- tions. Boundary conditions, in particular, for the three- layered structure „metal/piezoelectric/metal“ are: equality to zero of normal component of a stress tensor on the „metal/vacuum“ interface; equality of normal component of a stress tensor, equality of displacement vectors and equality to zero of the electric potential wave on the „metal/piezoelectric“ interface [16]: velocities for bulk and surface acoustic waves. Thanks to these properties the considered materials are actively ap- τ(1) 3 j x3 =d1 = 0; plied when developing various acousto-electronic devices. Aluminum and molybdenum in the form of the thin (1) (2) τ = τ 3 j 3 j ; φ(2) = 0 x3 =h x3 =h . . ; U (1) = U ( 2) ; x3 =h . . (2) metalizing film were used as materials for metal layers, τ(2) = τ(3) ; φ(2) = 0 ; U (2) = U (3) ; that is the metals which are often applied at electrodes production [13]. The choice of two metals with various 3 j τ(3) 3 j x3 =d2 = 0. x3 =d2 x3 =d2 acoustic impedances also allows to analyse how the ma- terial layer acoustic impedance influences the change of phase velocity of various elastic wave modes in an acous- tic line. 3 j x3 =0 Here d1, d2 and h - thickness of metal overlayer, sublayer and piezoelectric layer respectively. In this work the d1 = d2 condition is met, i.e. metal overlayer and sublayer thicknesses are equal. Substituting the elastic wave equation in the form of linear combinations of partial waves in boundary condi- tions: similar dependences for the layered structure Me/AlN/Me with orientation of the layer [100] (001) are presented in fig. 1, c. As a metal layer (Me) Al and Mo having various values of acoustic impedance Z = ρv were used. The acoustic impedances calculated for bulk longitudinal waves of all considered materials relate as U = åC(m)α(n)exp éi (k x + k (n) x - ωt )ù , Z < Z < Z < Z . The values of the material coni n i ë 1 1 3 3 û Al AlN ZnO Mo n φ = åC (m)α(n)exp éi (k x + k (n)x - ωt )ù , (3) stants for aluminum and molybdenum are taken in [13], ZnO - in [20], AlN - in [21]. The change range of phase 4 4 ë n 1 1 3 3 û velocities of the elastic wave running modes in the strucwhere the superscript n is equal to the number of partial waves in the corresponding layer, we will recieve a linear equation system - boundary conditions. Equality to zero of boundary conditions matrix continuant (2) which di- mension in this case is 20×20 elements allows to calculate phase velocities of the symmetric and antisymmetric modes of the Lamb wave and SH-wave. Variations of boundary conditions (2) define all types of elastic waves propagating in the three-layered plate. When using isotropic metals and piezoelectric materials of a symmetry group 6 mm for the considered layered structures Me/ZnO/Me and Me/AlN/Me, a set of equa- tions describing boundary conditions (2) is divided into two independent parts: a set of equations for the Lamb wave modes and a set of equations for SH-modes of the elastic wave. The computing experiment was carried out with the use of the software developed by authors. Analysis of mass load influence on waves velocities in three-layered plates. There are several determinations of piezoelectric structure sensitivity parameter to mass loading [17; 18]. In particular, S mass sensitivity of the multilayer resonator can be defined as the displacement of the relative resonator frequency normalized to the mass density surface [18]. However the disadvantage of this approach is that at great values of frequency f the changes of phase velocity at increase in thickness of a metal layer are smoothed. In this work the following formula was used for mass sensitivity of elastic wave determination: dr ç v ÷ S = 1 æ Dv ö , (4) è ø where Dv = v - vmet - the change of an elastic wave phase velocity when metallizating with a bulk density r ; d - the total thickness of metal layers. Velocities displace- ment is calculated concerning the elastic wave v velocity when metallizating with infinitely thin layer which does not change mechanical boundary conditions [19]. Let us note that this condition, as the computer ex- periment revealed, is carried out at a layer thickness less than 1 nanometer. Also the condition (4) levels the influ- ence of the piezoelectric effect and, therefore, parameter S will not depend on the piezoelectric layer type. The schematic representation of the studied layered structures is presented in fig. 1, а. Phase velocities disper- sion dependences of the Lamb wave fundamental modes and fundamental and first modes of SH-wave on value of the h´f parameter (piezoelectric thickness ´ frequency) for the layered structure of Me/ZnO/Me type with orienta- tion of the layer [100] (001) are presented in fig. 1, b, ture Me/ZnO/Me at metallization with infinitely thin lay- ers is ranging from values of a longitudinal wave phase velocity in ZnO crystal (v = 6084.14 m/s) or from zero (A0 mode) to the phase velocity of the Rayleigh wave in ZnO at the metallized surfaces ( v = 2677.6 m/s). It should be noted that due to „switching off“ of the piezoelectric effect at metallization of plate both surfaces the phase velocity of SH0 mode coincides with the velocity of a fast shear wave in ZnO crystal. In Me/AlN/Me structure phase velocities of elastic wave running modes also change from the longitudinal wave phase velocity values in AlN film (v = 10287.28 m/s) or from zero (mode A0) to the phase velocity of the Rayleigh wave in AlN at metallized surfaces (v = 5485 m/s). Unlike Me/ZnO/Me structure with the maximal value of the parameter h´f = 10000 m/s for calculation of Me/AlN/Me structure the range of change of the pa- rameter h´f to 20000 m/s was increased due to the fact that phase velocities values of all elastic wave modes significantly increased. Calculation of elastic wave phase velocities modes change was made at the ratio of a metal layer thickness to piezoelectric material thickness d/h = 0.001; 0.005; 0.01; 0.05; 0.1. Solid lines in fig. 1, b, c designated phase ve- locities modes of the Lamb elastic wave and SH-modes at metallization of a piezoelectric layer above and below by two infinitely metal layers which are not violating me- chanical boundary conditions are shown with solid lines in fig.1, b, c. Dashed lines in fig. 1, b, c correspond to the piezoelectric material layer d/h = 0.1. In absolute values the presented values of phase velocities can be received at piezoelectric material thickness layers h = 1 mkm and each of metal layers d = 0.1 mkm at frequency change from 1 MHz to 10 GHz (20 GHz for AlN). Metal use in a multilayer plate can considerably affect wave elastic energy redistribution due to acoustic proper- ties of layers and also an interference between elastic wave modes falling and reflected from borders of layers. For example, in Al/AlN/Al structure at values h´f increase and metal layers of maximal thickness (d/h = 0.1) use phase velocities of elastic wave modes tend to the value of phase velocity of the bulk shear wave in aluminum (fig. 1, c). The similar situation occurs also when used other metals in multilayer structure. At the same time it is noted that in Me/AlN/Me structure elastic wave modes velocities at identical type of metal and the parameter value h´f decrease more significantly than in Me/ZnO/Me structure. When comparing identical structures with metal layers different thickness it was revealed that the increase in metal layer thickness, as a rule, reduces the velocities of propagation of the first and the higher orders acoustic waves. In fig. 2 the graph of sensitivity S calculated by formula (4) for fundamental A0, S0 modes, the first modes A1, S1 of the Lamb wave and also fundamental and the two first modes of SH-wave depending on the parameter h´f in Al/ZnO/Al structure providing d/h = 0.001 is presented. It is possible to note that on modes dispersion dependences S0, A1, S1 of the Lamb wave the expressed S maxima are observed. The maximal value of S = 251 sq.m/kg for S0 mode at h´f = 2500 m/s , for A1 mode S = 248 sq.m/kg at h´f = 4700 m / s, for S1 mode S = 249 sq.m/kg at h´f = 7000 m / s (fig. 2). Fig. 1. Dispersive dependence of phase velocities on the product h´f for the layered structure Al/Piezoelectric/Al: a - the scheme of the layered structure; b - fundamental and first modes of Lamb and SH-wave in the structure Al/ZnO/Al; c - fundamental and first modes of Lamb and SH-wave in the structure Al/AlN/Al. Solid lines - mode velocities for infinitely thin metal layers, dash lines - mode velocities for metal layers thickness d/h = 0.1 Рис. 1. Дисперсионные зависимости фазовых скоростей от произведения h´f для слоистой структуры Al/пьезоэлектрик/Al: а - схема слоистой структуры; б - фундаментальные и первые моды волны Лэмба и SH-волны в структуре Al/ZnO/Al; в - фундаментальные и первые моды волны Лэмба и SH-волны в структуре Al/AlN/Al; сплошные линии - скорости мод при бесконечно тонких слоях металла, точечные линии - скорости при относительной толщине металла и слоя d/h = 0,1 Fig. 2. Dispersive dependence of the sensitivity S for Lamb and SH-modes of an elastic wave on the product h´f for the layered structure Al/ZnO/Al. Ratio of the metal layers and piezoelectric d/h = 0.001 Рис. 2. Дисперсионные зависимости чувствительности S мод упругой волны Лэмба и SH-волны от произведения h´f для слоистой структуры Al/ZnO/Al при отношении толщин металла и пьезоэлектрика d/h = 0,001 Therefore, at metallization of a plate from above and from below by two metal thin layers the modes of an elas- tic Lamb wave above the zero order (except for mode A0) have the sensitivity maxima occurring practically through identical intervals of the parameter h´f values. S value maxima are revealed only for the Lamb wave modes A1 and S1. They are explained by the condition of „resonance“ reflection of these waves from a piezoelectric boundary with metal layers realization, what leads to the emergence of local maxima in values of S (fig. 2). The graphs of sensitivity S fundamental mode of an antisym- metric wave A0 and SH-mode with shear-horizontal and horizontal polarization on the studied frequency range have no peculiarities. Only the apparent dependence of the module S decrease at the increase of the parameter h´f is characteristic for them. It is possible to note that the decrease (on module) of S values is connected with grad- ual mode А0 velocity increase at the values of the parame- ter h´f increase. At the same time the velocity of this mode remains the minimal of all considered elastic wave modes velocities. The decrease in S values for transversal and horizontal polarization modes is connected with the same gradual decrease of these modes velocities at the increase of the parameter h´f values. Further in this work only features of acting of S sensitivity curves for the Lamb elastic wave mode having maxima on the consid- ered range of the parameter h´f values are being de- scribed. In fig. 3 the graphs of S sensitivity of S0 fundamental mode and the first modes A1 и S1 of the Lamb wave de- pending on the metal layer material, the ratio of d/h and h´f parameter in Me/ZnO/Me structure are presented. The maximal values of S parameter for S0 mode are reached at h´f = 2400-2500 m / s and the ratio of d/h = 0.001 where condition h = 2λ/3 is satisfied, i. e. a piezoelectric layer thickness is equal to two thirds of the elastic wave S0 length (fig. 3, а, d). When using Al and Mo layers the value of parameter S = 251 sq.m/kg and S = 246 sq.m/kg, respectively. Essential decrease of S sensitivity for the considered elastic wave modes is particularly evident at d/h > 0.05. This decrease of S values at d/h increase is especially observed when using Mo - the metal with high acoustic impedance value. Numerical values of S maxima for both types of metal differ only at the considerable thicknesses of layers, for example, in Al/ZnO/Al structure at d/h = 0.1 the value of sensitivity of S = 200 sq.m/kg at h´f = 2300 m / s, and in Mo/ZnO/Mo structure S = 155 sq.m/kg at h´f = 2100 m / s (tab. 1). The maximal values of S parameter for A1 mode of the Lamb elastic wave are reached at h´f = 4700 m / s and d/h = 0.001 where the condition h = λ is satisfied, i. e. thickness of a piezoelectric layer is equal to A1 elastic wave length (fig. 3, b, e). The maximal values of S parameter for S1 mode of the Lamb elastic wave are reached at h´f = 7000 m / s and d/h = 0,001 where the condition h = 3λ/2 is satisfied, i. e. the piezoelectric layer thickness is equal to one and a half length of S1 elastic wave (fig. 3, c, f). At the same time it was also revealed that application of Mo layers leads to more considerable decrease of the S maximal sensitivity value than when using Al. According to the graphs in fig. 3 and to the data in tab. 1 it is noted that the increase in a metal layer thick- ness (d/h) always leads to wave frequency displacement in the S maximum to more low-frequency area. In fig. 4 the graphs of S sensitivity of the Lamb wave fundamental and first modes depending on the metal layer type, the ratio of d/h and the h´f parameter in Me/AlN/Me structure are shown. The maximal values of S parameter for S0 mode are reached at h´f = 5200 m / s and d/h = 0.001 on the condi- tion when h = 2λ/3, as well as in the previous structure for the similar mode (fig. 4, а, d). When using Al and Mo layers at S = 502 sq.m/kg and S = 493 sq.m/kg, respec- tively. Essential decrease of S sensitivity for the consid- ered elastic wave modes is shown at d/h > 0.05, in par- ticular when using Mo. S maxima numerical values for both types of metal differ only at the considerable thick- nesses of layers, for example, in Al/AlN/Al structure at d/h = 0.1 the sensitivity value of S = 432 sq.m/kg at h´f = 4800 m / s, and in Mo/AlN/Mo structure S = 239 sq.m/kg at h´f = 4200 m / s (tab. 2). The maximal values of S parameter for A1 mode of the Lamb elastic wave are reached at h´f = 9600 m / s and the ratio of d/h = 0.001 where h = λ condition is satisfied (fig. 4, b, e). The maximal values of S parameter for A1 mode of the Lamb elastic wave are reached at h´f = 13800-14000 m / s and d/h = 0,001 where h = 3λ/2 condition is satisfied (fig. 4, c, f). The values of S parameter in maxima exceed similar values by 1.5-2 times for Me/ZnO/Me structure, at the same time the increase in thickness of a metal layer also leads to the wave frequency displacement in maxima of S parameter values. At the same time unlike Me/ZnO/Me structure for the considered S0, A1, S1 modes of the Lamb elastic wave the essential increase in values of S sensitivity when using thick layers of Al is noted. In this case S at great values of the parameter h´f can be comparable or even exceed the values in the local maximum (fig. 4, a). At the same time it is also possible, for example, for A1 mode when the maximal value of sensitivity is noted at the moment of wave emergence, further the local minimum is observed and then S becomes almost dispersionless (fig. 4, b). In general the nature of parameter S changes in Me/AlN/Me structure is similar to changes of S in the three-layered plate of Me/ZnO/Me. Features of S sensitivity curves of Lamb and SH waves in Me/ZnO/Me and Me/AlN/Me structures. For both structures the characteristic feature is that only S0, A1 and S1 modes of the Lamb elastic wave possess the localized by maxima S values. Other modes - antisym- metric fundamental A0 with the lowest propagation veloc- ity and the mode with shear-horizontal polarization of SH0, SH1, SH2 - have no apparent features of sensitivity curves on the studied frequency range. Metal layers thicknesses increase results in the dis- placement of S sensitivity values maxima to the low- frequency area. At the same time S curve peaks amplitude and width with the increase in metal thickness for Me/ZnO/Me plate depend on the type of metal and the parameter h´f less than for Me/AlN/Me. When comparing various structures of Me/Piezoelectric/Me type it is revealed that sensitivity values the Lamb elastic wave modes at d/h = 0.001; 0.005; 0.01 differ insignificantly. The essential decrease of S sensitivity for the considered elastic wave modes at d/h > 0.05, in particular when using Mo, can limit metals layers of the considerable thickness appli- cation in sensor devices. Fig. 3. Dispersive dependence of the sensitivity for S0, A1 и S1 modes on the product h´f for the layered structure Me/ZnO/Me: a-c - aluninium; d-f - molybdenum. Ratio of the metal layers and piezoelectric d/h: 0.001, 0.005; 0.01; 0.05; 0.1 Рис. 3. Дисперсионные зависимости чувствительности мод S0, A1 и S1 от произведения h´f для слоистой структуры Me/ZnO/Me: а-в - алюминий; г-е - молибден. Отношение толщин металла и пьезоэлектрика d/h: 0,001, 0,005; 0,01; 0,05; 0,1 Тable 1 Максимальные значения чувствительности S в структуре Me/ZnO/Me The maximum sensitivity values of S in the Me/ZnO/Me structure d/h Mode S0 A1 S1 Metal Al Mo Al Mo Al Mo 0.001 S, m2/kg 251 246 248 247 249 248 h´f, m/s 2500 2400 4700 4700 7000 7000 0.005 S, m2/kg 248 244 246 242 248 243 h´f, m/s 2500 2400 4700 4700 6700 7000 0.01 S, m2kgг 246 240 244 237 245 237 h´f, m/s 2400 2400 4700 4600 6900 6900 0.05 S, m2/kg 225 199 231 190 233 177 h´f, m/s 2400 2300 4400 4400 6700 6500 0.1 S, m2/kg 200 155 212 141 227 124 h´f, m/s 2300 2100 4300 4100 6500 6300 Mean value of h´f in maximum S, m/s 2420 ± 104 2320 ± 162 4560 ± 114 4500 ± 149 6760 ± 114 6740 ± 188 Fig. 4. Dispersive dependence of the sensitivity for S0, A1 и S1 modes on the product h´f for the layered structure Me/AlN/Me: a-c - aluninium; d-f - molybdenum. Ratio of the metal layers and piezoelectric d/h: 0.001, 0.005; 0.01; 0.05; 0.1 Рис. 4. Дисперсионные зависимости чувствительности мод S0 , A1 и S1 от произведения h´f для слоистой структуры Me/AlN/Me: а-в - алюминий; г-е- молибден. Отношение толщин металла и пьезоэлектрика d/h: 0,001, 0,005; 0,01; 0,05; 0,1 Максимальные значения величины S в структуре Me/AlN/Me The maximum sensitivity values of S in the Me/AlN/Me structure Тable 2 d/h Mode S0 A1 S1 Metal Al Mo Al Mo Al Mo 0.001 S, m2/kg 502 493 329 327 323 321 h´f, m/s 5200 5200 9600 9600 14000 13800 0.005 S, m2kgг 500 489 327 318 321 312 h´f, m/s 5200 5000 9600 9400 13800 13600 0.01 S, m2kgг 497 479 324 306 319 302 h´f, m/s 5200 5000 9400 9200 13800 13400 0.05 S, m2kg 476 349 313 229 364 264 h´f, m/s 5000 4600 9000 7600 12000 11800 0.1 S, m2/kg 432 239 349 178 714 193 h´f, msс 4800 4200 7600 7400 12400 11200 Mean value of h´f in maximum S, m/s 5080 ± 222 4800 ± 497 9040 ± 1046 8640 ± 1309 13200 ± 1153 12760 ± 1465 The feature of the elastic wave S0 mode propaga- tion in Me/ZnO/Me layered structure is the emergence of S parameter negative values when using both Al and Mo metal layers, that is the value of the elastic wave phase velocity when metalizing, can increase (fig. 2, fig. 3, а, d). For Me/AlN/Me layered structure this effect is not ob- served. Conclusion. The influence of mass loading in the form of two metal layers (simulating the upper and lower electrode) on S sensitivity of elastic wave dispersion modes in piezoelectric layered structures Me/ZnO/Me and Me/AlN/Me change according to the frequency and d/h depends on the ratio of the metal layer acoustic imped- ance and acoustic line material. It is noted that the maxi- mal values of elastic wave modes sensitivity are reached at the layered Al/AlN/Al system configuration, that is in the system with low values of acoustic impedances of a bulk longitudinal wave for layers and piezoelectric acous- tic line materials. It is revealed that the maximal values of S sensitivity in all structures are observed in cases when a metal layer thickness is 2/3, 1, 3/2 of the corresponding Lamb elastic wave length. The results received during computer simulation can be useful when developing acoustic electronic devices on the Lamb waves and the SH-waves.

Sobre autores

O. Zolotova

Reshetnev Siberian State University of Science and Technology

Email: zolotova@sibsau.ru
31, Krasnoyarsky Rabochy Av., Krasnoyarsk, 660037, Russian Federation

S. Burkov

Siberian Federal University

79, Svobodny Av., Krasnoyarsk, 660041, Russian Federation

Bibliografia

Arquivos suplementares