Universality of mode symmetries in creating photonic Dirac cones

doi:10.1364/JOSAB.29.002770

« Show journal navigation

Universality of mode symmetries in creating photonic Dirac cones

Kazuaki Sakoda »View Author Affiliations

JOSA B, Vol. 29, Issue 10, pp. 2770-2778 (2012)
http://dx.doi.org/10.1364/JOSAB.29.002770

View Full Text Article

Acrobat PDF (661 KB)

Journal Search
Article Lookup

Article Tools

Citations

Alert me when this article is cited
Export Citation/Save

Article Navigation

▸ Article Outline
– Abstract
INTRODUCTION
THEORY
RESULTS AND DISCUSSION
CONCLUSION
– References and links

Article Tools
Full Text
Article Info
References
Cited By
Figures
Metrics
Download PDF

Viewing Equations in the
HTML Article

Optimal Viewing Recommendations

Full Text
Article Info
References (14)
Cited By
Figures (7)
Metrics

Abstract

The creation of photonic Dirac cones by accidental degeneracy in the Brillouin-zone center was recently reported for both metamaterials with localized electromagnetic resonant states and dielectric photonic crystals without well-defined resonance. Based on the anticipation that there should be a common physical origin in this phenomenon, we systematically examined the relation between mode symmetries and shapes of dispersion curves for both systems. The result strongly suggests the presence of universality of mode symmetries that enable the creation of photonic Dirac cones irrespective of the details of the sample structure.

1. INTRODUCTION

Photonic Dirac cones, or the isotropic linear dispersion relation around certain points in the Brillouin zone, have been attracting considerable interest during the last five years, since they offer various novel possibilities for physics and engineering. Haldane et al. [1

F. D. M. Haldane and S. Raghu, “Possible realization of directional optical waveguides in photonic crystals with broken time-reversal symmetry,” Phys. Rev. Lett. 100, 013904 (2008). [CrossRef]

S. Raghu and F. D. M. Haldane, “Analogs of quantum-Hall-effect edge states in photonic crystals,” Phys. Rev. A 78, 033834 (2008). [CrossRef]

] pointed out the presence of photonic Dirac cones on the Brillouin-zone boundary of two-dimensional triangular-lattice photonic crystals due to their structural symmetry and discussed unidirectional propagation of surface modes caused by time-reversal symmetry breaking due to the Faraday effect. Ochiai and Onoda extended the discussion to honeycomb-lattice photonic crystals [3

T. Ochiai and M. Onoda, “Photonic analog of graphene model and its extension: Dirac cone, symmetry, and edge states,” Phys. Rev. B 80, 155103 (2009). [CrossRef]

]. Zhang proposed optical simulation of Zitterbewegung, or trembling motion, in particle physics by propagating an optical pulse of the Dirac point frequency [4

X. Zhang, “Observing zitterbewegung for photons near the Dirac point of a two-dimensional photonic crystal,” Phys. Rev. Lett. 100, 113903 (2008). [CrossRef]

]. A pseudodiffusive transmission whose intensity was inversely proportional to the thickness of the photonic crystal was found by Sepkhanov et al. [5

R. A. Sepkhanov, Y. B. Bazaliy, and C. W. J. Beenakker, “Extremal transmission at the Dirac point of a photonic band structure,” Phys. Rev. A 75, 063813 (2007). [CrossRef]

] and numerically demonstrated by Diem et al. [6

M. Diem, T. Koschny, and C. M. Soukoulis, “Transmission in the vicinity of the Dirac point in hexagonal photonic crystals,” Physica B 405, 2990–2995 (2010). [CrossRef]

Recently, Nam et al. analyzed the wave propagation in a periodic metal-dielectric waveguide array, which may be regarded as a one-dimensional photonic crystal, to show the presence of the linear dispersion relation, or a one-dimensional Dirac cone, in the Brillouin-zone center by adjusting the coupling constants between the waveguides [7

S. H. Nam, A. J. Taylor, and A. Efimov, “Diabolical point and conical-like diffraction in periodic plasmonic nanostructures,” Opt. Express 18, 10120–10126 (2010). [CrossRef]

]. Huang et al. reported that Dirac cones can also be created in the Brillouin-zone center of two-dimensional dielectric photonic crystals by accidental degeneracy of two bands [8

X. Huang, Y. Lai, Z. H. Hang, H. Zheng, and C. T. Chan, “Dirac cones induced by accidental degeneracy in photonic crystals and zero-refractive-index materials,” Nat. Mater. 10, 582–586 (2011). [CrossRef]

]. They showed for square- and triangular-lattice photonic crystals that combinations of a nondegenerate band and a doubly degenerate band yield a Dirac cone together with an auxiliary quadratic dispersion surface. Because the Dirac point in the Brillouin-zone center is equivalent to a zero effective refractive index [8

], it has much potential for various applications such as scatter-free waveguides [9

M. Silveirinha and N. Engheta, “Tunneling of electromagnetic energy through subwavelength channels and bends using $ε$ -near-zero materials,” Phys. Rev. Lett. 97, 157403 (2006). [CrossRef]

] and lenses of arbitrary shapes [10

A. Alu, M. G. Silveirinha, A. Salandrino, and N. Engheta, “Epsilon-near-zero metamaterials and electromagnetic sources: Tailoring the radiation phase pattern,” Phys. Rev. B 75, 155410 (2007). [CrossRef]

On the other hand, we showed by tight-binding approximation and group theory that we can also create Dirac cones in the Brillouin-zone center of metamaterials, which are characterized by well-defined electromagnetic resonant states localized in their unit structures, by accidental degeneracy of two modes [11

K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef]

,12

K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]

]. We proved that the combination of

A_{1}

and

E

modes of square-lattice metamaterials and the combination of

A_{1 g}

and

T_{1 u}

modes of simple-cubic-lattice metamaterials create a Dirac cone with auxiliary quadratic dispersion surfaces [11

K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef]

]. We further proved that the combination of

E_{1}

and

E_{2}

modes of triangular-lattice metamaterials yields double Dirac cones, or a pair of identical Dirac cones, without an auxiliary quadratic dispersion surface [12

K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]

Because the presence of Dirac cones was found in two extreme systems, i.e., dielectric photonic crystals for which the extended-wave picture is appropriate and metamaterials with localized electromagnetic resonant states for which the tight-binding picture is appropriate, it is highly likely that there is a common physical origin in this phenomenon irrespective of the details of the sample structure.

In this paper, based on this anticipation, we systematically examine the relation between the symmetry of modes and the shapes of dispersion curves generated by accidental degeneracy for both dielectric photonic crystals and metamaterials with localized electromagnetic resonant states by numerical calculation and tight-binding approximation, respectively. We show that the two calculations give the same results, which strongly suggests the presence of universality of mode symmetries that enable the creation of photonic Dirac cones irrespective of the details of the sample structure, irrespective of the presence or absence of localized resonant states, and irrespective of the approximation used in the calculation.

This paper is organized as follows. In Section 2, the methods for calculating the dispersion curves of metamaterials and photonic crystals are briefly described. In Section 3, the shapes of the dispersion curves thus calculated are presented in the case of accidental degeneracy for various combinations of mode symmetries. A summary of the present study is given in Section 4. In Appendix A, we show the formulation of the tight-binding approximation with the magnetic fields of the resonant states of a single unit structure of metamaterials as basis functions. The method for deriving mutual relations among different electromagnetic transfer integrals is also described. The shapes of the dispersion curves thus calculated for metamaterials of the

C_{6 v}

and

C_{4 v}

symmetries are presented in Appendices B and C, respectively.

2. THEORY

For the analytical description of dispersion curves of periodic metamaterials by tight-binding approximation, electromagnetic transfer integrals defined by the following equation play the main role [11

K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef]

,12

K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]

L_{m}^{(i j)} \equiv \frac{1}{V} \int_{V} d r H^{(i) *} (r) \cdot L H^{(j)} (r - r_{m}),

(1)

where

{H^{(i)}}

(

i = 1, 2 \dots

) denote the magnetic field of resonant states of a unit cell,

{r_{m}}

are lattice vectors,

V

is the volume on which the periodic boundary condition is imposed, and

L

is the operator that defines the eigenvalue problem of the electromagnetic field:

L H_{k} \equiv \nabla \times [\frac{1}{ε (r)} \nabla \times H_{k}] = \frac{ω_{k}^{2}}{c^{2}} H_{k},

(2)

where

ε (r)

is the position-dependent dielectric constant,

ω_{k}

is the eigen angular frequency, and

c

is the light velocity in free space. The magnetic permeability is assumed to be unity, since we do not deal with magnetic materials.

We analytically solve the secular equation of the electromagnetic field for wave vector

{k}

in the vicinity of the Brillouin-zone center by taking into account

L_{m}^{(i j)}

’s for the nearest-neighbor lattice points [11

K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef]

,12

K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]

]. The spatial symmetry of the resonant states results in many relations among

L_{m}^{(i j)}

’s and simplifies the secular equation (see Appendix A). In the case of accidental degeneracy of

E_{1}

and

E_{2}

modes of triangular-lattice metamaterials with

C_{6 v}

(regular triangle) symmetry, for example, we can thus analytically solve the secular equation that is quartic in

ω^{2}

[12

K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]

]. For numerical calculations of the dispersion relation of photonic crystals, a conventional plane-wave method is used [13

K. Sakoda, Optical Properties of Photonic Crystals , 2nd ed. (Springer-Verlag, 2004).

3. RESULTS AND DISCUSSION

First, we examine the triangular lattice of

C_{6 v}

symmetry. In Table 1, the shapes of dispersion curves generated by accidental degeneracy on the

Γ

point of the Brillouin zone are summarized for metamaterials. As we mentioned previously, we obtain double Dirac cones by combining

E_{1}

and

E_{2}

modes [12

K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]

]. On the other hand, we found in the present study that (

A_{1}

E_{1}

), (

A_{2}

E_{1}

), and (

B_{2}

E_{2}

) combinations yield a Dirac cone with an auxiliary quadratic dispersion surface. For the rest of the mode combinations listed in Table 1, we obtain two or three different quadratic dispersion curves (see Appendix B for details).

Table 1. Types of Dispersion Curves Generated by Accidental Degeneracy of Two Modes (Mode 1 and Mode 2) for Triangular-Lattice Metamaterials of

C_{6 v}

Symmetry^a

Mode 1	Mode 2	Dispersion curves	Refs.
$A_{1}$	$B_{1}$	$Q + Q$	B.1
$A_{1}$	$B_{2}$	$Q + Q$	B.2
$A_{1}$	$E_{1}$	$D + Q$	B.3
$A_{1}$	$E_{2}$	$Q + Q + Q$	B.8
$A_{2}$	$E_{1}$	$D + Q$	B.5
$A_{2}$	$E_{2}$	$Q + Q + Q$	B.7
$B_{2}$	$E_{1}$	$Q + Q + Q$	B.6
$B_{2}$	$E_{2}$	$D + Q$	B.4
$E_{1}$	$E_{2}$	$D D$	[12 K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef] ]

^a Symmetries of the magnetic field of the two modes on the

Γ

point are given in the left column. Shapes of the dispersion curves are given in the middle column, where

D

D D

, and

Q

denote Dirac cone, double Dirac cones, and quadratic dispersion surface, respectively. B.1–B.8 and the number in brackets in the right column denote the section number of the appendix and the reference number, respectively, where the shapes of the dispersion curves are derived by tight-binding approximation.

For the photonic crystal, we assume a periodic triangular array of circular air cylinders in a uniform material with a dielectric constant of 12.6 (GaAs). First, we examine the structural parameters for materializing accidental degeneracy. Figure 1 shows the second to fifth lowest eigenfrequencies on the

Γ

point for transverse electric (TE) polarization with its electric field perpendicular to the cylinder axis. In this figure, the vertical axis denotes the eigen angular frequency

ω

normalized by the lattice constant of the photonic crystal

a

and the light velocity in free space, while the horizontal axis denotes the radius of the air cylinder

ρ

normalized by the lattice constant. The spatial symmetry of the magnetic field is also given to each mode.

Fig. 1. Normalized eigenfrequency (

ω a / 2 π c

) for the TE polarization on the

Γ

point of the photonic crystal composed of a regular triangular array of circular air cylinders in a medium with a dielectric constant of 12.6 (GaAs). The horizontal axis is the normalized radius of the air cylinders (

ρ / a

). The degeneracy points denoted by

a

e

correspond to the five panels in Fig. 2.

Figure 1 shows that there are six values of the normalized radius that materialize accidental degeneracy of eigenfrequencies. Dispersion curves for five of the six cases are given in Figs. 2(a)–2(e). They are calculated for

k

in the

Γ

-to-

M

and

Γ

-to-

K

directions in the Brillouin zone [see Fig. 2(f)]. The sixth case (normalized radius of 0.494) is the degeneracy of the

E_{1}

and

E_{2}

modes and gives a feature similar to Fig. 2(a).

Fig. 2. Dispersion curves of the triangular-lattice photonic crystal for TE polarization with a normalized radius of (a) 0.407, (b) 0.433, (c) 0.459, (d) 0.474, and (e) 0.492. (f) Brillouin zone of the triangular lattice. The frequencies of accidental degeneracy are circled in (a)–(e).

M / 10

and

K / 10

in (e) mean that the horizontal axis is magnified by 10 times.

Accidental degeneracy of the

E_{1}

and

E_{2}

modes results in double Dirac cones. As shown in Fig. 2(a), their dispersion in the vicinity of the

Γ

point is linear and isotropic and the slopes of the two dispersion lines are the same on the

Γ

point. For the combination of

A_{1}

and

E_{1}

modes [Fig. 2(b)] and that of

B_{2}

and

E_{2}

modes [Fig. 2(e)], accidental degeneracy results in the formation of a Dirac cone with an auxiliary quadratic dispersion surface. Two of the three dispersion curves are linear and isotropic in the vicinity of the

Γ

point, while the third curve is quadratic. Finally, for the combination of

B_{2}

and

E_{1}

modes [Fig. 2(c)] and that of

A_{1}

and

E_{2}

modes [Fig. 2(d)], accidental degeneracy results in three quadratic dispersion curves. All these features are the same as the case of triangular-lattice metamaterials listed in Table 1, which was analytically proved assuming the presence of localized resonant states in each unit structure.

In Fig. 3, structural parameters for realizing accidental degeneracy are examined for transverse magnetic (TM) polarization with its magnetic field perpendicular to the cylinder axis. The spatial symmetry of the magnetic field is also given to each mode.

Fig. 3. The second to fifth lowest normalized eigenfrequency (

ω a / 2 π c

) for the TM polarization on the

Γ

point of the photonic crystal composed of a regular triangular array of circular air cylinders in a dielectric medium with

ε = 12.6

(GaAs). The horizontal axis is the normalized radius of the air cylinders (

ρ / a

). The degeneracy points denoted by

a

and

b

correspond to the two panels in Fig. 4.

Let us make one remark here. The TM polarization is often referred to as

E

polarization, and its electric field, which is parallel to the cylinder axis, is discussed traditionally. Because the electric field is a genuine vector and the magnetic field is an axial vector, their transformation by symmetry operations is generally different. When we denote the symmetry operation and its matrix representation by

R

and

R

and the character of transformation of the electric (magnetic) field by

χ^{(E)}

(

χ^{(H)}

), we have the following relation [13

K. Sakoda, Optical Properties of Photonic Crystals , 2nd ed. (Springer-Verlag, 2004).

χ^{(H)} (R) = (\det R) χ^{(E)} (R) .

(3)

Because of this difference between

χ^{(E)}

and

χ^{(H)}

, the irreducible representation of the electric field is generally different from that of the magnetic field. In the present study, we always refer to the symmetry of the magnetic field to avoid confusion.

There are two values of the normalized radius that materialize accidental degeneracy for the TM polarization as shown in Fig. 3. The combination of

A_{2}

and

E_{1}

modes [Fig. 4(a)] results in the creation of a Dirac cone with an auxiliary quadratic dispersion surface, whereas the combination of

A_{2}

and

E_{2}

modes [Fig. 4(b)] results in the formation of three quadratic dispersion curves. Both features agree with the case of metamaterials listed in Table 1.

Fig. 4. Dispersion curves of the triangular-lattice photonic crystal for TM polarization with a normalized radius of (a) 0.445, (b) 0.483. The frequencies of accidental degeneracy are circled.

Now we proceed to the case of square lattices of

C_{4 v}

(regular square) symmetry. Table 2 summarizes the shapes of dispersion curves of periodic metamaterials generated by accidental degeneracy on the

Γ

point. As we mentioned previously, we obtain a Dirac cone with an auxiliary quadratic dispersion surface by combining

A_{1}

and

E

modes [11

K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef]

], whereas (

A_{1}

B_{1}

) and (

A_{1}

B_{2}

) combinations yield two quadratic curves [12

K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]

]. We further found in the present study that (

B_{1}

E

) and (

B_{2}

E

) combinations give a Dirac cone and an auxiliary quadratic dispersion surface and that (

A_{1}

A_{2}

) and (

A_{2}

B_{2}

) combinations yield two quadratic dispersion curves (see Appendix C for details).

Table 2. Types of Dispersion Curves Generated by Accidental Degeneracy of Two Modes (Mode 1 and Mode 2) for Square-Lattice Metamaterials of

C_{4 v}

Symmetry^a

Mode 1	Mode 2	Dispersion Curves	Refs.
$A_{1}$	$A_{2}$	$Q + Q$	C.2
$A_{1}$	$B_{1}$	$Q + Q$	[12 K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef] ], C.1
$A_{1}$	$B_{2}$	$Q + Q$	[12 K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef] ], C.3
$A_{2}$	$B_{2}$	$Q + Q$	C.4
$A_{1}$	$E$	$D + Q$	[11 K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef] ]
$B_{1}$	$E$	$D + Q$	C.5
$B_{2}$	$E$	$D + Q$	C.6

^a Symmetries of the magnetic field of the two modes on the

Γ

point are given in the left column. Shapes of the dispersion curves are given in the middle column, where

D

and

Q

denote Dirac cone and quadratic dispersion surface, respectively. C.1–C.6 and the number in brackets in the right column denote the section number of the appendix and the reference number, respectively, where the shapes of the dispersion curves are derived by tight-binding approximation.

On the other hand, for square-lattice photonic crystals of the

C_{4 v}

symmetry composed of circular air cylinders in a uniform material with a dielectric constant of 12.6, we found three values of the normalized radius that materialize accidental degeneracy on the

Γ

point. Dispersion curves for these three cases are shown in Fig. 5. Their shapes agree with the metamaterials listed in Table 2.

Fig. 5. (a) Brillouin zone of the square lattice of

C_{4 v}

symmetry. Dispersion curves of the square-lattice photonic crystal (b) for TE polarization with a normalized radius of 0.299 and for TM polarization with a normalized radius of (c) 0.421 and (d) 0.464. The frequencies of accidental degeneracy are circled.

M / 2

and

X / 2

in (c) mean that the horizontal axis is magnified by two times.

Thus, in all cases that we could examine numerically, dielectric photonic crystals and metamaterials with localized electromagnetic resonant states have the same relation between the mode symmetries and the shapes of the dispersion curves generated by accidental degeneracy in the vicinity of the

Γ

point. This fact strongly suggests the presence of universality in combinations of mode symmetries that enable the creation of photonic Dirac cones and double Dirac cones irrespective of the details of the sample structure. This knowledge is quite useful for designing specimens with Dirac cones to realize many applications described in Section 1. However, we should note that the universality has not been proved rigorously yet, since the present study is partly based on numerical calculation for a limited number of examples. Three-dimensional systems in particular were not analyzed, although we expect a similar universality of mode symmetries. So, the rigorous proof of the universality remains a future problem.

4. CONCLUSION

We analytically solved secular equations of periodic metamaterials of the

C_{6 v}

and

C_{4 v}

symmetries with localized electromagnetic resonant states by tight-binding approximation and obtained the shapes of dispersion curves in the vicinity of the Brillouin-zone center. Thus we clarified the combinations of mode symmetries that materialize Dirac cones and double Dirac cones by accidental degeneracy. We also examined dispersion curves of dielectric photonic crystals numerically by the plane-wave expansion method. By changing a structural parameter of assumed samples, we found 11 cases of accidental degeneracy on the

Γ

point and compared the shapes of dispersion curves with the metamaterials.

In all 11 cases, the relation between the mode symmetries and the shape of the dispersion curves was common to metamaterials with localized electromagnetic resonant states and dielectric photonic crystals without well-defined resonant states, which strongly suggests the presence of universality of mode symmetries that enable the creation of photonic Dirac cones irrespective of the details of the sample structure, irrespective of the presence or absence of localized resonant states, and irrespective of the approximation used in the calculation.

Appendices

APPENDIX A: FORMULATION

Maxwell’s wave equation for the magnetic field

H (r, t)

is given by

L H \equiv \nabla \times [\frac{1}{ε (r)} \nabla \times H] = - \frac{1}{c^{2}} \frac{\partial^{2}}{\partial t^{2}} H,

(A1)

where

ε (r)

is the periodic dielectric constant of the system. The magnetic permeability is assumed to be unity, since we do not deal with magnetic materials. We impose the periodic boundary condition on

H

and assume that

ε (r)

is real to make our problem well-defined [13

K. Sakoda, Optical Properties of Photonic Crystals , 2nd ed. (Springer-Verlag, 2004).

]. Thus, Eq. (A1) leads to an eigenvalue problem.

For a single metallic unit structure described by dielectric constant

ε_{s} (r)

, we assume the presence of localized resonant states, which characterize the optical response of the metamaterial. Thus, their magnetic fields satisfy the following eigenvalue equation:

L_{s} H^{(i)} (r) \equiv \nabla \times [\frac{1}{ε_{s}} \nabla \times H^{(i)} (r)] = \frac{ω_{i}^{2}}{c^{2}} H^{(i)} (r),

(A2)

where operator

L_{s}

is defined by the first equality of this equation,

i

is the suffix to distinguish the resonant states, and

ω_{i}

is the resonance angular frequency. Without loss of generality, we can assume that the eigenfunctions are normalized as follows:

\int_{V} d r H^{(i) *} (r) \cdot H^{(j)} (r) = V δ_{i j},

(A3)

where

V

is the volume on which we impose the periodic boundary condition.

We assume that the unit cell of the periodic metamaterial has a certain spatial symmetry and denote the group of its symmetry operations by

G

. In Appendices B and C, we deal with two cases of

G = C_{6 v}

and

G = C_{4 v}

. By definition, the transformation of the magnetic field by symmetry operation

R \in G

is given by

[R H] (r) \equiv R H (R^{- 1} r),

(A4)

where

R

is the matrix representation of

R

. For all

R \in G

, we can prove the following equation [12

K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]

L_{m}^{(i j)} = \frac{1}{V} \int_{V} d r [R H^{(i) *}] (r) \cdot L [R H^{(j)}] (r - R r_{m}),

(A5)

where

L_{m}^{(i j)}

is defined by Eq. (1). Because

H^{(i)}

is an eigenfunction of

L_{s}

, it is a basis of an irreducible representation of

G

, so

R H^{(i)}

can be expressed by a linear combination of eigenfunctions of the same irreducible representation [14

T. Inui, Y. Tanabe, and Y. Onodera, Group Theory and Its Applications in Physics (Springer, 1990).

]. Then, the right side of Eq. (A5) must be a linear combination of

{L_{m}^{(i j)}}

, so we obtain their mutual relations. We can thus derive many relations among

{L_{m}^{(i j)}}

for

C_{6 v}

and

C_{4 v}

, as will be shown in Appendices B and C.

Now, we solve the eigenvalue problem given by

L H_{k} (r) = \frac{ω_{k}^{2}}{c^{2}} H_{k} (r),

(A6)

by the tight-binding approximation. In Eq. (A6),

k

is a wave vector in the Brillouin zone. The secular equation has the following form [12

K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]

| B - \frac{ω_{k}^{2}}{c^{2}} I | = 0,

(A7)

where

I

is a unit tensor and elements of matrix

B

are given by

B_{i j} = \sum_{m} e^{i k \cdot r_{m}} L_{m}^{(i j)} .

(A8)

When we analyze the dispersion curves for small

k

around the Brillouin-zone center, the summation for

m

can be limited to the nearest-neighbor lattice points, since other terms give contributions of higher orders of

k

APPENDIX B: TRIANGULAR LATTICE

1. $A_{1} + B_{1}$

In this appendix, we examine the dispersion relation in the vicinity of the Brillouin-zone center for the triangular lattice of unit cells of the

C_{6 v}

symmetry. We denote the origin and the nearest-neighbor lattice points by 0 to 6, as shown in Fig. 6, where symmetry operations of the

C_{6 v}

point group are also shown.

Fig. 6. Numbering of the lattice points of the triangular lattice and symmetry operations of the

C_{6 v}

point group.

We start with the case of accidental degeneracy of an

A_{1}

mode and a

B_{1}

mode on the

Γ

point of the Brillouin zone. We denote the magnetic fields of the

A_{1}

and

B_{1}

modes by

H^{(1)}

and

H^{(2)}

, respectively. By using all

R \in C_{6 v}

in Eq. (A5), we can derive the following relations:

L_{0}^{(11)} \equiv M_{10}, L_{0}^{(22)} \equiv M_{20},

(B1)

L_{0}^{(12)} = L_{0}^{(21)} = 0,

(B2)

L_{1}^{(11)} = L_{2}^{(11)} = L_{3}^{(11)} = L_{4}^{(11)} = L_{5}^{(11)} = L_{6}^{(11)}, \equiv M_{11},

(B3)

L_{1}^{(22)} = L_{2}^{(22)} = L_{3}^{(22)} = L_{4}^{(22)} = L_{5}^{(22)} = L_{6}^{(22)}, \equiv M_{21},

(B4)

L_{1}^{(12)} = L_{3}^{(12)} = L_{5}^{(12)} = - L_{2}^{(12)} = - L_{4}^{(12)} = - L_{6}^{(12)} \equiv M_{1} .

(B5)

In addition to these relations, we can generally prove the following equation.

L_{m}^{(i j)} = L_{- m}^{(j i) *},

(B6)

where

- m

denotes the lattice point given by

- r_{m}

From Eq. (A8), we obtain

B_{11} = ξ_{Γ}^{(1)} + M_{11} [2 (\cos k_{x} a - 1) + 4 (\cos \frac{k_{x} a}{2} \cos \frac{\sqrt{3} k_{y} a}{2} - 1)],

(B7)

B_{22} = ξ_{Γ}^{(2)} + M_{21} [2 (\cos k_{x} a - 1) + 4 (\cos \frac{k_{x} a}{2} \cos \frac{\sqrt{3} k_{y} a}{2} - 1)],

(B8)

B_{12} = 2 i M_{1} (\sin k_{x} a - 2 \sin \frac{k_{x} a}{2} \cos \frac{\sqrt{3} k_{y} a}{2}),

(B9)

where

ξ_{Γ}^{(1)} \equiv M_{10} + 6 M_{11},

(B10)

ξ_{Γ}^{(2)} \equiv M_{20} + 6 M_{21},

(B11)

which are the square of the eigen angular frequencies at the

Γ

point divided by

c^{2}

for the

A_{1}

and

B_{1}

modes, respectively. Because the secular equation in the present case is a quadratic equation of

ω^{2}

, it can easily be solved. When we assume the accidental degeneracy of the two modes at the

Γ

point—that is,

ξ_{Γ}^{(1)} = ξ_{Γ}^{(2)} \equiv ξ_{Γ}

—the dispersion relation for small

k

ω \approx {\begin{cases} ω_{Γ} - \frac{3 a^{2} c^{2} k^{2} M_{11}}{2 ω_{Γ}}, \\ ω_{Γ} - \frac{3 a^{2} c^{2} k^{2} M_{21}}{2 ω_{Γ}}, \end{cases}

(B12)

where

k = \sqrt{k_{x}^{2} + k_{y}^{2}}

and

ω_{Γ} = c \sqrt{ξ_{Γ}}

M_{1}

gives higher-order corrections. So, the linear term is absent and the dispersion relations are quadratic in

k

in the vicinity of the

Γ

point.

2. $A_{1} + B_{2}$

For the combination of an

A_{1}

mode (mode 1) and a

B_{2}

mode (mode2), we can derive

L_{m}^{(12)} = 0 (m = 0, 1, \dots, 6)

(B13)

by examining Eq. (A5) for all

R \in C_{6 v}

. So, there is no interaction between the two modes and the dispersion relations are quadratic for small

k

3. $A_{1} + E_{1}$

E_{1}

modes are doubly degenerate. From group theory [14

T. Inui, Y. Tanabe, and Y. Onodera, Group Theory and Its Applications in Physics (Springer, 1990).

], we can assume without loss of generality that two eigenfunctions of the

E_{1}

mode transform like

x

and

y

. We denote the magnetic field of the

A_{1}

mode by

H^{(0)}

for this time and those of the

E_{1}

mode by

H^{(1)}

, which transforms like

x

, and

H^{(2)}

, which transforms like

y

. Then, we can derive the following relations from Eq. (A5):

L_{0}^{(00)} \equiv M_{01}, L_{0}^{(11)} \equiv L_{0}^{(22)} = M_{11},

(B14)

L_{1}^{(00)} = L_{2}^{(00)} = L_{3}^{(00)} = L_{4}^{(00)} = L_{5}^{(00)} = L_{6}^{(00)} \equiv M_{02},

(B15)

L_{1}^{(11)} = L_{4}^{(11)} \equiv M_{12}, L_{1}^{(22)} = L_{4}^{(22)} \equiv M_{13},

(B16)

L_{2}^{(11)} = L_{3}^{(11)} = L_{5}^{(11)} = L_{6}^{(11)} \equiv M_{14},

(B17)

L_{2}^{(22)} = L_{3}^{(22)} = L_{5}^{(22)} = L_{6}^{(22)} \equiv M_{15},

(B18)

L_{0}^{(12)} = L_{1}^{(12)} = L_{4}^{(12)} = 0,

(B19)

L_{2}^{(12)} = - L_{3}^{(12)} = L_{5}^{(12)} = - L_{6}^{(12)} \equiv M_{16},

(B20)

L_{0}^{(01)} = L_{0}^{(02)} = L_{1}^{(02)} = L_{4}^{(02)} = 0,

(B21)

L_{1}^{(01)} = - L_{4}^{(01)} \equiv M_{1},

(B22)

L_{2}^{(01)} = - L_{3}^{(01)} = - L_{5}^{(01)} = L_{6}^{(01)} \equiv M_{2},

(B23)

L_{2}^{(02)} = L_{3}^{(02)} = - L_{5}^{(02)} = - L_{6}^{(02)} \equiv M_{3} .

(B24)

We can further prove that

M_{14} = \frac{M_{12} + 3 M_{13}}{4},

(B25)

M_{15} = \frac{3 M_{12} + M_{13}}{4},

(B26)

M_{16} = \frac{\sqrt{3} (M_{12} - M_{13})}{4},

(B27)

M_{1} = 2 M_{2}, M_{3} = \sqrt{3} M_{2} .

(B28)

Then, we can obtain the elements of matrix

B

according to Eq. (A8). The secular equation is a cubic equation of

ω^{2}

. To obtain its analytical solutions, we need a lengthy calculation similar to that presented in [11

K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef]

] The solutions are, however, rather simple when we assume the accidental degeneracy of the

A_{1}

and

E_{1}

modes at the

Γ

point. The result is

ω \approx {\begin{cases} ω_{Γ} \pm \frac{3 a k c^{2} | M_{2} |}{ω_{Γ}} - \frac{a^{2} k^{2} M}{4 ω_{Γ}}, \\ ω_{Γ} - \frac{a^{2} k^{2} M}{4 ω_{Γ}}, \end{cases}

(B29)

where

M = M_{02} + M_{12} + M_{13}

. Two of the three dispersion curves are linear and isotropic with respect to

k

in the vicinity of the

Γ

point, and the third one is quadratic. So, we have a Dirac cone with an auxiliary quadratic dispersion surface.

4. $B_{2} + E_{2}$

E_{2}

modes are doubly degenerate. From group theory [14

T. Inui, Y. Tanabe, and Y. Onodera, Group Theory and Its Applications in Physics (Springer, 1990).

], we can assume without loss of generality that two eigenfunctions of the

E_{2}

mode transform like

2 x y

and

x^{2} - y^{2}

. We denote the magnetic field of the

B_{2}

mode by

H^{(0)}

and those of the

E_{2}

mode by

H^{(1)}

, which transforms like

2 x y

, and

H^{(2)}

, which transforms like

x^{2} - y^{2}

. Then, we can derive the relations among

L_{m}^{(i j)}

’s by using Eq. (A5). The results are completely the same as the case of

A_{1} + E_{1}

. So we have a Dirac cone with an auxiliary quadratic dispersion surface in the vicinity of the

Γ

point.

5. $A_{2} + E_{1}$

We denote the magnetic field of the

A_{2}

mode by

H^{(0)}

and those of the

E_{1}

mode by

H^{(1)}

, which transforms like

y

, and

H^{(2)}

, which transforms like

x

, for this time. Then, we can derive the relations among

L_{m}^{(i j)}

’s according to Eq. (A5). We obtain the same results as the case of

A_{1} + E_{1}

, except

M_{3} = - \sqrt{3} M_{2} .

(B30)

Since the dispersion relation, Eq. (B29), does not depend on the sign of

M_{3}

, we obtain the same dispersion, i.e., a Dirac cone with an auxiliary quadratic dispersion surface.

6. $B_{2} + E_{1}$

We denote the magnetic field of the

B_{2}

mode by

H^{(0)}

and those of the

E_{1}

mode by

H^{(1)}

, which transforms like

x

, and

H^{(2)}

, which transforms like

y

. Then, the relations among

L_{m}^{(i j)}

’s with

1 \leq i

j \leq 2

are the same as in Appendix B.3. For the rest of

L_{m}^{(i j)}

’s, we can derive the following relations:

L_{0}^{(00)} \equiv M_{01},

(B31)

L_{1}^{(00)} = L_{2}^{(00)} = L_{3}^{(00)} = L_{4}^{(00)} = L_{5}^{(00)} = L_{6}^{(00)} \equiv M_{02},

(B32)

L_{0}^{(01)} = L_{0}^{(02)} = L_{1}^{(01)} = L_{4}^{(01)} = 0,

(B33)

L_{2}^{(01)} = - L_{3}^{(01)} = L_{5}^{(01)} = - L_{6}^{(01)} \equiv M_{1},

(B34)

L_{1}^{(02)} = L_{4}^{(02)} \equiv M_{2},

(B35)

L_{2}^{(02)} = L_{3}^{(02)} = L_{5}^{(02)} = L_{6}^{(02)} \equiv M_{3} .

(B36)

We can further prove

M_{1} = - \sqrt{3} M_{3}, M_{2} = - 2 M_{3} .

(B37)

Then, we obtain the elements of matrix

B

in Eq. (A8) and solve the secular equation, Eq. (A7). After a simple but lengthy calculation, we obtain the dispersion relation as follows.

ω \approx {\begin{cases} ω_{Γ} - \frac{[M \pm \sqrt{3 ({| M_{1} |}^{2} + {| M_{16} |}^{2} \sin^{2} 2 ϕ)}] a^{2} c^{2} k^{2}}{4 ω_{Γ}}, \\ ω_{Γ} - \frac{a^{2} c^{2} k^{2} M}{4 ω_{Γ}}, \end{cases}

(B38)

where

M = M_{02} + M_{12} + M_{13}

and

ϕ

is the angle between vector

k

and the

x

axis. So, there are three quadratic dispersion surfaces.

7. $A_{2} + E_{2}$

We denote the magnetic field of the

A_{2}

mode by

H^{(0)}

and those of the

E_{2}

mode by

H^{(1)}

, which transforms like

x^{2} - y^{2}

, and

H^{(2)}

, which transforms like

2 x y

, for this time. Then, we can derive the relations among

L_{m}^{(i j)}

’s according to Eq. (A5). We obtain the same results as the case of

B_{2} + E_{1}

, except

M_{1} = \sqrt{3} M_{3} .

(B39)

Because the dispersion relation, Eq. (B38), does not depend on the sign of

M_{1}

, we obtain the same dispersion, i.e., three quadratic dispersion surfaces.

8. $A_{1} + E_{2}$

We denote the magnetic field of the

A_{1}

mode by

H^{(0)}

and those of the

E_{2}

mode by

H^{(1)}

, which transforms like

2 x y

, and

H^{(2)}

, which transforms like

x^{2} - y^{2}

. Then, the relations derived for

L_{m}^{(i j)}

’s are the same as the case of

B_{2} + E_{1}

. So, there are three quadratic dispersion surfaces in the vicinity of the

Γ

point.

APPENDIX C: SQUARE LATTICE

1. $A_{1} + B_{1}$

Next, we examine the square lattice of unit cells of the

C_{4 v}

symmetry. We denote the origin and the nearest-neighbor lattice points by 0 to 4, as shown in Fig. 7, where symmetry operations of the

C_{4 v}

point group are also shown.

Fig. 7. Numbering of the lattice points of the square lattice and symmetry operations of the

C_{4 v}

point group.

We start with the case of accidental degeneracy of an

A_{1}

mode and a

B_{1}

mode in the Brillouin-zone center. We denote the magnetic fields of the

A_{1}

and

B_{1}

modes by

H^{(1)}

and

H^{(2)}

, respectively. By using all

R \in C_{4 v}

in Eq. (A5), we can derive the following relations:

L_{0}^{(11)} \equiv M_{10}, L_{0}^{(22)} \equiv M_{20}, L_{0}^{(12)} = 0,

(C1)

L_{1}^{(11)} = L_{2}^{(11)} = L_{3}^{(11)} = L_{4}^{(11)} \equiv M_{11},

(C2)

L_{1}^{(22)} = L_{2}^{(22)} = L_{3}^{(22)} = L_{4}^{(22)} \equiv M_{22},

(C3)

L_{1}^{(12)} = - L_{2}^{(12)} = L_{3}^{(12)} = - L_{4}^{(12)} \equiv M_{1} .

(C4)

Then, we obtain the elements of matrix

B

and solve the secular equation. The dispersion relation thus obtained for small

k

in the vicinity of the

Γ

point is

ω \approx ω_{Γ} - \frac{a^{2} c^{2} k^{2}}{4 ω_{Γ}} [M_{11} + M_{21} \pm \sqrt{{(M_{11} - M_{21})}^{2} + 4 {| M_{1} |}^{2} \cos ϕ}],

(C5)

where

ϕ

is the angle between vector

k

and the

x

axis. So, there are two quadratic dispersion surfaces.

2. $A_{1} + A_{2}$

For this combination,

L_{0}^{(12)} = L_{1}^{(12)} = L_{2}^{(12)} = L_{3}^{(12)} = L_{4}^{(12)} = 0 .

(C6)

So, there is no interaction between the two modes and the dispersion relations are quadratic.

3. $A_{1} + B_{2}$

Equation (C6) holds for this combination, too. So, the dispersion relations are quadratic with respect to

k

in the vicinity of the

Γ

point.

4. $A_{2} + B_{2}$

We have completely the same relations for

L_{m}^{(i j)}

’s as

A_{1} + B_{1}

. So, there are two quadratic dispersion surfaces.

5. $B_{1} + E$

E

modes are doubly degenerate. From group theory [14

T. Inui, Y. Tanabe, and Y. Onodera, Group Theory and Its Applications in Physics (Springer, 1990).

], we can assume without loss of generality that two eigenfunctions of the

E

mode transform like

x

and

y

. We denote the magnetic field of the

B_{1}

mode by

H^{(0)}

and those of the

E

mode by

H^{(1)}

, which transforms like

x

, and

H^{(2)}

, which transforms like

y

. Then, we can derive the following relations by examining Eq. (A5) for all

R \in C_{4 v}

L_{0}^{(00)} \equiv M_{01}, L_{0}^{(11)} \equiv L_{0}^{(22)} = M_{11},

(C7)

L_{1}^{(00)} = L_{2}^{(00)} = L_{3}^{(00)} = L_{4}^{(00)} \equiv M_{02},

(C8)

L_{1}^{(11)} = L_{3}^{(11)} = L_{2}^{(22)} = L_{4}^{(22)} \equiv M_{12},

(C9)

L_{2}^{(11)} = L_{4}^{(11)} = L_{1}^{(22)} = L_{3}^{(22)} \equiv M_{13},

(C10)

L_{0}^{(12)} = L_{1}^{(12)} = L_{2}^{(12)} = L_{3}^{(12)} = L_{4}^{(12)} = 0,

(C11)

L_{0}^{(01)} = L_{0}^{(02)} = 0,

(C12)

L_{1}^{(01)} = - L_{3}^{(01)} = - L_{2}^{(02)} = L_{4}^{(02)} \equiv M_{1},

(C13)

L_{1}^{(02)} = L_{3}^{(02)} = L_{2}^{(01)} = L_{4}^{(01)} = 0 .

(C14)

Then, from Eq. (A8), elements of matrix

B

are

B_{00} = M_{01} + 2 M_{02} (\cos k_{x} a + \cos k_{y} a),

(C15)

B_{11} = M_{11} + 2 M_{12} \cos k_{x} a + 2 M_{13} \cos k_{y} a,

(C16)

B_{22} = M_{11} + 2 M_{13} \cos k_{x} a + 2 M_{12} \cos k_{y} a,

(C17)

B_{01} = 2 i M_{1} \sin k_{x} a,

(C18)

B_{02} = - 2 i M_{1} \sin k_{y} a,

(C19)

B_{12} = 0 .

(C20)

When we compare the structure of matrix

B

with that for the combination of

A_{1}

and

E

modes, Eqs. (9)–(14) of [11

K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef]

], they are the same with each other except the sign of

B_{02}

, which is equivalent to changing the sign of

k_{y}

. So, the dispersion relation of the

B_{1} + E

case is obtained by changing the sign of

k_{y}

in the dispersion relation of

A_{1} + E

. Because the structure of the metamaterial that we assume in this paper is symmetric about the

x

axis due to the

C_{4 v}

symmetry, the dispersion relation is invariant for the change of the sign of

k_{y}

. In conclusion, the shape of the dispersion curves is the same for the

B_{1} + E

case and the

A_{1} + E

case, which is a Dirac cone with an auxiliary quadratice dispersion surface.

6. $B_{2} + E$

Finally, we examine the combination of a

B_{2}

mode and an

E

mode. We denote the magnetic field of the

B_{2}

mode by

H^{(0)}

and those of the

E

mode by

H^{(1)}

, which transforms like

y

, and

H^{(2)}

, which transforms like

x

, for this time. Then, the relations among

L_{m}^{(i j)}

are the same as the

B_{1} + E

case except

L_{1}^{(01)} = - L_{3}^{(01)} = L_{2}^{(02)} = - L_{4}^{(02)} \equiv M_{1} .

(C21)

The structure of these relations are completely the same as the

A_{1} + E

case, Eqs. (27)–(32) of [11

K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef]

]. So, the accidental degeneracy of the

B_{2}

and

E

modes at the

Γ

point results in a Dirac cone with an auxiliary quadratice dispersion surface.

ACKNOWLEDGMENTS

This study was supported by a Grant-in-Aid for Scientific Research on Innovative Areas from the Japanese Ministry of Education, Culture, Sports and Technology (Grant No. 22109007).

REFERENCES

1.	F. D. M. Haldane and S. Raghu, “Possible realization of directional optical waveguides in photonic crystals with broken time-reversal symmetry,” Phys. Rev. Lett. 100, 013904 (2008). [CrossRef]
2.	S. Raghu and F. D. M. Haldane, “Analogs of quantum-Hall-effect edge states in photonic crystals,” Phys. Rev. A 78, 033834 (2008). [CrossRef]
3.	T. Ochiai and M. Onoda, “Photonic analog of graphene model and its extension: Dirac cone, symmetry, and edge states,” Phys. Rev. B 80, 155103 (2009). [CrossRef]
4.	X. Zhang, “Observing zitterbewegung for photons near the Dirac point of a two-dimensional photonic crystal,” Phys. Rev. Lett. 100, 113903 (2008). [CrossRef]
5.	R. A. Sepkhanov, Y. B. Bazaliy, and C. W. J. Beenakker, “Extremal transmission at the Dirac point of a photonic band structure,” Phys. Rev. A 75, 063813 (2007). [CrossRef]
6.	M. Diem, T. Koschny, and C. M. Soukoulis, “Transmission in the vicinity of the Dirac point in hexagonal photonic crystals,” Physica B 405, 2990–2995 (2010). [CrossRef]
7.	S. H. Nam, A. J. Taylor, and A. Efimov, “Diabolical point and conical-like diffraction in periodic plasmonic nanostructures,” Opt. Express 18, 10120–10126 (2010). [CrossRef]
8.	X. Huang, Y. Lai, Z. H. Hang, H. Zheng, and C. T. Chan, “Dirac cones induced by accidental degeneracy in photonic crystals and zero-refractive-index materials,” Nat. Mater. 10, 582–586 (2011). [CrossRef]
9.	M. Silveirinha and N. Engheta, “Tunneling of electromagnetic energy through subwavelength channels and bends using $ε$ -near-zero materials,” Phys. Rev. Lett. 97, 157403 (2006). [CrossRef]
10.	A. Alu, M. G. Silveirinha, A. Salandrino, and N. Engheta, “Epsilon-near-zero metamaterials and electromagnetic sources: Tailoring the radiation phase pattern,” Phys. Rev. B 75, 155410 (2007). [CrossRef]
11.	K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef]
12.	K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]
13.	K. Sakoda, Optical Properties of Photonic Crystals , 2nd ed. (Springer-Verlag, 2004).
14.	T. Inui, Y. Tanabe, and Y. Onodera, Group Theory and Its Applications in Physics (Springer, 1990).

OCIS Codes
(350.3618) Other areas of optics : Left-handed materials
(160.3918) Materials : Metamaterials
(160.5298) Materials : Photonic crystals

ToC Category:
Materials

History
Original Manuscript: July 18, 2012
Manuscript Accepted: August 8, 2012
Published: September 18, 2012

Citation
Kazuaki Sakoda, "Universality of mode symmetries in creating photonic Dirac cones," J. Opt. Soc. Am. B 29, 2770-2778 (2012)
http://www.opticsinfobase.org/josab/abstract.cfm?URI=josab-29-10-2770

Sort: Author | Year | Journal | Reset

References

F. D. M. Haldane and S. Raghu, “Possible realization of directional optical waveguides in photonic crystals with broken time-reversal symmetry,” Phys. Rev. Lett. 100, 013904 (2008). [CrossRef]
S. Raghu and F. D. M. Haldane, “Analogs of quantum-Hall-effect edge states in photonic crystals,” Phys. Rev. A 78, 033834 (2008). [CrossRef]
T. Ochiai and M. Onoda, “Photonic analog of graphene model and its extension: Dirac cone, symmetry, and edge states,” Phys. Rev. B 80, 155103 (2009). [CrossRef]
X. Zhang, “Observing zitterbewegung for photons near the Dirac point of a two-dimensional photonic crystal,” Phys. Rev. Lett. 100, 113903 (2008). [CrossRef]
R. A. Sepkhanov, Y. B. Bazaliy, and C. W. J. Beenakker, “Extremal transmission at the Dirac point of a photonic band structure,” Phys. Rev. A 75, 063813 (2007). [CrossRef]
M. Diem, T. Koschny, and C. M. Soukoulis, “Transmission in the vicinity of the Dirac point in hexagonal photonic crystals,” Physica B 405, 2990–2995 (2010). [CrossRef]
S. H. Nam, A. J. Taylor, and A. Efimov, “Diabolical point and conical-like diffraction in periodic plasmonic nanostructures,” Opt. Express 18, 10120–10126 (2010). [CrossRef]
X. Huang, Y. Lai, Z. H. Hang, H. Zheng, and C. T. Chan, “Dirac cones induced by accidental degeneracy in photonic crystals and zero-refractive-index materials,” Nat. Mater. 10, 582–586 (2011). [CrossRef]
M. Silveirinha and N. Engheta, “Tunneling of electromagnetic energy through subwavelength channels and bends using ε-near-zero materials,” Phys. Rev. Lett. 97, 157403 (2006). [CrossRef]
A. Alu, M. G. Silveirinha, A. Salandrino, and N. Engheta, “Epsilon-near-zero metamaterials and electromagnetic sources: Tailoring the radiation phase pattern,” Phys. Rev. B 75, 155410 (2007). [CrossRef]
K. Sakoda, “Dirac cone in two- and three-dimensional metamaterials,” Opt. Express 20, 3898–3917 (2012). [CrossRef]
K. Sakoda, “Double Dirac cones in triangular-lattice metamaterials,” Opt. Express 20, 9925–9939 (2012). [CrossRef]
K. Sakoda, Optical Properties of Photonic Crystals, 2nd ed. (Springer-Verlag, 2004).
T. Inui, Y. Tanabe, and Y. Onodera, Group Theory and Its Applications in Physics (Springer, 1990).

Alu, A.
Bazaliy, Y. B.
Beenakker, C. W. J.
Chan, C. T.
Diem, M.
Efimov, A.
Engheta, N.
Haldane, F. D. M.
Hang, Z. H.
Huang, X.
Inui, T.
Koschny, T.
Lai, Y.
Nam, S. H.
Ochiai, T.
Onoda, M.
Onodera, Y.
Raghu, S.
Sakoda, K.
Salandrino, A.
Sepkhanov, R. A.
Silveirinha, M.
Silveirinha, M. G.
Soukoulis, C. M.
Tanabe, Y.
Taylor, A. J.
Zhang, X.
Zheng, H.

Nat. Mater.

X. Huang, Y. Lai, Z. H. Hang, H. Zheng, and C. T. Chan, “Dirac cones induced by accidental degeneracy in photonic crystals and zero-refractive-index materials,” Nat. Mater. 10, 582–586 (2011). [CrossRef]

Opt. Express

Phys. Rev. A

S. Raghu and F. D. M. Haldane, “Analogs of quantum-Hall-effect edge states in photonic crystals,” Phys. Rev. A 78, 033834 (2008). [CrossRef]
R. A. Sepkhanov, Y. B. Bazaliy, and C. W. J. Beenakker, “Extremal transmission at the Dirac point of a photonic band structure,” Phys. Rev. A 75, 063813 (2007). [CrossRef]

Phys. Rev. B

T. Ochiai and M. Onoda, “Photonic analog of graphene model and its extension: Dirac cone, symmetry, and edge states,” Phys. Rev. B 80, 155103 (2009). [CrossRef]
A. Alu, M. G. Silveirinha, A. Salandrino, and N. Engheta, “Epsilon-near-zero metamaterials and electromagnetic sources: Tailoring the radiation phase pattern,” Phys. Rev. B 75, 155410 (2007). [CrossRef]

Phys. Rev. Lett.

X. Zhang, “Observing zitterbewegung for photons near the Dirac point of a two-dimensional photonic crystal,” Phys. Rev. Lett. 100, 113903 (2008). [CrossRef]
F. D. M. Haldane and S. Raghu, “Possible realization of directional optical waveguides in photonic crystals with broken time-reversal symmetry,” Phys. Rev. Lett. 100, 013904 (2008). [CrossRef]
M. Silveirinha and N. Engheta, “Tunneling of electromagnetic energy through subwavelength channels and bends using ε-near-zero materials,” Phys. Rev. Lett. 97, 157403 (2006). [CrossRef]

Physica B

M. Diem, T. Koschny, and C. M. Soukoulis, “Transmission in the vicinity of the Dirac point in hexagonal photonic crystals,” Physica B 405, 2990–2995 (2010). [CrossRef]

Other

K. Sakoda, Optical Properties of Photonic Crystals, 2nd ed. (Springer-Verlag, 2004).
T. Inui, Y. Tanabe, and Y. Onodera, Group Theory and Its Applications in Physics (Springer, 1990).

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

Click here to see a list of articles that cite this paper

Figures


Fig. 1.	Fig. 2.	Fig. 3.


Fig. 4.	Fig. 5.	Fig. 6.


Fig. 7.

« Previous Article | Next Article »

OSA is a member of CrossRef.

Select a Conference
Session or Paper #
Year

OpticsInfoBase

Journal of the Optical Society of America B

| OPTICAL PHYSICS

Universality of mode symmetries in creating photonic Dirac cones

Article Tools

Article Navigation

Abstract

1. INTRODUCTION

2. THEORY

3. RESULTS AND DISCUSSION

4. CONCLUSION

Appendices

APPENDIX A: FORMULATION

APPENDIX B: TRIANGULAR LATTICE

1. A1+ B1

2. A1+ B2

3. A1+ E1

4. B2+ E2

5. A2+ E1

6. B2+ E1

7. A2+ E2

8. A1+ E2

APPENDIX C: SQUARE LATTICE

1. A1+ B1

2. A1+ A2

3. A1+ B2

4. A2+ B2

5. B1+E

6. B2+E

ACKNOWLEDGMENTS

REFERENCES

References

Nat. Mater.

Opt. Express

Phys. Rev. A

Phys. Rev. B

Phys. Rev. Lett.

Physica B

Other

Cited By

Figures

1. $A_{1} + B_{1}$

2. $A_{1} + B_{2}$

3. $A_{1} + E_{1}$

4. $B_{2} + E_{2}$

5. $A_{2} + E_{1}$

6. $B_{2} + E_{1}$

7. $A_{2} + E_{2}$

8. $A_{1} + E_{2}$

1. $A_{1} + B_{1}$

2. $A_{1} + A_{2}$

3. $A_{1} + B_{2}$

4. $A_{2} + B_{2}$

5. $B_{1} + E$

6. $B_{2} + E$