Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Resonance modes in moiré photonic patterns for twistoptics

Open Access Open Access

Abstract

Twistronics has been studied for manipulating electronic properties through a twist angle in the formed moiré superlattices of two dimensional layer materials. In this paper, we study twistoptics for manipulating optical properties in twisted moiré photonic patterns without physical rotations. We describe a theoretic approach for the formation of single-layer twisted photonic pattern in square and triangular lattices through an interference of two sets of laser beams arranged in two cone geometries. The moiré period and the size of unit super-cell of moiré patterns are related to the twist angle that is calculated from the wavevector ratio of laser beams. The bright and dark regions in moiré photonic pattern in triangular lattices are reversible. We simulate E-field intensities and their cavity quality factors for resonance modes in moiré photonic pattern in square lattices. Due to the bandgap dislocation between the bright and dark regions, the resonance modes with very high quality-factors appears near bandgap edges for the moiré photonic pattern with a twist angle of 9.5 degrees. At the low frequency range, the resonance modes can be explained as Mie resonances. The cavity quality factor decreases for resonance modes when the twist angle is increased to 22.6 degrees.

© 2021 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

The moiré phenomenon occurs when repetitive structures are superposed against each other with a relative rotation angle [1]. Recently, twisted bilayers of two-dimensional (2D) materials have ignited significant interest due to the ability in experiments to control the relative twist angle between successive layers [259]. This approach was referred as “twistronics” [2,32,57] for manipulating their electronic properties through a twist angle in the formed moiré superlattices. Manipulating the twist angle between the two layers of 2D materials has resulted in a fine control of the electronic band structure [2], magic-angle flat-band superconductivity [37], the formation of moiré excitons [812], interlayer magnetism [13,39,58], topological edge states [8,1418], and correlated insulator phases [2022]. Research has progressed to moiré superlattices between two layers of van der Waals materials [10,13,18] including twisted bilayer transition metal dichalcogenides [4850].

Very recently, several groups have studied twistoptics: studies of optical properties in twisted moiré lattice [6062]. A moiré photonic lattice has been generated in a photorefractive crystal by a shallow modulation of the refractive index induced by two mutually twisted, periodic square lattices generated by laser interference (both sublattices interfere in one plane) [61,62]. The study of twistoptics is still in the very early stage [6062].

In this paper, we describe an approach for the formation of single-layer moiré photon pattern through an interference of two sets of interfering beams. The twist angle of the moiré interference pattern is determined by the wavevector ratio of these laser beams. We simulate the cavity quality factors for resonance modes in a wide frequency range for the moiré photonic pattern in square lattices with a twist angle of 9.5 degrees and in a certain frequency range for the one with a twist angle of 22.6 degrees. The high cavity quality factors are observed for resonance modes near photonic bandgap edges.

2. Theoretical approach for single-layer moiré photonic pattern for twistoptics

We introduce the moiré concept into the laser interference pattern in a square lattice. When laser beams (with wave-vectors of km and kn) interfere, the intensity of the interference pattern is modulated by a cosine function with a period and orientation that can be characterized by a wave-vector difference of (km-kn). Figure 1(a) shows an arrangement of wave-vectors (k1, k2, …, k8) for eight laser beams. (k1,k2,k3,k4) have the same magnitude and (k5,k6,k7,k8) also have the same magnitude that is much smaller than that for (k1,k2,k3,k4). The wave-vector difference (kn+1-kn) among (k1,k2,k3,k4) determines square photonic lattices with a small periodicity. As shown in Fig. 1(a), the direction of wavevector difference (k1-k8) is twisted from (k3-k8) by a rotation angle α. Thus the formed square lattices by (k1-k6), (k2-k7), (k3-k8), and (k4-k5) is twisted by a rotation angle α from those by (k1-k8), (k2-k5), (k3-k6), and (k4-k7).

 figure: Fig. 1.

Fig. 1. (a) Schematic of arrangement of 8 wave-vectors of interfering beams for the generation of a twist angle α; (b) Schematic of optical setup for 8 interfering laser beams; (c-d) moiré interference patterns with a twist angle of 7.2 degrees at different z-location=0 and 0.27Pz, respectively. The unit super-cell, moiré period Pmoiré, lattice period Psmall, bright and dark regions are indicated in figures. (e) x-z cross section of moiré interference pattern.

Download Full Size | PDF

In experimental setup, these wave-vectors (k1, k2, …, k8) can be achieved through the configuration of interfering angles [63,64]. As shown in Fig. 1(b), one set of interfering beams are arranged in a cone with a cone angle of 2β1 with their wave vectors of: {k1…4}={k(sin(β1)cos(45+n×90)°, sin(β1)sin(45+n×90)°, cos(β1)), n=0,1,2,3. The other set of interfering beams are arranged in a cone with a much smaller cone angle of 2β2 than 2β1 with wave vectors of {k5…8}={k(sin(β2)cos(45+n×90)°, sin(β2)sin(45+n×90)°, cos(β2)), n=0,1,2,3. In this case, the twist angle α in Fig. 1(a) is determined by Eq. (1):

$$\mathrm{\alpha } = 2\arctan \frac{{\textrm{sin}({{\beta_2}} )}}{{\textrm{sin}({\; {\beta_1}} )}}$$

The interference intensity of these 8 beams, I(r), is calculated by Eq. (2):

$$I(r )= \left\langle{\mathop \sum \limits_{i = 1}^\textrm{N} E_i^2({r,t} )}\right\rangle+ \mathop \sum \limits_{i < j}^\textrm{N} {E_i}{E_j}{e_i} \cdot {e_j}cos[{({{k_j} - {k_i}} )\cdot r + ({{\delta_j} - {\delta_i}} )} ].\; $$
where e is the electric field polarization, k is the wave vector, E is the electric field, and δ is the initial phase. The interference occurs among all beams. As shown in Fig. 2(c), the small period Psmall equals approximately 2π/(2ksin(β1)sin45). The moiré period Pmoiré can be approximately estimated from the cone angle of 2β2 and twist angle α by Eq. (3):
$${P_{moiré}} \cong \frac{{2\pi }}{{2ksin({{\beta_2}} )sin45}} = \frac{{2\pi }}{{2k \times \tan \left( {\frac{\mathrm{\alpha }}{2}} \right)\sin ({{\beta_1}} )sin45}} = \frac{{{P_{small}}}}{{\tan \left( {\frac{\mathrm{\alpha }}{2}} \right)}}$$

For bilayer graphene with a twist angle of θ, the moiré period equals a/[2sin(θ/2)], where a is the lattice constant of graphene [20]. From both Eq. (3) and the equation in above line for the moiré period for the twisted bilayer graphene, the moiré period and the size of unit super-cell increase with decreasing twist angle (rotational angle).

 figure: Fig. 2.

Fig. 2. (a) Schematic of the arrangement of 12 wave-vectors of interfering beams for the generation of a twist angle α indicted by dashed blue line and dotted purple line; (b-c) moiré interference patterns with a twist angle of 9.9 degrees at z-location=0 and 0.55Pz, respectively. The moiré period Pmoiré, and bright regions are indicated. Dashed red circles indicate the interface between two distinct pattern regions.

Download Full Size | PDF

The moiré period Pmoiré, small period Psmall (lattice period) and unit super-cell are indicated in Fig. 1(c) for a moiré pattern with a twist angle of 7.2 degrees. In Eq. (2), there is an interference intensity modulation in z-direction by cos(k(cos(β1)-cos(β2))z). The formed moiré photonic lattice has a feature distance between bright zones of Pz=λ/(cos(β1)-cos(β2)) in z-direction as indicated in Fig. 1(e). At location z=0 as indicated by the solid white line in Fig. 1(e), the pattern in Fig. 1(c) has a bright region at the center of the solid red square and dark region at the corners of the square, corresponding to AA and AB stacking region in twisted bilayer 2D layered materials [26], respectively. Figure 1(d) shows a moiré interference pattern at a location z=0.27Pz as indicated by the dashed white line in Fig. 1(e). Dual lattices appear at the central region inside the unit super-cell in the figure. The bright and dark regions have been described as almost-uniform and graded regions and the photonic patterns have been assigned as graded photonic super-crystals [63,64]. A pattern at a different z-location can be selected for the study of twistoptics. As shown in Fig. 1(e), we can get the clearest moiré pattern at z=0.

We can also generate moiré photonic patterns in triangular lattice. As shown in Fig. 2(a), one set of interfering beams have six beams and are arranged in a cone with a cone angle of 2β1 with their wave vectors of: {k1…6}={k(sin(β1)cos(n×60)°, sin(β1)sin(n×60)°, cos(β1)), n=0,1,..,5. The other six interfering beams are arranged in a cone with a much smaller cone angle of 2β2 than 2β1 with wave vectors of {k7…12}={k(sin(β2)cos(n×60)°, sin(β2)sin(n×60)°, cos(β2)), n=0,1,…5. The wave-vector difference (kn+1-kn) among (k1,…,k6) determines triangular photonic lattices with a small periodicity Psmall ${\cong} $ 2π/(2ksin(β1)sin60). As indicated by dashed blue lines in Fig. 2(a), the direction of vector difference (k1-k9) is twisted from (k4-k9) by a rotation angle α. The direction of vector difference (k1-k8) is also twisted from (k4-k8) by the same rotation angle α as indicated by dashed purple lines. Thus the formed triangular lattices by (k1-k9), (k2-k10), (k3-k11), (k4-k12), (k5-k7) and (k6-k8) is twisted by a rotation angle α from those by (k4-k9), (k5-k10), (k6-k11), (k1-k12) and (k2-k7). Similarly, the formed triangular lattices by (k1-k8), (k2-k9), (k3-k10), (k4-k11), (k5-k12) and (k6-k7) is twisted by the same rotation angle α from those by (k4-k8), (k5-k9), (k6-k10), (k1-k11) and (k2-k12). α=a1+a2 as shown in Fig. 2(a). These angles can be determined by following equations:

$${\alpha _1} = \arctan \frac{{\textrm{ksin}({{\beta_2}} )\sin ({60} )}}{{\textrm{ksin}({\; {\beta_1}} )- \textrm{ksin}({{\beta_2}} )\cos ({60} )}},\; $$
$$\; {\alpha _2} = \arctan \frac{{\textrm{ksin}({{\beta_2}} )\sin ({60} )}}{{\textrm{ksin}({\; {\beta_1}} )+ \textrm{ksin}({{\beta_2}} )\cos ({60} )}},$$
$$\alpha = {\alpha _1} + {\alpha _2} \cong 2\arctan \frac{{\textrm{ksin}({{\beta_2}} )\sin ({60} )}}{{\textrm{ksin}({\; {\beta_1}} )}},\; $$
$$\tan \left( {\frac{\alpha }{2}} \right) = \frac{{{p_{small}}}}{{{p_{moiré}}}}\sin ({60} ),$$

Figure 2(b) and 2(c) show moiré interference patterns with a twist angle of 9.9 degrees at z-location=0 and 0.55Pz, respectively. These two patterns are different due to the modulation from k(cos(β1)-cos(β2)z in the cosine function in Eq. (2). The pattern inside the dashed red circle in Fig. 2(b) has motif of cylinders in a triangular lattice orientated in directions of 0, 60 and 120 degrees. Outside the dashed red circle, the motifs are mostly in triangular shape in a triangular lattice orientated in 30, 90, and 150 degrees. In Fig. 2(c), the triangular lattices inside and outside the dashed red circle have a same orientation. There are air ring structures surrounding the cylinders inside the dashed red circle while there are only cylinder patterns outside. Thus there are two distinct photonic materials inside and outside the dashed red circles in Fig. 2(b) and 2(c). After half Pz, the bright and dark regions are switched. The bright regions are inside the dashed red circles in Fig. 2(b) while the bright regions are outside in Fig. 2(c).

3. Cavity quality factors and resonance modes in moiré photonic lattices with different twisted angles

In this section, we simulate cavity quality factors and electric field (E-field) intensity distributions for resonance modes in silicon moiré photonic pattern using the harmonic inversion function included in the MIT MEEP software [65,66]. The interference intensity function, I(r), in Eq. (2), is replaced by a binary Si/air structures by comparing I(r) with a threshold intensity, Ith, that is set to be 28% of maximum of I(r). A step function is used: ɛ(r) = 1 (for air) when I < Ith, and ɛ(r) = 12 (for Si) when I > Ith. Simulated E-field intensity distributions as shown in Fig. 3(a-l), is performed for the moiré photonic pattern with a twist angle of 9.5 degrees, a thickness of 2Psmall, and a unit cell size of 12Psmall×12Psmall, the same size as the one used for photonic band gap calculation [63,67]. The dashed red squares indicate the unit super-cell. The bright (graded) and dark (almost-uniform) regions are indicated by the dashed blue and green squares, respectively. The frequency Psmall/λ and Q-factor are labelled for each figure. At frequency Psmall/λ= (0.576, 0.531, 0.503, 0.491, 0.478, 0.432, 0.380, 0.343, 0.333, 0.267, 0.255, and 0.238), we obtain Q-factor=(291256, 402479, 32690, 863601, 521476, 254402, 417278, 13306, 1491234, 69204, 26000, and 167471), respectively.

 figure: Fig. 3.

Fig. 3. Simulated electric field intensity for resonance modes in moiré photonic lattice with a twist angle of 9.5°. Frequency Psmall/λ and Q-factor are labelled in each figure (a-l).

Download Full Size | PDF

Figure 4 plots Q-factors at different frequency for resonance modes in moiré Si photonic pattern in Fig. 3 with a twist angle of 9.5 degrees. The red rectangle in Fig. 4 indicates the location of photonic band gap for the bright (graded) region in Fig. 3 while the purple rectangle indicates the common photonic band gap in both bright (graded) and dark (almost-uniform) regions [63,67]. Within the frequency range indicated by these rectangles, there is no strong resonance mode in Fig. 4. At the top edge of band gap at Psmall/λ = 0.478, and bottom edge at 0.432, the resonance modes show up just outside the dashed blue squares in Fig. 3(e) and 3(f). A resonance mode in Fig. 3(g) appears just inside the blue square at Psmall/λ = 0.478, near the edge of photonic band gap in dark region. In Fig. 3(h), the beam propagates in the dark region and has weak resonance. The strong modes in Fig. 3(i) can be considered as a resonance at a point-cavity. At higher frequencies, resonance modes cover both bright (graded) and dark (almost-uniform) regions in Fig. 3(a-d). The resonance modes look like square shape orientated 45° and 0° relative to the horizontal direction. In Fig. 3(c), the resonance modes separate and appear in both bright and dark regions, however, the oscillation direction is rotated by 45 degrees in Fig. 3(d).

 figure: Fig. 4.

Fig. 4. Calculated Q-factors at different frequencies for resonance modes in Si moiré photonic pattern in Fig. 3 with a twist angle of 9.5 degrees. The red and purple rectangles indicate the locations of simulated photonic band gaps.

Download Full Size | PDF

Resonances in Fig. 3(j-l) occur at Psmall/λ between 0.267 and 0.238. Assuming the cylinder diameter D in the moiré Si photonic pattern is half of Psmall, we can calculate:

$$\frac{{\pi D\; \sqrt {12} }}{\mathrm{\lambda }} = 1.45\sim 1.29,\; \textrm{for}\; \frac{{{P_{small}}}}{\mathrm{\lambda }} = 0.267\sim 0.238\; ,\; $$
$$\frac{{2D\; \sqrt {12} }}{\mathrm{\lambda }} = 0.82\sim 0.92,\; \textrm{for}\; \frac{{{P_{small}}}}{\mathrm{\lambda }} = 0.267\sim 0.238\; ,$$

From Eq. (8) and (9), the ratio of the size of Si cylinder over the wavelength is near 1, thus these resonances in Fig. 3(j-l) can be assigned to Mie resonances [68].

Figure 5 shows simulated E-field intensity distributions for the Si moiré photonic pattern with a twist angle of 22.6 degrees and a thickness of 2Psmall. The dashed blue squares indicate the unit super-cell. The frequency Psmall/λ and Q-factor are labelled in each figure. At frequency Psmall/λ= (0.333, 0.343, and 0.355), we obtain Q-factor=(1273, 8131 and 37877), respectively. These Q-factors in Fig. 5 is much smaller than these in Fig. 3(h-i). It is due to the coupling of resonance modes in the photonic pattern with a small unit super-cell, i.e. 5Psmall×5Psmall in Fig. 5. We have not performed the simulation of Q-factors in other frequency ranges in Fig. 3. Further simulations are needed in order to obtain full knowledge of twistoptics. From Fig. 5(a), the resonance modes at the corners of the blue square (dark regions) oscillate in a direction 45 degrees relative to the horizontal, couples to weaker modes in the top and bottom dark regions with a vertical oscillations, and to weaker modes in the left and right dark regions with a horizontal oscillations. The resonance modes in Fig. 5(b) have already appeared in the bright regions outside the blue square during the simulation decay times. At a frequency of Psmall/λ=0.355 in Fig. 5(c), the mode is localized inside the blue square with a weaker coupling to bright regions outside and has a higher Q-factor compared with these modes in Fig. 5(b).

 figure: Fig. 5.

Fig. 5. Simulated electric field intensity for resonance modes at the dark regions (a), bright regions (b) and highly localized mode at the bright region (c) in Si moiré photonic pattern with a twist angle of 22.6 degrees. The unit super-cell is indicated by dashed blue squares.

Download Full Size | PDF

4. Discussion

The moiré photonic structure can be fabricated by either E-beam lithography or holographic lithography on a substrate of silicon-on-insulator. When E-beam lithography is used, the laser interference pattern from Eq. (2) is used as an input for the E-beam pattern writing program.

So far we picked up symmetric patterns with a desired moiré period Pmoiré. For Pmoiré=16, 12 and 5Psmall, the corresponding pattern has a twist angle of 7.2, 9.5 and 22.6 degrees, respectively. Further studies can include these asymmetric patterns generated through any twist angles. There are many moiré photonic patterns after combining different twist angles and z-locations. To address this challenge of huge work load, the machine learning method [69] can be applied to perform the research work on twistoptics.

5. Conclusion

We have presented an approach for the generation of moiré photonic patterns with different twist angles through interference of two sets of interfering laser beams, without physical rotations. The moiré photonic patterns have been demonstrated in square and triangular lattices. Resonance modes and cavity Q-factors have been simulated in these Si moiré photonic patterns with twist angles of 9.5 and 22.6 degrees. High Q-factors have appeared for these resonance modes near the edge of the photonic band gaps. The presented approach can lead toward a study of twistoptics for the manipulating of optical properties in moiré patterns with different twist angles.

Funding

National Science Foundation (1661842).

Acknowledgments

This work is supported by research grants from the U.S. National Science Foundation under Grant No. CMMI-1661842.

Disclosures

The authors declare no conflict of interest related to this paper.

References

1. I. Amidror, The Theory of the Moiré Phenomenon: Volume I: Periodic Layers, 2nd ed. (Springer2009).

2. S. Carr, D. Massatt, S. Fang, P. Cazeaux, M. Luskin, and E. Kaxiras, “Twistronics: Manipulating the electronic properties of two-dimensional layered structures through their twist angle,” Phys. Rev. B 95(7), 075420 (2017). [CrossRef]  

3. Y. Cao, V. Fatemi, S. Fang, K. Watanabe, T. Taniguchi, E. Kaxiras, and P. Jarillo-Herrero, “Unconventional superconductivity in magic-angle graphene superlattices,” Nature 556(7699), 43–50 (2018). [CrossRef]  

4. M. Yankowitz, S. Chen, H. Polshyn, Y. Zhang, K. Watanabe, T. Taniguchi, D. Graf, A. F. Young, and C. R. Dean, “Tuning superconductivity in twisted bilayer graphene,” Science 363(6431), 1059–1064 (2019). [CrossRef]  

5. H. S. Arora, R. Polski, Y. Zhang, A. Thomson, Y. Choi, H. Kim, Z. Lin, I. Z. Wilson, X. Xu, J.-H. Chun, K. Watanabe, T. Taniguchi, J. Alicea, and S. Nadj-Perge, “Superconductivity without insulating states in twisted bilayer graphene stabilized by monolayer WSe2,” Nature 583(7816), 379–384 (2020). [CrossRef]  

6. L. Balents, C. R. Dean, D. K. Efetov, and A. F. Young, “Superconductivity and strong correlations in moiré flat bands,” Nat. Phys. 16(7), 725–733 (2020). [CrossRef]  

7. J. Gonzalez and T. Stauber, “Kohn-luttinger superconductivity in twisted bilayer graphene,” Phys. Rev. Lett. 122(2), 026801 (2019). [CrossRef]  

8. F. Wu, T. Lovorn, and A. H. MacDonald, “Topological exciton bands in moiré heterojunctions,” Phys. Rev. Lett. 118(14), 147401 (2017). [CrossRef]  

9. H. Yu, G.-B. Liu, J. Tang, X. Xu, and W. Yao, “Moiré excitons: From programmable quantum emitter arrays to spin-orbit–coupled artificial lattices,” Sci. Adv. 3(11), e1701696 (2017). [CrossRef]  

10. Kha Tran, Galan Moody, Fengcheng Wu, Xiaobo Lu, Junho Choi, Kyounghwan Kim, Amritesh Rai, Daniel A. Sanchez, Jiamin Quan, Akshay Singh, Jacob Embley, André Zepeda, Marshall Campbell, Travis Autry, Takashi Taniguchi, Kenji Watanabe, Nanshu Lu, Sanjay K. Banerjee, Kevin L. Silverman, Suenne Kim, Emanuel Tutuc, Li Yang, Allan H. Macdonald, and Xiaoqin Li, “Evidence for moiré excitons in van der Waals heterostructures,” Nature 567(7746), 71–75 (2019). [CrossRef]  

11. X. Lu, X. Li, and L. Yang, “Modulated interlayer exciton properties in a two-dimensional moiré crystal,” Phys. Rev. B 100(15), 155416 (2019). [CrossRef]  

12. Y. Bai, L. Zhou, J. Wang, W. Wu, L. Mcgilly, D. Halbertal, C. Fan, B. Lo, F. Liu, J. Ardelean, P. Rivera, N. R. Finney, X. Yang, D. N. Basov, W. Yao, X. Xu, J. Hone, A. Pasupathy, and X. Zhu, “Excitons in strain-induced one-dimensional moiré potentials at transition metal dichalcogenide heterojunctions,” Nat. Mater. 19(10), 1068–1073 (2020). [CrossRef]  

13. W. Chen, Z. Sun, Z. Wang, L. Gu, X. Xu, S. Wu, and C. Gao, “Direct observation of van der Waals stacking–dependent interlayer magnetism,” Science 366(6468), 983–987 (2019). [CrossRef]  

14. S. Huang, K. Kim, D. K. Efimkin, T. Lovorn, T. Taniguchi, K. Watanabe, A. H. MacDonald, E. Tutuc, and B. J. LeRoy, “Topologically protected helical states in minimally twisted bilayer graphene,” Phys. Rev. Lett. 121(3), 037702 (2018). [CrossRef]  

15. F. Wu, T. Lovorn, E. Tutuc, I. Martin, and A. MacDonald, “Topological insulators in twisted transition metal dichalcogenide homobilayers,” Phys. Rev. Lett. 122(8), 086402 (2019). [CrossRef]  

16. Shaowen Chen, Minhao He, Ya-Hui Zhang, Valerie Hsieh, Zaiyao Fei, K. Watanabe, T. Taniguchi, David H. Cobden, Xiaodong Xu, Cory R. Dean, and Matthew Yankowitz, “Electrically tunable correlated and topological states in twisted monolayer-bilayer graphene,” Nat. Phys. 17, 374–380 (2021). [CrossRef]  

17. G. Hu, Q. Ou, G. Si, Y. Wu, J. Wu, Z. Dai, A. Krasnok, Y. Mazor, Q. Zhang, Q. Bao, C.-W. Qiu, and A. Alù, “Topological polaritons and photonic magic angles in twisted α-MoO3 bilayers,” Nature 582(7811), 209–213 (2020). [CrossRef]  

18. Q. Tong, H. Yu, Q. Zhu, Y. Wang, X. Xu, and W. Yao, “Topological mosaics in moiré superlattices of van der Waals heterobilayers,” Nat. Phys. 13(4), 356–362 (2017). [CrossRef]  

19. S. S. Sunku, G. X. Ni, B. Y. Jiang, H. Yoo, A. Sternbach, A. S. McLeod, T. Stauber, L. Xiong, T. Taniguchi, K. Watanabe, P. Kim, M. M. Fogler, and D. N. Basov, “Photonic crystals for nano-light in moiré graphene superlattices,” Science 362(6419), 1153–1156 (2018). [CrossRef]  

20. Y. Cao, V. Fatemi, A. Demir, S. Fang, S. L. Tomarken, J. Y. Luo, J. D. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, and P. Jarillo-Herrero, “Correlated insulator behaviour at half-filling in magic-angle graphene superlattices,” Nature 556(7699), 80–84 (2018). [CrossRef]  

21. G. W. Burg, J. Zhu, T. Taniguchi, K. Watanabe, A. H. MacDonald, and E. Tutuc, “Correlated insulating states in twisted double bilayer graphene,” Phys. Rev. Lett. 123(19), 197702 (2019). [CrossRef]  

22. Y. Cao, J. Luo, V. Fatemi, S. Fang, J. Sanchez-Yamagishi, K. Watanabe, T. Taniguchi, E. Kaxiras, and P. Jarillo-Herrero, “Superlattice-induced insulating states and valley-protected orbits in twisted bilayer graphene,” Phys. Rev. Lett. 117(11), 116804 (2016). [CrossRef]  

23. E. Koren, I. Leven, E. Lörtscher, A. Knoll, O. Hod, and U. Duerig, “Coherent commensurate electronic states at the interface between misoriented graphene layers,” Nat. Nanotechnol. 11(9), 752–757 (2016). [CrossRef]  

24. A. V. Rozhkov, A. Sboychakov, A. Rakhmanov, and F. Nori, “Electronic properties of graphene-based bilayer systems,” Phys. Rep. 648, 1–104 (2016). [CrossRef]  

25. K. Uchida, S. Furuya, J.-I. Iwata, and A. Oshiyama, “Atomic corrugation and electron localization due to moiré patterns in twisted bilayer graphenes,” Phys. Rev. B 90(15), 155451 (2014). [CrossRef]  

26. C. R. Woods, L. Britnell, A. Eckmann, R. S. Ma, J. C. Lu, H. M. Guo, X. Lin, G. L. Yu, Y. Cao, R. V. Gorbachev, A. V. Kretinin, J. Park, L. A. Ponomarenko, M. I. Katsnelson, Y. N. Gornostyrev, K. Watanabe, T. Taniguchi, C. Casiraghi, H.-J. Gao, A. K. Geim, and K. S. Novoselov, “Commensurate–incommensurate transition in graphene on hexagonal boron nitride,” Nat. Phys. 10(6), 451–456 (2014). [CrossRef]  

27. A. Ebnonnasir, B. Narayanan, S. Kodambaka, and C. V. Ciobanu, “Tunable MoS2 bandgap in MoS2-graphene heterostructures,” Appl. Phys. Lett. 105(3), 031603 (2014). [CrossRef]  

28. Z. Zhu, S. Carr, D. Massatt, M. Luskin, and E. Kaxiras, “Twisted Trilayer Graphene: A Precisely Tunable Platform for Correlated Electrons,” Phys. Rev. Lett. 125(11), 116404 (2020). [CrossRef]  

29. S. Carr, S. Fang, and E. Kaxiras, “Electronic-structure methods for twisted moiré layers,” Nat. Rev. Mater. 5(10), 748–763 (2020). [CrossRef]  

30. G. A. Tritsaris, S. Carr, Z. Zhu, Y. Xie, S. B. Torrisi, J. Tang, M. Mattheakis, D. Larson, and E. Kaxiras, “Electronic structure calculations of twisted multi-layer graphene superlattices,” 2D Mater. 7(3), 035028 (2020). [CrossRef]  

31. D. Waters, Y. Nie, F. Lüpke, Y. Pan, S. Fölsch, Y.-C. Lin, B. Jariwala, K. Zhang, C. Wang, H. Lv, D. X. Kyeongjae Cho, J. A. Robinson, and R. M. Feenstra, “Flat bands and mechanical deformation effects in the moiré superlattice of MoS2-WSe2 heterobilayers,” ACS Nano 14(6), 7564–7573 (2020). [CrossRef]  

32. P. Kang, “Indirect-to-direct bandgap transition in bilayer inse: roles of twistronics,” 2D Mater. 7(2), 021002 (2020). [CrossRef]  

33. F. J. Culchac, R. Del Grande, R. B. Capaz, L. Chico, and E. S. Morell, “Flat bands and gaps in twisted double bilayer graphene,” Nanoscale 12(8), 5014–5020 (2020). [CrossRef]  

34. Z. Zhang, Y. Wang, K. Watanabe, T. Taniguchi, K. Ueno, E. Tutuc, and B. J. LeRoy, “Flat bands in twisted bilayer transition metal dichalco- genides,” Nat. Phys. 16, 1093–1096 (2020). [CrossRef]  

35. A. Kerelsky, L. J. McGilly, D. M. Kennes, L. Xian, M. Yankowitz, S. Chen, K. Watanabe, T. Taniguchi, J. Hone, C. Dean, A. Rubio, and A. N. Pasupathy, “Maximized electron interactions at the magic angle in twisted bilayer graphene,” Nature 572(7767), 95–100 (2019). [CrossRef]  

36. F. Hu, S. R. Das, Y. Luan, T.-F. Chung, Y. Chen, and Z. Fei, “Real-space imaging of the tailored plasmons in twisted bilayer graphene,” Phys. Rev. Lett. 119(24), 247402 (2017). [CrossRef]  

37. S. Carr, S. Fang, P. Jarillo-Herrero, and E. Kaxiras, “Pressure depen- dence of the magic twist angle in graphene superlattices,” Phys. Rev. B 98(8), 085144 (2018). [CrossRef]  

38. Y. Luo, R. Engelke, M. Mattheakis, M. Tamagnone, S. Carr, K. Watan- abe, T. Taniguchi, E. Kaxiras, P. Kim, and W. L. Wilson, “In situ nanoscale imaging of moiré superlattices in twisted van der Waals heterostructures,” Nat. Commun. 11(1), 4209 (2020). [CrossRef]  

39. D. A. Bahamon, G. Gómez-Santos, and T. Stauber, “Emergent magnetic texture in driven twisted bilayer graphene,” Nanoscale 12, 15383–15392 (2020). [CrossRef]  

40. G. Abbas, Y. Li, H. Wang, W.-X. Zhang, C. Wang, and H. Zhang, “Recent advances in twisted structures of flatland materials and crafting moiré superlattices,” Adv. Funct. Mater. 30(36), 2000878 (2020). [CrossRef]  

41. Z. Zheng, F. Sun, W. Huang, J. Jiang, R. Zhan, Y. Ke, H. Chen, and S. Deng, “Phonon polaritons in twisted double-layers of hyperbolic van der waals crystals,” Nano Lett. 20(7), 5301–5308 (2020). [CrossRef]  

42. G. Hu, A. Krasnok, Y. Mazor, C.-W. Qiu, and A. Alù, “Moiré hyperbolic metasurfaces,” Nano Lett. 20(5), 3217–3224 (2020). [CrossRef]  

43. P. Sutter, R. Ibragimova, H.-P. Komsa, B. A. Parkinson, and E. Sutter, “Self-organized twist-heterostructures via aligned van der waals epitaxy and solid-state transformations,” Nat. Commun. 10(1), 5528 (2019). [CrossRef]  

44. Y. Jiang, X. Lai, K. Watanabe, T. Taniguchi, K. Haule, J. Mao, and E. Y. Andrei, “Charge order and broken rotational symmetry in magic-angle twisted bilayer graphene,” Nature 573(7772), 91–95 (2019). [CrossRef]  

45. P. Sutter, S. Wimer, and E. Sutter, “Chiral twisted van der Waals nanowires,” Nature 570(7761), 354–357 (2019). [CrossRef]  

46. F. Guinea and N. R. Walet, “Continuum models for twisted bilayer graphene: effect of lattice deformation and hopping parameters,” Phys. Rev. B 99(20), 205134 (2019). [CrossRef]  

47. L. Wang, S. Zihlmann, M.-H. Liu, P. Makk, K. Watanabe, T. Taniguchi, A. Baumgartner, and C. Schonenberger, “New generation of moiré superlattices in doubly aligned hBN/graphene/hBN heterostructures,” Nano Lett. 19(4), 2371–2376 (2019). [CrossRef]  

48. Y. Li and M. Koshino, “Twist-angle dependence of the proximity spin- orbit coupling in graphene on transition-metal dichalcogenides,” Phys. Rev. B 99(7), 075438 (2019). [CrossRef]  

49. M. H. Naik and M. Jain, “Ultraflatbands and shear solitons in moiré patterns of twisted bilayer transition metal dichalcogenides,” Phys. Rev. Lett. 121(26), 266401 (2018). [CrossRef]  

50. S. Zhu and H. T. Johnson, “Moiré-templated strain patterning in transition-metal dichalcogenides and application in twisted bilayer MoS2,” Nanoscale 10(44), 20689–20701 (2018). [CrossRef]  

51. T. Stauber, T. Low, and G. Gómez-Santos, “Linear response of twisted bilayer graphene: Continuum versus tight-binding models,” Phys. Rev. B 98(19), 195414 (2018). [CrossRef]  

52. R. Ribeiro-Palau, C. Zhang, K. Watanabe, T. Taniguchi, J. Hone, and C. R. Dean, “Twistable electronics with dynamically rotatable het- erostructures,” Science 361(6403), 690–693 (2018). [CrossRef]  

53. N. R. Finney, M. Yankowitz, L. Muraleetharan, K. Watanabe, T. Taniguchi, C. R. Dean, and J. Hone, “Tunable crystal symmetry in graphene–boron nitride heterostructures with coexisting moiré superlattices,” Nat. Nanotechnol. 14(11), 1029–1034 (2019). [CrossRef]  

54. S. S. Sunku, A. S. Mcleod, T. Stauber, H. Yoo, D. Halbertal, G. Ni, A. Sternbach, B.-Y. Jiang, T. Taniguchi, K. Watanabe, P. Kim, M. M. Fogler, and A. D. N. Basov, “Nano-photocurrent mapping of local electronic structure in twisted bilayer graphene,” Nano Lett. 20(5), 2958–2964 (2020). [CrossRef]  

55. L. J. McGilly, A. Kerelsky, N. R. Finney, K. Shapovalov, E.-M. Shih, A. Ghiotto, Y. Zeng, S. L. Moore, W. Wu, Y. Bai, K. Watanabe, T. Taniguchi, M. Stengel, L. Zhou, J. Hone, X.-Y. Zhu, D. N. Basov, C. Dean, C. E. Dreyer, and A. N. Pasupathy, “Seeing Moiré Superlattices,” Nat. Nanotechnol. 15(7), 580–584 (2020). [CrossRef]  

56. M. He, Y. Li, J. Cai, Y. Liu, K. Watanabe, T. Taniguchi, X. Xu, and M. Yankowitz, “Symmetry breaking in twisted double bilayer graphene,” Nat. Phys. 17, 26–302020. [CrossRef]  

57. T. Salamon, A. Celi, R. W. Chhajlany, I. Frérot, M. Lewenstein, L. Tarruell, and D. Rakshit, “Simulating twistronics without a twist,” Phys. Rev. Lett. 125(3), 030504 (2020). [CrossRef]  

58. A. L. Sharpe, E. J. Fox, A. W. Barnard, J. Finney, K. Watanabe, T. Taniguchi, M. Kastner, and D. Goldhaber-Gordon, “Emergent ferromagnetism near three-quarters filling in twisted bilayer graphene,” Science 365(6453), 605–608 (2019). [CrossRef]  

59. A. González-Tudela and J. I. Cirac, “Cold atoms in twisted-bilayer optical potentials,” Phys. Rev. A 100(5), 053604 (2019). [CrossRef]  

60. Kaiyuan Yao, Nathan R. Finney, Jin Zhang, Samuel L. Moore, Lede Xian, Nicolas Tancogne-Dejean, Fang Liu, Jenny Ardelean, Xinyi Xu, Dorri Halbertal, K. Watanabe, T. Taniguchi, Hector Ochoa, Ana Asenjo-Garcia, Xiaoyang Zhu, D. N. Basov, Angel Rubio, Cory R. Dean, James Hone, and P. James Schuck, “Nonlinear twistoptics at symmetry-broken interfaces,” arXiv:2006.13802, 2020.

61. Q. Fu, P. Wang, C. Huang, Y. V. Kartashov, L. Torner, V. V. Konotop, and F. Ye, “Optical soliton formation controlled by angle twisting in photonic moiré lattices,” Nat. Photonics 14(11), 663–668 (2020). [CrossRef]  

62. P. Wang, Y. Zheng, X. Chen, C. Huang, Y. V. Kartashov, L. Torner, V. V. Konotop, and F. Ye, “Localization and delocalization of light in photonic moiré lattices,” Nature 577(7788), 42–46 (2020). [CrossRef]  

63. D. Lowell, J. Lutkenhaus, D. George, U. Philipose, B. Chen, and Y. Lin, “Simultaneous direct holographic fabrication of photonic cavity and graded photonic lattice with dual periodicity, dual basis, and dual symmetry,” Opt. Express 25(13), 14444–14452 (2017). [CrossRef]  

64. D. Lowell, S. Hassan, M. Adewole, U. Philipose, B. Chen, and Y. Lin, “Holographic fabrication of graded photonic super-crystals using an integrated spatial light modulator and reflective optical element laser projection system,” Appl. Opt. 56(36), 9888 (2017). [CrossRef]  

65. V. A. Mandelshtam and H. S. Taylor, “Harmonic inversion of time signals,” J. Chem. Phys. 107(17), 6756–6769 (1997). [CrossRef]  

66. A. F. Oskooi, D. Roundy, M. Ibanescu, P. Bermel, J. D. Joannopoulos, and S. G. Johnson, “MEEP: a flexible free-software package for electromagnetic simulations by the FDTD method,” Comput. Phys. Commun. 181(3), 687–702 (2010). [CrossRef]  

67. S. Hassan, O. Sale, K. Alnasser, N. Hurley, H. Zhang, U. Philipose, and Y. Lin, “Broadband light-matter interaction due to resonance cavities in graded photonic super-crystals,” OSA Continuum 2(11), 3272–3280 (2019). [CrossRef]  

68. S. Lepeshov and Y. Kivshar, “Near-Field Coupling Effects in Mie-Resonant Photonic Structures and All-Dielectric Metasurfaces,” ACS Photonics 5(7), 2888–2894 (2018). [CrossRef]  

69. S. An, C. Fowler, B. Zheng, M. Y. Shalaginov, H. Tang, H. Li, L. Zhou, J. Ding, A. M. Agarwal, C. Rivero-Baleine, K. A. Richardson, T. Gu, J. Hu, and H. Zhang, “A Deep Learning Approach for Objective-Driven All-Dielectric Metasurface Design,” ACS Photonics 6(12), 3196–3207 (2019). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (5)

Fig. 1.
Fig. 1. (a) Schematic of arrangement of 8 wave-vectors of interfering beams for the generation of a twist angle α; (b) Schematic of optical setup for 8 interfering laser beams; (c-d) moiré interference patterns with a twist angle of 7.2 degrees at different z-location=0 and 0.27Pz, respectively. The unit super-cell, moiré period Pmoiré, lattice period Psmall, bright and dark regions are indicated in figures. (e) x-z cross section of moiré interference pattern.
Fig. 2.
Fig. 2. (a) Schematic of the arrangement of 12 wave-vectors of interfering beams for the generation of a twist angle α indicted by dashed blue line and dotted purple line; (b-c) moiré interference patterns with a twist angle of 9.9 degrees at z-location=0 and 0.55Pz, respectively. The moiré period Pmoiré, and bright regions are indicated. Dashed red circles indicate the interface between two distinct pattern regions.
Fig. 3.
Fig. 3. Simulated electric field intensity for resonance modes in moiré photonic lattice with a twist angle of 9.5°. Frequency Psmall/λ and Q-factor are labelled in each figure (a-l).
Fig. 4.
Fig. 4. Calculated Q-factors at different frequencies for resonance modes in Si moiré photonic pattern in Fig. 3 with a twist angle of 9.5 degrees. The red and purple rectangles indicate the locations of simulated photonic band gaps.
Fig. 5.
Fig. 5. Simulated electric field intensity for resonance modes at the dark regions (a), bright regions (b) and highly localized mode at the bright region (c) in Si moiré photonic pattern with a twist angle of 22.6 degrees. The unit super-cell is indicated by dashed blue squares.

Equations (9)

Equations on this page are rendered with MathJax. Learn more.

α = 2 arctan sin ( β 2 ) sin ( β 1 )
I ( r ) = i = 1 N E i 2 ( r , t ) + i < j N E i E j e i e j c o s [ ( k j k i ) r + ( δ j δ i ) ] .
P m o i r é 2 π 2 k s i n ( β 2 ) s i n 45 = 2 π 2 k × tan ( α 2 ) sin ( β 1 ) s i n 45 = P s m a l l tan ( α 2 )
α 1 = arctan ksin ( β 2 ) sin ( 60 ) ksin ( β 1 ) ksin ( β 2 ) cos ( 60 ) ,
α 2 = arctan ksin ( β 2 ) sin ( 60 ) ksin ( β 1 ) + ksin ( β 2 ) cos ( 60 ) ,
α = α 1 + α 2 2 arctan ksin ( β 2 ) sin ( 60 ) ksin ( β 1 ) ,
tan ( α 2 ) = p s m a l l p m o i r é sin ( 60 ) ,
π D 12 λ = 1.45 1.29 , for P s m a l l λ = 0.267 0.238 ,
2 D 12 λ = 0.82 0.92 , for P s m a l l λ = 0.267 0.238 ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.