To begin, Eq. (1) is simplified using the curl of curl identity: Display Formula
(2)We assume that there are no free charges, which demands that Display FormulaEquation (2) now simplifies to Display Formula
(3)Expressing Eq. (3) in terms of second order central finite differences gives the following: Display Formula
(4)where and are the electric field vectors along and , and and are the physical distances between the mesh elements along the and axis, respectively. Equation (4) is fully decoupled; that is, changes in happen only along (-axis), and changes in happen only along (-axis). This allows it to be rewritten as Display Formula
(5)To solve this eigenvalue problem, a mesh is generated along a unit cell of the PC, with values corresponding to mesh elements. Figure 5 shows an example of the mesh in which the electric field at each position is labeled . In the finite differences equation, the field at a given position, , corresponds to one position, . The and terms correspond to adjacent mesh cells in the lattice. Elements within the unit cell are shown as bounded by solid lines and those within adjacent cells are represented by dashed lines. The 1 further discusses the mesh elements of the unit cell, as well as the required mesh size at which calculated bands converge to a common value.