_{1}

^{*}

This theoretical paper analytically predicts the existence of new surface wave in propagation direction [101] in the cubic piezoelectromagnetics. The solution for the velocity of the new wave is given in an explicit form. Such new wave possesses one real component and two purely imaginary components. This corresponds to a leaky acoustic SH-wave. However, in this case the real component does not participate in the complete displacements. As a result, the new wave can represent the new shear-horizontal surface acoustic wave (SH-SAW) for suitable boundary conditions. For the me-chanically free surface, several combinations of the following electrical and magnetic boundary conditions were used: electrically closed, electrically open, magnetically closed, magnetically open surface. This new SH-SAW can propagate with the speed slightly larger than that for the SH bulk acoustic wave coupled with both the electrical and magnetic potentials. The existence conditions for the new SH-SAW were also discussed. They can be very complicated and depend on all the material parameters. It was also discussed that the results can be true for the left-handed metamaterials. It is thought that the new SH-SAW can be produced by electromagnetic acoustic transducers (EMATs) because generation of SH-SAWs is feasible with the EMATs. This can be a problem for experimentalists working in the research arenas such as acoustooptics, photonics, and acoustooptoelectronics.

There is a continuous interest in the various theoretical and experimental investigations of the magnetoelectric (ME) effect in (composite) materials for development of smart materials in the microwave technology. To start a short historical review of works on the magnetoelectric effect, it is possible first of all to mention several classical works [1-5] published in 1970s which have originally studied composite materials and the ME effect. Indeed, piezoelectrics and piezomagnetics can be bonded together in different ways to form composites that exhibit the ME effect. The ME effect can lead to the following: An electrical signal caused by the presence of the piezoelectric phase of such composites can be obtained as a result of the application of a magnetic field affecting the composite piezomagnetic phase. On the other hand, the composite piezomagnetic phase can magnetically respond because of the application of an electrical field to the ME composite possessing the piezoelectric phase.

It is thought that it is unnecessary to review all the theoretical and experimental works concerning the ME effect and creation and characterization of new two-phase ME composites. One can readily find thousands of papers on the subject. Therefore, it is suitable to give only several review papers on the subject. Indeed, the reader can find that the recent reviews cited in Refs. [6-11] contain enough citations. For the reader who would like to know more about the subject, it is strongly recommended to use the following several additional review papers cited in Refs. [12-15]. According to review paper [

Concerning the theoretical investigations of the propagation problems of shear-horizontal surface acoustic waves (SH-SAWs), it is essential to mention the excellent work [

It is thought that there is a lack of works on existence of leaky acoustic waves propagating in piezoelectromagnetics. The leaky acoustic waves, also called the leaky surface acoustic waves (LSAWs), are also very important in acoustics of solids. Generally, some LSAWs are numerically studied in appropriate piezoelectrics and they are longitudinal, for example, see Refs. [26-31]. The works on the LSAW studies frequently use promising noncubic and nontransversely-isotropic piezoelectric monocrystals such as Quartz, Langasite, LiTaO_{3}, etc. The longitudinal LSAWs propagating on suitable cuts of piezoelectric monocrystals are at least three-partial. It is well-known that the SH-SAWs such as the surface BGwave are two partial and can propagate only on suitable cuts of piezoelectrics or piezomagnetics. It was also mentioned in Ref. [

However, an existence of any SH-LSAW in addition to the SH-SAW was not demonstrated, probably, due to the fact that they are two-partial in pure piezoelectrics. As soon as the wave propagation in piezoelectromagnetics is treated, the SH-waves can become three-partial. So, it is expected that some LSAW propagation can be found in piezoelectromagnetics. The main purpose of this report is to analytically demonstrate a possible existence of new SH-wave propagating in the cubic piezoelectromagnetics on the same cuts in addition to the seven new SH-SAWs discovered in the recent book cited in Ref. [

Following the theoretical description of SH-SAW propagation in the cubic piezoelectromagnetics described in book [_{1}- and x_{2}-axes lie in the interface plane and the x_{3}-axis is directed perpendicular to both the axes. The SH-wave propagate along the x_{1}-axis and are polarized along the x_{2}-axis. In the case of the SH-SAW, the waves must damp towards the depth of the composite material, namely along the negative values of the x_{3}-axis. However, this is not obligatory in the case of the SH-LSAW because at least one partial component of three must have non-SAW behavior. Indeed, the plane wave approximation is used in the theory of SAW and LSAW propagation.

Using only nonzero material constants, it is possible to write the coupled equations of motion for the studied case. They read

In the first equation written above, t stands for time and U_{2} is the corresponding component of the mechanical displacement directed along the x_{2}-axis. Also, φ and ψ are the electrical and magnetic potentials, respectively, defined by E = –¶φ/¶x_{3} and H = –¶ψ/¶x_{3} where E and H are the electrical and magnetic fields. In the equations written above, the non-zero material constants are as follows: C_{44} = C_{66} = C (elastic stiffness constants), e_{16} = –e_{34} = e (piezoelectric constants), h_{16} = –h_{34} = h (piezomagnetic coefficients), ε_{11} = ε_{33} = ε (dielectric permittivity coefficients), μ_{11} = μ_{33} = μ (magnetic permeability coefficients), α_{11} = α_{33} = α (electromagnetic constants), and r is the mass density.

The plane wave solution for the mechanical displacement U_{2}, electrical potential φ, and magnetic potential ψ can be written in the following form:

In equation (4), the displacements are defined by U_{2} = U, U_{4} = φ, U_{5} = ψ, and are the initial amplitudes. Also, the directional cosines are n_{1} = 1, n_{2} = 0, and n_{3} ≡ n_{3}. V_{ph} and j stand for the phase velocity and the imaginary unity, respectively. The wavenumber k in Equation (4) is coupled with the components {k_{1}, k_{2}, k_{3}} of the wavevector k as follows: {k_{1}, k_{2}, k_{3}} = k {n_{1}, n_{2}, n_{3}}.

Exploiting the solutions in Equation (4) for equations from (1) to (3), the coupled equations of motion can be written in a tensor form. Using the corresponding components of the symmetric GL-tensor in the well-known modified Green-Christoffel equation [

(5)

where U^{0} = U_{2}^{0}, φ^{0} = U_{4}^{0}, and ψ^{0} = U_{5}^{0} represent three components of the eigenvector. Note that the readers themselves can obtain the explicit forms for the GL-tensor components in Equation (5) using plane wave solutions (4) for equations from (1) to (3).

It is obvious that the three-component eigenvector should be nonzero for each eigenvalue n_{3} = k_{3}/k. Suitable eigenvalues n_{3} can be obtained when the following matrix determinant of the system of homogeneous Equations (5) becomes equal to zero:

It is useful to rewrite the characteristic determinant in Equation (6) in the following form:

(7)

where m = 1 + sqr(n_{3}) and V_{t}_{4} = sqrt(C/r) is the speed of the bulk acoustic wave (BAW) uncoupled with both the electrical and magnetic potentials.

Analyzing the first row or the first column of the matrix determinant in Equation (7), it is natural to treat the following formula for the new wave velocity:

where V_{tem} = V_{t}_{4 }sqrt(1 + K_{em}^{2}) is the speed of the SH-BAW coupled with both the electrical and magnetic potentials because K_{em}^{2} represents the coefficient of the magnetoelectromechanical coupling (CMEMC). It reads

Therefore, the left-hand side in Equation (7) can be written in the form of two factors. The transformed equation can be written as follows:

The first factor in Equation (10) reveals the following eigenvalue n_{3}:

In the studied case of the wave propagation problem, the first eigenvalue in expression (11) can be real. Using it for Equation (5), one can obtain the corresponding eigenvector components. It is apparent that the simplest and most convenient eigenvector components for the first eigenvalue n_{3}^{(1)} in Equation (11) are

The second factor representing the determinant in Equation (10) can reveal the rest two eigenvalues. Expanding this determinant, one can get the following secular equation:

It is clear that Equation (13) also has two factors on the left-hand side. Therefore, the equality in Equation (13) is true when either the first or second factor equals to zero. In the first case of m = 0, the following eigenvalue can be obtained:

Note that in the wave propagation problem, the sign for the eigenvalue in relation (14) is chosen negative in order to cope with wave motions which must damp towards the depth of the bulk piezoelectromagnetics. This is like the SH-SAW propagation problem. Using the eigenvalue n_{3}^{(2)}, one can check that Equation (5) equals to zero with the following eigenvector components:

The equality to zero of the second factor in Equation (13) is possible as soon as the eigenvalue n_{3} becomes the following function of the CMEMC:

It is crucial to state that it is clearly seen in Equation (16) that for K_{em}^{2} = 1 one gets n_{3}^{(3)} = 0. This solidly demonstrates that K_{em}^{2} = 1 results in V_{new} = V_{tem}, see expression (8). Therefore, this case is unique and demonstrates some association with the SH-BAW solution. Also, it is apparent that the value of n_{3}^{(3)} becomes real as soon as K_{em}^{2} > 1 in expression (16). Therefore, it is necessary to cope with K_{em}^{2} < 1 for simplicity. The case of K_{em}^{2} = 0 gives n_{3}^{(3)} = n_{3}^{(2)} = –j. This is the case of two equal eigenvectors. This situation always leads to zero value of the boundary-condition determinant and results in zero values of complete displacements. This is like the situation occurred in the cubic piezoelectrics [_{tem} [_{em}^{2} = 0 is a unique case, but it is not possible to solidly say that this situation corresponds to an SH-SAW solution.

Utilizing the third eigenvalue defined by relation (16) for Equation (5), one can find that two possible sets of the eigenvector components can exist in this case. The first set of the eigenvector components is expressed as follows:

Also, the second set of the eigenvector components can be written as follows:

The first problem of the finding of the new SH-wave velocity, eigenvalues, and corresponding eigenvectors was resolved above. The speed of the new SH-wave is given by expression (8). Three eigenvalues n_{3} were found in the explicit forms. For the new SH-wave propagation in the cubic piezoelectromagnetics, the first eigenvalue is real and defined by expression (11), but the second and third eigenvalues are purely imaginary and defined by expressions (14) and (16), respectively. Indeed, all the eigenvalues possess the corresponding eigenvectors and the third eigenvalue has two different sets of the eigenvector components.

The second problem is the determination of the existence conditions of the new SH-wave propagation. For this purpose, it is necessary to consider the mechanical, electrical, and magnetic boundary conditions. Al’shits, Darinskii, and Lothe [_{3} = 0 must vanish, following books [17,20]. The electrical and magnetic boundary conditions such as the electrically closed surface (φ = 0) and the magnetically open surface (ψ = 0) can be also applied. These boundary conditions are those which reveal the surface BGM-wave speed in the transversely isotropic piezoelectromagnetic materials [16,17] and the cubic piezoelectromagnetics [_{32}(x_{3} = 0) = 0 at the interface x_{3} = 0 [

Following the theoretical treatments done in book [

In Equations (19), F_{1}, F_{2}, and F_{3} are the weight factors which must be determined. They are very important parameters. The complete mechanical displacement U^{Σ}, complete electrical potential φ^{Σ}, and complete magnetic potential ψ^{Σ} depend on them and can be written in the plane wave forms as follows:

where x_{3} < 0 and V_{new} is defined by expression (8).

To determine the values of the weight factors, it is necessary to use in Equation (19) the eigenvalues and the corresponding eigenvector components. Exploiting them and the first set of the eigenvector components defined by expression (17), one can obtain the explicit form for the first third-order boundary-condition determinant (BCD3) of the coefficient matrix in Equation (19). It is as simple as follows:

Expanding it, one can obtain the following secular equation, in which only material constants are involved because the phase velocity V_{ph} was determined in Equation (8):

It is clearly seen in Equation (24) that the left-hand side is formed by two factors. Therefore, there are two possibilities to satisfy the equality in Equation (24) when either the first or second factor equals to zero. The equality to zero of the first factor leads to the following existence condition for the new wave propagation in the cubic piezoelectromagnetics:

For this condition, it is possible to find the explicit forms of the weight factors. Using Equations (19) and (23), one can check that the simplest forms are expressed as follows:

In Equation (24), the equality to zero of the second factor representing the CMEMC K_{em}^{2} defined by relation (9) is inappropriate here. This is true because n_{3}^{(3)}(K_{em}^{2} = 0) = n_{3}^{(2)}, see formulae (14) and (16). In this case there are two equal eigenvalues which give the same sets of the eigenvector components. This situation was also mentioned in the previous section. As a result, F_{2} = –F_{3} together with F_{1} = 0 will zero the complete displacements defined by formulae from (20) to (22). So, this possibility is excluded.

Therefore, the complete displacements defined by formulae from (20) to (22) represent those corresponding to the new wave propagation. These new wave can be called surface electromagnetic wave due to F_{1} = 0 in expression (26) and U^{0(3)} = 0 (see expression (27) below) in expression (17) for condition (25). They then read:

where x_{3} < 0 means that the new wave must damp towards the depth of the cubic piezoelectromagnetics. In expressions from (27) to (29), the speed of the new wave V_{new} is defined by expression (8). Using definition (17) for U^{0(3)} and existence condition (25), it is flagrantly seen in Equation (27) that the complete mechanical displacement U^{Σ} is equal to zero.

This fact can be understood by the way that the wave propagation specifics exhibits a compensative character, namely U^{0(3)} = (hε – eα)K_{em}^{2} = 0 (17) due to hε = eα (25). This can mean that this slow wave is truly acoustic and the piezoelectromagnetic properties can compensate the mechanical ones resulting in zero value of the mechanical displacement during the wave propagation. Indeed, this slow wave propagates with the speed slightly above the SH-BAW speed V_{tem}. It is well-known that acoustic wave speeds including the V_{tem} and V_{new} are approximately five orders slower than the speed of the electromagnetic wave propagating in a bulk solid. It is defined by the following relation: V^{2} = (εµ)^{–1}. This extreme slowness of the new wave illustrates that some connection with the mechanical displacement is conserved. Therefore, the new wave can be called the surface acoustic magnetoelectric wave to illuminate that this wave relates to an acoustic branch, but not an optic branch.

This new electromagnetic wave can be also called the surface acoustic-phonon polariton (SAPP) to distinguish from the surface optic-phonon polaritons (SOPPs). The SOPPs [

Utilizing the second set of the eigenvector components defined by expression (18) for the third eigenvalue, the second BCD3 reads:

Expanding this BCD3, the following second existence condition of the wave propagation in the cubic piezoelectromagnetics can be revealed:

Therefore, the explicit forms of the weight factors are as follows:

For this second case, the complete displacements of the new wave can be written in the same forms defined by formulae from (27) to (29), where the weight factor F_{3} from definitions (32) and the second set of the eigenvector components given in expression (18) must be used. It is thought that it is vital to treat the other possible sets of the boundary conditions for comparison.

Using this set of the mechanical, electrical, and magnetic boundary conditions, it is possible to write the following homogeneous equations [

The matrix BCD3 representing a number can be then written in the following simplified form, using Equation (17):

The BCD3 in Equation (34) can be readily reduced to the following determinant:

It is clearly seen in the first row of determinant (35) that the determinant equals to zero when (25). This coincides with the result obtained in Subsection 2.1. Using the second row of the determinant, the explicit forms of the weight factors are written as follows:

The weight factors are used in equations from (27) to (29), too. It is also necessary to mention the second solution for Equation (35). This is as follows:

This second solution (37) looks like unreal or it is hard to realize it. However, it is thought that it is useful to record any solution to have a more complete picture of the problem of wave propagation.

Using eigenvector components (18) for Equations (33), one can get the following BCD3:

It is obvious that one can transform determinant (38) by the same way, see the transformations of determinant (35) written above. It was found that determinant (38) equals to zero when (31). This also coincides with the result obtained in Subsection 2.1. There is also the second solution which coincides with expression (37). For the case of the first solution for Equation (38), the explicit forms of the weight factors read:

The values of weight factors (39) are also utilized in equations from (27) to (29). Therefore, it is possible to state that the cases of the boundary conditions used in this and the previous subsections demonstrate the same results.

This is the case of the mechanically free (σ_{32} = 0), electrically closed (φ = 0), and magnetically closed (B = 0) surface. It is obvious that this is the case when the electrical and magnetic boundary conditions are used from Subsections 2.1 and 2.2, respectively. Therefore, the reader can readily form the corresponding three homogeneous equations in the matrix form (this form can be also found in book [

where the nondimensional value of K_{α}^{2} was introduced in the recent book cited in Ref. [_{α}^{2} = eh/(Cα).

Equation (40) represents the existence condition for the new wave propagation in the case studied in this subsection. The new wave propagates with the speed defined by relation (8). It is evident that Equation (40) is very complicated. For small values of α, the value of K_{α}^{2} can be very large. Therefore, it is necessary to cope with small values of the material constant h in order to have appropriate values of K_{α}^{2}, namely 0 < K_{α}^{2} < 1 similar to 0 < K_{em}^{2} < 1. It is clear that for small values of α and h, the value of K_{em}^{2} is completely defined by the piezoelectric properties. This can be the case of the dominant piezoelectric phase. However, it is thought that this dominance is not obligatory in the case of a large value of α. Indeed, the value of K_{em}^{2} can be large, but less than unity. Consequently, equality (40) can be fulfilled when the following occurs: K_{α}^{2} → K_{em}^{2} → K_{e}^{2}, where K_{e}^{2} is called the coefficient of the electromechanical coupling (CEMC). It represents the well-known non-dimensional characteristic for a pure piezoelectrics and can also define the piezoelectric phase of a piezoelectromagnetics. It is defined by the following relation:

Also, it is possible to introduce the weight factors for this case in non-dimensional forms. They are written as follows:

Therefore, it is possible to state that existence condition (40) can give non-zero values of the complete mechanical displacement U^{Σ} in Equation (27). This means that for the mechanically free, electrically closed, and magnetically closed surface, the new surface wave can represent the new SH-SAW coupled with both the electrical and magnetic potentials.

Using Equation (18), it is possible to form the second BCD3 for determination of the existence conditions of the new wave propagation for this set of the boundary conditions. It was found that three solutions (existence conditions) can exist in this case. The first existence condition coincides with that revealed in formula (31):. The second solution such as K_{em}^{2} = 0 is inappropriate because it gives two equal eigenvalues n_{3}^{(2)} = n_{3}^{(3)} = –j. This was discussed above. The third existence condition is complicated. This is a limit case and gives K_{m}^{2} → K_{em}^{2} → 1. However, K_{em}^{2} = 1 is also unsuitable because it gives V_{new} = V_{tem}. This was also discussed above. The limit condition reads:

where K_{m}^{2} is called the coefficient of the magnetomechanical coupling (CMMC). It represents the well-known non-dimensional characteristic for a pure piezomagnetics and can also define the piezomagnetic phase of a piezoelectromagnetics. It is defined by

The weight factors can be then written in nondimensional forms as follows:

In this case of the mechanically free (σ_{32} = 0), electrically open (D = 0), and magnetically open (ψ = 0) surface, two different sets of eigenvector components (17) and (18) must be also applied. It is obvious that this is the case when the electrical and magnetic boundary conditions can be borrowed from Subsections 2.2 and 2.1, respectively. Indeed, the reader can also form three homogeneous equations in the matrix form and expand the BCD3 of the coefficient matrix by the way demonstrated in Subsections 2.1 and 2.2. Using definition (17) to form the first BCD3 and expanding the BCD3, one can also obtain three limit existence conditions. All of them are unsuitable similar to those for the second BCD3 from Subsection 2.3. They read: (25), K_{em}^{2} = 0, and K_{e}^{2} → K_{em}^{2} → 1. Indeed, K_{em}^{2} = 1 is inapt because it gives the case of V_{new} = V_{tem} discussed above and results from the following limit condition:

Therefore, the weight factors are as follows:

Using definition (18) to form the second BCD3, one can find single suitable existence condition because the BCD3 can reduce to the following equality after expansion:

Existence condition (48) looks like existence condition (40) for the new wave propagation in the case of Subsection 2.3. In the case of this subsection, the new wave also propagates with the speed defined by relation (8). Indeed, for small values of α, the value of K_{α}^{2} can be very large. Therefore, it is necessary to cope with small values of the material constant e in order to get the following: 0 < K_{α}^{2} < 1 similar to 0 < K_{em}^{2} < 1. It is lucid that for very small values of α and e, the value of K_{em}^{2} is completely defined by the piezomagnetic properties: The piezomagnetic phase can be dominant and this dominance is not obligatory in the case of a large value of α. Indeed, the value of K_{em}^{2} can be large, but less than unity. Consequently, Equality (48) can be also fulfilled when the following occurs: K_{α}^{2} → K_{em}^{2} → K_{m}^{2}. It is expected that in the case of large values of α and e, Equality (48) with K_{α}^{2} < 1 can be also true for suitable values of all the material constants of the cubic piezoelectromagnetics. Finally, it is indispensable to mention that the weight factors for this case are those introduced in Equation (42) in the nondimensional forms.

These theoretical investigations carried out above soundly demonstrated that an additional surface wave can propagate in direction [

Concerning the propagation of the new wave in the cubic piezoelectromagnetics, it is possible to briefly discuss the revealed existence conditions given by formulae (25) and (31) in the previous section. These obtained existence conditions are true for the applied boundary conditions of Subsection 2.1 such as the normal component of the stress tensor must vanish, φ = 0 and ψ = 0. Using the other boundary conditions in subsections from 2.2 and 2.4, one can find the other existence conditions which can be different from obtained Equalities (25) and (31). Also, the suitable existence conditions obtained in Subsections 2.3 and 2.4 can result in the fact that the new surface wave can become the new SH-SAW. It is worth noting that the wave speed defined by Formula (8) is the same for each suitable existence condition resulting in the possible propagation of the new surface electromagnetic wave or the new SH-SAW.

Indeed, existence Conditions (25) and (31) completely depend on the material constants of the cubic piezoelectromagnetics. The first condition defined by Equality (25), namely demonstrates that this case can be realized with suitable cubic piezoelectromagnetics which can have as large as possible values of the piezoelectric constant e and the electromagnetic constant α. Their typical values are as follows: from ~0.1 to ~10 C/m^{2} and by about 10^{–12} s/m, respectively. On the other hand, the values of the piezomagnetic coefficient h and dielectric permittivity constant ε must be small enough. The typical values of h are from ~1 to ~10^{3} Tesla and those for ε are ~10^{–10 F/m. It is possible to say that here the piezoelectric phase must be dominant. It is obvious that it is necessary to use suitable cubic piezoelectromagnetics with a very weak piezomagnetic phase for experimental investigations. It is possible that suitable values of h must be from ~10–3} to ~10^{–5} Tesla or even less. It is thought that such small values of h can be also achieved in composites when a piezoelectric-phase matrix is used in which a weak piezomagnetic phase is properly added. Therefore, the experimental techniques for measurements of the material constant h must be improved because experimentalists frequently write zero instead of such small values of h.

Indeed, the values of all the material constants depend on the value of the applied magnetic or electric field. Therefore, one can deal with many parameters and it is possible to suggest that some material constants can depend stronger than the others. This can be different for different piezoelectromagnetics. Therefore, it is possible that a class of suitable piezoelectromagnetics can be formed. This activity can represent intensive experimental investigations in the future for decades. It is thought that some technical devices using such surface wave can be extremely sensitive (supersensors) because any infinitesimal change in the applied magnetic (or electric) field can cause a propagation problem for such surface wave. Therefore, the wave can vanish. As soon as propagation Condition (25) is restored, the wave can then propagate anew.

The second condition defined by Formula (31), namely also couples four material constants such as the piezoelectric e, piezomagnetic h, magnetic permeability µ, and electromagnetic α. Because the value of α is always very small, the value of h must be as large as possible and the values of e and µ must be as small as possible. It is obvious that the value of µ is restricted by the value of µ_{0} for a vacuum. However, it is thought that there is no restriction to minimize the value of the piezoelectric constant e down to suitable very small values. Indeed, this class of cubic piezoelectromagnetics is for those with a dominant piezomagnetic phase. Therefore, it is very important to account measured very small value of the piezoelectric constant e, but not write zero instead. It is possible that such native cubic piezoelectromagnetics can exist. Also, it is thought that artificial composites can be readily created when a strong cubic piezomagnetics can be used as the suitable matrix to solute a very weak piezoelectric phase.

Also, it is possible to mention the other artificial materials which are well-known. They are metamaterials. These materials possess µ < 0 and ε < 0 resulting in εµ > 0. It is apparent that for such metamaterials, the value of the electromagnetic constants α should have a negative sign to satisfy the conditions written in Formulae (25) and (31). This theoretical work does not have a purpose to solidly demonstrate that Conditions (25) and (31) can be also true for some suitable metamaterials. Indeed, this is possible. Note that the metamaterials are new material and they are extensively studied concerning various applications. However, it is possible to review some recent papers concerning a variety of investigations of the metamaterials. Refs. [37,38] have reported their studies of left-handed artificial materials (metamaterials) in the frequency region from 1 THz to 100 THz, and even above [

Also, it is necessary to state that this problem of the surface wave propagation discussed above in this work relates to the cubic piezoelectromagnetics (two-phase materials with the cubic symmetry) but not to the pure piezomagnetics or pure piezoelectrics possessing the cubic symmetry. These single-phase materials were theoretically studied in Refs. [

This report represents the theoretical description of the new wave propagation in direction [_{em}^{2} < 1. The existence conditions for the new wave propagation were revealed. These conditions can be complicated and are coupled with the values of the material constants of the cubic piezoelectromagnetics. It is wellknown that the value of the electromagnetic constant α is very small. Therefore, to satisfy some existence conditions, the two-phase (composite) material must possess the material properties resulting in the domination of the piezoelectric phase with a significantly weaker piezomagnetic phase or the domination of the piezomagnetic phase with a significantly weaker piezoelectric phase. Indeed, if the suitable existence conditions obtained in subsections from 2.1 to 2.4 can be realized, the new wave can propagate. Using suitable Condition (40) or (48), it is expected that the new wave can also represent the new SH-SAW. It is also expected that any infinitesimal perturbation of the medium surface along the new wave propagation way can cause a dramatic attenuation of such wave. This can be used in creation of various technical devices, for instance, supersensors.

The author thanks to the referee and the members of the Editorial Board for a large interest in my theoretical work.