**Ψ**existing in its bulk while its surface is charge free and current free. An alternative model is obtained by setting

**Ψ≡0**and incorporating surface charge and current densities characterized by an admittance

**γ**. Analysis of planewave reflection and refraction due to a topological-insulator half space reveals that the parameters

**Ψ**and

**γ**arise identically in the reflection and transmission coefficients, implying that the two classical models cannot be distinguished on the basis of any scattering scenario. However, as

**Ψ**disappears from the Maxwell equations applicable to any region occupied by the topological insulator, and because surface states exist on topological insulators as protected conducting states, the alternative model must be chosen.

## 1.

## Introduction

The discovery of topological insulators^{1} has prompted researchers in classical optics^{2}3.^{–}^{4} to examine electromagnetic scattering due to bound objects made of these materials, exemplified by chalcogenides such as ${\mathrm{Bi}}_{2}{\mathrm{Se}}_{3}$, ${\mathrm{Bi}}_{2}{\mathrm{Te}}_{3}$, and ${\mathrm{Sb}}_{2}{\mathrm{Te}}_{3}$. As a topological insulator is considered to be an isotropic material, its frequency-domain constitutive relations are formulated to contain a magnetoelectric pseudoscalar (denoted by $\mathrm{\Psi}$ here) in addition to the permittivity scalar $\u03f5$ and the permeability scalar $\mu $. The surface of the topological insulator is assumed to be charge free and current free, and the scattering problem can then be solved by following textbook techniques.^{5}

Yet, according to condensed-matter theory, surface states exist on topological insulators as protected conducting states,^{2} and the characteristic electromagnetic responses of these materials are due to those surface states. Should then a topological insulator’s optical response be modeled as due solely to either

I. the bulk constitutive parameter $\mathrm{\Psi}$ with the surface of the topological insulator being charge free and current free, or

II. a surface parameter (denoted by $\gamma $ here) that quantifies the charge density and current density on the surface of the topological insulator?

The topological insulator possesses the permittivity $\u03f5$ and permeability $\mu $ in both models. Both $\mathrm{\Psi}$ and $\gamma $ are admittances, and whereas the magnetoelectric constitutive parameter $\mathrm{\Psi}$ mediates between $\mathbf{D}$ and $\mathbf{B}$ as well as between $\mathbf{H}$ and $\mathbf{E}$ throughout the topological insulator, $\gamma $ is meaningful only on the surface of that material.

This communication is devoted to a comparison of models I and II, through the fundamental boundary-value problem of reflection and refraction of a plane wave. This problem is described and solved in Sec. 2 for both models. Section 3 contains a comparative discussion of the two models. Vectors are underlined. An $\mathrm{exp}(-i\omega t)$ dependence on time $t$ is implicit, with $\omega $ as the angular frequency and $i=\sqrt{-1}$.

## 2.

## Fundamental Boundary-Value Problem

Suppose that all space is divided into two mutually disjoint half spaces ${\mathcal{V}}_{\text{out}}=\{(x,y,z):z<0\}$ and ${\mathcal{V}}_{\text{in}}=\{(x,y,z):z>0\}$ separated by the surface $\mathcal{S}=\{(x,y,z):z=0\}$. We need to solve the frequency-domain macroscopic Maxwell equations:

## (1)

$$\begin{array}{l}\nabla \xb7\mathbf{B}(\mathbf{r},\omega )=0\\ \nabla \times \mathbf{E}(\mathbf{r},\omega )-i\omega \mathbf{B}(\mathbf{r},\omega )=\mathbf{0}\\ \nabla \xb7\mathbf{D}(\mathbf{r},\omega )=0\\ \nabla \times \mathbf{H}(\mathbf{r},\omega )+i\omega \mathbf{D}(\mathbf{r},\omega )=\mathbf{0}\end{array}\}$$Let the half space ${\mathcal{V}}_{\text{out}}$ be vacuous so that the constitutive equations

## (2)

$$\mathbf{D}(\mathbf{r},\omega )={\u03f5}_{0}\mathbf{E}(\mathbf{r},\omega ),\phantom{\rule[-0.0ex]{1em}{0.0ex}}\mathbf{H}(\mathbf{r},\omega )={\mu}_{0}^{-1}\mathbf{B}(\mathbf{r},\omega ),\phantom{\rule[-0.0ex]{1em}{0.0ex}}\mathbf{r}\in {\mathcal{V}}_{\text{out}},$$## (3)

$$\begin{array}{l}\nabla \xb7\mathbf{B}(\mathbf{r},\omega )=0\\ \nabla \times \mathbf{E}(\mathbf{r},\omega )-i\omega \mathbf{B}(\mathbf{r},\omega )=\mathbf{0}\\ \nabla \xb7\mathbf{E}(\mathbf{r},\omega )=0\\ \nabla \times \mathbf{B}(\mathbf{r},\omega )+i\omega {\mu}_{0}{\u03f5}_{0}\mathbf{E}(\mathbf{r},\omega )=\mathbf{0}\end{array}\},\phantom{\rule[-0.0ex]{1em}{0.0ex}}\mathbf{r}\in {\mathcal{V}}_{\text{out}}$$The frequency-domain constitutive relations of the material occupying ${\mathcal{V}}_{\text{in}}$ are

## (4)

$$\begin{array}{l}\mathbf{D}(\mathbf{r},\omega )=\u03f5(\omega )\mathbf{E}(\mathbf{r},\omega )+\mathrm{\Psi}(\omega )\mathbf{B}(\mathbf{r},\omega )\\ \mathbf{H}(\mathbf{r},\omega )={\mu}^{-1}(\omega )\mathbf{B}(\mathbf{r},\omega )-\mathrm{\Psi}(\omega )\mathbf{E}(\mathbf{r},\omega )\end{array}\},\phantom{\rule[-0.0ex]{1em}{0.0ex}}\mathbf{r}\in {\mathcal{V}}_{\text{in}},$$## (5)

$$\begin{array}{l}\nabla \xb7\mathbf{B}(\mathbf{r},\omega )=0\\ \nabla \times \mathbf{E}(\mathbf{r},\omega )-i\omega \mathbf{B}(\mathbf{r},\omega )=\mathbf{0}\\ \nabla \xb7\mathbf{E}(\mathbf{r},\omega )=0\\ \nabla \times \mathbf{B}(\mathbf{r},\omega )+i\omega \mu (\omega )\u03f5(\omega )\mathbf{E}(\mathbf{r},\omega )=\mathbf{0}\end{array}\},\phantom{\rule[-0.0ex]{1em}{0.0ex}}\mathbf{r}\in {\mathcal{V}}_{\text{in}}.$$When solving an electromagnetic boundary-value problem, it is common to use the boundary conditions:

## (6)

$$\begin{array}{l}\widehat{\mathbf{n}}({\mathbf{r}}_{\mathrm{S}})\xb7[{\mathbf{B}}_{\text{out}}({\mathbf{r}}_{\mathrm{S}},\omega )-{\mathbf{B}}_{\text{in}}({\mathbf{r}}_{\mathrm{S}},\omega )]=0\\ \widehat{\mathbf{n}}({\mathbf{r}}_{\mathrm{S}})\times [{\mathbf{E}}_{\text{out}}({\mathbf{r}}_{\mathrm{S}},\omega )-{\mathbf{E}}_{\text{in}}({\mathbf{r}}_{\mathrm{S}},\omega )]=\mathbf{0}\\ \widehat{\mathbf{n}}({\mathbf{r}}_{\mathrm{S}})\xb7[{\mathbf{D}}_{\text{out}}({\mathbf{r}}_{\mathrm{S}},\omega )-{\mathbf{D}}_{\text{in}}({\mathbf{r}}_{\mathrm{S}},\omega )]={\rho}_{s}({\mathbf{r}}_{\mathrm{S}},\omega )\\ \widehat{\mathbf{n}}({\mathbf{r}}_{\mathrm{S}})\times [{\mathbf{H}}_{\text{out}}({\mathbf{r}}_{\mathrm{S}},\omega )-{\mathbf{H}}_{\text{in}}({\mathbf{r}}_{\mathrm{S}},\omega )]={\mathbf{J}}_{s}({\mathbf{r}}_{\mathrm{S}},\omega )\end{array}\},\phantom{\rule[-0.0ex]{1em}{0.0ex}}{\mathbf{r}}_{\mathrm{S}}\in \mathcal{S},$$## (7)

$$\begin{array}{l}{\rho}_{s}({\mathbf{r}}_{\mathrm{S}},\omega )=\gamma (\omega )\widehat{\mathbf{n}}({\mathbf{r}}_{\mathrm{S}})\xb7{\mathbf{B}}_{\text{in}}({\mathbf{r}}_{\mathrm{S}},\omega )\\ {\mathbf{J}}_{s}({\mathbf{r}}_{\mathrm{S}},\omega )=-\gamma (\omega )\widehat{\mathbf{n}}({\mathbf{r}}_{\mathrm{S}})\times {\mathbf{E}}_{\text{in}}({\mathbf{r}}_{\mathrm{S}},\omega )\end{array}\},\phantom{\rule[-0.0ex]{1em}{0.0ex}}{\mathbf{r}}_{\mathrm{S}}\in \mathcal{S},$$Let an arbitrarily polarized plane wave in ${\mathcal{V}}_{\text{out}}$ be incident on $\mathcal{S}$. Then the primitive field phasors in ${\mathcal{V}}_{\text{out}}$ can be written as

## (8)

$$\begin{array}{l}\mathbf{E}(\mathbf{r},\omega )=\{[{a}_{s}{\widehat{\mathbf{u}}}_{y}+{a}_{p}(-{\widehat{\mathbf{u}}}_{x}{\tau}_{0}+{\widehat{\mathbf{u}}}_{z}\kappa )/{k}_{0}]\mathrm{exp}(i{\tau}_{0}z)\\ +[{r}_{s}{\widehat{\mathbf{u}}}_{y}+{r}_{p}({\widehat{\mathbf{u}}}_{x}{\tau}_{0}+{\widehat{\mathbf{u}}}_{z}\kappa )/{k}_{0}]\mathrm{exp}(-i{\tau}_{0}z)\}\mathrm{exp}(i\kappa x)\\ \mathbf{B}(\mathbf{r},\omega )=\frac{{k}_{0}}{\omega}\{[-{a}_{p}{\widehat{\mathbf{u}}}_{y}+{a}_{s}(-{\widehat{\mathbf{u}}}_{x}{\tau}_{0}+{\widehat{\mathbf{u}}}_{z}\kappa )/{k}_{0}]\mathrm{exp}(i{\tau}_{0}z)\\ +[-{r}_{p}{\widehat{\mathbf{u}}}_{y}+{r}_{s}({\widehat{\mathbf{u}}}_{x}{\tau}_{0}+{\widehat{\mathbf{u}}}_{z}\kappa )/{k}_{0}]\mathrm{exp}(-i{\tau}_{0}z)\}\mathrm{exp}(i\kappa x)\end{array}\},\phantom{\rule[-0.0ex]{1em}{0.0ex}}\mathbf{r}\in {\mathcal{V}}_{\text{out}},$$The primitive field phasors in ${\mathcal{V}}_{\text{in}}$ are given as

## (9)

$$\begin{array}{l}\mathbf{E}(\mathbf{r},\omega )=[{t}_{s}{\widehat{\mathbf{u}}}_{y}+{t}_{p}(-{\widehat{\mathbf{u}}}_{x}\tau +{\widehat{\mathbf{u}}}_{z}\kappa )/k]\mathrm{exp}(i\tau z)\mathrm{exp}(i\kappa x)\\ \mathbf{B}(\mathbf{r},\omega )=\frac{k}{\omega}[-{t}_{p}{\widehat{\mathbf{u}}}_{y}+{t}_{s}(-{\widehat{\mathbf{u}}}_{x}\tau +{\widehat{\mathbf{u}}}_{z}\kappa )/k]\mathrm{exp}(i\tau z)\mathrm{exp}(i\kappa x)\end{array}\},\phantom{\rule[-0.0ex]{1em}{0.0ex}}\mathbf{r}\in {\mathcal{V}}_{\text{in}},$$The foregoing expressions were substituted into Eqs. (2)_{2}, (4)_{2}, (6)_{2,4} and (7)_{2} to determine ${r}_{s}$, ${r}_{p}$, ${t}_{s}$, and ${t}_{p}$ in terms of ${a}_{s}$ and ${a}_{p}$. Thus,

## (10)

$${r}_{s}=\frac{[({\eta}_{r}-{\delta}_{r})(1+{\eta}_{r}{\delta}_{r})-{(G{\eta}_{0})}^{2}{\eta}_{r}^{2}{\delta}_{r}]{a}_{s}+2\text{\hspace{0.17em}\hspace{0.17em}}G{\eta}_{0}{\eta}_{r}^{2}{\delta}_{r}{a}_{p}}{({\eta}_{r}+{\delta}_{r})(1+{\eta}_{r}{\delta}_{r})+{(G{\eta}_{0})}^{2}{\eta}_{r}^{2}{\delta}_{r}},$$## (11)

$${r}_{p}=\frac{[({\eta}_{r}+{\delta}_{r})(1-{\eta}_{r}{\delta}_{r})+{(G{\eta}_{0})}^{2}{\eta}_{r}^{2}{\delta}_{r}]{a}_{p}+2\text{\hspace{0.17em}\hspace{0.17em}}G{\eta}_{0}{\eta}_{r}^{2}{\delta}_{r}{a}_{s}}{({\eta}_{r}+{\delta}_{r})(1+{\eta}_{r}{\delta}_{r})+{(G{\eta}_{0})}^{2}{\eta}_{r}^{2}{\delta}_{r}},$$## (12)

$${t}_{s}=\frac{2{\eta}_{r}(1+{\eta}_{r}{\delta}_{r}){a}_{s}+2\text{\hspace{0.17em}\hspace{0.17em}}G{\eta}_{0}{\eta}_{r}^{2}{\delta}_{r}{a}_{p}}{({\eta}_{r}+{\delta}_{r})(1+{\eta}_{r}{\delta}_{r})+{(G{\eta}_{0})}^{2}{\eta}_{r}^{2}{\delta}_{r}},$$## (13)

$${t}_{p}=\frac{2{\eta}_{r}({\eta}_{r}+{\delta}_{r}){a}_{p}-2\text{\hspace{0.17em}\hspace{0.17em}}G{\eta}_{0}{\eta}_{r}^{2}{a}_{s}}{({\eta}_{r}+{\delta}_{r})(1+{\eta}_{r}{\delta}_{r})+{(G{\eta}_{0})}^{2}{\eta}_{r}^{2}{\delta}_{r}},$$## (14)

$$G=\mathrm{\Psi}+\gamma ,\phantom{\rule[-0.0ex]{1em}{0.0ex}}{\eta}_{0}=\sqrt{\frac{{\mu}_{0}}{{\u03f5}_{0}}},\phantom{\rule[-0.0ex]{1em}{0.0ex}}{\delta}_{r}=\frac{\tau /k}{{\tau}_{0}/{k}_{0}},\phantom{\rule[-0.0ex]{1em}{0.0ex}}{\eta}_{r}=\sqrt{\frac{{\u03f5}_{0}\mu}{\u03f5{\mu}_{0}}}.$$We have verified that Eqs. (10) through (13) satisfy Eqs. (6)_{1,3} and (7)_{1}. Moreover, Eqs. (10) through (13) simplify to the standard results:^{6}^{,}^{7}

## (15)

$$\begin{array}{rr}{r}_{s}={a}_{s}\frac{{\eta}_{r}-{\delta}_{r}}{{\eta}_{r}+{\delta}_{r}},& {r}_{p}={a}_{p}\frac{1-{\eta}_{r}{\delta}_{r}}{1+{\eta}_{r}{\delta}_{r}}\\ {t}_{s}={a}_{s}\frac{2{\eta}_{r}}{{\eta}_{r}+{\delta}_{r}},& {t}_{p}={a}_{p}\frac{2{\eta}_{r}}{1+{\eta}_{r}{\delta}_{r}}\end{array}\}$$## 3.

## Discussion and Conclusion

Equations (10) through (13) can be recast in matrix form as

## (16)

$$\left[\begin{array}{c}{r}_{s}\\ {r}_{p}\end{array}\right]=\left[\begin{array}{cc}{r}_{ss}& {r}_{sp}\\ {r}_{ps}& {r}_{pp}\end{array}\right]\left[\begin{array}{c}{a}_{s}\\ {a}_{p}\end{array}\right],\phantom{\rule[-0.0ex]{1em}{0.0ex}}\left[\begin{array}{c}{t}_{s}\\ {t}_{p}\end{array}\right]=\left[\begin{array}{cc}{t}_{ss}& {t}_{sp}\\ {t}_{ps}& {t}_{pp}\end{array}\right]\left[\begin{array}{c}{a}_{s}\\ {a}_{p}\end{array}\right].$$Equations (10) through (13) do not contain $\mathrm{\Psi}$ and $\gamma $ separately, but their sum $G$ instead. Thus, measurements of the reflection coefficients ${r}_{s}$ and ${r}_{p}$ (or the transmission coefficients ${t}_{s}$ and ${t}_{p}$, if at all possible) cannot be used to discriminate between models I ($\gamma =0$) and II ($\mathrm{\Psi}=0$). Equations (6)_{4} and (7)_{2} together make it clear that measurements of the reflection and transmission coefficients of a slab made of a topological insulator cannot be used to discriminate between the two models; not only that, the solution of every scattering problem will depend on $G$, not on $\mathrm{\Psi}$ alone or $\gamma $ alone.

This impasse can be resolved by realizing that surface states exist on topological insulators as protected conducting states, and the characteristic behavior of these materials is due to those surface states. Furthermore, $\mathrm{\Psi}$ vanishes from the Maxwell equations (5) applicable to ${\mathcal{V}}_{\text{in}}$ occupied by the topological insulator; indeed, $\mathrm{\Psi}$ would vanish even if the topological insulator were bianisotropic.^{8} For both of these reasons, we must choose model II, which also satisfies the Post constraint $\mathrm{\Psi}\equiv 0$.^{9}

As the material occupying ${\mathcal{V}}_{\text{in}}$ is isotropic and achiral, cross-polarized reflection in this problem has been taken to arise from the Lorentz nonreciprocity inherent in Eq. (4).^{10} But now we see that surface states described by Eq. (7) by themselves are capable of yielding cross-polarized reflection, which is, therefore, not an indication of Lorentz nonreciprocity.

## Acknowledgments

A.L. is grateful to the Charles Godfrey Binder Endowment at Penn State for ongoing support of his research. T.G.M. acknowledges the support of EPSRC grant EP/M018075/1.

## References

## Biography

**Akhlesh Lakhtakia** received his degrees from the Banaras Hindu University and the University of Utah. He is the Charles Godfrey Binder Professor of engineering science and mechanics at the Pennsylvania State University. His current research interests include nanotechnology, bioreplication, forensic science, solar-energy harvesting, surface multiplasmonics, metamaterials, mimumes, and sculptured thin films. He is a fellow of OSA, SPIE, IoP, AAAS, APS, and IEEE. He received the 2010 SPIE Technical Achievement Award.

**Tom G. Mackay** is a reader in applied mathematics at the University of Edinburgh and an adjunct professor in the Department of Engineering Science and Mechanics at the Pennsylvania State University. He is a graduate of the Universities of Edinburgh, Glasgow, and Strathclyde, and a fellow of the Institute of Physics (UK) and SPIE. His research interests include the electromagnetic theory of novel and complex materials, including homogenized composite materials.