First Evidence of Surface SH-Wave Propagation in Cubic Piezomagnetics

This theoretical work provides with results of characteristics calculation of the ultrasonic surface Zakharenko waves (USZWs) existing in piezomagnetic cubic monocrystals of class m3m that can be readily used for non-destructive testing. The piezomagnetic waves propagate in direction [101] corresponding to relatively easy magnetization for the following piezomagnetics: Galfenol, Terfenol-D, and CoFe2O4 with cubic structures. The phase velocities of the USZW-waves and the coefficient of magnetomechanical coupling (CMMC) K were calculated for the crystals. It was found that the coefficient K for piezomagnetics with Km > 1/3 and Km >> 1/3 is about 8% to 9%, where K = 2 (VUSZW,o – VUSZW,c)/VUSZW,o and Km = h/(Cμ). Knowledge of piezomagnetic properties of cubic crystals makes possible the use of them in new products utilizing the phenomenon called the magnetoelectric effect. Also, this study is useful for possible application of cubic piezomagnetics in composite structures consisting of piezoelectric and (or) piezomagnetic materials and in the microwave technology. This broadens choice of possible piezomagnetic materials for utilization in various technical devices.


Introduction
Ferroic materials can display a spontaneous magnetization (ferromagnetic), polarization (ferroelectric) and strain (ferroelastic, shape memory alloy).Materials possessing two "ferro" properties simultaneously are called "biferroics".There are three kinds of biferroic materials: electroelastic, magnetoelastic, and magnetoelectric materials.The magnetoelectric materials (composites) can simultaneously possess ferromagnetic and ferroelectric properties.Composite materials or structures consisting of piezoelectric and piezomagnetic phases are able to facilitate the conversion of energy between electric and magnetic fields.Such phenomenon is called magnetoelectric (ME) effect.The ME effect of the composites was first reported by van Suchtelen [1].The ME effect was then studied by van den Boomgaard [2], van Run et al. [3], and van den Boomgaard et al. [4] for BaTiO 3 / CoFe 2 O 4 composites.Possible applications of ME materials include magnetic-electric energy converting components; solid state non-volatile memory, multi-state memory, which can find application in quantum computing area; and electrical/optical polarization components, which can find applications in communication, light computing, and solid state memories based on spintronics [5,6].Also, this research topic can be found in a recent review paper [7] by Fiebig.Note that the ME effect was originally predicted by P. Curie [8] in 1894.
Surface acoustic wave (SAW) devices are widely used in numerous branches of science and technology and their investigation especially in the case of interconnected physical fields is an important and actively developing branch of research and applications.About 120 years ago the first type of SAW was described by Lord Rayleigh [9] in connection with the problem of earthquakes.Bleustein [10] and Gulyaev [11] theoretically predicted that a pure shear-horizontal (SH) surface wave can be guided by the free surface of a piezoelectric half-space.Later, Maerfeld and Tournois [12] investigated SH acoustic waves guided by the interface of two half-spaces, and Danicki [13] described the acoustic waves which can be guided by an embedded conducting plane in the electro-elastic materials of class 6 mm.It is thought that the simplest type of SH-SAW represents Love waves [14] in layered systems, consisting of an isotropic layer and an isotropic substrate.Recent developments in physics and technology made possible to construct new magneto-electro-elastic materials which demonstrate interconnection between magnetic, electric, and elastic fields [15,16].When the electric field was connected with the elastic one (piezoelectric materials) it brought up new and unexpected possibilities for science and technology.Connecting the magnetic field with the electric and elastic ones in magneto-electro-elastic materials suggests a range of new possibilities.Alshits et al. [17] conducted a qualitative investigation on the existence of surface waves in half-infinite anisotropic elastic media with piezoelectric, piezomagnetic, and ME effects.Soh and Liu [18] gave the existence conditions of interfacial SH-waves in a piezoelectric/piezomagnetic bi-material.Although much attention has been recently concentrated on magneto-electro-elastic materials and several dynamic problems have been solved by Hu and Li [19], Li [20], Chen, et al., [21], etc., the investigation of SAW propagation in transversely isotropic magneto-electro-elastic materials [22][23][24][25] of class 6 mm is currently an open actual subject.In addition to the transversely isotropic materials there are piezoelectric and piezomagnetic crystals with cubic symmetries, in which the surface BG-waves [10,11] cannot exist, according to [26] by Gulyaev and Hickernell.That is also true for the interfacial MTwaves [12].However, new SH-SAWs called the ultrasonic surface Zakharenko waves (USZWs) were recently discovered in [27], which can propagate in cubic piezoelectrics.That is also true for cubic piezomagnetics.The purpose of this paper is to show the first evidence of the USZW existence in cubic piezomagnetics.It is possible to introduce piezomagnetic cubic crystals which are today highly-called for different technical devices.
The Cobalt Ferrite CoFe 2 O 4 (CoO•Fe 2 O 3 or CFO) has the spinel ferrite structure with a stoichiometry denoted by AB 2 O 4 , where A is a divalent transition metal ion or a mixture of a trivalent and a monovalent metal ion.A spinel (see Figure 1 and [28]) is a cubic structure in which O 2-ions form an fcc-lattice.It requires eight Formula equivalent AB 2 O 4 to form a repeating unit cell which contains 32 O 2-ions, 64 tetragonal sites (A sites) and 32 octahedral sites (B sites).However, only 8 tetragonal and 16 octahedral sites are occupied so that the ratio of A-atoms/B-atoms is 1:2.CFO is known to have highly anisotropic magnetic properties and has been widely used in magnetic media and microwave device applications.
In 1971, the U.S. Navy developed an alloy of rare metals Tb and Dy with iron to create Terfenol-D (Tb x Dy 1-x Fe y ).The name "Terfenol-D" comes from a combination of TER for terbium, FE for iron, NOL for Naval Ordnance Laboratory, and D for dysprosium.Recently an iron/gallium alloy termed Galfenol has also been discovered at the Naval Surface Warfare Center (former NOL), where the magnetostrictive material Terfenol-D was discovered.It was found that magnetostriction in iron/gallium alloys peaks at a volume fraction of~17% gallium.Unlike Terfenol-D representing a giant magnetostrictive material useable in practical operating conditions, FeGa-alloys are tough and not toxic, and can be machined and used without any special handling in devices where Terfenol-D may fracture.The cubic structure of Galfenol is shown in Figure 2. Terfenol-D has been commercially available since the late 1980's and has since grown into an international industry.The Terfenol-D rods can be readily manufactured using the freestand-zone-melt process resulting in a cubic crystal structure [29] featuring a large magnetic anisotropy.Also, FeGa-alloys appear promising in filling the role as a mechanically robust material with substantial magnetostrictive capability [30].Additional benefits of FeGa-alloys include substantially lower material costs than rare-earth based magnetostrictive alternatives.For example, the Terfenol-D (Tb 27 Dy 73 Fe 195 ) raw material cost is approximately $ 0.50/g and current production techniques utilize carefully controlled crystal growth processes [31].For comparison, the raw material cost for Fe 81 Ga 19 is approximately $ 0.08/g [31].Indeed, Terfenol-D can be used in magnetic transducers.It is wellknown that magnetic transducers are increasingly considered as actuators and sensors for numerous aerospace, aeronautic, automotive, industrial, and biomedical applications, for example see [32][33][34].Recent progress in studying the cubic structure piezomagnetic crystals Galfenol and Terfenol-D of class m3m for utilization in transducers and sensors can be additionally found in [35][36][37][38].Section 2 theoretically describes SAW propagation in piezomagnetic cubic crystals.

Theory
The rectangular coordinate system (x 1 , x 2 , x 3 ) is shown in Figure 3.The suitable crystal cut and propagation direction was obtained by 45 0 -rotation around the x 2 -axis directed along the lowest odd-order symmetry axis of a piezomagnetic cubic crystal.Therefore, the piezomagnetic waves with the anti-plane polarization propagate along the crystal surface in direction [101].It is noted that the x 2 -axis is perpendicular to the sagittal plane.Theoretical description of wave propagation in cubic piezomagnetics is similar to that for cubic piezoelectrics, for example see [27,39,40].The constitutive equations for a piezomagnetic material can be expressed in terms of the strains τ and the magnetic field H. Strains are related to mechanical displacements: τ ij = (∂U i /∂x j + ∂U j /∂x i )/2 [41].The governing mechanical equilibrium is ∂σ ij /∂x j = 0 and the governing magnetostatic equilibrium is ∂B i /∂x i = 0 where σ ij and B i are the stress tensor and magnetic flux, respectively.A piezomagnetic medium possesses the elastic C ijkl and piezomagnetic h ijk coefficients as well as the magnetic permeability coefficients μ ij and the medium density .
It is necessary to write governing equations for the linear case.Constitutive relations are written as follows: in which σ ij and τ ij are the stress and strain tensors, respectively; B m and H m are components of the magnetic induction (i.e.magnetic flux) and magnetic field: H m = -∂ψ/∂x m , where ψ is the magnetic potential.The indices i, j, k, l, m, and n run from 1 to 3. According to Voigt's usual notation, C ijkl , h ijm , and μ mn can be written as 6 × 6, 3 × 6, and 3 × 3 matrices standing for the elastic, piezomagnetic, and magnetic tensors, respectively.They are thermodynamically defined as follows: where the elastic constants C ijkl are defined at constant magnetic field and the magnetic constants μ mn are defined at constant strain.
In the quasi-static approximation, the equations of motion of an elastic medium and magnetostatics are In Equation ( 6), U i denote the mechanical displacement components; t is time.Using Equations ( 1), ( 2), (6) and ψ = U 4 , one can write the coupled equations of motion for a piezomagnetic medium in the following form: Solutions of homogeneous partial differential Equations (7) of the second order are found in the following plane wave form: where the index i runs from 1 to 3. U i 0 and ψ are initial amplitudes; j = (-1) 1/2 .k s r denotes the scalar multiplication of two vectors and ω is the angular frequency.{k 1 , k 2 , k 3 } = k{n 1 , n 2 , n 3 } are the components of the wavevector k s , {x 1 , x 2 , x 3 } are the components of the real space vector r, and {n 1 , n 2 , n 3 } are the directional cosines.
The coupled equations of motion can be readily written in the following simplified form, leaving only equations for waves with polarization perpendicular to the sagittal plane as well as non-zero components of the material ten- sors for the studied direction [101] of wave propagation: In Equation ( 8), the mechanical displacement component U 2 is directed along the x 2 -axis, see Figure 3: where the phase velocity is defined as  8).The dependence of the normalized piezomagnetic constants on the cubic crystal cuts and propagation directions is shown in Figure 4.
Substituting the mechanical displacement U 2 and magnetic potential ψ = U 4 of Equation ( 9) into Equation ( 8), the equations of motion can readily be written in the wellknown tensor form, using corresponding GL-components in the Green-Christoffel equation, (GL rw -δ rw ρV ph ) U r = 0 [39,41]:   and 2 44 3 with n 3 = k 3 /k.In the GL-equation, r and w run from 1 to 4, δ rw is the Kronecker delta for r < 4 and w < 4, δ 44 = 0, U r = {U 1 , U 2 , U 3 , ψ}.Therefore, the following system of two homogeneous equations for pure SH-waves can be written as: In Equation (10), the directional cosines are defined as follows: n 1 ≡ 1, n 2 ≡ 0 and n 3 = n 3 .The velocity V t4 is also defined as follows: V t4 = (C 44 /ρ) 1/2 .Setting the matrix determinant equal to zero in Equation ( 10), the suitable phase velocity V ph satisfying boundary conditions discussed in the following section and four polynomial roots n 3 (p) (V ph ), as well as the functions U 2 0 (V ph ) and φ 0 (V ph ), can then be found.For example, the functions can be taken in the following form: 0 42 GL   and . It is noted that for the piezomagnetic surface Bleustein-Gulyaev waves, the GL 24 and GL 42 components are as follows: Substituting m = 1 + n 3 2 and expanding the matrix determinant in Equation ( 10), the following polynomial can be obtained from Equation (10): ) of which two roots can be found as follows giving four polynomial roots of Equation ( 10) The polynomial roots n 3 (1,2,3,4) actually represent the eigenvalues.For each eigenvalue n 3 (n) , the corresponding eigenvector can be found in the following form: (U 2 (n) , ψ (n) ) where the index n runs from 1 to 4. The coefficient K m 2 in Equations ( 11) and ( 12) is the static coefficient of the magnetomechanical coupling (CMMC): It is also noted that the speed V tm of the bulk SH-wave is given by the following Formula: Analyzing the roots for propagation direction [101] in Equations ( 12) and ( 13), it can be found that all complex roots will be calculated when the expression under the square root in Equation ( 12) is negative.That fulfills for velocities V ph being lower than some velocity V Km obtained solving the following Equation from (12): and defined by the following formula: It is clearly seen in Equation ( 17) that the factor A Km Degree and direction Piezomagnetic constants, a.u.15) is also shown in Figure 5.The function ) 1/2 giving the following equality: V Km = V tm (see Figure 5).It is noted that only complex polynomial roots can exist when V ph < V Km .The CMMC K 0 2 = 1/3 is readily found by substituting the velocity V tm from Equation (15) instead of the phase velocity V ph in Equation (16).Note that when K m 2 < K 0 2 there are all imaginary roots for V ph > V Km , but a big K m 2 > K 0 2 gives real roots when V ph > V Km .It is also noted that only complex or imaginary roots with negative imaginary parts are chosen in order to cope with wave damping towards the depth of a crystal corresponding to negative values of the x 3 -axis shown in Figure 3.
Table 1 lists experimental data for Galfenol (83% of Fe and 17% of Ga) obtained in [42,43].It is clearly seen that the static CMMC K m 2 for cubic Galfenol can be as high as~0.42that is larger than 1/3, see also Table 2.The Galfenol density ρ was taken from [31], which can depend on material synthesis.For example, [31] gives a Galfenol density of~5900 kg/m 3 for a typical sintered disk that is only 74% of maximum possible density.It is noted that the material density of Fe is~7848 kg/m 3 (see [44]).According to [31,42] as well as many Galfenol Workshops at the University of Maryland, the shear elastic constant c 44 = c for Galfenol relatively slightly depends on both temperature and concentration of ga in fe.The elastic constant c can have values from~ 120 GPa to~135 gPa.According to [31], the static CMMC K m 2 for cubic Terfenol-D can be usually about two times larger than that for Galfenol under the same conditions.However, it is interesting to compare some characteristics of Galfenol and Terfenol-D.Therefore, Terfenol-D with a very big K m 2 ~1.2 was chosen from [45] for numerical calculations suggesting that Terfenol-d has the cubic symmetry.The piezomagnetics CoFe 2 O 4 [16,46] is also listed in Table 2 for comparison, suggesting the cubic symmetry for the crystal.
Since the experimental data listed in table 1 is given for different experimental conditions, it is possible to calculate the single piezomagnetic constant h = h 14 for cubic Galfenol, using known values of K m and C, as well as the measured values of magnetic permeability μ given in Table 1.The strain τ and magnetic flux B can be also written based on the following linear constitutive piezomagnetic equations: where s H is the compliance at constant magnetic field H and μ σ is the magnetic permeability at a constant stress σ.Note that Equation ( 19) incorporates the direct magnetostrictive effect in which magnetic flux is generated by stresses in the piezomagnetic material.The piezomagnetic constants q and q * are thermodynamically defined as follows, assuming q = q * : The piezomagnetic coefficient q symbolizes the axial magnetostriction coefficient and represents the change in strain per change in magnetic field at a constant stress.The parameter q* represents the change in magnetic induction per change in stress at a constant magnetic field.
The short circuit elastic modulus E B represents the stiffest-material condition, which occurs when all available magnetic energy has been transduced into elastic potential energy.When energy is transferred from the elastic to the magnetic regime the effective modulus decreases to the value of E H .It has been shown by invoking linear transduction considerations that these two elastic moduli are related by the magnetomechanical coupling coefficient K m 2 (0 ≤ K m 2 ≤ 1) as follows [47]: In a similar manner, the intrinsic or uncoupled magnetic permeability μ σ of Equation ( 19) can be reduced to a value corresponding to the constant strain permeability μ ε because of the energy conversion from the magnetic to the elastic regime.The permeabilities are related by It is emphasized that expressions ( 21) and ( 22) are only valid for linear systems and significant corrections are necessary in order to extend these formulae to the full nonlinear regime.However, they are highly useful to illustrate the exchange mechanisms taking place in magnetostrictive transduction.
The static CMMC K m 2 can be also written as follows: Therefore, the single non-zero piezomagnetic constant h can be calculated with the following formula, using Equations ( 14), (22), and (23): Table 2 lists the piezomagnetic constant h for cubic Galfenol corresponding to the experimentally measured values of the magnetic permeability μ σ at constant stress σ.
For the free space, Laplace's equation of type ψ = 0 is written in the following form: (k 1 2 + k 3 2 ) ψ 0 = 0 where μ 0 is the magnetic constant for the free space.The magnetic potential for the free space can be written as follows: ψ 0 = F (0) exp (-k 1 x 3 ) exp [j (k 1 x 1 -t)].The potential ψ 0 must decrease according to any increase in the coordinate x 3 (see Figure 3).It is now possible to use the material constants in mechanical and magnetic boundary conditions.

Boundary Conditions for Surface SH-waves
The mechanical boundary conditions of a magnetoelastic material which occupies the half-space x 3 < 0 must be satisfied together with the magnetic conditions of a magnetically closed (B 3 = 0) or magnetically open (ψ = 0) surface.The realization of each of these boundary conditions is described in [17].There is the single mechanical boundary condition for the normal component of the stress tensor: σ 32 = 0 at x 3 = 0, where Also, there are the magnetic boundary conditions: continuity of the normal component B 3 of the magnetic flux at x 3 = 0 being the interface between a vacuum (B 3 f ) and the crystal surface (B 3 = B 3 f ) where and continuity of the magnetic potential U 4 = ψ at x 3 = 0 (ψ = ψ f ) where Therefore, using Equations ( 25)-( 27), two homogeneous equations for the case of a magnetically closed surface can be readily written as follows: For a magnetically open surface, the following homogeneous equations can also be written in the corresponding form: (2) Therefore, the corresponding boundary-condition determinants (BCDs) of matrices in Equations ( 28) and (29) for both cases represent complex numbers.
The complete mechanical displacement U 2 Σ and magnetic potential ψ Σ = U 4 Σ can be written in the plane wave form as follows: The corresponding weight functions F (1) and F (2) are found from Equations ( 28) and (29), which can give the same eigenvectors (U 2 0(1) , ψ 0 (1) ) and (U 2 0(2) , ψ 0(2) ) for two equal eigenvalues n 3 (1) = n 3 (2) , and hence (1)   (2) . It is obvious that for this case the weight factors (1)  F  (2) will zero the complete mechanical displacement U 2 Σ and magnetic potential ψ Σ in Equations ( 28) and ( 29) giving "latent" characteristics in Equation (30).On the other hand, unequal eigenvalues n 3 (1) and n 3 (2) give different eigenvectors (U 2 0(1) , ψ 0(1) ) and (U 2 0(2) , ψ 0(2) ).It is thought that the V ph -solutions with the "latent" characteristics must be experimentally verified due to the fact that some elements of crystals symmetry (screw axis or glide reflection) can be broken near the surface.The same relates to simple reflections and axis if the surface restricts the crystal in such a way that it breaks some simple elements of the crystal class group.It is also noted that V phsolutions cannot be found when the surface Bleustein-Gulyaev waves are studied that is an additional difference for finding the BG-waves and USZW-waves.

Results and Discussions
Acoustic wave propagation along the surface of piezoelectric material in direction [101] as an instability problem can sustain the surface SH-waves in cubic piezoelectrics that was recently discovered in [27].Concerning piezomagnetics, it is thought that piezomagnetic SH-SAW can also propagate in direction [101] in cubic piezomagnetics, for instance, crystals of class m3m.In-deed, all cubic crystals of the class possess no piezoelectricity according to [48,49].Note that concerning transversely-isotropic materials, solutions for the surface BG-waves in piezomagnetics can be written in an identical manner to the wave solutions in piezoelectrics [50].It is also noted that the velocity V BGmo of piezomagnetic surface BG-wave propagating along the magnetically open surface of transversely-isotropic piezomagnetics (treating the simplest case) can be found with the following well-known formula: The velocity V BGmc of piezomagnetic BG-wave with the magnetically closed surface can be calculated with the following formula: Note that Formula (31) for the velocity V BGmo is not applicable for the case of cubic piezomagnetics.For example, formula (31) gives V BGmo ~4760.093893m/s for Galfenol with the largest value of K m 2 > 1/3 in the first row in Table 3.This value of V BGmo is close to that of V tm , but not to the value of V Km (see Table 3) representing the right result for cubic piezomagnetics with K m 2 > 1/3.The same there is for Terfenol-D with K m 2 ~1 in Table 2, for which the incorrect velocity calculated with Formula (31) is V BGmo ~1179.199962m/s.
For cubic piezomagnetics with K m 2 < 1/3, Formula (31) also gives values of V BGmo being closer to corresponding values of V tm than true values of V USZW,o .For example, the velocity V USZW,o for Galfenol in the fourth row in Table 3 with the smallest value of K m 2 < 1/3 (see also Table 2) is V USZW,o ~4328.763803m/s.However, Formula (31) Table 3.The velocities v Km , V t4 , and V tm (all in m/s) as well as the velocities V USZW,o and V USZW,c of the ultrasonic surface Zakharenko waves for [101] propagation direction in the piezomagnetic cubic crystals.The last column gives values of the K 2 calculated with Formula (33).All the Galfenols correspond to those listed in Tables 1 and 2 Material (symmetry group)  3 is already very close to the velocity V BGmo ~2985.08072m/s, which is also very close to the velocity V tm .Note that the relative magnetic permeability μ/μ 0 in Formula (31) for piezomagnetics is about one order larger than the corresponding relative dielectric constant ε/ε 0 for piezoelectrics for the same K m 2 = K e 2 , where K e 2 is the static coefficient of electromechanical coupling (CEMC), and ε and ε 0 are the single electric constants for a piezoelectrics and free space, respectively.That can result in mm/s-difference for piezomagnetics compared with cm/s-difference for piezoelectrics between the corresponding BG-wave velocity and bulk wave velocity V t .That can result in significantly larger penetration depth of surface BG and USZW-waves for piezomagnetics towards negative values of the x 3 -axis in Figure 3.It is also noted that the velocity V BGmc for a magnetically closed surface can be readily calculated with Formula (32), with which the velocity of interfacial Maerfeld-Tournois waves [12] is also calculated.The interfacial Maerfeld-Tournois waves can propagate along the interface of two transversely-isotropic half-spaces of class 6 mm.Concerning cubic piezomagnetics, it is expected that some interfacial waves can also propagate that can be further researched.
Using the velocities V USZW,o and V USZW,c for the USZWwaves, it is possible to evaluate the coefficient of magnetomechanical coupling (CMMC) K 2 for the piezomagnetic cubic crystals listed in Table 3.Indeed, it is thought that it is possible to use the well-known formula that is used for piezoelectrics.Therefore, the coefficient K 2 can be evaluated with the following formula: For Galfenol in Tables 2 and 3, the evaluated value of K 2 is as high as 2.26% for the smallest value of K m 2 .The value of K 2 increases with increase in K m 2 to the biggest value of 8.25%.This value of 8.25% for the piezomagnetic Galfenol in the first row of table 3 can be compared with that for Terfenol-D in the Table: K 2 (Terfenol-D)~ 9.13%.It is clearly seen in the second and third tables that a big value of K m 2 ~1.2 for Terfenol-D being about three times larger than the biggest value of K m 2 for Galfenol does not give a significant rise in the coefficient K 2 .The rise is smaller than 1%.That can be explained using the following fact: the velocity V USZW,o for cubic piezomagnetics with a giant coefficient K m 2 is situated slightly below the speed V Km , but not slightly below the speed V tm (see also Figure 5).That results in lower coef-ficient K 2 .Therefore, it can be concluded that the coefficient K 2 for cubic piezomagnetics with a giant K m 2 is about 8% to 9% and cannot be larger than 10%.It is thought that the coefficient K 2 for the surface BG-waves in transversely-isotropic piezomagnetics with a giant K m 2 can be larger than 10%, according to Formulae (31)- (33) and Figure 8 shows the dependence of the V ph on the coefficient K m for the Fe-17Ga alloy.Note that any change in the material density will change both the velocities V t4 and V tm (V t4 ), but not the coefficients K m 2 and K 2 .

Conclusions
This paper addresses the theoretical investigations of piezomagnetic acoustic waves propagating along the surface of cubic crystals.The ultrasonic surface Zakharenko waves (USZWs) can propagate in cubic piezomagnetics when direction , where K m 2 is the static coefficient of magnetomechanical coupling.The simulations have shown that for a cubic crystal with a large K m 2 > 1/3 (even K m 2 > 1 for Terfenol-D) the coefficient K 2 is as big as 8%-9% and not larger than 10%.It is noted that a big K m 2 > 1 does not give a significant rise to the coefficient K 2 for the case of K m 2 ~1/3.That is due to the fact that velocity solutions for large values of K m 2 > 1/3 are found just below the velocity V Km , but not just below the velocity V tm .This theoretical study of cubic piezomagnetics allows use of suitable surfaces of piezomagnetic cubic crystals in non-destructive testing and evaluation that was also recently done for cubic piezoelectrics.

Figure 1 .Figure 2 .
Figure 1.The normal spinel structure of Cobalt Ferrite CoFe 2 O 4 (CoO•Fe 2 O 3 ) with a stoichiometry denoted by AB 2 O 4 .Only the lower half of the unit cell is shown for projection on (001)

Figure 4 .
Figure 4.The dependence of the normalized piezomagnetic constants h 14 = h 36 , h 16 , and h 34 on the propagation directions with the Euler angles {0 o , θ, 0 o } for piezomagnetic cubic crystals, where there is rotation around the x 2 -axis

Figure 5 . 2 ) 1 / 2 2 ) for small values of the K m 2 Table 1 .
Figure 5.The dependence of both the factor A Km and the function (1 + K m 2 ) 1/2 on the values of K m 2 from 0 up to 3. The insertion shows the A Km (K m 2 ) for small values of the K m 2

Figure 5 .
The sample behaviors of the boundarycondition determinant (BCD3) around the phase velocity solutions for the piezomagnetic monocrystals CoFe 2 O 4 and Galfenol are shown in Figures6 and 7, respectively.

Figure 6 . 4 PhaseFigure 7 .VelocityFigure 8 .
Figure 6.The behavior of the boundary-condition determinant (BCD3) for the piezomagnetic monocrystal CoFe 2 O 4 is the wavenumber in direction of wave propagation).For piezomagnetic cubic crystals, it is possible to write C 44 = C 66 = C and μ 11 = μ 33 = μ.It is noted that many propaga- tion directions perpendicular to direction [010] can exist.It is also possible to cut a cubic crystal in order to study wave propagation in direction [101] with h 14 = h 36 = 0 and the non-zero piezomagnetic constants {-h 16 , h 34 } that is shown in Equation (

Table 2 ,
% BGmo ~4328.783915m/s that is closer to V tm ~ 4328.798060m/s.Therefore, the differences Δ BG = V BGmo -V tm and Δ USZW = V USZW,o -V tm are~0.014m/s and~0.034m/s, respectively.For comparison, for the piezomagnetic CoFe 2 O 4 with the smallest value of K m 2 ~4% from the true velocity V USZW,o in Table