Next Article in Journal
On K-Banhatti, Revan Indices and Entropy Measures of MgO(111) Nanosheets via Linear Regression
Previous Article in Journal
A Low-Inertia and High-Stiffness Cable-Driven Biped Robot: Design, Modeling, and Control
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Influence of the Schottky Junction on the Propagation Characteristics of Shear Horizontal Waves in a Piezoelectric Semiconductor Semi-Infinite Medium

1
School of Mathematics and Physics, University of Science and Technology Beijing, Beijing 100083, China
2
Department of Applied Mechanics, University of Science and Technology Beijing, Beijing 100083, China
3
School of Civil and Resource Engineering, University of Science and Technology Beijing, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(4), 560; https://doi.org/10.3390/math12040560
Submission received: 26 December 2023 / Revised: 22 January 2024 / Accepted: 5 February 2024 / Published: 13 February 2024
(This article belongs to the Special Issue Mathematical Analysis in Materials Science and Spectral Theory)

Abstract

:
In this paper, a theoretical model of the propagation of a shear horizontal wave in a piezoelectric semiconductor semi-infinite medium is established using the optimized spectral method. First, the basic equations of the piezoelectric semiconductor semi-infinite medium are derived with the consideration of biased electric fields. Then, considering the propagation of a shear horizontal wave in the piezoelectric semiconductor semi-infinite medium, two equivalent mathematical models are established. In the first mathematical model, the Schottky junction is theoretically treated as an electrically imperfect interface, and an interface characteristic length is utilized to describe the interface effect of the Schottky junction. To legitimately confirm the interface characteristic length, a second mathematical model is established, in which the Schottky junction is theoretically treated as an electrical gradient layer. Finally, the dispersion and attenuation curves of shear horizontal waves are numerically calculated using these two mathematical models to discuss the influence of the Schottky junction on the dispersion and attenuation characteristics of shear horizontal waves. Utilizing the equivalence of these two mathematical models and the above numerical results, the numerical value of the interface characteristic length is reliably legitimately confirmed; this value is independent of the thickness of the upper metal layer, the doping concentration of the lower n-type piezoelectric semiconductor substrate, and biasing electric fields. Only the biasing electric field parallel to the Schottky junction can provide an evident influence on the attenuation characteristics of shear horizontal waves and enhance the interface effect of the Schottky junction. Since the second mathematical model is also a validation of our previous mathematical model established through the state transfer equation method, some numerical results calculated using these two mathematical models are compared with those obtained using the previous method to verify the correctness and superiority of the research work presented in this paper. Since these two mathematical models can better calculate the dispersion and attenuation curves of high-frequency waves in micro- and nano-scale piezoelectric semiconductor materials, the establishment of mathematical models and the revelation of physical mechanisms are fundamental to the analysis and optimization of micro-scale resonators, energy harvesters, and amplifications.

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 O x y plane. The z -direction is the polarization direction. The thickness of the upper grounded metal layer H 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 y axis. The geometrical thickness of the Schottky junction is y n , 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 E ¯ x and E ¯ y , respectively. Consider the propagation of the shear horizontal wave in this PSC substrate in the O x y plane along the x 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 ‘ U ’ and ‘ L ’, 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
σ i j L = c i j k l L u l , k L + e k i j L ϕ , k L , D k L = e k i j L u j , i L ε k l L ϕ , l L ,
J k L = q n ¯ L μ k l L ϕ , l L + q n L μ k l L E ¯ l + q d k l L n , l L
where u l L , ϕ L , n L , σ i j L , D k L , and J k L 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, c i j k l L , e k i j L , ε k l L , μ k l L , and d k l L are the elastic, piezoelectric, dielectric, electron carrier migration, and electron carrier diffusion parameter tensors, respectively. The carrier charge is q = 1.602 × 10 19 C . The steady electron carrier concentration is n ¯ L , which is a function of the space coordinate y and the vertical biasing electric field E ¯ y mathematically
n ¯ L = N D L e x p μ y y L d y y L q N D L e y y z L 2 / c y z y z L + ε y y L 1 2 y 2 y n y 1 2 y n 2 2 E ¯ y y + y n y n y 0 N D L < y y n
where N D L are the donor concentration i.e., the doping concentration, in the lower n-type PSC substrate. The geometrical thickness of the Schottky junction is
y n = 2 e y y z L 2 / c y z y z L + ε y y L W U W L q 2 N D L
where W U and W L 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 x and y and time coordinate t mathematically
u x L = u y L = 0 , u z L = u z L x , y , t , ϕ L = ϕ L x , y , t , n L = n L x , y , t
Inserting Equation (4) into Equation (1) leads to the Cauchy stress components
σ z x L = σ x z L = c x z x z L u z L x + e y y z L ϕ L x , σ z y L = σ y z L = c y z y z L u z L y + e y y z L ϕ L y
electric displacement components
D x L = e x x z L u z L x ε x x L ϕ L x , D y L = e y y z L u z L y ε y y L ϕ L y
and electron carrier current density components
J x L = q n ¯ L μ x x L ϕ L x + q n L μ x x L E ¯ x + q d x x L n L x ,
J y L = q n ¯ L μ y y L ϕ L y + q n L μ y y L E ¯ y + q d y y L n L y
With the consideration of the piezoelectric and semiconductor effects, the governing equations of the lower n-type PSC substrate [11,12] are
σ i j , i L = ρ L 2 u j L t 2 , D k , k L = q n L , J k , k L = q n L t
where ρ L 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
c y z y z L 2 x 2 + 2 y 2 u z L + e y y z L 2 x 2 + 2 y 2 ϕ L = ρ 2 u z L t 2 e y y z L 2 x 2 + 2 y 2 u z L ε y y L 2 x 2 + 2 y 2 ϕ L = q n L n ¯ L μ y y L ϕ L y n ¯ L μ y y L 2 x 2 + 2 y 2 ϕ L + μ y y L E ¯ x x + μ y y L E ¯ y y n L + d y y L 2 x 2 + 2 y 2 n L = n L t
Without the piezoelectric and semiconductor effects, under the supposition of small deformations, the constitutive equations of the upper metal layer [28] are
σ i j U = c i j k l U u l , k U
where u j U , σ i j U , and c i j k l U 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 u j U is also a function of space coordinates x and y and time coordinate t mathematically
u x U = u y U = 0 , u z U = u z U x , y , t
Inserting Equation (11) into Equation (10) leads to the Cauchy stress components of the upper metal layer
σ z x U = σ x z U = c x z x z U u z U x , σ z y U = σ y z U = c y z y z U u z U y
The governing equation of the upper metal layer [28] is
σ i j , i U = ρ U 2 u j U t 2
where ρ U is the mass density of the upper metal layer. Inserting Equation (12) into Equation (13) leads to
c y z y z U 2 x 2 + 2 y 2 u z U = ρ U 2 u z U t 2
On the surface y = H , consider the surface condition of the upper metal layer with respect to the Cauchy stress [28]
σ y z U y = H = 0
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
u z L x , y , t , ϕ L x , y , t , n L x , y , t = U z L y , U ϕ L y , U n L y e x p k L x ω L t , < y 0
where k L and ω L 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 Q m are used to expand these three undetermined amplitudes U z L , U ϕ L , and U n L , which are all functions of the space coordinate y mathematically
U z L y = m = 0 z m L R m L y m = 0 M z m L R m L y ,
U ϕ L y = m = 0 ϕ m L R m L y m = 0 M ϕ m L R m L y ,
U n L y = m = 0 n m L R m L y m = 0 M n m L R m L y , < y 0
where z m L , ϕ m L and n m L are expansion coefficients and
R m L y = e y 2 Q m y , R m L = R l L = 0 , Q m y ¯ = e y ¯ m ! d m d y ¯ m y ¯ m e y ¯
Inserting Equations (16) and (17) into Equation (9), multiplying both sides by R l L y and taking the definite integrals lead to
ω 2 ρ L m = 0 M z m L 0 R m L y R l L y d y + c y z y z L m = 0 M z m L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y + e y y z L m = 0 M ϕ m L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y = k 2 c y z y z L m = 0 M z m L 0 R m L y R l L y d y + k 2 e y y z L m = 0 M ϕ m L 0 R m L y R l L y d y
e y y z L m = 0 M z m L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y ε y y L m = 0 M ϕ m L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y + q m = 0 M n m L 0 R m L y R l L y d y = k 2 e y y z L m = 0 M z m L 0 R m L y R l L y d y k 2 ε y y L m = 0 M ϕ m L 0 R m L y R l L y d y
i ω m = 0 M n m L 0 R m L y R l L y d y N D L μ y y L m = 0 M ϕ m L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y + d y y L m = 0 M n m L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y + i k E ¯ x μ y y L m = 0 M n m L 0 R m L y R l L y d y + E ¯ y μ y y L m = 0 M n m L R m L 0 R l L 0 0 R m L y R l L y d y = N D L μ y y L k 2 m = 0 M ϕ m L 0 R m L y R l L y d y + d y y L k 2 m = 0 M n m L 0 R m L y R l L y d y
According to Equation (11), assume the mathematical expressions of the displacement of the shear horizontal wave in the upper metal layer as
u z U x , y , t = U z U y e x p i k U x ω U t , 0 y H
where k U and ω U 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 P m are used to expand the undetermined amplitude U z U , which is also a function of the space coordinate y mathematically
U z U y = m = 0 z m U R m U y m = 0 M z m U R m U y , 0 y H
where z m U is expansion coefficients, m = 0,1 , 2 , , M . M is an integer.
R m U y = P m 2 y H H , P m y ¯ = d m 2 m m ! d y ¯ m y ¯ 2 1 m
Inserting Equations (22) and (23) into Equation (14), multiplying both sides by R l U y and taking the definite integrals leads to
ρ U ω 2 m = 0 M z m U 0 H R m U y R l U y d y + c y z y z U m = 0 M z m U R m U y R l U y 0 H R m U y R l U y 0 H + 0 H R m U y R l U y d y = c y z y z U k 2 m = 0 M z m U 0 H R m U y R l U y d y
Considering the geometrical thickness of the Schottky junction is much smaller than that of the upper metal layer, i.e., y n H , 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
n i σ i j L , n i D i L y = 0 n i σ i j U , 0 y = 0 + = 0 , D l , l I
where D l I is the electric displacement vector of the electrically imperfect interface,
D l I = e l i j I u j , i I ε l j I ϕ , j I
and e l i j I and ε l j I are the interface piezoelectric and dielectric parameter tensors, respectively. The interface boundary condition to the displacement and electric potential at the boundary y = 0 is
u z L , ϕ L y = 0 = u z U , 0 y = 0 + = u z I , ϕ I
where u z I and ϕ I are the displacement vector and electric potential of the electrically imperfect interface, respectively. Inserting the interface normal vector n i = 0,1 , 0 into Equation (26) leads to
σ y z L , D y L y = 0 σ y z U , 0 y = 0 + = 0 , D x , x I
Inserting Equation (28) into Equation (27) leads to
D x I = e y y z I u z , x I
and
e y y z I = f e L e y y z L
where f e L is the interface characteristic length related to the interface piezoelectric parameter. Inserting Equations (30) and (31) into Equations (28) and (29) leads to
u z L , ϕ L , σ y z L y = 0 = u z U , 0 , σ y z U y = 0 +
e y y z L u z L y ε y y L ϕ L y y = 0 = f e L e y y z L 2 u z L x 2 y = 0
Inserting Equations (16), (17), (22), and (23) into Equations (15), (32) and (33) leads to
m = 0 M z m U R m U H = 0
m = 0 M z m U R m U 0 = m = 0 M z m L R m L 0
m = 0 M ϕ m L R m L 0 = 0
c y z y z L m = 0 M z m L R m L 0 + e y y z L m = 0 M ϕ m L R m L 0 = c y z y z U m = 0 M z m U R m U 0
e y y z L m = 0 M z m L R m L 0 ε y y L m = 0 M ϕ m L R m L 0 = k 2 f e L e y y z L m = 0 M z m L R m L 0
and
k U = k L = k ,   ω U = ω L = ω
Inserting Equations (34) and (35) into Equation (25) leads to
ρ U ω 2 m = 0 M z m U 0 H R m U y R l U y d y
c y z y z U m = 0 M z m U R m U 0 R l U 0 + R m U H R l U H 0 H R m U y R l U y d y
+ c y z y z U m = 0 M z m L R m L 0 R l U 0 = k 2 c y z y z U m = 0 M z m L 0 H R m U y R l U y d y
Inserting Equations (36)–(38) into Equations (19)–(21) leads to
ω 2 ρ L m = 0 M z m L 0 R m L y R l L y d y c y z y z L m = 0 M z m L R m L 0 R l L 0 0 R m L y R l L y d y + e y y z L m = 0 M ϕ m L 0 R m L y R l L y d y + c y z y z U m = 0 M z m U R m U 0 R l L 0 = k 2 c y z y z L m = 0 M z m L 0 R m L y R l L y d y + k 2 e y y z L m = 0 M ϕ m L 0 R m L y R l L y d y
e y y z L m = 0 M z m L 0 R m L y R l L y d y R m L 0 R l L 0 ε y y L m = 0 M ϕ m L 0 R m L y R l L y d y + q m = 0 M n m L 0 R m L y R l L y d y = k 2 e y y z L m = 0 M z m L 0 R m L y R l L y d y k 2 ε y y L m = 0 M ϕ m L 0 R m L y R l L y d y + k 2 f e L e y y z L m = 0 M z m L R m L 0 R l L 0
i ω m = 0 M n m L 0 R m L y R l L y d y N D L μ y y L m = 0 M ϕ m L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y + d y y L m = 0 M n m L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y + E ¯ y μ y y L m = 0 M n m L R m L 0 R l L 0 0 R m L y R l L y d y + i k E ¯ x μ y y L m = 0 M n m L 0 R m L y R l L y d y = N D L μ y y L k 2 m = 0 M ϕ m L 0 R m L y R l L y d y + d y y L k 2 m = 0 M n m L 0 R m L y R l L y d y
Combining Equations (40)–(43) and letting l changes from 0 to M leads to
ω 2 A P + i ω B P + C P + i k F P = k 2 D P
where the unknown expansion coefficients vector is
P = z 0 U , z 1 U , , z M U , z 0 L , z 1 L , , z M L , ϕ 0 L , ϕ 1 L , , ϕ M L , n 0 L , n 1 L , , n M L T
The explicit expressions of coefficient matrices A , B , C , D , and F 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
Q = k P
Inserting Equation (46) into Equation (44) leads to
ω 2 A P + i ω B P + C P = i k F P + k D Q
Combining Equations (46) and (47) leads to
O E ω 2 A + i ω B + C O P Q = k E O i F D P Q
where O and E are the null and unit matrices, respectively. Equation (48) is a generalized eigenvalue problem for the apparent wave number component k . 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 A , B , C , D , and F 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 f e L should be legitimately confirmed. To legitimately confirm the interface characteristic length f e L , 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 y and the vertical biasing electric field E ¯ y mathematically, as shown in Equation (2). Therefore, Equation (21) are modified as
i ω m = 0 M n m L y n R m L y R l L y d y μ y y L N D L m = 0 M ϕ m L y n R m L y R l L y d y + d y y L m = 0 M n m L R m L 0 R l L 0 R m L 0 R l L 0 + y n R m L y R l L y d y d y y L m = 0 M n m L y n 0 R m L y R l L y d y + μ y y L m = 0 M n m L y n i k E ¯ x R m L y + E ¯ y R m L y R l L y d y = k 2 μ y y L N D L m = 0 M ϕ m L y n R m L y R l L y d y + k 2 d y y L m = 0 M n m L y n R m L y R l L y d y
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 n ¯ L 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 y = 0 , consider the following interface condition [10]
u z U , 0,0 , σ y z U y = 0 + = u z L , ϕ L , n L , σ y z L y = 0
Inserting Equations (16), (17), (22) and (23) into Equations (52) leads to
m = 0 M z m L R m L 0 = m = 0 M z m U R m U 0
m = 0 M ϕ m L R m L 0 = 0
m = 0 M n m L R m L 0 = 0
c y z y z L m = 0 M z m L R m L 0 + e y y z L m = 0 M ϕ m L R m L 0 = c y z y z U m = 0 M z m U R m U 0
Inserting Equations (52)–(54) into Equations (19), (20) and (49) leads to
ω 2 ρ L m = 0 M z m L 0 R m L y R l L y d y c y z y z L m = 0 M z m L R m L 0 R l L 0 0 R m L y R l L y d y + e y y z L m = 0 M ϕ m L 0 R m L y R l L y d y + c y z y z U m = 0 M z m U R m U 0 R l L 0 = k 2 c y z y z L m = 0 M z m L 0 R m L y R l L y d y + k 2 e y y z L m = 0 M ϕ m L 0 R m L y R l L y d y
e y y z L m = 0 M z m L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y ε y y L m = 0 M ϕ m L R m L 0 R l L 0 + 0 R m L y R l L y d y + q m = 0 M n m L 0 R m L y R l L y d y = k 2 e y y z L m = 0 M z m L 0 R m L y R l L y d y k 2 ε y y L m = 0 M ϕ m L 0 R m L y R l L y d y
i ω m = 0 M n m L y n R m L y R l L y d y μ y y L N D L m = 0 M ϕ m L y n R m L y R l L y d y + d y y L m = 0 M n m L R m L 0 R l L 0 + 0 R m L y R l L y d y d y y L m = 0 M n m L y n 0 R m L y R l L y d y + μ y y L m = 0 M n m L y n i k E ¯ x R m L y + E ¯ y R m L y R l L y d y = k 2 μ y y L N D L m = 0 M ϕ m L y n R m L y R l L y d y + k 2 d y y L m = 0 M n m L y n R m L y R l L y d y
Combining Equations (40) and (55)–(57) and letting l changes from 0 to M leads to
ω 2 A ¯ P + i ω B ¯ P + C ¯ P + i k F ¯ P = k 2 D ¯ P a = 1 ,
where the unknown expansion coefficients vector is
P = z 0 U , z 1 U , , z M U , z 0 L , z 1 L , , z M L , ϕ 0 L , ϕ 1 L , , ϕ M L , n 0 L , n 1 L , , n M L T
The explicit expressions of coefficient matrices A ¯ , B ¯ , C ¯ , D ¯ , and F ¯ 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
Q = k P
Inserting Equation (60) into Equation (58) leads to
ω 2 A ¯ P + i ω B ¯ P + C ¯ P = i k F ¯ P + k D ¯ Q
Combining Equations (60) and (61) leads to
O E ω 2 A ¯ + i ω B ¯ + C ¯ O P Q = k E O i F ¯ D ¯ P Q
Equation (62) is also a generalized eigenvalue problem for the apparent wave number component k without the interface characteristic length f e L . Considering these two mathematical models established in Section 3 and Section 4 are equivalent, the interface characteristic length f e L 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 f e L 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 f e L , the thickness of the upper Au layer is taken as 1   μ m , 5   μ m or 10   μ m , respectively, and the doping concentration of the lower n-type ZnO substrate is 1 × 10 21   m 3 , 5 × 10 21   m 3 , or 10 × 10 21   m 3 , 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 ω ¯ = ω H / c y z y z L / ρ L . The vertical axis of the dispersion curve is the wave speed of shear horizontal waves c = r e a l ω / k . The vertical axis of the attenuation curve is the dimensionless attenuation parameter k ¯ = i m a g k H .
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 1   μ m and the doping concentration of lower n-type ZnO substrate is 1 × 10 21   m 3 , 5 × 10 21   m 3 , or 10 × 10 21   m 3 , 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 1475   m / s and 0   m / s , 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 2747   m / s and 0   m / s , respectively, when ω ¯ 0 . 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 < y y n . 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 1   μ m , 5   μ m , or 10   μ m , respectively, and the doping concentration of lower n-type ZnO substrate is 1 × 10 21   m 3 or 10 × 10 21   m 3 , 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
f e L = 0.014774258   m
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 f e L 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 y = 0 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 y = 0 + . 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 f e L 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 10 to the power of 7 , 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 f e L 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,
γ x = μ y y L E ¯ x ρ L / c y z y z L , γ y = μ y y L E ¯ y ρ L / c y z y z L
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 E ¯ x and E ¯ y 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 E ¯ x makes the attenuation curves of shear horizontal waves shift towards higher values, which means that the existence of E ¯ x can greatly enhance the attenuation characteristics of shear horizontal waves. But the biasing electric field E ¯ y has negligible influence on the attenuation characteristic of shear horizontal waves. This is because the direction of the biased electric field E ¯ x 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 1   μ m , 5   μ m , or 10   μ m , respectively, and the doping concentration of lower n-type ZnO substrate is 1 × 10 21 m 3 or 10 × 10 21   m 3 , 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 f e L 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.

6. Concluding Remarks

The influence of the Schottky junction on the dispersion and attenuation characteristics of shear horizontal waves in an n-type PSC substrate with a metal layer is the primary concerns of the present work. According to the surface and interface piezoelectric theories, a new interface effect theoretical model of the Schottky junction is derived based on the physical properties of the Schottky junction. An interface characteristic length that related to the interface piezoelectric parameter is utilized to express the electrically imperfect interface boundary conditions of the Schottky junction. To legitimately confirm the interface characteristic length, the Schottky junction is treated as an electrical gradient layer. Two equivalent mathematical models are established through the optimized spectral method. From the dispersion and attenuation curves of shear horizontal waves and the interface characteristic length calculated and confirmed through these two equivalent mathematical models, the following conclusions can be drawn:
  • Compared with our previous numerical results, the more accurate, wider dimensionless angular frequency intervals and higher orders of the dispersion and attenuation curves of shear horizontal waves can be calculated through these two equivalent mathematical models, which embodies the value of this paper’s research work.
  • The second mathematical model is more accurate, but the computational stability and engineering practicality of the first mathematical model is better than that of the second mathematical model.
  • Only the biasing electric field parallel to the Schottky junction can provide evident influence on the attenuation characteristic of shear horizontal waves and enhance the influence of the Schottky junction; adjusting the biasing electric field can regulate attenuation characteristics of shear horizontal waves and optimize the acoustic characteristics of energy harvesters formed by PSC semi-infinite medium
Since the influence of the Schottky junction mainly concentrates on the attenuation characteristic of the shear horizontal wave, the discussion in this paper provides theoretical guidance for the design and manufacture of SAW resonators, energy harvesters, and acoustic wave amplifications [29,30]. Through the first mathematical model and the interface characteristic length, we might consider the propagation characteristics of shear horizontal waves in more complex PSC materials in subsequent work, such as the PSC layered composite structures containing multi-layer Schottky junctions.

Author Contributions

Conceptualization, X.G.; methodology, X.G.; software, X.G.; validation, X.G.; formal analysis, Y.W.; investigation, Y.W.; resources, C.X.; data curation, C.X.; writing—original draft preparation, Z.W.; writing—review and editing, Z.W.; visualization, C.D.; supervision, C.D.; All authors have read and agreed to the published version of the manuscript.

Funding

The whole research is supported by the National Natural Science Foundation of China (12202039) and the Interdisciplinary Research Project for Young Teachers of USTB (Fundamental Research Funds for the Central Universities) (FRF-IDRY-21-006).

Data Availability Statement

The data that support the findings of this study are available on request from the corresponding author, Xiao Guo ([email protected]), upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest in this work. The authors declare that they do not have any commercial or associative interest that represents a conflict of interest in connection with the work submitted.

Appendix A

To obtain the geometrical thickness and steady electron carrier concentration of the Schottky junction, consider the constitutive equations of the Schottky junction
σ y z L y = c y z y z L u z , y L y + e y y z L ϕ , y L y
D y L y = e y y z L u z , y L y ε y y L ϕ , y L y
J y L y = q n ¯ L y μ y y L d ϕ L y d y + q n ¯ L y μ y y L E ¯ y + q d y y L d n ¯ L y d y , y n y 0
Consider the Cauchy stress tensor in the Schottky junction is zero
σ y z L y = 0 , y n y 0
Inserting Equation (A1) into Equation (A4) leads to
u z , y L y = e y y z L c y z y z L ϕ , y L y , y n y 0
Inserting Equation (A5) into Equation (A2) leads to
D y L y = e y y z L 2 c y z y z L + ε y y L ϕ , y L y , y n y 0
Consider the governing equations of the Schottky junction
D y , y L = q N D L , y n y 0
where N D L are the donor concentration i.e., the doping concentration, in the lower n-type PSC substrate. Inserting Equation (A7) into Equation (A8) leads to
ϕ , y y L y = q N D L e y y z L 2 / c y z y z L + ε y y L , y n y 0
The mathematical solution of ϕ L can be obtained from Equation (A8)
ϕ L y = q N D L e y y z L 2 / c y z y z L + ε y y L 1 2 y 2 + A y + B , y n y 0
where A and B are two unknown constants. To determine the expression of A , the potential boundary conditions at the lower boundary of the Schottky junction are introduced
ϕ , y L y n = E ¯ y
Combining Equations (A9) and (A10) leads to
ϕ L y = q N D L e y y z L 2 / c y z y z L + ε y y L 1 2 y 2 y n y + B E ¯ y y , y n y 0
The contact potential condition between the upper metal layer and lower n-type PSC substrate is introduced
ϕ L y n ϕ L 0 = W U W L q + E ¯ y y n
Inserting Equation (A1) into Equation (A12) leads to the geometrical thickness of the Schottky junction
y n = 2 e y y z L 2 / c y z y z L + ε y y L W U W L q 2 N D L
Consider the electron carrier current density tensor in the Schottky junction is zero
J y L y = 0 , y n y 0
Inserting Equations (A3) and (A11) into Equation (A14) leads to
n ¯ L y = C e x p μ y y L d y y L q N D L e y y z L 2 / c y z y z L + ε y y L 1 2 y 2 y n y 2 E ¯ y y , y n y 0
The boundary conditions of steady-state carrier concentration are introduced
n ¯ L y n = N D L
Combining Equations (A15) and (A16), and μ y y L / d y y L = q / k B T leads to the steady electron carrier concentration of the Schottky junction
n ¯ L = N D L e x p q k B T q N D L e y y z L 2 / c y z y z L + ε y y L 1 2 y 2 y n y 1 2 y n 2 2 E ¯ y y + y n y n y 0 N D L < y y n
The explicit expressions of coefficient matrices A , B , C , D , and F in Equation (44) are
A = A 11 l m O O O O A 22 l m O O O O O O O O O O , B = O O O O O O O O O O O O O O O B 44 l m ,
C = C 11 l m C 12 l m O O C 21 l m C 22 l m C 23 l m O O C 32 l m C 33 l m C 34 l m O O C 43 l m C 44 l m , D = D 11 l m O O O O D 22 l m D 23 l m O O D 32 l m D 33 l m O O O D 43 l m D 44 l m ,
F = O O O O O O O O O O O O O O O F 44 l m
where
A 11 l m = ρ U 0 H R m U y R l U y d y , A 22 l m = ρ L 0 R m L y R l L y d y ,
B 44 l m = 0 R m L y R l L y d y ,
C 11 l m = c y z y z U R m U 0 R l U 0 + R m U H R l U H 0 H R m U y R l U y d y ,
C 12 l m = c y z y z U R m L 0 R l U 0 , C 21 l m = c y z y z U R m U 0 R l L 0 ,
C 22 l m = c y z y z L R m L 0 R l L 0 0 R m L y R l L y d y ,
C 23 l m = e y y z L 0 R m L y R l L y d y ,
C 32 l m = e y y z L 0 R m L y R l L y d y R m L 0 R l L 0 ,
C 33 l m = ε y y L 0 R m L y R l L y d y , C 34 l m = q 0 R m L y R l L y d y ,
C 43 l m = N D L μ y y L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y ,
C 44 l m = d y y L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y
+ E ¯ y μ y y L R m L 0 R l L 0 0 R m L y R l L y d y ,
D 11 l m = c y z y z U 0 H R m U y R l U y d y , D 22 l m = c y z y z L 0 R m L y R l L y d y ,
D 23 l m = e y y z L 0 R m L y R l L y d y ,
D 32 l m = e y y z L 0 R m L y R l L y d y + f e L e y y z L R m L 0 R l L 0 ,
D 33 l m = ε y y L 0 R m L y R l L y d y , D 43 l m = N D L μ y y L 0 R m L y R l L y d y
D 44 l m = d y y L 0 R m L y R l L y d y , F 44 l m = E ¯ x μ y y L 0 R m L y R l L y d y
and
R m L 0 = 1 , R m L 0 = 1 2 + m , 0 R m L y R l L y d y = δ m , l ,
0 R m L y R l L y d y = 1 2 δ m , l + h 1 m l 1 ,
0 R m L y R l L y d y = 1 4 δ m , l + m l h 1 m l 1 ,
R m U H = 1 , R m U 0 = ( 1 ) m ,   R m U H = m ( m + 1 ) H ,
R m U 0 = ( 1 ) m + 1 m ( m + 1 ) H , 0 H R m U y R l U y d y = H 2 m + 1 δ m , l ,
0 H R m U y R l U y d y = k = 0 n 2 δ 2 k 1 , l m = 2 n k = 0 n 2 δ 2 k , l m = 2 n + 1 ,
0 H R m U y R l U y d y
= 2 m m + 1 h 2 min m , l · m i n m , l + 1 h m + l 0   m o d   2 0 m + l 1   m o d   2 ,
δ m , l = 0 m l 1 m = l , h 1 y = 0 y < 0 1 y 0
The explicit expressions of coefficient matrices A ¯ , B ¯ , C ¯ , D ¯ and F ¯ in Equation (58) are
A ¯ = A ¯ 11 l m O O O O A ¯ 22 l m O O O O O O O O O O , B ¯ = O O O O O O O O O O O O O O O B ¯ 44 l m ,
C ¯ = C ¯ 11 l m C ¯ 12 l m O O C ¯ 21 l m C ¯ 22 l m C ¯ 23 l m O O C ¯ 32 l m C ¯ 33 l m C ¯ 34 l m O O C ¯ 43 l m C ¯ 44 l m , D ¯ = D ¯ 11 l m O O O O D ¯ 22 l m D ¯ 23 l m O O D ¯ 32 l m D ¯ 33 l m O O O D ¯ 43 l m D ¯ 44 l m ,
F ¯ = O O O O O O O O O O O O O O O F ¯ 44 l m
where
A ¯ 11 l m = ρ U 0 H R m U y R l U y d y , A ¯ 22 l m = ρ L 0 R m L y R l L y d y ,
B ¯ 44 l m = y n R m y R l ( y ) d y ,
C ¯ 11 l m = c y z y z U R m U 0 R l U 0 + R m U H R l U H 0 H R m U y R l U y d y ,
C ¯ 12 l m = c y z y z U R m L 0 R l U 0 , C ¯ 21 l m = c y z y z U R m U 0 R l L 0 ,
C ¯ 22 l m = c y z y z L R m L 0 R l L 0 0 R m L y R l L y d y ,
C ¯ 23 l m = e y y z L 0 R m L y R l L y d y ,
C ¯ 32 l m = e y y z L R m L 0 R l L 0 R m L 0 R l L 0 + 0 R m L y R l L y d y ,
C ¯ 33 l m = ε y y L R m L 0 R l L 0 + 0 R m L y R l L y d y ,
C ¯ 34 l m = q 0 R m L y R l L y d y ,
C ¯ 43 l m = μ y y L N D L y n R m L y R l L y d y , C ¯ 44 l m = d y y L R m L 0 R l L 0 + 0 R m L y R l L y d y d y y L y n 0 R m L y R l L y d y + μ y y L E ¯ y y n R m L y R l y d y ,
D ¯ 11 l m = c y z y z U 0 H R m U y R l U y d y , D ¯ 22 l m = c y z y z L 0 R m L y R l L y d y ,
D ¯ 23 l m = e y y z L 0 R m L y R l L y d y , D ¯ 32 l m = e y y z L 0 R m L y R l L y d y ,
D ¯ 33 l m = ε y y L 0 R m L y R l L y d y , D ¯ 43 l m = μ y y L N D L y n R m L y R l L y d y ,
D ¯ 44 l m = d y y L y n R m L y R l L y d y , F ¯ 44 l m = μ y y L E ¯ x y n R m L y R l L y d y
Since the second mathematical model is also a validation of our previous mathematical model, Table A1 summarizes the current work with previous work.
Table A1. The summarization of the current work with previous work.
Table A1. The summarization of the current work with previous work.
The Current WorkReference [16]
The First Mathematical ModelThe Second Mathematical Model
The constitutive equationsThe lower n-type PSC substrateWith the consideration of biasing electric fieldsWithout the consideration of biasing electric fields
The upper metal layerIdentical
The governing equationsThe lower n-type PSC substrateIdentical
The upper metal layerIdentical
The surface condition of the upper metal layerIdentical
Theoretical treatment of the Schottky junctionAn electrically imperfect interfaceAn electrical gradient layer
Theoretical modeling methodThe spectral methodThe optimized spectral methodThe state transfer equation method
The numerical calculation of dispersion and attenuation curvesGeneralized eigenvalue problemsSearching for complex roots of complex transcendental equations

References

  1. Cao, X.; Hu, S.; Liu, J.; Shi, J. Generalized Rayleigh surface waves in a piezoelectric semiconductor substrate. Meccanica 2019, 54, 271–281. [Google Scholar] [CrossRef]
  2. Liu, Y.; Lin, S.; Li, Y.; Li, C.; Liang, Y. Numerical investigation of shear horizontal wave in layered composite piezoelectric structures using the SIGA-PML approach. Compos. Part B 2019, 158, 230–238. [Google Scholar] [CrossRef]
  3. Singh, A.K.; Kumari, R.; Ray, A.; Chattopadhyay, A. Love-type waves in a piezoelectric-viscoelastic bimaterial composite structure due to an impulsive point source. Int. J. Mech. Sci. 2019, 152, 613–629. [Google Scholar] [CrossRef]
  4. Goyal, R.; Kumar, S. Quantifying viscoelastic, piezoelectric and couple stress effects on Love-type wave propagation. Smart Mater. Struct. 2019, 28, 105021. [Google Scholar] [CrossRef]
  5. Yang, G.; Zhang, M.; Hu, J.; Huang, B.; Xu, M.; Zhu, Z.; Du, J. Effects of initial stress on band gap of Love waves in a layered domain inverted phononic crystal structure. Ultrasonics 2020, 106, 106145. [Google Scholar] [CrossRef] [PubMed]
  6. Luo, Y.; Zhang, C.; Chen, W.; Yang, J. Piezotronic Effect of a Thin Film With Elastic and Piezoelectric Semiconductor Layers Under a Static Flexural Loading. J. Appl. Mech. 2019, 86, 051003. [Google Scholar] [CrossRef]
  7. Liang, Y.; Hu, Y. Influence of Doping Concentration on the Outputs of a Bent ZnO Nanowire. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2019, 66, 1793–1797. [Google Scholar] [CrossRef]
  8. Qu, Y.; Jin, F.; Yang, J. Torsion of a piezoelectric semiconductor rod of cubic crystals with consideration of warping and in-plane shear of its rectangular cross section. Mech. Mater. 2022, 172, 104407. [Google Scholar] [CrossRef]
  9. He, J.; Ye, W. Fully-coupled piezoelectric field in the semiconductor heterostructure with an ellipsoidal anisotropic inhomogeneity. Mech. Mater. 2023, 187, 104837. [Google Scholar] [CrossRef]
  10. Tian, R.; Nie, G.; Liu, J.; Pan, E.; Wang, Y. On Rayleigh waves in a piezoelectric semiconductor thin film over an elastic half-space. Int. J. Mech. Sci. 2021, 204, 106565. [Google Scholar] [CrossRef]
  11. Jiao, F.; Wei, P.; Zhou, X.; Zhou, X. The dispersion and attenuation of the multi-physical fields coupled waves in a piezoelectric semiconductor. Ultrasonics 2019, 92, 68–78. [Google Scholar] [CrossRef]
  12. Jiao, F.; Wei, P.; Zhou, Y.; Zhou, X. Wave propagation through a piezoelectric semiconductor slab sandwiched by two piezoelectric half-spaces. Eur. J. Mech./A Solids 2019, 75, 70–81. [Google Scholar] [CrossRef]
  13. Tian, R.; Liu, J.; Pan, E.; Wang, Y. SH waves in multilayered piezoelectric semiconductor plates with imperfect interfaces. Eur. J. Mech./A Solids 2020, 81, 103961. [Google Scholar] [CrossRef]
  14. Al-Ahmadi, N.A. Schottky barrier inhomogeneities at the interface of different epitaxial layer thicknesses of n-GaAs/Ti/Au/Si: Al0.33Ga0.67As. Heliyon 2020, 6, e04852. [Google Scholar] [CrossRef]
  15. Yang, M.M.; Luo, Z.D.; Mi, Z.; Zhao, J.; Sharel P., E.; Alexe, M. Piezoelectric and pyroelectric effects induced by interface polar symmetry. Nature 2020, 377, 584. [Google Scholar] [CrossRef]
  16. Xu, C.; Wei, P.; Wei, Z.; Guo, X. Shear horizontal wave in a piezoelectric semiconductor substrate covered with a metal layer with consideration of Schottky junction effects. Appl. Math. Model. 2022, 109, 509–518. [Google Scholar] [CrossRef]
  17. Pang, Y.; Feng, W.; Liu, J.; Zhang, C. SH wave propagation in a piezoelectric/piezomagnetic plate with an imperfect magnetoelectroelastic interface. Waves Random Complex Media 2019, 29, 580–594. [Google Scholar] [CrossRef]
  18. Song, H.X.; Ke, L.L.; Su, J.; Yang, J.; Kitipornchai, S.; Wang, Y.S. Surface effect on the contact problem of a piezoelectric half-plane. Int. J. Solids Struct. 2020, 185–186, 380–393. [Google Scholar] [CrossRef]
  19. Wang, X.; Li, P.; Jin, F. A generalized dynamic model of nanoscale surface acoustic wave sensors and its applications in Love wave propagation and shear-horizontal vibration. Appl. Math. Model. 2019, 75, 101–115. [Google Scholar] [CrossRef]
  20. Zhu, F.; Pan, E.; Qian, Z.; Wang, Y. Dispersion curves, mode shapes, stresses and energies of SH and Lamb waves in layered elastic nanoplates with surface/interface effect. Int. J. Eng. Sci. 2019, 142, 170–184. [Google Scholar] [CrossRef]
  21. Enzevaee, C.; Shodja, H.M. Torsional surface wave propagation in a transversely isotropic FG substrate with piezoelectric over-layer within surface/interface theory. Acta Mech. 2020, 231, 2203–2216. [Google Scholar] [CrossRef]
  22. Guo, X.; Wei, P.; Xu, M.; Lan, M. Dispersion relations of anti-plane elastic waves in micro-scale one dimensional piezoelectric semiconductor phononic crystals with the consideration of interface effect. Mech. Mater. 2021, 161, 104000. [Google Scholar] [CrossRef]
  23. Zhang, Q.; Li, M.; Zhao, M. Dynamic analysis of a piezoelectric semiconductor nanoplate with surface effect. Mater. Today Commun. 2022, 33, 104406. [Google Scholar] [CrossRef]
  24. Li, D.; Li, S.; Zhang, C.; Chen, W. Propagation characteristics of the shear horizontal wave in piezoelectric semiconductor nanoplates incorporating surface effect. Int. J. Mech. Sci. 2023, 247, 108201. [Google Scholar] [CrossRef]
  25. Karpfinger, F.; Gurevich, B.; Bakulin, A. Modeling of wave dispersion along cylindrical structures using the spectral method. J. Acoust. Soc. Am. 2008, 124, 859. [Google Scholar] [CrossRef] [PubMed]
  26. Li, K.; Jing, S.; Yu, J.; Zhang, X.; Zhang, B. The Complex Rayleigh Waves in a Functionally Graded Piezoelectric Half-Space: An Improvement of the Laguerre Polynomial Approach. Materials 2020, 13, 2320. [Google Scholar] [CrossRef]
  27. Li, K.; Jing, S.; Yu, J.; Zhang, B. Complex Rayleigh Waves in Nonhomogeneous Magneto-Electro-Elastic Half-Spaces. Materials 2021, 14, 1011. [Google Scholar] [CrossRef]
  28. Pang, Y.; Xu, C.; Ge, T.; Feng, W. Shear horizontal wave propagation along a periodic metal grating surface of a magneto-electro-elastic substrate. J. Appl. Phys. 2019, 125, 165104. [Google Scholar] [CrossRef]
  29. Patel, R.; Adhikari, M.S.; Boolchandani, D. Active area optimization of film bulk acoustic resonator for improving performance parameters. Electron. Lett. 2020, 56, 1191–1194. [Google Scholar] [CrossRef]
  30. Hu, Y.; Dong, B.; Lei, L.; Wang, Z.; Ruan, S. High-order harmonic film bulk acoustic resonator based on a polymer reflector. Sensors 2022, 22, 7439. [Google Scholar] [CrossRef]
Figure 1. Propagation of the shear horizontal wave in a PSC semi-infinite medium.
Figure 1. Propagation of the shear horizontal wave in a PSC semi-infinite medium.
Mathematics 12 00560 g001
Figure 2. The dispersion and attenuation curves of shear horizontal waves calculated by the second mathematical model without the consideration of biasing electric fields. (a) The dispersion curves H = 1   μ m ; (b) The attenuation curves H = 1   μ m .
Figure 2. The dispersion and attenuation curves of shear horizontal waves calculated by the second mathematical model without the consideration of biasing electric fields. (a) The dispersion curves H = 1   μ m ; (b) The attenuation curves H = 1   μ m .
Mathematics 12 00560 g002
Figure 3. The comparisons of attenuation curves of shear horizontal waves calculated by the first and second mathematical models without the consideration of biasing electric fields. (a) The attenuation curves H = 1   μ m ,   N D L = 1 × 10 21   m 3 ; (b) The attenuation curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ; (c) The attenuation curves H = 5   μ m ,   N D L = 1 × 10 21   m 3 ; (d) The attenuation curves H = 5   μ m ,   N D L = 10 × 10 21     m 3 ; (e) The attenuation curves H = 10   μ m ,   N D L = 1 × 10 21   m 3 ; (f) The attenuation curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 .
Figure 3. The comparisons of attenuation curves of shear horizontal waves calculated by the first and second mathematical models without the consideration of biasing electric fields. (a) The attenuation curves H = 1   μ m ,   N D L = 1 × 10 21   m 3 ; (b) The attenuation curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ; (c) The attenuation curves H = 5   μ m ,   N D L = 1 × 10 21   m 3 ; (d) The attenuation curves H = 5   μ m ,   N D L = 10 × 10 21     m 3 ; (e) The attenuation curves H = 10   μ m ,   N D L = 1 × 10 21   m 3 ; (f) The attenuation curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 .
Mathematics 12 00560 g003
Figure 4. The dispersion and attenuation curves of shear horizontal waves with the variation of the dimensionless biasing electric field γ x . (a) The dispersion curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (b) The attenuation curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (c) The dispersion curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (d) The attenuation curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (e) The dispersion curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (f) The attenuation curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 .
Figure 4. The dispersion and attenuation curves of shear horizontal waves with the variation of the dimensionless biasing electric field γ x . (a) The dispersion curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (b) The attenuation curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (c) The dispersion curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (d) The attenuation curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (e) The dispersion curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (f) The attenuation curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 .
Mathematics 12 00560 g004aMathematics 12 00560 g004b
Figure 5. The dispersion and attenuation curves of shear horizontal waves with the variation of the dimensionless biasing electric field γ y . (a) The dispersion curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (b) The attenuation curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 , γ x = 0 ; (c) The dispersion curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (d) The attenuation curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (e) The dispersion curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (f) The attenuation curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 .
Figure 5. The dispersion and attenuation curves of shear horizontal waves with the variation of the dimensionless biasing electric field γ y . (a) The dispersion curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (b) The attenuation curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 , γ x = 0 ; (c) The dispersion curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (d) The attenuation curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (e) The dispersion curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (f) The attenuation curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 .
Mathematics 12 00560 g005aMathematics 12 00560 g005b
Figure 6. The comparisons of attenuation curves of shear horizontal waves calculated by the first and second mathematical models with the consideration of biasing electric fields. (a) The attenuation curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (b) The attenuation curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (c) The attenuation curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (d) The attenuation curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (e) The attenuation curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (f) The attenuation curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 .
Figure 6. The comparisons of attenuation curves of shear horizontal waves calculated by the first and second mathematical models with the consideration of biasing electric fields. (a) The attenuation curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (b) The attenuation curves H = 1   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (c) The attenuation curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (d) The attenuation curves H = 5   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 ; (e) The attenuation curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ y = 0 ; (f) The attenuation curves H = 10   μ m ,   N D L = 10 × 10 21   m 3 ,   γ x = 0 .
Mathematics 12 00560 g006aMathematics 12 00560 g006b
Table 1. The physical constants of Au and n-type ZnO [10].
Table 1. The physical constants of Au and n-type ZnO [10].
c y z y z ρ e y y z ε y y μ y y d y y W
ZnO 43 5700 0.48 7.61 0.01 2.6 4.45
Au 42 19,300 0 0 0 0 5.1
Note: cyzyz: GPa, ρ: kgm−3, eyyz: cm−2, εyy: 10−11 CV−1m−1, μyy: m2V−1s−1, dyy: 10−4m2s−1, W: eV.
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

Guo, X.; Wang, Y.; Xu, C.; Wei, Z.; Ding, C. Influence of the Schottky Junction on the Propagation Characteristics of Shear Horizontal Waves in a Piezoelectric Semiconductor Semi-Infinite Medium. Mathematics 2024, 12, 560. https://doi.org/10.3390/math12040560

AMA Style

Guo X, Wang Y, Xu C, Wei Z, Ding C. Influence of the Schottky Junction on the Propagation Characteristics of Shear Horizontal Waves in a Piezoelectric Semiconductor Semi-Infinite Medium. Mathematics. 2024; 12(4):560. https://doi.org/10.3390/math12040560

Chicago/Turabian Style

Guo, Xiao, Yilin Wang, Chunyu Xu, Zibo Wei, and Chenxi Ding. 2024. "Influence of the Schottky Junction on the Propagation Characteristics of Shear Horizontal Waves in a Piezoelectric Semiconductor Semi-Infinite Medium" Mathematics 12, no. 4: 560. https://doi.org/10.3390/math12040560

APA Style

Guo, X., Wang, Y., Xu, C., Wei, Z., & Ding, C. (2024). Influence of the Schottky Junction on the Propagation Characteristics of Shear Horizontal Waves in a Piezoelectric Semiconductor Semi-Infinite Medium. Mathematics, 12(4), 560. https://doi.org/10.3390/math12040560

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