1. Introduction
Surface acoustic waves (SAWs) are elastic waves that propagate on a solid surface, including Rayleigh waves propagating on the surface of a substrate and decaying along the depth direction of a semi-infinite medium [
1,
2], Love waves propagating in a layer over a substrate [
3,
4,
5], and so on. As SAWs have been widely used in nondestructive testing, biological experiments, sensors, and seismology, studying the propagation of SAWs on the surface of a solid medium or in a layer is of great importance. During the propagation of SAWs in piezoelectric materials, the piezoelectric effect leads to variation in coupled mechanical displacement and electric potential. Resulting from the semiconductor effect, the variation in electric potential can lead to the migration and diffusion of hole and electron carriers in some practical piezoelectric materials, i.e., piezoelectric semiconductors (PSCs). Due to piezoelectric and semiconductor effects, PSCs have been used to make devices for energy harvesting and acoustic wave amplification. Based on the classical piezoelectric basic equations with respect to stress and electric displacement, some additional equations concerning charge carrier current density are theoretically added [
6,
7,
8,
9], which have been utilized to investigate the propagation of generalized Rayleigh SAWs in a PSC half-space [
1] and a PSC thin film over an elastic half-space [
10], bulk waves in a PSC slab sandwiched by two piezoelectric half-spaces [
11] and a PSC half-space [
12], and shear-horizontal (SH) waves in single-layered and multilayered PSC plates [
13].
As a result of the different doping methods used, some practical PSCs can be either p-type or n-type PSCs, e.g., zinc oxide (ZnO), while some practical PSCs can only be n-type PSCs, e.g., gallium nitride (GaN). When an n-type PSC and metal are bonded together, a Schottky junction appears inevitably at the interface between these two bonded solids [
14,
15]. Compared with the physical properties of the bulk material, the Schottky junction often has special physical properties. For instance, the steady electron carrier concentrations inside the Schottky junction can vary dramatically, which may lead to drastic changes in the electric potential and electron carrier concentration due to the piezoelectric effect and the semiconductor effect. The influence of the Schottky junction on the propagation characteristics of SAWs should also be considered with the decrease in the geometric dimensions of the metal and PSC employed. Compared with the above literature, the discussion of the influence of Schottky junctions on the propagation characteristics of a shear horizontal wave is relatively lacking; this motivated our previous research work [
16], in which the dispersion and attenuation characteristics of a shear horizontal wave in an n-type PSC substrate with a metal layer were studied. The Schottky junction that appeared at the interface between the PSC substrate and the metal layer was considered an electrical gradient layer.
On the other hand, in order to describe the propagation characteristics of SAWs in micro or nanometer-scale materials, the surface and interface effect theory has been established, in which the surface or interface can be regarded as an imperfect interface that leads to a jump in the wave motion quantities between the two sides of the surface or interface [
17,
18]. With the consideration of the piezoelectric effect, surface or interface imperfect boundary conditions, governing equations, and interface characteristic lengths were introduced, and then some surface and interface piezoelectric theories were also established and utilized to investigate the propagation characteristics of elastic waves in micro- or nanometer-scale piezoelectric materials [
19]. Based on Hamilton’s principle and a variational approach, Wang et al. established a generalized dynamic model to depict the wave propagation properties in SAW nano-devices [
19]. Utilizing surface elasticity, surface piezoelectricity, and surface permittivity, the surface effect was considered to investigate the Love wave propagation in a typical surface acoustic wave device composed of a piezoelectric ceramic transducer film and an aluminum substrate and the shear horizontal vibration of a piezoelectric plate. Based on the surface/interface elasticity (effect), Zhu et al. derived formulations and carried out the calculation of the wave motion characteristics of SH and Lamb waves in layered nanostructures [
20]. It was found that the surface/interface effect can be significant when the wavenumber is large and the thickness of the layer is small, as these characteristics affect the dispersion curves, mode shapes, and kinetic and strain energy on the wave. Making use of the dual variable and position method, Tian et al. investigated the propagation characteristics of SH waves in multilayered PSC plates with imperfect interfaces and discussed the effect of stacking sequences and imperfect interfaces [
13]. It was observed that the mechanical imperfect interface can increase the real part of the wavenumber and decrease the imaginary part of the wavenumber for the 1st and 2nd modes, while the electric or electric current imperfect interface has nearly no influence on the dispersion and attenuation curves. Under the mathematical framework of surface/interface elasticity theory, Enzevaee and Shodja formulated the propagation of the torsional surface waves in a medium consisting of a functionally graded substrate bonded to a thin piezoelectric over-layer [
21]. It is observed that the surface/interface effects result in a reduction in the phase velocity and the phase velocity reduces by increasing surface/interface parameters for all modes at a specific wave number. Based on the transfer matrix method and the Bloch theorem, Guo et al. studied the influence of the interface effect on the dispersion relation of anti-plane elastic waves in nano-scale one dimensional piezoelectric and PSC phononic crystal [
22]. It is found that homo- and hetero-junctions can be approximately treated as rigid imperfect interfaces, which can enhance the width of band-gaps of anti-plane elastic waves, especially for high-frequency short-wavelength anti-plane elastic waves. Based on surface effect of nanostructures and first-order shear deformation theory, Zhang et al. developed a PSC nanoplate model to investigate the influence of surface effect on the natural characteristics and electromechanical responses of the nanoplate subjected to dynamic forces [
23]. It is found that surface effect causes the stiffness to be larger, the natural frequencies to be larger, and the piezoelectric effect to be stronger of the nanoplate. Based on the Gurtin-Murdoch surface model, Li et al. conduct an analysis of horizontally polarized shear (SH) waves propagating in an infinite n-type PSC plate with nano-thickness and elaborate the influence of the surface effect, initial electron concentration, and thickness of the PSC plate on the propagation characteristics of SH waves, such as frequency spectrum, phase velocity, and cut-off frequency [
24]. It is found that the propagation characteristics of SH waves are size-dependent, and the macroscopic effective shear rigidity is softened at the nanoscale due to surface effect, which thus leads to the reduction of the frequency and phase velocity of SH waves.
Since the geometrical thickness of the Schottky junction is common at the micro-scale which is much smaller than that of the body materials, it can also be treated as an electrically imperfect interface without geometrical thickness but with piezoelectric and semiconductor effects. Based on the physical properties of the Schottky junction and the existing interface piezoelectric theories, the theoretical model of the interface effect of the Schottky junction can also be established to better investigate the influence of the Schottky junction on the propagation characteristics of SAW in micro or nanometer-scale PSC materials, which motivates this paper’s work. First, considering the coupling mechanical displacement, electric potential, and charge carrier perturbation, it is mathematically derived that the basic equations of a PSC semi-infinite medium in
Section 2, which consists of a lower n-type PSC substrate, an upper grounded metal layer, and a Schottky junction that appears at the interface between these two regions. Two biasing electric fields are further considered in the basic equations of the lower n-type PSC substrate, which provides evident influence on the physical properties of the Schottky junction. Then, through the spectral method [
25,
26,
27], two equivalent mathematical models are established to investigate the dispersion and attenuation characteristics of the shear horizontal wave that propagates in this PSC semi-infinite medium, instead of the state transfer equation method. Searching for complex roots of complex transcendental equations can be avoided, and the numerical calculation of the dispersion and attenuation curves of the shear horizontal wave is transformed into solving the eigenvalues. To improve the efficiency of the above numerical calculation, the calculating intervals of definite integrals are adjusted reasonably. In the first mathematical model established in
Section 3, the Schottky junction is considered as a two-dimensional electrically imperfect interface without geometrical thickness but with piezoelectric and semiconductor characteristics. To legitimately confirm the interface characteristic length that describes the electrically imperfect interface boundary conditions of the Schottky junction, the second mathematical model is established in
Section 4, in which the Schottky junction is considered as an electrical gradient layer. Therefore, the second mathematical model is also a validation of our previous mathematical model. Finally, based on the above two equivalent mathematical models, the numerical results of the dispersion and attenuation curves of shear horizontal waves are provided in
Section 5. The numerical value and variation law of the interface characteristic length are also legitimately confirmed and studied. Some numerical results are compared our previous numerical results to verify the correctness and superiority of the research work in this paper. After that, the conclusion remarks are given in
Section 6.
2. Problem Formulation and Basic Equations
Consider a PSC semi-infinite medium, formed by a lower n-type PSC substrate and an upper grounded metal layer, as shown in
Figure 1. The upper grounded metal layer is isotropic. The lower n-type PSC substrate is transversely isotropic in the
plane. The
-direction is the polarization direction. The thickness of the upper grounded metal layer
is at the micro-scale. Both the upper metal layer and lower PSC substrate are homogeneous. Therefore, a Schottky junction appears at the interface between these two regions, in which the steady electron carrier concentration is inhomogeneous along the
axis. The geometrical thickness of the Schottky junction is
, which is much smaller than that of the upper metal layer. The horizontal and vertical biasing electric fields applied to the PSC semi-infinite medium are
and
, respectively. Consider the propagation of the shear horizontal wave in this PSC substrate in the
plane along the
axis. The Schottky junction and biasing electric fields provide evident influence on the dispersion and attenuation characteristics of the shear horizontal wave. To establish the mathematical models of the shear horizontal wave that propagates in this piezoelectric semiconductor semi-infinite medium and investigate the influence of the Schottky junction on the dispersion and attenuation characteristics of the shear horizontal wave, the basic equations of the upper metal layer and lower PSC substrate are mathematically derived in this section. To distinguish these two regions in the latter formulation, the physical parameters of these two regions are indicated by the superscript ‘
’ and ‘
’, respectively.
Under the supposition of small deformations and the quasi-static electric field approximation and with the consideration of biasing electric fields, the constitutive equations of the lower n-type PSC substrate [
11,
12] are
where
,
,
,
,
, and
are the displacement vector, electric potential scalar, perturbations of electron carrier concentration, Cauchy stress tensor, electric displacement vector and electron carrier current density vector, respectively. With the consideration of the piezoelectric and semiconductor effects,
,
,
,
, and
are the elastic, piezoelectric, dielectric, electron carrier migration, and electron carrier diffusion parameter tensors, respectively. The carrier charge is
. The steady electron carrier concentration is
, which is a function of the space coordinate
and the vertical biasing electric field
mathematically
where
are the donor concentration i.e., the doping concentration, in the lower n-type PSC substrate. The geometrical thickness of the Schottky junction is
where
and
are the work functions for the metal and n-type PSC substrate, respectively. The detailed derivation of the above two physical quantities is shown in the
Appendix A. Based on plane strain assumption, the displacements, electric potential, and carrier concentration perturbation of the shear horizontal wave are all functions of space coordinates
and
and time coordinate
mathematically
Inserting Equation (4) into Equation (1) leads to the Cauchy stress components
electric displacement components
and electron carrier current density components
With the consideration of the piezoelectric and semiconductor effects, the governing equations of the lower n-type PSC substrate [
11,
12] are
where
is the mass density. Inserting Equations (5)–(7) into Equation (8) leads to the governing equations expressed by the displacements, electric potential and carrier concentration perturbations
Without the piezoelectric and semiconductor effects, under the supposition of small deformations, the constitutive equations of the upper metal layer [
28] are
where
,
, and
are the displacement vector, Cauchy stress tensor, and elastic parameter tensor of the upper metal layer, respectively. Based on the plane strain assumption, the displacement component
is also a function of space coordinates
and
and time coordinate
mathematically
Inserting Equation (11) into Equation (10) leads to the Cauchy stress components of the upper metal layer
The governing equation of the upper metal layer [
28] is
where
is the mass density of the upper metal layer. Inserting Equation (12) into Equation (13) leads to
On the surface
, consider the surface condition of the upper metal layer with respect to the Cauchy stress [
28]
Equations (9), (14) and (15) provide the basic equations of the PSC semi-infinite medium, which can be utilized to establish the mathematical model of the shear horizontal wave that propagates in this PSC semi-infinite medium through the spectral method. However, as the core part of the PSC semi-infinite medium, the basic equations of the Schottky junction are not fully given. Therefore, when the Schottky junction is treated as an electrically imperfect interface or an electrical gradient layer, two different groups of basic equations can be given and two equivalent mathematical models can be established in the next two sections. In order to improve the computational efficiency of these two mathematical models, we reasonably modify some calculating intervals of definite integrals in the spectral method.
3. The First Mathematical Model
According to Equation (4), assume the mathematical expressions of the displacement, electric potential, and carrier concentration perturbation of the shear horizontal wave in the lower n-type PSC substrate as
where
and
are the apparent wave number component and angular frequency of the shear horizontal wave in the lower n-type PSC substrate, respectively. Based on the spectral method [
25,
26,
27], the Laguerre polynomials
are used to expand these three undetermined amplitudes
,
, and
, which are all functions of the space coordinate
mathematically
where
,
and
are expansion coefficients and
Inserting Equations (16) and (17) into Equation (9), multiplying both sides by
and taking the definite integrals lead to
According to Equation (11), assume the mathematical expressions of the displacement of the shear horizontal wave in the upper metal layer as
where
and
are the apparent wave number component and angular frequency of the shear horizontal wave in the upper metal layer, respectively. Based on the spectral method [
25,
26,
27], the Legendre polynomials
are used to expand the undetermined amplitude
, which is also a function of the space coordinate
mathematically
where
is expansion coefficients,
.
is an integer.
Inserting Equations (22) and (23) into Equation (14), multiplying both sides by
and taking the definite integrals leads to
Considering the geometrical thickness of the Schottky junction is much smaller than that of the upper metal layer, i.e.,
, the Schottky junction is treated as a two-dimensional electrically imperfect interface without geometrical thickness but with piezoelectric and semiconductor effects in this mathematical model. The interface boundary condition to the Cauchy stress and electric displacement at the Schottky junction is
where
is the electric displacement vector of the electrically imperfect interface,
and
and
are the interface piezoelectric and dielectric parameter tensors, respectively. The interface boundary condition to the displacement and electric potential at the boundary
is
where
and
are the displacement vector and electric potential of the electrically imperfect interface, respectively. Inserting the interface normal vector
into Equation (26) leads to
Inserting Equation (28) into Equation (27) leads to
and
where
is the interface characteristic length related to the interface piezoelectric parameter. Inserting Equations (30) and (31) into Equations (28) and (29) leads to
Inserting Equations (16), (17), (22), and (23) into Equations (15), (32) and (33) leads to
and
Inserting Equations (34) and (35) into Equation (25) leads to
Inserting Equations (36)–(38) into Equations (19)–(21) leads to
Combining Equations (40)–(43) and letting
changes from
to
leads to
where the unknown expansion coefficients vector is
The explicit expressions of coefficient matrices
,
,
,
, and
in Equation (44) are given in
Appendix A. In order to obtain the dispersion and attenuation characteristics of the shear horizontal wave, define another vector as
Inserting Equation (46) into Equation (44) leads to
Combining Equations (46) and (47) leads to
where
and
are the null and unit matrices, respectively. Equation (48) is a generalized eigenvalue problem for the apparent wave number component
. Given the numerical values of the angular frequency of the shear horizontal wave
, the dispersion and attenuation curves of the shear horizontal wave can be obtained. Searching for complex roots of complex transcendental equations in our previous mathematical model established through the state transfer equation method can be avoided mathematically. Since many elements of matrices
,
,
,
, and
are related to the definite integrals of Legendre and Laguerre polynomials, the analytical results of the definite integrals of the above two polynomials are given in the
Appendix A of this paper. Using these analytic results, the generalized eigenvalue problem in Equation (48) can be calculated efficiently. However, to obtain the dispersion and attenuation curves of the shear horizontal wave, the interface characteristic length
should be legitimately confirmed. To legitimately confirm the interface characteristic length
, another equivalent mathematical model is established in
Section 4, i.e., the second mathematical model.
4. The Second Mathematical Model
In this mathematical model, the Schottky junction is theoretically treated as an electrical gradient layer, in which the steady electron carrier concentration is a function of the space coordinate
and the vertical biasing electric field
mathematically, as shown in Equation (2). Therefore, Equation (21) are modified as
It is observed that the integral range of Equation (49) is different from that of Equation (21), i.e., the optimized spectral method. The mathematical reason why we set the integral range of Equation (49) this result is that the steady electron carrier concentration
is constant in this integral range according to Equation (2), and such mathematical operation can simplify subsequent numerical calculations of definite integrals and improve the computational efficiency of the spectral method. At the boundary
, consider the following interface condition [
10]
Inserting Equations (16), (17), (22) and (23) into Equations (52) leads to
Inserting Equations (52)–(54) into Equations (19), (20) and (49) leads to
Combining Equations (40) and (55)–(57) and letting
changes from
to
leads to
where the unknown expansion coefficients vector is
The explicit expressions of coefficient matrices
,
,
,
, and
in Equation (59) are given in
Appendix A. In order to obtain the dispersion and attenuation characteristics of the shear horizontal wave, define another one vector as
Inserting Equation (60) into Equation (58) leads to
Combining Equations (60) and (61) leads to
Equation (62) is also a generalized eigenvalue problem for the apparent wave number component
without the interface characteristic length
. Considering these two mathematical models established in
Section 3 and
Section 4 are equivalent, the interface characteristic length
can be inversely confirmed through the comparison of dispersion and attenuation curves numerically calculated using these two mathematical models in the next section.
5. Numerical Results and Discussion
In this section, the dispersion and attenuation curves of shear horizontal waves are numerically calculated using the mathematical models established in
Section 3 and
Section 4. Meanwhile, the interface characteristic length
are legitimately confirmed through the comparison of dispersion and attenuation curves calculated using these two equivalent mathematical models. The upper metal layer and lower n-type PSC substrate are defined as Au and ZnO, respectively. The constitutive parameters of these two materials in the Cartesian coordinate system are contained in
Table 1 [
10]. To study the variation law of the interface characteristic length
, the thickness of the upper Au layer is taken as
,
or
, respectively, and the doping concentration of the lower n-type ZnO substrate is
,
, or
, respectively. Since the numerical calculation of the dispersion and attenuation curves of shear horizontal waves is transformed into solving the eigenvalues, many eigenvalues that do not belong to the dispersion and attenuation curves of shear horizontal waves are also obtained. To solve the above computational complexity of the spectral method, we set the above numerical calculation examples to compare with our previous research work and to quickly refine the eigenvalues that belong to the dispersion and attenuation curves of shear horizontal waves from the original calculation dates. The horizontal axes of the dispersion and attenuation curves are the dimensionless angular frequency
. The vertical axis of the dispersion curve is the wave speed of shear horizontal waves
. The vertical axis of the attenuation curve is the dimensionless attenuation parameter
.
To verify the correctness and superiority of the research work in this paper,
Figure 2 shows the dispersion and attenuation curves of shear horizontal waves calculated by the second mathematical model without the consideration of biasing electric fields, when the thickness of the upper Au layer is
and the doping concentration of lower n-type ZnO substrate is
,
, or
, respectively. Similar to
Figure 2 and
Figure 3 in Reference [
16], it is also observed from
Figure 2a,b that the phase velocities and dimensionless attenuation parameters of the first and second orders shear horizontal wave tend to
and
, respectively, with the increase of the dimensionless angular frequency; the phase velocities and dimensionless attenuation parameters of the first and second orders shear horizontal wave both gradually increase with the increase of the dimensionless angular frequency. At different doping concentrations, the dispersion curves of shear horizontal waves coincide approximately, which means that the Schottky junction has negligible influence on the dispersion characteristic of shear horizontal waves; but the attenuation curves of shear horizontal waves change obviously, which means that the Schottky junction has evident influence on the attenuation characteristic of shear horizontal wave. Therefore, the Schottky junction can be theoretically treated as a two-dimensional electrically imperfect interface. However, different from
Figure 2 shown in Reference [
16], the phase velocities and dimensionless attenuation parameters of the first order shear horizontal waves don’t tend to
and
, respectively, when
. The reason for these numerical calculation results is that the second mathematical model is established using the optimized spectral method. On the one hand, searching for complex roots of complex transcendental equations is transformed into eigenvalue calculation problems; the initial values which are closely related to the complex roots and need to be given in advance are no longer given. On the other hand, the numerical calculation of the transfer matrix of Schottky junction through the truncated Magnus series in Reference [
16] is also avoided as the integral interval of Equation (49) is
. All the above mathematical operations can reduce the numerical calculation errors and improve the computational efficiency. The robustness of the second mathematical model is better than that of the mathematical model in Reference [
16]. So, we can calculate and plot the wider dimensionless angular frequency intervals and higher orders of the dispersion and attenuation curves of shear horizontal waves in this paper.
Since the Schottky junction has evident influence on the attenuation characteristic of shear horizontal wave, the comparisons of the attenuation curves of shear horizontal waves are utilized to legitimately confirm the numerical value of the interface piezoelectric parameter.
Figure 3 shows the comparisons of the attenuation curves of shear horizontal waves calculated by the first and second mathematical models without the consideration of biasing electric fields, when the thickness of the upper Au layer is
,
, or
, respectively, and the doping concentration of lower n-type ZnO substrate is
or
, respectively. The discrete numerical points are calculated using the first mathematical model, in which the interface characteristic length related to the interface piezoelectric parameter is
Because the Schottky junction is theoretically treated as a two-dimensional electrically imperfect interface, there is only one interface characteristic length that needs to be determined. This provides great convenience for the numerical calculations in this section. It is noted from Equation (63) that the numerical value of
is negative. This numerical calculation result can be explained by Equations (32) that the sum of the electric displacement of the shear horizontal wave at the boundary
and the gradient of the interface electric displacement is exactly equal to 0, which is the electric displacement of the shear horizontal wave at the boundary
. The continuous curves are also calculated using the second mathematical model. Since these two mathematical models are equivalent, all the discrete numerical points fall on the continuous curves approximately. It is further observed from
Figure 3 that the numerical value of
is independent of the thickness of the upper Au layer and the doping concentration of the lower n-type ZnO substrate. This is because the thickness of the Schottky junction is usually on the order of
to the power of
, which is much smaller than the thickness of the upper Au layer. Due to the upper metal layer being a grounded Au, the electric potential in the whole upper Au layer is zero, which causes the interface electric potential to be zero as well and there is no interface boundary condition related to the semiconductor effect, as shown in Equation (32). Since the interface characteristic length
is a constant, the first mathematical model is more practical, and can be utilized to investigate the interaction between high-frequency elastic waves and Schottky junctions.
To investigate the influence of biasing electric fields, the following dimensionless biasing electric fields are introduced, namely,
Figure 4 and
Figure 5 show the dispersion and attenuation curves of shear horizontal waves with the variation of the above two dimensionless biasing electric fields. It is observed from
Figure 4 and
Figure 5 that these two biasing electric fields
and
have almost no influence on the dispersion characteristic of shear horizontal waves. This numerical calculation result can be explained by Equation (7) that these two biasing electric fields are only related to the carrier’s current density components, and the electrical energy accounts for a relatively low proportion of the total energy of shear horizontal waves. While the increase of
makes the attenuation curves of shear horizontal waves shift towards higher values, which means that the existence of
can greatly enhance the attenuation characteristics of shear horizontal waves. But the biasing electric field
has negligible influence on the attenuation characteristic of shear horizontal waves. This is because the direction of the biased electric field
is the same as the propagation direction of shear horizontal waves, as shown in
Figure 1.
Figure 6 shows the comparisons of attenuation curves of shear horizontal waves calculated by the first and second mathematical models with the consideration of biasing electric fields, when the thickness of the upper Au layer is
,
, or
, respectively, and the doping concentration of lower n-type ZnO substrate is
or
, respectively. The discrete numerical points and continuous curves are also calculated using the first and second mathematical models, respectively. It is also observed from
Figure 6 that all the discrete numerical points fall on the continuous curves approximately, which means that the numerical value of
in Equation (63) is independent of biasing electric fields. This numerical calculation result can also be explained by Equation (32) that there is no biasing electric field in the electrically imperfect interface boundary condition.