## Topological nature of bound states in the radiation continuum

Bo Zhen, 1, ∗ Chia Wei Hsu, 1, 2, † Ling Lu, 1 A. Douglas Stone, 3 and Marin Soljaˇ ci´ c 1

1 Department of Physics, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, Massachusetts 02139, USA 2 Department of Physics, Harvard University, 17 Oxford Street, Cambridge, Massachusetts 02138, USA 3 Department of Applied Physics, Post Office Box 208284, Yale University, New Haven, CT 06520, USA (Dated: August 4, 2014)

∗ E-mail: bozhen@mit.edu

† These authors contributed equally to this work.

Bound states in the continuum (BICs) are unusual solutions of wave equations describing light or matter: they are discrete and spatially bounded, but exist at the same energy as a continuum of states which propagate to infinity. Until recently, BICs were constructed through fine-tuning parameters in the wave equation 1-5 or exploiting the separability of the wave equation due to symmetry 6-8 . More recently, BICs that that are both robust and not symmetryprotected ('accidental') have been predicted 9-13 and experimentally realized 13 in periodic structures; the simplest such system is a periodic dielectric slab 13 , which also has symmetry-protected BICs. Here we show that both types of BICs in such systems are vortex centers in the polarization direction of far-field radiation. The robustness of these BICs is due to the existence of conserved and quantized topological charges, defined by the number of times the polarization vectors wind around the vortex centers. Such charges can only be generated or annihilated by making large changes in the system parameters, and then only according to strict rules, which we derive and test numerically. Our results imply that laser emission based on such states will generate vector beams 14 .

These BICs exist in photonics crystal slabs, which are dielectric slabs with a periodic modulation of refractive index at the wavelength scale 15 surrounded by air or a low-index medium. Photonic crystal slabs are used in many applications including surface-emitting lasers 16 , photovoltaics 17 , LEDs 18 , and biosensing 19,20 . The periodic modulation alters the dispersion relation of light in the slabs and gives rise to photonic bands, analogous to electronic band structures in solids. In general light can escape from the surface of the slab and propagate to the far-field, but a portion of the photonic band states (those 'below the light line') are perfectly confined to the slab by the generalized form of total internal reflection 15 . In contrast, modes above the light line appear generically as resonances with finite lifetimes due to their coupling to the continuum of extended modes 21 . It has been known for some time that bound states with infinite lifetimes also exist above the light line at isolated highsymmetry wavevectors. This type of BIC arises from the symmetry mismatch between their mode profiles inside the photonic crystal and those of the external propagating modes (radiation continuum) 8,21-23 . Recently, a new type of BIC, unprotected by symmetry, has been found (both theoretically and experimentally) to exist at arbitrary wavevectors in bands above the light line in photonic crystal slabs 13 . These BICs occur 'accidentally', when the

relevant couplings to the continuum (see below) all vanish simultaneously. However their existence does not require fine-tuning of system parameters; small changes in parameters simply shift the position of these special points along the band diagram. An intuitive understanding of why such BICs exist and are robust was previously lacking. Recently, an explanation based on accidental triangular symmetry of the radiating fields was proposed 24 but does not explain the robustness of these BICs and their occurrences in TE-like bands.

We now show that both types of BICs in photonic crystal slabs are vortex centers in the polarization direction of the far-field radiation of the slabs. Using the Bloch theorem for photonic crystals 15 , we write the electric field of a resonance as E k ( ρ , z ) = e i k · ρ u k ( ρ , z ), where k = k x x ˆ + k y y ˆ is the two-dimensional wave vector, ρ = xx ˆ + yy ˆ is the in-plane coordinate, u k is a periodic function in ρ , and z is the normal direction to the slab. While the fields inside the slab are periodically modulated, outside the slab each state consists of propagating plane waves and/or evanescent waves that decay exponentially away from the surface. For states above the light line (resonances), and wavelengths below the diffraction limit, the only non-zero propagating-wave amplitudes are the zero-order (constant in-plane) Fourier coefficients of u k , given by c k ( ) = c x ( k )ˆ + x c y ( k )ˆ (Fig. 1a). y Here, c x ( k ) = ˆ x · 〈 u k 〉 , c y ( k ) = ˆ y · 〈 u k 〉 , and the brackets denote spatial average over one unit cell on any horizontal plane outside the slab. Note that c k ( ) is the projection of 〈 u k 〉 onto the xy plane; it points in the polarization direction of the resonance in the far field, so we refer to c k ( ) as the 'polarization vector'.

A resonance turns into a BIC when the outgoing power is zero, which happens if and only if c x = c y = 0. In general, c x and c y are both complex functions of k , and varying the wave vector components ( k , k x y ) is not sufficient to guarantee a solution where c x = c y = 0. However, when the system is invariant under the operation C T z 2 , implying that glyph[epsilon1] ( x, y, z ) = glyph[epsilon1] ∗ ( -x, -y, z ), we show that c x and c y can be chosen to be real numbers simultaneously; in other words, the far field is linearly polarized (see Supplementary Information, here C z 2 is 180 ◦ rotation operator around z axis, and T is the time reversal operator). When the system also has up-down mirror symmetry ( σ z ), the outgoing waves on one side of the slab determine those on the other; for such systems, BICs are stable because they correspond to the intersections between the nodal line of c x and the nodal line of c y in the k x -k y plane. Such a nodal intersection naturally causes a vortex in the polarization vector field centered on the BIC, as illustrated in Fig. 1b, for the simplest case. Along the nodal line of c x (or

c y ), the direction of c k ( ) is along the y axis (or x axis), as illustrated in Fig. 1b. As one encircles the nodal intersection (BIC) in the k x -k y plane each component of the polarization vector flips sign as its nodal line is crossed so as to create a net circulation of ± 2 π in the polarization field. At the nodal intersection the polarization direction becomes undefined, since at the BIC there is zero emission into the far-field. Conversely one could say that BICs cannot radiate because there is no way to assign a far-field polarization that is consistent with neighbouring k points. Thus robust BICs are only possible when there is vorticity in the polarization field.

Vortices are characterized by their topological charges. Here, the topological charge ( q ) carried by a BIC is defined as:

$$q = \frac { 1 } { 2 \pi } \oint _ { C } d { \mathbf k } \cdot \nabla _ { \mathbf k } \phi ( { \mathbf k } ), \quad q \in \mathbb { Z }$$

which describes how many times the polarization vector winds around the BIC. Here, φ ( k ) = arg[ c x ( k ) + ic y ( k )] is the angle of the polarization vector, and C is a closed simple path in k space that goes around the BIC in the counterclockwise direction.The fields u k are chosen to be smooth functions of k , so φ ( k ) is differentiable in k along the path. The polarization vector has to come back to itself after the closed loop, so the overall angle change must be an integer multiple of 2 π , and q must be an integer. Fig. 1c shows examples of how the polarization vector winds around a BIC with charge q = +1 and also around a BIC with charge q = -1 along a loop C marked by 1 → → → → 2 3 4 1. Similar definitions of winding numbers as in Eq. 1 can be found in describing topological defects 25 of continuous two-dimensional spins, dislocations in crystals, and quantized vortices in helium II 26 . This formalism describing polarization vortices is also closely related to Berry phases in describing adiabatic changes of polarization of light 27 and Dirac cones in gaphenes 28 .

The far-field pattern at a definite k point by itself does not reflect the vorticity of polarization around a BIC, but laser emission centered on such a BIC will. Laser emission always has a finite width in k -space and this wave-packet will be centered on the BIC; hence it will consist of a superposition of plane waves from the neighborhood of the BIC, leading naturally to a spatial twist in the polarization for the outgoing beam. Such beams have been studied previously, and are known as vector beams 14 , although their connection with BICs does not appear to have been realized. The number of twists in the polarization direction is known as the order number of the vector beam, and we now see that it is given by

the topological charge carried by the BIC. Note that these vector beams are different from optical vortices, which usually have a fixed polarization direction.

In the example of Fig. 2, we show the topological charges of BICs for a structure that has been experimentally realized in ref. 13. In this example, there are five BICs on the lowestfrequency TM-like band (Fig. 2a). We obtain polarization vectors c k ( ) from finite-difference time-domain (FDTD) calculations, which reveal five vortices with topological charges of ± 1 at these five k points (Fig. 2b). As discussed above, the BICs and their topological charges can also be identified from the nodal-line crossings and the gray-scale colors of c x and c y (Fig. 2c).

The winding number of polarization vector along a closed path is given by the sum of the topological charges carried by all BICs enclosed within this path 25 . When system parameters vary continuously, the winding number defined on this path remains invariant, unless there are BICs crossing the boundary. Therefore, topological charge is a conserved quantity. This conservation rule leads to consequences/restrictions on behaviors of the BICs (Supplementary Information). For example, as long as the system retains C T z 2 and σ z symmetries, a BIC can only be destroyed through annihilation with another BIC of the exact opposite charge, or through bringing it outside of the continuum (below the light line).

The conservation of topological charges allows us to predict and understand the behaviors of BICs when the parameters of the system are varied over a wide range, as we now illustrate. First, consider the lowest-frequency TM-like mode (TM 1 band) of a 1D-periodic structure in air shown in Fig. 3a. This grating consists of a periodic array of dielectric bars with periodicity of a , width w = 0 45 . a , and refractive index n = 1 45. . Its calculated band structure is shown in Fig. 3b. When the thickness of the grating is h = 1 50 . a , there are two BICs on the k x axis, as indicated by the radiative quality factor of the resonances (Fig. 3c). The polarization vector c k ( ), also shown in Fig. 3c, characterizes both BICs as carrying charges q = +1. When the grating thickness is decreased to h/a = 1 43 (all other . parameters fixed), the two BICs move towards the center of the Brillouin zone, meet at the Γ point, and deflect onto the k y axis (Fig. 3d). This is inevitable due to the conservation of the topological charges: annihilation cannot happen between two BICs of the same charge.

Annihilation of BICs is only possible when charges of opposite signs are present. This can be seen in the lowest-frequency TE-like band of the same structure (Fig. 3e,f). When

h/a = 1 04, there are two off-. Γ BICs with charge -1 and a BIC with charge +1 at the Γ point (Fig. 3e). As h/a decreases, the two -1 charges move to the center and eventually annihilate with the +1 charge, leaving only one BIC with charge q = -1 (Fig. 3f).

Generation of BICs is also restricted by charge conservation, and can be understood as the reverse process of annihilation. We provide an example by considering the lowest-frequency TE-like mode in a photonic crystal slab of n = 3 6 with a square lattice of cylindrical air holes . of diameter d = 0 5 . a (Fig. 4a). As the slab thickness increases, BICs are generated at the Γ point. Each time, four pairs of BICs with exact opposite charges are generated, consistent with charge conservation and C 4 v symmetry of the structure. With further increase of the slab thickness, the eight BICs move outward along high-symmetry lines and eventually go outside of the continuum (fall below the light line).

Although the examples discussed so far only show topological charges of ± 1, other values of charges can be found in higher-frequency bands of the PC or in structures with higher rotational symmetry. For example, Fig. S2 in Supplementary Information shows a stable BIC of charge -2 at the Γ point arising from the double degeneracy of nodal lines caused by the C 6 v symmetry of the system.

The symmetries of the system also restrict the possible values of topological charges, since the nodal curves must respect the point symmetry. BICs at k points related by in-plane point-group symmetries of the system (mirror reflections and rotations) have the same topological charges (Supplementary Information). Using this fact, one can calculate all possible topological charges of BICs at high symmetry k points. For a system with C n symmetry, the possible topological charges at the Γ point on a singly-degenerate band are given in Table S1 (Supplementary information). This is consistent with all examples in this paper. This table can be used to predict the charges in other systems of interest and to design high order vector beams.

## Conclusions

We have demonstrated that BICs in photonic crystal slabs are associated with vortices in the polarization field and explained their robustness in terms of conserved topological charges. We derive the symmetries that constrain these charges and explain their generation, evolution and annihilation. We conjecture that all robust BICs 6-13 will correspond to vortices in an appropriate parameter space. Our finding connects electromagnetic BICs to a wide range of physical phenomena including Berry phases around Dirac points 28 , topo-

logical defects 25 , and general vortex physics 26 . Optical BICs in photonic crystals have a wealth of applications. Lasing action can naturally occur at BIC states where the quality factor diverges. The angular (wavevector) tunablity of the BICs makes them great candidates for on-chip beam-steering 29 . Furthermore, photonic crystal lasers through BICs are naturally vector beams 30,31 , which are important for particle accelerations, optical trapping and stimulated emission depletion microscopy.

## Acknowledgments

The authors thank Chong Wang, Scott Skirlo, David Liu, Fan Wang, Nicholas Rivera, Dr. Ido Kaminer, Dr. Homer Reid, Prof. Xiaogang Wen, Prof. Steven G. Johnson, and Prof. John D. Joannopoulos for helpful discussions. This work was partly supported by the Army Research Office through the Institute for Soldier Nanotechnologies under contract no. W911NF-07-D0004. B.Z., L.L., and M.S. were partly supported by S3TEC, an Energy Frontier Research Center funded by the US Department of Energy under grant no. desc0001299. L.L. was supported in part by the Materials Research Science and Engineering Center of the National Science Foundation (award no. DMR-0819762). A.D.S. was partly supported by NSF grant DMR-1307632.

## Author contributions

All authors discussed the results and made critical contributions to the work.

## Additional information

Supplementary information is available in the online version of the paper. Reprints and permissions information is available online at www.nature.com/reprints. Correspondence and requests for materials should be addressed to B.Z.

## Competing financial interests

The authors declare no competing financial interests.

## Supplementary information

## A. Symmetry requirements for stable BICs

Here, we give the proof that stable BICs at arbitrary k points can be found when the system is invariant under C T z 2 and σ z operators, and that stable BICs at C z 2 -invariant k points can be found when the system has C z 2 symmetry. Here, C z 2 means 180 ◦ rotation around z axis, and T means the time reversal operator. The schematics of the symmetry requirement is summarized in Fig. S1.

In region I , systems are invariant under the symmetry operator C T z 2 , namely glyph[epsilon1] glyph[star] ( x, y, z ) = glyph[epsilon1] ( -x, -y, z ). Let u k be an eigenfunction of the master operator 15 Θ = k 1 glyph[epsilon1] ( ∇ + i k ) × ∇ ( + i k ) × , and recall that k here only has x and y components since we are considering a slab structure that does not have translational symmetry in z . A short derivation shows that at any k point, u k ( r ) and C z 2 u ∗ k ( C z 2 r ) are both eigenfunctions of Θ k ( r ) with the same eigenvalue, so they must differ at most by a phase factor,

$$u _ { k } ( r ) & = e ^ { i \theta _ { k } } C _ { 2 } ^ { z } u _ { k } ^ { * } ( C _ { 2 } ^ { z } r ) \\ & = e ^ { i \theta _ { k } } ( - u _ { k } ^ { x * }, - u _ { k } ^ { y * }, u _ { k } ^ { z * } ) | _ { ( - x, - y, z ) }$$

Here θ k is an arbitrary phase factor. Meanwhile, we are free to multiply u k with any phase factor, and it remains a valid eigenfunction. For our purpose here, we explicitly choose the phase factor of u k such that e iθ k = -1 for all k . With this choice, we can average over x and y to get c k ( ) = c ∗ ( k ) for all k . That is, the polarization vector c k ( ) is purely real.

Using the fact that systems in region I also have the up-down mirror symmetry σ z , namely glyph[epsilon1] ( x, y, z ) = glyph[epsilon1] ( x, y, -z ), we can link the radiation loss above and below the photonic crystal slab denoted by c ↑ and c ↓ . At any k point, u k ( r ) and σ z u k ( σ z r ) are both eigenfunctions of Θ ( ) with the same eigenvalue, so k r

$$u _ { k } ( { \mathbf r } ) & = e ^ { i \theta _ { k } } \sigma _ { z } u _ { k } ( \sigma _ { z } { \mathbf r } ) \\ & = e ^ { i \theta _ { k } } ( u _ { k } ^ { x }, u _ { k } ^ { y }, - u _ { k } ^ { z } ) | _ { ( x, y, - z ) }$$

with θ k being an arbitrary phase factor (not to be confused with the one in Eq. (S.1)). Since σ 2 z = 1, we can apply Eq. (S.2) twice to show that e iθ k = ± 1. Averaging over x and y , we see that c ↑ = ± c ↓ .

After using these two symmetries, the number of independent real variables in all radiation coefficients c ↑ ↓ , x,y has been reduced from 8 to 2. Given that the number of independent tuning parameters is also 2: ( k , k x y ), we are able to get stable BICs. Note that the combination of C T z 2 and σ z is just one sufficient condition for stable BICs in photonic crystal slabs. There might be other different choices of symmetries. For example, PT and σ z is equivalent to C T z 2 and σ z , where P is the inversion operator. Also, the requirement of σ z is not necessary when there is leakage to one direction only (such as BICs on the surface of a photonic bandgap structure 12 ).

In region II , stable BICs at C z 2 -invariant k points can be found. Systems in this region have C z 2 symmetry, namely glyph[epsilon1] ( x, y, z ) = glyph[epsilon1] ( -x, -y, z ). k points are C z 2 -invariant when -k = k + G , with G being a reciprocal lattice vector. A short derivation shows that at any k point, u k ( r ) and C z 2 u -k ( C z 2 r ) are both eigenfunctions of Θ k ( r ) with the same eigenvalue, so

$$u _ { k } ( r ) = e ^ { i \theta _ { k } } C _ { 2 } ^ { z } u _ { - k } ( C _ { 2 } ^ { z } r ),$$

with θ k being an arbitrary phase factor (not to be confused with the two phase factors above). At these high-symmetry k points, using Bloch theorem we know: u -k = u k + G = u k , so we can apply Eq. (S.3) twice to get e iθ k = ± 1. When this factor is +1, we can average over x and y to see that c k ( ) = 0, corresponding to a BIC at this C z 2 -invariant k point.

In region III , both kinds of BICs can be found, where C z 2 , T and σ z are all present. All our numerical examples are within this region to make it easier to understand the relation and interaction between different types of BICs.

## B. Consequences of topological charge conservation

Since topological charge is a conserved quantity, there are a few consequences and restriction on the evolution of BICs. First, BICs are stable as long as the system retains required symmetries; however, perturbations that break these two required symmetries eliminate the existence of BICs. When C T z 2 symmetry is broken, the coefficients ( c x and c y ) require complex components, meaning the radiation becomes elliptically polarized instead of linearly polarized. When σ z symmetry is broken, the coefficients c ↑ ↓ , x,y are still real numbers, but radiation towards the top and towards the bottom become separate degrees of freedom and so they do not vanish simultaneously in general. Second, when BICs collide into each other

in the moment space, the sum of all topological charges they carry remains the same before and after the collision.

## C. Example of charge -2

We consider the lowest-frequency TE-like mode of a photonic crystal slab with a hexagonal lattice of cylindrical air holes (shown in Fig. S2a). The refractive index of the slab is n = 1 5; the air-hole diameter is 0 5 . . a ; and the thickness of the slab is 0.5a, where a is the lattice constant. This system has C z 6 symmetry. Normalized lifetime plot indicates a BIC at the center of the Brillouin zone shown in Fig. S2b. The polarization vector field characterizes the BIC carrying charge -2 shown in Fig. S2c. Charge -2 can also be understood from the double degeneracy of both nodal lines of c x (green) and c y (red), shown in the inset of Fig. S2c. All four nodal lines are pinned at Γ point stabilized by the C 6 symmetry.

## D. BICs related by point group symmetries have the same topological charges

Here, we prove that when the structure has a certain in-plane point group symmetry R (namely, glyph[epsilon1] ( r ) = glyph[epsilon1] ( R r ); R can be a combination of rotation and reflection on the x y -plane) and when the band has no degeneracy, a BIC at k indicates there is another BIC at R k with the same topological charge. The assumption here is that the eigenfunctions u k at different k points already have their phases chosen to ensure the reality of c k ( ), and the signs of u k at different k points have been chosen such that u k is continuous with respect to k (so that a small change in k leads to a small change in u k ).

We start by relating the eigenfunction at k and the eigenfunction at R k . Let u k be an eigenfunction of operator Θ k . Since the system is invariant under transformation R , we know ˆ O R u k is an eigenfunction of Θ R k , so in the absence of degeneracy, we can write ˆ O R u k = α k u R k , where α k is some number. The number α k must have unit magnitude (due to the normalization of u k and u R k ) and must be real-valued (because c k ( ) is real-valued), so it can only take on discrete values of ± 1. Also, α k must be a continuous function of k since u k is continuous with respect to k . Since α k is both discrete-valued and continuous, it must be a constant. Then, we may denote this constant with its value at the Γ point, as α k = α Γ . Note that R Γ = Γ , so we can determine coefficient α Γ using the mode profile:

ˆ O R u Γ = α Γ u Γ . In conclusion, we have u R k = α O Γ ˆ R u k .

Now we consider how the angle φ ( R k ) is related to φ ( k ). The vector field u k transforms under the rotation operator as ( ˆ O R u k )( r ) = R u k ( R -1 r ), so averaging over x and y we get 〈 ˆ O R u k 〉 = R〈 u k 〉 . Let P be the operator that projects a 3D vector onto the x y -plane, namely P r = r -( r · ˆ)ˆ; it commutes with z z R , since it does not alter the x or y component. Then c k ( ) = P 〈 u k 〉 , and

$$c ( \mathcal { R } k ) = P \langle \mathbf u _ { \mathcal { R } k } \rangle = P \langle \alpha _ { \Gamma } \hat { O } _ { \mathcal { R } } \mathbf u _ { k } \rangle = \alpha _ { \Gamma } \mathcal { R } P \langle \mathbf u _ { k } \rangle = \alpha _ { \Gamma } \mathcal { R } c ( k ).$$

So, the polarization vector at the the transformed k point is simply the original polarization vector transformed and times ± 1. So, the angle of the polarization vector only changes by a constant in the case of proper rotations (where det R = 1); in the case of improper rotations (where det R = -1), it also changes sign. So, in general, we can write

$$\phi ( \mathcal { R } { \mathbf k } ) = ( \det \mathcal { R } ) \phi ( { \mathbf k } ) + c$$

with c being a constant depending on R and α Γ . It follows that ∇ R k φ ( R k ) = (det R R∇ ) k φ ( k ), so the topological charge at R k is

$$\text{gical charge at } \mathcal { R } k & \text{is} \\ q _ { \mathcal { R } k } & = \frac { 1 } { 2 \pi } \oint _ { C _ { \mathcal { R } k } } \nabla _ { k ^ { \prime \prime } } \phi ( k ^ { \prime \prime } ) \cdot d k ^ { \prime \prime } \\ & = \frac { 1 } { 2 \pi } \oint _ { \mathcal { R } ^ { - 1 } C _ { \mathcal { R } k } } \nabla _ { \mathcal { R } k ^ { \prime } } \phi ( \mathcal { R } k ^ { \prime } ) \cdot \mathcal { R } d k ^ { \prime } \\ & = \frac { 1 } { 2 \pi } ( \det \mathcal { R } ) \oint _ { C _ { \mathcal { R } k } } \nabla _ { \mathcal { R } k ^ { \prime } } \phi ( \mathcal { R } k ^ { \prime } ) \cdot \mathcal { R } d k ^ { \prime } \\ & = \frac { 1 } { 2 \pi } ( \det \mathcal { R } ) ^ { 2 } \oint _ { C _ { \mathcal { R } } } \mathcal { R } \nabla _ { k ^ { \prime } } \phi ( k ^ { \prime } ) \cdot \mathcal { R } d k ^ { \prime } \\ & = q _ { k }, \\ \\ - \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

where C R k is a closed simple path that is centered on R k and loops in counterclockwise direction, R -1 C R k is this loop transformed by R -1 (which centers on k in counterclockwise direction if R is a proper rotation, or in clockwise direction if R is improper), and C k is this transformed loop traversed in counterclockwise direction.

In conclusion, we have proven that if a system has certain point group symmetry R , then the topological charges carried by the BIC at k and at R k on a singly degenerate band have to be the same. This conclusion agrees with all examples in Figs. 2-4.

## E. Allowed charges at Γ in systems with different symmetries

Allowed topological charges at high symmetry k points can be determined by the field eigenvalues of the rotational symmetry of a system. For systems with m -fold rotational symmetry, we can first determine the relationship between polarization direction at wavevector k and at rotated wavevector R k ( φ ( k ) and φ ( R k )) using Eq. (S.4). Since the wavevector gets back to its original point if applying this rotation m -times: R m k = k , we can then apply this relationship m times and get how many times the polarization vector rotates around the center of the Brillouin zone. From there, we categoraize all possible charges allowed at Γ as shown in Table S1. Allowed charges depend on two factors. The first one is which symmetry representation the band belongs to. The second one is the degeneracy of nodal lines at Γ, because more nodal lines intersecting at the same point usually leads to more oscillations in color and thus higher topological charges. This factor is reflected by the integer number n , depending on the number of equivalent Γ points at this frequency 30 . Note that only singly degenerate bands are considered in this Letter, having no crossing with other bands in the bandstructures, as can be seen in Table S1. Further research directions may include BICs on degenerate bands, as well as the search of BICs with higher-order and potentially fractional topological charges.

TABLE I. Allowed stable topological charges at Γ for singly degenerate bands. A B ( ) corresponds to modes of different representations of the symmetry operator 32 . Note that only singly degenerate representations of symmetry operators are included in here.

| Symmetries   | Representation   | Charges   | Allowed n     | Allowed charges   |
|--------------|------------------|-----------|---------------|-------------------|
| C 2          | A                | ± 1 +2 n  | 0             | ± 1               |
| C 2          | B                | 0 +2 n    | 0             | 0                 |
| C 3          | A                | 1 +3 n    | 0 , ± 1 , ... | +1,+4,-2,...      |
| C 4          | A                | 1 +4 n    | 0 , ± 1 , ... | +1,+5,-3,...      |
| C 4          | B                | - 1 +4 n  | 0 , ± 1 , ... | -1,-5,+3,...      |
| C 6          | A                | 1 +6 n    | 0 , ± 1 , ... | +1,+7,-5,...      |
| C 6          | B                | - 2 +6 n  | 0 , ± 1 , ... | -2,+4,-8,...      |

## a Field decomposition

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0237v1-figure-3c8422a1-8d35-486b-b5fc-735f6e30f754.jpg)

FIG. 1. Stable bound states in the continuum (BICs) as vortex centers of polarization vectors. a , Schematics of radiation field decomposition for resonances of a slab structure. The spatially-averaged Bloch part of the electric field 〈 uk 〉 is projected onto the x y -plane as the polarization vector c = ( c x , c y ). A resonance turns into a BIC if and only if c x = c y = 0. b , Schematic illustration for the nodal lines of c x (green) and of c y (red) in a region of k space near a BIC. The direction of vector c (shown in arrows) becomes undefined at the nodal line crossing, where a BIC is found. c , Two possible configurations of the polarization field near a BIC. Along a closed loop in k -space containing a BIC (loop goes in counterclockwise direction, 1 → → → 2 3 4), the polarization vector either rotates by angle 2 π (denoted by topological charge q = +1) or rotates by angle -2 π (denoted by topological charge q = -1). Different regions of the k space are colored in four gray-scale colors according to the signs of c x and c y . In this way, a BIC happens where all four gray-scale colors meet, and charge q = +1 corresponds to the color changing from white to black along the counterclockwise loop C , and charge q = -1 corresponds to the color changing from black to white.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0237v1-figure-f8f5cd09-c4ff-479a-bfab-772d031c7f67.jpg)

FIG. 2. Characterization of BICs using topological charges. a , Calculated radiative quality factor Q of the TM 1 band on a square-lattice photonic crystal slab (as in ref. 13), plotted in the first Brillouin zone. Five BICs can be seen. b , Directions of the polarization vector field reveal vortices with topological charges of ± 1 at each of the five k points. The area shaded in blue indicates modes below the lightline and thus bounded by total internal reflection. c , Nodal lines and gray-scale colors of the polarization vector fields (same coloring scheme as in Fig. 1c).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0237v1-figure-51ab4794-ca44-4da4-bd28-9d181fbb7997.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0237v1-figure-076aea9d-54eb-476e-b413-b90f720b21ec.jpg)

FIG. 3. Evolution of BICs and conservation of topological charges. a , Schematic drawing of a photonic crystal slab with one-dimensional periodicity in x and is infinitely long in y . b , Calculated TM-like band structure along k x axis and along k y axis. The area shaded in yellow indicates the light cone, where there is a continuum of radiation modes in the surrounding medium. c, d , An example showing topological charges with the same sign bouncing off each other. As the slab thickness h decreases, the two BICs with charge +1 move along the k x axis, meet at the origin, and then deflect onto the k y axis. This can be understood from the conservation of topological charges or from the evolution of nodal lines. e, f , An example of topological charge annihilation happening on the lowest-frequency TE-like band of the same structure with different slab thicknesses. As the slab thickness h decreases, two BICs with charge -1 meet with a BIC with charge +1 at the origin. These three BICs annihilate to yield one BIC with charge -1, as governed by charge conservation.

- [1] von Neumann, J. &amp; Wigner, E. Uber merkw¨rdige diskrete eigenwerte. ¨ u Phys. Z. 30 , 465-467 (1929).
- [2] Friedrich, H. &amp; Wintgen, D. Interfering resonances and bound states in the continuum. Phys. Rev. A 32 , 3231-3242 (1985).
- [3] Bulgakov, E. N. &amp; Sadreev, A. F. Bound states in the continuum in photonic waveguides inspired by defects. Phys. Rev. B 78 , 075105 (2008).
- [4] Corrielli, G., Della Valle, G., Crespi, A., Osellame, R. &amp; Longhi, S. Observation of surface states with algebraic localization. Phys. Rev. Lett. 111 , 220403 (2013).
- [5] Weimann, S. et al. Compact surface Fano states embedded in the continuum of waveguide arrays. Phys. Rev. Lett. 111 , 240403 (2013).
- [6] Evans, D. V., Levitin, M. &amp; Vassiliev, D. Existence theorems for trapped modes. J. Fluid Mech. 261 , 21-31 (1994). URL http://journals.cambridge.org/article\_ S0022112094000236 .
- [7] Plotnik, Y. et al. Experimental observation of optical bound states in the continuum. Phys. Rev. Lett. 107 , 183901 (2011).
- [8] Lee, J. et al. Observation and differentiation of unique highQ optical resonances near zero wave vector in macroscopic photonic crystal slabs. Phys. Rev. Lett. 109 , 067401 (2012).

- [9] Porter, R. &amp; Evans, D. Embedded rayleigh-bloch surface waves along periodic rectangular arrays. Wave Motion 43 , 29 - 50 (2005).
- [10] Marinica, D. C., Borisov, A. G. &amp; Shabanov, S. V. Bound states in the continuum in photonics. Phys. Rev. Lett. 100 , 183902 (2008).
- [11] Liu, V., Povinelli, M. &amp; Fan, S. Resonance-enhanced optical forces between coupled photonic crystal slabs. Opt. Express 17 , 21897-21909 (2009).
- [12] Hsu, C. W. et al. Bloch surface eigenstates within the radiation continuum. Light: Science &amp; Applications 2 , e84 (2013).
- [13] Hsu, C. W. et al. Observation of trapped light within the radiation continuum. Nature 499 , 188-191 (2013).
- [14] Zhan, Q. Cylindrical vector beams: from mathematical concepts to applications. Adv. Opt. Photon. 1 , 1-57 (2009).
- [15] Joannopoulos, J. D., Johnson, S. G., Winn, J. N. &amp; Meade, R. D. Photonic Crystals: Molding the Flow of Light (Princeton University Press, 2008), 2 edn.
- [16] Hirose, K. et al. Watt-class high-power, high-beam-quality photonic-crystal lasers. Nature Photon. 8 , 406-411 (2014).
- [17] Ko, D.-H. et al. Photonic crystal geometry for organic solar cells. Nano Letters 9 , 2742-2746 (2009).
- [18] Wierer, J. J., David, A. &amp; Megens, M. M. III-nitride photonic-crystal light-emitting diodes with high extraction efficiency. Nature Photonics 3 , 163-169 (2009).
- [19] Ganesh, N. et al. Enhanced fluorescence emission from quantum dots on a photonic crystal surface. Nature Nanotech. 2 , 515-520 (2007).
- [20] Yanik, A. A. et al. Seeing protein monolayers with naked eye through plasmonic fano resonances. Proc. Natl. Acad. Sci. U.S.A. 108 , 11784-11789 (2011).
- [21] Fan, S. &amp; Joannopoulos, J. D. Analysis of guided resonances in photonic crystal slabs. Phys. Rev. B 65 , 235112 (2002).
- [22] Pacradouni, V. et al. Photonic band structure of dielectric membranes periodically textured in two dimensions. Phys. Rev. B 62 , 4204-4207 (2000).
- [23] Ochiai, T. &amp; Sakoda, K. Dispersion relation and optical transmittance of a hexagonal photonic crystal slab. Phys. Rev. B 63 , 125107 (2001).

- [24] Yang, Y., Peng, C., Liang, Y., Li, Z. &amp; Noda, S. Analytical perspective for bound states in the continuum in photonic crystal slabs. Phys. Rev. Lett. 113 , 037401 (2014).
- [25] Mermin, N. D. The topological theory of defects in ordered media. Reviews of Modern Physics 51 , 591 (1979).
- [26] Donnelly, R. J. Quantized vortices in helium II , vol. 2 (Cambridge University Press, 1991).
- [27] Berry, M. The adiabatic phase and pancharatnam's phase for polarized light. Journal of Modern Optics 34 , 1401-1407 (1987).
- [28] Ma˜ nes, J. L., Guinea, F. &amp; Vozmediano, M. A. Existence and topological stability of fermi points in multilayered graphene. Physical Review B 75 , 155424 (2007).
- [29] Kurosaka, Y. et al. On-chip beam-steering photonic-crystal lasers. Nature Photonics 4 , 447450 (2010).
- [30] Iwahashi, S. et al. Higher-order vector beams produced by photonic-crystal lasers. Opt. Express 19 , 11963-11968 (2011).
- [31] Kitamura, K., Sakai, K., Takayama, N., Nishimoto, M. &amp; Noda, S. Focusing properties of vector vortex beams emitted by photonic-crystal lasers. Opt. Lett. 37 , 2421-2423 (2012).
- [32] Sakoda, K. Optical properties of photonic crystals , vol. 80 (Springer, 2005).

FIG. 4. Generation of BICs. a , Schematic drawing of a photonic crystal slab with twodimensional periodicity. b , Generation of BICs on the TE 1 band when the slab thickness h is increased. Each time, four pairs of BICs with charges ± 1 are generated simultaneously, consistent with the charge conservation and C 4 v symmetry. Insets show the locations of BICs in the k space and their corresponding topological charges for h/a = 1 0 1 2 1 35 1 8 . , . , . , . , and 2 4. As the slab thick-. ness increases, the BICs move outward and eventually fall below the light line into the area shaded in dark blue.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0237v1-figure-5b7017fa-1e6f-4398-b799-518ad4a8ee92.jpg)

FIG. S1. Symmetry requirements for BICs. Systems in the blue circle are invariant under operators C T z 2 and σ z , where stable BICs at arbitrary wavevectors can be found. In the red circle, where C z 2 is a symmetry of the system, robust BICs can be found at high-symmetry wavevector points. Here, high-symmetry wavevectors mean C z 2 -invariant ones, while arbitrary wavectors are not necessarily C z 2 -invariant. In the overlapping area (region III ), both types BICs can be found. All numerical examples in this Letter are within region III .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0237v1-figure-8a501513-0f7e-44c3-b7ef-a8a263e2d790.jpg)

FIG. S2. Stable BIC with topological charge -2. a , Schematic drawing of the photonic crystal slab. b , Q plotted in the first Brillouin zone, showing a BIC at the Γ point. c , Polarization vector field characterizes the BIC with a stable topological charge of -2, as can be shown from double degeneracies of both nodal lines.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0237v1-figure-e621e52f-85ed-437a-8cf0-dd659ec0c941.jpg)
Bo Zhen Chia Wei Hsu Ling Lu A. Doug Stone Marin Soljacic
Topological nature of bound states in the radiation continuum
Bound states in the continuum (BICs) are unusual solutions of wave equations
describing light or matter: they are discrete and spatially bounded, but exist
at the same energy as a continuum of states which propagate to infinity. Until
recently, BICs were constructed through fine-tuning parameters in the wave
equation or exploiting the separability of the wave equation due to symmetry.
More recently, BICs that that are both robust and not symmetry-protected
(accidental) have been predicted and experimentally realized in periodic
structures; the simplest such system is a periodic dielectric slab, which also
has symmetry-protected BICs. Here we show that both types of BICs in such
systems are vortex centers in the polarization direction of far-field
radiation. The robustness of these BICs is due to the existence of conserved
and quantized topological charges, defined by the number of times the
polarization vectors wind around the vortex centers. Such charges can only be
generated or annihilated by making large changes in the system parameters, and
then only according to strict rules, which we derive and test numerically. Our
results imply that laser emission based on such states will generate vector
beams.