Open access

The Eigen Theory of Electromagnetic Waves in Complex Media

Written By

Shaohua Guo

Submitted: 08 October 2010 Published: 05 July 2011

DOI: 10.5772/16313

From the Edited Volume

Behaviour of Electromagnetic Waves in Different Media and Structures

Edited by Ali Akdagli

Chapter metrics overview

3,219 Chapter Downloads

View Full Metrics

1. Introduction

Since J. C. Maxwell presented the electromagnetic field equations in 1873, the existence of electromagnetic waves has been verified in various medium (Kong, 1986; Monk, 2003). But except for Helmholtz’s equation of electromagnetic waves in isotropic media, the laws of propagation of electromagnetic waves in anisotropic media are not clear to us yet. For example, how many electromagnetic waves are there in anisotropic media? How fast can these electromagnetic waves propagate? Where are propagation direction and polarization direction of the electromagnetic waves? What are the space patterns of these waves? Although many research works were made in trying to deduce the equations of electromagnetic waves in anisotropic media based on the Maxwell’s equation (Yakhno, 2005, 2006; Cohen, 2002; Haba, 2004), the explicit equations of electromagnetic waves in anisotropic media could not be obtained because the dielectric permittivity matrix and magnetic permeability matrix were all included in these equations, so that only local behaviour of electromagnetic waves, for example, in a certain plane or along a certain direction, can be studied.

On the other hand, it is a natural fact that electric and magnetic fields interact with each other in classical electromagnetics. Therefore, even if most of material studies deal with the properties due to dielectric polarisation, magneitc materials are also capable of producing quite interesting electro-magnetic effects (Lindellm et al., 1994). From the bi-anisotropic point of view, magnetic materials can be treated as a subclass of magnetoelectric materials. The linear constitutive relations linking the electric and magnetic fields to the electric and magnetic displacements contain four dyadics, three of which have direct magnetic contents. The magnetoelectric coupling has both theoretical and practical significance in solid state physics and materials science. Though first predicted by Pierre Curie, magnetoeletric coupling was originally through to be forbidden because it violates time-reversal symmetry, until Laudau and Lifshitz (Landau & Lifshitz, 1960) pointed out that time reversal is not a symmetry operation in some magnetic crystal. Based on this argument, Dzyaloshinskii (Dzyaloshinskii, 1959) predicted that magnetoelectric effect should occur in antiferromagnetic crystal Cr2O3, which was verified experimentally by Astrov (Astrov, 1960). Since then the magnetoelectric coupling has been observed in single-phase materials where simultaneous electric and magnetic ordering coexists, and in two-phase composites where the participating phase are pizoelectric and piezomagnetic (Bracke & Van Vliet,1981; Van Run et al., 1974). Agyei and Birman (Agyei & Birman, 1990) carried out a detailed analysis of the linear magnetoelectric effect, which showed that the effect should occur not only in some magnetic but also in some electric crystals. Pradhan (Pradhan, 1993) showed that an electric charge placed in a magnetoelectric medium becomes a source of induced magnetic field with non-zero divergence of volume integral. Magnetoelectric effect in two-phase composites has been analyzed by Harshe et al. ( Harshe et al., 1993), Nan (Nan, 1994) and Benveniste (Benveniste, 1995). Broadband transducers based on magnetoelectric effect have also been developed (Bracke & Van Vliet, 1981). Although the development mentioned above, no great progress in the theories of electromagnetic waves in bi-anisotropic media because of the difficulties in deal with the bi-coupling in electric field and magnetic one of the Maxell’s equation and the bi-anisotropic constitutive equation by classical electromagnetic theory.

Recently there is a growing interest modeling and analysis of Maxwell’s equations (Lee & Madsen, 1990; Monk, 1992; Jin et al., 1999). However, most work is restricted to simple medium such as air in the free space. On the other hand, we notice that lossy and dispersive media are ubiquitous, for example human tissue, water, soil, snow, ice, plasma, optical fibers and radar-absorbing materials. Hence the study of how electromagnetic wave interacts with dispersive media becomes very important. Some concrete applications include geophysical probing and subsurface studied of the moon and other planets (Bui et al., 1991), High power and ultra-wide-band radar systems, in which it is necessary to model ultra-wide-band electromagnetic pulse propagation through plasmas (Dvorak & Dudley, 1995), ground penetrating radar detection of buried objects in soil media (liu & Fan, 1999). The Debye medium plays an important role in electromagnetic wave interactions with biological and water-based substances (Gandhi & Furse, 1997). Until 1990, some paper on modeling of wave propagation in dispersive media started making their appearance in computational electromagnetics community. However, the published papers on modeling of dispersive media are exclusively restricted to the finite-difference time-domain methods and the finite element methods (Li & Chen, 2006; Lu et al., 2004). To our best knowledge, there exist only few works in the literature, which studied the theoretical model for the Maxwell’s equation in the complex anisotropic dispersive media, and no explicit equations of electromagnetic waves in anisotropic dispersive media can be obtained due to the limitations of classical electromagnetic theory.

Chiral materials have been recently an interesting subject. In a chiral medium, an electric or magnetic excitation will produce simultaneously both electric and magnetic polarizations. On the other hand, the chiral medium is an object that cannot be brought into congruence with its mirror image by translation and rotation. Chirality is common in a variety of naturally occurring and man-made objects. From an operation point of view, chirality is introduced into the classical Maxwell equations by the Drude-Born-Fedorov relative constitutive relations in which the electric and magnetic fields are coupled via a new materials parameter (Lakhtakia, 1994; Lindell et al., 1994), the chirality parameter. These constitutive relations are chosen because they are symmetric under time reversality and duality transformations. In a homogeneous isotropic chiral medium the electromagnetic fields are composed of left-circularly polarized (LCP) and right- circularly polarized (RCP) components (Jaggard et al., 1979; Athanasiadis & Giotopoulos, 2003), which have different wave numbers and independent directions of propagation. Whenever an electromagnetic wave (LCP, RCP or a linear combination of them) is incident upon a chiral scatterer, then the scattered field is composed of both LCP and RCP components and therefore both LCP and RCP far-field patterns are derived. Hence, in the vector problem we need to specify two directions of propagation and two polarizations. In recent years, chiral materials have been increasingly studied and there is a growing literature covering both their applications and the theoretical investigation of their properties. It will be noticed that the works dealing with wave phenomena in chiral materials have been mainly concerned with the study of time-harmonic waves which lead to frequency domain studies (Lakhtakia et al., 1989; Athanasiadis et al., 2003).

In this chapter, the idea of standard spaces is used to deal with the Maxwell’s electromagnetic equation (Guo, 2009, 2009, 2010, 2010, 2010). By this method, the classical Maxwell’s equation under the geometric presentation can be transformed into the eigen Maxwell’s equation under the physical presentation. The former is in the form of vector and the latter is in the form of scalar. Through inducing the modal constitutive equations of complex media, such as anisotropic media, bi-anisotropic media, lossy media, dissipative media, and chiral media, a set of modal equations of electromagnetic waves for all of those media are obtained, each of which shows the existence of electromagnetic sub-waves, meanwhile its propagation velocity, propagation direction, polarization direction and space pattern can be completely determined by the modal equations.This chapter will make introductions of the eigen theory to reader in details. Several novel theoretical results were discussed in the different parts of this chapter.

Advertisement

2. Standard spaces of electromagnetic media

In anisotropic electromagnetic media, the dielectric permittivity and magnetic permeability are tensors instead of scalars. The constitutive relations are expressed as follows

D = ε E ,     B = μ H E1

Rewriting Eq.(1) in form of scalar, we have

D i = ε i j E j ,    B i = μ i j H j E2

where the dielectric permittivity matrix ε and the magnetic permeability matrix μ are usually symmetric ones, and the elements of the matrixes have a close relationship with the selection of reference coordinate. Suppose that if the reference coordinates is selected along principal axis of electrically or magnetically anisotropic media, the elements at non-diagonal of these matrixes turn to be zero. Therefore, equations (1) and (2) are called the constitutive equations of electromagnetic media under the geometric presentation. Now we intend to get rid of effects of geometric coordinate on the constitutive equations, and establish a set of coordinate-independent constitutive equations of electromagnetic media under physical presentation. For this purpose, we solve the following problems of eigen-value of matrixes.

( ε λ I ) ϕ = 0 ,     ( μ γ I ) φ = 0 E3

where λ i ( i = 1 , 2 , 3 ) and γ i ( i = 1 , 2 , 3 ) are respectively eigen dielectric permittivity and eigen magnetic permeability, which are constants of coordinate-independent. ϕ i ( i = 1 , 2 , 3 ) and φ i ( i = 1 , 2 , 3 ) are respectively eigen electric vector and eigen magnetic vector, which show the electrically principal direction and magnetically principal direction of anisotropic media, and are all coordinate-dependent. We call these vectors as standard spaces. Thus, the matrix of dielectric permittivity and magnetic permeability can be spectrally decomposed as follows

ε = Φ Λ Φ Τ ,     μ = Ψ Π Ψ Τ E4

where Λ = d i a g [ λ 1 , λ 2 , λ 3 ] and Π = d i a g [ γ 1 , γ 2 , γ 3 ] are the matrix of eigen dielectric permittivity and eigen magnetic permeability, respectively. Φ = { ϕ 1 , ϕ 2 , ϕ 3 } and Ψ = { φ 1 , φ 2 , φ 3 } are respectively the modal matrix of electric media and magnetic media, which are both orthogonal and positive definite matrixes, and satisfy Φ T Φ = I , Ψ T Ψ = I .

Projecting the electromagnetic physical qualities of the geometric presentation, such as the electric field intensity vector E , magnetic field intensity vector H , magnetic flux density vector B and electric displacement vector D into the standard spaces of the physical presentation, we get

D * = Φ Τ D , E * = Φ Τ E E5
B * = Ψ Τ B , H * = Ψ Τ H E6

Rewriting Eqs.(5) and (6) in the form of scalar, we have

D i * = ϕ i Τ D i = 1 , 2 , 3 , E i * = ϕ i Τ E i = 1 , 2 , 3 E7
B i * = φ i Τ B i = 1 , 2 , 3 , H i * = φ i Τ H i = 1 , 2 , 3 E8

These are the electromagnetic physical qualities under the physical presentation.

Substituting Eq. (4) into Eq. (1) respectively, and using Eqs.(5) and (6) yield

D i * = λ i E i * i = 1 , 2 , 3 E9
B i * = γ i H i * i = 1 , 2 , 3 E10

The above equations are just the modal constitutive equations in the form of scalar.

Advertisement

3. Eigen expression of Maxwell’s equation

The classical Maxwell’s equations in passive region can be written as

× H = t D , × E = t B E11

Now we rewrite the equations in the form of matrix as follows

[ 0 z y z 0 x y x 0 ] { H 1 H 2 H 3 } = t { D 1 D 2 D 3 } E12

or

[ Δ ] { H } = t { D } E13
[ 0 z y z 0 x y x 0 ] { E 1 E 2 E 3 } = t { B 1 B 2 B 3 } E14

or

[ Δ ] { E } = t { B } E15

where [ Δ ] is defined as the matrix of electric and magnetic operators.

Substituting Eq. (1) into Eqs. (13) and (15) respectively, we have

[ Δ ] { H } = t [ ε ] { E } E16
[ Δ ] { E } = t [ μ ] { H } E17

Substituting Eq. (16) into (17) or Eq. (17) into (16), yield

[ ] { H } = t 2 [ μ ] [ ε ] { H } E18
[ ] { E } = t 2 [ μ ] [ ε ] { E } E19

where [ ] = [ Δ ] [ Δ ] is defined as the matrix of electromagnetic operators as follows

[ ] = [ ( z 2 + y 2 ) x y 2 x z 2 y x 2 ( x 2 + z 2 ) y z 2 z x 2 z y 2 ( x 2 + y 2 ) ] E20

In another way, substituting Eqs. (5) and (6) into Eqs. (13) and (15), respectively, we have

[ Δ ] [ Φ ] { E * } = t [ Ψ ] { B * } E21
[ Δ ] [ Ψ ] { H * } = t [ Φ ] { D * } E22

Rewriting the above in indicial notation, we get

{ Δ i * } E i * = t { φ i } B i * i = 1 , 2 , 3 E23
{ Δ i * } H i * = t { ϕ i } D i * i = 1 , 2 , 3 E24

where, Δ i * is the electromagnetic intensity operator, and i th row of [ Δ * ] = [ Δ ] [ Φ ] .

Advertisement

4. Electromagnetic waves in anisotropic media

4.1. Electrically anisotropic media

In anisotropic dielectrics, the dielectric permittivity is a tensor, while the magnetic permeability is a scalar. So Eqs. (18) and (19) can be written as follows

[ ] { H } = t 2 μ 0 [ ε ] { H } E25
[ ] { E } = t 2 μ 0 [ ε ] { E } E26

Substituting Eqs. (4) - (6) into Eqs. (25) and (26), we have

[ * ] { H * } = t 2 μ 0 [ Λ ] { H * } E27
[ * ] { E * } = t 2 μ 0 [ Λ ] { E * } E28

where [ * ] = [ Φ ] T [ ] [ Φ ] is defined as the eigen matrix of electromagnetic operators under the standard spaces. We can note from Appendix A that it is a diagonal matrix. Thus Eqs. (27) and (28) can be uncoupled in the form of scalar

i * H i * + μ 0 λ i t 2 H i * = 0 i = 1 , 2 , 3 E29
i * E i * + μ 0 λ i t 2 E i * = 0 i = 1 , 2 , 3 E30

Eqs.(29) and (30) are the modal equations of electromagnetic waves in anisotropic dielectrics.

4.2. Magnetically anisotropic media

In anisotropic magnetics, the magnetic permeability is a tensor, while the dielectric permittivity is a scalar. So Eqs. (18) and (19) can be written as follows

[ ] { H } = t 2 ε 0 [ μ ] { H } E31
[ ] { E } = t 2 ε 0 [ μ ] { E } E32

Substituting Eqs. (4) - (6) into Eqs. (31) and (32), we have

[ * ] { H * } = t 2 ε 0 [ Π ] { H * } E33
[ * ] { E * } = t 2 ε 0 [ Π ] { E * } E34

where [ * ] = [ Ψ ] T [ ] [ Ψ ] is defined as the eigen matrix of electromagnetic operators under the standard spaces. We can also note from Appendix A that it is a diagonal matrix. Thus Eqs. (33) and (34) can be uncoupled in the form of scalar

i * H i * + ε 0 γ i t 2 H i * = 0 i = 1 , 2 , 3 E35
i * E i * + ε 0 γ i t 2 E i * = 0 i = 1 , 2 , 3 E36

Eqs.(35) and (36) are the modal equations of electromagnetic waves in anisotropic magnetics.

Advertisement

5. Electromagnetic waves in bi-anisotropic media

5.1. Bi-anisotropic constitutive equations

The constitutive equations of bi-anisotropic media are the following (Lindellm & Sihvola, 1994; Laudau & Lifshitz, 1960)

D = ε E + ξ H E37
B = ξ E + μ H E38

where ξ is the matrix of magneto-electric parameter, and a symmetric one.

Substituting Eqs. (5) and (6) into Eqs. (37) and (38), respectively, and multiplying them with the transpose of modal matrix in the left, we have

Φ T D = Φ T ε Φ E * + Φ T ξ Ψ H * E39
Ψ Τ B = Ψ Τ ξ Φ E * + Ψ Τ μ Ψ H * E40

Let G = Φ T ξ Ψ = Ψ T ξ Φ , that is a coupled magneto-electric matrix, and using Eq. (4), we have

D * = Λ E * + G H * E41
B * = G E * + Π H * E42

Rewriting the above in indicial notation, we get

D i * = λ i E i * + g i j H j * i = 1 , 2 , 3 j = 1 , 2 , 3 E43
B i * = γ i H i * + g i j E j * i = 1 , 2 , 3 j = 1 , 2 , 3 E44

Eqs. (43) and (44) are just the modal constitutive equations for bi-anisotropic media.

5.2. Eigen equations of electromagnetic waves in bi-anisotropic media

Substituting Eqs. (43) and (44) into Eqs. (23) and (24), respectively, we have

{ Δ i * } E i * = t { φ i } ( γ i H i * + g i j E j * ) E45
{ Δ i * } H i * = t { ϕ i } ( λ i E i * + g i j H j * ) E46

From them, we can get

( { Δ i * } t { ϕ i } g i j δ i j ) T ( { Δ i * } + t { φ i } g i j δ i j ) E i * = t 2 { φ i } { ϕ i } T λ i γ i E i * E47
( { Δ i * } + t { φ i } g i j δ i j ) T ( { Δ i * } t { ϕ i } g i j δ i j ) H i * = t 2 { ϕ i } { φ i } T γ i λ i H i * E48

The above can also be written as the standard form of waves

i * E i * + t { Δ i * } T ( { φ } { ϕ } ) g i i E i * + t 2 { ϕ } T { φ } ( λ i γ i g i i 2 ) E i * = 0 i = 1 , 2 , 3 E49
i * H i * + t { Δ i * } T ( { φ } { ϕ } ) g i i H i * + t 2 { ϕ } T { φ } ( λ i γ i g i i 2 ) H i * = 0 i = 1 , 2 , 3 E50

where, i * = { Δ i * } T { Δ i * } is the electromagnetic operator. Eqs.(49) and (50) are just equations of electric field and magnetic field for bi-anisotropic media.

5.3. Applications

5.3.1. Bi-isotropic media

The constitutive equations of bi-isotropic media are the following

D = [ ε 0 0 0 ε 0 0 0 ε ] E + [ ξ 0 0 0 ξ 0 0 0 ξ ] H E51
B = [ ξ 0 0 0 ξ 0 0 0 ξ ] E + [ μ 0 0 0 μ 0 0 0 μ ] H E52

The eigen values and eigen vectors of those matrix are the following

Λ = d i a g [ ε , ε , ε ] , Π = d i a g [ μ , μ , μ ] E53
Φ = Ψ = [ 1 0 0 0 1 0 0 0 1 ] E54

We can see from the above equations that there is only one eigen-space in isotropic medium, which is a triple-degenerate one, and the space structure is the following

W = W 1 ( 3 ) [ ϕ 1 , ϕ 2 , ϕ 3 ] E55
ϕ 1 * = φ 1 * = 3 3 [ 1 , 1 , 1 ] T E56

Then the eigen-qualities and eigen-operators of bi-isotropic medium are respectively shown as belows

E 1 * = ϕ 1 * T E = 3 3 ( E 1 + E 2 + E 3 ) E57
1 * = ( x 2 + y 2 + z 2 ) , g 11 = ξ E58

So, the equation of electromagnetic wave in bi-isotropic medium becomes

( x 2 + y 2 + z 2 ) E 1 * = ( μ ε ξ 2 ) t 2 E 1 * E59

the velocity of electromagnetic wave is

c ( 1 ) = 1 μ ε ξ 2 E60

5.3.2. Dzyaloshinskii’s bi-anisotropic media

Dzyaloshinskii’s constitutive equations of bi-anisotropic media are the following

D = [ ε 0 0 0 ε 0 0 0 ε z ] E + [ ξ 0 0 0 ξ 0 0 0 ξ z ] H E61
B = [ ξ 0 0 0 ξ 0 0 0 ξ z ] E + [ μ 0 0 0 μ 0 0 0 μ z ] H E62

The eigen values and eigen vectors of those matrix are the following

Λ = d i a g [ ε , ε , ε z ] , Π = d i a g [ μ , μ , μ z ] E63
Φ = Ψ = [ 1 0 0 0 1 0 0 0 1 ] E64

We can see from the above equations that there are two eigen-spaces in Dzyaloshinskii’s bi-anisotropic medium, in which one is a binary-degenerate one, the space structure is the following

W = W 1 ( 2 ) [ ϕ 1 , ϕ 2 ] W 2 1 [ ϕ 3 ] E65

Then the eigen-qualities and eigen-operators of Dzyaloshinskii’s bi-anisotropic medium are respectively shown as belows

E 2 * = ϕ 2 T E = E 3 , | E 1 * | = ( E ϕ 2 T E 2 * ) T ( E ϕ 2 T E 2 * ) = E 1 2 + E 2 2 E66
1 * = ( x 2 + y 2 + 2 z 2 2 x y 2 ) , 2 * = ( x 2 + y 2 ) , g 11 = ξ , g 22 = ξ z E67

So, the equations of electromagnetic wave in Dzyaloshinskii’s bi-anisotropic medium become

( x 2 + y 2 + 2 z 2 2 x y 2 ) E 1 2 + E 2 2 = ( μ ε ξ 2 ) t 2 E 1 2 + E 2 2 E68
( x 2 + y 2 ) E 3 = ( μ z ε z ξ z 2 ) t 2 E 3 E69

the velocities of electromagnetic wave are

c ( 1 ) = 1 μ ε ξ 2 E70
c ( 2 ) = 1 μ z ε z ξ z 2 E71

It is seen both from bi-isotropic media and Dzyaloshinskii’s bi-anisotropic medium that the electromagnetic waves in bi-anisotropic medium will go faster duo to the bi-coupling between electric field and magnetic one.

Advertisement

6. Electromagnetic waves in lossy media

6.1. The constitutive equation of lossy media

The constitutive equation of lossy media is the following

D = ε E + t σ d E d τ d τ E72

It is equivalent to the following differential constitutive equation

D ˙ = ε E ˙ + σ E E73

Let

D e = ε E , D ˙ d = σ E E74

Eq.(73) can be written as

D ˙ = D ˙ e + D ˙ d E75

or

t { D } = ( [ ε ] t + [ σ ] ) { E } E76

Using Eq.(5), the above becomes

t { D * } = ( [ Φ ] T [ ε ] [ Φ ] t + [ Φ ] T [ σ ] [ Φ ] ) { E * } E77

According to Appendix B and Eq.(77), we have

t { D * } = ( [ Λ ] t + [ Γ ] ) { E * } E78

Rewriting the above in indicial notation, we get

t D i * = ( λ i t + η i ) E i * E79

Eq.(79) is just the modal constitutive equations for lossy media.

6.2. Eigen equations of electromagnetic waves in lossy media

Substituting Eqs. (10) and (79) into Eqs. (23) and (24), respectively, we have

{ Δ i * } E i * = t { φ i } γ i H i * i = 1 , 2 , 3 E80
{ Δ i * } H i * = { ϕ i } ( λ i t + η i ) E i * i = 1 , 2 , 3 E81

From them, we can get

i * E i * + t t ξ i γ i λ i E i * + t ξ i γ i η i E i * = 0 i = 1 , 2 , 3 E82
i * H i * + t t ξ i γ i λ i H i * + t ξ i γ i η i H i * = 0 i = 1 , 2 , 3 E83

where ξ i = { ϕ i * } T { φ i * } . Eqs.(82) and (83) are just equations of electric field and magnetic field for bi-anisotropic media.

6.3. Applications

In this section, we discuss the propagation laws of electromagnetic waves in an isotropic lossy medium. The material tensors in Eqs.(1) and (72) are represented by the following matrices

ε = [ ε 1 1 0 0 0 ε 1 1 0 0 0 ε 1 1 ] , μ = [ μ 11 0 0 0 μ 11 0 0 0 μ 11 ] , σ = [ σ 11 0 0 0 σ 11 0 0 0 σ 11 ] E84

The eigen values and eigen vectors of those matrix are the following

Λ = d i a g [ ε 11 , ε 11 , ε 11 ] , Π = d i a g [ μ 11 , μ 11 , μ 11 ] , Γ = d i a g [ σ 11 , σ 11 , σ 11 ] E85
Φ = Ψ = Θ = [ 1 0 0 0 1 0 0 0 1 ] E86

We can see from the above equations that there is only one eigen-space in an isotropic lossy medium, which is a triple-degenerate one, and the space structure is the following.

W m a g = W 1 ( 3 ) [ ϕ 1 , ϕ 2 , ϕ 3 ] , W e l e = W 1 ( 3 ) [ φ 1 , φ 2 , φ 3 ] E87

where, ϕ 1 * = 3 3 { 1 , 1 , 1 } T , φ 1 * = 3 3 { 1 , 1 , 1 } T , ξ 1 = 1 .

Then the eigen-qualities and eigen-operators of an isotropic lossy media are respectively shown as follows

E 1 * = 3 3 ( E 1 + E 2 + E 3 ) E88
H 1 * = 3 3 ( H 1 + H 2 + H 3 ) E89
1 * = 1 3 [ ( x 2 + y 2 + z 2 ) ] E90

So, the equation of electromagnetic wave in lossy media becomes

( x 2 + y 2 + z 2 ) E 1 * = 1 c 2 t 2 E 1 * + 1 τ 2 t E 1 * E91

Rewriting it in the component form, we have

( x 2 + y 2 + z 2 ) E 1 = 1 c 2 t 2 E 1 + 1 τ 2 t E 1 E92
( x 2 + y 2 + z 2 ) E 2 = 1 c 2 t 2 E 2 + 1 τ 2 t E 2 E93
( x 2 + y 2 + z 2 ) E 3 = 1 c 2 t 2 E 3 + 1 τ 2 t E 3 E94

where, c is the velocity of electromagnetic wave, τ is the lossy coefficient of electromagnetic wave

c = 1 μ 11 ε 11 ,    τ = 1 μ 11 σ 11 E95

Now, we discuss the the propagation laws of a plane electromagnetic wave in x-axis. In this time, Eq. (92) becomes

x 2 E 1 = 1 c 2 t 2 E 1 + 1 τ 2 t E 1 E96

Let the solution of Eq. (96) is as follows

E 1 = A exp [ i ( k x ω t ) ] E97

Substituting the above into Eq. (96), we have

k 2 = ω 2 c 2 + i ω τ 2 E98

From Eq.(96), we can get

k ¯ = k 1 + i k 2 E99

where

k 1 = ω c [ 1 + ( 1 + c 4 ω 2 τ 4 ) 1 2 2 ] 1 2 , k 2 = ω c [ 1 + ( 1 + c 4 ω 2 τ 4 ) 1 2 2 ] 1 2 E100

Then, the solutions of electromagnetic waves are the following

E 1 = A e k 2 x e i ( k 1 x ω t ) = A ¯ e i ( k 1 x ω t ) E101

It is an attenuated sub-waves.

Advertisement

7. Electromagnetic waves in dispersive media

7.1. The constitutive equation of dispersive media

The general constitutive equations of dispersive media are the following

D = ε E + ε 1 E ˙ + ε 2 E ¨ + E102
B = μ H + μ 1 H ˙ + μ 2 H ¨ + E103

where ε i , ( i = 1 , 2 , ) and μ i , ( i = 1 , 2 , ) are the higher order dielectric permittivity matrix and the magnetic permeability matrix respectively, and all symmtric ones.

Substituting Eqs. (5) and (6) into Eqs. (101) and (102), respectively, and multiplying them with the transpose of modal matrix in the left, we have

D * = Φ T ε Φ E * + Φ T ε 1 Φ E ˙ + Φ T ε 2 Φ E ¨ + E104
B * = Ψ T μ Ψ H + Ψ T μ 1 Ψ H ˙ + Ψ T μ 2 Ψ H ¨ + E105

It can be proved that there exist same standard spaces for various order electric and magnetic fields in the condition close to the thermodynamic equilibrium. Then, we have

D i * = λ i E i * + λ i ( 1 ) E ˙ i * + λ i ( 2 ) E ¨ i * + E106
B i * = γ i H i * + γ i ( 1 ) H ˙ i * + γ i ( 2 ) H ¨ i * + E107

Eqs. (105) and (106) are just the modal constitutive equations for the general dispersive media.

7.2. Eigen equations of electromagnetic waves in dispersive media

Substituting Eqs. (105) and (106) into Eqs. (23) and (24), respectively, we have

{ Δ i * } E i * = t { φ i } ( γ i H i * + γ i ( 1 ) H ˙ i * + γ i ( 2 ) H ¨ i * + ) E108
{ Δ i * } H i * = t { ϕ i } ( λ i E i * + λ i ( 1 ) E ˙ i * + λ i ( 2 ) E ¨ i * + ) E109

From them, we can get

1 ξ i i * E i * + γ i λ i t t E i * + ( γ i λ i ( 1 ) + γ i ( 1 ) λ i ) t t t E i * + ( γ i λ i ( 2 ) + γ i ( 1 ) λ i ( 1 ) + γ i ( 2 ) λ i ) t t t t E i * + = 0 E110
1 ξ i i * H i * + γ i λ i t t H i * + ( γ i λ i ( 1 ) + γ i ( 1 ) λ i ) t t t H i * + ( γ i λ i ( 2 ) + γ i ( 1 ) λ i ( 1 ) + γ i ( 2 ) λ i ) t t t t H i * + = 0 E111

Eqs.(109) and (110) are just equations of electric field and magnetic field for general dispersive media.

7.3. Applications

In this section, we discuss the propagation laws of electromagnetic waves in an one-order dispersive medium. The material tensors in Eqs.(101) and (102) are represented by the following matrices

ε = [ ε 1 1 0 0 0 ε 1 1 0 0 0 ε 1 1 ] , μ = [ μ 11 0 0 0 μ 11 0 0 0 μ 11 ] , ε 1 = [ ε 1 1 0 0 0 ε 1 1 0 0 0 ε 1 1 ] , μ 1 = [ μ 11 0 0 0 μ 11 0 0 0 μ 11 ] E112

The eigen values and eigen vectors of those matrix are the following

Λ = d i a g [ ε 11 , ε 11 , ε 11 ] , Λ 1 = d i a g [ ε 11 , ε 11 , ε 11 ] E113
Π = d i a g [ μ 11 , μ 11 , μ 11 ] , Π 1 = d i a g [ μ 11 , μ 11 , μ 11 ] E114
Φ = Ψ = [ 1 0 0 0 1 0 0 0 1 ] E115

We can see from the above equations that there is only one eigen-space in isotropic one-order dispersive medium, which is a triple-degenerate one, and the space structure is the following

W m a g = W 1 ( 3 ) [ ϕ 1 , ϕ 2 , ϕ 3 ] , W e l e = W 1 ( 3 ) [ φ 1 , φ 2 , φ 3 ] E116

where, ϕ 1 * = 3 3 { 1 , 1 , 1 } T , φ 1 * = 3 3 { 1 , 1 , 1 } T , ξ 1 = 1 . Thus the eigen-qualities and eigen-operators of isotropic one-order dispersive medium are known as same as Eqs. (88) – (90). The equations of electromagnetic wave in one-order dispersive medium become

( x 2 + y 2 + z 2 ) E 1 * = 1 c 2 t t E 1 * + ( μ 11 ε 11 ( 1 ) + μ 11 ( 1 ) ε 11 ) t t t E 1 * E117
( x 2 + y 2 + z 2 ) H 1 * = 1 c 2 t t H 1 * + ( μ 11 ε 11 ( 1 ) + μ 11 ( 1 ) ε 11 ) t t t H 1 * E118

in which c is the velocity of electromagnetic wave

c = 1 μ 11 ε 11 E119

Now, we discuss the propagation laws of a plane electromagnetic wave in x-axis. In this time, Eq.(116) becomes

x 2 E 1 = 1 c 2 t t E 1 + ( μ 11 ε 11 ( 1 ) + μ 11 ( 1 ) ε 11 ) t t t E 1 E120

Let

E 1 = A exp [ i ( k x ω t ) ] E121

Substituting the above into Eq.(119), we have

k 2 = ω 2 c 2 i ω 3 ( μ 11 ε 11 ( 1 ) + μ 11 ( 1 ) ε 11 ) E122

From the above, we can get

k ¯ = k 1 + i k 2 E123

where

k 1 = ω c [ 1 + ( 1 + c 4 ( μ 11 ε 11 ( 1 ) + μ 11 ( 1 ) ε 11 ) 2 ω 2 ) 1 2 2 ] 1 2 , k 2 = ω c [ 1 + ( 1 + c 4 ( μ 11 ε 11 ( 1 ) + μ 11 ( 1 ) ε 11 ) 2 ω 2 ) 1 2 2 ] 1 2 E124

Then, the solutions of electromagnetic waves are

E 1 = A e k 2 x e i ( k 1 x ω t ) = A ¯ e i ( k 1 x ω t ) E125

It is an attenuated sub-waves.

Advertisement

8. Electromagnetic waves in chiral media

8.1. The constitutive equation of chiral media

The constitutive equations of chiral media are the following

D = ε E - χ t H E126
B = χ t E + μ H E127

where χ is the matrix of chirality parameter, and a symmtric one.

Substituting Eqs. (5) and (6) into Eqs. (124) and (125), respectively, and multiplying them with the transpose of modal matrix in the left, we have

D * = Φ T ε Φ E * - Φ T χ Ψ t H * E128
B * = Ψ T χ Φ t E * + Ψ T μ Ψ H * E129
Let Γ = Ψ T χ Φ , that is a coupled chiral matrix, and using Eq. (4), we have
D * = Λ E * - Γ T t H * E130
B * = Γ T t E * + Π H * E131

For most chiral, Γ = d i a g [ ς 1 , ς 2 , ς 3 ] . Then we have

D i * = λ i E i * ς i t H i * E132
B i * = ς i t E i * + γ i H i * E133

Eqs.(130) and (131) are just the modal constitutive equations for anisotropic chiral media.

8.2. Eigen equations of electromagnetic waves in chiral media

Substituting Eqs. (130) and (131) into Eqs. (23) and (24), respectively, we have

{ Δ i * } E i * = t { φ i } ( ς i t E i * + γ i H i * ) i = 1 , 2 , 3 E134
{ Δ i * } H i * = t { ϕ i } ( λ i E i * ς i t H i * ) i = 1 , 2 , 3 E135

From them, we can get

i * E i * + ξ i ς i 2 t t t t E i * + ( 2 ς i i * + ξ i λ i γ i ) t t E i * = 0 i = 1 , 2 , 3 E136
i * H i * + ξ i ς i 2 t t t t H i * + ( 2 ς i i * + ξ i λ i γ i ) t t H i * = 0 i = 1 , 2 , 3 E137

where ξ i = { ϕ i } T { φ i } = 1 , i = { Δ i * } T { φ i } . Eqs.(134) and (135) are just equations of electric field and magnetic field for chiral media.

8.3. Applications

In this section, we discuss the propagation laws of electromagnetic waves in an isotropic chiral medium. The material tensors in Eqs.(124) and (125) are represented by the following matrices

ε = [ ε 1 1 0 0 0 ε 1 1 0 0 0 ε 1 1 ] , μ = [ μ 11 0 0 0 μ 11 0 0 0 μ 11 ] , χ = [ χ 11 0 0 0 χ 11 0 0 0 χ 11 ] E138

The eigen values and eigen vectors of those matrix are the following

Λ = d i a g [ ε 11 , ε 11 , ε 11 ] , Π = d i a g [ μ 11 , μ 11 , μ 11 ] , Γ = d i a g [ χ 11 , χ 11 , χ 11 ] E139
Φ = Ψ = [ 1 0 0 0 1 0 0 0 1 ] E140

We can see from the above equations that there is only one eigen-space in isotropic medium, which is a triple-degenerate one, and the space structure is the following

W m a g = W 1 ( 3 ) [ ϕ 1 , ϕ 2 , ϕ 3 ] , W e l e = W 1 ( 3 ) [ φ 1 , φ 2 , φ 3 ] E141

where, ϕ 1 * = 3 3 { 1 , 1 , 1 } T , φ 1 * = 3 3 { 1 , 1 , 1 } T , ξ 1 = 1 .

Then the eigen-qualities and eigen-operators of isotropic chiral medium are respectively shown as follows

E 1 * = 3 3 ( E 1 + E 2 + E 3 ) , H 1 * = 3 3 ( H 1 + H 2 + H 3 ) E142
1 * = 1 3 [ ( x 2 + y 2 + z 2 ) ] , 1 * = 3 3 ( x + y + z ) E143

So, the equations of electromagnetic wave in isotropic chiral medium become

( x 2 + y 2 + z 2 ) E 1 * = χ 11 2 t t t t E 1 * + [ 2 3 χ 11 ( x + y + z ) + 1 c 2 ] t t E 1 * E144
( x 2 + y 2 + z 2 ) H 1 * = χ 11 2 t t t t H 1 * + [ 2 3 χ 11 ( x + y + z ) + 1 c 2 ] t t H 1 * E145

where, c is the velocity of electromagnetic wave

c = 1 μ 11 ε 11 E146

Now, we discuss the propagation laws of a plane electromagnetic wave in x-axis. In this time, Eq.(142) becomes

x 2 E 1 = χ 11 2 t t t t E 1 + [ 2 3 χ 11 x + 1 c 2 ] t t E 1 E147
Let
E 1 = A exp [ i ( k x ω t ) ] E148

Substituting the above into Eq.(145), we have

k 2 = ( 1 c 2 + i 2 3 χ 11 k ) ω 2 χ 11 2 ω 4 E149
or
k 2 i 2 3 χ 11 ω 2 k + ( χ 11 2 ω 2 μ 11 ε 11 ) ω 2 = 0 E150
  1. when

    ω 2 μ 11 ε 11 χ 11 2 0 E151

    By Eq.(148), we have

    k 1 = k 1 + i k 1 = i χ 11 ω 2 [ 3 + ( x + i y ) ] E152
    k 2 = k 2 + i k 2 = i χ 11 ω 2 [ 3 ( x + i y ) ] E153

    From them, we can get

    k 1 = χ 11 ω 2 y , k 1 = ( 3 + x ) χ 11 ω 2 E154
    k 2 = χ 11 ω 2 y , k 2 = ( 3 x ) χ 11 ω 2 E155

    where, x = 2 2 3 + [ 9 + ( μ 11 ε 11 χ 11 2 ω 2 ) 2 ω 4 ] 1 2 ,

    y = 2 2 3 + [ 9 + ( μ 11 ε 11 χ 11 2 ω 2 ) 2 ω 4 ] 1 2 E156
    .

    Then, the solution of electromagnetic waves is the following

    E 1 = A 1 e k 1 x e i ( k 1 x ω t ) + A 2 e k 2 x e i ( k 2 x ω t ) E157

    It is composed of two attenuated sub-waves.

  2. when

    ω 2 μ 11 ε 11 χ 11 2 0 E158

    By Eq.(148), we have

    k 1 = k 1 + i k 1 = i χ 11 ω 2 [ 3 + 3 + 1 ω 2 ( ω 2 μ 11 ε 11 χ 11 2 ) ] E159
    k 2 = k 2 + i k 2 = i χ 11 ω 2 [ 3 3 + 1 ω 2 ( ω 2 μ 11 ε 11 χ 11 2 ) ] E160

    where

    k 1 = 0 , k 1 = χ 11 ω 2 [ 3 + 3 + 1 ω 2 ( ω 2 μ 11 ε 11 χ 11 2 ) ] E161
    k 2 = χ 11 ω 2 [ 3 + 1 ω 2 ( ω 2 μ 11 ε 11 χ 11 2 ) 3 ] , k 2 = 0 E162

    Then, the solution of electromagnetic waves is the following

    E 1 = A 1 e k 1 x e i ω t + A 2 e i ( k 2 x ω t ) E163

    It is seen that there only exists an electromagnetic sub-wave in opposite direction.

  3. when

    ω 2 μ 11 ε 11 χ 11 2 = 0 E164

    By Eq.(148), we have

    k = k + i k = i 2 3 χ 11 ω 2 E165

    where

    k = 0 , k = 2 3 χ 11 ω 2 E166

    Then, the solution of electromagnetic waves is the following

    E = A e k x e i ω t E167

    No electromagnetic sub-waves exist now.

Advertisement

9. Conclusion

In this chapter, we construct the standard spaces under the physical presentation by solving the eigen-value problem of the matrixes of dielectric permittivity and magnetic permeability, in which we get the eigen dielectric permittivity and eigen magnetic permeability, and the corresponding eigen vectors. The former are coordinate-independent and the latter are coordinate-dependent. Because the eigen vectors show the principal directions of electromagnetic media, they can be used as the standard spaces. Based on the spaces, we get the modal equations of electromagnetic waves for anisotropic media, bi-anisotropic media, dispersive medium and chiral medium, respectively, by converting the classical Maxwell’s vector equation to the eigen Maxwell’s scalar equation, each of which shows the existence of an electromagnetic sub-wave, and its propagation velocity, propagation direction, polarization direction and space pattern are completely determined in the equations. Several novel results are obtained for anisotropic media. For example, there is only one kind of electromagnetic wave in isotropic crystal, which is identical with the classical result; there are two kinds of electromagnetic waves in uniaxial crystal; three kinds of electromagnetic waves in biaxial crystal and three kinds of distorted electromagnetic waves in monoclinic crystal. Also for bi-anisotropic media, there exist two electromagnetic waves in Dzyaloshinskii’s bi-anisotropic media, and the electromagnetic waves in bi-anisotropic medium will go faster duo to the bi-coupling between electric field and magnetic one. For isotropic dispersive medium, the electromagnetic wave is an attenuated sub-waves. And for chiral medium, there exist different propagating states of electromagnetivc waves in different frequency band, for example, in low frequency band, the electromagnetic waves are composed of two attenuated sub-waves, in high frequency band, there only exists an electromagnetic sub-wave in opposite direction, and in the critical point, no electromagnetic can propagate. All of these new theoretical results need to be proved by experiments in the future.

Advertisement

10. Appendix A: Proof of the eigenmode of electromagnetic operator matrix

The Maxwell’s equation of anisotropic dielectrics is the following

[ ] { H } = t 2 μ 0 [ ε ] { H } E168

Using the representation transform relationship Eq. (6), we have

[ ] [ Φ ] { H * } = t 2 μ 0 [ ε ] [ Φ ] { H * } E169

Substituting the spectral decomposition matrix of dielectric permittivity Eq. (4) into above, we have

[ ] [ Φ ] { H * } = t 2 μ 0 [ Φ ] [ Λ ] { H * } E170

Comparing the both sides of above equation, we can get

[ ] [ Φ ] = t 2 μ 0 [ Φ ] [ Λ ] E171

Multiplying the both sides of above with the transpose of modal matrix in the left, we have

[ Φ ] T [ ] [ Φ ] = t 2 μ 0 [ Λ ] E172

It is seen that the right side above is a diagonal matrix, which shows that the electromagnetic operators matrix can also be spectrally decomposed in standard spaces, then we get

[ * ] = t 2 μ 0 [ Λ ] E173

Rewriting above in the form of scalar, we have

i * = t 2 μ 0 λ i E174

11. Appendix B: Spectrally decomposition of lossy matrix

The Helmholtz’s free energy of electromagnetic system with lossy property is the following

ψ = ψ ( D , B , D d ) = 1 2 ( D D d ) ε 1 ( D D d ) + 1 2 B μ 1 B E175

Differentiating the above with lossy variable, and using Eq.(1), we have

R = ψ D d = ( D D d ) ε 1 = E E176

According to the Onsager’s principle, for the process of closing to equilibrium, the rate of lossy variable is proportion to the driving force, that is

ψ D i d + β i j d D j d d t = 0 E177

where, β i j = β j i is the general friction coefficient. Rewritting the above in matrix form, we have

{ R } = [ B ] { D ˙ d } E178

Projecting the lossy electric displacement vector D into the standard spaces of the physical presentation, we get

{ D d } = a { φ } E179

Using Eq.(B2), we have

{ R } = ω a ˙ { φ } E180

Substituting Eqs. (B5) and (B6) into Eq. (B4), the condition of non-zero solution to α is the following

( [ B ] ω [ I ] ) { φ } = 0 E181

It is seen that the general friction coefficient matrix can also be spectrally decomposed in standard spaces, so we have

[ B ] = [ Φ ] [ Ω ] [ Φ ] T E182

Comparing Eq. (B4) with Eq. (74), it is known that we can also spectrally decompose the lossy matrix in standard spaces

[ σ ] = [ Φ ] T [ Γ ] [ Φ ] E183

where,

Γ = Ω 1 = d i a g [ η 1 , η 2 , η 3 ] E184
.

References

  1. 1. Agyei A. Birman J. 1990 On the linear magnetoelectric effect. Journal of Physics: Condensed Matter, 2 13 3007 3020 , 0953-8984
  2. 2. Astrov D. 1960 The magnetoelectric effect in antiferromagnetics. Soviet Physics, Journal of Experimental and Theoretical Physics, 38 11 984 985 , 0038-5646
  3. 3. Athanasiadis C. Giotopoulos S. 2003 The Atkinson-Wilcox expansion theorem for electromagnetic chiral waves. Applied Mathematics Letters, 16 5 675 681 , 0893-9659
  4. 4. Athanasiadis C. Roach G. Stratis I. 2003 A time domain analysis of wave motions in chiral materials. Mathematische Nachrichten, 250 1 3 16 , 1522-2616
  5. 5. Benveniste Y. 1995 Magnetoelectric effect in fibrous composites with piezoelectric and piezomagnetic phases. Physical Review B 51 9 424 427 , 1098-0121
  6. 6. Bracke L. Van Vliet M. 1981 A broadband magneto-electric transducer using a composite material. International Journal of Electronics, 51 2 255 262 , 0020-7217
  7. 7. Bui M. Stuchly S. Costache G. 1991 Propagation of transients in dispersive dielectric media. IEEE Transactions on Microwave Theory and Techniques, 39 7 1165 1171 , 0018-9480
  8. 8. Cohen G. 2002 Higher-order numerical methods for transient wave equations, Springer, 354041598 Berlin
  9. 9. Dzyaloshinskii I. 1959 On the magneto-electrical effect in antiferromagnets. Soviet Physics, Journal of Experimental and Theoretical Physics, 10 3 628 669 , 0038-5646
  10. 10. Dvorak S. Dudley D. 1995 Propagation of ultra-wide-band electromagnetic pulses through dispersive media. IEEE Transactions on Electromagnetic Compatibility, 37 1 192 200 , 0018-9375
  11. 11. Gandhi O. Furse C. 1997 Currents induced in the human body for exposure to ultrawideband electromagnetic pulses. IEEE Transactions on Electromagnetic Compatibility, 39 2 174 180 , 0018-9375
  12. 12. Guo S. 2009 An eigen theory of electromagnetic waves based on the standard spaces. International Journal of Engineering Science, 47 3 405 412 , 0020-7225
  13. 13. Guo S. 2009 An Eigen theory of Static Electromagnetic Field for Anisotropic Media. Applied Mathematics and Mechanics, 30 5 643 652 , 0253-4827
  14. 14. Guo S. 2010 A second-order eigen theory for static electromagnetic fields. Journalof Modern Physics, 1 2 100 107 , 2153-1196
  15. 15. Guo S. 2010 An eigen theory of electro-magnetic acoustic waves in magnetoelectroelastic media. Acta Mechanica, 211 2 173 180 , 0001-5970
  16. 16. Guo S. 2010 An fully dynamical theory of piezoelectromagnetic waves. Acta Mech., 215 3 335 344 , 0001-5970
  17. 17. Haba Z. 2004 Green functions and propagation of waves in strongly inhomogeneous media. Journal of Physis A: Mathematical and General, 37 9 9295 9302 , 0305-4470
  18. 18. Harshe G. Dougherty J. Newnham R. 1993 Theoretical modelling of multilayer magnetoelectric composites. International Journal of Applied Electromagnetic Matter, 4 2 145 159 , 1383-5416
  19. 19. Jaggard D. Mickelson A. Papas C. 1979 On electromagnetic waves in chiral media. Applied Physics A: Materials Science and Processing, 18 2 211 216 , 0947-8396
  20. 20. Jin J. Zunoubi M. Donepudi K. Chew W. 1999 Frequency-domain and time-domain finite element solution of Maxwell’s equations using spectral Lanczos decomposition method. Computer Methods in Applied Mechanics and Engineering, 169 2 279 296 , 0045-7825
  21. 21. Kong J. 1986 Electromagnetic wave theory, John Wiley & Sons, 978-0-96681-439-2 New York
  22. 22. Lakhtakia A. 1994 Beltrami fields in chiral media. World Scientific, 9-81021-403-0
  23. 23. Lakhtakia A. Varadan V. K. Varadan V. V. 1989 Time-harmonic electromagnetic fields in chiral media. Springer-Verlag, 3-54051-317-5
  24. 24. Landau L. Lifshitz I. 1960 Electrodynamics of Continuous Media. Pergamon Press, 0-44486-165-3
  25. 25. Lee R. Madsen N. 1990 A mixed finite element formulation for Maxwell’s equations in the time domain. Journal Computational Physics, 88 2 284 304 , 0021-9991
  26. 26. Lindell I. Sihvola A. Tretyakov S. Viitanen A. 1994 Electromagnetic Waves in Chiral and Bi-isotropic Media. Artech House, ISBN, Boston
  27. 27. Lindell I. Sihvola A. Tretyakov S. Viitanen A. 1994 Electromagnetic Waves in Chiral and Bi-isotropic Media., Artech House, 1-59829-320-6
  28. 28. Li J. Chen Y. 2006 Analysis of a time-domain finite element method for 3D Maxwell’s equations in dispersive media. Computer Methods in Applied Mechanics and Engineering, 195 11 4220 4229 , 0045-7825
  29. 29. Liu Q. Fan G. 1999 Simulation of GPR in dispersive media using a frequency-dependent PSTD algorithm. IEEE Transactions on Geoscience and Remote Sensing, 37 5 2317 2324 , 0196-2892
  30. 30. Lu T. Zhang P. Cai W. 2004 Discontinuous Galerkin methods for dispersive and lossy Maxwell’s equations and PML boundaryconditions. Journal Computational Physics, 200 4 549 580 , 0021-9991
  31. 31. Monk P. 1992 A comparison of three mixed methods for time-dependent Maxwell’s equations. SIAM Journal of Science Statistics Computer, 13 10 1097 1122 , 0960-3174
  32. 32. Monk P. 2003 Finite Element methods for Maxwell’s equations, Clarendon Press, 0-19850-888-3
  33. 33. Pradhan T. 1993 Magnetic charge induced by electric charge in a magneto-electric medium. Europhysics Letters, 21 7 971 973 , 0295-5075
  34. 34. Van Run A. Terrell D. Scholing J. 1974 An in situ grown eutectic magnetoelectric composite material. Journal of Materials Science, 9 11 1710 1714 , 0022-2461
  35. 35. Yakhno V. Yakhno T. Kasap M. 2006 A novel approach for modeling and simulation of electromagnetic waves in anisotropic dielectrics. Internatioal Journal of Solids and Structures, 43 12 6261 6276 , 0020-7683

Written By

Shaohua Guo

Submitted: 08 October 2010 Published: 05 July 2011