**Radiation and Propagation of Waves in Magnetic Materials with Helicoidal Magnetic Structure**

Igor V. Bychkov, Dmitry A. Kuzmin and Vladimir G. Shavrov

Additional information is available at the end of the chapter

http://dx.doi.org/10.5772/64014

#### **Abstract**

In this chapter, we are shortly reviewing some problems of electromagnetic and acoustic wave propagation and radiation in the magnets with helicoidal spin structure. We show the band structure of the coupled wave spectrum in the materials. The band gap width depends on the spiral angle (or, equivalently, on external magnetic field value). Interaction of spin and electromagnetic waves leads to opening the gap in spinelectromagnetic dispersion. This gap leads to opacity window in reflection spectrum of spiral magnet plate. The opacity window closes at phase transition into collinear ferromagnetic state and reaches a maximum at simple spiral state. At the frequencies near band gap boundaries, the rotation of polarization plane of propagating electro‐ magnetic wave is observed. Account of interaction of spin and electromagnetic waves with acoustic subsystem leads to opening the gap in spin-acoustic spectrum. This gap leads to some features in electromagnetic reflectance spectrum and to rotation of acoustic wave polarization plane, i.e. to acoustic Faraday effect. We also show the possibility of acoustic and electromagnetic wave radiation by helicoidal magnets at phase transition into collinear ferromagnetic state. Some features of electromagnetic waves generation by spiral magnets placed in homogeneous magnetic field with harmonical time-dependence are also discussed.

**Keywords:** spiral magnets, helicoidal magnetic structure, electromagnetic waves, spin waves, acoustic waves, Faraday effect, acoustic Faraday effect

© 2016 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## **1. Introduction**

The magnets with spiral (helicoidal, or modulated) magnetic structure have been under investigation for a long time [1, 2]. These materials include the rare-earth metals, compounds based on them, and some of the compounds based on transition metals. Such materials may be semiconductors, as well as dielectrics. The reasons of such ordering may be the following: competing exchange interaction of atoms in the first and the second coordination spheres, asymmetric exchange, relativistic interaction (Dzyaloshinskii-Moriya interaction). The main characteristic of modulated structures is a modulation wave number *q* = 2π/*L*, where *L* is the modulation period. The modulation period of such structures is much greater than the lattice constant of the crystal. Modulation period of helicoidal magnets usually changes continuous‐ lywiththechangeoftemperature.Itmaybeincommensurablewiththecrystal'slatticeconstants. This is the reason why modulated, or long-periodic, magnetic structures are frequently called "incommensurable" ones.

Modulated magnetic structure leads to sufficient features of dynamical properties of helicoidal magnets comparing with the properties of collinear ferro- and anti-ferromagnets. Thus, the spectrum of spin waves has a band structure. It has no activation frequency when wave number of propagating wave is equal to the modulation wave number [2], i.e. *k* = *q*, while in collinear ferro-and anti-ferromagnets this situation is observed at *k* = 0. Magnetoelastic waves in ferromagnets with spiral structure have been studied in details, for example, in works [3–8]. Works [4–6] are devoted to investigation of spin and acoustic oscillations interaction in spiral magnets without taking into account the spontaneous deformations in the ground state. These deformations were consistently considered in work [8]. Interaction of spin and electromagnetic waves in magnetic dielectrics with collinear magnetic structure has been investigated in [9]. Peculiarities of spin waves and helicons hybridization have been considered in [10]. The result of this work has been generalized for the case of arbitrary propagation direction in [11]. Later, coupled spin-Alfven waves in ferromagnetic metals [12], and coupled plasmon-spin-electro‐ magnetic waves in ferromagnetic semiconductors and metals with different anisotropy types [13] have been studied. These investigations show that interaction of the spin subsystem with the electromagnetic waves usually leads to the shift of activation frequency of quasi-spin branch: an additional term of electromagnetic nature, namely magnetostatic frequency, appears. This interaction may also lead to reduce the phase speed of electromagnetic wave.

Modulated magnetic structure may manifest itself in the modulation of electromagnetic characteristics of the magnet. Phenomenological approach has been successfully applied to investigate the magneto optical effects in the medium with the simple spiral structure [14]. Some features of light propagation have been predicted in works [15, 16]. Authors of these works have not taken into account the dynamics of magnetic subsystem. Such approach is valid for optical frequencies, when magnetic subsystem does not keep pace with changes of field's magnetic component. As for investigation of electrodynamical processes at lower frequencies, taking into account dynamics of magnetic subsystem may lead to qualitatively new effects.

Coupled spin-electromagnetic and spin-acoustic-electromagnetic waves have been investi‐ gated in detail for magnets with the simple spiral structure in [17, 18]. Wave hybridization in such magnets leads to opening the gaps in spectrums of electromagnetic waves, as well as in acoustic ones. Due to the band structure of the spectrum, some features in electromagnetic and acoustic reflectance may be observed. All these features should have a place in magnets with conical spiral ordering as well.

**1. Introduction**

162 Magnetic Materials

"incommensurable" ones.

new effects.

The magnets with spiral (helicoidal, or modulated) magnetic structure have been under investigation for a long time [1, 2]. These materials include the rare-earth metals, compounds based on them, and some of the compounds based on transition metals. Such materials may be semiconductors, as well as dielectrics. The reasons of such ordering may be the following: competing exchange interaction of atoms in the first and the second coordination spheres, asymmetric exchange, relativistic interaction (Dzyaloshinskii-Moriya interaction). The main characteristic of modulated structures is a modulation wave number *q* = 2π/*L*, where *L* is the modulation period. The modulation period of such structures is much greater than the lattice constant of the crystal. Modulation period of helicoidal magnets usually changes continuous‐ lywiththechangeoftemperature.Itmaybeincommensurablewiththecrystal'slatticeconstants. This is the reason why modulated, or long-periodic, magnetic structures are frequently called

Modulated magnetic structure leads to sufficient features of dynamical properties of helicoidal magnets comparing with the properties of collinear ferro- and anti-ferromagnets. Thus, the spectrum of spin waves has a band structure. It has no activation frequency when wave number of propagating wave is equal to the modulation wave number [2], i.e. *k* = *q*, while in collinear ferro-and anti-ferromagnets this situation is observed at *k* = 0. Magnetoelastic waves in ferromagnets with spiral structure have been studied in details, for example, in works [3–8]. Works [4–6] are devoted to investigation of spin and acoustic oscillations interaction in spiral magnets without taking into account the spontaneous deformations in the ground state. These deformations were consistently considered in work [8]. Interaction of spin and electromagnetic waves in magnetic dielectrics with collinear magnetic structure has been investigated in [9]. Peculiarities of spin waves and helicons hybridization have been considered in [10]. The result of this work has been generalized for the case of arbitrary propagation direction in [11]. Later, coupled spin-Alfven waves in ferromagnetic metals [12], and coupled plasmon-spin-electro‐ magnetic waves in ferromagnetic semiconductors and metals with different anisotropy types [13] have been studied. These investigations show that interaction of the spin subsystem with the electromagnetic waves usually leads to the shift of activation frequency of quasi-spin branch: an additional term of electromagnetic nature, namely magnetostatic frequency, appears. This interaction may also lead to reduce the phase speed of electromagnetic wave.

Modulated magnetic structure may manifest itself in the modulation of electromagnetic characteristics of the magnet. Phenomenological approach has been successfully applied to investigate the magneto optical effects in the medium with the simple spiral structure [14]. Some features of light propagation have been predicted in works [15, 16]. Authors of these works have not taken into account the dynamics of magnetic subsystem. Such approach is valid for optical frequencies, when magnetic subsystem does not keep pace with changes of field's magnetic component. As for investigation of electrodynamical processes at lower frequencies, taking into account dynamics of magnetic subsystem may lead to qualitatively A special kind of dynamic effects in materials with a modulated magnetic structure can occur during phase transitions, accompanied by appearance/disappearance of the modulated structure. Similar effects are well known for periodically distributed sources of different nature and they usually consist of different wave radiation. Thus, nanostructured metallic films may be sources of radiation when they are irradiated by laser pulses [19, 20]. Theoretical model for describing this effect has been proposed in [21]. Terahertz wave radiation may occur when magnons are excited by femtosecond laser pulses in anti-ferromagnets [22], or by infrared pulses in diluted magnetic semiconductors [23]. In helicoidal magnets, the spin structure is reorganized during phase transition. This may also lead to radiation of electromagnetic waves and sound. Recently [24], we have observed experimentally electromagnetic waves radiation by Heusler alloy during magneto-structural phase transition, possible mechanisms and theory of this effect have been discussed as well. We have supposed that the main contribution is caused by moving dislocations with dipole moment. Here, we show that in spiral magnets during orientational magnetic phase transition, electromagnetic radiation may be comparable to and even be more than in Heusler alloys.

Nowadays, scientific interest on the investigation of helicoidal magnets is caused by the fact that in such magnets, spontaneous polarization may appear frequently due to specific magnetic ordering. Such materials, combining both electric and magnetic properties, are called multiferroics. First, in 2003 such effect was observed in TbMnO3 [25], then in 2004 in TbMn2O5 [26], and many other materials. Experiments clearly show the coupling between spontaneous polarization and modulated magnetic ordering. Phenomenologically, coupling of polarization and non-uniform magnetic ordering may be described by energy term of third order *PM*∂*M*, which is allowed by symmetry. This coupling is linear according to polarization. Thus, even weak interaction of such type will lead to electric polarization at corresponding magnetic arrangement. It this chapter we will not consider multiferroic materials, but quali‐ tatively all results are valid even in case of multiferroics. Details of peculiarities of wave interaction in helicoidal multiferroics may be found, for example, in our previous works [27– 29].

This chapter consists of introduction, two sections and conclusions. In the first section, we review some peculiarities of electromagnetic and acoustic wave propagation in helical magnets. The results cover the coupled wave spectrum, electromagnetic reflection from halfspace and plate of magnet with the spiral magnetic structure, Faraday rotation, and acoustic Faraday rotation (i.e. rotation of polarization plane of acoustic wave). In the second section, we consider the possibility of electromagnetic and acoustic wave radiation in spiral magnets during the phase transition. An infinitely fast phase transition and the linear change of the transverse magnetization component are considered as model processes.

## **2. Influence of electromagnetic-spin-acoustic interaction on wave propagation**

This section is mainly based on results of our previous investigations [30–33].

Let us consider the magnet with spiral magnetisation distribution at the ground state. Such magnetic order may be realized if free energy of the crystal has the following form:

$$F = \frac{\alpha}{2} \left(\frac{\partial M\_{\iota}}{\partial \mathbf{x}\_{\iota}}\right)^{2} + \frac{\mathfrak{B}\_{1}}{2} M\_{z}^{2} + \frac{\mathfrak{B}\_{2}}{2} M\_{z}^{4} - M\_{z} H + b\_{ij\text{in}} M\_{i} M\_{j} u\_{\text{lm}} + c\_{ij\text{in}} u\_{\text{ij}} u\_{\text{lm}} + F\_{\text{in}},\tag{1}$$

where **M** = (*Mx, My, Mz*) is the magnetization of the crystal, α is the exchange constant, β1 and β<sup>2</sup> are anisotropy constants, *bijlm* are magnetostriction constants, *cijlk* are elasticity constants, *i*, *j*, *l*, *k* = *x*, *y*, *z*. The sum is assumed where indexes are repeated. The term *Fin* corresponds to inhomogeneous magnetization distribution at the ground state. It may have an exchange nature (such structures we will call "exchange spiral"): *Fin* <sup>=</sup>*γ*(∂<sup>2</sup>*Mi* /∂ *xi* 2 )/2. Another reason may be relativistic interaction (such situation we will call "relativistic spiral"): *Fin* = α1**M***rot***M**. In (1) we suppose that external magnetic field is directed along *z*-axis, which corresponds to anisotropy axis.

For concreteness, we will assume that crystal with exchange spiral structure is from hexagonal crystal system (*D*3*<sup>d</sup>*, *C*3*<sup>v</sup>*,*D*6, *D*6*<sup>h</sup>*). Magnets with relativistic spiral structure will be assumed as isotropic ones. Their elastic properties may be defined by two characteristics (Lamé constants λ and μ), and magnetostriction will be characterized by one "effective" constant *b*.

Spiral magnetization distribution may be described by the following coordinate dependence of magnetization components:

$$M\_x = M\_0 \sin\theta \cos qz,\newline M\_y = M\_0 \sin\theta \sin qz,\newline M\_z = M\_0 \cos\theta.\newline \tag{2}$$

In (2), *q* = 2π/*L* is the spiral wave number, *L* is the spiral period, θ is the spiral angle, *M*0 is the magnetization amplitude. In case of θ = 0, one will have a ferromagnetic ground state, while θ = π/2 corresponds to simple spiral structure. When 0 < θ < π/2, the "ferromagnetic spiral" phase is realized.

Minimizing the free energy (1) with appropriate expression for *Fin*, taking into account ground state magnetization distribution of spiral type (2), we will get the conditions of the phase stability. For crystals with exchange spiral, these conditions are the following: γ > 0, α < 0, *q* = (−α/2γ)1/2. For relativistic spiral, one will have α1 ≠ 0, α >0, *q* = α1/α. In both cases an external magnetic field is coupled with the spiral angle:

$$M\_0 \cos \Theta \Big[ \tilde{\mathfrak{B}}\_1 + h\_{mc} + \left( \tilde{\mathfrak{B}}\_2 - h\_{mc} \Big/ M\_0^2 \right) M\_0^2 \cos^2 \Theta + \alpha q^2 + \tilde{\Delta} \Big] + H = 0, 1$$

**2. Influence of electromagnetic-spin-acoustic interaction on wave**

magnetic order may be realized if free energy of the crystal has the following form:

*F M M M H b MMu c uu F*

nature (such structures we will call "exchange spiral"): *Fin* <sup>=</sup>*γ*(∂<sup>2</sup>*Mi*

<sup>a</sup> æ ö ¶ b b = + + -+ + + ç ÷ ¶ è ø

Let us consider the magnet with spiral magnetisation distribution at the ground state. Such

1 2 2 4 , 2 22

where **M** = (*Mx, My, Mz*) is the magnetization of the crystal, α is the exchange constant, β1 and β<sup>2</sup> are anisotropy constants, *bijlm* are magnetostriction constants, *cijlk* are elasticity constants, *i*, *j*, *l*, *k* = *x*, *y*, *z*. The sum is assumed where indexes are repeated. The term *Fin* corresponds to inhomogeneous magnetization distribution at the ground state. It may have an exchange

be relativistic interaction (such situation we will call "relativistic spiral"): *Fin* = α1**M***rot***M**. In (1) we suppose that external magnetic field is directed along *z*-axis, which corresponds to

For concreteness, we will assume that crystal with exchange spiral structure is from hexagonal crystal system (*D*3*<sup>d</sup>*, *C*3*<sup>v</sup>*,*D*6, *D*6*<sup>h</sup>*). Magnets with relativistic spiral structure will be assumed as isotropic ones. Their elastic properties may be defined by two characteristics (Lamé constants

Spiral magnetization distribution may be described by the following coordinate dependence

In (2), *q* = 2π/*L* is the spiral wave number, *L* is the spiral period, θ is the spiral angle, *M*0 is the magnetization amplitude. In case of θ = 0, one will have a ferromagnetic ground state, while θ = π/2 corresponds to simple spiral structure. When 0 < θ < π/2, the "ferromagnetic spiral"

Minimizing the free energy (1) with appropriate expression for *Fin*, taking into account ground state magnetization distribution of spiral type (2), we will get the conditions of the phase stability. For crystals with exchange spiral, these conditions are the following: γ > 0, α < 0, *q* = (−α/2γ)1/2. For relativistic spiral, one will have α1 ≠ 0, α >0, *q* = α1/α. In both cases an external

0 00 sin cos , sin sin , cos . *M M qz M M qz M M x yz* =q =q = q (2)

λ and μ), and magnetostriction will be characterized by one "effective" constant *b*.

*z z z ijlm i j lm ijlm ij lm in*

/∂ *xi* 2

)/2. Another reason may

(1)

This section is mainly based on results of our previous investigations [30–33].

**propagation**

164 Magnetic Materials

anisotropy axis.

phase is realized.

of magnetization components:

magnetic field is coupled with the spiral angle:

2

*i*

*M*

*x*

*i*

where β˜ <sup>1</sup> and β˜ 2 are anisotropy constants renormalized by magnetostriction [28], and *h me* =(*b*<sup>11</sup> −*b*12)<sup>2</sup>*M*<sup>0</sup> <sup>2</sup> / (*c*<sup>11</sup> <sup>−</sup>*c*12), *Δ*˜ <sup>=</sup>γ*<sup>q</sup>* <sup>4</sup> for magnets with exchange spiral structure, while *<sup>h</sup> me* <sup>=</sup>*<sup>b</sup>* <sup>2</sup> *M*<sup>0</sup> <sup>2</sup> / <sup>2</sup>*μ*, *Δ*˜ <sup>=</sup> <sup>−</sup>2*α*1*q* for relativistic ones.

There is a critical magnetic field value *Hcr*, when the phase transition into collinear ferromag‐ netic state occurs. Diagonal components of the equilibrium stress tensor are homogeneous and depend on the spiral angle, while the components *uxz* and *uyz* depend on both spiral angle and *z*-coordinate, *uxy* = 0.

To obtain the spectrum of coupled spin, acoustic and electromagnetic waves, one should use the Landau-Lifshitz and Maxwell equations with equation of motion of the elastic medium. Simplifying such system of equations by small oscillations method, linearizing and taking the circular components, one will obtain the following:

( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 2 2 4 4 2 2 1 44 0 2 2 2 0 44 33 31 0 0 0 1 1 cos sin sin cos 2 2 2 1 3 sin ( ) sin 2 2 sin 2 2 sin 2 2 sin cos , *k me me kq z z z m k m k mk q m k q igb M k u k <sup>i</sup> gM b k q u k q ig b b M k q u k q gM h k q gM h k* ± ± ± ± ± ± ± ± é ù ±w = q w + w q + w q q - ê ú ë û é ù -w q + -q - ê ú ë û - q± - - q± + +q - q m m m m m ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 2 2 44 0 0 <sup>1</sup> sin 2 sin 2 4 <sup>1</sup> sin , <sup>2</sup> *m k <sup>z</sup> k qmkq mkq k q <sup>i</sup> gb M k q u k q k q u k q gM h k q h k q* - + - - + + - + + w = q w - -w + - é ù ë û - q - --+ + + é ù ë û + q +- - é ù ë û ( ) ( ) ( ) 2 22 44 0 sin cos , *t z <sup>i</sup> s k u k kb M m k q m k* ± ± é ù w- = é ù q +q ë û ë û <sup>r</sup> m ( ) ( ) ( ) ( ) ( ) ( ) ( ) 2 22 33 31 0 2 22 2 2 cos , 4, 4. *l z z z z s k u k i b b kM m k kv h k m k h k m k* ± ± é ù w - =- - q r ë û é ù w - = -w p =- p ë û

This system of equations should be complemented by the condition of magnetization vector permanence |**M**| =*const*, which for Fourier components has a form:

$$
\sin\Theta \left[ m\_- \left( k - q \right) + m\_+ \left( k + q \right) \right] + 2 \cos\Theta m\_z \left( k \right) = 0.
$$

In all equations above we have used the following notations: **h**(*k*), **m**(*k*), **u**(*k*) are Fourier components of corresponding vector **h**, **m**, **u**; *v* =*c* / *ε* is electromagnetic wave speed in the medium; ε is permittivity of the crystal; *ω*2*<sup>k</sup>* <sup>±</sup> <sup>=</sup>*ω*<sup>20</sup> <sup>+</sup> *gM*0*<sup>L</sup>* <sup>⊥</sup> ±(*k*), *ω*1*<sup>k</sup>* <sup>=</sup>*ω*<sup>10</sup> <sup>+</sup> *gM*0sin<sup>2</sup> *θL* <sup>∥</sup>(*k*)=0, *L* <sup>⊥</sup> ± (*k*)= −*α*(*q* <sup>2</sup> −*k* <sup>2</sup> )−*γ*(*q* <sup>4</sup> −*k* <sup>4</sup> ) <sup>+</sup> <sup>2</sup>*α*1(*<sup>q</sup>* <sup>∓</sup> *<sup>k</sup>*), *<sup>L</sup>* ||(*k*)= <sup>−</sup>*α*(*<sup>q</sup>* <sup>2</sup> <sup>−</sup>*<sup>k</sup>* <sup>2</sup> )−*γ*(*q* <sup>4</sup> −*k* <sup>4</sup> ) + 2*α*1*q*, *st* <sup>2</sup> <sup>=</sup>*c*<sup>44</sup> / *<sup>ρ</sup>*, *sl* 2 = *с*<sup>33</sup> / *ρ*, *ωme*<sup>4</sup> = *gM*0*h me*<sup>4</sup> = *gb*<sup>44</sup> <sup>2</sup>*M*<sup>0</sup> <sup>3</sup> / *<sup>c</sup>*44, *ω*<sup>10</sup> <sup>=</sup> *gM*<sup>0</sup> *<sup>h</sup> me*<sup>4</sup> <sup>−</sup>sin<sup>2</sup> *θ*(*β*˜ <sup>1</sup> +(*β*˜ <sup>2</sup> + 2*β*2)*M*<sup>0</sup> 2 cos<sup>2</sup> *<sup>θ</sup>* <sup>+</sup> *<sup>h</sup> me*sin<sup>2</sup> *θ*) .

Dispersion equation may be obtained from this system of equations. Solution of dispersion equation gives the propagation constant *k* for the waves with opposite circular polarization. We should note that due to the symmetry of the problem, the propagation of the wave with fixed circular polarization along *z*-axis is similar to the propagation of the wave with opposite circular polarization in inverse direction. Detailed investigation shows that waves with opposite circular polarizations have different propagation constants. This situation is valid for electromagnetic waves, as well as for acoustic ones. The difference in propagation constants will lead to different speed of waves with opposite circular polarization, hence, the rotation of polarization plane of linearly polarized wave. Rotation angle may be calculated as *Δφ* =(*k*<sup>+</sup> −*k*−)*l* / 2, where *k*+ and *k*− are propagation constants of the waves with corresponding circular polarization, and *l* is the distance, which the wave goes in the medium.

For investigation of electromagnetic wave reflection from the plate of magnet with spiral magnetic structure, one should add the boundary conditions: *Hx*,*<sup>y</sup>* (*out*) =*Hx*,*<sup>y</sup>* (*in*) , *Ex*,*<sup>y</sup>* (*out*) =*Ex*,*<sup>y</sup>* (*in*) , *Dz* (*out*) =*Dz* (*in*) , *∂m<sup>i</sup>* /*∂z* = 0, and σ*jz* = 0. Indexes (*in*) and (*out*) denote the fields inside and outside the magnet. We suppose that the outer medium is a vacuum. When one investigates the reflection of electromagnetic waves from semi-infinite structure, these boundary conditions should be applied to the single boundary.

For numerical estimations we will use following constant values [2]: *bij* ∼ 20 erg/(Oe × cm<sup>4</sup> ), ρ ∼ 10 g/cm3 , *ν<sup>t</sup>* ∼ 3 × 105 cm/s, *ν<sup>l</sup>* ∼ 5 × 105 cm/s, *M*<sup>0</sup> ∼ 500 Oe. The period of the structure for relativistic interaction-caused spiral magnet is usually much larger that in case of exchangecaused spiral structures. For example, Fe*x*Co1−*<sup>x</sup>*Si alloys, which symmetry allows Dzyaloshinskii-Moriya interaction, for *x* = 0.3 in spiral state has a modulation period *L* = 230 nm (*q* ∼ 3 × 105 cm−1) [34]. Other examples of the magnets with Dzyaloshinskii-Moriya interaction-caused spiral structures are FeGe (*L* = 70 nm, *q* ∼ 8 × 10<sup>5</sup> cm−1) [35], and MnSi (*L* = 18 nm, *q* ∼ 3 × 10<sup>6</sup> cm−1) [36]. Different modulated states exist in erbium single crystal due to the competing exchange interaction. In conical state, wave number of the structure is 5*c*\* /21 (*c*\* = 2π/*c* is inverse lattice constant, *c* = 0.56 nm is lattice constant) [37], i.e. *q* ∼ 3 × 10<sup>7</sup> cm−1. We will use *q* ∼ 10<sup>5</sup> cm−1 and *q* ∼ 108 cm−1, consequently.

First of all, let us discuss the influence of spin-electromagnetic interaction on wave propaga‐ tion. For such purposes, the problem may be simplified. We may neglect the elasticity and magnetostriction terms in free energy (1). In such case, dispersion equation may be calculated in form:

Radiation and Propagation of Waves in Magnetic Materials with Helicoidal Magnetic Structure http://dx.doi.org/10.5772/64014 167

( ) (( ) )( ) ( ( ) ) (( ( ) ) )( ) (( ) )(( ( ) ) ) 2 2 22 2 22 2 2 1 2 22 2 2 2 2 2 2 2 1 2 2 4 2 22 2 2 2 2 1 2 2 sin 2 cos 2 2 2 sin 2 cos sin 2 0. *k MM kq M M M k M kq M vk vk vk q vk q v k vk q* + - - + - é ù w w - + w - W -w w W q- w q ´ ë û é ù w w- - + ê ú - + w - - W -w w W q+ w q ë û -W q w - W - w w w - - W - w w =

Here, *ω<sup>M</sup>* =4*πgM*0 is magnetostatic frequency, Ω1*<sup>M</sup>* =Ω1*<sup>k</sup>* <sup>−</sup>*<sup>q</sup>* + *ω<sup>M</sup>* , Ω*<sup>k</sup>* <sup>−</sup>2*<sup>q</sup>* <sup>−</sup> <sup>=</sup> *gM*<sup>0</sup> (*α*(*<sup>q</sup>* <sup>2</sup> <sup>−</sup>(*<sup>k</sup>* <sup>−</sup>2*q*)2) <sup>+</sup> *<sup>γ</sup>*(*<sup>q</sup>* <sup>4</sup> <sup>−</sup>(*<sup>k</sup>* <sup>−</sup>2*q*)4) <sup>−</sup>2*α*<sup>1</sup> (*k* −*q*)), *Ω<sup>k</sup>* <sup>+</sup> <sup>=</sup> *gM*0(*α*(*<sup>q</sup>* <sup>2</sup> <sup>−</sup>*<sup>k</sup>* 2) <sup>+</sup> *<sup>γ</sup>*(*<sup>q</sup>* <sup>4</sup> <sup>−</sup>*<sup>k</sup>* 4) <sup>+</sup> <sup>2</sup>*α*<sup>1</sup> (*k* −*q*)), Ω1*<sup>k</sup>* <sup>−</sup>*<sup>q</sup>* = *gM*0(*β*<sup>1</sup> + 3*β*2*M*<sup>0</sup> 2 cos<sup>2</sup> *<sup>θ</sup>* <sup>−</sup>*α*(*<sup>q</sup>* <sup>2</sup> <sup>−</sup>(*<sup>k</sup>* <sup>−</sup>*q*)2) <sup>+</sup> *<sup>γ</sup>*(*<sup>q</sup>* <sup>4</sup> <sup>−</sup>(*<sup>k</sup>* <sup>−</sup>*q*)4) <sup>+</sup> <sup>2</sup>*α*1*q*). For relativistic spiral, one should put α1 = 0, while in exchange spiral case γ = 0.

Results of calculations are shown in **Figure 1**.

sin q -+ + + q = é ù *mkq mkq mk* - + ( ) ( ) 2cos 0. *<sup>z</sup>* ( ) ë û

In all equations above we have used the following notations: **h**(*k*), **m**(*k*), **u**(*k*) are Fourier components of corresponding vector **h**, **m**, **u**; *v* =*c* / *ε* is electromagnetic wave speed in the

) <sup>+</sup> <sup>2</sup>*α*1(*<sup>q</sup>* <sup>∓</sup> *<sup>k</sup>*), *<sup>L</sup>* ||(*k*)= <sup>−</sup>*α*(*<sup>q</sup>* <sup>2</sup> <sup>−</sup>*<sup>k</sup>* <sup>2</sup>

Dispersion equation may be obtained from this system of equations. Solution of dispersion equation gives the propagation constant *k* for the waves with opposite circular polarization. We should note that due to the symmetry of the problem, the propagation of the wave with fixed circular polarization along *z*-axis is similar to the propagation of the wave with opposite circular polarization in inverse direction. Detailed investigation shows that waves with opposite circular polarizations have different propagation constants. This situation is valid for electromagnetic waves, as well as for acoustic ones. The difference in propagation constants will lead to different speed of waves with opposite circular polarization, hence, the rotation of polarization plane of linearly polarized wave. Rotation angle may be calculated as *Δφ* =(*k*<sup>+</sup> −*k*−)*l* / 2, where *k*+ and *k*− are propagation constants of the waves with corresponding

For investigation of electromagnetic wave reflection from the plate of magnet with spiral magnetic structure, one should add the boundary conditions:

fields inside and outside the magnet. We suppose that the outer medium is a vacuum. When one investigates the reflection of electromagnetic waves from semi-infinite structure, these

For numerical estimations we will use following constant values [2]: *bij* ∼ 20 erg/(Oe × cm<sup>4</sup>

the competing exchange interaction. In conical state, wave number of the structure is 5*c*\*

= 2π/*c* is inverse lattice constant, *c* = 0.56 nm is lattice constant) [37], i.e. *q* ∼ 3 × 10<sup>7</sup>

cm−1, consequently.

First of all, let us discuss the influence of spin-electromagnetic interaction on wave propaga‐ tion. For such purposes, the problem may be simplified. We may neglect the elasticity and magnetostriction terms in free energy (1). In such case, dispersion equation may be calculated

relativistic interaction-caused spiral magnet is usually much larger that in case of exchangecaused spiral structures. For example, Fe*x*Co1−*<sup>x</sup>*Si alloys, which symmetry allows Dzyaloshinskii-Moriya interaction, for *x* = 0.3 in spiral state has a modulation period *L* = 230 nm (*q* ∼ 3 × 105 cm−1) [34]. Other examples of the magnets with Dzyaloshinskii-Moriya

, *ν<sup>t</sup>* ∼ 3 × 105 cm/s, *ν<sup>l</sup>* ∼ 5 × 105 cm/s, *M*<sup>0</sup> ∼ 500 Oe. The period of the structure for

cm−1) [36]. Different modulated states exist in erbium single crystal due to

<sup>3</sup> / *<sup>c</sup>*44, *ω*<sup>10</sup> <sup>=</sup> *gM*<sup>0</sup> *<sup>h</sup> me*<sup>4</sup> <sup>−</sup>sin<sup>2</sup>

circular polarization, and *l* is the distance, which the wave goes in the medium.

boundary conditions should be applied to the single boundary.

interaction-caused spiral structures are FeGe (*L* = 70 nm, *q* ∼ 8 × 10<sup>5</sup>

cm−1 and *q* ∼ 108

<sup>±</sup> <sup>=</sup>*ω*<sup>20</sup> <sup>+</sup> *gM*0*<sup>L</sup>* <sup>⊥</sup>

±(*k*), *ω*1*<sup>k</sup>* <sup>=</sup>*ω*<sup>10</sup> <sup>+</sup> *gM*0sin<sup>2</sup>

2 cos<sup>2</sup>

) + 2*α*1*q*, *st*

*<sup>θ</sup>* <sup>+</sup> *<sup>h</sup> me*sin<sup>2</sup>

)−*γ*(*q* <sup>4</sup> −*k* <sup>4</sup>

/*∂z* = 0, and σ*jz* = 0. Indexes (*in*) and (*out*) denote the

*θ*(*β*˜ <sup>1</sup> +(*β*˜ <sup>2</sup> + 2*β*2)*M*<sup>0</sup>

*θL* <sup>∥</sup>(*k*)=0,

2 =

), ρ

/21 (*c*\*

cm−1. We

cm−1) [35], and MnSi (*L* =

<sup>2</sup> <sup>=</sup>*c*<sup>44</sup> / *<sup>ρ</sup>*, *sl*

*θ*) .

medium; ε is permittivity of the crystal; *ω*2*<sup>k</sup>*

<sup>2</sup>*M*<sup>0</sup>

)−*γ*(*q* <sup>4</sup> −*k* <sup>4</sup>

*L* <sup>⊥</sup> ±

166 Magnetic Materials

*Hx*,*<sup>y</sup>* (*out*) =*Hx*,*<sup>y</sup>* (*in*) , *Ex*,*<sup>y</sup>* (*out*) =*Ex*,*<sup>y</sup>* (*in*) , *Dz* (*out*) =*Dz* (*in*) , *∂m<sup>i</sup>*

∼ 10 g/cm3

18 nm, *q* ∼ 3 × 10<sup>6</sup>

will use *q* ∼ 10<sup>5</sup>

in form:

(*k*)= −*α*(*q* <sup>2</sup> −*k* <sup>2</sup>

*с*<sup>33</sup> / *ρ*, *ωme*<sup>4</sup> = *gM*0*h me*<sup>4</sup> = *gb*<sup>44</sup>

**Figure 1a** shows the dispersion of coupled spin-electromagnetic waves in spiral magnets with different spiral angles (or external magnetic field values). One can see that the spectrum has a band structure. The band gap is observed near the point where non-interacting dispersion curves are crossed. This band gap appears due to the resonant interaction of spin and electro‐ magnetic waves in a magnet. From **Figure 1a** we can see that with decrease of angle, the gap shifts toward lower frequencies and its width decreases. Calculations show that in case of exchange spiral, a zone of opacity is much narrower than in case of relativistic spiral, and lies at higher frequencies. It should be also noted that the magnitude of the interaction of spin and electromagnetic waves depends on the spiral angle.

**Figure 1b** shows the frequency dependence of the reflectance of electromagnetic wave from the plate of spiral magnet with width *l* = 100 μm for different spiral angles. One can see that reflectance is near to unit at frequencies corresponding to the band gap of the spectrum of coupled waves, i.e. opacity window is observed. This window exists only for spiral states with θ ≠ 0, while it disappears in the collinear ferromagnetic state.

The results of calculating polarization plane rotation are shown in **Figure 1c** and **d**. One can see that rotation angle depends on external magnetic field almost linearly at frequencies upper and lower the band gap. The maximum values of rotation correspond to the band gap boundaries.

Let us now investigate the effect of acoustic subsystem on wave propagation in the magnet. For such purposes we should solve all equations for exact free energy expression (1). The details of this investigation may be found in work [33].

Calculations of the spectrum of coupled spin-electromagnetic-acoustic waves in spiral magnets show that despite the gap corresponding to spin-electromagnetic interaction, an additional gap corresponding to spin-acoustic interaction is opened. This gap reaches a maximal width at ferromagnetic spiral state with θ = π/4. The spectrum of coupled waves near *k* = 0 is shown in **Figure 2a**.

**Figure 1.** Dispersion of coupled spin-electromagnetic waves in spiral magnets with different spiral angles (a). On the inset the small *k* region is shown. Frequency dependence of the reflectance of electromagnetic wave from the plate of spiral magnet with width *l* = 100 μm for different spiral angles (b). Electromagnetic waves polarization plane rotation: near upper frequency of band gap (c), and near lower frequency of band gap (d).

**Figure 2b** shows results of investigation of electromagnetic waves reflection from semi-infinite spiral magnet. One can see that besides the opacity window corresponding to spin-electro‐ magnetic interaction, there are some features at the frequencies of spin-acoustic interaction. This effect corresponds to resonant interaction of electromagnetic, spin and acoustic waves. The width of such acoustically caused electromagnetic opacity window depends on the spiral angle (or, equivalently, on external magnetic field value), as well.

The gap opening in spin-acoustic spectrum leads to the possibility of polarization plane rotation for acoustic waves, i.e. acoustic Faraday effect. The results of calculation of rotation angle at the sample distance of 1 cm for different frequencies are shown in **Figure 2c**. The breaks on the curves correspond to band gap in the spectrum, which moves with the change of the external magnetic field. One can see that maximal rotation of the polarization plane is observed near the band gap.

The results we have discussed in this section show that spiral magnets may be perspective materials for electromagnetic and acoustic wave manipulation (such as filtering, polarizing, and phase shifting).

Radiation and Propagation of Waves in Magnetic Materials with Helicoidal Magnetic Structure http://dx.doi.org/10.5772/64014 169

**Figure 2.** Spectrum of coupled waves (a) (left – schematic spectrum, right – calculated spectrum near *k* = 0). Electro‐ magnetic waves reflection from semi-infinite spiral magnet (b). On the inset the feature associated with the acoustic gap is shown. Acoustic Faraday effect (c).

## **3. Radiation of electromagnetic and acoustic waves by helicoidal magnets at phase transitions**

In this section we partly follow our previous works [38, 39].

**Figure 1.** Dispersion of coupled spin-electromagnetic waves in spiral magnets with different spiral angles (a). On the inset the small *k* region is shown. Frequency dependence of the reflectance of electromagnetic wave from the plate of spiral magnet with width *l* = 100 μm for different spiral angles (b). Electromagnetic waves polarization plane rotation:

**Figure 2b** shows results of investigation of electromagnetic waves reflection from semi-infinite spiral magnet. One can see that besides the opacity window corresponding to spin-electro‐ magnetic interaction, there are some features at the frequencies of spin-acoustic interaction. This effect corresponds to resonant interaction of electromagnetic, spin and acoustic waves. The width of such acoustically caused electromagnetic opacity window depends on the spiral

The gap opening in spin-acoustic spectrum leads to the possibility of polarization plane rotation for acoustic waves, i.e. acoustic Faraday effect. The results of calculation of rotation angle at the sample distance of 1 cm for different frequencies are shown in **Figure 2c**. The breaks on the curves correspond to band gap in the spectrum, which moves with the change of the external magnetic field. One can see that maximal rotation of the polarization plane is observed

The results we have discussed in this section show that spiral magnets may be perspective materials for electromagnetic and acoustic wave manipulation (such as filtering, polarizing,

near upper frequency of band gap (c), and near lower frequency of band gap (d).

angle (or, equivalently, on external magnetic field value), as well.

near the band gap.

168 Magnetic Materials

and phase shifting).

In context of macroscopic electrodynamics, electromagnetic wave radiation may be described by the following equation:

$$
\Delta \mathbf{E} - \frac{1}{c^2} \frac{\partial^2 \mathbf{E}}{\partial t^2} = \frac{4\pi}{c^2} \left[ \frac{\partial}{\partial t} \left( \mathbf{j} + \frac{\partial \mathbf{P}}{\partial t} + c \operatorname{rot} \mathbf{M} \right) + c^2 \operatorname{grad} \operatorname{div} \left( \mathbf{P} \right) \right] \tag{3}
$$

In (3) **E** is electric field strength, **P** and **M** are polarization and magnetization of the medium, consequently, **j** is the current density, *c* is speed of light in the vacuum. Assuming magnetic material to be non-conductive with dielectric constant ε, Eq. (3) may be simplified:

$$
\Delta \mathbf{E} - \frac{\varepsilon}{c^2} \frac{\partial^2 \mathbf{E}}{\partial t^2} = \frac{4\pi}{c^2} \frac{\partial \mathbf{j\_M}}{\partial t} \tag{4}
$$

In (4) magnetic current density has been introduced: **jM** = *c*·rot(**M**). For spiral magnets with magnetization distribution given by (2), one may calculate magnetic current density **jM** = (*j*0cos(*qz*), *j*0sin(*qz*), 0), *j*0 = −*cqM*0sin(θ). So, we may conclude that any process with changing the spiral wave number *q*, the amplitude of magnetization vector *M*0, or the spiral angle θ will be accompanied by electromagnetic radiation.

The form of radiation Eq. (4) allows one to see that time-varying non-uniform magnetization distribution is equivalent to time-varying distributed conductive currents. Generally, this equation should be solved together with the magnetization motion equation (usually, Landau-Lifshitz equation with additional relaxation terms). This system of equations can be solved only numerically. However, it is known that the frequency of precession of the magnetization vector in ferromagnetic materials (ferromagnetic resonance frequency) is about 10–100 GHz, and the relaxation time of magnetization to equilibrium state is of the order of few microsec‐ onds at usual conditions. In this case, if we are interested in processes with characteristic frequencies far from resonance, the motion of magnetization can be neglected, and one can assume that magnetization has always ground state value (2). This makes reasonable the assumption that only spiral angle changes during magnetic field induced phase transition.

For investigation of sound waves generation, we should solve equation of motion for elastic medium:

$$
\rho \Box^2 u\_i / \partial t^2 = \partial \sigma\_y / \partial \mathfrak{x}\_j \,, \sigma\_y = \partial F / \partial u\_y \, . \tag{5}
$$

In (5), *F* denotes free energy density defined by Eq. (1) for spiral magnets. Eq. (5), similarly to the case of electromagnetic waves, may be transformed into wave equations with non-zero source functions, which in general may depend on both coordinate and time. It is convenient to introduce circular components (*u*, *M*)±= (*u*, *M*)*x* ± *i*(*u*, *M*)*y*, transversal and longitudinal sound velocities *vt* = (2*c*44/)1/2 and *vl* = (*c*33/)1/2, consequently. Calculations show that only transverse sound will be excited. We will have

$$\frac{\partial^2 u\_\pm}{\partial t^2} - \nu\_t^2 \frac{\partial^2 u\_\pm}{\partial \boldsymbol{\varpi}^2} = \pm \frac{iqb\_{44}M\_0^2}{\rho} \sin(2\theta) \exp(\pm iq\boldsymbol{\varpi}).\tag{6}$$

Let's consider an infinite crystal. At initial time, it will be assumed that there was no electro‐ magnetic radiation, i.e., we should solve the Eq. (4) or (6) with zero initial conditions. In this case, solutions of Eq. (4) or (6) can be obtained from the Duhamel's principle [40].

Here, we would not solve the problem of phase transition kinetics. We will consider only two simplest cases: an infinitely fast phase transition and linear variation of source functions. Comparing Eqs. (6) and (4) one may note that source function for radiation of electromagnetic waves is proportional to time derivative of sin[θ(*t*)], while in case of acoustic waves it is defined by time dependence of sin[2θ(*t*)]. Thus, maximal radiation of electromagnetic waves may be expected at phase transition from simple spiral state to collinear ferromagnetic one, while maximal acoustic radiation should be observed at phase transition from "ferromagnetic spiral" phase with θ = π/4 to collinear ferromagnetic state.

2 22 2

be accompanied by electromagnetic radiation.

medium:

170 Magnetic Materials

velocities *vt*

= (2*c*44/)1/2 and *vl*

sound will be excited. We will have

4 *ct c t* e¶ p ¶ D- =

**<sup>E</sup> jM <sup>E</sup>** (4)

¶ ¶

In (4) magnetic current density has been introduced: **jM** = *c*·rot(**M**). For spiral magnets with magnetization distribution given by (2), one may calculate magnetic current density **jM** = (*j*0cos(*qz*), *j*0sin(*qz*), 0), *j*0 = −*cqM*0sin(θ). So, we may conclude that any process with changing the spiral wave number *q*, the amplitude of magnetization vector *M*0, or the spiral angle θ will

The form of radiation Eq. (4) allows one to see that time-varying non-uniform magnetization distribution is equivalent to time-varying distributed conductive currents. Generally, this equation should be solved together with the magnetization motion equation (usually, Landau-Lifshitz equation with additional relaxation terms). This system of equations can be solved only numerically. However, it is known that the frequency of precession of the magnetization vector in ferromagnetic materials (ferromagnetic resonance frequency) is about 10–100 GHz, and the relaxation time of magnetization to equilibrium state is of the order of few microsec‐ onds at usual conditions. In this case, if we are interested in processes with characteristic frequencies far from resonance, the motion of magnetization can be neglected, and one can assume that magnetization has always ground state value (2). This makes reasonable the assumption that only spiral angle changes during magnetic field induced phase transition.

For investigation of sound waves generation, we should solve equation of motion for elastic

In (5), *F* denotes free energy density defined by Eq. (1) for spiral magnets. Eq. (5), similarly to the case of electromagnetic waves, may be transformed into wave equations with non-zero source functions, which in general may depend on both coordinate and time. It is convenient to introduce circular components (*u*, *M*)±= (*u*, *M*)*x* ± *i*(*u*, *M*)*y*, transversal and longitudinal sound

( ) ( ) 2 2 <sup>2</sup>

± ± ¶ ¶ - =± q ±

case, solutions of Eq. (4) or (6) can be obtained from the Duhamel's principle [40].

2 2 sin 2 exp . *<sup>t</sup>*

*v iqz*

Let's consider an infinite crystal. At initial time, it will be assumed that there was no electro‐ magnetic radiation, i.e., we should solve the Eq. (4) or (6) with zero initial conditions. In this

2 44 0

*u u iqb M*

¶ ¶ r

*t z*

2 2 , . *<sup>i</sup> ij j ij ij* r¶ ¶ = ¶s ¶ s = ¶ ¶ *u t x Fu* (5)

= (*c*33/)1/2, consequently. Calculations show that only transverse

(6)

Let us assume that sin[θ(*t*)] = Θ(-*t*); Θ(*t*) is the step-like Heaviside function; ∂Θ(-*t*)/∂*t* = −δ(*t*); δ(*t*) is the Dirac's delta function. This situation corresponds to instantaneous actuation of magnetic field with the value *Hcr* at time *t* = 0, and collapse of the spin spiral. Such a process is quite unreal since magnetization requires finite relaxation time to get to the final collinear ferromagnetic state. However, despite the unreality of the case, its analysis allows us to estimate an upper limit of emitted wave amplitude. Calculations show that electromagnetic field has a form:

$$E\_{\pm}\left(z,t\right) = E\_{x}\left(z,t\right) \pm iE\_{y}\left(z,t\right) = \mp 4\pi iM\_{o}\varepsilon^{-l/2}\sin\left(q\gamma t\right)\exp\left[\mp iqz\right] \tag{7}$$

In case of linear time dependence of sin[θ(*t*)]: sin[θ(*t*)] = *Pt*, at *t* ≤ *P*−1, *P* = *const*, we will have:

$$E\_{\pm}\left(z,t\right) = \pm 4\pi i M\_{\circ} P\left(qc\right)^{-1} \left[1 - \cos\left(q\nu t\right)\right] \exp\left[\mp iqz\right], t \le P^{-1} \tag{8}$$

It is seen that in contrast to the fast actuation of the magnetic field, in addition to the standing wave, a constant component of the electric field appears as well. The amplitude and the constant component are directly proportional to the speed of the state change *P*.

Solutions given by expressions (7) and (8) are superposition of waves travelling to the left and to the right, which, due to unbound periodic source, give a standing wave by adding each other. In the real material, multiple reflections of the excited waves from the boundaries and some part of electromagnetic energy will be emitted beyond the material. If refractive index of the environment is the same as for crystal (i.e. mediums are ideally conjugated) with thickness *d*, two electromagnetic wave pulses with wave number *q* and frequency *qv* will be emitted from each sample's side.

Maximum sound waves generation will take place at phase transition from θ = π/4 to collinear ferromagnetic state. Mathematically, it may be expressed as sin[2θ(*t*)] = Θ(*t*). In such case, solution may be easily calculated:

$$\ln u\_{\pm}(z,t) = \pm i \, b\_{44} M\_0^2 \exp\left[\pm iqz\right] \left[1 - \cos\left(q\mathbf{v}\_i t\right)\right] \left\langle \mathbf{p} q \mathbf{v}\_i^2 \right\rangle \tag{9}$$

Let us consider now the model process, when sin[2θ(*t*)] = η*t*, *t* < η−1 (i.e. linear time-dependence of source function). Calculations show that the following sound oscillations will be excited:

$$u\_{\pm}\left(z,t\right) = \pm \eta i b\_{44} M\_{0}^{2} \mathfrak{p}^{-1} q^{-2} \nu\_{\prime}^{-3} \left[q \nu\_{\prime} t - \sin\left(q \nu\_{\prime} t\right)\right] \exp\left[\pm iqz\right], t < \eta^{-1}. \tag{10}$$

The frequency of standing electromagnetic waves in formulas (7) and (8) is ω = *qv* ∼ 1015 − 1016 s−1. At "slow" change of the state, in mode determined by formula (8), the signal amplitude is proportional to *P*. For *P* ∼ 107 s−1, the amplitude of emitted waves is small: *E* ∼ 1012 CGSE (energy density is about 10−11 J/cm3 ). For ultrafast magnetization reversal processes, the time of magnetization switch is usually about 10−12 s, or *P* ∼ 1012 s−1, and we have *E* ∼ 10<sup>2</sup> CGSE (energy density is about 10−3 J/cm3 ). The maximal radiation corresponding to the situation of the instantaneous change of the state has amplitude of electric field *E* ∼ 10<sup>3</sup> CGSE (energy density is about 0.1 J/cm<sup>3</sup> ). This radiation will be emitted like a pulse. Characteristic time length of this pulse is *d*/*v* + *P*−1. For sample with *d* ∼ 1 μm, we will have *d*/*v* ∼ 10−14 s. This time is much smaller than usual magnetization switching time, so in real experiments, characteristic time length of the pulse should be comparable with switching time for *d*≪*vP*−1. Despite low energy density of the pulse, it may have valuable power density. In case of *P* = 1012 s−1 this value may reach about 1 GW/cm<sup>3</sup> .

For acoustic wave generation in case of infinitely fast phase transition, we will have the oscillations frequency ω = *qvt* . For Dzyaloshinskii-Moriya interaction-caused spiral magnets the linear frequency *f* = ω/(2π) ∼ 5 GHz, and amplitude is *u*<sup>±</sup> ∼ 10−5 cm. In exchange-caused case *f* ∼ 5 THz, *u*<sup>±</sup> ∼ 10−8 cm. Both cases correspond to hypersound oscillations. In case of "slow" phase transition sound oscillations with the same frequency are excited. Amplitude of excited oscillations linearly depends on the "phase transition speed" η. We may calculate that *u*<sup>±</sup> ∼ 10−21 × η cm and *u*± ∼ 10−27 × η cm for Dzyaloshinskii-Moriya and exchange interaction-caused spiral magnets, respectively. For example, ultrafast magnetization reversal processes have time of magnetization switch of about 10−12 s, or η ∼ 1012 s−1. The amplitudes will have values *u*<sup>±</sup> ∼ 10−9 cm and *u*<sup>±</sup> ∼ 10−15 cm for Dzyaloshinskii-Moriya and exchange interaction-caused spiral magnets, consequently. One can see that in case of exchange spiral structure, the amplitude is negligibly small. In experiments it is possible to detect the pulse of sound with characteristic time length of about τ ∼ *d*/*vt* , where *d* is the sample size. For sample size *d* ∼ 1 mm we will have time length τ ∼ 10−6 s.

The abovementioned results make one believe that a similar effect may be observed when the spiral magnet is placed in homogeneous periodic magnetic field, i.e. when the spiral angle periodically changes near its equilibrium value. For investigation of this problem, let us consider the spiral magnet in time-dependent external magnetic field *Hz*(*t*) = *H*0 + *h*0(*t*), where *H*0 is a constant magnetic field defining the ground state and equilibrium value of spiral angle, and *h*0(*t*) = *h*0exp[*i*Ω*t*] is a periodic component. For simplicity, we will assume that the periodic component of the field has small amplitude. In such context, magnetization may be described by static component **M0**, which is described by Eq. (2) with the spiral angle defined by *H*0, and by dynamical component *m*(*z*,*t*). Dynamical component may be obtained from Landau-Lifshitz equation. We will neglect relaxation processes for simplicity. Oscillated components may be assumed in form of Fourier series *m*(*z*,*t*) = Σ*mn*(*t*)exp[−*inqz*]. Calculations show that longitudinal periodic magnetic field will excite only harmonics of *mx* and *my* with *n* = ±1. We will have the following:

Let us consider now the model process, when sin[2θ(*t*)] = η*t*, *t* < η−1 (i.e. linear time-dependence of source function). Calculations show that the following sound oscillations will be excited:

The frequency of standing electromagnetic waves in formulas (7) and (8) is ω = *qv* ∼ 1015 − 1016 s−1. At "slow" change of the state, in mode determined by formula (8), the signal amplitude is proportional to *P*. For *P* ∼ 107 s−1, the amplitude of emitted waves is small: *E* ∼ 1012 CGSE

of magnetization switch is usually about 10−12 s, or *P* ∼ 1012 s−1, and we have *E* ∼ 10<sup>2</sup> CGSE

of this pulse is *d*/*v* + *P*−1. For sample with *d* ∼ 1 μm, we will have *d*/*v* ∼ 10−14 s. This time is much smaller than usual magnetization switching time, so in real experiments, characteristic time length of the pulse should be comparable with switching time for *d*≪*vP*−1. Despite low energy density of the pulse, it may have valuable power density. In case of *P* = 1012 s−1 this value may

For acoustic wave generation in case of infinitely fast phase transition, we will have the

the linear frequency *f* = ω/(2π) ∼ 5 GHz, and amplitude is *u*<sup>±</sup> ∼ 10−5 cm. In exchange-caused case *f* ∼ 5 THz, *u*<sup>±</sup> ∼ 10−8 cm. Both cases correspond to hypersound oscillations. In case of "slow" phase transition sound oscillations with the same frequency are excited. Amplitude of excited oscillations linearly depends on the "phase transition speed" η. We may calculate that *u*<sup>±</sup> ∼ 10−21 × η cm and *u*± ∼ 10−27 × η cm for Dzyaloshinskii-Moriya and exchange interaction-caused spiral magnets, respectively. For example, ultrafast magnetization reversal processes have time of magnetization switch of about 10−12 s, or η ∼ 1012 s−1. The amplitudes will have values *u*<sup>±</sup> ∼ 10−9 cm and *u*<sup>±</sup> ∼ 10−15 cm for Dzyaloshinskii-Moriya and exchange interaction-caused spiral magnets, consequently. One can see that in case of exchange spiral structure, the amplitude is negligibly small. In experiments it is possible to detect the pulse of sound with

The abovementioned results make one believe that a similar effect may be observed when the spiral magnet is placed in homogeneous periodic magnetic field, i.e. when the spiral angle periodically changes near its equilibrium value. For investigation of this problem, let us consider the spiral magnet in time-dependent external magnetic field *Hz*(*t*) = *H*0 + *h*0(*t*), where *H*0 is a constant magnetic field defining the ground state and equilibrium value of spiral angle, and *h*0(*t*) = *h*0exp[*i*Ω*t*] is a periodic component. For simplicity, we will assume that the periodic component of the field has small amplitude. In such context, magnetization may be described by static component **M0**, which is described by Eq. (2) with the spiral angle defined by *H*0, and

the instantaneous change of the state has amplitude of electric field *E* ∼ 10<sup>3</sup>

<sup>±</sup> = ±h r é ù - ± <h ë û (10)

). For ultrafast magnetization reversal processes, the time

). The maximal radiation corresponding to the situation of

). This radiation will be emitted like a pulse. Characteristic time length

. For Dzyaloshinskii-Moriya interaction-caused spiral magnets

, where *d* is the sample size. For sample size *d* ∼ 1

CGSE (energy

( ) ( ) [ ] <sup>2123</sup> <sup>1</sup> 44 0 , sin exp , . *tt t u z t ib M q v qv t qv t iqz t* --- -

(energy density is about 10−11 J/cm3

(energy density is about 10−3 J/cm3

.

characteristic time length of about τ ∼ *d*/*vt*

mm we will have time length τ ∼ 10−6 s.

density is about 0.1 J/cm<sup>3</sup>

172 Magnetic Materials

reach about 1 GW/cm<sup>3</sup>

oscillations frequency ω = *qvt*

$$\begin{split} \frac{\partial^2 m\_{\pm}}{\partial z \partial t} &= \frac{-\text{g}q \, \Omega M\_{0} \, \sin \theta \, h\_{0z} e^{q\Omega}}{\Omega\_{\text{res}}^2 - \Omega^2} \Big\{ \pm 2 \, \text{g}\alpha\_{l} q M\_{0} \, \cos \theta e^{\mp qz} + \left( \Omega \pm \, \text{g}\alpha\_{l} q M\_{0} \cos \theta \right) e^{\pm qz} \Big\}; \\ m\_{\pm} &= m\_{\pm} \pm i m\_{\pm}; \, \Omega\_{\text{res}}^2 = \text{S} \alpha\_{l}^2 q^2 q^2 M\_{0}^2 \cos^2 \theta. \end{split} \tag{11}$$

One can see that the source function defined by Eq. (11) has a resonance at Ω = Ω*res*. This is due to resonant behaviour of corresponding susceptibility component. In real materials, due to unavoidable damping processes, the source function will have a finite value even at resonant frequency. We will take into account these processes formally by putting Ω + *i*δ instead of Ω, where δ is a formal damping parameter.

Solving Eq. (4) with use of Green's function formalism for magnet's size *d*, we will have the following stationary electric field amplitudes at the left boundary of the sample:

$$\begin{split} E\_{\pm}\left(0, t\right) &= \frac{i2\pi gM\_{o}\sin\Theta h\_{0\pm}e^{i\Omega t}e^{-i\Omega d\eta}}{\Omega\_{\text{rx}}^{2} - \left(\Omega + i\gamma\right)^{2}} \times \\ &\times \left\{ \frac{\Omega + i\gamma \pm \text{g}\alpha\_{1}qM\_{0}\cos\Theta}{\Omega \mp q\gamma} \Bigg[ q\upsilon\left(e^{\alpha\omega\mu\upsilon} - 1\right) \pm \Omega \left(1 - e^{\pm\omega\mu}\right) \right] \mp \\ &\mp \frac{2\text{g}\alpha\_{1}qM\_{0}\cos\Theta}{\Omega \pm q\upsilon} \Bigg[ q\upsilon\left(e^{\alpha\omega\mu\upsilon} - 1\right) \mp \Omega \left(1 - e^{\mp\omega\mu}\right) \Bigg] \end{split} \tag{12}$$

One may calculate the energy of emitted waves *W*± = |*E*±|2 /8π. Dependencies of this energy from the frequency and from the equilibrium spiral angle are shown in **Figure 3**.

One can see that intensities of left- and right-polarized emitted waves differ from each other. When 0 < θ < π/2 (i.e. when the static component of magnetization is tilted along spiral axis) the energy of negatively polarized waves is greater than that of positively polarized ones. When θ = π/2, energies of waves for both polarizations are equal. When π/2 < θ < π, the energy of positively polarized waves prevails. The increase of damping parameter leads to decrease of emitted wave energy, especially near the resonance. When the damping parameter becomes much greater than the resonant frequency, the energy of emitted waves slowly varies with the frequency, equal for waves of both polarizations, and has a maximum at θ = π/2.

**Figure 3.** Energy of emitted electromagnetic waves versus frequency of exciting field and the static value of spiral angle.

Similarly, one may consider acoustic wave excitation in spiral magnets in homogeneous harmonical magnetic field. Wave equations may be written as follows:

$$\begin{split} \frac{\partial^2 u\_\pm}{\partial \varpi^2} - \frac{\ddot{u}\_\pm}{\nu\_\times} &= -\frac{b\_{44}}{c\_{44}} \Big( M\_{\pm 0} \frac{\partial m\_\pm}{\partial \varpi} + m\_\pm \frac{\partial M\_{\pm 0}}{\partial \varpi} \Big); \\ \frac{\partial^2 u\_\pm}{\partial \varpi^2} - \frac{\ddot{u}\_\pm}{\nu\_l} &= -\frac{2b\_{13}}{c\_{33}} \Big( M\_{\pm 0} \frac{\partial m\_\pm}{\partial \varpi} + m\_\times \frac{\partial M\_{\times 0}}{\partial \varpi} + M\_{\pm 0} \frac{\partial m\_\wp}{\partial \varpi} + m\_\wp \frac{\partial M\_{\wp 0}}{\partial \varpi} \Big). \end{split} \tag{13}$$

Both longitudinal and transversal acoustic waves may be excited. Transverse magnetization oscillations will excite both longitudinal and transverse sound, while longitudinal magnetization oscillations affect only transverse sound. Similar to calculations made in [41], for simplicity we will suppose that magnetization oscillations are uniform, i.e. ∂*mi* (*z*, *t*)/ ∂ *z* ≈*mi*<sup>0</sup> *δ*(*z*)−*δ*(*z* −*d*) exp(−*iωt*), here δ(*z*) is Dirac's delta-function, *d* is the sample size. Wave equations may be solved by Green's function formalism. Green's functions for longitudinal and transverse acoustic waves *Gt,l*(*z*, *z*0) consist of two terms. The first one (*i*/ 2*kt,l*)exp(*ikt,l*|*z*–*z*0|) corresponds to the sound excited in point *z*<sup>0</sup> and propagating in +*z*direction. The second one (*i*/2*kt,l*)exp(*ikt,l*|*z*+*z*0|) corresponds to the wave excited in the same point, but propagating in opposite direction. In both cases, *kt,l* = *ω/vt,l*. Complete Green's function will be as follows:

Radiation and Propagation of Waves in Magnetic Materials with Helicoidal Magnetic Structure http://dx.doi.org/10.5772/64014 175

$$G\_{t,l}\left(z, z\_0\right) = ik\_{t,l}^{-1} \exp\left(-ik\_{t,l}z\right)\cos\left(k\_{t,l}z\_0\right). \tag{14}$$

Calculations show that amplitudes of excited waves are:

**Figure 3.** Energy of emitted electromagnetic waves versus frequency of exciting field and the static value of spiral

harmonical magnetic field. Wave equations may be written as follows:

44 0

*M m*

*z z*

*z z x x y y*

¶ æ ö ¶ ¶ ¶ ¶ - =- + + + ç ÷ ¶ ¶¶ ¶¶ è ø

*u u b mM m M*

*zv c zz zz*

2 0 0

Similarly, one may consider acoustic wave excitation in spiral magnets in homogeneous

13 0 0

*xx y y*

*MmMm*

Both longitudinal and transversal acoustic waves may be excited. Transverse magnetization oscillations will excite both longitudinal and transverse sound, while longitudinal magnetization oscillations affect only transverse sound. Similar to calculations made in [41], for simplicity we will suppose that magnetization oscillations are uniform, i.e.

(*z*, *t*)/ ∂ *z* ≈*mi*<sup>0</sup> *δ*(*z*)−*δ*(*z* −*d*) exp(−*iωt*), here δ(*z*) is Dirac's delta-function, *d* is the sample size. Wave equations may be solved by Green's function formalism. Green's functions for longitudinal and transverse acoustic waves *Gt,l*(*z*, *z*0) consist of two terms. The first one (*i*/ 2*kt,l*)exp(*ikt,l*|*z*–*z*0|) corresponds to the sound excited in point *z*<sup>0</sup> and propagating in +*z*direction. The second one (*i*/2*kt,l*)exp(*ikt,l*|*z*+*z*0|) corresponds to the wave excited in the same point, but propagating in opposite direction. In both cases, *kt,l* = *ω/vt,l*. Complete Green's

<sup>2</sup> .

;

(13)

angle.

174 Magnetic Materials

∂*mi*

2

2

function will be as follows:

2 0 44

*t*

&&

&&

*l*

33

*uu b m M*

± ± <sup>±</sup> <sup>±</sup> ¶ ¶ æ ö ¶ - =- + ç ÷ ¶ è ø ¶ ¶

*zv c z z*

$$\begin{aligned} \boldsymbol{u}\_{\pm} &= -\frac{i\boldsymbol{b}\_{44}M\_{0}}{k\_{\circ}c\_{44}} \Bigg[ \begin{aligned} &m\_{\circ 0}\cos\Theta \Big[ &1-\cos\big(k\_{i}d\big) \Big] \pm \\ &\pm iq \frac{m\_{\circ 0}\sin\Theta}{q^{2}-k\_{i}^{2}} \Bigg[ \begin{aligned} &k\_{i}\exp\big(\pm iqd\big)\sin\big(k\_{i}d\big) \mp \\ &\mp iq \left(1-\exp\big(\pm iqd\big)\cos\big(k\_{i}d\big)\Big) \end{aligned} \right] \Bigg] \end{aligned} \tag{15}$$
 
$$\boldsymbol{u}\_{\pm} = -\frac{4ik\_{\restriction 3}M\_{0}\sin\Theta}{k\_{\restriction l}c\_{33}} \begin{aligned} &\left[m\_{\simeq 0}\frac{k\_{l}}{q^{2}-k\_{l}^{2}}\left[q\sin\operatorname{qd}\sin k\_{i}d - k\_{i}\cos qd\cos k\_{i}d\right] - \\ &-m\_{\geqslant 0}\sin qd \Big[ \cos k\_{i}d - \frac{q\sin k\_{i}d}{q+k\_{l}} \Bigg] \end{aligned} \tag{15}$$

**Figure 4.** Fragments of acoustic waves energy density dependencies from frequency and sample size for different po‐ larization of excited field. Wave number of spiral is *q* = 105 cm−1.

Energy density of excited waves is *W* = *c*33(*uzz*) 2 /2 + 2*c*44[(*uxz*) 2 +(*uyz*) 2 ]. Fragments of its depend‐ encies from frequency and sample size for different polarization of excited field are shown in **Figure 4**. One can see that these dependencies have a complex oscillating behaviour. Oscilla‐ tions are caused by interference of sound waves propagating from opposite sample bounda‐ ries. Another reason of oscillations is periodicity of source function. Eq. (15) have resonances at *kt,l* = *q*. At microwave frequencies these resonances may be observed in magnets with relativistic spiral structure (*q*∼ 105 cm−1, *ωres*∼*qvt,l*∼ 1010 rad/s), while in magnets with exchange spiral structure, resonance frequencies correspond to terahertz frequency range (*q* ∼ 10<sup>8</sup> cm−1, *ωres* ∼ *qvt,l* ∼ 1013 rad/s).

## **4. Conclusions**

We have reviewed some features of electromagnetic and acoustic wave propagation in magnets with helicoidal spin order. There is an opacity window in reflectance spectrum of electromagnetic waves from the plate of spiral magnet. The width of this window has a maximal value at phase of simple spiral, and decreases with decrease of spiral angle (or, equivalently, with increase of external magnetic field). Field dependencies of polarization plane rotation angle have a resonant behaviour. Maximum rotation may be observed near boundaries of band gaps for both electromagnetic and acoustic waves. Both electromagnetic and acoustic waves may be radiated by spiral magnets during the phase transition. The amplitude of radiated waves strictly depends on the speed of phase transition.

In general, the results show the possibility of practical applications of spiral magnets for electromagnetic and acoustic wave manipulation and generation.

## **Acknowledgements**

The work was supported partly by Russian Science Foundation Grant # 14-22-00279 and Russian Foundation for Basic Researches Grant # 15-07-99654.

## **Author details**

Igor V. Bychkov1\*, Dmitry A. Kuzmin1 and Vladimir G. Shavrov2

\*Address all correspondence to: bychkov@csu.ru

1 Chelyabinsk State University, Chelyabinsk, Russian Federation

2 Kotel'nikov Institute of Radio Engineering and Electronics of RAS, Moscow, Russian Fed‐ eration

## **References**

Energy density of excited waves is *W* = *c*33(*uzz*)

*ωres* ∼ *qvt,l* ∼ 1013 rad/s).

**Acknowledgements**

**Author details**

eration

Igor V. Bychkov1\*, Dmitry A. Kuzmin1

\*Address all correspondence to: bychkov@csu.ru

**4. Conclusions**

176 Magnetic Materials

2

encies from frequency and sample size for different polarization of excited field are shown in **Figure 4**. One can see that these dependencies have a complex oscillating behaviour. Oscilla‐ tions are caused by interference of sound waves propagating from opposite sample bounda‐ ries. Another reason of oscillations is periodicity of source function. Eq. (15) have resonances at *kt,l* = *q*. At microwave frequencies these resonances may be observed in magnets with relativistic spiral structure (*q*∼ 105 cm−1, *ωres*∼*qvt,l*∼ 1010 rad/s), while in magnets with exchange spiral structure, resonance frequencies correspond to terahertz frequency range (*q* ∼ 10<sup>8</sup>

We have reviewed some features of electromagnetic and acoustic wave propagation in magnets with helicoidal spin order. There is an opacity window in reflectance spectrum of electromagnetic waves from the plate of spiral magnet. The width of this window has a maximal value at phase of simple spiral, and decreases with decrease of spiral angle (or, equivalently, with increase of external magnetic field). Field dependencies of polarization plane rotation angle have a resonant behaviour. Maximum rotation may be observed near boundaries of band gaps for both electromagnetic and acoustic waves. Both electromagnetic and acoustic waves may be radiated by spiral magnets during the phase transition. The

In general, the results show the possibility of practical applications of spiral magnets for

The work was supported partly by Russian Science Foundation Grant # 14-22-00279 and

2 Kotel'nikov Institute of Radio Engineering and Electronics of RAS, Moscow, Russian Fed‐

and Vladimir G. Shavrov2

amplitude of radiated waves strictly depends on the speed of phase transition.

electromagnetic and acoustic wave manipulation and generation.

Russian Foundation for Basic Researches Grant # 15-07-99654.

1 Chelyabinsk State University, Chelyabinsk, Russian Federation

/2 + 2*c*44[(*uxz*)

2 +(*uyz*) 2

]. Fragments of its depend‐

cm−1,


magnetic fields. IEEE Transactions on Magnetics. 2013;49(8):4695–4698. DOI: 10.1109/ TMAG.2013.2261055

[29] Buchelnikov V.D., Bychkov I.V., Kuzmin D.A., Lamekhov S.J., Shavrov V.G. Reflecting electromagnetic waves from a surface of TbMnO3 with sinusoidal antiferromagnetic structure. Bulletin of the Russian Academy of Sciences: Physics. 2013;77(9):1120–1122. DOI: 10.3103/S1062873813090062

[16] Sementsov D.I. Features of light propagation in helical magnetic structures. Optics and

[17] Manzhos I.V., Chupis I.E. Electromagnetic-spin waves in crystals with a simple spiral

[18] Buchelnikov V.D., Bychkov I.V., Shavrov V.G. Coupled magnetoelastic and electro‐ magnetic waves in uniaxial crystals having spiral magnetic structure. Journal of Magnetism and Magnetic Materials. 1993;118(1–2):169–174. DOI:

[19] Welsh G.H., Wynne K. Generation of ultrafast terahertz radiation pulses on metallic nanostructured surfaces. Optics Express. 2009;17(4):2470–2480. DOI: 10.1364/OE.

[20] Welsh G.H., Hunt N.T., Wynne K. Terahertz-pulse emission through laser excitation of surface plasmons in a metal grating. Physical Review Letters. 2007;98:026803. DOI:

[21] Gladun A.D., Leiman V.G., Arsenin A.V. On the mechanism of generation of terahertz electromagnetic radiation upon irradiation of a nanostructured metal surface by femtosecond laser pulses. Quantum Electronics. 2007;37(12):1166–1168. DOI: 10.1070/

[22] Nishitani J., Kozuki K., Nagashima T., Hangyo M. Terahertz radiation from coherent antiferromagnetic magnons excited by femtosecond laser pulses. Applied Physics

[23] Rungsawang R., Perez F., Oustinov D., Gómez J., Kolkovsky V., Karczewski G., et al. Terahertz radiation from magnetic excitations in diluted magnetic semiconductors. Physical Review Letters. 2013;110(17):177203. DOI: 10.1103/PhysRevLett.110.177203

[24] Bychkov I., Kuzmin D., Kalenov D., Kamantsev A., Koledov V., Kuchin D., et al. Electromagnetic waves generation in Ni2.14Mn0.81GaFe0.05 Heusler alloy at structur‐ al phase transition. Acta Physica Polonica A. 2015;127(2):588–590. DOI: 10.12693/

[25] Kimura T., Goto T., Shintani H., Ishizaka K., Arima T., Tokura Y. Magnetic control of ferroelectric polarization. Nature. 2003;426:55–58. DOI: 10.1038/nature02018

[26] Hur N., Park S., Sharma P.A., Ahn J.S., Guha S., Cheong S.W. Electric polarization reversal and memory in a multiferroic material induced by magnetic fields. Nature.

[27] Bychkov I.V., Kuzmin D.A., Lamekhov S.J., Shavrov V.G. Magnetoelectric susceptibil‐ ity tensor of multiferroic TbMnO3 with cycloidal antiferromagnetic structure in external

[28] Bychkov I.V., Kuzmin D.A., Shavrov V.G. Spectrum of coupled waves in orthorhombic multiferroics with cycloidal antiferromagnetic structure in external electric and

field. Journal of Applied Physics. 2013;113:17C726. DOI: 10.1063/1.4798820

magnetic structure. Low Temperature Physics (in Russian). 1988;14:606.

Spectroscopy (in Russian). 1981;56:37.

10.1016/0304-8853(93)90173-Y

10.1103/PhysRevLett.98.026803

QE2007v037n12ABEH013618

APhysPolA.127.588

Letters. 2010;96(22):221906. DOI: 10.1063/1.3436635

2004;429:392–395. DOI: 10.1038/nature02572

17.002470

178 Magnetic Materials


induced phase transition. Journal of Magnetism and Magnetic Materials. 2016. DOI: 10.1016/j.jmmm.2016.05.060

