Access the full text.

Sign up today, get DeepDyve free for 14 days.

Acoustics
, Volume 1 (4) – Sep 20, 2019

/lp/multidisciplinary-digital-publishing-institute/on-discovery-of-the-twelfth-and-thirteenth-new-nondispersive-sh-saws-0Qz8kZ04BT

- Publisher
- Multidisciplinary Digital Publishing Institute
- Copyright
- © 1996-2021 MDPI (Basel, Switzerland) unless otherwise stated Disclaimer The statements, opinions and data contained in the journals are solely those of the individual authors and contributors and not of the publisher and the editor(s). MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Terms and Conditions Privacy Policy
- ISSN
- 2624-599X
- DOI
- 10.3390/acoustics1040044
- Publisher site
- See Article on Publisher Site

acoustics Article On Discovery of the Twelfth and Thirteenth New Nondispersive SH-SAWs in 6 mm Magnetoelectroelastics Aleksey Anatolievich Zakharenko International Institute of Zakharenko Waves (IIZWs), Krasnoyarsk 660014, Russia; aazaaz@inbox.ru Received: 20 May 2019; Accepted: 18 September 2019; Published: 20 September 2019 Abstract: This report acquaints the reader with an extra two new shear-horizontal surface acoustic waves (SH-SAWs). These new SH-SAWs can propagate along the free surface of the transversely isotropic (6 mm) magnetoelectroelastic materials. These (composite) materials can simultaneously possess the piezoelectric, piezomagnetic, and magnetoelectric eects. Some competition among these eects can lead to suitable solutions found for the following three possible coupling mechanisms: e – h", e – h, " – . Here, the mechanically free interface between the solid and a vacuum was considered. This report discovers the twelfth (thirteenth) new SH-SAW for the magnetically closed (electrically open) case and continuity of both the normal component of the electrical (magnetic) displacement and the electrical (magnetic) potential when the coupling mechanism e – h" (e – h) works. The propagation velocities were obtained in explicit forms that take into account the contribution of the vacuum material parameters. The discovered waves were then graphically studied for the purpose of disclosing the dissipation phenomenon (the propagation velocity becomes imaginary) caused by the coupling with the vacuum properties. The obtained results can be useful for further investigations of interfacial and plate SH-waves, constitution of technical devices, nondestructive testing and evaluation, and application of some gravitational phenomena. Keywords: transversely isotropic piezoelectromagnetics; magnetoelectric eect; new nondispersive SH-waves; synergetics; complex systems 1. Introduction According to their book [1], Landau and Lifshits have predicted the existence of the linear magnetoelectric (ME) eect in 1956. Landau has demonstrated the problem of the existence of the ME eect in solids for Dsyaloshinskii. The latter has theoretically found and reported in a 1959 paper [2] that antiferromagnetic Cr O can possess the ME eect. One year later, Astrov [3] experimentally 2 3 conﬁrmed the existence of this eect in Cr O . The ME eect can be found in a class of ME solids called 2 3 piezoelectromagnetics (PEMs), also known as magnetoelectroelastics (MEEs). These continuous media can possess both the piezoelectric and piezomagnetic eects resulting in the existence of the ME eect. The ME solids relate to smart materials because any change in the magnetic subsystem can cause some change in the electrical subsystem via the mechanical subsystem, or vice versa. Up to today, various ME materials (composites) were found (created) and they were then investigated both experimentally and theoretically. However, only a few are suitable for commercial use. Some recent reviews on the ME materials and their applications are listed in References [4–8]. The reader can also ﬁnd about one hundred reviews on the subject in the literature because the smart ME materials can be used for dierent applications including spintronics, i.e., the future electronics without free charge carriers. The linear ME eect characterized by the nonsymmetric second-rank tensor of the electromagnetic constants is quite small. Therefore, creation of two-phase piezoelectric-piezomagnetic composites Acoustics 2019, 1, 749–762; doi:10.3390/acoustics1040044 www.mdpi.com/journal/acoustics Acoustics 2019, 1 750 is commercially required in order to enhance the value of the electromagnetic constant and to have samples with stronger ME eect. Concerning the ME composites, Scott has stated in his recent review [7] that the outstanding PEM composite PZT–Terfenol-D is competitive with some respects to SQUID superconducting devices [9] because such multiferroic composites have already accomplished the level of commercial technical devices, for instance, weak magnetic ﬁeld sensors. Also, it is worth noting that some created composite devices [10] can cost as much as a one cent coin and have linear dimensions signiﬁcantly smaller than the coin radius. The ME eect of the composites can be even several orders stronger than that for the ME monocrystals. Concerning the ME monocrystals, a recent review by Kimura [6] has stated that the unique monocrystal Sr Co Fe O was discovered in 2010. 3 2 24 41 This Z-type hexagonal ferrite actually possesses the realizable ME eect appropriate for practical and commercial employments. The piezoelectric, piezomagnetic, and magnetoelectric properties can be found simultaneously in piezoelectromagnetic composites or metamaterials. The metamaterials cannot be found in nature and represent special class of magnetoelectric composite materials. Extensive investigations of their properties begun at the beginning of the 21st century [11]. A handbook [12] combines review articles on physics and applications of metamaterials by authors from many international research centers. Such smart materials are already used today in medicine [13] and optical nanoantennas [14]. Also, smart piezoelectromagnetic composites can be used instead of conventional piezoelectric materials. The exploitation of piezoelectromagnetics instead of piezoelectrics is more preferable to experimentally generate anti-plane polarized (surface) acoustic waves with noncontact methods. For this purpose, the electromagnetic acoustic transducers (EMATs) as the famous noncontact method [15–17] can be used. Today there is a single review paper written by the author and published in 2013 [8]. This review touches upon the acoustic waves’ propagation with the anti-plane polarization along the free surface of hexagonal (6 mm) PEM solids. The propagation of these shear-horizontal surface acoustic waves (SH-SAWs) is studied when the ME eect is taken into account. Also, extra four new SH-SAWs were recently discovered [18,19]. Some of the new waves studied [20] can have a dramatic dependence on the ME eect. These four new SH-SAWs must be added to the known seven new SH-SAWs discovered by the author in a 2010 book [21] and three SH-SAWs by Melkumyan introduced in his 2007 paper [22]. The latter three SH-SAWs were also studied in [23,24]. They are naturally called the piezoelectric exchange surface Melkumyan (PEESM) wave, the piezomagnetic exchange surface Melkumyan (PMESM) wave, and the surface Bleustein-Gulyaev-Melkumyan (BGM) wave for convenience. The latter Melkumyan SH-SAW was called the BGM wave in order to have an association with the well-known BG-wave [25,26] existing in purely piezoelectric (or purely piezomagnetic) materials of class 6 mm. However, Bleustein and Gulyaev did not study any wave propagation in a PEM continuum possessing the ME eect. Note that according to Zakharenko, A.A. [8], theoretical work [21,27] has also conﬁrmed the existence of the aforementioned three Melkumyan SH-SAWs. In 2003, Chapter 7 of dissertational work [28] by Darinskii has discussed the wave existence in piezoelectric-piezomagnetic solids (piezoelectromagnetics). Chapter 7 was written following the theoretical investigations published in 1992 [29] and 1994 [30], respectively. Following the discussions [28], it is expected that at least two SAWs, namely one in-plane polarized wave and one anti-plane polarized wave can propagate in piezoelectromagnetics. This is similar to the wave propagation in the pure piezoelectrics or the pure piezomagnetics. Darinskii has also assumed that more than two SAWs can propagate in such media possessing the aforementioned three eects due to a competition among the eects. However, he has written down that such possibility was not recorded. As an example, a propagation velocity was obtained in an explicit form for the case of the mechanically clamped surface. This ﬁnal result depends on the coecient of the magnetoelectromechanical coupling (CMEMC) when the dependence on the electromagnetic constant was neglected, i.e., the ME eect was not taken into account. It is obvious that the problem of wave propagation in such smart piezoelectromagnetic materials relates to synergetics that studies various complex systems. Acoustics 2019, 1 751 In comparison, the theory developed by the author of this report utilizes the CMEMC depending on the electromagnetic constant . It is obvious that this material constantly participates in three possible CMEMC coupling mechanisms such as e – h", e – h, and " – . Besides the mass density % there are the following material constants for a piezoelectromagnetics: the elastic stiness constant C, piezoelectric constant e, piezomagnetic coecient h, dielectric permittivity coecient ", magnetic permeability coecient , and aforementioned electromagnetic constant . It is obvious that two ﬁrst aforementioned CMEMC coupling mechanisms represent some exchange mechanisms and do not pertain entirely to one of three eects mentioned above. It is possible to assume that only the third mechanism is relevant to the ME eect. Exploiting one of the coupling mechanisms, suitable solutions can be found. This was discussed in reference [31]. This theoretical report investigates the SH-wave propagation in the transversely isotropic (6 mm) PEM solids and adds the twelfth and thirteenth new SH-SAWs to the family of already known eleven SH-SAWs brieﬂy discussed above. To obtain the explicit forms for the new SH-wave velocities is the main purpose of the following three sections. 2. Theory Leading to an Extra Two New Results Consider a bulk PEM monocrystal or composite possessing the hexagonal symmetry of class 6 mm. Let us study the wave propagation in such transversely isotropic ME medium exploiting the rectangular coordinate system {x , x , x } also known as the Cartesian coordinate system invented 1 2 3 by Descartes, the famous mathematician. First of all, it is necessary to choose proper propagation direction in order to deal with the pure SH-wave coupled with both the electrical (') and magnetic ( ) potentials. The ﬁtting propagation direction for the 6 mm solid is described in reference [32]. For this purpose, the propagation direction is managed along the x -axis and perpendicular to the 6-fold symmetry axis. The x -axis is directed along the 6-fold symmetry axis and the x -axis is managed 2 3 along the normal to the free surface of the ME material. The coordinate beginning is situated at the interface between the solid and a vacuum and the x -axis negative values are managed towards the depth of the piezoelectromagnetics. It is very important to state that the problem of acoustic wave propagation coupled with both the electrical and magnetic potentials is treated. Therefore, the quasi-static approximation must be used here because the speed of any acoustic wave is ﬁve orders slower than the speed of the electromagnetic wave in a vacuum or solid. As a result of the chosen propagation direction, the dierential form of the coupled equations of motion [21,33–37] written in the common form can be separated into two independent sets of equations of motion. The ﬁrst set is for the purely mechanical wave with the in-plane polarization coupling two mechanical displacements U and U along the x -axis and x -axis, respectively. Besides, the second 1 3 1 3 set is for the mechanical SH-wave possessing the mechanical displacement U coupled with both the electrical and magnetic potentials. This theoretical report has an interest in the study of the second case. So, the dierential form of the coupled equations of motion representing the partial dierential equations of the second order can be composed as follows: 2 2 2 2 2 2 2 > @ @ @ @ @ U @ U @ U > C + + e + + h + = , > 2 2 2 2 2 2 2 > @x @x @x @x @x @x @t > 1 3 1 3 1 3 > 2 2 2 2 2 2 < @ @ @ @ @ U @ U e + " + + = 0, > . (1) 2 2 2 2 2 2 @x @x @x @x @x @x 3 3 3 1 1 1 2 2 2 2 > 2 2 @ @ @ @ > @ U @ U > h + + + = 0 2 2 2 2 2 2 @x @x @x @x @x @x 1 3 1 3 1 3 In coupled Equation (1), the material parameters such as %, C, e, h, ", , and are deﬁned in the previous section, the mechanical displacement along the x -axis is U = U , and t is time. It is natural to 2 2 write down the solutions of Equation (1) in the following plane wave form: U = U exp[j(k x + k x + k x !t)]. (2) I 1 1 2 2 3 3 I Acoustics 2019, 1 752 In Equation (2), the index I can be equal to 2, 4, 5. It is necessary to keep in mind that the electrical potential is ' = U and the magnetic potential is = U . The unknown coecients U , 4 5 0 0 1/2 ' , and must be determined. The imaginary unity is deﬁned by j = (– 1) and ! is the angular frequency. The wavevector components are written as follows: {k ,k ,k } = k{n ,n ,n }, where k is the 1 2 3 1 2 3 wavenumber in the propagation direction because the directional cosines are n = 1, n = 0, and n n . 1 2 3 3 A substitution of solutions (2) into Equation (1) leads to the tensor form of the coupled equations of motion known as the Green-Christoel equation. The tensor form then reads: 0 1 0 1 0 1 B C C m V /V em hm U 0 B CB C B C ph t4 B CB C B C B CB C B C B CB C B C B CB C B C B CB C = B 0 C. (3) B em "m m CB C B C B CB C B C B C@ A @ A @ A hm m m 1/2 In Equations’ set (3), V = !/k stands for the phase velocity, V = (C/) denotes the purely ph t 4 mechanical bulk acoustic wave (BAW) with the shear-horizontal (SH) polarization, and m = 1 + n . All the suitable eigenvalues n must be determined and the corresponding eigenvector components U , 0 0 ' , and must be also obtained by resolving the set of three homogeneous Equations (3). For this purpose, an expansion of the determinant of the coecient matrix in Equation (3) leads to a sixth order polynomial in one indeterminate n . Therefore, six suitable eigenvalues n must be found. They read 3 3 as follows: (1,2) (3,4) n = n = j, (4) 3 3 (5,6) n = j 1 V /V , (5) tem ph where 1/2 V = C/ 1 + K . (6) tem em In Equation (6), V represents the velocity of the SH-BAW coupled with both the electrical and tem magnetic potentials via the following coecient of the magnetoelectromechanical coupling (CMEMC): 2 2 e + "h 2eh e(e h) h(e h") eM hM 2 1 K = = = . (7) em 2 2 CM C(" ) C(" ) In Equation (7), it is convenient to distinguish the following three CMEMC coupling mechanisms discussed in reference [31]: M = e h", (8) M = e h, (9) M = " . (10) For the study of the surface wave propagation in the solid, it is necessary to leave only three suitable polynomial roots of six to deal with the wave damping towards the depth of the piezoelectromagnetics. For the coordinate system described above at the beginning of this section, the suitable values of n must have a negative sign. Let us use the ﬁrst, third, and ﬁfth eigenvalues as the suitable ones. Therefore, it is necessary to ﬁnd now the eigenvector components by substituting the suitable eigenvalues into Equation (3). The problem of the suitability of the eigenvectors was also recently discussed in reference [38]. The reader can do the same and ﬁnd the following forms of the ﬁrst set of the eigenvector components: 0 1 0 1 0 1 0(1) 0(3) B U C B U C B 0 C B C B C B C B C B C B C B C B C B C 0(1) 0(3) B C B C B C B C = B C = B C, (11) B C B C B C B C B C B C @ A @ A @ A 0(1) 0(3) " Acoustics 2019, 1 753 0 1 0 1 0 1 0 1 0(5) 2 B M /C C B U C B M /CK C 1 B M C B C 3 1 em B C B C B C B C B C B C B C 1 M B C B C B C 1 B C B 2 2 C 0(5) 2 B C B C B C B K K C B C = B eh/CK C = = B M C, (12) B em C 2 B C B em C B C 2 B C 2 B C B C B C B C K CK M @ A @ A @ A 0(5) 2 2 @ A 3 em 2 2 em " + e /CK M " K K em 1 em e because (e h")(e h) M M 1 2 2 2 K K = = , (13) em C(" ) CM (e h") 2 2 K K = = , (14) em e C"(" ) C"M eh eh K = = , (15) C C K = . (16) C" In eigenvector components (12), the nondimensional parameter K is deﬁned by Equation (15) and combines the terms with the electromagnetic constant in the CMEMC. Also, the other nondimensional parameter K called the coecient of the electromechanical coupling (CEMC) is deﬁned by Equation (16). The reader can also ﬁnd the common factors before the eigenvector components in Equation (12). It is necessary to leave them because the dimension of corresponding eigenvector components of dierent eigenvectors must be the same and they can play a crucial role for some ﬁnal results of obtaining propagation velocity. Indeed, Equation (12) demonstrates three dierent forms of the eigenvector. However, they are actually identical. It is thought that it is convenient to utilize the ﬁrst form when dierent determinants of the boundary conditions are constructed. They will be used in the following sections. The second form of the eigenvector in Equation (12) uses the nondimensional parameters 2 2 2 K , K , and K deﬁned above that can demonstrate a physical sense of the components. The third e em form already utilizes the CMEMC coupling mechanisms introduced after the CMEMC K (7). Other em dierent forms of the eigenvector components are discussed in reference [38] stating that forms (12) are preferable. It is possible to brieﬂy note here that eigenvectors (11) are in general uncertain. To resolve this uncertainty, the nonzero eigenvector components (11) can be connected with the corresponding eigenvector components (12) because eigenvector (12) is certain. So, one has to be familiar that there is the connection among the corresponding eigenvector components via the ﬁrst CMEMC coupling mechanism. This reads as follows: 0(1) 0(1) 0(3) 0(3) 0(5) 0(5) e + h = e + h = e + h = e h" = M . (17) There is also the second set of the eigenvector components for the same eigenvalues deﬁned by Equations (4) and (5). These components can be introduced as follows: 0 1 0 1 0 1 0(1) 0(3) B U C B U C B 0 C B C B C B C B C B C B C B C B C B C 0(1) 0(3) B C B C B C B C = B C = B C, (18) B C B C B C B C B C B C @ A @ A @ A 0(1) 0(3) 0 1 0 1 0 1 0 1 0(5) 2 B M /C C B U C B M /CK C 2 B M C 2 B C 3 B C B em C B C B C B C B C B C 1 B C M B C B C 2 2 2 B C 0(5) 2 2 B C B C B C B C B K K C B C = B h /CK C = = B M C, (19) B em m C em 2 B C B C B C 2 B C 2 B C B C B C B C @ A @ A K CK M @ A 0(5) 2 em @ A em 3 2 2 + eh/CK M K K em em where (e h) 2 2 K K = = , (20) em m C(" ) CM K = . (21) C Acoustics 2019, 1 754 In eigenvector components (19) there is the nondimensional parameter K called the coecient of the magnetomechanical coupling (CMMC) deﬁned by Equation (21). The connection of the corresponding eigenvector components is plausible through the following second CMEMC coupling mechanism: 0(1) 0(1) 0(3) 0(3) 0(5) 0(5) e + h = e + h = e + h = e h = M . (22) Using Equation (2) and the found eigenvalues with the eigenvector components, the complete S S parameters such as the complete mechanical displacement (U = U ), complete electrical potential S S S S (U = ), and complete magnetic potential (U = ) can be naturally written as follows: 4 5 0(p) (p) S (p) U = F U exp [jk(x + n x V t)], (23) p=1,3,5 1 3 ph I I 3 where the index I runs its values of 2, 4, 5. It is possible to ﬁnally state that the studied case relates to the three-partial wave problem because (1) (3) (5) each complete parameter (23) depends on three terms with the weight factors F , F , and F . Let us (1) (3) (5) further use F , F , and F instead of F , F , and F , respectively. Also, it is natural to further 1 2 3 exploit F = F + F because there are two identical eigenvalues (4) that give two identical eigenvectors 1 2 (11) or (18). Exploiting complete parameters (23), various determinants of the boundary conditions can be composed. The possible forms of such determinants depend on the mechanical, electrical, and magnetic boundary conditions at the interface between a vacuum and the piezoelectromagnetic solid. Several decades ago Al’shits, Darinskii, and Lothe [29] perfectly described possible boundary conditions for the treated case of the wave propagation along the surface of the piezoelectromagnetics. The used mechanical boundary condition at the interface relates to the normal component of the stress tensor: = 0. The electrical boundary conditions can be the continuity of the normal component of the electrical displacement D , continuity of the electrical potential ', D = 0, or ' = 0. Besides, 3 3 the possible magnetic boundary conditions are the continuity of the normal component of the magnetic displacement B , continuity of the magnetic potential , B = 0, or = 0. It is reasonable that all the 3 3 boundary conditions will not be given in this report in their explicit ﬁnal forms because the reader can ﬁnd them in book [21] and open access publication [18,19]. Also, the reader must be familiar with the following vacuum parameters: the dielectric permittivity 7 2 12 8 constant is " = 10 /(4C ) = 8.854187817 10 (F/m), where C = 2.99782458 10 (m/s) is the speed 0 L L of light in a vacuum. The magnetic permeability constant is = 4 10 (H/m) = 12.5663706144 7 f 10 (H/m). The constant " is the proportionality coecient between the electric induction D and f f f the electric ﬁeld E : D = " E , where the superscript “f ” is used for the free space (vacuum) and the f f electric ﬁeld components can be deﬁned as follows: E = @' /x . Similarly, the constant is the i i 0 proportionality coecient between the vacuum magnetic induction B and the vacuum magnetic ﬁeld f f f f f H : B = H , where the magnetic ﬁeld components can be deﬁned as follows: H = @ /x . So, 0 i i f f Laplace’s equations such as D' = 0 and D = 0 can be used for the potentials in a vacuum. It is also required that both the potentials must exponentially vanish in a vacuum far from the free surface of the piezoelectromagnetics. Utilizing the vacuum parameters, let’s treat some additional new cases that were not recorded in book [21] and reference [18,19,22]. This is the main purpose of the following two sections. 3. The New Results for the Case of B = 0 and Continuity of Both D and ' 3 3 This is the case of the continuity of the electric displacement D . Therefore, the vacuum electric constant " introduced in the previous section must be included in the further calculations. Using the ﬁrst set of the eigenvector components given in Equations (11) and (12), the coupled three equations for this case (Equations (188), (189), and (190) in book [21]) are then rewritten as follows: Acoustics 2019, 1 755 2 2 (" + " )CK F + b(" + " )C 1 + K F = 0, 0 em 0 em 3 2 2 2 2 CK F + eh CK F = 0, (24) em em : 2 2 " + " CK F + 0F = 0, 0 3 em where b = 1 V /V . (25) ph tem It is necessary here to state that the set of Equations (24) already embodies three consistent equations in two unknowns because the third equation is the main equation and represents a sum of the rest two. It is unnecessary to acquaint the reader with all the relatively complicated mathematics leading to ﬁnal Equations’ set (24). The reader can read book [21] because the case of this section was treated in Chapter XIII of the book. Therefore, Equations’ set (24) epitomizes Equations (188), (189), and (190) in the book chapter when the vacuum electric constant " is taken into account. Therefore, the resulting propagation velocity exempliﬁes the velocity V of the ﬁfth new new SH-SAW discovered in book [21]: h i 1/2 V = V 1 b . (26) tem new5 n5 In Equation (26), the parameter b already includes the vacuum electric constant " that is missing n 5 0 in the book. Therefore, this parameter can be expressed as follows: 2 2 K K em b = . (27) n5 (" + " ) 1 + K em Setting " = 0, this parameter can be simpliﬁed and written similar to the result obtained in book [21] and studied in reference [23]: 2 2 2 2 K K K K em EM em b = = , (28) n5 2 2 2 1 + K V 1 + K em em 1/2 where V = (") and V = 1/ stand for the electromagnetic wave speed and the exchange EM speed, respectively. Also, it is crucial to state that ﬁnal results (27) and (28) relate to the third coupling mechanism of the CMEMC deﬁned by Equation (10). One can also check that result (28) is also true for the following electrical and magnetic boundary conditions: ' = 0 and = 0 that also supports the propagation of the surface Bleustein-Gulyaev-Melkumyan wave mentioned above in Introduction. It is obvious that zero electromagnetic constant, = 0, provides b = 0 in Equations (27) and (28) and V = V n 5 new 5 tem occurs in Equation (26). This means that the magnetoelectric eect is vital for the existence of the ﬁfth new SH-SAW. Indeed, as soon as <> 0, the ﬁfth new SH-SAW can propagate. This is similar to the eighth and tenth new SH-SAWs discovered in reference [18,19], respectively, and studied in reference [20]. However, the existence conditions for the latter two new SH-SAWs are more complicated. These SH-SAWs cannot propagate in some PEM (composite) solids of class 6 mm [20] and require large enough value of the electromagnetic constant . Generally, the value of the is very small for real piezoelectromagnetics that are available today for commercial aims. Next, it is possible to treat a new possible case that was not recorded in previous work [18,19,21]. This new case pertains to the ﬁrst coupling mechanism (8) of the CMEMC. Consequently, one has to ﬁttingly reform Equations’ set (24). To get the suitable new set of three consistent equations, the following corresponding factors must be employed for three Equations (24): [e h(" + " )]/(" + " ), e/, 0 0 andh(" + " )/ " + " . In the reformed Equations (24) with the factors, the ﬁrst equation is 0 0 already the main equation that represents a sum of the rest two. As a result, a successive subtraction of the last two equations from the ﬁrst one leads to the following equation for propagation velocity of the twelfth new SH-SAW: Acoustics 2019, 1 756 h i 1/2 V = V 1 b , (29) new12 tem n12 where 2 2 2 2 eV K K K K em EM em b = = . (30) n12 2 2 2 e h(" + " ) 1 + K (e h" )V h 1 + K 0 0 em em EM If one does not want to use the vacuum electric constant " in calculations setting " = 0, 0 0 Equation (30) reduces to the following form: 2 2 2 2 eV e K K K K em EM em b = = . (31) n12 2 2 2 e h" 1 + K eV h 1 + K em em EM It is indispensable to brieﬂy discuss that in this case, = 0 does not ban the SH-SAW existence. 2 2 Using Equation (7) for K and (15) for K , the reader can check this statement. Indeed, parameter em b ( = 0) reduces to the following form: 2 2 e K b ( = 0) = ! . (32) n12 2 2 2 2 " =0 0 1 + K + K ( ) C " + " 1 + K + K e m e m 4. The New Results for the Case of D = 0 and Continuity of Both B and 3 3 For this case, the continuity condition for the magnetic displacement B allows one to include the vacuum magnetic constant . Using the second set of the eigenvector components given in Equations (18) and (19), the three coupled equations (namely, Equations (185), (186), and (187) from book [21]) can be rewritten in the following forms: 2 2 "( + )CK F + b"( + )C 1 + K F = 0, > 0 0 3 em em h i 2 2 "( + ) CK F + 0F = 0, (33) 0 em 3 : 2 2 2 2 CK F eh CK F = 0. em em In Equations’ set (33), the second equation represents the main equation that can be formed by a sum of the ﬁrst and third equations. Taking into account the vacuum parameter , the velocity V new 0 5 of the ﬁfth new SH-SAW deﬁned by Equation (26) depends on the parameter b . This parameter b n 5 n 5 already includes the vacuum magnetic constant instead of " . So, the parameter b reads: 0 0 5 2 2 K K em b = . (34) n5 "( + ) 1 + K em To obtain the extra new result, it is necessary to multiply the three equations in set (33) by the following factors: [e( + ) h]/"( + ), e( + )/ " + " , and h/, respectively. 0 0 0 0 This transformation of three Equations (33) sets the ﬁrst equation as the main equation. Therefore, the second and third transformed equations must be successively subtracted from the ﬁrst. This leads to the existence of the thirteenth new result, namely the velocity V of the thirteenth new SH-SAW new 13 that can be written as follows: h i 1/2 V = V 1 b , (35) new13 tem n13 where 2 2 2 2 h"V K K K K EM em em b = = . (36) n13 2 2 2 e( + ) h ( ) 0 1 + K e + e h "V 1 + K em em EM If the reader does not take into account the small vacuum parameter in (36), the parameter b 0 n 13 reduces to the following form: Acoustics 2019, 1 757 2 2 2 2 h"V h K K K K em EM em b = = . (37) n13 2 2 2 e h 1 + K e h"V 1 + K em em EM It is also possible to inspect the limit case for the electromagnetic constant = 0. With = 0 in Equations (36) and (37), the parameter b reduces to the following form: n 13 2 2 ( ) b = 0 = ! . (38) n13 2 2 2 2 =0 0 1 + K + K C( + ) 1 + K + K e m e m Thus, it is possible to conclude that the magnetoelectric eect does not play a vital role in the existence of the 12th and 13th new SH-SAWs because = 0 does not lead to the lack of the instability of the SH-BAW: V , V due to b ( = 0) , 0 in Equation (32) and V , V due to b ( = 0) new 12 tem n 12 new 13 tem n 13 , 0 in Equation (38). It is obvious that some graphical investigations of the new results obtained in both this section and the previous sections are required. Also, it is possible to demonstrate the inﬂuence of the inclusion of the corresponding vacuum parameter (" or ) in the calculations. This is the main 0 0 purpose of the following section. 5. Numerical Results and Discussion It is necessary to numerically investigate the theoretical results obtained in the previous two sections. The obtained explicit forms for the new SH-SAW velocities are quite complicated. Therefore, they can be studied numerically with discussion in this section below. Obtained Equations (29) and (35) in the third and fourth sections, respectively, are ﬁnal and have suitable forms for further numerical study. For the calculations, the material parameters can be borrowed from reference [39–42]. Let us use the material parameters for two piezoelectromagnetic materials (composites BaTiO –CoFe O and PZT-5H–Terfenol-D) given in references [39,40]. The PZT-5H–Terfenol-D 3 2 4 10 2 2 10 material parameters are: C = 1.45 10 (N/m ), e = 8.5 [C/m ], h = 83.8 (T), " = 75.0 10 6 2 3 (F/m), = 2.61 10 (N/A ), = 8500 (kg/m ). This composite is known as one of the strongest piezoelectromagnetics. The composite material BaTiO –CoFe O possesses signiﬁcantly weaker 3 2 4 piezoelectromagnetic properties. Therefore, the reader has a contrast for comparison and better understanding of the problem of the SH-wave propagation in the piezoelectromagnetic materials. 10 2 2 The BaTiO –CoFe O material parameters are: C = 4.4 10 (N/m ), e = 5.8 (C/m ), h = 275.0 (T), " = 3 2 4 10 6 2 3 56.4 10 (F/m), = 81.0 10 (N/A ), = 5730 (kg/m ). First of all, it is natural to contemplate the velocity V of the ﬁfth new SH-SAW given by new 5 Equation (26) and discovered in book [21] several years ago. It depends on the parameter b deﬁned by Equations (27), (28), and (34). Equations (27) and (34) contain the constants " and , respectively. 0 0 Figure 1 shows the dependence of the normalized velocity V /V on the electromagnetic constant new tem , namely on the normalized parameter /". It is obvious that the value of V /V must be new 5 tem smaller than unity because one deals here with one of the surface acoustic waves. Also, the value of the 2 2 /" must be conﬁned between zero and unity to satisfy the limitation condition such as < " [4,5]. The insertion in Figure 1 demonstrates that at = 0 (i.e., /" = 0) there is V = V and then for new tem the larger values of the /", the SH-SAW speed V is slightly slower than the SH-BAW speed new 5 2 2 V . However there is V = V anew at /" ~ 0.1565 for BaTiO –CoFe O and at /" ~ 0.2793 tem new tem 5 3 2 4 for PZT-5H–Terfenol-D symbolizing the fact that there is no instability of the SH-BAW at such large enough values of /". Indeed, V = V can occur when b changes its sign at some certain value new tem n 5 5 of /" > 0 given above for both the studied solids. One can also read paper [23] on the study of the ﬁfth new SH-SAW, in which the parameter D = V – V was used instead of V /V . For ! tem new new tem 5 5 " there is V !1 due to K ! 1 and therefore, V >> V occurs resulting in V /V ! 0. tem tem new 5 new 5 tem em For the twelfth new SH-SAW deﬁned by Equations (29), (30), and (31) there can occur more complicated dependence of the normalized velocity V /V on the parameter /". This is shown new 12 tem in Figure 2. For the case of result (29) with (30) when the parameter " is included in the calculations, 0 Acoustics 2019, 2 757 numerical study. For the calculations, the material parameters can be borrowed from reference [39– 42]. Let us use the material parameters for two piezoelectromagnetic materials (composites BaTiO3– CoFe2O4 and PZT-5H–Terfenol-D) given in references [39,40]. The PZT-5H–Terfenol-D material 10 2 2 −10 −6 parameters are: C = 1.45 × 10 (N/m ), e = 8.5 [C/m ], h = 83.8 (T), ε = 75.0 × 10 (F/m), μ = 2.61 × 10 2 3 (N/A ), ρ = 8500 (kg/m ). This composite is known as one of the strongest piezoelectromagnetics. The composite material BaTiO3–CoFe2O4 possesses significantly weaker piezoelectromagnetic properties. Therefore, the reader has a contrast for comparison and better understanding of the problem of the SH-wave propagation in the piezoelectromagnetic materials. The BaTiO3–CoFe2O4 material 10 2 2 −10 −6 parameters are: C = 4.4 × 10 (N/m ), e = 5.8 (C/m ), h = 275.0 (T), ε = 56.4 × 10 (F/m), μ = 81.0 × 10 2 3 (N/A ), ρ = 5730 (kg/m ). First of all, it is natural to contemplate the velocity Vnew5 of the fifth new SH-SAW given by Acoustics 2019, 1 758 Equation (26) and discovered in book [21] several years ago. It depends on the parameter bn5 defined by Equations (27), (28), and (34). Equations (27) and (34) contain the constants ε0 and μ0, respectively. Figure 1 shows the dependence of the normalized velocity Vnew5/Vtem on the electromagnetic constant 2 2 α, namely on the normalized parameter α /εμ. It is obvious that the value of Vnew5/Vtem must be smaller the velocity V touches the SH-BAW speed V at the same nonzero value of /" that is given new 12 tem than unity because one deals here with one of the surface acoustic waves. Also, the value of the α /εμ in the context above for the case of V = V . After that right away, the speed V rapidly new tem new 5 12 must be confined between zero and unity to satisfy the limitation condition such as α < εμ [4,5]. The decreases down to zero and then even becomes imaginary, illuminating the dissipation phenomenon. insertion in Figure 1 demonstrates that at α = 0 (i.e. α /εμ = 0) there is Vnew5 = Vtem and then for the 2 2 The velocity V stays imaginary for the values of /" in the narrow -range between /" ~ largernew values of the α /εμ, the SH-SAW speed Vnew5 is slightly slower than the SH-BAW speed Vtem. 2 2 However there is Vnew5 = Vtem anew at α /εμ ~ 0.1565 for BaTiO3–CoFe2O4 and at α /εμ ~ 0.2793 for PZT- 0.1569 and ~0.1571 for BaTiO –CoFe O and between ~ 0.27975 and ~ 0.28046 for PZT-5H–Terfenol-D. 3 2 4 5H–Terfenol-D symbolizing the fact that there is no instability of the SH-BAW at such large enough These are quite narrow -ranges for the existence of the dissipation phenomenon due to taking into 2 2 values of α /εμ. Indeed, Vnew5 = Vtem can occur when bn5 changes its sign at some certain value of α /εμ account the vacuum electric constant " . It is expected that a narrow -range for the dissipation > 0 given above for both the studied solids. One can also read paper [23] on the study of the fifth new phenomenon can exist even for a very weak PEM material and this -range can be naturally shifted SH-SAW, in which the parameter Δ = Vtem – Vnew5 was used instead of Vnew5/Vtem. For α → εμ there is Vtem → ∞ due to 𝐾 →∞ and therefore, Vtem >> Vnew5 occurs resulting in Vnew5/Vtem → 0. towards smaller values of the normalized parameter /". Acoustics 2019, 2 758 2 2 velocity Vnew12 stays imaginary for the values of α /εμ in the narrow α-range between α /εμ ~ 0.1569 and ~ 0.1571 for BaTiO3–CoFe2O4 and between ~ 0.27975 and ~ 0.28046 for PZT-5H–Terfenol-D. These Figure 1. The inﬂuence of the inclusion of the vacuum parameters " and in the dependence of the Figure 1. The influence of the inclusion of the vacuum parameters ε0 and μ0 in the dependence of the 0 0 are quite narrow α-ranges for the existence of the dissipation phenomenon due to taking into account normalized velocity Vnew5/Vtem on the electromagnetic constant α, i.e. on the nondimensional parameter normalized velocity V /V on the electromagnetic constant , i.e., on the nondimensional parameter the vacuum electric con new stem tant ε0. It is expected that a narrow α-range for the dissipation phenomenon α /εμ. With Equation (26) using (28), the gray and black lines are calculated for BaTiO3–CoFe2O4 and can exist even for a very weak PEM material and this α-range can be naturally shifted towards smaller /". With Equation (26) using (28), the gray and black lines are calculated for BaTiO –CoFe O and 3 2 4 PZT-5H–Terfenol-D, respectively. The dotted lines are calculated by using (27) with ε0 and (34) with values of the normalized parameter α /εμ. PZT-5H–Terfenol-D, respectively. The dotted lines are calculated by using (27) with " and (34) with . 0 0 μ0. Note that all three cases (27), (28), and (34) coincide which is clearly seen in the inclusion. Note that all three cases (27), (28), and (34) coincide which is clearly seen in the inclusion. For the twelfth new SH-SAW defined by Equations (29), (30), and (31) there can occur more complicated dependence of the normalized velocity Vnew12/Vtem on the parameter α /εμ. This is shown in figure 2. For the case of result (29) with (30) when the parameter ε0 is included in the calculations, the velocity Vnew12 touches the SH-BAW speed Vtem at the same nonzero value of α /εμ that is given in the context above for the case of Vnew5 = Vtem. After that right away, the speed Vnew12 rapidly decreases down to zero and then even becomes imaginary, illuminating the dissipation phenomenon. The Figure 2. The normalized velocity V /V versus the nondimensional parameter 2 /". The gray new tem Figure 2. The normalized velocity V 12 new12/Vtem versus the nondimensional parameter α /εμ. The gray and black and black lines are line for s are for BaTiO BaTiO –CoFe 3–CoFe O 2O and 4 and PZT PZT-5H -5H–T –Terfenol-D, respectively erfenol-D, respectively . Using. Eq Using uationEquation (29), (29), 3 2 4 the dotted and solid lines are for the cases with the parameter ε0 in Equation (30) and without ε0 in the dotted and solid lines are for the cases with the parameter " in Equation (30) and without " in 0 0 Equation (31), respectively. The insertion shows the dissipation phenomenon by the dotted black line Equation (31), respectively. The insertion shows the dissipation phenomenon by the dotted black line at α /εμ ~ 0.28 for PZT-5H–Terfenol-D. at /" ~ 0.28 for PZT-5H–Terfenol-D. Using the dotted black line for PZT-5H–Terfenol-D shown in the insertion of Figure 2, it is Using pthe ossib dotted le to note th black at Vline new12 = for VtemPZT occurs w -5H–T hen erfenol-D bn12 = 0, i.e. b shown n12 defined in by Equ the insertion ation (30) chan of Figur ges its sign e 2, it .is possible Indeed, the value of bn12 rapidly approaches an infinity and also quickly returns from an infinity to note that V = V occurs when b = 0, i.e., b deﬁned by Equation (30) changes its sign. new 12 tem n 12 n 12 already with the changed sign. It is clearly seen in Figure 2 that these peculiarities of Vnew12 = Vtem and Indeed, the value of b rapidly approaches an inﬁnity and also quickly returns from an inﬁnity the imaginary velocity Vnew12 within the narrow α-range are absent for the case of Equation (29) with already with the changed sign. It is clearly seen in Figure 2 that these peculiarities of V = V new 12 tem (31). Also, both cases (30) and (31) have the imaginary values of the velocity Vnew12 at some large values 2 2 2 2 and the imaginary of α /εμ due velocity to large vV alues of within (bn12) > 1: the α /εμ narr > ~ 0.96 ow 826 -range for Baar TiO e3– absent CoFe2O4for and the α /εμ case > ~ of 0.85 Equation 933 (29) new for PZT-5H–Terfenol-D. Note that for the case without ε0 shown by the solid lines in Figure 2, this with (31). Also, both cases (30) and (31) have the imaginary values of the velocity V at some new 12 dissipation commences at slightly larger value of α /εμ for PZT-5H–Terfenol-D. Such small 2 2 2 large values of /" due to large values of (b ) > 1: /" > ~ 0.96826 for BaTiO –CoFe O and 12 3 2 4 discrepancy was not recorded for the weaker PEM composite BaTiO3–CoFe2O4. It is also necessary to mention that the dissipation phenomenon can exist at very small values of 2 2 α /εμ that are smaller than some threshold value (α /εμ)th for the 8th and 10th new SH-SAWs th th discovered in references [18,19], respectively. Some propagation peculiarities of the 8 and 10 new 2 2 SH-SAWs were studied in reference [20]. If α /εμ > (α /εμ)th, the 8th and 10th new SH-SAWs can propagate. The value of the electromagnetic constant α is very small for real piezoelectromagnetic 2 2 monocrystals and even composite materials. Therefore, the existence condition of α /εμ > (α /εμ)th for the 8th and 10th new SH-SAWs is extremely important. Also, the existence of the dissipation 2 2 phenomenon with the imaginary propagation velocity at a narrow α-range, 0 < α /εμ < (α /εμ)th, is th caused by the corresponding vacuum parameter, ε0 or μ0. Theoretical work [20] also states that the 8 and 10th new SH-SAWs are apt for constitution of technical devices to study (to sense) the magnetoelectric effect. Figure 3 shows the behavior of the normalized velocity Vnew13/Vtem in dependence on the nondimensional parameter α /εμ for composite materials BaTiO3–CoFe2O4 and PZT-5H–Terfenol-D. The figure insertion illuminates the dependence of Vnew13(α) for small values of α /εμ. It is clearly seen Acoustics 2019, 1 759 /" > ~ 0.85933 for PZT-5H–Terfenol-D. Note that for the case without " shown by the solid lines in Figure 2, this dissipation commences at slightly larger value of /" for PZT-5H–Terfenol-D. Such small discrepancy was not recorded for the weaker PEM composite BaTiO –CoFe O . 3 2 4 It is also necessary to mention that the dissipation phenomenon can exist at very small values 2 2 of /" that are smaller than some threshold value ( /") for the 8th and 10th new SH-SAWs th discovered in references [18,19], respectively. Some propagation peculiarities of the 8th and 10th new 2 2 SH-SAWs were studied in reference [20]. If /" > ( /") , the 8th and 10th new SH-SAWs can th propagate. The value of the electromagnetic constant is very small for real piezoelectromagnetic 2 2 monocrystals and even composite materials. Therefore, the existence condition of /" > ( /") th for the 8th and 10th new SH-SAWs is extremely important. Also, the existence of the dissipation 2 2 phenomenon with the imaginary propagation velocity at a narrow -range, 0 < /" < ( /") , th is caused by the corresponding vacuum parameter, " or . Theoretical work [20] also states that 0 0 the 8th and 10th new SH-SAWs are apt for constitution of technical devices to study (to sense) the magnetoelectric eect. Figure 3 shows the behavior of the normalized velocity V /V in dependence on the new 13 tem nondimensional parameter /" for composite materials BaTiO –CoFe O and PZT-5H–Terfenol-D. Acoustics 2019, 2 3 2 4 759 The ﬁgure insertion illuminates the dependence of V () for small values of /". It is clearly new 13 for one of the strongest PEM composites such as PZT-5H–Terfenol-D that the inclusion of the vacuum seen for one of the strongest PEM composites such as PZT-5H–Terfenol-D that the inclusion of the parameter μ0 in the calculations with Equation (35) and (36) can result in disappearance of the vacuum parameter in the calculations with Equations 2 (35) and (36) can result in disappearance of dissipation phenomenon existing at the large values of α /εμ ~ 1. For the significantly weaker PEM the dissipation composite phenomenon such as BaTiO existing 3–CoFe2O at 4, the this phenomen large values on cannot e of /" xist~ . A 1.lso, For the inser the signiﬁcantly tion clearly weaker demonstrates that at the very small values of α /εμ, the speed Vnew13(α → 0) is slower for the case of PEM composite such as BaTiO –CoFe O , this phenomenon cannot exist. Also, the insertion clearly 3 2 4 Equations (35) and (37) without the vacuum parameter μ0 in comparison with the case of Equations demonstrates that at the very small values of /", the speed V (! 0) is slower for the case of new 13 th (35) and (36). This means that the penetration depth for the 13 new SH-SAW is larger for the latter Equations (35) and (37) without the vacuum parameter in comparison with the case of Equations (35) case. 0 and (36). This means that the penetration depth for the 13th new SH-SAW is larger for the latter case. Figure 3. The normalized velocity V /V versus the nondimensional parameter 2 /". The gray Figure 3. The normalized velocity new V 13 new13/V tem tem versus the nondimensional parameter α /εμ. The gray and black lines are for BaTiO3–CoFe2O4 and PZT-5H–Terfenol-D, respectively. The dotted lines and black lines are for BaTiO –CoFe O and PZT-5H–Terfenol-D, respectively. The dotted lines 3 2 4 correspond to the calculations with Equation (35) when the parameter μ0 is included in Equation (36). correspond to the calculations with Equation (35) when the parameter is included in Equation (36). The solid lines are for case (37) without the parameter μ0. The insertion shows the dependence for The solid lines are for case (37) without the parameter . The insertion shows the dependence for small values of α /εμ. small values of /". It is well-known that experimentally measured values of the electromagnetic constant α are very It is well-known that experimentally measured values of the electromagnetic constant are very small, α << εμ [4] and [5]. This means that the aforementioned dissipation phenomenon cannot be reachable for all known PEM monocrystals and composites. Indeed, this research is primary and small, << " [4] and [5]. This means that the aforementioned dissipation phenomenon cannot be allows the further theoretical investigations of various SH-waves in PEM plates. In general, the plate reachable for all known PEM monocrystals and composites. Indeed, this research is primary and waves have the following peculiarity: at large values of the plate thickness the speed of the plate SH- allows the further theoretical investigations of various SH-waves in PEM plates. In general, the plate wave must approach the speed of the corresponding SH-SAW. Therefore, it is helpful to know the waves have the following peculiarity: at large values of the plate thickness the speed of the plate characteristics of SAWs before development of any theory of SH-wave propagation in plates. It is well-known that plates are used for further miniaturization of various technical devices with a higher SH-wave must approach the speed of the corresponding SH-SAW. Therefore, it is helpful to know the level of integration. Also, various SH-waves such as SH-SAWs and plate SH-waves can be characteristics of SAWs before development of any theory of SH-wave propagation in plates. It is significantly more sensitive to various chemicals. This fact is called for implementation in sensor well-known that plates are used for further miniaturization of various technical devices with a higher technologies. These SH-waves can be also used in nondestructive testing and evaluation of the certain level of intsur egration. faces of variou Also,s PE various M solid SH-waves s. such as SH-SAWs and plate SH-waves can be signiﬁcantly Finally, it is necessary here to mention the recently developed theory concerning the coexistence of the electromagnetic and gravitational phenomena [43,44]. It is obvious that the obtained Equation (29) and (35) for the calculation of the 12th and 13th new SH-SAWs can be also true for the case when gravitational phenomena are exploited instead of the electromagnetic phenomena. Indeed, it is possible to use the material parameters corresponding to the purely gravitational case. In this case, –1/2 speeds (29) and (35) will depend on the gravitational wave speed VGC = (γη) , where γ and η are the gravitic constant and cogravitic constant, respectively, instead of the electromagnetic wave speed VEM –1/2 = (εμ) . For the purely electromagnetic case, it is also possible to state that speeds (29) and (35) –1/2 2 2 depend on the speed of light in a vacuum, CL = (ε0μ0) because ε0 = 1/(μ0CL ) and μ0 = 1/(ε0CL ). In the purely gravitational case, there are the vacuum gravitic constant γ0 = 1/(η0CL ) and the vacuum Acoustics 2019, 1 760 more sensitive to various chemicals. This fact is called for implementation in sensor technologies. These SH-waves can be also used in nondestructive testing and evaluation of the certain surfaces of various PEM solids. Finally, it is necessary here to mention the recently developed theory concerning the coexistence of the electromagnetic and gravitational phenomena [43,44]. It is obvious that the obtained Equations (29) and (35) for the calculation of the 12th and 13th new SH-SAWs can be also true for the case when gravitational phenomena are exploited instead of the electromagnetic phenomena. Indeed, it is possible to use the material parameters corresponding to the purely gravitational case. In this case, speeds 1/2 (29) and (35) will depend on the gravitational wave speed V = ( ) , where and are the GC gravitic constant and cogravitic constant, respectively, instead of the electromagnetic wave speed 1/2 V = (") . For the purely electromagnetic case, it is also possible to state that speeds (29) and EM 1/2 2 (35) depend on the speed of light in a vacuum, C = (" ) because " = 1/( C ) and = L 0 0 0 0 L 0 2 2 1/(" C ). In the purely gravitational case, there are the vacuum gravitic constant = 1/( C ) and the 0 L 0 0 L 2 1/2 vacuum cogravitic constant = 1/( C ) due to C = ( ) . However, it is expected that sound 0 0 L L 0 0 experimental veriﬁcations of some gravitational phenomena can be released only in several decades. For instance, the theoretically predicted existence of gravitational waves by Albert Einstein in 1916 [45] was veriﬁed only in 2016 [46] with very expensive experiments carried out by an international group consisting of over 1000 researchers. This happens due to the weakness of the gravitational phenomena. However, these phenomena are extremely important. Therefore, it is expected that it is preferable for the ﬁrst time to deal with the electromagnetic phenomena regarding the experimental proof. In the case of the acoustic wave propagation in solids, taking into account both the gravitational and the electromagnetic phenomena can result in the dependence of the acoustic wave velocity on both the speeds V and V as well as on the speeds of the new fast waves that can propagate in both solids EM GC and a vacuum. In both continuous media the propagation speeds of the new fast waves can be thirteen orders faster than the speed of light. The propagation of these new fast waves [47,48] in a vacuum represents a great interest due to a possibility to develop the instant interplanetary communication [47]. For this purpose, it is necessary to focus on solutions of many theoretical, mathematical, experimental, and engineering problems that can lead to the development of perfect infrastructures for the instant interplanetary communication. 6. Conclusions This theoretical study has obtained the explicit forms for the propagation velocities of the 12th and 13th new SH-SAWs that can propagate in the transversely isotropic piezoelectromagnetics of symmetry class 6 mm. The discovered new SH-waves correspond to dierent electrical and magnetic boundary conditions for the electrical and magnetic potentials and the normal components of the electrical and magnetic inductions. This investigation also has an interest in the record of the inﬂuence of the inclusion of the corresponding vacuum parameter, " or , on the propagation velocity of the 0 0 new SH-SAWs. For the ﬁfth new SH-SAW discovered in book [21], it was found that the inﬂuence of the inclusion of the " or is not signiﬁcant and these vacuum parameters can be neglected in 0 0 calculations. However, this is not true for the 12th and 13th new SH-SAWs discovered in this paper. It was found that with the vacuum electric constant " for the 12th new SH-SAW there can be found one extra narrow -range for /" << 1, in which the propagation velocity becomes imaginary, i.e., there is some dissipation. This dissipation phenomenon can also exist at large values of the nondimensional parameter /" < ~ 1. For the 13th new SH-SAW there is the other picture: the presence of the vacuum magnetic constant can result in the disappearance of the dissipation phenomenon for the values of the /" being just below unity. This was found only for PZT-5H–Terfenol-D representing one of the strongest PEM solids. The obtained theoretical results can be useful for the further investigations of the SH-wave propagation in PEM plates, nondestructive testing and evaluation, constitution of various technical devices based on PEM SH-SAWs. Also, some gravitational Acoustics 2019, 1 761 phenomena discussed in recently developed theory [43] and [44] can be applied as a supplementary to the electromagnetic properties or instead of them. Funding: This research received no external funding. Conﬂicts of Interest: The author declares no conﬂict of interest. References 1. Landau, L.D.; Lifshits, E.M. Course of Theoretical Physics. Electrodynamics of Continuous Media, 1st ed.; Gostekhizdat: Moscow, Russia, 1957; Volume 8, p. 532. 2. Dzyaloshinskii, I.E. On the magneto-electrical eects in antiferromagnetics. Sov. Phys. J. Exp. Theor. Phys. 1960, 10, 628–629. 3. Astrov, D.N. The magnetoelectric eect in antiferromagnets. Sov. Phys. J. Exp. Theor. Phys. 1960, 11, 708–709. 4. Fiebig, M. Revival of the magnetoelectric eect. J. Phys. D Appl. Phys. 2005, 38, R123–R152. [CrossRef] 5. Özgür, Ü.; Alivov, Y.; Morkoç, H. Microwave ferrites, Part 2: Passive components and electrical tuning. J. Mater. Sci. Mater. Electron. 2009, 20, 911–952. [CrossRef] 6. Kimura, T. Magnetoelectric hexaferrites. Annu. Rev. Condens. Matter Phys. 2012, 3, 93–110. [CrossRef] 7. Scott, J.F. Prospects for ferroelectrics: 2012–2022. ISRN Mater. Sci. 2013, 2013, 24. [CrossRef] 8. Zakharenko, A.A. Piezoelectromagnetic SH-SAWs: A review. Can. J. Pure Appl. Sci. 2013, 7, 2227–2240. [CrossRef] 9. Nan, C.W.; Bichurin, M.I.; Dong, S.; Viehland, D.; Srinivasan, G. Multiferroic magnetoelectric composites: Historical perspective, status, and future directions. J. Appl. Phys. 2008, 103, 031101. [CrossRef] 10. Israel, C.; Mathur, N.D.; Scott, J.F. A one-cent room-temperature magnetoelectric sensor. Nat. Mater. 2008, 7, 93–94. [CrossRef] [PubMed] 11. Vendik, I.B.; Vendik, O.G. Metamaterials and their utilization in techniques of ultrahigh frequencies: A review. J. Tech. Phys. 2013, 83, 3–28. 12. Capolino, F. (Ed.) Phenomena and Theory of Metamaterials; CRC Press: Boca Raton, FL, USA, 2009; Volume 1, p. 926. 13. Radu, X.; Garray, D.; Craeye, C. Toward a wire medium endoscope for MRI imaging. Metamaterials 2009, 3, 90–99. [CrossRef] 14. Filonov, D.S.; Krasnok, A.E.; Slobozhanyuk, A.P.; Kapitanova, P.V.; Nenasheva, E.A.; Kivshar, Y.S.; Belov, P.A. Experimental veriﬁcation of the concept of all-dielectric nanoantennas. Appl. Phys. Lett. 2012, 100, 2011131. [CrossRef] 15. Thompson, R.B. Physical principles of measurements with EMAT transducers. In Physical Acoustics; Mason, W.P., Thurston, R.N., Eds.; Academic Press: New York, NY, USA, 1990; Volume 19, pp. 157–200. 16. Hirao, M.; Ogi, H. EMATs for Science and Industry: Noncontacting Ultrasonic Measurements; Kuwer Academic: Boston, MA, USA, 2003. 17. Ribichini, R.; Cegla, F.; Nagy, P.B.; Cawley, P. Quantitative modeling of the transduction of electromagnetic acoustic transducers operating on ferromagnetic media. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 2010, 57, 2808–2817. [CrossRef] [PubMed] 18. Zakharenko, A.A. New nondispersive SH-SAWs guided by the surface of piezoelectromagnetics. Can. J. Pure Appl. Sci. 2013, 7, 2557–2570. [CrossRef] 19. Zakharenko, A.A. A study of new nondispersive SH-SAWs in magnetoelectroelastic medium of symmetry class 6 mm. Open J. Acoust. 2015, 5, 95–111. [CrossRef] 20. Zakharenko, A.A. Dramatic inﬂuence of the magnetoelectric eect on the existence of the new SH-SAWs propagating in magnetoelectroelastic composites. Open J. Acoust. 2015, 5, 73–87. [CrossRef] 21. Zakharenko, A.A. Propagation of Seven New SH-SAWs in Piezoelectromagnetics of Class 6 mm; Lambert Academic Publishing (LAP): Riga, Latvia, 2010; p. 84. 22. Melkumyan, A. Twelve shear surface waves guided by clamped/free boundaries in magneto-electro-elastic materials. Int. J. Solids Struct. 2007, 44, 3594–3599. [CrossRef] 23. Zakharenko, A.A. On wave characteristics of piezoelectromagnetics. Pramana J. Phys. 2012, 79, 275–285. [CrossRef] 24. Zakharenko, A.A. Analytical investigation of surface wave characteristics of piezoelectromagnetics of class 6 mm. ISRN Appl. Math. 2011, 2011, 8. [CrossRef] Acoustics 2019, 1 762 25. Gulyaev, Y.V. Electroacoustic surface waves in solids. Sov. Phys. J. Exp. Theor. Phys. Lett. 1969, 9, 37–38. 26. Bleustein, J.L. A new surface wave in piezoelectric materials. Appl. Phys. Lett. 1968, 13, 412–413. [CrossRef] 27. Wei, W.-Y.; Liu, J.-X.; Fang, D.-N. Existence of shear horizontal surface waves in a magneto-electro-elastic material. Chin. Phys. Lett. 2009, 26, 3. [CrossRef] 28. Darinskii, A.N. Waves in Finite Homogeneous and Layered Inhomogeneous Media with Arbitrary Anisotropy. Ph.D. Thesis, Institute of Crystallography, Russian Academy of Sciences, Moscow, Russia, 2003; pp. 127–148. (In Russian) 29. Al’shits, V.I.; Darinskii, A.N.; Lothe, J. On the existence of surface waves in half-inﬁnite anisotropic elastic media with piezoelectric and piezomagnetic properties. Wave Motion 1992, 16, 265–283. [CrossRef] 30. Al’shits, V.I.; Darinskii, A.N.; Lothe, J.; Lyubimov, V.N. Surface acoustic waves in piezocrystals: An example of surface wave existence with clamped boundary. Wave Motion 1994, 19, 113–123. [CrossRef] 31. Zakharenko, A.A. Peculiarities study of acoustic waves’ propagation in piezoelectromagnetic (composite) materials. Can. J. Pure Appl. Sci. 2013, 7, 2459–2461. [CrossRef] 32. Gulyaev, Y.V. Review of shear surface acoustic waves in solids. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 1998, 45, 935–938. [CrossRef] [PubMed] 33. Zakharenko, A.A. Thirty Two New SH-Waves Propagating in PEM Plates of Class 6 mm; Lambert Academic Publishing (LAP): Riga, Latvia, 2012; p. 162. 34. Pan, E. Exact solution for simply supported and multilayered magneto-electro-elastic plates. J. Appl. Mech. 2001, 68, 608–618. [CrossRef] 35. Annigeri, A.R.; Ganesan, N.; Swarnamani, S. Free vibrations of simply supported layered and multiphase magneto-electro-elastic cylindrical shells. Smart Mater. Struct. 2006, 15, 459–467. [CrossRef] 36. Zakharenko, A.A. Twenty Two New Interfacial SH-Waves in Dissimilar PEMs; Lambert Academic Publishing (LAP): Riga, Latvia, 2012; p. 148. [CrossRef] 37. Berlincourt, D.A.; Curran, D.R.; Jae, H. Piezoelectric and piezomagnetic materials and their function in transducers. In Physical Acoustics; Mason, W.P., Ed.; Academic Press: New York, NY, USA, 1964; Volume 1A, pp. 169–270. 38. Zakharenko, A.A. Some problems of ﬁnding of eigenvalues and eigenvectors for SH-wave propagation in transversely isotropic piezoelectromagnetics. Can. J. Pure Appl. Sci. 2014, 8, 2783–2787. [CrossRef] 39. Zakharenko, A.A. Investigation of SH-wave fundamental modes in piezoelectromagnetic plate: Electrically closed and magnetically closed boundary conditions. Open J. Acoust. 2014, 4, 90–97. [CrossRef] 40. Zakharenko, A.A. Consideration of SH-wave fundamental modes in piezoelectromagnetic plate: Electrically open and magnetically open boundary conditions. Waves Random Complex Media 2013, 23, 373–382. [CrossRef] 41. Wang, B.-L.; Mai, Y.-W. Applicability of the crack-face electromagnetic boundary conditions for fracture of magnetoelectroelastic materials. Int. J. Solids Struct. 2007, 44, 387–398. [CrossRef] 42. Liu, T.J.-C.; Chue, C.-H. On the singularities in a bimaterial magneto-electro-elastic composite wedge under antiplane deformation. Compos. Struct. 2006, 72, 254–265. [CrossRef] 43. Zakharenko, A.A. On piezogravitocogravitoelectromagnetic shear-horizontal acoustic waves. Can. J. Pure Appl. Sci. 2016, 10, 4011–4028. [CrossRef] 44. Zakharenko, A.A. The problem of ﬁnding of eigenvectors for 4P-SH-SAW propagation in 6 mm media. Can. J. Pure Appl. Sci. 2017, 11, 4103–4119. [CrossRef] 45. Einstein, A. Die Grundlage der allgemeinen Relativitätstheorie. Ann. Phys. 1916, 4, 284. [CrossRef] 46. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Abernathy, M.R.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; et al. Observation of gravitational waves from a binary black hole merger. Phys. Rev. Lett. 2016, 116, 16. [CrossRef] [PubMed] 47. Zakharenko, A.A. On necessity of development of instant interplanetary telecommunication based on some gravitational phenomena for remote medical diagnostics and treatment. Can. J. Pure Appl. Sci. 2018, 12, 4481–4487. [CrossRef] 48. Zakharenko, A.A. On separation of exchange terms for four-potential acoustic SH-wave case with dependence on gravitational parameters. Hadron. J. 2018, 41, 349–370. [CrossRef] © 2019 by the author. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Acoustics – Multidisciplinary Digital Publishing Institute

**Published: ** Sep 20, 2019

**Keywords: **transversely isotropic piezoelectromagnetics; magnetoelectric effect; new nondispersive SH-waves; synergetics; complex systems

Loading...

You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!

Read and print from thousands of top scholarly journals.

System error. Please try again!

Already have an account? Log in

Bookmark this article. You can see your Bookmarks on your DeepDyve Library.

To save an article, **log in** first, or **sign up** for a DeepDyve account if you don’t already have one.

Copy and paste the desired citation format or use the link below to download a file formatted for EndNote

Access the full text.

Sign up today, get DeepDyve free for 14 days.

All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.