The finite element method applied to the study of two-dimensional photonic crystals

Imanol Andonegui,Angel J. Garcia-Adeva

DOI: https://doi.org/10.1364/OE.21.004072

2012-10-29

Abstract:Calculations of the photonic band structure, transmission coefficients, and quality factors of various two-dimensional, periodic and aperiodic, dielectric photonic crystals by using the finite element method (FEM) are reported. The fundamental equations governing the propagation of electromagnetic waves in inhomogeneous media are revisited together with the boundary conditions required for each of the performed calculations. A detailed account of the eigenvalue and harmonic propagation analysis of the electromagnetic problem is reported for several periodic and finite-length structures. It is found that this method reproduces quite well previous results for these lattices obtained with the standard plane wave method with regards to the eigenvalue analysis (photonic band structure calculations). However, in contrast with frequency methods, the finite element method easily allows one to study the time-harmonic propagation of electromagnetic fields and, thus, to calculate the transmission coefficient of finite clusters in a natural way. Moreover, the advantages of using this real space method for structures of arbitrary complexity are also discussed. In addition, point defect cluster quality factor calculations are reported by means of FEM and they are compared to the ones obtained with the FDTD and \texttt{Harminv} methods. As a result, FEM comes out as an effective, stable, robust, and rigorous tool to study light propagation and confinement in both periodic and aperiodic dielectric photonic crystals.

Mesoscale and Nanoscale Physics,Optics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to study the photonic band structure, transmission coefficient and quality factor of two - dimensional photonic crystals through the finite element method (FEM). Specifically, the author hopes to calculate and analyze the following problems through FEM:

1. **Photonic band structure**: Calculate the photonic band structures of different two - dimensional periodic and non - periodic dielectric photonic crystals. This includes a detailed eigenvalue analysis of photonic crystals with square and triangular lattices.

2. **Transmission coefficient**: Study the transmission coefficient of finite - length photonic crystal clusters, especially for the natural calculation method of time - harmonic electromagnetic field propagation.

3. **Quality factor**: Calculate the quality factor of point - defect clusters and compare it with the results of the FDTD and Harminv methods.

4. **Applicability to arbitrary complex structures**: Explore the advantages of FEM in dealing with arbitrary complex geometric structures, especially its performance when dealing with frequency - dependent dielectric functions and discontinuous dielectric functions.

5. **Research on disordered photonic crystals**: Evaluate the feasibility of FEM in studying disordered photonic crystals, especially for cases of random displacement or component disorder caused by the manufacturing process itself.

### Specific problem description

- **Photonic band structure**: Calculate the photonic band structure through FEM and verify whether its results are consistent with the traditional plane - wave method.

\[
\nabla^2 E_z(\vec{r})+k_0^2 \epsilon_r(\vec{r}) E_z(\vec{r}) = 0
\]

where \(E_z(\vec{r})\) is the z - component of the electric field at position \(\vec{r}\), \(\epsilon_r(\vec{r})\) is the relative dielectric constant of the photonic crystal, \(k_0=\frac{\omega}{c}\), \(\omega\) is the angular frequency of the incident electric field, and \(c\) is the speed of light in free space.

- **Transmission coefficient**: Study the transmission coefficient of finite - length photonic crystal clusters, especially the setting of boundary conditions for TE and TM polarization cases.

\[
\hat{n} \times \vec{H}=\vec{0} \quad \text{(TE polarization)}
\]
\[
\hat{n} \times \vec{E}=\vec{0} \quad \text{(TM polarization)}
\]

- **Quality factor**: Calculate the quality factor of point - defect clusters and compare it with the results of the FDTD and Harminv methods.

### Conclusion

The main purpose of the paper is to evaluate the effectiveness and superiority of FEM in studying two - dimensional photonic crystals, especially its advantages in dealing with complex geometric structures and disordered systems. By comparing with traditional methods, the author proves that FEM is an effective, stable, robust and rigorous tool that can be used to study light propagation and confinement in photonic crystals.