Flat-band localization and self-collimation of light in photonic crystals

We investigate the optical properties of a photonic crystal (PC) composed of a quasi-one-dimensional flat-band lattice array through finite-difference time-domain simulations. The photonic bands contain flat bands (FBs) at specific frequencies, which correspond to compact localized states as a consequence of destructive interference. The FBs are shown to be nondispersive along the Г → X line, prohibiting optical transmission with incident light in x direction. On the other hand, the photonic band for the FB frequency is found to be dispersive along the Г → Y line, resulting in nonzero optical transmission. Such anisotropic optical response of the PC due to the FB localization of light in a single direction only results in a self-collimation of light propagation throughout the PC at the FB frequency.

In this paper, we investigate the photonic band structure and optical properties of a 2D PC, in the form of an array of quasi-1D flat-band lattices, using finite-difference time-domain (FDTD) simulations. We demonstrate that non-dispersive modes are found to be only along the Г → X line in the Brillouin zone, having dispersive nature along the Г → Y line. The CLSs are shown to be odd symmetric with respect to the transverse direction of the lattice. Contrary to the isotropic nature of the full-k flat bands, the anisotropic flat-band modes of our 2D PC allow the suppression of light propagation in the x direction. Notably, for a specific frequency corresponding to the CLSs, we reveal a self-collimation of light propagation for any focused excitation on our 2D PC, thanks to the nondispersive band in one direction only.

Photonic Bands and Transmittance Spectra
Here, we consider a 2D PC structure where infinitely long dielectric nanopillars are arranged periodically to form a specific lattice, as displayed in Fig. 1(a). Such a lattice structure is described by an array of quasi-1D tunable diamond (TD) lattices. Let us notice that the dielectric nanopillars in the array are infinitely long in the vertical direction, so that we may regard it as a 2D PC, assuming homogeneous electromagnetic fields along the z direction.
The unit cell of the 2D PC contains three nanopillars, which can be split into two sublattices: majority (A and B) and minority (C), as depicted in Fig. 1(a). All dielectric nanopillars are identical with radius r, and the period of the TD lattice is a. In this work, we suppose ε = 12.9 for the dielectric nanopillars, which indicates dielectric constant of GaAs as a widely used material in near-infrared ranges.
In order to determine the presence of the FBs, we investigate the optical properties of the 2D PC by numerical simulations in both frequency and time domains, using freely available MPB and MEEP software packages, respectively 37,38 . Numerical results from FDTD simulations, that we compute here, are convinced by using adequate resolution, 32 × 32 pixels per a unit cell. (see Supplementary Information) Fig. 1(b) shows the photonic band structure of the PC calculated in the 2D Brillouin zone, for TE polarization with light propagation in x − y plane (E x = E y = 0 but E z ≠ 0). The PC exhibits a wide band gap from f ≃ 0.5 to 0.67, a narrow band gap from f ≃ 0.8 to 0.84, and so on. Here, we focus on the photonic bands along the Г → X line in momentum space, which corresponds to the one-dimensional nature of the quasi-1D TD lattice laid along the x direction in real space. Indeed, one can see that there are two FBs at specific frequencies f ≃ 0.5 and 1.08, as shown in Fig. 2(a,b). These FBs, which correspond to the third and tenth bands, are found to be flat only along the Г → X line. Indeed, we achieve strongly localized electromagnetic waves in A and B, resulting in zero field intensity in C. (see lower pannels in Fig. 2(a,b)).
Meanwhile, the presence of the one-directional FB along Г → X line can be theoretically confirmed through tight-binding approach mimicking our PC. Complicated interactions of electromagnetic fields in the PC are well described with tight-binding method by taking into account higher-order hoppings. Beyond the nearest-neighbor hopping, our tight-binding approach to the FB is in good agreement with the numerically computed photonic band results. Leaving the details for the theoretical frameworks in Supplementary Information, let us emphasize that the FB along Г → X line can be secured even when higher-order hoppings are considered.
Since the FB reflects the existence of CLSs, propagation of light is expected to be prohibited at the FB frequencies in PCs. In order to examine the role of FBs for light propagation through the PC, we investigate the optical responses of the PC to a plane wave with TE polarization. We now consider a finite length of the PC in the x direction with a periodic boundary condition in the y direction. Twenty-one unit cells along the x direction are taken into account in the FDTD simulation, a number believed to be large enough to reflect the photonic band structure that we found in the infinite PC. We put a plane wave source sufficiently far away from one facet of the finite PC, and calculate the energy flux of the transmitted wave throughout the PC to yield optical transmittance as a ratio of energy fluxes between incident and transmitted waves. The rightmost and leftmost pannels of Fig. 1(b) show the transmittance spectra T Γ→X and T Γ→Y as functions of the incident plane wave frequency in the x and y directions, respectively. T Γ→X exhibit photonic band gaps in the specific frequency ranges where transmittance drops to zero. www.nature.com/scientificreports www.nature.com/scientificreports/ Overall, the photonic band gaps are well reflected in the transmittance spectra. However, we observe that certain photonic bands along the Г → X line seem to have no contribution to T Γ→X . As displayed by the shaded regions in Fig. 1(b), the third, seventh, eighth, and tenth bands seem to not be involved in T Γ→X . The reason for these inactive photonic bands will be discussed in the next section.
On the other hand, let us consider optical transmittance T Γ→Y through the 2D PC with a plane wave source propagating along the y direction. Because the unit cell considered in this study contains three sublattices that are anisotropic (i.e., the unit cell has no symmetry with respect to 90° rotation), the PC exhibits an anisotropic photonic band structure ( Fig. 1(b)). Overall, T Γ→Y is similar to T Γ→X , but a distinguishing feature is that the third band along Г → Y now contributes to optical transmittance. As highlighted in Fig. 2(b), the third band still seems to be flat over Г → Y. However, by investigating the related field intensities along the Г → Y line, it is revealed to be a dispersive band, as displayed in the lower panel of Fig. 2(c), so that the electromagnetic wave can be transferred from A and B to C, and vice versa.

Field Distributions and Compact Localized states
The optical property analysis results of the 2D PC in this study introduce two important questions: (i) how the FBs differ from the dispersive bands, and (ii) why certain photonic bands, including the FBs, have no contribution to the optical transmittance. To answer these questions, we analyze field distribution in the unit cell of the PC for the resulting photonic bands.
We first compare the electric field distributions of the flat and dispersive bands. As shown in Fig. 3, both dispersive and FBs exhibit a distinguishing distributions of field intensity |E| 2 . The second band, which is found to be dispersive in Fig. 1(b), indeed has a state that spreads over the unit cell, whereas the third band supports CLSs as |E| 2 is strongly localized in A and B nanopillars only. Such a localization of the electric field results from destructive interference in the C nanopillar as the electric fields in A and B have the opposite sign. On the other hand, the dispersive band comprises the same sign electric fields in A and B, causing constructive interference in C. It is worth mentioning that the dispersive and FBs can be characterized by even and odd parities in the y direction of the unit cell, respectively.
As aforementioned, two FBs were found in the given frequency range f = 0 to 1.1, corresponding to the third and tenth bands. Like the third band, the tenth band also exhibits a strong localization of electric field intensity. By comparing the electric fields in Fig. 3(d,e and g,h), one can easily see that the tenth FB is attributed to odd parity with respect to the y direction like the third band. Further, the tenth band occurs as a consequence of destructive interference between excited modes in A and B nanopillars, contrary to fundamental modes for the third band.
At this time, it should be noted that our electromagnetic wave source is provided as a plane wave propagating along the x direction, in other words, the plane wave source is characterized as even parity with respect to the y direction. This even-parity plane wave source does not excite the FBs, so that the FBs (the third and tenth bands) are not reflected in the transmittance spectra, as shown in Fig. 1. This is because the FBs are dark modes with respect to the plane-wave source. Similarly, the seventh and eighth bands are also dark, as highlighted by color-filled symbols in Fig. 4(a). By investigating their electric field patterns, we reveal that both the seventh and eighth bands indeed have odd parity with respect to the y direction (see Fig. 4). Here, let us notice that there are slight mismatches between the shaded regions in Fig. 4(a,b), which are because of finite-resolution errors. Additionally, let us discuss the seventh band for which even though the excited modes in A and B are out-of-phase seemingly to produce destructive interference in C, the band is dispersive and its field distribution is changed along the Г → X line. (see details in Supplementary Information).

self-collimation of Light propagation
Due to the anisotropy of the 2D PC, the transmittance spectra with an incident wave along the y direction are different from those for the x direction (see Fig. 1). In particular, it is worthwhile to note that T Γ→Y is nonzero at frequencies f ≃ 0.5 and 1.08, corresponding to the FBs along the Г → X line. In other words, we expect anisotropic light propagation throughout the 2D PC at the FB frequencies.
Now, we consider point sources inside the 2D PC instead of a plane wave source outside the PC. The point source radiates light isotropically, with a frequency precisely set at f = 0.501107, which is the frequency of the fundamental FB (third band). In order to examine the effects of parity, we place two monochromatic point sources at the centers of A and B nanopillars by changing their phase difference. Here, we stimulate a finite size (21 × 21 unit cells), and for simplicity, observe light propagation before the radiated waves are reflected from the boundaries of the PC. Figure 5 shows the field intensities with in-phase and out-of-phase point sources on A and B after some time steps in the FDTD simulations. As expected, for the out-of-phase point sources, the odd-parity waves in A and B nanopillars correspond to the FB shown in Fig. 2(a), resulting in prohibited light propagation in the x direction. Therefore, we observe an interesting optical phenomenon where the radiated light from the point sources is strongly collimated into the y direction, as exhibited in Fig. 5(a). Such a phenomenon is called "self-collimation"light propagation is allowed in only one direction because of vanishing group velocity in the perpendicular directions [39][40][41] . On the other hand, for the in-phase point sources, there is no photonic mode correlated with the even-parity waves in A and B, resulting in marginal effects in terms of light propagation (see very low field intensity in Fig. 5(b)). It is noteworthy that the collimated photon beam induced by the FB is expected to be extraordinary than the conventional self-collimation effects 42,43 , since the leakage into the perpendicular direction to the collimation is strongly suppressed by the zero group velocity of the FB mode.
As the existence of FBs has been believed to be a consequence of destructive interference, not only opposite phases but also equal amplitudes of the electromagnetic waves are therefore required. Despite this, the self-collimation effect remains observable even with a single point source, as displayed in Fig. 5(c). Such a robustness of the FB-induced collimation effect is understood as follows. The out-of-phase and in-phase excitations in use for Fig. 5(a,b) are represented as vector φ a = (1, −1)√2 and φ b = (1, 1)√2; the former corresponds to the FB mode whereas the latter does not. On the other hand, the single-source excitation in Fig. 5(c) is expressed as a vector composed of the linear combination of φ a and φ b , i.e., φ c = A 1 φ a + A 2 φ b . One can easily see that φ a is always dominant, since φ b barely affects light propagation in the PC. Consequently, the self-collimation effect in the FB PC is expected to be observable with arbitrary excitations in A and B.
Since the FB is formed by destructive interference for excited mode as well as fundamental modes, it is also expected to observe the self-collimation effect at the excited-mode-induced FB frequency, i.e., f = 1.08. Despite the existence of the FB mode at f = 1.08, the self-collimation of light propagation is not evidently observed, since the light sources excite not only the dispersive mode in the y direction but also dispersive modes along other directions (X → M and M → Γ). Details are shown in Supplementary Information. Lastly, let us discuss the fact that the operation frequency of the self-collimation effect should be finely tuned at the flat-band frequency (here, we set f = 0.501107). If one slightly increases/decreases the frequency of incident light, the self-collimation effect suddenly disappears since there is no photonic band just above/blow the FB. Such an extremely narrow bandwidth of the self-collimation effect enables our 2D PC to be promising for sensing slight changes in dielectric constant of backgrounds, leading to potential sensing applications like humidity sensors, gas sensors, DNA or polymer detectors. www.nature.com/scientificreports www.nature.com/scientificreports/

Conclusions
We have investigated the photonic band structure and optical properties of a 2D PCs, which is constructed by a periodic array of quasi-1D TD lattices. It has been shown that the 2D PC exhibits flat bands only along the Г → X line in the Brillouin zone, reflecting the existence of compact localized states as a consequence of destructive interference between the odd-parity waves in the A and B nanopillars. Thus, we can conclude that the FBs of the TD lattice remain even if the TD lattices are not isolated. We have also found that FBs occur not only for the fundamental modes but also for excited modes in the nanopillars. Further, we have learned that the FBs are not involved in optical transmittance through the PC because of their nondispersive nature. Moreover, by changing the direction of the incident plane wave, the anisotropy of the photonic bands and transmittance spectra were investigated. This anisotropic property of the 2D PC leads to an intriguing optical phenomenon-the self-collimation of light propagation, resulting from the existence of the nondispersive FB.
Flat band-induced strong localization acquired in disorder-free PCs can be widely employed to send photons over long distances with optical fiber communication technology. Also, our finding of a self-collimation effect at FB frequency may be practically applicable to photonics and metamaterials research. The strong collimation effect induced by the FB mode provides benefits to optical interferometry seeking highly focused light, and also the single-frequency operation of the FB-induced self-collimation effect may beget a possibility of sensing applications such as refractive-index sensors.