# Magnetodielectric and Metalomagnetic 1D Photonic Crystals Homogenization:ε-μ Local Behavior

Authors: J. I. Rodríguez Mora; J. Flores Méndez; B. Zenteno Mateo; F. Severiano Carrillo; R. C. Ambrosio Lázaro
Journal: Advances in Materials Science and Engineering
 (2016)

---

## Abstract

A theory for calculating the effective optic response of photonic crystals with metallic and magnetic inclusions is reported, for the case when the wavelength of the electromagnetic fields is much larger than the lattice constant. The theory is valid for any type of Bravais lattice and arbitrary form of inclusions in the unitary cell. An equations system is obtained for macroscopic magnetic field and magnetic induction components expanding microscopic electromagnetic fields in Bloch waves. Permittivity and permeability effective tensors are obtained comparing the equations system with an anisotropic nonlocal homogenous medium. In comparison with other homogenization theories, this work uses only two tensors: nonlocal permeability and permittivity. The proposal showed here is based on the use of permeability equations, which are exact and very simple. We present the explicit form of these tensors in the case of binary 1D photonic crystals.


---

## Body

## 1. Introduction

The material properties and their potential applications research are very important in the industry and scientific areas. However, all these impacts could not be possible without the great advances developed in the last 100 years in materials science. Thanks to metallurgy, ceramics, and plastics advances, engineers have been able to modify a wide range of mechanical and electrical properties in materials. Semiconductor physics advances have allowed adjusting the conductive properties of certain materials, initiating the electronics revolution. Scientists have discovered high-temperature superconductors with new alloys and ceramics. It is impossible to cite all applications derived from materials science. In recent years, scientists have a new challenge, to study artificial materials and their optical properties. A similar challenge is to study the mechanical properties of artificial materials [1, 2]. A great advance in this area is known as nanophotonics, which tries to solve the information transmission lossless problem through light. Continuous search for new technologies has brought surprising properties to a new materials class known as photonic crystals (PF) [3]. These artificial materials can serve as the new raw material for photonic circuits age. These applications include both submicron sized laser sources as conducting channels (optical waveguides) and logical components (amplifiers, transistors photonic, etc.). Essentially, a photonic crystal is a material with a periodic arrangement composed of multiple elements distributed which scatters light in a coherent and cooperative way. This produces, a similar mechanism that have the electrons in semiconductors, a forbidden energy range for photons propagation, in this case a photonic forbidden gap (gap), which allows control light propagation. By combination and structuration of these elements, photonic circuits can be designed similar to electronics with cables, switches, splitters, and so forth.If the wavelength of the light is greater than the photonic crystal lattice constant, the properties can be described by means of homogeneous medium effective parameters. These parameters can explain the optical properties of photonic crystals in a simple and clear manner. The calculation of these parameters is one of the great interesting theoretical problems.Homogenization theory of photonic crystals with metallic and magnetic inclusions is presented. Explicit formulas for permittivity and not local permeability tensors calculation are obtained and the application of these relations to calculate the effective parameters of a dielectric-magnetic 1D periodicity photonic crystal is illustrated. Finally, the theoretical results are compared with respect to theories based on bianisotropic response.


## 2. Material Characterization

### 2.1. Lattice Homogenization with Metalomagnetic Inclusions

Let us consider a 1D PC composed of metallic and magnetic inclusions. A 1D binary crystal (composed of two materials) representation is shown in Figure1.Figure 1
Superlattice or 1D binary photonic crystal. Periodicity inZ direction.The metalomagnetic isotropic inclusions properties of the photonic crystal will be fully described by two material parameters: the permittivity and permeability, which are position functions; therefore,(1)εz=ε0+iωσ^z,μzσ^z=σz-iωε0χz,having the same periodicity as the 1D PC. Here, σ(z) and χ(z) are the usual conductivity and electric susceptibility for isotropic materials in the unit cell. Note that, for metals, σ^=σ, for dielectrics, σ^=-iωχ, and, for air, σ^=0 [4].According to these assumptions, you can study electromagnetic fields behavior in PC at structural level using “microscopic” Maxwell’s equations:(2)∇×hr=jr-iωε0er,(3)∇·er=-iε0ω∇·jr,(4)∇×er=iωbr,(5)∇·br=0.Furthermore, current density j and magnetic induction b are related to magnetic field h and electric field e by materials equations:(6)jr=σ^zer,br=μzhr.Given the 1D PC periodicity, the generalized conductivity and magnetic permeability can be expanded in Fourier series:(7)σ^z=∑Gzσ^GzeiGzz,μz=∑GzμGzeiGzz,Gz being the 1D lattice reciprocal vectors.


## 2.1. Lattice Homogenization with Metalomagnetic Inclusions

Let us consider a 1D PC composed of metallic and magnetic inclusions. A 1D binary crystal (composed of two materials) representation is shown in Figure1.Figure 1
Superlattice or 1D binary photonic crystal. Periodicity inZ direction.The metalomagnetic isotropic inclusions properties of the photonic crystal will be fully described by two material parameters: the permittivity and permeability, which are position functions; therefore,(1)εz=ε0+iωσ^z,μzσ^z=σz-iωε0χz,having the same periodicity as the 1D PC. Here, σ(z) and χ(z) are the usual conductivity and electric susceptibility for isotropic materials in the unit cell. Note that, for metals, σ^=σ, for dielectrics, σ^=-iωχ, and, for air, σ^=0 [4].According to these assumptions, you can study electromagnetic fields behavior in PC at structural level using “microscopic” Maxwell’s equations:(2)∇×hr=jr-iωε0er,(3)∇·er=-iε0ω∇·jr,(4)∇×er=iωbr,(5)∇·br=0.Furthermore, current density j and magnetic induction b are related to magnetic field h and electric field e by materials equations:(6)jr=σ^zer,br=μzhr.Given the 1D PC periodicity, the generalized conductivity and magnetic permeability can be expanded in Fourier series:(7)σ^z=∑Gzσ^GzeiGzz,μz=∑GzμGzeiGzz,Gz being the 1D lattice reciprocal vectors.


## 3. Average Fields

Thee(r) and b(r) “microscopic” fields within the periodic structure obey the Bloch theorem. Therefore, fields are obtained in Bloch wave form:(8)er=eik·r∑GzeGzeiGzz=eik·reGz=0+eik·r∑Gz′⁡eGzeiGzz,br=eik·r∑GzbGzeiGzz=eik·rbGz=0+eik·r∑Gz′⁡bGzeiGzz,where ∑′ indicates that Gz=0 term is excluded from the sum.For macroscopic electromagnetic fields, the corresponding microscopic fields (8) are averaged over a length L much greater than the lattice constant a but much smaller than the wavelength 2π/k (a≪L≪2π/k). This averaging procedure smooths the rapid oscillations of the microscopic fields within the interval of length L. Thus, in the case of wavelengths much larger than the lattice constant a,(9)ka≪1,and macroscopic fields are described as plane waves:(10)Er≡er=eGz=0eik·r=E0eik·r,Br≡br=bGz=0eik·r=B0eik·r.Another equivalent way to define the macroscopic fields is to eliminate the Gz≠0 terms in (8) expansions. This definition is correct when inequality (9) is accomplished and the Gz=0 component in such expansions proves to be a smooth function in variations over distances much larger than the lattice constant.


## 4. Permittivity and Permeability Effective

In a homogeneous magnetic medium, the constitutive equations that describe the relationship of the displacementD and the magnetic induction B vectors to the electric field E and magnetic field H can be written as(11)D=ε⃡·E,B=μ⃡·H,where ε⃡ and μ⃡ are the permittivity and permeability tensors. Maxwell’s equations in this case are written in the form(12)∇×H=-iωε⃡·E,(13)∇·D=0,(14)∇×E=iωB,(15)∇·B=0.Removing the magnetic field H in (12) and (14) and using (11), the wave equation is obtained in electric field E terms:(16)k×μ⃡-1·k×E0+ω2ε⃡·E0=0.Equivalently, the electric field E can be eliminated from (12) and (14), and write the wave equation for magnetic field H:(17)k×ε⃡-1·k×H0+ω2μ⃡·H0=0.To obtain photonic crystal effective tensors, (16) and (17) must be compared with macroscopic fields exact equations that satisfy E0=e(0) and H0=h(0) amplitudes. Equations (6)–(8) are substituted in (2) and (4) and expansions are used:(18)1εz=∑GzηGze-iGzz,1μz=∑GzξGze-iGzz.Removing the magnetic field components, an electric field components equations system is obtained:(19)∑Gz′D⃡Ek;Gz,Gz′·eGz′=0,where(20)D⃡Ek;Gz,Gz′=k+Gzz^×ξGz-Gz′k+Gz′z^×I⃡+iωσ^Gz-Gz′+ε0ω2I⃡.Here I⃡ is the unitary dyadic. From equations with Gz≠0 (19), e(Gz) components can be expressed in e(0) component terms. Substituting in the Gz=0 equation, equations system (19) takes the form(21)D⃡E-1k;Gz=0,Gz′=0-1·eGz=0=0.Here, D⃡E-10,0 is a 3 × 3 block, obtained from the inverse of the matrix (20). Similarly, the magnetic field e(Gz) components equations system is obtained:(22)∑Gz′D⃡Hk;Gz,Gz′·hGz′=0,where(23)D⃡Hk;Gz,Gz′=k+Gzz^×ηGz-Gz′k+Gz′z^×I⃡+ω2μGz-Gz′I⃡.Eliminating Gz nonzero components in (22), an h(0) homogeneous algebraic equation is obtained:(24)D⃡H-1k;Gz=0,Gz′=0-1·hGz=0=0,where D⃡H-10,0 is a 3×3 block, from the inverse of the matrix (23).


## 5. Effective Parameters

Equating (16) and (17) with (21) and (24), respectively, a system of two algebraic equations that define the effective tensors ε⃡ and μ⃡ is obtained:(25)ε⃡effk=1ω2D⃡E-1k;0,0-1-k×μ⃡eff·k×I⃡,μ⃡effk=1ω2D⃡H-1k;0,0-1-k×ε⃡eff·k×I⃡.In the long wavelengths limit case (k=0), the effective tensors are given explicitly by(26)ε⃡eff0=1ω2D⃡E-10;0,0-1,μ⃡eff0=1ω2D⃡H-10;0,0-1.


## 6. Results

The results obtained with the theory developed in the previous section are described in this section. The effective permittivity and permeability for a 1D photonic crystal composed of alternating layers of dielectric and magnetic materials in the direction ofz are calculated. The set of parameters for the calculation are: εb=12.25 and μb=1 for dielectric matrix (host) and εa=13 and μa=8 for isotropic ferrite layer (inclusion) [5].In Figures2 and 3, effective permeability and permittivity tensors perpendicular components εeff=εxx=εyy and μeff=μxx=μyy behavior is showed versus magnetic inclusion filling fraction (f=da/a, where da is the magnetic layer thickness).Figure 2
Effective permittivity versus filling fraction.Figure 3
Effective permeability versus filling fraction graphic representation.Local approach was used in calculations; that is,k=0 (see (26) formulas). εeff permittivity and μeff permeability calculated are filling fraction f functions and are in agreement with classical quasistatic results: εeff=εaf+εb(1-f), μeff=μaf+μb(1-f).The formulas obtained in the previous section are validated with these results and open the possibility for future studies of electromagnetic waves propagation in 1D PC nonlocal effects. In such case, it will be necessary to include in the formulas obtained in this work (25) the dispersion relation k(ω). Considering the condition ka≪1 (9) gives nonlocal effective permittivity and permeability corrections.


## 7. Conclusions

The analytical relationships for 1D periodical photonic crystal effective permittivity and permeability tensors with periodicitya were obtained. Demonstrating that effective ε-μ material response of photonic crystal is not local (25), thus the photonic crystal homogenized effective parameters depend on wave vector. We point out that theoretical formalism developed in this work uses only two effective tensors: permittivity and permeability (25), unlike the bianisotropic response [4–7], which uses four effective tensors (permittivity, permeability, and two magnetoelectric tensors). Obviously, permittivity and permeability definitions used here and the bianisotropic approach [4–7] are different for finite wave vectors. However, in the local case (k=0) the definitions are in agreement. Unlike the previous works [5–11], in this limit, our formulas provide the dependence of the effective parameters in photonic crystals upon the filling fraction of the inclusions.



---

*Source: 7451941-2016-05-12.xml*