Next Article in Journal
Influence of Weather Conditions and Mechanical Reclamation on Molding Sand with Alkali-Phenolic Binder for Manganese Cast Steel
Next Article in Special Issue
Features of the Defect Structure and Luminescence of Nominally Pure Lithium Niobate Crystals Produced Using Different Technologies
Previous Article in Journal
Experimental Investigation of the Spall Propagation Mechanism in Bearing Raceways
Previous Article in Special Issue
Silver Itaconate as Single-Source Precursor of Nanocomposites for the Analysis of Chloride Ions
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Materials
Volume 16
Issue 1
10.3390/ma16010069
materials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of SAW Temperature Properties in KTiOPO4 Single Crystal

by
Rinat Taziev
1 and
Victor Atuchin
1,2,3,4,*
1
Laboratory of Optical Materials and Structures, Institute of Semiconductor Physics, SB RAS, 630090 Novosibirsk, Russia
2
Research and Development Department, Kemerovo State University, 650000 Kemerovo, Russia
3
Department of Industrial Machinery Design, Novosibirsk State Technical University, 630073 Novosibirsk, Russia
4
R&D Center “Advanced Electronic Technologies”, Tomsk State University, 634034 Tomsk, Russia
*
Author to whom correspondence should be addressed.
Materials 2023, 16(1), 69; https://doi.org/10.3390/ma16010069
Submission received: 11 November 2022 / Revised: 19 December 2022 / Accepted: 19 December 2022 / Published: 21 December 2022
(This article belongs to the Special Issue Advances in Crystal Material: Design, Characterization, and Applications)
Browse Figures
Versions Notes

Abstract

:
The surface acoustic wave (SAW) properties of potassium titanyl phosphate (KTiOPO4, KTP) single crystal were evaluated by numerical methods. The phase velocity, electromechanical coupling coefficient, power flow deflection angle, and temperature coefficient of delay (TCD) were determined for different crystal cuts of KTP. It was shown that SAW has the electromechanical coupling coefficient of 0.59% and the TCD of 62 ppm/°C on the Z-cut and wave propagation direction along the crystal X + 70°-axis. For the Z-cut and wave propagation direction along the X-axis, the pseudo-surface wave (PSAW) is characterized by the coupling coefficient of 0.46% and the TCD value of 57 ppm/°C. The Bleustein–Gulyaev (BG) wave has the TCD value of 35 ppm/°C and 41 ppm/°C on the Y- and X-cuts of KTP, respectively.

1. Introduction

Potassium titanyl phosphate KTiOPO4 (KTP) is a ferroelectric crystal commonly used for the frequency doubling of Nd3+-doped laser sources. KTP is of great significance in the field of nonlinear optics due to its large nonlinear coefficients and a wide transparency range [1,2,3,4,5]. KTP crystals are used in the devices for the conversion of IR radiation to the THz wave and THz radiation to the THz wave [5]. The generation of THz wave is a parametric process in KTP, and this material has great advantages in comparison with LiNbO3:MgO crystals. In the KTP structure, wide-range substitutions are possible in cation and anion sublattices, and a great diversity of crystals is presently known in the KTP family [6,7,8,9,10,11,12]. Accordingly, a significant tuning of physical properties can be reached via a solid solution formation. It is known that KTP-type crystals are characterized by high thermal and chemical stabilities [6,9,13,14]. Presently, large-size monodomain KTP crystals of high optical quality and high electrical resistivity are available for industrial applications [15,16,17].
As a multifunctional material, KTP possesses piezoelectric properties, which can be useful for SAW device applications. The surface acoustic wave excitation by an interdigital transducer (IDT) in KTP was first experimentally investigated in [18]. A full set of elastic, piezoelectric, and dielectric constants of KTP was established in [19,20]. It was shown that the pseudo-surface wave and the Bleustein−Gulyaev (BG) wave are excited on Z- and Y-cuts, respectively, for the wave propagation direction along the X-axis in KTP [21]. The pseudo-surface acoustic wave is similar to the Rayleigh wave because it has no transverse displacement component. BG waves, due to their high sensitivity to fluid viscosity, are of considerable interest for the use in biomedical applications [22,23,24,25,26]. For the BG wave propagation, the displacement of material particles occurs only parallel to the substrate surface (there is no normal component), and it prevents its penetration into a liquid on the surface. Therefore, the BG wave cannot be used for microfluidic excitation, and it makes them suitable for measuring a liquid viscosity. The surface acoustic wave is widely used in the temperature, pressure sensors, and microfluidic applications [27,28,29,30,31,32,33]. Generally, the SAW has three components of mechanical displacement on the crystal surface: the first one is parallel to the direction of wave propagation, the second one is perpendicular to the direction of propagation, and the third one is normal to the substrate surface. Typically, the magnitude of the SAW amplitude reaches from a fraction to several nanometers, and it decreases exponentially with the depth into the crystal. In the layer with a thickness of about one wavelength, almost all the wave energy is concentrated, and in such a way that the SAW is strongly localized near the substrate surface. In addition, since the wave has both longitudinal and transverse displacement components, it can effectively “sense” the properties of an adjacent medium (gas, liquid, solid films, small particles, etc.). For microfluidic applications, the normal displacement component is important, since the Rayleigh wave penetrates a liquid drop on the substrate surface, thus, transforming into a leaky wave. The leaky wave creates the acoustic pressure in a liquid medium, which can cause the liquid drop movement called acoustic streaming, that forms the basis for many microfluidic devices [34,35,36,37,38,39,40,41,42,43]. To study the effects of the temperature change on the propagation of SAW, pseudo-SAW and BG waves, the dependence of material constants on temperature should be determined. This allows the temperature sensor creation on the base of SAW in many microfluidic devices.
Recently, the temperature dependence of the piezoelectric constants of KTP was determined with the use of single sample by the method of resonant ultrasound spectroscopy (RUS) [44]. It utilizes a wide resonance spectrum of bulk acoustic eigen modes excited in one cubic sample, containing all the necessary information on KTP material constants. However, at present, the characteristics of the surface, pseudo-surface and Bleustein–Gulyaev waves and their temperature properties are not studied in detail numerically or experimentally on KTP.
This work is aimed at the numerical study of the temperature properties of surface acoustic waves in KTP. For generality, it is also necessary to investigate such important SAW parameters as the phase velocity, electromechanical coupling coefficient, and power flow deflection angle on X-, Z-, and Z-rotated cuts of KTP. The paper is organized as follows: (1) The contour plots of the main SAW characteristics (phase velocity, electromechanical coupling coefficient, power flow deflection angle, and temperature coefficient of delay) are simulated by changing Euler angles. (2) For selected Z- and X-cuts of the KTP crystal, the SAW parameters are investigated in detail to find the SAW propagation directions which can be of potential interest for SAW device applications. (3) FEM simulations are performed in the time domain to investigate the selected SAW, pseudo-SAW, and Bleustein–Gulyaev wave propagation orientations in KTP.

2. Materials and Methods

KTP belongs to the orthorhombic crystal system (mm2 point group), and it has 17 independent material constants necessary to simulate its acoustic properties: nine elastic constants, five piezoelectric, and three dielectric constants. In our numerical simulations, all necessary material constants of the KTP crystal were taken from [19,20]. All these constants (elastic stiffness tensor CEij, piezoelectric stress tensor eij, and dielectric constants εSij) can be defined as the matrix elements [19,20]:
C E = ( 1 . 70 0 . 333 0 . 40 0 0 0 0 . 333 1 . 74 0 . 35 0 0 0 0 . 40 0 . 35 1 . 495 0 0 0 0 0 0 0 . 5919 0 0 0 0 0 0 0 . 5454 0 0 0 0 0 0 0 . 4313 ) , 10 11 N / m 2
e = ( 0 0 0 0 0.332 0 0 0 0 0.403 0 0 0.5428 0.707 1.636 0 0 0 ) , C / m 2
ε S / ε 0 = ( 11.44 0 0 0 11.47 0 0 0 15.47 ) ,   ε 0 = 8.854 * 10 − 12 F / m
The KTP crystal mass density was assumed to be 3.024 g/cm3 [19]. For aluminium electrodes, the isotropic elastic constants are taken to be as C11 = 1.11 × 1011 N/m2, C12 = 0.61 × 1011 N/m2, and C66 = 0.25 × 1011 N/m2, and the density is ρ = 2.695 g/cm3. The SAW parameter computations were performed using our FEM/BEM program (V.1.3, R. Taziev, Novosibirsk, Russia), described in [45]. Our program allows calculating SAW parameters, leaky and bulk wave properties, and simulate effective permittivity, matrix Green function, and the conductance of a finite number IDT electrode structure. In the simulation, the crystal cut is determined by two Euler angles α and μ, and the propagation direction vector is determined by the third Euler angle θ. For angle α = 0°, the crystal cut angle μ and wave propagation direction θ changed in the range from 0 to 90° with the step of 2° (determining Euler angles is provided in Figure 1e).
The SAW properties in piezoelectric materials are described by the system of wave equations [46,47]:
ρ ∂ 2 u j ∂ t 2 − C i j k l ∂ 2 u k ∂ x i ∂ x l − e k i j ∂ 2 φ ∂ x i ∂ x k = 0 , e i k l ∂ 2 u k ∂ x i ∂ x l − ε i k ∂ 2 φ ∂ x i ∂ x k = 0 ,   i , j , k , l = 1 , 2 , 3
here φ and uj are the electric potential and mechanical displacement vector, respectively, and ρ is the crystal density. The coordinate system for the SAW problem is taken with the x3 outward normal to stress-free substrate surface, and x1 is chosen as a wave propagation direction vector. To solve system (1), it is assumed that the solution is a linear combination of partial waves, each of which propagates with the same wave vector k and velocity V along axis x1, does not change in direction x2 and is attenuated in direction-x3 along the substrate depth. The partial waves are taken to be (uj, φ)n = (αj, α4)nexp(ikbnx3)exp(ik(x1 − Vt) with the slowness parameter bn to be chosen so that the partial wave satisfies the wave Equation (1) for a given value of V, where k is the propagation vector. By substituting the solution into Equation (1), eight roots bn for a given velocity V, which are complex-conjugate pairs, are determined. Only four roots bn satisfy the condition of vanishing of each partial wave at large substrate depths. Finally, the unknown velocity V is determined from the boundary conditions at the crystal surface: the vanishing of the normal components of elastic stresses at the free substrate surface and the continuity of the electric displacement or the equality to zero of the electric potential at the crystal surface. A detailed efficient iterative method of searching the velocity is proposed and described in [46,47].
An important SAW parameter is an electromechanical coupling coefficient for piezoelectric materials. It is defined as the percentage difference in the velocity between free surface (Vfree) and stress-free surface coated with an infinitely thin perfect conductor (Vshorted) [47,48]:
K2/2 = (Vfree − Vshorted)/Vfree
This value determines the magnitude of SAW wave generated by the IDT structure.
Another important SAW parameter is the power flow deflection angle (or beam steering angle). The power flow angle is the angle between the power flow vector and the propagation direction vector (phase velocity vector). The power deflection angle ϕ can be found to be as a derivative of phase velocity, with respect to the propagation direction angle θ [47,48]:
tan ϕ   =   1 V d V d θ
It is preferable if the angle value ϕ is small or equal to zero for SAW devices [49].
The temperature behavior of some SAW device characteristics is dependent on the temperature change in the delay time between transducers. For SAW resonators and narrow-band SAW filters, it is necessary to have the center device frequency which has low sensitivity to the substrate temperature change [49]. Oppositely, in temperature sensors, it is desirable to have a large device frequency sensitivity to the substrate temperature change with a good linearity [50]. The SAW device parameters sensitivity to the temperature variation are determined by the temperature coefficient of delay (TCD) [48]:
T C D     =   1 τ d τ d T   =   α − 1 V d V d T
where τ(T) = L(T)/V(T) is the delay time between input and output IDT, L(T) is the distance between two IDT, V(T) is the phase velocity, T is the substrate temperature, and α is the thermal expansion coefficient along the wave propagation direction in the crystal. The delay time τ(T) is calculated using formula τ(T) = L0(1 + α(T – 20 °C))/V(T), where L0 is the distance between input and output IDT at room temperature 20 °C, α is the thermal expansion coefficient along the wave propagation direction, and V(T) is the phase velocity at temperature T. The phase velocity V(T) is calculated using the material data of elastic, piezoelectric and dielectric constants tabulated in the temperature range from 20 to 140 °C with a step of 20 °C [44]. The temperature dependence of the crystal density is determined by relation ρ(T) = ρ0[1 − (α11 + α22 + α33)(T − 20 °C)], where ρ0 is the density at room temperature (20 °C). In the simulation (see Figure 1d), the average TCD value of wave was calculated by the formula:
T C D   =   1 τ d τ d T   ≈   τ ( 140   ° C ) − τ ( 20   ° C ) τ ( 20   ° C ) 1 120   ° C
The thermal expansion coefficients of KTP were taken from [19]: α11 = 6.8 × 10−6/°C, α22 = 9.6 × 10−6/°C, and α33 = −1.3 × 10−6/°C.

3. Results and Discussion

The SAW-simulated parameters versus two Euler angles μ and θ are presented in Figure 1. It can be seen that the phase velocity changes in the range from 3620 to 4040 m/s. The maximum value of SAW K2/2 value is about 0.6%. The SAW piezoelectric coupling coefficients, which are equal to the values between 0.25% and 0.56%, are in the region determined by the Euler angles μ and θ: μ = 0°~60° and θ = 10°~75° (see Figure 1b). The maximum SAW power flow deflection angle value is about 4° in the determined region of K2/2 (see Figure 1c). The SAW TCD, which values are in the wide range from 35 to 105 ppm/°C is shown in Figure 1d. It should be noted that the small of SAW TCD (35 ppm/°C) values are in the region of small electromechanical coupling coefficients. Large SAW TCD values in KTP are comparable to the values in the well-known LiNbO3 single crystal [48]. As can be seen, KTP does not have orientations with zero SAW TCD values.
The SAW velocity, ΔV/V, power flow deflection angle, and TCD versus propagation direction on the Z-cut of KTP are presented in Figure 2. The maximum ΔV/V value is equal to 0.59% along the X + 70°-axis wave propagation direction in the crystal (see Figure 2b). For this propagation direction, the SAW has the phase velocity of 4047 m/s, zero beam steering angle and the TCD of 61 ppm/°C. The dark rectangle in Figure 2a denotes the pseudo-surface wave solution [47]. This wave arises from the leaky wave branch, when the leaky wave propagates along the crystal X-axis. It has an in-plane displacement structure, without the transverse displacement component that causes the leaky wave attenuation. As soon as the wave propagation direction deviates from the crystal X-axis, there is an arising attenuation due to the appearance of the transverse bulk partial component. The pseudo-SAW has the following parameters: V = 3967 m/s, ΔV/V = 0.46% and TCD = 57 ppm/°C. The calculated pseudo-surface wave velocity is in a good agreement with the measured value of ≈3940 m/s [21]. It should be noted that the SAW degenerates into the nonpiezoelectric shear-horizontal (SH) bulk wave along the X-axis propagation direction in the crystal. The calculated SH bulk wave has the velocity value of 3776 m/s.
The SAW velocity ΔV/V, power flow deflection angle and TCD on the X-cut versus wave propagation direction θ on the crystal are displayed in Figure 3. For the X-cut and wave propagation direction along the crystal Z-axis, the SAW has the following parameters: V = 3927 m/s, ΔV/V = 0.13%, and TCD = 72 ppm/°C. The dark circle in Figure 3a denotes the Bleustein–Gulyaev wave solution. The TCD value equals 41 ppm/°C for the BG wave propagation direction along the Y-axis. In the same propagation direction, there exists the nonpiezoactive surface wave, which has a normal u3 and longitudinal u1 displacement structure similar to those in Rayleigh wave. The main SAW properties for different cuts and wave propagation directions are presented on Table 1.
The relative delay time as a function of substrate temperature with discrete steps of 20 °C for the SAW and pseudo-surface wave on the Z-cut of KTP are shown in Figure 4. As can be seen in Figure 4, the relative delay time slightly deviates from the monotonic linear function of the crystal temperature, thus the approximated Formula (6) can be used for the wave TCD estimation.
The effective permittivity versus slowness and the IDT conductance as a function of frequency for the Z-cut and wave propagation direction along the crystal X + 70°-axis are shown in Figure 5. The transducer had 100 Al electrodes with the period of 20 μm and the electrode thickness of 0.2 μm. The transducer aperture was 1600 μm. As can be seen, there is a leaky wave response near the SAW pole. The leaky wave has a phase velocity of about 4173 m/s, an attenuation of 0.14 dB per wavelength, and ΔV/V ≈ 0.072%. The electromechanical coupling coefficient of the leaky wave is eight times smaller than it is for the surface wave. For this reason, its parasitic contribution to the output signal is very small, as can be seen in Figure 5b for the SAW transducer conductance value.
The effective permittivity versus slowness and the IDT conductance as a function of frequency for Z-cut and the pseudo-surface wave propagation direction along the X-axis on the KTP crystal are illustrated in Figure 6. As can be seen in Figure 6b, there is only the pseudo-SAW excitation on the KTP crystal surface. The IDT had 100 Al electrodes with the period of 20 μm and electrode thickness 0.2 μm. The transducer aperture was 1600 μm.
The normalized absolute values of SAW and pseudo-SAW displacements generated by the IDT consisting of five electrode pairs with periodicity 20 μm in the uniform IDT are illustrated in Figure 7. The transducer is located at the center on the top of Figure 7a,b. The generation source was a harmonic voltage signal with the duration of two periods. For the Z-cut and the wave propagation direction along the crystal X + 70°-axis, the period was taken equal to 1/f0 = 1/200 MHz (see Figure 5b). For the pseudo-surface wave propagation direction along the X-axis, it was taken equal to 1/f0 = 1/197 MHz (Figure 6b). Here, f0 is the frequency of the maximum IDT conductance value in Figure 5 and Figure 6, respectively. It is shown in Figure 6a,b that the excitation of SAW and pseudo-SAW by the IDT flow in a similar way without the large difference in the bulk wave generation on the crystal surface. It can be seen that the IDT spreads a small amount of transducer energy into the bulk acoustic waves, thus propagating into the crystal depths. Both of these propagation orientations are suitable for SAW device applications. The time-domain FEM simulations were performed by the FlexPDE6 software (PDE Solutions Inc., 9408 E. Holman Rd.pokane Valley, WA, USA) [51].
The IDT conductance as a function of frequency for the BG wave propagation direction along the X-axis on the crystal Y-cut is shown in Figure 8. The transducer has 100 Ti electrodes with the period of 16 μm and electrode thickness 0.06 μm and acoustic aperture 1250 μm (Figure 8a). In the simulation, the isotropic elastic constants of Ti electrodes are assumed to be C11 = 1.57 × 1011 N/m2, C12= 0.75 × 1011 N/m2, and C66 = 0.41 × 1011 N/m2, and the density is ρ = 4.50 g/cm3. The frequency of the maximum of IDT conductance amplitude (≈268 MHz) is in good agreement with the measurements (≈266 MHz) [21]. In Figure 8b, the transducer has 100 Al electrodes with period 20 μm, and the electrode thickness varies from 0.4 to 2 μm. The transducer aperture was 1600 μm. The strong resonance behavior of IDT conductance as a function of electrode thickness was observed. It may be very useful for high-sensitivity liquid viscosity sensor applications. It should be noted that the BG wave has a very small electromechanical coupling coefficient value for infinitely thin electrodes [51]. As the electrode thickness increases, the IDT conductance and electromechanical coupling coefficient of the BG wave also increases [52,53].
The dependence of shear-wave displacement penetration into the KTP crystal depth, as a function of IDT electrode thickness in infinite periodic structure, is shown in Figure 9. It should be noted that the resonance frequencies obtained by FlexPDE are in the excellent agreement with the results obtained by our FEM/BEM program [45] (Figure 8b). The larger the electrode thickness is, the stronger the Bleustein–Gulayev wave localization to the crystal near-surface region is. It can be seen that most of the BG wave energy is accumulated inside the electrode as the electrode thickness increases. A combined (mixed) wave propagation in the IDT electrode structure can be considered as a superposition of the BG wave and the local shear oscillation mode within the electrode region. This phenomenon is indicated by an increase in the asymmetric conductivity behavior with the increasing electrode thickness (see Figure 8b). As noted above, most of the acoustic energy of this mixed wave is stored on the high-amplitude shear polarization mode confined within the electrode [54,55,56,57,58]. For a large electrode thickness, each IDT electrode oscillates as an elastic resonator coupled to its neighbours by evanescent waves [57,58,59].
In Figure 10 are the results of the time-domain FEM modeling of BG wave excitation by five electrode pairs at the center of uniform periodic structure of 41 electrodes with the 20 μm periodicity on the Y-cut and the wave propagation along the X-axis of KTP. The transducer is located at the center on top of Figure 10. The excitation voltage source is the harmonic pulse with a duration of two periods. For the electrode thickness of 0.08 μm, the period was taken equal to 1/f0 = 1/198 MHz, and for the electrode thickness of 2 μm, it was taken equal to 1/f0 = 1/190.8 MHz, where f0 is taken as the frequency of the IDT conductance maximum in Figure 8b. As it can be seen in Figure 10c, the BG wave displacement is highly localized at the crystal surface for the electrode thickness of h= 2 μm (h/λ = 10%). For this electrode thickness, there was a drastic change in the BG wave field structure. The wave is strongly localized in the near-surface region where the electrodes are located. A significant part of the wave energy is in the inner region of the electrodes, and that leads to a strong surface wave slowdown. A short BG wave pulse is strongly slowed down in the transducer structure and the reflection from the electrodes in the transducer induces oscillatory movements inside the uniform periodic structure, radiating its energy by small portions at the ends of the periodic structure. In this case, the BG wave is strongly coupled with the electrode oscillations and becomes as if trapped in the periodic structure of the electrodes and is periodically reflected from them.
It should be noted that such SAW parameters as phase velocity, electromechanical coupling coefficient, power flow deflection angle, and TCD can be measured using standard equipment from the 1970s [48,49,60]. However, a new technique is needed for the investigation of SAW short impulses [61,62,63,64]. A conventional method to study the SAW characteristics is the Laser Doppler Vibrometry (LDV) to measure out-of-plane displacement of wave [65]. However, the BG-type surface acoustic waves have no out-of-plane displacement components within the wave aperture; therefore, new methods should be developed for their detection [66].
In summary, detailed studies of the temperature properties of surface acoustic waves in the KTP crystal showed that there are no cuts and wave propagation directions in which the wave delay time changes weakly with the temperature change, similar to the ST-cut of alpha-quartz (Figure 1d). The SAW TCD values are in the range from 35 to 105 ppm/°C in the KTP crystal, which are of the same order of magnitude as for the well-known LiTaO3 and LiNbO3 single crystals. KTP possesses sufficiently high SAW coupling coefficients on the Z-cut for the use in SAW devices. For the Z-cut and the wave propagation direction along the X + 70°-axis, it is about 0.59%, and the SAW TCD is equal to 61 ppm/°C. For the Z-cut and the wave propagation direction along the X-axis, the pseudo-SAW has the electromechanical coupling coefficient of 0.46% and the TCD of 57 ppm/°C (Table 1). The pseudo-surface wave and SAW have zero power deflection angle value. The calculated pseudo-surface wave velocity value is in a good agreement with the experimental data [21]. In both KTP crystal cuts, spurious bulk acoustic waves are not generated by the IDT on the crystal surface, and these features make KTP attractive for SAW device applications.
As for the BG waves, their important properties are listed in Table 1. Their electromechanical coupling coefficient is vanishingly small, but it increases rapidly as the transducer electrode thickness increases (Figure 8b). For the Y-cut and the propagation direction along the X-axis, the BG wave, with the increasing electrode thickness in the periodic structure, is strongly localized in the KTP crystal near-surface region. It leads to a combined (mixed) wave propagation is considered as a superposition of BG wave and the local shear oscillation mode within the electrode region. That BG wave behavior can be useful for sensor applications.

4. Conclusions

The SAW phase velocity, electromechanical coupling coefficient, power flow deflection angle, and temperature coefficient of delay (TCD) were calculated for different crystal cuts and wave propagation directions in KTP. The isolines of the SAW temperature properties (SAW delay temperature coefficient), which were presented as a function of two Euler angles, show that the SAW has its TCD in the range from 35 to 105 ppm/°C in KTP. These TCD values are comparable to those of well-known LiTaO3 and LiNbO3 crystals. It was shown that KTP does not have orientations with zero SAW TCD values. Analysis of contour plots of SAW parameters (velocity, electromechanical coupling coefficient, power flow deflection angle, and temperature coefficient of delay (TCD)) showed that SAW has the maximum SAW piezoelectric coupling coefficient value (0.59%) on the Z-cut and the wave propagation direction along the X + 70°-axis. The SAW has the TCD of 61 ppm/°C. For the Z-cut and the wave propagation direction along the X-axis, the pseudo-surface wave has the high coupling coefficient of 0.46%, and its TCD is equal to 57 ppm/°C. The BG wave has TCD values of 35 and 41 ppm/°C on the Y- and X-cuts of KTP, respectively.
The calculated velocities of the pseudo-SAW and BG waves in KTP were in excellent agreement with the experimental data [21]. For the BG wave on the Y-cut and X-axis wave propagation direction, it was shown that, for an electrode thickness ratio to the period of 10%, the BG wave is strongly coupled with electrode oscillations, and it leads to a wave slowdown and strong localization. The short pulse of the BG wave is strongly slowed down inside the transducer structure, and the reflection from the electrodes induces the pulse to make oscillatory movements inside the homogeneous periodic structure, radiating its energy in small portions at the ends of the periodic structure.

Author Contributions

Conceptualization, V.A.; Methodology, R.T.; Software, R.T.; Formal analysis, R.T.; Investigation, R.T. and V.A.; Writing—original draft, R.T.; Writing—review & editing, V.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fan, T.Y.; Huang, C.E.; Hu, B.Q.; Eckardt, R.C.; Fan, Y.X.; Byer, R.L.; Feigelson, R.S. Second harmonic generation and accurate index of refraction measurements in flux-grown KTiOP04. Appl. Opt. 1987, 26, 2390–2394. [Google Scholar] [CrossRef] [PubMed]
  2. Vanherzeele, H.; Bierlein, J.D. Magnitude of the nonlinear-optical coefficients of KTiOPO4. Opt. Lett. 1992, 17, 982–984. [Google Scholar] [CrossRef] [PubMed]
  3. Hansson, G.; Karlsson, H.; Wang, S.; Laurell, F. Transmission Measurements in KTP and Isomorphic Compounds. Appl. Opt. 2000, 39, 5058–5069. [Google Scholar] [CrossRef] [PubMed]
  4. Wakui, K.; Takahashi, H.; Furusawa, A.; Sasaki, M. Photon subtracted squeezed states generated with periodically poled KTiOPO4. Opt. Express 2007, 15, 3568–3574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Mamrashev, A.; Nikolaev, N.; Antsygin, V.; Andreev, Y.; Lanskii, G.; Meshalkin, A. Optical Properties of KTP Crystals and Their Potential for Terahertz Generation. Crystals 2018, 8, 310. [Google Scholar] [CrossRef] [Green Version]
  6. Hagerman, M.E.; Poeppelmeier, K.R. Review of the Structure and Processing-Defect-Property Relationships of Potassi-um Titanyl Phosphate: A Strategy for Novel Thin-Film Photonic Devices. Chem. Mater. 1995, 7, 602–621. [Google Scholar] [CrossRef]
  7. Savatinova, I.; Savova, I.; Liarokapis, E.; Ziling, C.C.; Atuchin, V.V.; Armenise, M.N.; Passaro, V.M.N. A comparative analysis of Rb:KTP and Cs:KTP optical waveguides. J. Phys. D Appl. Phys. 1998, 31, 1667–1672. [Google Scholar] [CrossRef]
  8. Atuchin, V.V.; Alekseeva, O.A.; Kesler, V.G.; Pokrovsky, L.D.; Sorokina, N.I.; Voronkova, V.I. Chemical shifts of atomic core levels and structure of K1−xTi1−xSbxOPO4, x=0–0.23, solid solutions. J. Solid State Chem. 2006, 179, 2349–2355. [Google Scholar] [CrossRef]
  9. Atuchin, V.V.; Kesler, V.G.; Meng, G.; Lin, Z.S. The electronic structure of RbTiOPO4 and the effects of the A-site cation substitution in KTiOPO4-family crystals. J. Phys. Condens. Matter 2012, 24, 405503. [Google Scholar] [CrossRef]
  10. Dong, X.; Huang, L.; Hu, C.; Zeng, H.; Lin, Z.; Wang, X.; Ok, K.M.; Zou, G. CsSbF2SO4: An Excellent Ultraviolet Nonlinear Optical Sulfate with a KTiOPO4 (KTP)-Type Structure. Angew. Chem. Int. Ed. 2019, 58, 6528–6534. [Google Scholar] [CrossRef]
  11. Yang, F.; Wang, L.; Huang, L.; Zou, G. The study of structure evolvement of KTiOPO4 family and their nonlinear optical properties. Co-ord. Chem. Rev. 2020, 423, 213491. [Google Scholar] [CrossRef]
  12. Jiang, C.; Jiang, X.; Wu, C.; Huang, Z.; Lin, Z.; Humphrey, M.G.; Zhang, C. Isoreticular Design of KTiOPO4-like Deep-Ultraviolet Transparent Materials Exhibiting Strong Second-Harmonic Generation. J. Am. Chem. Soc. 2022, 144, 20394–20399. [Google Scholar] [CrossRef]
  13. Atuchin, V.V.; Kesler, V.G.; Maklakova, N.Y.; Pokrovsky, L.D.; Semenenko, V.N. Study of KTiOPO4 surface by X-ray photoelectron spectroscopy and reflection high-energy electron diffraction. Surf. Interface Anal. 2002, 34, 320–323. [Google Scholar] [CrossRef]
  14. Atuchin, V.; Maklakova, N.Y.; Pokrovsky, L.; Semenenko, V. Restoration of KTiOPO4 surface by annealing. Opt. Mater. 2003, 23, 363–367. [Google Scholar] [CrossRef]
  15. Rusov, V.A.; Gorchakov, A.V.; Belostotsky, A.L.; Sapojnikov, V.K.; Kaplun, A.B.; Meshalkin, A.B.; Shapovalov, V.N. KTiOPO4 electro-optical devices with extinction ratio up to 1:1000. Proc. SPIE 2001, 4350, 117–119. [Google Scholar] [CrossRef]
  16. Wang, X.; Yuan, X.; Li, W.; Qi, J.; Wang, S.; Shen, D. Flux growth of large potassium titanyl phosphate crystals and their electro-optical applications. J. Cryst. Growth 2002, 237–239, 672–676. [Google Scholar] [CrossRef]
  17. Roth, M.; Tseitlin, M. Growth of large size high optical quality KTP-type crystals. J. Cryst. Growth 2010, 312, 1059–1064. [Google Scholar] [CrossRef]
  18. Buritskii, K.S.; Dianov, E.M.; Kiselev, A.B.; Maslov, V.A.; Shcherbakov, E.A. Excitation of surface acoustic waves in Rb:KTP. Electron. Lett. 1991, 27, 1896–1897. [Google Scholar] [CrossRef]
  19. Chu, D.K.T.; Bierlein, J.D.; Hunsperger, R.G. Piezoelectric and acoustic properties of potassium titanyl phosphate (KTP) and its isomorphs. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1992, 39, 683–687. [Google Scholar] [CrossRef]
  20. Chu, D.K.T.; Hsiung, H.; Cheng, L.K.; Bierlein, J.D. Curie temperatures and dielectric properties of doped and undoped KTiOPO4 and isomorphs. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1993, 40, 819–824. [Google Scholar] [CrossRef]
  21. Chu, D.K.T.; Bierlein, J.D. Surface acoustic wave and Bleustein–Gulyaev wave generation in KTiOPO4crystals. Appl. Phys. Lett. 1993, 63, 2041–2043. [Google Scholar] [CrossRef]
  22. El Baroudi, A.; Le Pommellec, J.Y. Bleustein–Gulyaev waves in a finite piezoelectric material loaded with a viscoelastic fluid. Wave Motion 2021, 101, 102695. [Google Scholar] [CrossRef]
  23. Elhady, A.; Abdel-Rahman, E. Measurement of the electric permittivity using Bleustein–Gulyaev wave sensor. J. Micromech. Microeng. 2022, 32, 034004. [Google Scholar] [CrossRef]
  24. Elhady, A.; Basha, M.; Abdel-Rahman, E.M. Analysis of tunable Bleustein–Gulyaev permittivity sensors. J. Appl. Phys. 2021, 129, 164501. [Google Scholar] [CrossRef]
  25. Cui, W.; Mu, L.; Duan, X.; Pang, W.; Reed, M.A. Trapping of sub-100 nm nanoparticles using gigahertz acoustofluidic tweezers for biosensing applications. Nanoscale 2019, 11, 14625–14634. [Google Scholar] [CrossRef] [PubMed]
  26. Huang, Y.; Das, P.K.; Bhethanabotla, V.R. Surface acoustic waves in biosensing applications. Sens. Actuators Rep. 2021, 3, 100041. [Google Scholar] [CrossRef]
  27. Zu, H.; Wu, H.; Wang, Q.-M. High-Temperature Piezoelectric Crystals for Acoustic Wave Sensor Applications. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 2016, 63, 486–505. [Google Scholar] [CrossRef]
  28. Huang, J.; Cegla, F.; Wickenden, A.; Coomber, M. Simultaneous Measurements of Temperature and Viscosity for Viscous Fluids Using an Ultrasonic Waveguide. Sensors 2021, 21, 5543. [Google Scholar] [CrossRef]
  29. Giurgiutiu, V.; Xu, B.; Liu, W. Development and Testing of High-temperature Piezoelectric Wafer Active Sensors for Extreme Environments. Struct. Health Monit. 2010, 9, 513–525. [Google Scholar] [CrossRef]
  30. Taziev, R.M.; Kolosovsky, E.A.; Kozlov, A.S. Deformation-sensitive cuts for surface acoustic waves in α-quartz. In Proceedings of the 1993 IEEE International Frequency Control Symposium, Salt Lake City, UT, USA, 2–4 June 1993. [Google Scholar] [CrossRef]
  31. Taziev, R.M. Stress, temperature and pressure behavior of SAW on langasite plates. In Proceedings of the 2001 IEEE International Frequncy Control Symposium and PDA Exhibition (Cat. No.01CH37218), Seattle, WA, USA, 8 June 2001. [Google Scholar] [CrossRef]
  32. Go, D.B.; Atashbar, M.Z.; Ramshani, Z.; Chang, H.-C. Surface acoustic wave devices for chemical sensing and microfluidics: A review and perspective. Anal. Methods 2017, 9, 4112–4134. [Google Scholar] [CrossRef]
  33. Wu, G.; Kong, L.; Fan, M.; Yu, F.; Jiang, C.; Cheng, X.; Zhao, X. Electro-elastic features of YBa3(PO4)3 and YbBa3(PO4)3 crystals with pure face-shear mode for acoustic wave sensor applications. J. Materiomics 2021, 7, 828–836. [Google Scholar] [CrossRef]
  34. Caliendo, C.; Hamidullah, M. Guided acoustic wave sensors for liquid environments. J. Phys. D: Appl. Phys. 2019, 52, 153001. [Google Scholar] [CrossRef] [Green Version]
  35. Collins, D.J.; Ma, Z.; Ai, Y. Highly Localized Acoustic Streaming and Size-Selective Submicrometer Particle Concentration Using High Frequency Microscale Focused Acoustic Fields. Anal. Chem. 2016, 88, 5513–5522. [Google Scholar] [CrossRef] [PubMed]
  36. Jin, S.; Ye, G.; Cao, N.; Liu, X.; Dai, L.; Wang, P.; Wang, T.; Wei, X. Acoustics-Controlled Microdroplet and Microbubble Fusion and Its Application in the Synthesis of Hydrogel Microspheres. Langmuir 2022, 38, 12602–12609. [Google Scholar] [CrossRef] [PubMed]
  37. Athanassiadis, A.G.; Ma, Z.; Moreno-Gomez, N.; Melde, K.; Choi, E.; Goyal, R.; Fischer, P. Ultrasound-Responsive Systems as Components for Smart Materials. Chem. Rev. 2021, 122, 5165–5208. [Google Scholar] [CrossRef]
  38. Tayebi, M.; Yang, D.; Collins, D.J.; Ai, Y. Deterministic Sorting of Submicrometer Particles and Extracellular Vesicles Using a Combined Electric and Acoustic Field. Nano Lett. 2021, 21, 6835–6842. [Google Scholar] [CrossRef]
  39. Zhou, J.; Habibi, R.; Akbaridoust, F.; Neild, A.; Nosrati, R. Paper-Based Acoustofluidics for Separating Particles and Cells. Anal. Chem. 2020, 92, 8569–8578. [Google Scholar] [CrossRef]
  40. Wu, Z.; Cai, H.; Ao, Z.; Nunez, A.; Liu, H.; Bondesson, M.; Guo, S.; Guo, F. A Digital Acoustofluidic Pump Powered by Localized Fluid-Substrate Interactions. Anal. Chem. 2019, 91, 7097–7103. [Google Scholar] [CrossRef]
  41. Reyes, C.; Fu, L.; Suthanthiraraj, P.P.A.; Owens, C.E.; Shields, C.W., IV; López, G.P.; Charbonneau, P.; Wiley, B.J. The Limits of Primary Radiation Forces in Bulk Acoustic Standing Waves for Concentrating Nanoparticles. Part. Part. Syst. Charact. 2018, 35, 1700470. [Google Scholar] [CrossRef]
  42. Xu, X.; Wang, G.; Jiao, W.; Ji, W.; Jiang, M.; Zhang, X. Multi-level sorting of nanoparticles on multi-step optical waveguide splitter. Opt. Express 2018, 26, 29262–29271. [Google Scholar] [CrossRef]
  43. Ota, N.; Yalikun, Y.; Suzuki, T.; Lee, S.W.; Hosokawa, Y.; Goda, K.; Tanaka, Y. Enhancement in acoustic focusing of micro and nanoparticles by thinning a microfluidic device. R. Soc. Open Sci. 2019, 6, 181776. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Zhang, Y.; Tang, L.; Ji, N.; Liu, G.; Wang, J.; Jiang, H.; Cao, W. Temperature dependence of full set tensor properties of KTiOPO4 single crystal measured from one sample. J. Appl. Phys. 2016, 119, 124104. [Google Scholar] [CrossRef] [Green Version]
  45. Taziev, R.M. FEM/BEM for Simulation of LSAW Devices. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 2007, 54, 2060–2069. [Google Scholar] [CrossRef]
  46. Campbell, J.J.; Jones, W.R. A method for estimating optimal crystal cuts and propagation directions for excitation of piezoelectric surface waves. IEEE Trans. Sonics Ultrason. 1968, 15, 209–217. [Google Scholar] [CrossRef]
  47. Farnell, G.W. SAW Propagation in Piezoelectric Solids; Wave Electronics; Wilkinson, C.D.W., Collins, J.H., Dunsmuir, R., Masotti, L., De La Rue, R.M., Eds.; Elsevier Scientific Publishing Company: Amsterdam, The Netherlands, 1976; pp. 1–24. [Google Scholar]
  48. Slobodnik, A.J. Surface acoustic waves and SAW materials. Proc. IEEE 1976, 64, 581–595. [Google Scholar] [CrossRef]
  49. Morgan, D.P. Surface Acoustic Wave Filters: With Applications to Electronic Communications and Signal Processing; Academic Press: Amsterdam, The Netherlands, 2007; p. 448. ISBN 9780123725370. [Google Scholar]
  50. Xu, H.; Jin, H.; Dong, S.; Chen, J.; Song, X.; Xuan, W.; Shi, L.; Huang, S.; Zhang, P.; Luo, J. A langasite surface acoustic wave wide-range temperature sensor with excellent linearity and high sensitivity. AIP Adv. 2021, 11, 015143. [Google Scholar] [CrossRef]
  51. FlexPDE6. Available online: https://www.pdesolutions.com (accessed on 30 July 2016).
  52. Bright, V.M.; Hunt, W.D. Analysis of Bleustein–Gulyaev wave propagation under thin periodic metal electrodes. J. Appl. Phys. 1991, 70, 594–602. [Google Scholar] [CrossRef]
  53. Soluch, W. Calculation of Bleustein-Gulyaev wave parameters in KTiOPO4 crystal. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 1995, 42, 977. [Google Scholar] [CrossRef]
  54. Baghai-Wadji, A.-R.; Plessky, V.P.; Simonian, A.V. Propagation of Rayleigh type SAW in periodic structures with resonating elements. Sov. Phys. Acoust. 1992, 38, 442–449. [Google Scholar]
  55. Garova, E.A.; Maradudin, A.A.; Mayer, A.P. Interaction of Rayleigh waves with randomly distributed oscillators on the surface. Phys. Rev. B 1999, 59, 13291–13296. [Google Scholar] [CrossRef]
  56. Wegert, H.; Mayer, E.A.; Reindl, L.M.; Ruile, W.; Mayer, A.P. Interaction of SAWs with resonating structures on the surface. In Proceedings of the 2010 IEEE International Ultrasonics Symposium, San Diego, CA, USA, 11–14 October 2010; pp. 185–188. [Google Scholar] [CrossRef]
  57. Baghai-Wadji, A.R.; Maradudin, A.A. Shear horizontal surface acoustic waves on large amplitude gratings. Appl. Phys. Lett. 1991, 59, 1841–1843. [Google Scholar] [CrossRef]
  58. Maradudin, A.A.; Ryan, P.; McGurn, A.R. Shear horizontal acoustic surface shape resonances. Phys. Rev. B 1988, 38, 3068–3074. [Google Scholar] [CrossRef]
  59. Socié, L.; Benchabane, S.; Robert, L.; Khelif, A.; Laude, V. Surface acoustic wave guiding in a diffractionless high aspect ratio transducer. Appl. Phys. Lett. 2013, 102, 113508. [Google Scholar] [CrossRef]
  60. Oliner, A.A. Acoustic Surface Waves; Topics in Applied Physics; Springer: Berlin/Heidelberg, Germany, 1978; Volume 24, p. 334. [Google Scholar] [CrossRef]
  61. Frazer, T.D.; Knobloch, J.L.; Hoogeboom-Pot, K.M.; Nardi, D.; Chao, W.; Falcone, R.W.; Murnane, M.M.; Kapteyn, H.C.; Hernandez-Charpak, J.N. Engineering Nanoscale Thermal Transport: Size- and Spacing-Dependent Cooling of Nanostructures. Phys. Rev. Appl. 2019, 11, 024042. [Google Scholar] [CrossRef] [Green Version]
  62. Ronchi, A.; Sterzi, A.; Gandolfi, M.; Belarouci, A.; Giannetti, C.; Del Fatti, N.; Banfi, F.; Ferrini, G. Discrimination of nano-objects via cluster analysis techniques applied to time-resolved thermo-acoustic microscopy. Ultrasonics 2021, 114, 106403. [Google Scholar] [CrossRef]
  63. Nardi, D.; Travagliati, M.; Siemens, M.E.; Li, Q.; Murnane, M.M.; Kapteyn, H.C.; Ferrini, G.; Parmigiani, F.; Banfi, F. Probing Thermomechanics at the Nanoscale: Impulsively Excited Pseudosurface Acoustic Waves in Hypersonic Phononic Crystals. Nano Lett. 2011, 11, 4126–4133. [Google Scholar] [CrossRef]
  64. Gandolfi, M.; Peli, S.; Diego, M.; Danesi, S.; Giannetti, C.; Alessandri, I.; Zannier, V.; Demontis, V.; Rocci, M.; Beltram, F.; et al. Ultrafast Photoacoustic Nanometrology of InAs Nanowires Mechanical Properties. J. Phys. Chem. C 2022, 126, 6361–6372. [Google Scholar] [CrossRef]
  65. Smagin, N.; Djoumi, L.; Herth, E.; Vanotti, M.; Fall, D.; Blondeau-Patissier, V.; Duquennoy, M.; Ouaftouh, M. Fast time-domain laser Doppler vibrometry characterization of surface acoustic waves devices. Sens. Actuators A Phys. 2017, 264, 96–106. [Google Scholar] [CrossRef]
  66. Elhady, A.; Abdel-Rahman, E.M. Characterization of Shear Horizontal Waves Using a 1D Laser Doppler Vibrometer. Sensors 2021, 21, 2467. [Google Scholar] [CrossRef]
Materials 16 00069 g001 550
Figure 1. Isolines of SAW parameters as a function of two Euler angles μ and θ in KTP: (a) phase velocity (m/s); (b) K2/2 (%); (c) power deflection angle ϕ (degrees); (d) TCD (ppm/°C); and (e) determining the Euler angle (α, μ, θ) rotation of the Cartesian coordinate system (x1,x2,x3) from the initial crystallophysic orthogonal coordinate axes (X, Y, Z) of the crystal. For Euler angles (0,0,0), the Cartesian coordinate system (x1,x2,x3) coincides with the crystallophysic orthogonal coordinate axes (X, Y, Z) of the crystal.
Figure 1. Isolines of SAW parameters as a function of two Euler angles μ and θ in KTP: (a) phase velocity (m/s); (b) K2/2 (%); (c) power deflection angle ϕ (degrees); (d) TCD (ppm/°C); and (e) determining the Euler angle (α, μ, θ) rotation of the Cartesian coordinate system (x1,x2,x3) from the initial crystallophysic orthogonal coordinate axes (X, Y, Z) of the crystal. For Euler angles (0,0,0), the Cartesian coordinate system (x1,x2,x3) coincides with the crystallophysic orthogonal coordinate axes (X, Y, Z) of the crystal.
Materials 16 00069 g001
Materials 16 00069 g002 550
Figure 2. The propagation direction dependence of SAW parameters from the X-axis on the Z-cut of KTP: (a) phase velocity; (b) ΔV/V; (c) power flow deflection angle (degrees); and (d) TCD. The dark rectangle denotes the pseudo-surface wave solution.
Figure 2. The propagation direction dependence of SAW parameters from the X-axis on the Z-cut of KTP: (a) phase velocity; (b) ΔV/V; (c) power flow deflection angle (degrees); and (d) TCD. The dark rectangle denotes the pseudo-surface wave solution.
Materials 16 00069 g002
Materials 16 00069 g003 550
Figure 3. SAW propagation direction dependences from the Y-axis on the X-cut of KTP: (a) phase velocity; (b) ΔV/V; (c) power flow deflection angle (degrees); and (d) TCD. The dark circle denotes the BG wave solution.
Figure 3. SAW propagation direction dependences from the Y-axis on the X-cut of KTP: (a) phase velocity; (b) ΔV/V; (c) power flow deflection angle (degrees); and (d) TCD. The dark circle denotes the BG wave solution.
Materials 16 00069 g003
Materials 16 00069 g004 550
Figure 4. Relative delay time as a function of substrate temperature: (a) for the SAW propagation along the X + 70°-axis and (b) for the pseudo-surface wave propagation along the X-axis on the Z-cut of KTP.
Figure 4. Relative delay time as a function of substrate temperature: (a) for the SAW propagation along the X + 70°-axis and (b) for the pseudo-surface wave propagation along the X-axis on the Z-cut of KTP.
Materials 16 00069 g004
Materials 16 00069 g005 550
Figure 5. (a) Effective permittivity εeff/ε0 versus dimensionless slowness VSAW/V and (b) surface wave IDT conductance response versus frequency for the uniform IDT having 100 Al electrodes with period λ = 20 μm and electrode thickness h/λ = 1% on the Z-cut and the wave propagation direction along the crystal X + 70°-axis.
Figure 5. (a) Effective permittivity εeff/ε0 versus dimensionless slowness VSAW/V and (b) surface wave IDT conductance response versus frequency for the uniform IDT having 100 Al electrodes with period λ = 20 μm and electrode thickness h/λ = 1% on the Z-cut and the wave propagation direction along the crystal X + 70°-axis.
Materials 16 00069 g005
Materials 16 00069 g006 550
Figure 6. (a) Dependence of effective permittivity εeff/ε0 versus dimensionless slowness VShear/V and (b) IDT conductance as a function of frequency for a homogeneous IDT having 100 Al electrodes with period λ = 20 μm and electrode thickness h/λ = 1% on the Z-cut and the pseudo-surface wave propagation direction along the crystal X-axis.
Figure 6. (a) Dependence of effective permittivity εeff/ε0 versus dimensionless slowness VShear/V and (b) IDT conductance as a function of frequency for a homogeneous IDT having 100 Al electrodes with period λ = 20 μm and electrode thickness h/λ = 1% on the Z-cut and the pseudo-surface wave propagation direction along the crystal X-axis.
Materials 16 00069 g006
Materials 16 00069 g007 550
Figure 7. Normalized absolute value of acoustic wave displacement excited by the IDT consisting of five electrode pairs with period λ = 20 μm and thickness h/λ = 1%: (a) for the Z-cut and the pseudo-surface wave propagation direction along the X-axis, and (b) for the Z-cut and the SAW propagation direction along the X + 70°-axis on KTP.
Figure 7. Normalized absolute value of acoustic wave displacement excited by the IDT consisting of five electrode pairs with period λ = 20 μm and thickness h/λ = 1%: (a) for the Z-cut and the pseudo-surface wave propagation direction along the X-axis, and (b) for the Z-cut and the SAW propagation direction along the X + 70°-axis on KTP.
Materials 16 00069 g007
Materials 16 00069 g008 550
Figure 8. (a) Bleustein–Gulyaev wave IDT conductance response versus frequency for uniform IDT with 100 Ti electrodes with period λ = 16 μm and electrode thickness 60 nm; (b) for the uniform IDT with 100 Al electrodes with period λ = 20 μm and different thicknesses on the Y-cut of KTP and the wave propagation direction along the crystal X-axis.
Figure 8. (a) Bleustein–Gulyaev wave IDT conductance response versus frequency for uniform IDT with 100 Ti electrodes with period λ = 16 μm and electrode thickness 60 nm; (b) for the uniform IDT with 100 Al electrodes with period λ = 20 μm and different thicknesses on the Y-cut of KTP and the wave propagation direction along the crystal X-axis.
Materials 16 00069 g008
Materials 16 00069 g009 550
Figure 9. Spatial distribution of normalized amplitude U2 into the substrate depth for different IDT electrode thicknesses on the Y-cut and the X-axis wave propagation direction in KTP (calculated by FlexPDE [51]). One period of IDT containing the infinite periodic electrode structure with period λ = 20 μm is shown: (a) the electrode thickness h/λ = 4% and resonance frequency f0 = 210.6 MHz (resonance frequency f0 = 210.0 MHz calculated by our FEM/BEM program [45]); (b) the electrode thickness h/λ = 6% and resonance frequency f0 = 205.6 MHz (resonance frequency f0 = 205.7 MHz calculated by our FEM/BEM program); (c) the electrode thickness h/λ = 10% and resonance frequency f0 = 190.8 MHz (resonance frequency f0 = 190.8 MHz calculated by our FEM/BEM program).
Figure 9. Spatial distribution of normalized amplitude U2 into the substrate depth for different IDT electrode thicknesses on the Y-cut and the X-axis wave propagation direction in KTP (calculated by FlexPDE [51]). One period of IDT containing the infinite periodic electrode structure with period λ = 20 μm is shown: (a) the electrode thickness h/λ = 4% and resonance frequency f0 = 210.6 MHz (resonance frequency f0 = 210.0 MHz calculated by our FEM/BEM program [45]); (b) the electrode thickness h/λ = 6% and resonance frequency f0 = 205.6 MHz (resonance frequency f0 = 205.7 MHz calculated by our FEM/BEM program); (c) the electrode thickness h/λ = 10% and resonance frequency f0 = 190.8 MHz (resonance frequency f0 = 190.8 MHz calculated by our FEM/BEM program).
Materials 16 00069 g009
Materials 16 00069 g010 550
Figure 10. The time dependence of the BG wave excitation by a short harmonic voltage in a homogeneous IDT structure on the Y-cut and the wave propagation direction along the crystal X-axis. The IDT consists of 5 electrode pairs in a periodic structure of 41 electrodes with period λ = 20 μm and various electrode thicknesses h/λ: (a) the spatial distribution of electric potential on the electrodes at the excitation time, (b,c) the amplitude of shear wave displacement distribution for electrode thicknesses, h = 0.08 μm and h = 0.2 μm, respectively. Here, τ is the duration time of wave excitation by the transducer. The FlexPDE software is used for FEM simulations [51].
Figure 10. The time dependence of the BG wave excitation by a short harmonic voltage in a homogeneous IDT structure on the Y-cut and the wave propagation direction along the crystal X-axis. The IDT consists of 5 electrode pairs in a periodic structure of 41 electrodes with period λ = 20 μm and various electrode thicknesses h/λ: (a) the spatial distribution of electric potential on the electrodes at the excitation time, (b,c) the amplitude of shear wave displacement distribution for electrode thicknesses, h = 0.08 μm and h = 0.2 μm, respectively. Here, τ is the duration time of wave excitation by the transducer. The FlexPDE software is used for FEM simulations [51].
Materials 16 00069 g010
Table
Table 1. SAW parameters on the different single KTP crystal cuts.
Table 1. SAW parameters on the different single KTP crystal cuts.
Crystal Cut and Wave Propagation DirectionWave Type Velocity, m/sΔV/V, %TCD, ppm/°C
Y,XBGW42890.01835
X,YBGW44840.03741
Z,XPSAW39670.4657
Z,YPSAW41130.6662
Z,X + 70°SAW40470.5961
X,ZSAW39270.1372
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Taziev, R.; Atuchin, V. Analysis of SAW Temperature Properties in KTiOPO4 Single Crystal. Materials 2023, 16, 69. https://doi.org/10.3390/ma16010069

AMA Style

Taziev R, Atuchin V. Analysis of SAW Temperature Properties in KTiOPO4 Single Crystal. Materials. 2023; 16(1):69. https://doi.org/10.3390/ma16010069

Chicago/Turabian Style

Taziev, Rinat, and Victor Atuchin. 2023. "Analysis of SAW Temperature Properties in KTiOPO4 Single Crystal" Materials 16, no. 1: 69. https://doi.org/10.3390/ma16010069

APA Style

Taziev, R., & Atuchin, V. (2023). Analysis of SAW Temperature Properties in KTiOPO4 Single Crystal. Materials, 16(1), 69. https://doi.org/10.3390/ma16010069

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