1 Introduction
Non-Hermitian systems are constructed by introducing gainand-loss distributions [
1,
2] or nonreciprocal interactions [
3,
4], illustrating a good deal of unusual physics [
5]. Nonreciprocal (anisotropic) coupling, is characterized by unbalanced couplings
κab ≠
κba between two lattice sites
a and
b [
4]. Where
means that the mode amplitude undergoes gain or loss while couple between lattice sites
a and
b [
6,
7]. Systems that break Lorentz reciprocity are nonreciprocal and prevent light from retracing the forward path [
8,
9]. Nonreciprocity exists in topologically protected unidirectional edge states of topological photonics [
10,
11]. Introducing nonreciprocal coupling into non-Hermitian topological photonics leads to intriguing phenomena [
12], including non-Hermitian skin effect [
3,
4,
13], higher-order exceptional points (EPs) [
6], revised bulk-boundary correspondence [
14], and new definitions of topological invariants [
15].
Nonreciprocal systems were initially based on magneto-optical materials [
16]. Recently, several approaches have been developed to generate nonreciprocity, including parity-time-symmetric nonlinear cavities [
17], use of energy loss [
18], spatial–temporal modulation [
19,
20], and metamaterials [
21]. However, these principles for implementing nonreciprocal interlayer coupling have practical difficulties, particularly in topological photonic systems because they may lose original topological properties after adding non-reciprocity [
22].
In this paper, we provide a scheme for realizing the non-reciprocal interlayer coupling system by constructing on-site gain/loss in bilayer non-Hermitian topological systems. We reveal similarity transformations between nonreciprocal interlayer coupling and on-site gain/loss in the one-dimensional bilayer Su–Schrieffer–Heeger (SSH) model and two-dimensional bilayer C
6v topological photonic crystal (PC). The similarity transformations reveal that novel behaviors like delocalization [
23], skin effect [
4,
24], and breakdown of the conventional bulk-boundary correspondence [
14] are generic non-Hermitian phenomena not tied to a specific microscopic provenance of the non-Hermiticity [
5]. The topological number of the bilayer nonreciprocal interlayer coupling system, defined using a gauge-smoothed Wilson loop, is equal to that of the bilayer on-site gain-or-loss system. Topological phase transitions and parity-time-phase transitions of the non-Hermitian topological states occur by modulating the strength of nonreciprocal interlayer coupling or on-site gain/loss quantity. These results have great potential applications in reconfigurable laser arrays [
23,
25–
27], and for studying non-Hermitian topological physics, such as non-Hermitian band topology [
14].
2 Nonreciprocity-induced topological phase transition
The bilayer non-Hermitian SSH model [
28] is constructed by stacked nonreciprocal interlayer coupling photonic waveguide arrays (Fig. 1a). The alternating distance between in-layer nearest-neighbor waveguide determines
t1 (short hopping) and
t2 (long hopping) [
29–
31]. Following coupled-mode theory under tight-binding approximation and applying Fourier transformation [
6,
30], the Bloch Hamiltonian of the unit cell (black dotted box) under periodic boundary conditions (PBCs) is
where K is Bloch wave vector. Hm is the Hamiltonian of monolayer SSH model. H12 and H21 are nonreciprocal interlayer coupling matrices. See Appendix A for the forms of Hm, H12, and H21. We apply a similarity transformation to ∶
where and are two-by-two identity matrix and Pauli matrix, and is a four-by-four identity matrix. The Bloch Hamiltonian of bilayer on-site gain-and-loss SSH model is obtained as shown in Appendix A. κ is isotropic interlayer hopping (IIH) (gray arrows in Fig. 1b). The gain (and loss) strengths in gain (and lossy) waveguides are γ (Fig. 1b) (See Appendix B for generalized derivations with arbitrary gain/loss). When K = 0, the eight periodic-boundary-condition eigenvalues are
Figures 1c–j compare the bulk bands given by and . When γ = 0, there are four intersections (red pentagrams) (Fig. 1c). When γ = 0.5, an intersection becomes two EPs (magenta and pink crosses), which move away from each other along the first Brillouin zone (FBZ) as γ increases, while the complex energy region expands from the intersection to both sides until the edge of the FBZ (Fig. 1d). Then all EPs begin to move toward K = 0, while the complex energy region expands from edge to center of the FBZ (Fig. 1e). The real part of bands (RPBs) moves toward zero energy, and the central gap closes when (Fig. 1f). The RPBs approach and form two central degenerate points (DPs) when (black crosses) (Fig. 1g). When γ = t1−κ, two EPs merge into one EP at K = 0 (Fig. 1h). The RPBs form a central DP when (Fig. 1i). When γ = t1+κ, the other two EPs merge into one EP at K = 0 (Fig. 1j).
Hermitian SSH model is topologically nontrivial for
t1 <
t2 and trivial for
t1 >
t2. However, the bilayer structure can be changed from trivial phase to topologically nontrivial phase by increasing non-Hermitian quantities. The topological number of the bilayer non-Hermitian system is defined by the winding number [
32,
33], which can be calculated using a gauges-moothed Wilson loop [
33,
34]:
where and are the mth (m = 1, …, 8) right and left eigenstates of . Ki (i = 0, 1, 2,…, N) is discrete Bloch wave vector, and −π = K0 < K1 < K2 < … < KN = π, where N is a large integer number. Given the relation between the mth right and left eigenstates of and (See Appendix C for the deduced process.):
The Hermitian conjugate form of Eq. (6) is
By multiplying Eq. (7) with Eq. (5), we get , through which the gain-and-loss system is topologically equivalent to the nonreciprocal interlayer coupling system. has chiral symmetry , where . Em is the mth band, sorting the RPBs in ascending order. En and Em satisfy chiral symmetry for En = −Em*. Ej forms EPs with Em. For each m, we add and , and obtain four numbers characterizing topological number of the system (Fig. 2a). The four numbers change from 0 mod 1 to ±0.5 mod 1 at γc = 1.697, indicating a topological phase transition. The bilayer non-Hermitian system is topologically nontrivial when γ > γc and trivial when γ < γc, where is the topological phase boundary in the phase diagram (Fig. 2b).
Under open boundary conditions (OBCs), the waveguide array in each layer has 40 waveguides. Figure 2c compares real parts (RPs) and imaginary parts (IPs) of the open-boundary-condition
E −
γ relation given by
and
using the same parameters as above. The energy bands in Fig. 2c are colored according to the ratio,
ρ1, of the sum of field intensity (SoFI) of the four sites at the boundaries of the bilayer chain to SoFI of all sites.
ρ1 describes the degree of locality of the edge states' normalized field distributions. When
γ < 1.697, there are only bulk states. When
γ > 1.697, the RPs of the eigenvalues of the two pairs of degenerate topological edge states are close to zero, and the IPs of the eigenvalues are opposite. Figure 2d compares the normalized field distributions of four representative edge states of
and
when
γ = 3. The edge states of
are localized at boundaries of the SSH chain in the first (second) layer if the IPs of corresponding eigenvalues are positive (negative). The transformation matrix between
and
is
.
is shown in Fig. 2e, where
and
are the eigenstates of
and
. The bilayer non-Hermitian SSH model defined by waveguide arrays can be fabricated inside glasses using femtosecond-laser direct writing techniques [
35–
37]. A re-exposure technique can be applied to introduce point scatterers inside waveguides, making the system be non-Hermitian [
38].
3 Nonreciprocity-induced topological interface states
Given the C
6v PC with six sites per unit cell [
39], a topologically trivial or nontrivial bandgap is opened when intercell (
t1) and intracell (
t2) nearest-neighbor couplings are not equal [
39,
40]. With PBCs (OBCs) applied in the
x (
y) direction, the bilayer supercell of C
6v topologically nontrivial PC (
t1 >
t2) with zigzag-type domain walls [
41–
43] consists of 40 unit cells along
y direction per layer. The non-Hermitian domain walls are constructed by nonreciprocal interlayer coupling (Fig. 3a) and on-site gain–loss (Fig. 3b).
Using tight-binding approximation and Bloch theorem of periodic lattice, the Hamiltonian of the bilayer nonreciprocal interlayer coupling supercell is [
44]
where Hmono is the Hamiltonian of monolayer supercell without gain or loss, and κ denotes IIH. After the similarity transformation is applied to with , the Hamiltonian of the bilayer on-site gain-and-loss supercell is
where HGL and HLG are the Hamiltonians of the first layer and second layer with non-Hermitian domain walls, respectively. is identity matrix (N = 240). See Appendix D for the forms of Hmono, HGL, and HLG.
Figure 3c compares the periodic-boundary-condition E − γ relation given by and using parameters t1 = 5, t2 = 1, κ = 0.05, and K = 0. The eight eigenvalues whose RPs vary with γ are indicated in blue, and the corresponding IPs are indicated in red, as is the case for the bilayer non-Hermitian SSH model in Appendix A. Figures 3d − i compare parts of the projected bands given by and with different γ. Eight bands, which are new topological interface states localized at the bilayer non-Hermitian domain walls, appear in the bandgap. The projected bands of the eight DITISs (Fig. 4d − i) are similar to the bulk bands in Fig. 2e–j. However, the non-Hermitian domain walls cannot result in any new states in the topologically trivial PC (See Fig. 8 of Appendix D).
With OBCs applied in the x and y directions, the two-dimensional bilayer finite-size C6v topologically nontrivial PC with non-Hermitian domain walls consists of 10 (20) unit cells along the x (y) direction per layer. The Hamiltonians of the finite-size bilayer non-Hermitian domain walls constructed by nonreciprocal interlayer coupling and on-site gain–loss are and . Figure 4a compares the E − γ relation given by and . The energy bands are colored according to the ratio ρ2 of SoFI of the four sites on the boundaries of domain walls to SoFI of all sites. The eigenvalues whose RPs do not vary with γ are indicated in blue. The eigenvalues whose RPs vary with γ are indicated in purple, which are non-Hermitian DITISs localized at the bilayer non-Hermitian domain walls.
When
γ > 4.808, the normalized field distributions of two pairs of the degenerate edge of interface states (EOISs) of
are localized at the boundaries of non-Hermitian domain walls in the first (second) layer if the IPs of corresponding eigenvalues are positive (negative), and the RPs of corresponding eigenvalues of EOISs are close to zero (Fig. 4b, c). Figure 4d compares the normalized field distributions of two representative EOISs when
γ = 6. The normalized field distributions of EOISs of
are localized at the boundaries of the non-Hermitian domain walls in the first and second layer simultaneously. The transformation matrix between
and
is
.
is shown in Fig. 4e, where
and
are the eigenstates of
and
. The above results indicate that the topological numbers of bilayer non-Hermitian C
6v-typed DITISs can be defined as is the case for the bilayer non-Hermitian SSH model. The two-dimensional photonic systems can be experimentally realized at microwave frequencies. The photonic crystal platform is based on commercial alumina ceramics (Al
2O
3) with bandgap at microwave frequencies [
45,
46]. Al
2O
3 doped with chromium dioxide can introduce losses [
47], so the non-Hermitian control is achieved by doping or not doping chromium with Al
2O
3. The non-uniform dissipation distribution can be equivalent to the case of gain–loss distribution [
30,
48].
4 Conclusion
We have proposed a universal method to equivalently implement nonreciprocal interlayer coupling using on-site gain/loss in one-dimensional and two-dimensional bilayer topological systems through similarity transformation. The similarity transformation provides a convenient tool for understanding and implementing the non-Hermitian skin effect, especially in three-dimensional topological systems. The topological number of the bilayer nonreciprocal interlayer coupling system, which is defined using the gauge-smoothed Wilson loop, can be proved to be equal to the bilayer onsite gain-and-loss system. Topological phase transitions and parity-time-phase transitions of non-Hermitian topological states occur as a result of modulating the strength of nonreciprocal interlayer coupling or on-site gain/loss quantity. The topological origin of DITISs in the C
6v-typed domain wall can be understood via the bilayer non-Hermitian SSH model because they have the same form of transformation matrices, as is the case for the
E −
γ relation and eigenstate characteristics under both PBCs and OBCs. Our results offer new perspectives for studying non-Hermitian topological photonics and manipulating non-Hermitian topological states in bilayer non-Hermitian topological systems. We focused here on a photonic crystal for electromagnetic waves, but a similar lattice design may be applied to other bosonic systems, such as acoustic and mechanical structures [
23,
24]. The design principles should be generalizable to various frequencies including radio frequency [
49], microwave frequencies [
50], and optical frequencies [
13].