High contrast elliptic operators in honeycomb structures

TL;DR

This paper analyzes the band structure of high contrast elliptic operators with honeycomb symmetry, revealing the existence and behavior of Dirac points which are crucial for understanding phenomena in photonic and electronic materials.

Contribution

It extends the analysis of Dirac points and band structures to high contrast elliptic operators with honeycomb symmetry, including asymptotic expansions and error bounds.

Findings

Existence of Dirac points in low and high energy bands.

Asymptotic expansions of eigenpairs and Dirac velocity in high contrast limit.

Differences between high contrast elliptic operators and strong binding Schrödinger regimes.

Abstract

We study the band structure of self-adjoint elliptic operators $\mathbb{A}_g= -\nabla \cdot \sigma_{g} \nabla$, where $\sigma_g$ has the symmetries of a honeycomb tiling of $\mathbb{R}^2$. We focus on the case where $\sigma_{g}$ is a real-valued scalar: $\sigma_{g}=1$ within identical, disjoint "inclusions", centered at vertices of a honeycomb lattice, and $\sigma_{g}=g \gg1 $ (high contrast) in the complement of the inclusion set (bulk). Such operators govern, e.g. transverse electric (TE) modes in photonic crystal media consisting of high dielectric constant inclusions (semi-conductor pillars) within a homogeneous lower contrast bulk (air), a configuration used in many physical studies. Our approach, which is based on monotonicity properties of the associated energy form, extends to a class of high contrast elliptic operators that model heterogeneous and anisotropic honeycomb media. Our results concern the global behavior of dispersion surfaces, and the existence of conical crossings (Dirac points) occurring in the lowest two energy bands as well as in bands arbitrarily high in the spectrum. Dirac points are the source of important phenomena in fundamental and applied physics, e.g. graphene and its artificial analogues, and topological insulators. The key hypotheses are the non-vanishing of the Dirac (Fermi) velocity $v_D(g)$, verified numerically, and a spectral isolation condition, verified analytically in many configurations. Asymptotic expansions, to any order in $g^{-1}$, of Dirac point eigenpairs and $v_D(g)$ are derived with error bounds. Our study illuminates differences between the high contrast behavior of $\mathbb{A}_g$ and the corresponding strong binding regime for Schroedinger operators.