Numerical investigation of the flat band Bloch modes in a 2D photonic crystal with Dirac cones

A numerical method combining complex-k band calculations and absorbing boundary conditions for Bloch waves is presented. We use this method to study photonic crystals with Dirac cones. We demonstrate that the photonic crystal behaves as a zero-index medium when excited at normal incidence, but that the zero-index behavior is lost at oblique incidence due to excitation of modes on the flat band. We also investigate the formation of monomodal and multimodal cavity resonances inside the photonic crystals, and the physical origins of their different line-shape features.