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

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.


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.

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 Figure 1.
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, having the same periodicity as the 1D PC.Here, () and () are the usual conductivity and electric susceptibility for isotropic materials in the unit cell.Note that, for metals, σ = , for dielectrics, σ = −, 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: ∇ × e (r) = b (r) , Furthermore, current density j and magnetic induction b are related to magnetic field h and electric field e by materials equations: Given the 1D PC periodicity, the generalized conductivity and magnetic permeability can be expanded in Fourier series: being the 1D lattice reciprocal vectors.

Average Fields
The e(r) and b(r) "microscopic" fields within the periodic structure obey the Bloch theorem.Therefore, fields are obtained in Bloch wave form: where ∑  indicates that   = 0 term is excluded from the sum.
For macroscopic electromagnetic fields, the corresponding microscopic fields (8) are averaged over a length  much greater than the lattice constant  but much smaller than the wavelength 2/k ( ≪  ≪ 2/k).This averaging procedure smooths the rapid oscillations of the microscopic fields within the interval of length .Thus, in the case of wavelengths much larger than the lattice constant , Advances in Materials Science and Engineering 3 and macroscopic fields are described as plane waves: Another equivalent way to define the macroscopic fields is to eliminate the   ̸ = 0 terms in (8) expansions.This definition is correct when inequality ( 9) is accomplished and the   = 0 component in such expansions proves to be a smooth function in variations over distances much larger than the lattice constant.

Permittivity and Permeability Effective
In a homogeneous magnetic medium, the constitutive equations that describe the relationship of the displacement D and the magnetic induction B vectors to the electric field E and magnetic field H can be written as where ← →  and ← →  are the permittivity and permeability tensors.Maxwell's equations in this case are written in the form Removing the magnetic field H in (12) and (14) and using (11), the wave equation is obtained in electric field E terms: Equivalently, the electric field E can be eliminated from ( 12) and ( 14), and write the wave equation for magnetic field H: To obtain photonic crystal effective tensors, ( 16) and ( 17) must be compared with macroscopic fields exact equations that satisfy E 0 = e(0) and H 0 = h(0) amplitudes.Equations ( 6)-( 8) are substituted in ( 2) and ( 4) and expansions are used: Removing the magnetic field components, an electric field components equations system is obtained: where Here ← → I is the unitary dyadic.From equations with   ̸ = 0 (19), e(  ) components can be expressed in e(0) component terms.Substituting in the   = 0 equation, equations system (19) takes the form Here,

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: Advances in Materials Science and Engineering In the long wavelengths limit case ( = 0), the effective tensors are given explicitly by

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 of  are calculated.The set of parameters for the calculation are:   = 12.25 and   = 1 for dielectric matrix (host) and   = 13 and   = 8 for isotropic ferrite layer (inclusion) [5].
In Figures 2 and 3, effective permeability and permittivity tensors perpendicular components  eff =   =   and  eff =   =   behavior is showed versus magnetic inclusion filling fraction ( =   /, where   is the magnetic layer thickness).
Local approach was used in calculations; that is, k = 0 (see (26) formulas). eff permittivity and  eff permeability calculated are filling fraction  functions and are in agreement with classical quasistatic results: 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  ≪ 1 (9) gives nonlocal effective permittivity and permeability corrections.

Conclusions
The analytical relationships for 1D periodical photonic crystal effective permittivity and permeability tensors with periodicity  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][5][6][7], which uses four effective tensors (permittivity, permeability, and two magnetoelectric tensors).Obviously, permittivity and permeability definitions used here and the bianisotropic approach [4][5][6][7] are different for finite wave vectors.However, in the local case (k = 0) the definitions are in agreement.Unlike the previous works [5][6][7][8][9][10][11], in this limit, our formulas provide the dependence of the effective parameters in photonic crystals upon the filling fraction of the inclusions.