On Propagation Problems of New Surface Wave in Cubic Piezoelectromagnetics

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 mechanically 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.


Introduction
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][2][3][4][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][7][8][9][10][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][13][14][15].According to review paper [9] by Fiebig, such two-phase ME composites are potential candidates for utilization as various acoustic devices, magnetic-field probes, and hydrophones, in electronic packaging and medical ultrasonic imaging, also as sensors and actuators.However, all these review papers did not mention achievements in the theory of the wave propagation in the piezoelectromagnetic composite materials possessing the ME effect.
Concerning the theoretical investigations of the propagation problems of shear-horizontal surface acoustic waves (SH-SAWs), it is essential to mention the excellent work [16] published in 2007 by Arman Melkumyan.In his theoretical work, he discovered twelve new SH-SAWs propagating on the free surface of the transversely isotropic piezoelectromagnetics.Treating the wave propagation in these piezoelectromagnetics, Ref. [17] has also introduced seven new SH-SAWs and also confirmed that the SH-SAW called the surface Bleustein-Gulyaev-Mel-kumyan (BGM) wave can propagate in such piezoelectromagnetics.It is worth noticing that the surface BGMwave in piezoelectromagnetics is analogous to the wellknown surface Bleustein-Gulyaev (BG) wave [18,19] propagating in the transversely isotropic piezoelectrics.Moreover, recent book [20] has analytically found that the BGM-wave can also propagate in the cubic piezoelectromagnetics.This book published in 2011 has also discovered seven new SH-SAWs propagating on the surface of the cubic piezoelectromagnetics and discussed the main differences between the wave propagations in the cubic piezoelectromagnetics and the transversely isotropic piezoelectric composite materials.Note that the wave propagations in cubic piezomagnetics [21] and cubic piezoelectrics [22] are also different from those in the transversely isotropic materials.Also, Ref. [23] has stated that SH-SAWs can easily be produced by electromagnetic acoustic transducers (EMATs).According to books [24,25], the EMATs offer a series of advantages over traditional piezoelectric transducers.So, the EMATs can be used for measurements of SH-SAW characteristics when the wave propagations in the piezoelectrics, piezomagnetics, or piezoelectromagnetics are studied.
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][27][28][29][30][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. [27] that the LSAWs can be used as one of two possible ways to resolve the problem of movement of mobile communication equipments towards higher frequency range.Indeed, the LSAW speed should be higher than the SAW one, using the same crystal cut and propagation direction.
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. [20].It is also noted that the experimental investigations of cubic piezoelectromagnetics are performed not frequently compared with those of the transversely isotropic materials.Therefore, phase velocity calculations for the new SH-wave propagating in concrete cubic piezoelectromagnetics will be not carried out in this report.Note that one can find a proposal of several new cubic piezoelectromagnetic composite materials in book [20].The following section provides the theoretical description of the problem.

Theory
Following the theoretical description of SH-SAW propagation in the cubic piezoelectromagnetics described in book [20], it is possible to briefly introduce the problem focusing on the SH-wave.First of all, it is necessary to state that in this case, the piezoelectromagnetic acoustic waves are managed along direction [101] like that in the cubic piezoelectrics [22].The theory uses the usual quasistatic approximations, the equations of motion, electrostatics and magnetostatics like those described in books [17,20].The coordinate beginning is positioned at the interface between the cubic piezoelectromagnetics and a vacuum.The rectangular coordinate x 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, φ are the electrical and magnetic potentials, respectively, defined by E = -φ/x 3 and H = -ψ/x 3 where E and H are th and ψ e 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  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 0 2,4,5

U
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 nary unity, respectively.The wavenumber k in E (4 k 1 (5) where U 0 = U 2 0 , φ 0 = U 4 0 , and ψ 0 = U 0 represent three components of the eigenvector imagiquation ) is coupled with the components {k 1 , k 2 , k 3 } of the wavevector k as follows: { , k 2 , 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 [32] for brevity, three homogeneous equations can be represented in the following matrix form: tain the explic Note th the readers themselves can ob it forms for the GL-tensor components in Equation (5) using plane wave solutions (4) for equations from (1) to (3).
alue n 3 = k at It is obvious that the three-component eigenvector should be nonzero for each eigenv 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: where m = 1 + sqr(n 3 ) and V t4 = sqrt(C/) 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: ) is the speed of the SH-BAW coupled with both the electrical and 2 ).It reads tion can be written as follows: Therefore, the left-hand side in Equation ( 7) can be written in the form of two factors.The transformed equa- 0) reveals the following eigenvalue n 3 : In the studied case of the wave propagation problem, th nts.It is apparent that the simplest and most convenient eigenvector ei (1)   The first factor in Equation (1 e first eigenvalue in expression ( 11) can be real.Using it for Equation (5), one can obtain the corresponding eigenvector compone components for the first genvalue n 3 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: propagation problem genvalue n 3 (2) , one can check that ze (14) 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 .Using the ei-Equation ( 5) equals to ro 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: (8).Therefore, this case is unique and demonstrates some association with the SH Also, it is apparent that the value of n 3 (3) becomes real as so (16) 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 -BAW solution.
on 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 [22], cubic piezomagnetics [21], and cubic piezoelectromagnetics below the SH-BAW velocity V tem [20].Therefore, it is possible to state that K 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: 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 in the cubic piezoelectromagnetics, the first eigenvalue is real and defined by expression (11), but the second and th blem of wave propagation in piezoel (17) 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 propagation ird 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 [33] have theoretically studied the two-phase materials and discussed possible boundary conditions for the pro ectromagnetics.In this case of the wave propagation, it is also possible to use the mechanical boundary condition such as the normal component of the stress tensor at the solid-vacuum interface x 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 [20].The other possible electrical and magnetic boundary conditions are as follows: electrically open surface (D = 0) and the magnetically closed surface (B = 0) where D and B are the normal components of the electrical and magnetic displacements.The mechanical boundary condition used in all the theoretical treatments below is applied to the normal component of the stress tensor, namely σ 32 (x 3 = 0) = 0 at the interface x 3 = 0 [20] between the crystal surface and a vacuum.This is the case called the mechanically free surface.It is possible to consider possible boundary conditions.

σ 32 = 0, φ = 0, and ψ = 0
Following the theoretical treatments done in book [20], the utilization of these boundary conditions results in the following system of three homogeneous equations written in the matrix form: 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 a corresponding eigenvector components.Exploiting them and the first set of the eigenvector components by expression (17), one can obtain the explicit form for the first third-order boundary-condition determinan D3) of the coefficient matrix in Equation (19).It nd the defined t (BCis 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 th side is formed by two factors.Therefore, there are two po electromagnetics: e left-hand ssibilities 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 piezo 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 (9) is i n the previous section.As a result, F 2 = -F 3 together with F 1 = 0 will zero the complet ments defined by formulae from (20)  nappropriate here.This is true because n 3 (K em = 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 i e displace-).So, this ssibility 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 dam wards the depth of the cubic piezoelectromagnetics.In expressions from ( 27) to ( 29), the speed of the n V new is defined by expression (8).Using definit 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 wa propagation specifics exhibits a compensative character, s well-known that acoustic w rol p toew wave ion (17) ve 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 i ave 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 [34] are well-known and the SOPPs on crystal substrates have applications in microscopy, biosensing, and photonics.The SOPPs represent electromagnetic surface modes formed by the strong coupling of light and optical phonons in polar crystals, and are generally excited using IR or THz radiation.Generation and cont of the SOPPs are essential for realizing novel applications in microscopy, data storage, thermal emission, or in the field of metamaterials.However, a comparatively little attention is still focused on the SOPPs despite certain advantages including their ability to be generated in a wide spectral range, from IR to THz wavelengths at the surface of a large variety of semiconductors, insulators, and ferroelectrics [34].It was also mentioned in Ref. [35] that surface electromagnetic effects could enhance the efficiency of numerous physical and chemical processes [36], as these effects lead to an increase of the electromagnetic fields at the surface, giving rise to an improved experimental sensitivity.
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: re the weight factor F 3 from definitions (32) and the seco tor components given in expressio is y conditions for comparison.

σ 32 = 0, D = 0, and B = 0
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), whe nd set of the eigenvecn (18) must be used.It thought that it is vital to treat the other possible sets of the boundar Using this set of the mechanical, electrical, and magnetic boundary conditions, it is possible to write the following homogeneous equations [20]: n CU e h n eU ) 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 h e . 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: to mention the second solution for Equation (35).This is as follows: it is us s (33), one can get the following BCD3: (36) The weight factors are used in equations from (27) to (29), too.It is also necessary This second solution (37) looks like unreal or it is hard to realize it.However, it is thought that eful to record any solution to have a more complete picture of the problem of wave propagation.
Using eigenvector components (18) for Equation It is obvious that one can transform determinant (38) by the same way, see the transformations of determinant Copyright © 2012 SciRes.OJA A. A. ZAKHARENKO 110 (35) written above.It was found that determinant (38) equals to zero when h e    ed in Sub (31).This also coincides with the result obtain section 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 results.echan m the corresponding three homogeneous equations in the matrix form (this form al same 2.3.σ 32 = 0, φ = 0, and B = 0 This is the case of the m ically 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 for can be so found in book [20]) similar to Equations ( 19) and (33).Also, the reader can use definitions ( 17) and ( 18) in order to get two different BCDs of the coefficient matrix.Using expression (17), it is possible to expand the corresponding first BCD3 and to obtain the following equation: value of K em 2 is completely defined by the piezoelectric properties.This can be the case of the dominan electric phase.However, it is thought that this dominance is not obligatory in the case of a large value of α.Indeed, the follo where the nondimensional value of K α 2 was introduced in the recent book cited in Ref. [20], K α 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 t piezothe value of K em 2 can be large, but less than unity.Consequently, equality (40) can be fulfilled when wing occurs: 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 e the new wave propagation for this set of the boundary co incides with that revealed in formula ( 31 = -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: 2.4.σ 32 = 0, D = 0, and ψ = 0 In this case the mechanically free (σ = 0), electrically op s 2.2 and 2.1, respectively.Indeed, the reader can also form three homogene-of 32 en (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 Subsection ous equations in the matrix form and of the coefficient matrix by the Subsections 2.1 and 2.2.Using definition (17) to form th ce conditions.All of them are unsuitable similar to those for the second Subsection 2.3.They read: expand the BCD3 way demonstrated in e first BCD3 and expanding the BCD3, one can also obtain three limit existen BCD3 from h e    (25), K em 1 is inapt because it sed above and results 2 = 0, and K e 2 → K em 2 → 1.Indeed, K em 2 = gives the case of V new = V tem discus 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, propagates with the speed defined by re for small values of α, the value of K 2 can be very large.Th es of the material constant e in order to get the following: 1 similar to 0 < K 2 < 1.It is lucid that for v 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 fo occurs: It is expected that in the ca of large values of α and e, Equality (48) with K 2 < 1 can be propagate in the transversely isotropic piezoelectromagnetic materials.It is worth noting that as pply the theory developed in this work the new wave also lation (8).Indeed, α erefore, it is necessary to cope with small valu 0 < K α 2 < ery small em 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 llowing se α 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.

Discussion
These theoretical investigations carried out above soundly demonstrated that an additional surface wave can propagate in direction [101] in the cubic piezoelectromagnetics.This new wave has a unique propagation characteristic such as follows: Its speed is higher than the SH-BAW speed.This fact can relate this type of wave to leaky sur-face SH-waves.However, it is actually a surface wave.As a result, it is possible to state that in the cubic piezoelectromagnetics, nine different surface waves (including such SH-SAW as the surface BGM-wave) can propagate in the same suitable propagation direction, using different boundary conditions.For comparison, ten different SH-SAWs under application of the same possible boundary conditions can soon as one will a to the transversely isotropic piezoelectromagnetics to discover any SH-LSAW, one can possibly find that in such materials the suitable propagation speed is equal to zero.This can mean that such SH-SAW cannot be found in the transversely isotropic piezoelectromagnetics.However, it is possible that some types of leaky waves can exist in the materials.
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 h e    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 ma-trix is used in 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 which Formula (31), namely h e    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 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 piezomagnetic weak pie em of the experimental repo s can be used as the suitable matrix to solute a very zoelectric 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 [37].The main probl rts in Refs.[37][38][39] is the fact that they do not provide the complete set of material constants for investigated unique composites.Ref. [40] reported a study of piezoelectric-piezomagnetic multilayer with ε < 0 and μ < 0 in which dielectric polariton and magnetic polariton can be simultaneously created.Also, it was recently stated that phonon-polaritons in a form of a band-like structure can exist in piezomagnetic superlattices (PMSL) with periodically up and down polarized domain structures [41] in which the piezomagnetic coefficients are periodically modulated.Also, some theoretical approaches exist which use material constants with real and imaginary parts (complex numbers) to describe wave propagation in multi-layered structures [42].Some works concerning the investigations of different metamaterials can be also found in Refs.[43,44].It is noted that Ref. [44] deals with three-dimensional bulk metamaterials.Therefore, it is thought that SH-BAWs and SH-SAWs can be investigated in bulk metamaterials.
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.[21] and [22], correspondingly.That is true because the piezoelectric constant e vanishes in the case of pure piezomagnetics.As a result, existence Condition (25)

 
 cannot be fulfilled because of h = 0 and e ≠ 0 in the case of pure piezoelectrics or h ≠ 0 and e = 0 in the case of pure piezomagnetics.Indeed, one can note that it is possible to require α = 0 for both the cases mentioned above to fulfill Equalities (25) and (31).However, one will cope here with zeroes on both sides of Equalities (25) and (31).For that reason, it is necessary to work with non-zero values of the material constants: h ≠ 0, e ≠ 0, and α ≠ 0.

Conclusion
This report represents the theoretical description of the new wave propagation in direction [101] in the cubic piezoelectromagnetics.The new wave can propagate with the speed higher than that for the SH-BAW coupled with both the electrical and magnetic potentials when 0 < K 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 condi-tions, 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 ficantly weaker piezoelectric phase.able existence conditions obtained in c Composit s," Journal of Materials S p. 1710-1714.76, pp.727-728.doi:10.1080/00150197608236711phase with a signi Indeed, if the suit 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.

2 e
):xistence conditions of nditions.It was found that three solutions (existence conditions) can exist in this case.The first existence conbecause it gives two equal eigenval n = 0 is inapues n (