## Abstract

The Cell Method, a new efficient numerical method suitable for working with periodic structures having anisotropic inhomogeneous media with curved shapes, is proposed in order to calculate the band gap of 2D photonic crystals for in-plane propagation of TM and TE waves. Moreover some numerical comparisons with other numerical methods will be provided.

## 1. Introduction

In the study of the properties of photonic band gap materials it is important to find if some bands of frequencies exist (photonic band gap) where the propagation of the electromagnetic waves is prohibited in all the directions. In order to find these frequencies a lot of methods have been proposed up to now. The Plane Wave Expansion Method (PWE)  works with inhomogeneous media having curved shapes and deals with large not sparse hermitian matrices. The Finite Elements Method (FEM) works with inhomogeneous media having curved shapes as well and a generalized eigenvalues problem with large sparse (quasi band diagonal) matrices has to be solved. The Finite Difference Method (FD) works with inhomogeneous media without curved shapes and deals with large sparse (quasi band diagonal) hermitian matrices  in non generalized eigenvalues problems. In this paper we will use a new method, The Cell Method (CM) , that works with inhomogeneous media having curved shapes, like FEM, and deals with large sparse (quasi band diagonal) hermitian matrices in non generalized eigenvalues problems, like FD.

## 2. The Cell Method

The Cell Method is based on the followings tools:

• Two kinds of oriented space elements (space cells): the space elements endowed with an inner orientation that we will call primal space cells (primal points P, lines L, surfaces S and volumes V) and the space elements endowed with an outer orientation that we will call dual space cells (dual points , lines , surfaces and volumes ) (Fig.1a).
• Two oriented cell complexes and a relationship of duality between them. A space complex, synonymous of three dimensional grid, is a structured collection of points, lines, surfaces and volumes (cells).The primal complex is made by inner oriented space cells, the dual complex is made by outer oriented space cells and the relationship of duality means that we need a biunivocal matching respectively between the primal points P and the dual volumes , the primal lines L and the dual surfaces , the primal surfaces S and the dual lines and between the primal volumes V and the dual points . Fig. 1. (a)Inner and outer orientations of the space cells in the 3D space (b)Primal and dual complexes and global variables associated with their primal and dual space cells respectively.

• The use of the global (integral) variables to represent physical variables and their physically coherent association with the oriented space cells (Fig.1b):

Electric voltage V[L](= Ld l), Magnetic flux Φ[S](= S B·d s), Magnetic voltage F[](= H·d l), Electric flux ψ[](= D·d s), Electric charge content Qc[](= ρdv), Electric current I[](= J·ds).

The electromagnetic laws, that are links between global variables, can be divided into two categories, according to the type of variables that are linked:

• Topological equations (Field equations).These equations link global variables associated with space cells belonging to the same kind of complex (either the primal or the dual one).Since they do not concern metric concepts, they can be enforced on the cell complexes in a discrete not approximated form by using appropriate incidence matrices. If we denote with G, C, D and = DT, = CT, = -GT the incidence matrices  related to the primal and dual cell complexes respectively, which are the discrete counterparts of the differential operators gradient, curl and divergence, the topological equations in the frequency domain can be expressed as follows:
• Faraday-Neumann law Ca = -jωΦ
• Magnetic Gauss law = 0
• Maxwell-Ampére law C̃F = ψ + I
• Electric Gauss law D̃ψ = Q c

where V, Φ, F, ψ, I, Q c are scalar arrays.

• Constitutive relations.These equations link global variables associated with space cells belonging to different kind of cell complexes (i.e. the electric constitutive relation links the electric voltage V, associated with the primal lines L, and the electric flux ψ, associated with the dual surfaces ).They concern metric concepts and material parameters and they can be enforced on the cell complexes only approximatively in a discrete form by using suitable constitutive matrices:
• Electric constitutive relation ψ = M ε V
• Magnetic constitutive relation Φ = M μ F

If an orthogonality occurs (L and S, as in the Fig.1b), then the matrices M ε,Mμ can be built diagonal  and assure, at the same time, the convergence of the method. In general, both for orthogonal or non-orthogonal grids, it is possible to build the matrices M ε,M μ by the Microcell Interpolation Scheme (MIS) . Fig. 2. Examples of 2D photonic crystals and periodicity of the unit cells in the xy-plane (a) for rectangular unit cells (b) for quadrangular unit cells. (c) Dependence of some voltages on the boundary of a rectangular unit cell by the periodic boundary condition.

## 3. Analysis of 2d periodic structures with the Cell Method

The 2D photonic crystals are structures which have a transverse periodicity (xy-plane) and a longitudinal uniformity (z-axis).In order to study their properties it is enough to deal with a unit cell, the atomic part of the periodic structure (Fig.2), with appropriate periodic boundary conditions .In case of in-plane propagation, which is the propagation only in the transverse plane, it is possible to analyze separately the TMz modes from the TEz modes. Fig. 3. TMz case Fig. 4. Examples of 2Dprimal and dual grids (a) Triangular primal cell complex, barycentric dual cell complex for the TMz cases (b) Barycentric primal cell complex, triangular dual cell complex for the TEz cases

#### 3.1 Two dimensional TMz case

In CM for the TMz case we have to work on very simple 3D cell complexes that are, in particular, 2D extrude cell complexes along the z-axes (Fig.3). This fact allows us to work like in a pure 2D case (Fig.4a) by projecting the 3D cell complexes on the xy-plane, where:

• The 2D unit cell is discretized by a triangular primal cell complex and a barycentric dual cell complex (Fig.4a).
• The association global variables → space cells is:

Electric voltage V → Primal point P, Magnetic flux Φ → Primal line L Magnetic voltage F → Dual line , Electric flux ψ → Dual surface

• The topological equations are:

Faraday-Neumann law: GV = -Φ, Maxwell-Ampére law: C̃F = ψ.

• The constitutive relations are:

Electric: ψ = M ε V, Magnetic: F = M ν Φ.

From the previous equations and since = -G T in the 2D cases, we obtain the matricial equation:

$GTMνGV=ω2MεV$

that is the Helmholtz equation in CM for the TMz case, where V is the array of the N(=number of nodes=number of primal points P) electric voltages. If we build the matrices M ν and M ε with the MIS, for triangular grids the G T M ν G is a symmetric positive semidefinite band diagonal matrix and the M ε is a diagonal positive definite matrix. For a rectangular unit cell of sizes [a, b] (Fig.2a) the appropriate boundary condition is :

$V(x+a,y+b)=e−jβxa−jβybVxy$

where βx and βy are the phase constants in the x and y directions. Since some voltages on the boundary are dependent of other voltages, because the periodic conditions (Fig.2c), it is possible to reduce the dimension of the problem by means of the following relation:

$V=T(βx,βy)V′$

where V′ is a vector of the independent voltages, T(βx,βy) is a sparse matrix whose entries T ij are 1 if V i = Vj, either e -xa or e -yb or e -jβxa-jβyb according to the boundary condition between V i and Vj, 0 otherwise. Thus the Helmholtz equation with the periodic boundary conditions reduces to:

$AV′=ω2BV′$

where A = T H G T M ν GT is a hermitian semidefinite positive matrix and B = T H M ε T is a diagonal real positive definite matrix. This generalized eigenvalues problem can be transformed into a simple eigenvalues problem:

$MTMV”=ω2V”$

where D = √B, M TM = D -1 AD -1 and V’’ = DV′.The matrix M TM is still sparse hermitian positive semidefinite and, by an appropriate renumbering of the nodes, the structure of this matrix becomes quasi band diagonal. Thus, in order to find its eigenvalues, can be used the same efficient method developed for the FD method in .

#### 3.2 Two dimensional TEz case

In CM for the TEz case it is possible to do the same considerations as the TMz cases by swaping the role of the primal and dual cell complexes (Fig.4b):

• The 2D unit cell is discretized by a barycentric primal cell complex and a triangular dual cell complex.
• The association global variables → space cells is:

Electric voltage V → Primal line L, Magnetic flux Φ → Primal surface S, Magnetic voltage F → Dual point , Electric flux ψ → Dual line .

• The topological equations are:

Maxwell-Ampére law: G̃F = ψ, Faraday-Neumann law: Ca = - Φ.

• The constitutive relations are:

Electric: V = M ε-1 ψ, Magnetic: Φ = M μ F.

From the previous equations and since T = C in the 2D cases, we obtain the matricial equation:

$G˜TMε−1G˜V=ω2MμF$

that is the Helmholtz equation in CM for the TEz modes, where F is the array of the N(=number of nodes=number of dual points ) magnetic voltages. For a rectangular unit cell of sizes [a, b], by the boundary condition:

$F(x+a,y+b)=e−jβxa−jβybFxy$

a similar eigenvalue problem like Eq.(5) has to be solved:

$MTMF”=ω2F”$

## 4. Numerical results and discussion

In order to evaluate the reliability of CM, some tests have been performed. The first test was the calculus of the photonic band structure for TM and TE waves in a square lattice (lattice constant a) of dielectric rods (r = 12, radius=0.2a). The results, compared with PWE , show the existence of a TM photonic band gap (Fig.5a). In (Fig.5b) the normalized frequencies of the first TM mode in M and the second TM mode in X, for the first test, calculated by CM and FEM, are compared for different number of grid points. The second test was the calculus of the photonic band structure for TM and TE waves in a triangular lattice (lattice constant a) of air columns (radius=0.48a) in a dielectric substrate (r = 13).The results, in excellent agreement with  where a PWE has been used, exhibit a complete band gap (Fig.6a). The third test was the calculus of the photonic band structure for TM and TE waves in a square lattice of anisotropic dielectric rods (xx = yy = 23.04, zz = 38.44, Filling Fraction=0.4) (Fig.6b). The results, in good agreement with  where a PWE has been used, exhibit a complete band gap between (0.219–0.254) 2πc/a.In conclusion, the good agreement of the CM results with the results of other methods shows the effectiveness of the present approach. Moreover, in the test 1, CM exhibits more accuracy than FEM for the same number of grid points.

## 5. Conclusion

In this paper an efficient and accurate method for the analysys of the band gap of 2D photonic crystals for in-plane propagation of TM and TE modes with anisotropic media has been proposed and some comparisons with the results of other methodologies have been provided. Further work is however necessary to extend the method to the out of plane propagation cases. Fig. 5. (a) Results of the test 1 (b) A comparison of some CM and FEM results in the test 1. Fig. 6. (a) Results of the test 2 (b) Results of the test 3.

1. J.D. Joannopulos, R.D. Meade, and J.N. Winn, Photonic Crystals (Princeton University Press, 1995).

2. H.Y.D. Yang,”Finite Difference Analysis of 2-D Photonic Crystals,” IEEE Trans. Microwave Theory Tech. 442688–2695 (1996). [CrossRef]

3. E. Tonti, ”Finite Formulation of the Electromagnetic Field,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001) 1–44.

4. M. Marrone, ”Computational Aspects of Cell Method in Electrodynamics,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001), 317–356.

5. S.G. Johnson, S. Fan, P.R. Villeneuve, J.D. Joannopulos, and L.A. Kolodziejski, ”Guided modes in photonic crystal slabs,” Phys. Rev. B 3, 5751–5758 (1999) [CrossRef]

6. Z.Y. Li, B.Y. Gu, and G.Z. Yang, ”Improvement of absolute band gaps in 2Dphotonic crystals by anisotropy in dielectricity,” Eur.Phys. J. B 11, 65–73 (1999)

7. M. Clemens and T. Weiland, ”Discrete Electromagnetism with the Finite Integration Technique,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001), 65–87.

### References

• View by:
• |
• |
• |

1. J.D. Joannopulos, R.D. Meade, and J.N. Winn, Photonic Crystals (Princeton University Press, 1995).
2. H.Y.D. Yang,”Finite Difference Analysis of 2-D Photonic Crystals,” IEEE Trans. Microwave Theory Tech. 442688–2695 (1996).
[Crossref]
3. E. Tonti, ”Finite Formulation of the Electromagnetic Field,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001) 1–44.
4. M. Marrone, ”Computational Aspects of Cell Method in Electrodynamics,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001), 317–356.
5. S.G. Johnson, S. Fan, P.R. Villeneuve, J.D. Joannopulos, and L.A. Kolodziejski, ”Guided modes in photonic crystal slabs,” Phys. Rev. B 3, 5751–5758 (1999)
[Crossref]
6. Z.Y. Li, B.Y. Gu, and G.Z. Yang, ”Improvement of absolute band gaps in 2Dphotonic crystals by anisotropy in dielectricity,” Eur.Phys. J. B 11, 65–73 (1999)
7. M. Clemens and T. Weiland, ”Discrete Electromagnetism with the Finite Integration Technique,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001), 65–87.

#### 1999 (2)

S.G. Johnson, S. Fan, P.R. Villeneuve, J.D. Joannopulos, and L.A. Kolodziejski, ”Guided modes in photonic crystal slabs,” Phys. Rev. B 3, 5751–5758 (1999)
[Crossref]

Z.Y. Li, B.Y. Gu, and G.Z. Yang, ”Improvement of absolute band gaps in 2Dphotonic crystals by anisotropy in dielectricity,” Eur.Phys. J. B 11, 65–73 (1999)

#### 1996 (1)

H.Y.D. Yang,”Finite Difference Analysis of 2-D Photonic Crystals,” IEEE Trans. Microwave Theory Tech. 442688–2695 (1996).
[Crossref]

#### Clemens, M.

M. Clemens and T. Weiland, ”Discrete Electromagnetism with the Finite Integration Technique,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001), 65–87.

#### Fan, S.

S.G. Johnson, S. Fan, P.R. Villeneuve, J.D. Joannopulos, and L.A. Kolodziejski, ”Guided modes in photonic crystal slabs,” Phys. Rev. B 3, 5751–5758 (1999)
[Crossref]

#### Gu, B.Y.

Z.Y. Li, B.Y. Gu, and G.Z. Yang, ”Improvement of absolute band gaps in 2Dphotonic crystals by anisotropy in dielectricity,” Eur.Phys. J. B 11, 65–73 (1999)

#### Joannopulos, J.D.

S.G. Johnson, S. Fan, P.R. Villeneuve, J.D. Joannopulos, and L.A. Kolodziejski, ”Guided modes in photonic crystal slabs,” Phys. Rev. B 3, 5751–5758 (1999)
[Crossref]

J.D. Joannopulos, R.D. Meade, and J.N. Winn, Photonic Crystals (Princeton University Press, 1995).

#### Johnson, S.G.

S.G. Johnson, S. Fan, P.R. Villeneuve, J.D. Joannopulos, and L.A. Kolodziejski, ”Guided modes in photonic crystal slabs,” Phys. Rev. B 3, 5751–5758 (1999)
[Crossref]

#### Kolodziejski, L.A.

S.G. Johnson, S. Fan, P.R. Villeneuve, J.D. Joannopulos, and L.A. Kolodziejski, ”Guided modes in photonic crystal slabs,” Phys. Rev. B 3, 5751–5758 (1999)
[Crossref]

#### Li, Z.Y.

Z.Y. Li, B.Y. Gu, and G.Z. Yang, ”Improvement of absolute band gaps in 2Dphotonic crystals by anisotropy in dielectricity,” Eur.Phys. J. B 11, 65–73 (1999)

#### Marrone, M.

M. Marrone, ”Computational Aspects of Cell Method in Electrodynamics,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001), 317–356.

J.D. Joannopulos, R.D. Meade, and J.N. Winn, Photonic Crystals (Princeton University Press, 1995).

#### Tonti, E.

E. Tonti, ”Finite Formulation of the Electromagnetic Field,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001) 1–44.

#### Villeneuve, P.R.

S.G. Johnson, S. Fan, P.R. Villeneuve, J.D. Joannopulos, and L.A. Kolodziejski, ”Guided modes in photonic crystal slabs,” Phys. Rev. B 3, 5751–5758 (1999)
[Crossref]

#### Weiland, T.

M. Clemens and T. Weiland, ”Discrete Electromagnetism with the Finite Integration Technique,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001), 65–87.

#### Winn, J.N.

J.D. Joannopulos, R.D. Meade, and J.N. Winn, Photonic Crystals (Princeton University Press, 1995).

#### Yang, G.Z.

Z.Y. Li, B.Y. Gu, and G.Z. Yang, ”Improvement of absolute band gaps in 2Dphotonic crystals by anisotropy in dielectricity,” Eur.Phys. J. B 11, 65–73 (1999)

#### Yang, H.Y.D.

H.Y.D. Yang,”Finite Difference Analysis of 2-D Photonic Crystals,” IEEE Trans. Microwave Theory Tech. 442688–2695 (1996).
[Crossref]

#### Eur.Phys. J. B (1)

Z.Y. Li, B.Y. Gu, and G.Z. Yang, ”Improvement of absolute band gaps in 2Dphotonic crystals by anisotropy in dielectricity,” Eur.Phys. J. B 11, 65–73 (1999)

#### IEEE Trans. Microwave Theory Tech. (1)

H.Y.D. Yang,”Finite Difference Analysis of 2-D Photonic Crystals,” IEEE Trans. Microwave Theory Tech. 442688–2695 (1996).
[Crossref]

#### Phys. Rev. B (1)

S.G. Johnson, S. Fan, P.R. Villeneuve, J.D. Joannopulos, and L.A. Kolodziejski, ”Guided modes in photonic crystal slabs,” Phys. Rev. B 3, 5751–5758 (1999)
[Crossref]

#### Other (4)

J.D. Joannopulos, R.D. Meade, and J.N. Winn, Photonic Crystals (Princeton University Press, 1995).

E. Tonti, ”Finite Formulation of the Electromagnetic Field,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001) 1–44.

M. Marrone, ”Computational Aspects of Cell Method in Electrodynamics,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001), 317–356.

M. Clemens and T. Weiland, ”Discrete Electromagnetism with the Finite Integration Technique,” in Geometric Methods for Computational Electromagnetics, PIER32, F.L. Teixeira and J.A. Kong, ed.(EMW Publishing2001), 65–87.

### Cited By

OSA participates in Crossref's Cited-By Linking service. Citing articles from OSA journals and other participating publishers are listed here.

### Figures (6)

Fig. 1. (a)Inner and outer orientations of the space cells in the 3D space (b)Primal and dual complexes and global variables associated with their primal and dual space cells respectively.
Fig. 2. Examples of 2D photonic crystals and periodicity of the unit cells in the xy-plane (a) for rectangular unit cells (b) for quadrangular unit cells. (c) Dependence of some voltages on the boundary of a rectangular unit cell by the periodic boundary condition.
Fig. 4. Examples of 2Dprimal and dual grids (a) Triangular primal cell complex, barycentric dual cell complex for the TMz cases (b) Barycentric primal cell complex, triangular dual cell complex for the TEz cases
Fig. 5. (a) Results of the test 1 (b) A comparison of some CM and FEM results in the test 1.
Fig. 6. (a) Results of the test 2 (b) Results of the test 3.

### Equations (8)

$G T M ν GV = ω 2 M ε V$
$V ( x + a , y + b ) = e − j β x a − j β y b V x y$
$V = T ( β x , β y ) V ′$
$AV ′ = ω 2 BV ′$
$M TM V ” = ω 2 V ”$
$G ˜ T M ε − 1 G ˜ V = ω 2 M μ F$
$F ( x + a , y + b ) = e − j β x a − j β y b F x y$
$M TM F ” = ω 2 F ”$