1 Introduction
Photonic crystal (PC) cavities that can effectively confine light within small volume for extended durations have attracted significant attention owing to its potential applications in photonic integrated circuits [
1−
6]. Recent progress on topological photonics opens new routes for designing photonic cavity with strong robustness [
7−
13]. For example, by introducing the topological one-way edge states in a closed concave boundary, whispering-gallery-mode cavities featuring nonreciprocity and geometry independence are designed to produce nonreciprocal lasing [
14−
17]. However, the requirement of the time-reversal symmetry breaking hinders its extension in optical application. To overcome this limitation, researchers have leveraged the photonic pseudospin (such as orbital and valley) degree of freedoms to preserve the time-reversal symmetry. This approach has led to the discovery of the photonic topological insulator based whispering-gallery-mode cavities [
18−
20] and topology-controlled photonic cavities [
21−
23]. In addition, the higher-order band topology [
9,
24−
26], which goes beyond the conventional bulk-edge correspondences, enables the manipulation of the electromagnetic wave in lower dimensions, resulting in the proposal of the higher-order topology induced cavity [
27−
31].
On the other hand, photonic bound states, arising from the interplay between the topological defects (e.g., dislocation [
32,
33] and disclination [
34−
36]) and topological photonics are also considered promising candidates for photonic cavity [
37]. Generally, the underlying physics of these photonic bound states can be understood by solving the celebrated Dirac equations, where the topological defect acts as the Dirac mass term [
38,
39]. Depending on the nature of the spatially varying Dirac mass, the photonic bound states can be further classified as Jackiw−Rebbi and Jackiw−Rossi modes [
40,
41], corresponding to the kink (or monopole) and vortex solutions to the Dirac equations, respectively. For example, a photonic bound state induced by dislocation was demonstrated in a topological PC, fully explained by a model of sign switch of the mass term of the one-dimensional (1D) Dirac equation [
32]. Moreover, employing the strategy of Kekulé modulation, the Dirac vortex cavities are experimentally demonstrated in systems with spatial varying mass distribution [
42,
43], finding potential application on topological cavity surface emitting laser [
44,
45]. Nevertheless, most of existing proposals rely on the complex geometric designs or multiple-parameter fine-tuning.
To address the need for compact cavity design, we propose a concept of gradient dislocation, formed by placing two chunks of square-lattice PCs with same width but differ by a periodic number together. The gradient distribution of the relative displacement between upper and lower PCs is described by a 1D position-dependent mass system. Crucially, the signs of the Dirac mass on the two sides are opposite, giving rise to a mass domain-wall at the gradient dislocation. We demonstrate the Dirac cavity mode, as a zero-energy solution to the 1D Dirac equation, exhibiting the smallest modal area and a relatively high Q-factor. The robustness of the Dirac cavity mode and its possible application on the index sensing are also discussed.
2 Results and discussion
We start by introducing the concept of gradient dislocation. Considering a finite square-lattice PCs, a column of the PCs is removed from the upper-half of the square-lattice PCs. The missing lattice constant is then uniformly distributed to each column of the upper-half PCs, resulting in a gradient dislocation. From another perspective, the gradient dislocation is formed by adjoining two finite PCs with identical width but differing one period along the -direction. As shown in Fig.1(a), the upper-half PCs has periods along -direction (lattice constant is denoted as ) while the lower-half PCs has period along -direction. To maintain the same width, the lattice constant of the lower-half PCs gives . Obviously, the difference between upper-half and lower-half PCs can be almost ignored when the lattice periodic number is sufficiently large.
To elucidate the formation of Dirac cavity mode in a gradient dislocation, we select three interface configurations at left-/right-side, and middle of the gradient dislocation for consideration. For the middle area of the gradient dislocation, glide symmetry guarantees the emergence of a pair of gapless interface states [
33]. However, away from the middle area, the breaking of the glide symmetry leads to a pair of gapped interface states. Since the gradient dislocation is symmetric with respect to the center line, the gapped interface states exhibit opposite Dirac mass terms at left- and right-sides of the gradient dislocation. Consequently, a mass domain-wall forms at the center of the gradient dislocation [see Fig.1(b)]. According to the Jackiw−Rebbi theory, there will be a photonic bound states localized at the gradient dislocation center.
To realize the above scenarios, we design a square-lattice PCs comprising a single cylinder rod with permittivity in each unit cell. The rod diameter is , where refer to the lattice constant. As shown in Fig.2(a), we consider two primitive cell configurations: UC1, where the rod is split into two halves positioned at the upper and lower edges, and UC2, where the rod is split into four quarters positioned at the four corners. Throughout this work, we consider 2D transverse-magnetic harmonic modes and all the numerical simulations are carried out using the commercial software COMSOL Multiphysics based on finite element method. The photonic bands for UC1 and UC2 are shown in Fig.2(b). Here, we focus on the band gap between the first and the second bands, taking the bulk polarization (or the Wannier center) as the bulk topological invariant. Intuitively, the Wannier center for UC1 resides midway between the upper and lower edges, while for UC2 it resides at the corners of the unit cell (see Appendix A for details). We then consider an interface formed by adjoining UC1 and UC2, as depicted in Fig.2(c). In this configuration, the Wannier center in the lower half PC is shifted by half a lattice constant relative to that in the upper half PC. Eigenmode calculations in Fig.2(d) indicate that the interface formed by UC1 and UC2 supports a pair of gapless interface states, as a consequence of a synthetic Kramers theorem induced by the combination of the glide symmetry and time-reversal symmetry. From another perspective, a relative displacement parameter , i.e., , between upper and lower interface generates the glided-symmetric interface. Deviations from break glide symmetry, resulting in the gapped interface states [see the projected band structure of the interface states with in Fig.2(d)]. We emphasize that the Wannier center is ill-defined when the displacement parameter is not integer times of the half lattice constant. To this end, we employ a Dirac mass, which defined as the frequency gap between two projected interface bands and is proportional to the overlapping integral of the electromagnetic field of the two interface states at , to characterize the topological nature of the interface state gap. For example, the right panel of Fig.2(d) shows that the electric field patterns of the interface states at for the upper and lower bands with and . The overlapping integral results indicate a sign flip in the Dirac mass despite their identical projected band structure.
We then construct a finite square-lattice PCs with a size of , which is bounded with perfect matched layers, and remove a column of the upper-half PCs and distribute the missing lattice constant to each column of upper-half PCs uniformly [see Fig.2(e)]. For the th column of the upper-half PCs, the accumulated relative displacement PCs is , where refer to the relative displacement per column. Fig.2(f) shows the evolution of the band edges of the interface states versus relative displacement in Fig.2(f). It is seen that interface states exchange its frequency ordering upon crossing , indicating that the interface configured with and are of opposite Dirac masses. Furthermore, Fig.2(g) plots the eigen spectra of the interface states for five representative cases, namely, , and . The frequency gap first decreases from to , and then increase from to . The underlying physics of the gradient dislocation can be described by the 1D Dirac equation with a position dependent and sign switch mass distribution (see details in Appendix B). Notably, the Dirac mass has opposite signs on the left and right sides of the gradient dislocation center. According the Jackiw-Rebbi theory, this suggests an emergent photonic bound state localize at the center of the gradient dislocations. As expected, the eigen calculation reveal a photonic bound state — the Dirac cavity mode — with even-symmetry localized at the center of the gradient dislocation [see Fig.2(h)]. In addition to the Dirac cavity mode, we notice that there exist other bound modes in the bulk band gap, which are also the solutions of the Dirac equation and can be interpreted as standing wave solutions [also see in Fig.3(a)]. Furthermore, we also implement an eigen calculation for a gradient dislocation with 52 periods along -direction. Interestingly, the eigen mode in Fig.2(i) reveals a photonic bound states with odd symmetry emerging at the center of the gradient dislocation.
Since a larger periodic number
corresponds to a larger gradient structure with smaller relative displacement per column, the periodic number
plays a vital role in realizing high-performance Dirac cavity mode. Fig.3(a) shows the eigen spectra of the finite gradient dislocation versus periodic number. It is observed that the frequency of the Dirac cavity mode gradually increases and towards to pin at
, which corresponding to the frequency of degenerated interface states with
. Additionally, we also notice that other bound states within the bulk band gap correspond to standing wave solutions of the 1D Dirac equation and are identified as higher-order Dirac cavity modes (see details in Appendix B and Appendix C). We then pay special attention to the performance of Dirac cavity by studying its
Q-factor and modal area. As shown in Fig.3(b), the
Q-factor increase rapidly for
<40, and nearly saturates for
>40. The reason is that there are two decay mechanisms for the Dirac cavity mode: it can decay into the upper and lower PCs or it can decay along the gradient dislocation, which are characterized by the two different decay rates
and
[
1]. Since that
increases exponentially with
, while
is roughly
-independent, the overall
Q-factor
increases exponentially initially with
and then saturated at a certain
. For the Dirac cavity mode, the
Q-factor saturated at nearly
around
=40, while the higher-order Dirac mode I and II exhibit lower and higher saturated
Q-factor, respectively. For comparison, we calculate both the
Q-factor and the modal area of all eigen modes within the bulk band gap for
. The results in Fig.3(c) show that the Dirac cavity mode possesses a relatively high
Q-factor with lowest modal area, indicating superior performance for applications. Alternatively, we compute the distribution of topological cavity mode in momentum space by utilizing the Fourier transformation. Fig.3(d) shows that the Dirac cavity mode is mainly localized at the edge of the Brillouin zone
and lies outside of the line cone indicated by the dashed circle.
To further distinguish the Dirac cavity mode from higher-order Dirac mode, we consider a perturbed gradient dislocation by inserting a perfect electric conductor (PEC) into different positions of the gradient dislocation [see Fig.4(a)]. The distance between the PEC and the gradient dislocation center is denoted as . We implement the eigen calculation of the perturbed gradient dislocation with ranging from 0 to , and the results are presented in Fig.4(b). When the PEC is inserted at the center , the persevered mirror symmetry of gradient dislocation results in two degenerate Dirac cavity modes with odd and even symmetry, respectively [see Fig.4(c)]. Inserting the PEC away from the center (e.g., ) breaks mirror symmetry and significantly distorts the Dirac cavity mode [see the lower panel of Fig.4(b)]. Remarkably, the eigenfrequency of Dirac cavity mode remains unchanged [see the red points in Fig.4(b)]. In contrast, the eigenfrequencies of higher-order Dirac cavity modes shift, and their field distributions are disrupted. Fig.4(d) presents the Q-factor and the modal area of the Dirac cavity mode versus PEC position . It is seen that both the Q-factor and the modal area exhibit a sudden decrease when . This occurs because that the field profile of the Dirac cavity mode is largely distorted. As the PEC moves away from the center of gradient dislocation, both Q-factor and the modal area recover to their original levels as the Dirac cavity field profile restores.
Finally, we propose a coupled cavity-waveguide system by combining a gradient dislocation with a glide-symmetric waveguide. As shown in Fig.5(a), the gradient dislocation is formed by adjoining two finite PCs (orange and blue) of identical width but differing by one period. The boundary is surrounded by perfect matched layers. For convenience, we set the lattice constant nm. A glide-symmetric waveguide is introduced by shifting the lower PCs (green area) relative to the middle PCs (blue area) by half a lattice constant. The distance between the glide-symmetric waveguide and the gradient dislocation is denoted as . By tuning the distance , the coupling strength between the Dirac cavity and glide-symmetric waveguide can be effectively controlled. As shown in Fig.5(b), transmission dips confirm that the propagation of the interface states is strongly perturbed near the cavity eigenfrequencies. Utilizing the coupled cavity waveguide system, we further explore its potential application on refractive index sensing. Fig.5(c) presents the transmission spectra for refractive index change as the refractive index of the surrounding medium increases from to at an interval of . It is observed that the resonance dip redshifts with increasing of , indicating a strong dependence of the Dirac cavity mode frequency on the surrounding refractive index. Fig.5(d) plots the resonance wavelength shift () versus refractive index change () for different coupling strengths (controlled by ) between the Dirac cavity and the glide-symmetric waveguide. Linear fitting yields a refractive index sensitivity for , demonstrating the potential of the Dirac cavity for index sensing application.
3 Conclusion
In conclusion, we propose an ultracompact photonic cavity design based on gradient dislocation. The underlying physics is well described by 1D Dirac equation with a position-dependent mass term that switches sign. The resultant photonic bound states — Dirac cavity mode — emerges a kink solution to the Dirac equation, localized at the center of gradient dislocation. We demonstrate that the Dirac cavity mode exhibits an ultrasmall modal area with relatively high Q-factor compared with other higher-order Dirac cavity modes. We also discuss the robustness of Dirac cavity and its potential application in refractive index sensing. Our ultracompact cavity design is valuable for scalable and tunable integrated photonic devices and can be generalized to other classical wave systems.