Flat-band localization and self-collimation of light in photonic crystals
Nojoon Myoung Hee Chul Park Ajith Ramachandran Elefterios Lidorikis Jung-Wan Ryu
physics.optics




    < g r a p h i c s >




We investigate the optical properties of a photonic crystal 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, but dispersive along the Γ→ Y line. The FB localization of light in a single direction only results in a self-collimation of light propagation throughout the photonic crystal at the FB frequency.


Photonic crystals (PCs), artificial materials governed by periodic dielectric functions, have been thoroughly studied for applications covering diverse spheres of interest. One key feature of PCs is the ability to manipulate light propagation through them—properties of electromagnetic waves inside PCs can be influenced by designing a periodic distribution of the dielectric function. Similar to electronic bands in solids, photonic bands exhibit band gaps where photonic states are prohibited, resulting in the prohibition of light propagation in a certain frequency range. Owing to the existence of photonic band gaps, PCs have been widely exploited as waveguides,<cit.> solar cells,<cit.> lasesr<cit.> etc. In particular, one intriguing aspect of PCs is a localization of light, where a primary approach is to introduce a defect or imperfection to a PC that leads to localized modes at the defect, regarded as analogous to microcavity modes.<cit.>

Light localization in PCs has been an interesting topic in terms of both fundamentals and related applications. In 2D PCs composed of optical cavities, for instance, Anderson localization of light can be experimentally confirmed by introducing a random disorder to the PC;<cit.> notably, the stronger the localization, the less diffractive the light propagation through the disordered PC.<cit.> In order to obtain stronger light localization in disordered PCs, a larger level of disorder is required.<cit.> Observation of such strong light localization, moreover, requires confirmation through statistical means, i.e., ensemble average over many disordered PCs. Thus, it is not feasible to make use of Anderson localization in device applications, despite its non-diffusive features in terms of photon transport.

Recently, a novel type of localization has begun to attract attention, the so-called compact localized states (CLSs) proposed by lattice models.<cit.> In specific quasi-1D lattices, such as cross-stitch, tunable diamond, etc., dispersion relations exhibit flat bands (FBs) over the whole Brillouin zone, corresponding to CLSs, or in other words a delta-function-like localization in real space.<cit.> Similarly, FBs and CLSs are found in specific 2D lattices like Lieb and Kagome, among others.<cit.> The existence of CLSs has been experimentally demonstrated in PCs,<cit.> polaritonic systems,<cit.> and mechanical lattices.<cit.> Contrary to Anderson localization, there is no need for statistical tasks to demonstrate CLSs because they require no disorder. Indeed, it has been experimentally shown that flat-band modes are non-diffractive, as diffraction is prohibited via destructive interference.<cit.> However, since flat-band modes are a consequence of destructive interference between amplitudes residing on specific sites in PCs, the corresponding CLSs are expected to be easily perturbed by inhomogeneity, and thus require sophisticated fabrication and measurement techniques.

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 localization of light corresponding to the 1D flat-band mode can remain compact even in the presence of coupling between individual quasi-1D lattices. The CLSs are found to be odd symmetric with respect to the transverse direction of the lattice. The photonic band structure of the 2D PC exhibits an anisotropic nature, following anisotropic couplings between lattice sites. Notably, for a specific frequency corresponding to the CLSs, we reveal a self-collimation of light propagation with a given out-of-phase excitation, thanks to strong localization in one direction only. Understanding defect-free localization of light in PCs paves the way for enhancing performance of optical devices, telecommunications, and sensing applications.



§ RESULTS AND DISCUSSION






 §.§ 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. <ref>(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. <ref>(b). Note that coupling between A (or B) and C differs from that between A and B. All dielectric nanopillars are identical with radius r, and the period of the TD lattice is a. In our 2D PC, this quasi-1D TD lattice is repeated with the period a, forming a square unit cell with three sublattices.





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<cit.>. Figure <ref>(a) shows the photonic band structure for TE polarization (E_x=E_y=0 but E_z≠0) of the PC calculated in the 2D Brillouin zone. 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. <ref>(a) and (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. <ref>(a) and (b))  

Since the FB reflects the existence of CLSs, propagation of light is expected to be prohibited at the flat-band 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. <ref> 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. 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. <ref>, 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 subsection.

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. <ref>) where the symmetric point Y replaces X in Fig. <ref>(a). 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. <ref>(c), 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. <ref>(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 flat bands 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. <ref>,  both dispersive and FBs exhibit a distinguishing distributions of field intensity |E|^2. The second band, which is found to be dispersive in Fig. <ref>(a), 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. <ref>(d) and (e), and (g) and (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 the 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. <ref>. 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-filed symbols in Fig. <ref>(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. <ref>). Here, let us notice that there are slight mismatches between the shaded regions in Fig. <ref>(a) and (b), which are because of finite-resolution errors.



 §.§ 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. <ref>). 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 flat-band frequencies.

Now, we consider  point sources inside the 2D PC instead of a plane wave source outside the PC. The point source radiates isotropic light 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 <ref> 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. <ref>(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. <ref>(a). Such a phenomenon is called “self-collimation”—light propagation is allowed in one only direction because of vanishing group velocity in the perpendicular directions.<cit.> 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. <ref>(b)).

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 effects still remain observable even with a single point source, as displayed in Fig. <ref>(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 Figs. <ref>(a) and (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. <ref>(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.



§ CONCLUSIONS


We have investigated the photonic band structure and optical properties of a two-dimensional photonic crystals, which is constructed by a periodic array of quasi-1D tunable diamond 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 photonic and metamaterial research. 


This work was supported by Project IBS-R024-D1, the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (No. 2017076824), and research funds from Chosun University 2017.




