Robust Topological Edge States in C6 Photonic Crystals

1. Introduction

Photonic crystals (PCs) are systems whose electromagnetic features periodically modulate in the real space [1]. In particular, two-dimensional (2D) photonic crystals can present topological behavior by introducing a perturbation in the Hamiltonian of the system in order to obtain a non-zero Chern number [2,3]. 2D topological photonic crystals are currently in the scientific limelight not only for possessing tremendous technological potential but also for having opened several avenues of basic science exploration [4,5,6,7,8,9,10,11]. It is known that when the photonic crystal is in a nontrivial topological phase interesting phenomena can emerge. We can highlight photonic analogies of the quantum Hall effect [12,13], higher-order topological photonic crystals [14,15,16], and the quantum valley Hall effect [17,18,19,20].

It is known from the literature that phononic crystals in a triangular lattice, with $C_6$ point symmetry group, present a doubly degenerate Dirac cone at the $\Gamma$ point of the Brillouin zone. Moreover, the degenerate bands present $p-$ and $d-$ waves orbitals in the band structure of the transverse magnetic (TM) polarization [21,22]. However, if a perturbation is introduced in the crystal the double degenerescence is broken and a complete topological bandgap is opened in the band structure [23]. In addition, the edge-bulk correspondence guarantees that edge states will emerge inside the bandgap and they will be protected by the topology [24,25]. Then, the topological protection ensures that the edge states are localized and they are robust against defects, disorder and reflection [26,27].

In this work we investigate the propagation of electromagnetic waves, band structure and topological features of a two-dimensional topological photonic crystal composed of six dielectric cusped-oval-shaped (COS) rods. We induce a bandgap in the system by introducing a perturbation in the rods’ orientation angle, lifting the double degeneracy at the $\Gamma$ point of the Brillouin zone. In the following, we will show that two edge modes emerge in the induced bandgap. This paper is organized as follows. In Section 2 we introduce our system and explore its features. In Section 3 we introduce the perturbation and study the topological behavior associated with positive and negative perturbations. In Section 4 we study the emergence of edge states around the interface between the topological and trivial photonic crystals. The robustness of the edge states is addressed in Section 5. Finally, in Section 6 we give concluding remarks about this work.

2. The Photonic System

Consider a triangular lattice with six dielectric rods ($\epsilon =13$) per site surrounded by air. The unit cell is composed of dielectric COS rods, with $a=1\mu m$, as we can see in Figure 1. The COS rods are built taking into account two cylindrical rods of radius $r=0.17a$, which are shifted by a distance $d=\pm a/15$. Then, we consider the intersection area between the shifted cylindrical rods. Finally, the COS rods are placed in their locations which are distant $R=a/3$ from the unit cell’s center (see Figure 1), in order to obtain a 2D system with $C_6$ symmetry point group.

Each individual COS rod in the unit cell can be rotated around its respective center by the orientation angle ${\Phi }_i$, as illustrated in Figure 1. The rods present an anisotropic angular orientation which can vary spatially. Therefore, we may expect that an angular perturbation will lead the system to undergoes topological phase transitions from trivial to nontrivial domain [28]. A perturbation $\varphi$ is introduced in the orientation angle ${\Phi }_i$, so that we can write the orientation angle of the i-th COS rod as

$${\Phi }_i=(2i-1)\frac{\pi }{6}+{\varphi }_0+\varphi .$$

(1)

Here, $i=1,2,..,6$ is the rod index, ${\varphi }_0$ is the initial unperturbed angle and $\varphi$ is the angular perturbation introduced in our system. We remark that, because of the symmetry, the perturbation $\varphi$ has a period $\pi$, which means that the $\varphi$ and $\varphi +\pi$ induce the same topological behavior in the system. Here we assume ${\varphi }_0=\pi /4$ (see Figure 2(a)).

It is well known from the literature that triangular lattices with six “artificial atoms” have two 2D irreducible representations in the $C_6$ point symmetry group, which are associated with the symmetry of the triangular lattice [29]. As a consequence, a doubly degenerate Dirac cone appears at the Brillouin zone center [30]. The degenerate bands are pseudospin states which are related to $p_x$ ($p_y$) and $d_{xy}$ ($d_{x^2-y^2}$) orbitals, corresponding to odd and even parity in the real space, respectively [21]. The pseudospin states can be written as [28]

$$p_{\pm }=\frac{1}{\sqrt{2}}(p_x\pm i{p}_y)$$

(2)

and

$$d_{\pm }=\frac{1}{\sqrt{2}}(d_{x^2-y^2}\pm i{d}_{xy}).$$

(3)

Considering the unperturbed photonic crystal ($\varphi =0$), we calculate the band structure for the TM modes ($E_z,H_x,H_y\ne 0$) using the COMSOL Multiphysics software [31] which is based on the finite element method (FEM). The band structure is shown in Figure 2(a). We can observe a doubly degenerate Dirac cone, at the $\Gamma$ point, between the second and fifth bands, which is a consequence of the $C_6$ symmetry group of the system.

In Figure 2(b) we plot the electric field along the z direction ($E_z$), at the Dirac point, with $\omega a/\left(2\pi c\right)=0.437$. We found four states that are related to dipole and quadrupole modes. More specifically, $p_x$ and $p_y$ orbitals are dipole modes, while $d_{xy}$ and $d_{x^2-y^2}$ orbitals are quadrupole modes. In the next section we introduce a nonzero perturbation in order to lift the doubly degenerescence and induce a complete photonic bandgap in the band structure.

3. Topological Phase Transition

In Section 2 we have found a doubly degenerate Dirac cone at the $\Gamma$ point, as well as, orbitals p-like and d-like which are associated to the degenerate bands. Let us now study the consequences of considering a nonzero perturbation $\varphi$ in the rods’ orientation angle. In order to illustrate the effects of the angular perturbation, we evaluate $\omega a/\left(2\pi c\right)$vs$\varphi$ for the second, third, fourth, and fifth bands, as we can see in Figure 3. We can observe that as $\mid \varphi \mid$ increases from 0, the bandgap width monotonically increases reaching its maximum at $\varphi =-\pi /4$ and $\varphi =\pi /4$. On the contrary, as $\mid \varphi \mid$ increases from $\varphi =-\pi /4$ and $\varphi =\pi /4$, the bandgap width monotonically decreases until the doubly degenerate Dirac cone is recovered at $\varphi =-\pi /2$ and $\varphi =\pi /2$. This is a consequence, as mentioned before, that the angular perturbation $\varphi$ has a period $\pi$.

We next illustrate the opening of the bandgap in the band structure for two values of $\varphi$, one positive and other negative. Figure 4 shows the band structure corresponding to $\varphi =-\pi /4$ and $\varphi =\pi /4$. The perturbation opens a gap between $\omega a/\left(2\pi c\right)=0.432$ and $\omega a/\left(2\pi c\right)=0.4408$, for the positive case, and between $\omega a/\left(2\pi c\right)=0.431$ and $\omega a/\left(2\pi c\right)=0.4415$, for the negative case, corresponding to a gap–midgap ratio [1] of $\Delta \omega /{\omega }_m=0.0202$ and $\Delta \omega /{\omega }_m=0.0241$, respectively. One can observe that a complete bandgap is opened for both cases. Note that once the bandgap is opened, edge states can appear inside the gap as we will see later in this paper.

In order to study the topological behavior close to the $\Gamma$ point, we can write an effective Hamiltonian by using the $\mathbf{k}$.$\mathbf{p}$ perturbed model from which we can obtain the Chern number [32,33,34]

$$C_{\pm }=\pm \frac{1}{2}[sign\left(B\right)+sign\left(M\right)].$$

(4)

Here B is the diagonal term of the effective Hamiltonian close to the $\Gamma$ point which is essentially negative. Also, $M=({\omega }_d-{\omega }_p)/2$, where ${\omega }_d$ and ${\omega }_p$ are the eigenmodes of orbit d and orbit p, respectively [35]. The eigenmode ${\omega }_p$ is related to the double degenerate dipole states of $p_{\pm }$, while ${\omega }_d$ is related to the double degenerate quadrupole states of $d_{\pm }$[32]. If ${\omega }_p<{\omega }_d$ then $M>0$, hence $C_{\pm }=0$ and the photonic crystal is topologically trivial. However, if ${\omega }_p>{\omega }_d$ then $M<0$, hence $C_{\pm }=\pm 1$ and the photonic crystal is in a topological phase. Therefore, the inversion of the bands between the degenerate modes at the $\Gamma$ point leads to the topological phase transition [36].

For this work we consider the inversion of the bands that occurs between $\varphi =-\pi /4$ and $\varphi =\pi /4$ for the degenerate modes, as shown in Figure 5. We can observe that for the positive case the frequency of the dipole modes is lower than the frequency of the quadrupole modes, while for the negative case the frequency of the dipole modes is higher than the frequency of the quadrupole modes. Therefore, for $\varphi =\pi /4$ we obtain $C_{\pm }=0$, which corresponds to a trivial photonic crystal, and for the case $\varphi =-\pi /4$ we obtain $C_{\pm }=\pm 1$, which corresponds to a topological photonic crystal.

We can obtain important information about the topological behavior of the photonic bands from the electromagnetic (EM) energy density distribution in the real space. The general idea is that the EM energy density has peaks that are shifted towards the maximal localized Wyckoff points (WP) at which the Wannier functions (WF) of the system are centered. This is a consequence of the relationship between the Wilson-loop (WL) operator and the maximally-localized WF. The spectrum of the WL operator is a useful method for characterizing the topological phases of physical systems. First, we write the WL operator as a path-ordered exponential of the Barry fase which is defined by [37]

$${\mathcal{W}}_{mn}\left(l\right)=\mathcal{P}{e}^{-i{\int }_l{\mathbf{A}}_{mn}\left(\mathbf{k}\right).d\mathbf{l}}.$$

(5)

Here $\mathcal{P}$ is the path ordering operator, and ${\mathbf{A}}_{mn}$ is the Berry connection for $m=n$. Second, it is known from the literature that there is a connection between the WL operator and the WF. The WF, which are defined as a Fourier transformation of the Bloch states, can be written as [38]

$$w_{i\mathbf{R}}\left(\mathbf{r}\right)\propto {\int }_{BZ}{e}^{-i\mathbf{k}.\mathbf{R}}\sum _j{U}_{ij}^{\mathbf{k}}{\psi }_{j\mathbf{k}}\left(\mathbf{r}\right).$$

(6)

Here $U_{ij}^{\mathbf{k}}$ denotes the mixing matrix, which represents the mixing of the Bloch modes in the reciprocal space.

When we take into account the maximally-localized WF, the mixing matrix takes values to minimize the delocalization of the wave function according to the eigenvalue of the WL operator [39]. The sum of the phases of the operator’s eigenvalues provides a straight line in the Brillouin zone that corresponds to the expectation value of the projected position operator calculated over the maximally-localized WF [40]. For trivial systems, the WL spectrum does not present winds and the Chern number is zero. Moreover, the maximally-localized Wannier functions are exponentially localized in specific points of the real space. On the other hand, for topological systems the WL eigenvalues present winds and the maximally-localized WFs are not exponentially localized, but polynomially localized in the real space between consecutive unit cells [38]. Therefore, we can identify different topological phases by looking for the positions of the maxima of the WF in the real space.

As mentioned before, the peaks of the EM energy density are shifted towards the maximal localized WP at which the WF are centered. Therefore, the EM density energy can be a useful tool to probe the topological features of the system. We can obtain the EM energy density from the local density of states for a set of connected bands $\Lambda$ and considering the TM polarization as [38]

$$n_{\Lambda }\left(\mathbf{r}\right)=\frac{6}{\pi S}\sum _{n\in \Lambda }{\int }_{BZ}d\mathbf{k}\epsilon \left(\mathbf{r}\right)\mid {\mathbf{E}}_{\mathbf{k}n}\left(\mathbf{r}\right){\mid }^2.$$

(7)

Here S is the area of the unit cell in the real space and ${\mathbf{E}}_{\mathbf{k}n}\left(\mathbf{r}\right)$ is the electric field of the n-th band. It is possible to find the total EM energy density by summing all sets of bands, i.e., $n\left(\mathbf{r}\right)={\sum }_{\Lambda }{n}_{\Lambda }\left(\mathbf{r}\right)$. Moreover, the total energy density is written in terms of the electromagnetic field WF, i.e., ${\mathbf{E}}_{\mathbf{k}n}^w\left(\mathbf{r}\right)$ [41]:

$$n\left(\mathbf{r}\right)=\frac{6}{\pi S}\sum _{\Lambda }\sum _{n\in \Lambda }\sum _{\mathbf{R}}\epsilon \left(\mathbf{r}\right)\mid {\mathbf{E}}_{\mathbf{k}n}^w\left(\mathbf{r}\right){\mid }^2.$$

(8)

Since the WFs are linked to the WL, Eq. 8 allows us to indirectly identify the topological behavior by calculating the EM density, and study its maximal localization in the real space, without directly evaluating the WL [38].

For the photonic crystal considered in this work, we evaluate the EM energy density distribution for the positive and negative perturbation cases ($\varphi =\pi /4$ and $\varphi =-\pi /4$, respectively), as we can see in Figure 6. From Eq. 7 we can notice that the energy density depends on the permeability parameter $\epsilon \left(\mathbf{r}\right)$, which assume the values $\epsilon \left(\mathbf{r}\right)=1$ and $\epsilon \left(\mathbf{r}\right)=\epsilon$ for the background and rods, respectively. The difference between these two values of permeability makes the intensity of the EM density in the rods much higher than in the background. Thus, evaluating this quantity in the entire unit cell does not result in trustful data for the EM density energy localization in the background. In order to circumvent this problem, we separately evaluate the energy density for the rods and for the background. For the latter case, we consider $\epsilon \to 1$ for the rods.

Since the topological behavior of a bandgap can be defined as the sum of the topological behavior of the bands below that bandgap [8], we can focus on the set of bands below the bandgap at the $\Gamma$ point (see Figure 4). Figure 6 shows the EM energy density distribution in the unit cell and the maximal WP. In particular, results in the literature show that topological phases tend to present the associated EM energy density around the $3c$ WP in the edge of the unit cell [38]. From Figure 6(a) we conclude that the first band does not contribute to the topological features of the bandgap since it has a homogeneous energy distribution. Therefore, we can focus on the other bands below the bandgap, i.e., the second and third bands (see Figure 4(a)). As expected, we see that the maximal EM energy density in the rods is located in the regions around the $3c$ WP at the edge of the unit cell. The same behavior is observed for the background, but with the maximum EM energy density located right on the $3c$ WP. Therefore, we can infer that the negative perturbation leads our system to a nontrivial topological phase. On the other hand, in Figure 6(b) all the bands below the bandgap, i.e., the first, second and third bands (see Figure 4(b)) contribute to the topological features. We can observe that the maximum EM energy density is localized inside the rods but far from the maximal WP. Focusing on the background, we notice that the maximum is located between the rods, in the middle of the distance between the $1a$ WP and the edge of the unit cell. Thus, for the positive case, the perturbation leads the system to a trivial topological phase.

Let us quantify the localization illustrated in Figure 6. In order to do so, we set two lines: (i) one along the direction $({\mathbf{a}}_1+{\mathbf{a}}_2)$ and (ii) another along the direction $({\mathbf{a}}_1-{\mathbf{a}}_2)$. Next, we evaluate the EM energy density along those lines as shown in Figure 7. Figure 7(a) and (c) show the energy density $n_{\Lambda }/max\left[n_{\Lambda }\right]$ along the $({\mathbf{a}}_1-{\mathbf{a}}_2)$ and $({\mathbf{a}}_1+{\mathbf{a}}_2)$ directions, respectively, for $\varphi =-\pi /4$. On the other hand, Figure 7(b) and (d) show the energy density $n_{\Lambda }/max\left[n_{\Lambda }\right]$ along the $({\mathbf{a}}_1-{\mathbf{a}}_2)$ and $({\mathbf{a}}_1+{\mathbf{a}}_2)$ directions, respectively, for $\varphi =\pi /4$.

Analyzing Figure 7, we observe that bands 2 and 3, for the negative case, have the maximal EM energy density at the $3c$ WPs (red-solid line in Figure 7(a) and (c)), corresponding to polynomially localized WF. On the other hand, for the positive case, bands 1, 2, and 3 have the maximal EM energy density located in the region either between $1a$ and $2b$ WPs or between $1a$ and $3c$ WPs (red-solid line in Figure 7(b) and (d)), corresponding to exponentially localized WF. We can infer that the results about localization of the EM energy density reinforce our conclusion about the topological behavior for the negative perturbation case and the trivial behavior for the positive perturbation case. It is important to highlight that the results illustrated in Figure 7 are in complete agreement with the previous results obtained from the $\mathbf{k}.\mathbf{p}$ perturbed model, which we used to evaluate the Chern number.

4. Edge States

We have shown that the perturbation $\varphi$ opens a complete bandgap at $\Gamma$ point and we can identify two different phases, the topological phase for $\varphi <0$ and a trivial phase for $\varphi >0$. On the other hand, it is known from the literature that the edge-bulk correspondence guarantees that if we build a slab composed of two photonic crystals, with diferent topological invariants, i.e., Chern numbers, robust edge modes localized around the interface between the photonic crystals emerge inside the bandgap [34,42,43]. Those edge states are topologically protected and are robust against defects, disorder and allow transmission without any reflection, with no signifcant energetic loss [26,44,45]. Thus, in order to study the emergence of the edge states in our system, we built a supercell composed of 30 unit cells: 15 topological unit cells ($\varphi =-\pi /4$) and 15 trivial unit cells ($\varphi =\pi /4$). Next, we project the calculated band structure along the $\Gamma \to M$ direction as we can see in Figure 8.

From Figure 8 it is easy to identify two edge modes inside the gap and we notice that they emerge close to the $\Gamma$ point of the Brillouin zone. Those modes travel in opposite directions, and the traveling direction is reversed if we make $\mathbf{k}\to -\mathbf{k}$. In Figure 8(b) and Figure 8(c) we show the profile of $E_z$, the phase of $E_z$, and the Poynting vector for $\omega a/\left(2\pi c\right)=0.4349$ and $\omega a/\left(2\pi c\right)=0.4383$, respectively. We used $\mathbf{k}=k_x\widehat{x}=0.015k_0\widehat{x}$ with $k_0=4\pi /\left(\sqrt{3}a\right)$. Both modes are well localized at the interface and, by comparing $\omega a/\left(2\pi c\right)=0.4349$ and $\omega a/\left(2\pi c\right)=0.4383$, we see that the Poynting vectors have different directions, which is a confirmation of the pseudospin behavior of the edge states [18,46]. Moreover, focusing on the phase of the electric field $E_z$, we can identify that the mode with $\omega a/\left(2\pi c\right)=0.4349$ presents a clockwise polarization, while the mode with $\omega a/\left(2\pi c\right)=0.4383$ presents an anticlockwise polarization, corresponding to a pseudospin-down and pseudospin-up, respectively. Therefore, the pseudospin-up is associated with the interface state with group velocity and energy flux from the left to the right, while the pseudospin-down is associated with the interface state with group velocity and energy flux from the right to the left [46]. In the next section, we study the robustness and localization of the edge modes in our photonic system.

5. Robustiness of the Edge States

The topological protection guarantees that the propagating modes, associated with the edge states, have good robustness against defects, disorder and reflection at the interfaces [47,48,49]. So, in order to verify the edge states robustness, we build a $(35,20)a$ slab with a horizontal interface between topological ($\varphi =-\pi /4$) and trivial ($\varphi =\pi /4$) photonic crystals, respectively. We set a source of light on the left and a detector on the right, as shown in Figure 9(a). Then, we calculate the normalized electric field defined as $E_N=\sqrt{\mid E_z{\mid }^2}$, shown in Figure 9(b), for $\omega a/\left(2\pi c\right)=0.4353$. Next, we introduce small defects at the interface: (i) a small cavity removing a rod (Figure 9(c)), (ii) a bigger rod of size $d_0=0.3d$ (Figure 9(e)), (iii) a Ag rod (the yellow one in Figure 9(g)), and (iv) a disorder at the interface changing one negative perturbed unit cell for one positive perturbed unit cell (Figure 9(i)). We also introduce extensive defects: (i) a Z interface (Figure 10(a)), and an (ii) Omega (interface Figure 10(c)). The Normalized Electric Field $E_N$, corresponding to the defects mentioned above, is illustrated in Figure 9(b), Figure 9(d), Figure 9(f), Figure 9(h), Figure 9(j), Figure 10(b) and Figure 10(d), respectively. Besides the normalized electric field, the figures also show the Poynting vector (red arrows) around the defects, and around two points of the Z and Omega interfaces, as we can see in the zoom windows.

For the pristine interface case, as expected, the electric field and the Poynting vector are localized around the interface, while the energy flux is from the left to the right (see Figure 9(b)). Furthermore, comparing Figure 9(d) and Figure 9(b), we realize that there is no significant change in the electric field around the cavity’s position, and the Poynting vector is not captured by the cavity but just travels around it. Therefore, we can see in Figure 9(d) that the cavity created by removing a rod at the interface does not cause significant changes in the electrical field distribution and Poynting vector. Figure 9(f) corresponds to the bigger rod defect. Again, the presence of the defect does not cause changes in the edge mode, i.e., the electrical field distribution and Poynting vector are not affected. Despite the bigger rod, reflection does not occur and the energy flux is not affected by this defect at all. The electric field and the Poynting vector for the Ag rod defect are illustrated in Figure 9(h). We could expect major changes in this case because of the energy losses involved. However, as in the previous cases, light does not experience significant changes. In fact, the energy losses are much smaller than the transmittance, as we will see later. On the other hand, we can observe local changes in Figure 9(j). The disorder at the interface creates a different path for light. In this case the energy flux locally changes around the defect, but the global behavior does not change, i.e., the energy flux keeps flowing close to the interface and with no significant reflections. As the mode is localized around the interface, both the Poynting vector’s direction and electric field’s distribution deform to follow the new interface shape at the position of the defect. Thus, we realize that the electric field survives and remains localized around the interface despite any defect introduced in the interface. In addition, the Poynting vector just walks around the defect or ignores it, and the flux of energy remains unchanged. Much more impressive are the results for the extensive defects: the Z and Omega interfaces. For the Z interface case, the interface has two corners in which light faces two changes of direction that could cause reflections, but the Poynting vector and electric field just follow the interface’s contour and do not present any reflection in the corners (see Figure 10(b)). The Omega interface case is a more complex extensive defect. In fact, it has 6 corners which means that light faces six changes of direction! Despite the six corners, once again the Poynting vector and electric field follow the Omega interface’s contour and they do not present any reflection in the corners (see Figure 10(d)). Therefore, we can conclude that the edge states are robust against defects, disorders and reflection. This is guaranteed by the topological protection due to the Bulk-Edge correspondence. This behavior is agreement with previous works which studied edge modes in topological valley photonic crystals [50,51] and topological pseudospin photonic crystals [21,28,52,53].

In Figure 9 and Figure 10 is provided a very good qualitative piece of information on the topological protection and robustness of the edge modes. However, it is important to quantify the robustness of the edge modes. The quantitative information is provided by the calculation of the transmittance of light along the system. Thus, the calculated transmittance through the slabs schematized in Figure 9 and Figure 10 is plotted in Figure 11. We can observe in Figure 11 that transmittance changes very little when we introduce small defects like a cavity, a bigger rod, an Ag rod, and a disorder. In fact, the effect of the defects is indeed small, so that the transmittance remains around 1. It should be remarked that for the Ag rod defect we would expect some losses because of the metallic character of the defect. However, just minor changes occur in the transmittance which remains around 1. Let us discuss now the extensive defects. For the Z interface case we can observe a peak around 0.9, which is a 10% reduction in the transmittance in relation to the pristine interface. Despite this reduction in transmittance, the peak corresponding to the edge mode survives. A similar behavior is observed for the case of the Omega interface. Therefore, we can conclude that the transmittance peaks survive for the edge mode despite the small or extensive defects introduced in the interface. In short, our results show that the robustness of the interface mode is guaranteed against small and extensive defects in the interface, which means that light travels along the interface without changes in its energy flux, without reflection, and with minimal energy losses. Similar results were found in the literature with interface somehow modified: the peak can be reduced but the transmittance of light is at least $\ge 0.8$, which means that most part of light can travel through the considered system without reflection or absorption [11,54].

Before concluding, after studying the robustness of the edge modes, let us take look at the localization of the edge modes. In order to do so, we calculate $E_N/max\left[E_N\right]$ point-by-point on a line perpendicular to the interface (the blue line in Figure 12 and Figure 13). Therefore in Figure 12 is shown the intensity of $E_N/max\left[E_N\right]$ point-by-point along the line of the interface with no defect and with small defects: the interface with a cavity, the interface with a bigger rod, the interface with a disorder, and the interface with an Ag rod. The same is illustrated in Figure 13 for the extensive defects: the Z interface, and the Omega interface. We consider in all cases $\omega a/\left(2\pi c\right)=0.4353$. From Figure 12(a)-Figure 12(e) we can infer that the electric field is localized around the interface, regardless of the small defect considered, i.e., the electric field is near zero in rods far from the interface for both cases: the topological and trivial photonic crystals. This corresponds to a topological insulator behavior of our system because it does not present fields in the bulk and the electric field is different from zero only around the interface [46,50,55]. For the special cases of extensive defects, i.e., the Z and Omega interfaces, we set two perpendicular lines: (i) one at the interface before the first change of orientation, and (ii) another at the interface after the second change of orientation (see Figure 13(a) and Figure 13(b)). As for the small defect cases, the interface modes are well localized. Notice that the normalized electric field rapidly goes to zero when the profile moves away from the interface for both, the Z and Omega interfaces. In short, in all defect cases considered here, the localization of the modes allows the confinement of the light, and the interface works as a waveguide for the propagation of electromagnetic waves. Finally, all those features described here make this system a good candidate for topological wave guides and open up an opportunity for phototransport applications.

6. Conclusions

In this work we have proposed a two-dimensional photonic crystal composed of dielectric cusped-oval-shaped rods (COS rods), in a triangular lattice with $C_6$ point symmetry group, with the unit cell composed of six COS rods. The band structure of the system was obtained through the software COMSOL Multiphysics, which is based on the finite element method (FEM). Our results show that, as we introduced an angular perturbation $\varphi$ in the COS rods, a complete bandgap is opened in the band structure of the system. We have studied the topological phases of the system. It was found that for negative perturbations the photonic crystal is in a topological phase, while for positive perturbations the photonic crystal is in a trivial phase. The topological phase transition occurs for $\varphi =0$. We have also shown that the electromagnetic (EM) density energy distribution in the system can be used as a probe of the character of the band structure. This is because the EM density energy is proportional to the modulus square of the Wannier functions (WF) which are centered at the (maximal) localized Wyckoff positions (WP). Moreover, as a consequence of the perturbation, two edge states emerge inside the bandgap, which are localized around the interface and are topologically protected. Our results show that the edge modes are pseudospin modes traveling around the interface between the topological (negative perturbation) and trivial (positive perturbation) photonic crystals. They are topologically protected and robust against disorder and defects at the interface. For all defect cases considered in this work the transmittance is mildly affected. Therefore, the system studied is an excellent candidate for technological applications, once the flux of light can be controlled without any significant energetic loss or reflection. These features made the photonic crystal here proposed a good candidate to guide and confine light, with potential for phototransport applications.