The propagation of acoustic waves in a homogeneous isotropic semiconducting layer sandwiched between two homogeneous transversely isotropic piezoelectric halfspaces has been investigated. The mathematical model of the problem is depicted by a set of partial differential equations of motion, Gauss equation in piezoelectric material and electron diffusion equation in semiconductor along with the boundary conditions to be satisfied at the piezoelectric-semiconductor interfaces. The secular equations describing the symmetric and asymmetric modes of wave propagation have been derived in compact form after obtaining the analytical expressions for various field quantities that govern the wave motion. The complex secular equation has been solved numerically using functional interaction method along with irreducible cardano method. The computer simulated results are obtained with the help of MATLAB software for 6 mm cadmium selenide (<i>CdSe</i>) piezoelectric material and n-type silicon (<i>Si</i>) semiconductor in respect of dispersion curve, attenuation and specific loss factor of energy dissipation for symmetric (sym) and asymmetric (asym) modes of wave propagation. The study may find applications in non-destructive testing, resonators, waveguides etc.
Lord Rayleigh [
Kagan [
Wu and Zhu [
Sorokin [
Sandwich structures found applications in many fields such as spacecraft, aircraft, automobiles, boats and ships. Keeping in view the wide range applications and above referred work an attempt is made to study the acoustic waves in a sandwich structure consisting of semiconducting layer between two piezoelectric halfspaces. The phase velocity, attenuation and specific loss factor of energy dissipation for acoustic and optical modes of wave propagation in case of symmetric and asymmetric modes of wave propagation, in a sandwiched structure have been computed numerically from analytically developed model. The effect of variation of semiconducting layer thickness on the considered wave fields has also been taken in account to show the wave guide nature of the thin layer composite structure.
We consider an n-type homogeneous isotropic semiconductor layer of thickness sandwiched between two homogeneous transversely isotropic, piezoelectric halfspaces. We take the origin of Cartesian coordinate system at any point on the mid-plane of the semiconductor layer and z-axis along the thickness of the layer as shown in
We choose x-axis along the direction of wave propa-
gation in such a way that all particles on the line parallel to y-axis are equally displaced. Thus all field quantities are independent of y-coordinate.
The basic governing equations of motion and electron diffusion for the composite structure under study, in the absence of body forces and electric sources, are given below:
1) Semiconductor (n-type) elastic layer (Sharma et al. [
2) Piezoelectric medium (Sharma et al. [
where the notations, , ,
have been used. In the above equations the superposed dots on various quantities denote time differentiation and comma notation is used for spatial derivatives. Throughout this paper the superscripts and refer to piezoelectric and semiconductor materials, respectively. Here are Lamè’s parameters; is the density; is the elastodiffusive constants of electrons; is the diffusion coefficient of electron; and respectively are the life and relaxation times of the carrier field; and are the equilibrium and non-equilibrium values of electrons concentration; is the coefficient of linear thermal expansion of the semiconductor material. The quantities are flux-like constants and is the uniform temperature; and are displacement vectors for semiconductor and piezoelectric materials, respectively. The quantities, , andare the electric potential, density, elastic parameters and piezoelectric constants; and are the electric permittivity perpendicular and along the axis of symmetry of piezoelectric material, respectively.
The non-vanishing components of stress tensor, current density and electric displacement in both the media are:
Semiconductor:
Piezoelectric:
where and are the stress tensors. The quantities and respectively denote the current density and carrier density gradient in semiconducting layer; is the electric displacement vector of piezoelectric material and is the electronic charge.
Boundary ConditionsThe requirement of continuity of stresses, displacements, electric fields and current density at the interfaces of two media leads to the following interfacial boundary conditions
In order to simplify the model, we define following quantities
where is the characteristic frequency, and, are respectively, the longitudinal and shear wave velocities.
Introducing the quantities (9) in Equations (1) to (7), we obtain (on suppressing the primes for convenience)
In order to facilitate the solution in semiconductor layer, we introduce the scalar and vector point potential functions and through the relations
The equations (10)-(11) with the help of relations (17) provide us
The equation (19) corresponds to purely transverse waves in the semiconductor which get decoupled from rest of the motion and is not affected by the concentration of charge carrier fields.
We take harmonic wave solution of the form
where, here, and are the non-dimensional phase velocity, angular frequency and wave number, respectively, in which primes have been suppressed for convenience.
The use of solution (21) in equations (18)-(20) and (12)-(14) after straightforward algebraic reductions and simplifications leads to the formal solution satisfying the radiation conditions in semiconductor layer and piezoelectric as under:
for semiconductor layer
for piezoelectric halfspace we have
for piezoelectric halfspace we have where
Consequently, the expressions for stresses, current density and displacements in the semiconductor layer are obtained from equations (15) and (17) via (22) and (23) as:
where, , and are the unknowns to be determine.
Using the solution (24) and (25) in Equation (16) we obtain the stresses and electric displacement as:
for piezoelectric halfspace, and
for piezoelectric halfspace, where
.
Upon imposing the boundary conditions (8) we obtain a coupled system of twelve homogeneous simultaneous algebraic equations in twelve unknowns and which has a non-trivial solution if the determinant of the coefficients vanishes. This requirement leads to a dispersion relation of the type
where A, B, C and D are the matrices whose non-zero elements are defined as below:
The non-zero elements of matrix A are
The non-zero elements of matrix B are
, , ,
, , ,
, ,
, ,
, ,
, ,
, ,
The non-zero elements of matrix C can be obtained from corresponding non-zero elements of of matrix B by replacing hyperbolic tangent function with the hyperbolic cotangent function and non-zero elements of matrix D are same as that of corresponding non-zero elements of matrix A. After some algebraic reductions and simplifications the secular equation (37) can be rewritten as:
where
Here +1 stands for asymmetric mode of wave propagation and −1 for symmetric mode of wave propagation. The quantities are defined as
In general, wave number and hence the phase velocities of the waves are complex quantities, therefore the waves are attenuated in space. In order to solve the secular equation (38), we take
where, , and are real numbers.
Here, it may be noted that and respectively represent the phase velocity and attenuation coefficient of the waves. Using representation (41) in various relevant relations, the complex roots can be computed from Equations (27) with the help of Cardano method. The roots are further used to solve secular equation (38) to obtain phase velocity and attenuation coefficient of the surface waves by using function iteration numerical technique outlined below:
In general the secular equation (38) is of the form which on using representation (41) leads to a system of two real equations and . In order to apply functional iteration method we write and, where the functions and are selected in such a way that they satisfy the conditions
for all in the neighborhood of the root. If be an initial approximation to root, then we can construct the successive approximations according to the formulae
The sequence of approximations to the root will converge to the actual root provided lies in the neighborhood of the actual root. For initial value of, the roots are computed from equations (27) by using Cardano method for each value of the non-dimensional wave number for assigned frequency. The values of so obtained are then used in secular Equation (38) to obtain the current values of and each time which are further used to generate the sequence (43). The process is terminated as and when the condition, being arbitrarily small number to be selected at random to achieve the accuracy level, is satisfied. The procedure is continuously repeated for different values of non-dimensional wave number to obtain corresponding values of the phase velocity and attenuation coefficient. Thus, the real phase velocity and attenuation coefficient during the propagation of Rayleigh type waves in the composite structure under study can be computed from dispersion relation (38).
The specific loss is the direct method of defining the internal friction for a material. According to Kolsky [
In this section we present the discussion of the results obtained through computer simulations from the analytical developments in the previous sections for a sandwich structure consisting of n-type silicon layer and cadmium selenide halfspaces, whose material parameters and constants are defined as [
1) Piezoelectric halfspace:
, , , , , , , ,
,.
2) Semiconductor layer:
, , , , ,.
The numerical computations have been performed with the help of MATLAB software employing the procedure outlined in section 5 for the first two symmetric (sym) and asymmetric (asym) modes of wave propagation. Here the quantity denotes the non-dimensional wave number of acoustic waves traveling in the semiconducting layer sandwiched between piezoelectric halfspaces.
From
sat higher frequencies (short wavelengths), the semiconductor material layer behave like a thick slab and hence coupling between upper and lower interfaces is reduced as a result the properties of symmetric and asymmetric modes become almost similar. Under such situations the wave motion at the upper interface is not confined to the lower interface and displacement becomes localized near the respective boundaries and as a result the wave dispersion curves asymptotically approach to that of Rayleigh surface wave in case of acoustic modes and to shear wave velocity for optical modes in semiconductor material.
Moreover, the unbounded displacement field is characterized by the singularities of circular tangent functions. It is also observed that the coupling effect of various interacting fields decreases with increasing thickness of the semiconductor layer resulting in lower phase velocity as the wave travel along the interface (guided surface) with little disturbance to the semiconductor layer. The phase velocity of various modes at long wavelengths is noticed to be quite high because such waves penetrate deep in to the medium thereby creating significant disturbance in either of the material components because various interacting field become operative and contribute in increasing the magnitude of the phase velocity. Moreover, at long wavelengths the semiconductor layer behaves like a thin structure and hence its properties are greatly influenced by the substrate piezoelectric material under such situation through the interfaces which also attribute in high magnitude of phase velocity.
ation of both acoustic and optical modes decreases before it attains a steady state for. It is attributed to the fact that large wavelength waves interact with the medium in contact to the maximum extent in contrast to that at short wavelengths in which case these follow the surface without significantly disturbing the medium. Thus the medium offer a high resistance to wave propagation at long wavelengths as compared to that at short wavelengths. At large wavelengths optical modes of wave propagation possess small value of attenuation as compared to that of acoustic modes; however the trend gets reversed at small wavelengths. The optical modes both symmetric and asymmetric possess higher attenuation.
The increase of phase velocity with decreasing thickness of the semiconductor layer shows that the thin sandwiched layer acts like an acoustic waveguide. Acoustic waveguide confine the acoustic energy in the bulk of the
layer through multiple reflections from the top and bottom interface as wave propagates and it generates a particle displacement at both the interfaces to ensure the channelized flow of energy through the waveguide.
1) Functional iteration numerical technique along with Cardano method has been successfully used to solve the complex secular equation in order to compute the phase velocities, attenuation coefficients and specific loss factors of energy dissipation for different symmetric and asymmetric modes of wave propagation in a sandwich piezoelectric-semiconductor structure.
2) Optical modes possess higher phase velocities than that of acoustic modes in symmetric as well as asymmetric mode of wave propagation.
3) Phase velocity profiles show dispersive character and decrease with the increasing wave number to become smooth and stable at short wavelengths.
4) The attenuation and specific loss factor increase with increasing wave number initially and then decrease after attaining maximum value in each case.
5) At large wavelengths acoustic modes possess high attenuation and specific loss factor of energy dissipation and reverse trends are noticed at small wavelengths.
6) Phase velocity increases with decreasing thickness of semiconducting layer while attenuation and specific loss factor of energy dissipation increase with the thickness of the layer showing a waveguide nature of the thin sandwiched semiconductor layer.
7) The study may find applications in fabrication of micro-electromechanical surface acoustic wave devices.
The authors are thankful to the reviewer for his useful suggestions for the improvement of this work. The author (JNS) thankfully acknowledges the financial support provided by CSIR, New Delhi via project grant No: 025(0184) EMR-II.