1State Key Laboratory of Transient Optics and Photonics, Xi’an Institute of Optics and Precision Mechanics of Chinese Academy of Sciences, Xi’an 710119, China
2University of Chinese Academy of Sciences, Beijing 100049, China
Parity–time () symmetric lattices have been widely studied in controlling the flow of waves, and recently, moiré superlattices, connecting the periodic and non-periodic potentials, have been introduced for exploring unconventional physical properties in physics, while the combination of both and nonlinear waves therein remains unclear. Here, we report a theoretical survey of nonlinear wave localizations in symmetric moiré optical lattices, with the aim of revealing localized gap modes of different types and their stabilization mechanism. We uncover the formation, properties, and dynamics of fundamental and higher-order gap solitons as well as vortical ones with topological charge, all residing in the finite bandgaps of the underlying linear-Bloch wave spectrum. The stability regions of localized gap modes are inspected in two numerical ways: linear-stability analysis and direct perturbance simulations. Our results provide an insightful understanding of soliton physics in combined versatile platforms of symmetric systems and moiré patterns.
1. INTRODUCTION
The moiré pattern—a periodic pattern overlapping its copy with a relative twist—as a novel two-dimensional (2D) material has shown great unique physical properties in condensed matter physics, including tunable flatbands, superconductivity, and correlated insulator phases at twist (magic) angles in twisted double-layer graphene [1–4], leading to an emerging field called twistronics—manipulating electronic properties through the relative twist angle [5]. Recently, the studies of moiré patterns and physics have entered regimes of optics and photonics: particularly, reconfigurable photonic moiré lattices were created in 2D photorefractive media by optical induction [6], and localization–delocalization transition of light [7] and optical soliton formation induced by the twisting angle [8,9] were observed, respectively, in linear and nonlinear contexts; magic-angle lasers with unique confinement mechanisms were fabricated in nanostructured moiré superlattices [10]; and multifrequency soliton generation in quadratic nonlinear media with commensurate–incommensurate photonic moiré lattices was predicted [11], to name just a few examples.
In addition, moiré optical lattices have been proposed in the context of ultracold atoms, and the associated moiré physics therein is being revealed [12–14]. To be specific, simulating twisted bilayers is possible by using cold atoms in state-dependent optical lattices that show Dirac-like physics and a band narrowing feature, making them ideal candidates to observe similar physics (such as strongly correlated phenomena in condensed matters) with larger rotation angles [12]. Simulating twistronics without a twist, a highly tunable scheme that rules out a physical bilayer or twist to synthetically emulate twisted bilayer systems in the setting of ultracold atoms trapped in an optical lattice, was proposed [13]. Moiré physics including tunable flatbands and Larkin–Ovchinnikov superfluids in spin-twisted optical lattices (rather than bilayers) was investigated [14], and considering that electromagnetically induced regular optical (or photonic) lattices via atomic coherence in atomic ensembles are mature in experiments, we recently proposed a related scheme to create electromagnetically induced moiré optical lattices in a three-level coherent atomic gas (either hot or cold atoms) in the regime of electromagnetically induced transparency [15]. Consequently, we can safely conclude that optics and ultracold atoms open flexible and promising routes toward the realization of moiré optical lattices and associated moiré physics.
Moiré patterns bridge the gap between periodic structures and aperiodic ones, offering a new platform for studying nonlinear localization of light [8,9]. Particularly, conventional periodic structures such as photonic crystals and lattices in optics [16–18] and optical lattices in the context of ultracold atoms [19–21] exhibit finite photonic or atomic bandgaps [22], the precise control of which and the corresponding nonlinearity could result in the emergence of a new spatially self-localized state called gap solitons (GSs) under repulsive (defocusing) nonlinearity [23–32]. Experimentally, optical GSs have been confirmed in optical Bragg gratings [33] and nonlinear photonic crystals [34], and the creation of atomic GSs of Bose–Einstein condensates (BECs) in optical lattices [35]. It is worth noting that parity–time () symmetric lattices as an interesting periodic structure were heavily studied in optical [36–41] and matter-wave [42,43] media and beyond in past years, providing a fertile land for investigating nonlinear waves including GSs [44–49]. However, to the best of our knowledge, the combination of moiré patterns and symmetry has not yet been reported.
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
In this work, we address the formation, properties, and dynamical stability of matter-wave GSs in a BEC trapped with symmetric moiré optical lattices constituted of two 2D lattices with a twist (rotation) angle. Spatially localized nonlinear excitations of three kinds of coherent matter waves—fundamental GSs and high-order ones grouped as two fundamental modes, as well as gap vortices with topological charge—situated inside the first and second atomic bandgaps of the linear Bloch-wave spectrum are found, underlining the tunable flatbands of the twisting periodic structure and robust stability of the reconfigurable nonlinear localized matter-wave structures within the associated gaps. The stability and instability properties of the nonlinear localized modes are assessed by linear-stability analysis and direct perturbance simulations, and both show good agreement.
2. RESULTS AND DISCUSSION
A. Theoretical Model
1. Gross–Pitaevskii Equation
The dynamics of a BEC cloud in 2D symmetric moiré optical lattices is described by the Gross–Pitaevskii (GP) equation for the dimensionless macroscopic wave function [19–21,49,50]: where Laplacian , and ; the last term represents the nonlinearity where the repulsive atom–atom interaction controlled by Feshbach resonance is chosen. We stress that in the context of nonlinear optics, the corresponding nonlinear Schrödinger equation for the propagation dynamics of field amplitude is deduced readily by substituting the time with propagation distance in Eq. (1). The moiré lattice under study yields where is the strength, is the imaginary potential strength, and the strength contrast between the two sublattices is defined as , setting for discussion. Note that the potential Eq. (2) matches the requirement of symmetry, , and it reduces to the usual (non-) moiré lattice at , and to the conventional lattice at . The plane is related to the rotation plane with a twisting angle :
2. Linear Tunable Flatband Properties
Depicted in Fig. 1 is the linear Bloch spectrum of the moiré optical lattice [Eq. (2)] under Pythagorean angle , with the Pythagorean triples for natural numbers . One can see from Fig. 1(a) that the widths of the first and second finite gaps compress rapidly with increasing at a defined angle [i.e., ], and non-Hermitian degeneracy arises at an exceptional point (singularity) , at which there is no Bloch gap. With an increase of strength contrast , more and more flatbands appear, widening the first gap while splitting the second gap, according to Fig. 1(b). For the 2D square moiré lattice at Pythagorean angle, the associated first reduced Brillouin zone in reciprocal space is given in Fig. 1(c). Typical real and imaginary parts of such lattices are shown in Figs. 1(d) and 1(e) at , and in Figs. 1(g) and 1(h) at ; the corresponding linear bandgap structures are respectively shown in Figs. 1(f) and 1(i), where there are broad first and second gaps.
Figure 1.Bandgap structures for 2D symmetric moiré optical lattices at with increasing imaginary potential strength (a) and strength contrast (b). (c) First Brillouin zone in 2D reciprocal space. Contour plots of the real (d), (g) and imaginary (e), (h) parts of the lattice (blue, lattice potential minima; red, lattice potential maxima) at (d), (e) and (g), (h), and their corresponding bandgap diagrams (f), (i) at and in reduced zone representation. I and II in (f), (i) represent the first and second bandgaps.
The stationary solution at chemical potential of Eq. (1) is given by (or with the propagation constant, in the context of nonlinear optics), yielding
Since we are interested in GS solutions supported by 2D moiré lattices, their stability property is a key issue, which is evaluated by linear-stability analysis. Thus, we perturb the solution as , where is the unperturbed solution constructed from Eq. (4), and and are small perturbations under eigenvalue . Substituting it into Eq. (1) would lead to the linear eigenvalue problem:
One can see from Eqs. (5) and (6) that the solution is stable only when all real parts of the eigenvalues are zero []; it is unstable otherwise.
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
B. Nonlinear Localized Modes and Their Properties
1. Fundamental Gap Solitons
The typical nonlinear localized mode in moiré optical lattices is the fundamental mode, matter-wave GSs, populated within atomic finite gaps, the characteristic profile of which is displayed in Figs. 2(a)–2(c). It is observed that the real wave function, , resembles a bright GS, while its imaginary part, , takes the form of a dipole. The condensate population, which is the number of atoms, , as a function of the chemical potential for GSs in moiré lattices at , is shown in Fig. 2(d), displaying an “anti-Vakhitov–Kolokolov” (anti-VK) criterion, , a necessary but not sufficient condition for the stability of GSs in periodic structures with repulsive (defocusing) nonlinearity [29–32]. For a given inside the first gap, the dependency between number of atoms and strength contrast , , is obtained in Fig. 2(e), revealing a slight vibration of when altering ; such a feature provides a flexible opportunity to launch fundamental GSs in moiré optical lattices with changeable strength contrasts, in addition to the twisting angle . Figure 2(f) depicts versus imaginary potential strength at and (within the second finite gap), showing a decrease tendency. We emphasize that the fundamental GSs are very stable in the midst of the first and second gaps; they are unstable as long as they are excited near both edges of the Bloch bands, and the directly perturbed evolutions shown in the following will prove it.
Figure 2.Typical profile of a fundamental GS supported by the 2D symmetric moiré optical lattice at (a)–(c). Corresponding real (a) and imaginary (b) parts, and contour plot of the module (c). Condensate population, , as a function of chemical potential (d), strength contrast (e), and imaginary potential strength (f) at . Other parameters: , in (a)–(c). in (f).
Besides the fundamental mode reported in Fig. 2, moiré optical lattices can also support higher-order spatially localized gap modes that may be considered as composite structures of several fundamental GSs. Two examples of such higher-order modes, composed of two out-of-phase and in-phase fundamental GSs, are displayed in Figs. 3(a), 3(b) and Figs. 3(c), 3(d), respectively. Their contour plots are shown in the second line of Fig. 3, with an emphasis on their tilted placements structured by twisting optical lattices; conforming to the fundamental counterpart in Fig. 2(b), the imaginary sections of the composite GSs emerge always as a dipole mode. The depending relations of both higher-order GS modes are summed up in Fig. 3(i), demonstrating, once again, the empirical stability of anti-VK criterion [29–32].
Figure 3.Typical profiles of higher-order GSs grouped as two out-of-phase (a), (b) and in-phase (c), (d) fundamental GSs at ; corresponding contour plots are displayed in the second line (e)–(h). (i) Condensate population, , as a function of chemical potential at (black, out-of-phase mode; red dashed, in-phase mode). Other parameters: , for B and , for C.
It is instructive to see whether moiré optical lattices can sustain robust stable complex localized gap modes represented as gap vortices with topological charge (winding number) . Our numerical calculations demonstrate that it is possible to create stable gap vortices with in such a novel periodic structure, provided that they are excited within finite gaps. Typical profiles of 2D hollow gap vortices with are in the form of four fundamental GSs entangled with phase, according to Figs. 4(a) and 4(b), where the corresponding real and imaginary parts of the wave functions and the associated phase structures are included. Evidently, both real and imaginary parts show the similarity of having positive and negative values; and counterintuitively, the imaginary wave function for the vortex gap mode does not exhibit a dipole-like feature for each fundamental GS, in contrast to their fundamental counterparts and higher-order ones as depicted in Fig. 2(b) and Figs. 3(b), 3(d), 3(f), and 3(h). The non-dipole feature of the imaginary parts of the gap vortices may be explained by the unique and inherent property of localized vortical modes—the phase factor, which, together with the structural property of moiré optical lattices, determines the wave structures (both real and imaginary parts) of gap vortices. This can be clearly observed from the corresponding phase structures of the gap vortices displayed in the right columns of the first two lines of Fig. 4. Also, in Fig. 4(c), we obtain the curve for such kinds of gap vortices with topological charge , showing an increase relationship.
Figure 4.Profiles of gap vortices consisting of four fundamental GSs with vortex charge prepared within (a) and near the upper edge (b) of the first finite gap; corresponding condensate population, , as a function of chemical potential at (c). Panels for the top and center lines denote, respectively, contour plot of the module, real and imaginary parts, as well as the associated phase structure. Other parameters for gap vortices marked by points D and E: (a) , ; (b) , .
The stability properties of all 2D localized gap modes (fundamental and higher-order GSs, gap vortices) have been measured in linear-stability analysis and direct numerical simulations of the perturbed evolutions, using Eqs. (5), (6), and (1), as pointed out previously. Our findings suggest the common feature of spatially localized modes supported by moiré periodic structures, that is, the GSs, is stable within the middle portions of finite gaps and unstable close to the band edges. For the fundamental GSs in Figs. 5(a) and 5(b), higher-order modes in Figs. 5(c) and 5(d), and gap vortices at topological charge in Figs. 5(e) and 5(f), with the stable ones being prepared in gaps and unstable ones near the band edges, their directly perturbed evolutions are depicted in the second line of Fig. 5, which are in good agreement with their linear eigenvalue spectra produced by linear-stability analysis in the bottom line of Fig. 5. It is important to emphasize that during evolution, the unstable localized gap modes develop multiple side peaks that reduce the necessary number of atoms () for sustaining the original gap modes; by contrast, shapes (good coherence) of stable modes can always remain. It must be pointed that the real eigenvalues in Figs. 5(a), 5(c), and 5(e) for stable GSs are not exactly zero, because we use the truncated matrix limited by our insufficient memory and CPU power, confirming the requirement of large spatial grid points in dealing with moiré optical lattices (compared to that with conventional ones).
Figure 5.Profiles of fundamental GSs (a), (b), higher-order GSs (c), (d), and gap vortices with (e), (f) prepared within the first (a), (c), (e), (f) and second (b), (d) finite gaps. Corresponding perturbed evolutions and linear eigenvalue spectra obtained from linear-stability analysis are displayed in the second and third lines, respectively. Other parameters: (a) , ; (b) , ; (c) , ; (d) , ; (e) , ; (f) , . White noise with 10% of its soliton’s amplitude is applied for all.
5. Stability Regions for Stable Localized Gap Modes
Before proceeding to the next section, we stress that the stability regions of all the localized gap modes considered above are accumulated in Table 1. It is seen that there are wide stable regions within the first two finite gaps, no matter what kind of localized gap mode, providing compelling evidence that localized gap modes are robust enough in moiré optical lattices with flatbands. Another common feature of all these localized modes is that they are always populated at the lattice potential minima, which change with the change of the twist angle-dependent lattice structure.
Stability Regions (Characterized by μ) of Nonlinear Localized Modes within the First and Second Finite Gaps
Localized Modes
Stability Regions ()
Fundamental gap solitons
and
Out-of-phase gap solitons
and
In-phase gap solitons
and
Gap vortices
and
C. Experimental Consideration
Although our focus is merely on BECs loaded onto moiré optical lattices, the physical system under study could directly apply to the nonlinear optics context for describing light propagation in moiré photonic crystals and lattices, since both contexts share the same model [51]. The only difference is to replace the time by propagation direction in Eq. (1), and the chemical potential by propagation constant in Eq. (4); then the wave function is replaced by the electromagnetic field amplitude of laser pulse , and the number of atoms would be the soliton power . Considering the fact that optical lattices have been successfully fabricated in optics (and photonics) [36–41] and ultracold atoms [42,43], the optical lattices of moiré type could also be practically realized using current-state-of-the-art experimental technologies. Specifically, this is achievable in optics using optical waveguides of semiconductor structure with alternating -symmetric gain–loss modulation [52] or distributed-feedback optical structures with gain or loss regions [53]. In terms of atomic gases, this is accomplishable in a gas of multilevel atoms such as the atomic system working under electromagnetically induced transparency [15,42,43]. Particularly, GSs have been created in both contexts [33–35], and the spatially localized gap modes predicted here can thus be envisioned in both experimental platforms.
D. Materials and Methods
All the numerical results presented above obey the following numerical procedure: the stationary GS solution is first found from Eq. (4) via a modified squared-operator iteration method [54]; then its stability is measured by means of the linear-stability analysis [Eqs. (5) and (6)] in the Fourier collocation method [54], and direct perturbance simulations [Eq. (1)] using the fourth-order Runge–Kutta method in real time.
Sign up for Photonics Research TOC. Get the latest issue of Photonics Research delivered right to you!Sign up now
3. CONCLUSION
We have addressed an as yet unresolved issue of the excitations and stability of 2D GSs supported by symmetric moiré optical lattices that exhibit a tunable flatband feature. Three categories of GSs are found, which are fundamental GSs and higher-order ones as well as gap vortices with topological charge, residing in both the first and second finite bandgaps of the associated diffraction diagram, and their stability is evaluated in linear-stability analysis and direct perturbance simulations. The experimental platforms for observing them are discussed, and we envision that the localized gap modes in symmetric moiré periodic structures are within reach in contexts of nonlinear optics and atomic media.
[50] B. I. Mantsyzov, R. N. Kuzmin. Coherent interaction of light with a discrete periodic resonant medium. Zhurnal Eksperimentalnoi i Teoreticheskoi Fiziki, 91, 65-77(1986).
[51] D. E. Pelinovsky. Localization in Periodic Potential: From Schrödinger Operators to the Gross-Pitaevskii Equation(2011).