2.1. SAW Propagation in Layered Media
Surface acoustic waves can travel along the surface of an elastic half-space with a phase velocity which depends on the crystallographic cut of the propagating medium and on the propagation direction. In single-material half-space, the confinement mechanism of the SAWs depends on the presence of a stress-free surface. The propagation of the SAWs along the surface of a piezoelectric medium is excited and detected by metal interdigitated transducers (IDTs) which consists of spatially periodic thin-film metal electrodes placed in contact with the medium surface; the fingers periodicity represents the wavelength λ of the travelling acoustic mode. By applying an alternating voltage to the launching IDT, a periodic electric field is imposed on the piezoelectric substrate’s surface and a SAW is produced because of the generation of a periodic strain field. The generated SAW travels through a certain distance (the wave path) and reaches the output IDT which detects the electric voltage associated with the SAW. The IDTs parameters (such as the IDT finger width and spacing, fingers overlapping, number of finger pairs) affect the IDTs performance; the double electrode split IDT structure minimizes internal reflections (it is called non-reflective transducers because the internal reflection phenomenon is shifted at the double of the centre frequency) that otherwise complicate the IDT response [
9,
10,
11,
12]. The ability of the split IDT to generate harmonics of the fundamental mode strongly depends on the metallization ratio of the fingers [
13,
14]: for 0.5 metallization ratio, the split finger IDT allows highly efficiency at the third harmonic and hence high working frequency ranges which would otherwise require expensive lithographic techniques; the 5
th and 7
th harmonics are suppressed as the resulting wave would have the same piezoelectric potential at all IDT fingers, meaning that no electrical energy can be transferred in the respective surface wave mode by applying a voltage between the IDT electrodes.
If the propagating medium consists of a layer/substrate structure, it can sustain the propagation of multiple SAWs: the first mode is called Rayleigh wave while the higher order modes are called Sezawa, R3, R4, and so on [
1]. The condition required for this to happen is that the velocity of the transverse BAW of the half-space material is larger than that of the layer material. The number of SAW modes and their velocities depend on the layer thickness, h: for very thin film thickness (h/λ << 1), only the fundamental Rayleigh mode propagates with a velocity close to the SAW velocity of the substrate material (about 3450 m/s for fused quartz); by increasing the layer thickness (h/λ >> 1), the Rayleigh velocity asymptotically reaches the SAW velocity of the layer material (2644 m/s for c-ZnO). The velocity of higher order modes asymptotically reaches the velocity of BAW shear in the layer material as the layer thickness increases. The amplitude profile of the Rayleigh wave is predominantly confined in the layer and decays exponentially with the depth, while that of the higher order modes has an exponential tail in the substrate. The latter modes have a layer thickness-to-wavelength cut-off at which the phase velocity is equal to the substrate shear BAW velocity. Right at the cut-off, the SAW mode couples with bulk modes and shows a leaky nature, as the acoustic power flows into the bulk substrate, thus resulting in a large insertion loss.
2.2. FEM Analysis
The propagation of SAWs along ZnO/SiO
2 structure was investigated by using 2D FEM analysis since the acoustic field was unchanged along the transverse direction, due to the crystallographic symmetry of the c-ZnO and the isotropy of the fused silica substrate. Comsol Multiphysics 6.1 software was used in the eigenfrequency and frequency domain studies. The models use a piezoelectric multiphysics coupling node with solid mechanics and electrostatic interfaces. The unit cell (width λ = 80 µm) consists in the following domains, as shown in
Figure 1a: Air (2·λ thick), ZnO (from 1 to 80 µm thick), SiO
2 (5·λ thick) and PML (2·λ thick). Four Al electrodes (0.15 µm thick and λ/8 wide) were placed onto the surface of the ZnO layer to perform the frequency domain study, as shown in
Figure 1b: two neighbors electrodes were connected to a fixed electric potential (terminal type: voltage, 1 V) and the other two were grounded. Periodic boundary conditions (
) are applied to the lateral sides of the unit cell. The bottom boundary is fixed (
), while electrically free or shorted boundary condition on the ZnO free surface (
) is used to calculate the electromechanical coupling coefficient (as explained later).
The material constants were extracted from the material library of COMSOL. The ZnO has 0.01 permittivity loss and 0.002 mechanical loss; SiO2 has isotropic mechanical loss equal to 0.01. An extremely fine mesh (i.e., maximum automatically generated physics-defined triangular elements) was chosen for all the FEM simulations to get more accurate results.
Eigenfrequency and frequency domain studies were performed to identify the type of the acoustic modes experimentally detected based on a comparison of the experimental data with numerical calculations.
Figure 2 a-c shows, as an example, the solid displacement, the longitudinal and shear displacement components (u and v) of the LS wave for a ZnO layer 4 µm thick.
The fundamental and harmonic Rayleigh waves have elliptical polarization and travel along the surface of the propagating medium, up to a depth of about 1.5·λ from the surface. The LS waves are mostly longitudinally polarized, and their vertical displacement component increases while approaching the cut off. The LS waves leak energy from the layer into the substrate and becomes evanescent, thus exhibiting their leaky nature as demonstrated by the substantial vertical displacement component penetrating into the SiO2 substrate.
The phase velocity of the SAWs were calculated for different ZnO layer thicknesses and fixed wavelength (λ = 80 µm) by performing the sweep parameter eigenfrequency study: the modes velocity was calculated as v = λf, being λ = 80 µm the wavelength and f the eigenfrequency solution.
Figure 3 shows the Rayleigh and Sezawa dispersion curves; the black points are the experimentally measured data which will be discussed in the next paragraphs.
For very thin film thicknesses, only the Rayleigh mode propagates; the Sezawa mode appears as it exceeds its cut-off value (
= 0.68). With further increases in h/λ (not shown in
Figure 3), several unattenuated higher order Rayleigh modes appear at a succession of values of
. Above the transverse BAW velocity of the substrate, no modes are allowed to propagate unattenuated: by assuming non-zero dielectric loss and mechanical damping, the velocity dispersion curves of the Sezawa and higher order mode extend between the transverse and longitudinal BAW velocity in the substrate where they become leaky modes, broadened through coupling to the bulk modes of the substrate, into which they radiate. At small h/λ, the LS velocity is close to that of the longitudinal BAW in SiO
2 while in the high h/λ limit it tends to the velocity of the shear vertical BAW in the ZnO.
The electromechanical coupling coefficient K
2, which is a measure of the electrical-to-acoustic energy conversion efficiency by electrical means, was evaluated from eigenfrequency analysis by calculating the phase velocity of the modes with electrically free (v
f) and shorted (v
s) boundary condition onto the ZnO free surface (
) and using the following approximate formula:
.
Figure 4a shows the K
2 dispersion curves for the Rayleigh and Sezawa waves.
The admittance Y vs frequency curves were calculated by performing the frequency domain study: the waves propagation loss was calculated by applying the formula
, where
,
is the group velocity,
is the phase velocity,
[
15];
and
were evaluated from the real part of the Y vs. frequency curve as the abscissa of the peak and as the range of frequencies covered by the half-power bandwidth.
Figure 4b shows the propagation loss dispersion curves for the Rayleigh and Sezawa waves for different values of the isotropic loss factor of the SiO
2 substrate; the shadowed portion of the α curve highlights a low-loss region for the Sezawa mode.
At very low h/λ ratios, the acoustic penetration depth of the Rayleigh wave extends beyond the interface between the thin film and the substrate, making the substrate crucial for losses. With increasing the layer thickness, the Rayleigh wave experiences a decreasing propagation loss and an increasing K2 which asymptotically reaches that of the ZnO half-space while the propagation loss tends to zero since the wave travels predominantly along the ZnO layer thus it becomes unaffected by the presence of the lossy substrate (the wave suffers only the loss from the layer).
The Sezawa wave is less surfacy than the Rayleigh wave as it travels inside both the substrate and the layer. It has a very low K
2 and high propagation loss for very small h/λ since the acoustic penetration depth extends beyond the film/substrate interface. From
Figure 4a it can be noticed that, below the cut off
, the K
2 of the LS increases and reaches a peak (about 3.7%) at h/λ = 0.375; as h/λ decreases further, the K
2 goes to zero because the ZnO layer thins and
α increases. For h/λ larger than the cut off, the K
2 tends to zero and the loss asymptotically reaches the loss in the substrate since the Sezawa wave is mostly confined in the SiO
2 substrate which is non-piezoelectric and lossy.
The α dispersion curves of the Rayleigh and Sezawa waves were calculated for different values of the SiO2 isotropic loss factor ranging from 0.01 up to 0.05. The Rayleigh wave seems to be unaffected by the choice of the substrate loss which is consistent with the surface character of the wave. On the contrary, the Sezawa wave is affected by the SiO2 loss (since the wave mostly travels into the SiO2 substrate) but the shape of the dispersion curve remains unchanged for the different values of the parameter.
From the calculations carried out, it turned out that the Sezawa mode is characterized by a dispersion of leaky waves: the cut-off frequency is not a discrete point in the phase velocity vs. h/λ curves but is a long-period termination of a wave surface that approaches the bulk longitudinal velocity of the substrate with increasing propagation loss except in a small low-loss region whose existence is confirmed by the measured data. Since the low-loss conditions for LS waves propagating in ZnO/SiO2 structure can approach velocity of 5840 m/s (equal to almost double the Rayleigh wave velocity), operation frequencies of several hundreds of MHz are achievable, which are suitable to implement high frequency SAW resonators and sensors.