Table of Contents Author Guidelines Submit a Manuscript
Advances in Materials Science and Engineering
Volume 2016, Article ID 7451941, 5 pages
Research Article

Magnetodielectric and Metalomagnetic 1D Photonic Crystals Homogenization: Local Behavior

1Benemérita Universidad Autónoma de Puebla, Ciudad Universitaria, Boulevard Valsequillo Esquina, Avenida San Claudio s/n, Colonia San Manuel, Apartado Postal J-39, 72570 Puebla, PUE, Mexico
2División de Estudios de Posgrado e Investigación, Maestría en Ciencias en Ingeniería Mecánica, Instituto Tecnológico de Puebla, Avenida Tecnológico 420, Maravillas, 72220 Puebla, PUE, Mexico

Received 23 February 2016; Accepted 12 April 2016

Academic Editor: Francesco dell’Isola

Copyright © 2016 J. I. Rodríguez Mora et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


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.

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 Figure 1.

Figure 1: Superlattice or 1D binary photonic crystal. Periodicity in 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,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, [4].

According to these assumptions, you can study electromagnetic fields behavior in PC at structural level using “microscopic” Maxwell’s equations:Furthermore, current density and magnetic induction are related to magnetic field and electric field 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.

3. Average Fields

The and “microscopic” fields within the periodic structure obey the Bloch theorem. Therefore, fields are obtained in Bloch wave form:where indicates that 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 (). 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 ,and macroscopic fields are described as plane waves:Another equivalent way to define the macroscopic fields is to eliminate the terms in (8) expansions. This definition is correct when inequality (9) is accomplished and the 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 displacement and the magnetic induction vectors to the electric field and magnetic field can be written aswhere and are the permittivity and permeability tensors. Maxwell’s equations in this case are written in the formRemoving the magnetic field in (12) and (14) and using (11), the wave equation is obtained in electric field terms:Equivalently, the electric field can be eliminated from (12) and (14), and write the wave equation for magnetic field :To obtain photonic crystal effective tensors, (16) and (17) must be compared with macroscopic fields exact equations that satisfy and 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:whereHere is the unitary dyadic. From equations with (19), components can be expressed in component terms. Substituting in the equation, equations system (19) takes the formHere, is a 3 × 3 block, obtained from the inverse of the matrix (20). Similarly, the magnetic field components equations system is obtained:whereEliminating nonzero components in (22), an homogeneous algebraic equation is obtained:where is a 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:In the long wavelengths limit case (), the effective tensors are given explicitly by

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 of are calculated. The set of parameters for the calculation are: and for dielectric matrix (host) and and for isotropic ferrite layer (inclusion) [5].

In Figures 2 and 3, effective permeability and permittivity tensors perpendicular components and behavior is showed versus magnetic inclusion filling fraction (, where 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, (see (26) formulas). permittivity and 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 . Considering the condition (9) gives nonlocal effective permittivity and permeability corrections.

7. 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 [47], which uses four effective tensors (permittivity, permeability, and two magnetoelectric tensors). Obviously, permittivity and permeability definitions used here and the bianisotropic approach [47] are different for finite wave vectors. However, in the local case the definitions are in agreement. Unlike the previous works [511], in this limit, our formulas provide the dependence of the effective parameters in photonic crystals upon the filling fraction of the inclusions.

Competing Interests

The authors declare that they have no competing interests.


This work was supported by Faculty of Electronic Science (FCE), BUAP, Mexico.


  1. G. De Felice and N. Rizzi, “Homogenization for materials with microstructure,” Pressure Vessels and Piping Division of the American Society of Mechanical Engineers, vol. 369, pp. 33–38, 1997. View at Google Scholar
  2. A. Cecchi and N. L. Rizzi, “Heterogeneous elastic solids: a mixed homogenization-rigidification technique,” International Journal of Solids and Structures, vol. 38, no. 1, pp. 29–36, 2000. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  3. J. D. Joannopoulos, S. G. Johnson, J. N. Winn, and R. D. Meade, Photonic Crystals: Molding the Flow of Light, Princeton University Press, Princeton, NJ, USA, 2008.
  4. P. Halevi and F. Pérez-Rodríguez, “From photonic crystals (via homogenization) to metamaterials,” in Proceedings of the Complex Photonic Media, vol. 6320 of Proceedings of SPIE, San Diego, Calif, USA, August 2006. View at Publisher · View at Google Scholar
  5. E. Reyes-Ayona and P. Halevi, “Mean field theory of metallo-dielectric photonic crystals with magnetic components: the long-wavelength limit,” in Photonic Metamaterials, 66380G, vol. 6638 of Proceedings of SPIE, September 2007. View at Publisher · View at Google Scholar
  6. V. Cerdán-Ramírez, B. Zenteno-Mateo, M. P. Sampedro, M. A. Palomino-Ovando, B. Flores-Desirena, and F. Pérez-Rodríguez, “Anisotropy effects in homogenized magnetodielectric photonic crystals,” Journal of Applied Physics, vol. 106, Article ID 103520, 2009. View at Publisher · View at Google Scholar · View at Scopus
  7. J. A. Reyes-Avendaño, U. Algredo-Badillo, P. Halevi, and F. Pérez-Rodríguez, “From photonic crystals to metamaterials: the bianisotropic response,” New Journal of Physics, vol. 13, Article ID 073041, 2011. View at Publisher · View at Google Scholar · View at Scopus
  8. A. A. Krokhin and E. Reyes, “Homogenization of magnetodielectric photonic crystals,” Physical Review Letters, vol. 93, no. 2, Article ID 023904, 2004. View at Publisher · View at Google Scholar · View at Scopus
  9. E. Reyes, A. A. Krokhin, and J. Roberts, “Effective dielectric constants of photonic crystal of aligned anisotropic cylinders and the optical response of a periodic array of carbon nanotubes,” Physical Review B, vol. 72, no. 15, Article ID 155118, 2005. View at Publisher · View at Google Scholar
  10. A. A. Krokhin, E. Reyes, and L. Gumen, “Low-frequency index of refraction for a two-dimensional metallodielectric photonic crystal,” Physical Review B, vol. 75, no. 4, Article ID 045131, 2007. View at Publisher · View at Google Scholar
  11. G. P. Ortiz, B. E. Martínez-Zérega, B. S. Mendoza, and W. L. Mochán, “Effective optical response of metamaterials,” Physical Review B, vol. 79, no. 24, Article ID 245132, 9 pages, 2009. View at Publisher · View at Google Scholar · View at Scopus