Higher-order exceptional point and Landau-Zener Bloch oscillations in driven non-Hermitian photonic Lieb lattices

We propose a scheme to realize parity-time (PT) symmetric photonic Lieb lattices of ribbon shape and complex couplings, thereby demonstrating the higher-order exceptional point (EP) and Landau-Zener Bloch (LZB) oscillations in presence of a refractive index gradient. Quite diﬀerent from non-Hermitian ﬂatband lattices with on-site gain/loss, which undergo thresholdless PT symmetry breaking, the spectrum for such quasi-one-dimensional Lieb lattices has completely real values when the index gradient is applied perpendicular to the ribbon, and a triply degenerated (third-order) EP with coalesced eigenvalues and eigenvectors emerges only when the amplitude of gain/loss ratio reaches a certain threshold value. When the index gradient is applied parallel to the ribbon, the LZB oscillations exhibit intriguing characteristics including asymmetric energy transition and pseudo-Hermitian propagation as the ﬂatband is excited. Meanwhile, a secondary emission occurs each time when the oscillatory motion passes through the EP, leading to distinct energy distribution in the ﬂatband when a dispersive band is excited. Such novel phenomena may appear in other non-Hermitian ﬂatband systems. Our work may also bring insight and suggest a photonic platform to study the symmetry and topological characterization of higher-order EPs that may ﬁnd unique applications in for example enhancing sensitivity.

Introduction.-In the past decade, the physics of dispersionless flatbands has attracted considerable interest, partly due to their connection to various fascinating phenomena including the fractional quantum Hall effect [1][2][3], high-temperature superconductivity [4], and tunable strictly compact nonlinear excitations [5][6][7][8], among many others.Flatband phenomena also emerge in twisted bilayer graphene lattices that serve as a new platform for exploration of correlated phases [9][10][11].The flatband eigenmodes are typically compact localized states (CLSs) in real space due to local destructive interferences of the Bloch wave functions.Such CLSs have been observed in various systems, particularly in artificial photonic lattices [12][13][14][15].For instance, the first direct experimental observation of CLSs was realized in Lieb photonic lattices [12,13], which stimulated a surge of theoretical and experimental studies devoted to flatband physics -for a survey of the state of the art, see [16][17][18][19][20]. Efforts in realizing synthetic magnetic flux for light have allowed to observe Aharonov-Bohm caging in photonic rhombic lattices, an effect that originates from flatbands [21,22].Recently, investigations on flatbands have been extended to parity-time (PT) symmetric non-Hermitian systems [23][24][25].The existence of secondorder non-Hermitian degeneracies and the higher-order exceptional points (EPs) have led to many intriguing features such as unidirectional invisibility, topological phase transitions, PT-symmetric lasers, and even enhanced sensitivity [26][27][28][29][30][31][32][33].It has been shown that flatbands with completely real eigenvalues can appear in many non-Hermitian systems.For example, a bipartite sublattice symmetry enforces the existence of non-Hermitian flatbands, which are typically embedded in an auxiliary dispersive band -akin to the flatband bound states in the continuum [34].Defect states can emerge from a non-Hermitian flatband of photonic zero modes [35].In addition, flatbands can arise due to the interplay between synthetic magnetic flux and non-Hermiticity [36,37].Recently, non-Hermitian flatband lattices have been demonstrated in experiment using laserwritten photonic waveguide arrays with tailored loss distributions [38].
On the other hand, Bloch oscillations, the oscillatory motion of particles in a periodic potential driven by an external force, is one of the most striking wave dynamics in condensed matter physics [39].Landau-Zener tunneling between Bloch bands, a counteracting effect limiting Bloch oscillations, appears at the edges of the Brillouin zone where the band gap is shallow [40].While Bloch oscillations have been demonstrated in various wave systems including optical waves [41][42][43][44], typical examples of Landau-Zener Bloch (LZB) oscillations are better represented by their classical analogue observed in photonic waveguide arrays [45,46].Moreover, flatband Bloch oscillations together with triply Landau-Zener tunneling between the flatband and dispersive bands have also been proposed in optical lat-tices [47][48][49].Quite different from the Hermitian systems, the complex phase in non-Hermitian lattices indicates nonconserved energy, resulting in super-Bloch oscillations with amplification or damping, and even chiral Zener tunneling at the non-Bloch-band collapse point [50,51].In PTsymmetric photonic lattices, noticeable secondary emissions have been observed due to the degeneracy of the EPs [52,53].Very recently, Bloch oscillations and Landau-Zener tunneling have been investigated in a PT-symmetric flatband rhombic lattice by introducing on-site gain/loss [54].However, due to the band touching between flatband and dispersive bands [55,56], such a system does not support any unbroken PT phase or completely real spectra once the non-Hermitian parameters are introduced, as generally encountered in non-Hermitian flatband systems [57][58][59].Meanwhile, peculiar features like secondary emissions unique to the non-Hermitian systems remain undiscovered in flatband lattices as far as we know.Thus, a natural question arises: can we establish a non-Hermitian flatband system that will exhibit (nonzero) threshold PT symmetry breaking, and how would the LZB oscillation dynamics change when the interplay between the flatband, EP and the driving field comes to play the role in such a system?
In this work, we propose a scheme to establish quasi-onedimensional (1D) photonic PT-symmetric Lieb lattices (the Lieb ribbon), and to address the above questions by considering the LZB oscillations in the non-Hermitian flatband system.By introducing asymmetric complex-valued couplings and a refractive-index gradient perpendicular to the ribbon, we obtain the flatband lattice possessing the PTsymmetric phase and a triply degenerated higher-order EP.We show that, in such a non-Hermitian flatband system, PT symmetry breaking occurs only when the amplitude of gain/loss ratio exceeds a nonzero threshold.Moreover, an asymmetric energy distribution and pseudo-Hermitian propagation can take place for small non-Hermitian parameters.Importantly, we uncover that a wave-packet can become either amplified or damped during the LZB oscillations accompanied with strong secondary emissions in a non-Hermitian system.
A ribbon of Lieb photonic lattice composed of an array of periodically arranged evanescently coupled waveguides is shown in Fig. 1(a), which consists of five sites (A, B, C, D and E) per unit cell.Such a structure represents the quasi-1D form of two-dimensional Lieb lattices, which have attracted great interest in the past several years and have been employed to demonstrate fundamental flatband phenomena [12][13][14]60].Although in a simple flatband lattice geometry, the Lieb ribbon allows to clarify the effects of the refractive index gradient fields on the flatband CLSs and LZB oscillations in a PT-symmetric non-Hermitian system.The external driven fields are applied by applying linear refractive index gradients either parallel (x-gradient) and perpendicular (y-gradient) to the ribbon.The non-Hermiticity is introduced by employing asymmetric complex-valued couplings (t 1 and t 2 ).Such non-Hermitian couplings can be realized by embedding gain or loss media between adjacent waveguides or periodically modulating the on-site gain/loss amplitude ratio [34,61,62].The Hamiltonian of the system in the tight-binding approximation can be written as:  , c) Calculated real and imaginary parts of the spectrum as a function of g when no index gradient is present (i.e., ∆βx = 0 , ∆βy = 0).The red dashed circles represent the position of the EP, and the pink filled circles show the value in which the upmost and lowest dispersive bands become complex.(d) A flatband CLS occupies two plaquettes for g = 0 with colored and filled circles representing the sites of nonzero energies, and the characters right next to them show the corresponding amplitudes.For simplicity, we set t = 1, and the lattice period is 1. .
e n are the creation and annihilation operators in the n-th unit cell on the A, B, C, D and E sites, respectively.t, t 1 and t 2 are the coupling coefficients.We incorporate an arbitrary amount of PT-symmetric gain and loss into the system in a balanced way by setting t 1 = t (1 + ig), t 2 = t (1 − ig), where g is the non-Hermitian parameter [34,58].Here ∆β x and ∆β y denote the wave-number difference between adjacent waveguides and thus define the x-gradient and ygradient, respectively.The effective propagation constant (or the on-site energy) of the waveguides is determined by the refractive index gradient strength.We assume that the sites of the same unit cell have on-site energy difference ∆β y along y direction in the presence of y-gradient, while the effective propagation constant of the n-th unit cell is shifted by ∆β x compared to the (n+1)-th unit cell for the x-gradient.Meanwhile, sublattices A, C and E of the nth unit cell have the same effective propagation constant which is shifted by ∆β x /2 compared to the B (D) in the same unit cell.This arrangement can be described by the discrete T is the five-component wave function describing the field amplitude in unit cell n, and the momentum space Hamiltonian H k is given by ∆β y and Λ E = n∆β x − ∆β y are the gradients applied to the A, B, C, D and E sites, respectively.)Let us start from considering the case without the external field, i.e., ∆β y = 0, ∆β y = 0, and calculate the eigenvalues β as a function of g.The real and imaginary parts are depicted in Figs.1(b) and 1(c), respectively.The Hermitian regime g = 0 (where both t 1 and t 2 are real) reveals that the lattice has five bands distributed symmetrically with respect to the central flatband β = 0, which touches two dispersive bands at k = π as marked by red circles.The system undergoes thresholdless PT symmetry breaking, i.e., the central two dispersive bands become complex for any value of g = 0. Consequently, the band touching point is also the EP that representing the transition between the real and complex band regions.At the same time, band gap width of upmost and lowest dispersive bands decreases with the increasing of g.For g = √ 2/2, all the four dispersive bands are in the broken PT symmetry phase, as indicated by the pink filled circle.Nevertheless, it can be clearly found that the flatband β = 0 survives due to the chiral symmetry of the lattice [63,64].Moreover, different from uniform lattices (t 1 = t 2 ) where a CLS occupies four sites in a single plaquette [Fig.1(a)], the introduction of asymmetric couplings requires that a flatband CLS occupies at least two plaquettes.The irreducible CLS is sketched in Fig. 1(d) for t = 1, with B(D) and C sites have unequal amplitudes and opposite phase, ensuring destructive interference in the neighboring A(E) sites.Note that such a state becomes overlapping of two fundamental CLSs in Fig. 1(a) for g = 0 (t 1 = t 2 = t = 1).Next, we con-FIG.2. Same as in Figs.1(b-d) but with a y-gradient ∆βy = 0.3t (t = 1).Due to the presence of y-gradient field, two symmetric gaps open for g = 0 and a triply degenerated EP (marked by red filled circles) with nonzero value appears as g is increased to gc = 0.17.The new CLS occupies two plaquettes, where extra A and E sites are included, with opposite phases (represented by blue and green colors) and same amplitude equalling field strength ∆βy.
sider the case where a y-gradient is applied.Figures.2(a) and 2(b) show the typical spectrum for ∆β y = 0.3t (t = 1, and thus ∆β y = 0.3).The band degeneracy is lifted with two symmetric gaps emerge between the flatband and dispersive bands for g = 0.With the increase of the non-Hermitian parameter g, the band gap decreases and coincides with the PT symmetry breaking transition for critical value g c = 0.5 1 − ∆β 2 y − 1 − 4∆β 2 y = 0.17, at which the eigenmodes of the dispersive bands coalesce into a degenerated one.More importantly, the flatband still survives, and different from previously studied non-Hermitian systems, here the EP becomes triply degenerated in virtue of the existence of flatband.Accordingly, all the eigenvalues are completely real for g < g c , and a flatband lattice with nonzero threshold PT symmetry breaking is obtained.For g > g c the system turns into the symmetrybroken phase.As g is further increased to cross the value g = 0.5 1 − ∆β 2 y + 1 − 4∆β 2 y = 0.65, the highest and lowest bands merge and two pairs of complex-conjugates eigenvalues are observed.Furthermore, the CLS shown in Fig. 1(a) is no longer an exact solution of the lattice, as for ∆β y = 0 non-zero amplitudes appear in sites A and E as well as the y-gradient remodulates the amplitude in the central C site [Fig.2(c)].
To verify our analytical results, we numerically simulate the evolution dynamics of CLSs for g = 0.3 with propagation distance z=6 mm.We consider 260 unit cells of waveguide arrays arranged according to Fig. 1 and excite the central two plaquettes.Figures.3(a1) and 3(b1) show the intensity and phase of the input beam Ψ n similar to that of flatband CLSs shown in Figs.1(d) and 2(c), respectively.The first row of Fig. 3 shows the corresponding results of lattices without the refractive index gradient.It can be clearly seen that the initial excitation propagates without any diffraction due to the excitation of flatband modes [Fig.3(a3)].Furthermore, the initial intensity remains constant and consequently is conserved even though the system is in a PT symmetry broken phase.Moreover, the phase measurement obtained by interfering the output with an inclined plane wave further confirms that the initial phase structure is well preserved [Fig.3(a4)].Instead, if the input beam has equal phase, a discrete diffraction pattern is obtained [Fig.3(a2)] and the power is not conserved as shown in the inset.Similar results can be found for the case ∆β y = 0.For better visibility, we set ∆β y = 0.5t (A and E sites become more visible).The excitation [Fig.3(b1)] now contains two more sites and stays invariant during propagation [Fig.3(b3)].However, due to the presence of y-gradient, the discrete diffraction of equal phase input becomes asymmetrical [Fig.3(b2)].It should be noted that the system is in PT-symmetry unbroken phase as the critical value for the PT transition is set to g c = 0.43 in this case.Consequently, the quasi-energy Q = n Ψ (n) Ψ * (−n) is also conserved in Fig. 3(b2) [50].
FIG. 3. Evolution dynamics of the flatband CLSs in PTsymmetric photonic lattices (first row) without refractive index gradients ∆βx = 0, ∆βy = 0 and (second row) with a y-gradient ∆βy = 0.5t.First column: input beams.Insets show the phase structure, where white region represents π phase and gray area represents 0 phase.Second column: in-phase output with discrete diffraction appears.Insets show the evolution of quasienergy Q, which is conserved in (b2) as the system is operated below the PT symmetry threshold value.Third column: outof-phase output with each image normalized so that the total intensity is 1.Though the system is non-Hermitian, the original excitations do not lead to any diffraction and the total energy remains constant.Fourth column: corresponding zoom-in interferograms of the third column.We set t=250 mm −1 , g = 0.3, and propagation length z=6 mm.Now, we consider the case for ∆β x = 0 and investigate the flatband CLSs dynamics in the PT-symmetric Lieb lattice ribbon.In our platform, applying a constant refractive index gradient parallel to the lattice is similar to a constant dc electric field applied to a crystal, and thus resulting in Bloch oscillations of the wave-packet.Linear combinations of CLSs shown in Fig. 1(a) with a broad Gaussian distribution b n = −c n = d n = e −n 2 /2σ 2 (σ = 15) is set as an input excitation at normal incidence [Fig.4(a)].Different from the rhombic lattices where flatband and CLSs do not change in the presence of x-gradient [47,49], it is easy to find that the input CLSs no longer stay strictly localized for ∆β x = 0.However, the central band remains nearly flat when x-gradient is small, and the input still can be used to excite flatband [48].Figures . 4(b-e) depict the numerical results with ∆β x = 0.05t.(σ = 15) at normal incidence.Colored sites represent nonzero amplitudes with blue and green sites having π phase difference.(b) ∆βy = 0, g = 0.005, the input keeps compactness with slight and symmetric oscillations for non-Hermitian lattices in the absence of y-gradient.(c) ∆βy = 0.12t, g = 0.005, the total energy stays constant over many periods and the energy distributes asymmetrically for non-Hermitian lattices in the presence of y-gradient.(d) ∆βy = 0.12t, g = 0.015, the energy of the wave-packet is amplified during oscillations.(e) ∆βy = 0.12t, g = 0, initial state evolves in a symmetric way for Hermitian lattices in the presence of y-gradient for direct comparison.The oscillation period is l = 2π/∆βx and the propagation length L = 5l is divided into 2000 steps.For all figures, the input peak intensity is normalized to 1.
For the non-Hermitian case ∆β y = 0, g = 0.005, the input remains almost invariant with only sightly symmetric Bloch oscillations observed [Fig.4(b)].In this case, the spectrum is symmetric and the gap width is zero [see Fig. 1(b)], and the initial state ceases to be a linear combination of CLSs depicted in Fig. 1(d) with a part of the wave-packet splitting into dispersive bands.In addition, the system is PT symmetry broken and therefore the energy is obviously amplified.Noticeably different features can be found for the cases of ∆β y = 0, g = 0 [Figs.4(c) and 4(d)].More specifically, we observe sudden switches from an almost unaltered state into an oscillating pattern and Landau-Zener tunneling due to the triply interband transitions.The beam intensity starts to exhibit an asymmetric energy distribution with the increase of g.This asymmetry is likely attributed to the asymmetry of the energy spectrum, i.e., the spectrum become asymmetric with the simultaneous appearance of xgradient, y-gradient and non-Hermitian coupling, which is similar to the cases with magnetic field in the Hermitian systems [47].Note that, compared with the case in Fig. 2, the Wannier-Stark ladder spectrum is complex-valued even though g < g c (∆β y = 0.12t, and the critical value at EP is g c = 0.06) [50][51][52][53].This is because the external xgradient field breaks the PT symmetry of the Hamiltonian.Remarkably, we find that a pseudo-Hermitian propagation is restored for very small values of g (g g c ), where the average total power remains constant over many periods of Bloch oscillations, see Fig. 4(c) for g = 0.005.Meanwhile, due to the complex band structure, the wave-packet is amplified each time passing through the EP for bigger value of g, see Fig. 4(d) for g = 0.015.For comparison, we also investigate the Hermitian case in the presence of y-gradient in Fig. 4(e) (∆β y = 0, g = 0).In this case, the initial state is conserved and it evolves in a symmetric way, reflecting the band symmetry shown in Fig. 2(a) for g = 0.One of the most striking characteristics of non-Hermitian Bloch oscillations is the existence of the so-called secondary emission that occurs at the EP [52,53].The degeneracy at the EP causes a redistribution of the initial excitation.Therefore, each time the wave-packet passes through the EP, a new branch of Bloch oscillations appears.However, such a phenomenon is not obvious when the central flatband is excited.To reveal this effect, we simulate the Bloch oscillations of the dispersive bands.Different form the excitation in Fig. 4, the input beam has equal phase structure and all the sublattices including A and E are excited (a n = b n = c n = d n = e n = e −n 2 /2σ 2 ) [Fig.5(a)] .In the Hermitian case g = 0 [Fig.5(b)], this initial excitation yields sinusoidal energy oscillations together with only weak Landau-Zener tunneling between two dispersive bands.Moreover, no energy distributes in flatband.The situation nevertheless changes drastically in non-Hermitian regime, as for g = 0.015 the energy starts to evolve into the flatband [Fig.5(c)].Such an effect appears more clearly with the increase of g, as the secondary emissions occur each time when the oscillatory motion passes the triply degenerated EP, leading to distinct energy distribution into the flatband and strong energy tunneling between the flatband and dispersive bands [Fig.5(d)].Let us notice that in both cases Figs.5(c) and 5(d) the total energy decreases along the propagation.This can be changed so that the total energy becomes amplified along the propagation while preserving the oscillations observed in Fig. 5(d) simply by reversing the sign of the x-gradient field -as illustrated in Fig. 5(e).These results indicate that the wave-packet can become either amplified or damped during the LZB oscillations accompanied with strong secondary emissions, depending on the direction of the driving field (i.e., the applied index gradient as in our photonic lattices) in the non-Hermitian system [50].
We have theoretically and numerically investigated the LZB oscillations in flatband photonic Lieb ribbons with nonzero threshold PT symmetry breaking.
The PT symmetry is introduced by judiciously tailored couplings with asymmetric complex coupling coefficients, and the refractive-index gradient is applied acting as an external field.In such a system, the flatband can survive and a nonzero triply degenerated EP can be established in the presence of y-gradient.Furthermore, the LZB oscillations exhibit intriguing features such as asymmetric energy distribution, pseudo-Hermitian propagation, and secondary emissions at the EP.The proposed artificial lattices can also be realized and other physical platforms, where the presence of flatband in a PT-symmetric phase could lead to intriguing fundamental phenomena when disorder and nonlinearity are introduced.In future works, it could also be interesting to investigate such higher-order EPs in twodimensional systems and non-Hermitian topological transitions associated with them.

FIG. 1 .
FIG. 1.(a) Schematic diagram of a driven Lieb lattice ribbon formed by coupled optical waveguide arrays with each unit cell consisting of five sites (A, B, C, D and E).The coupling coefficients are t, t1 = t (1 + ig) and t2 = t (1 − ig), where g is the non-Hermitian parameter.Colored waveguides show the foursite CLS with equal amplitude but opposite phase of a uniform lattice (i.e., when t1 = t2), ensuring destructive interference in the neighboring A(E) sites.∆βx and ∆βy represent x-and y-gradient, respectively.(b, c) Calculated real and imaginary parts of the spectrum as a function of g when no index gradient is present (i.e., ∆βx = 0 , ∆βy = 0).The red dashed circles represent the position of the EP, and the pink filled circles show the value in which the upmost and lowest dispersive bands become complex.(d) A flatband CLS occupies two plaquettes for g = 0 with colored and filled circles representing the sites of nonzero energies, and the characters right next to them show the corresponding amplitudes.For simplicity, we set t = 1, and the lattice period is 1..

FIG. 4 .
FIG. 4. Flatband LZB oscillations with an x-gradient ∆βx = 0.05t.(a) An initial excitation beam formed by combining the CLSs with a broad Gaussian distribution (red line) bn = −cn = dn = e −n 2 /2σ 2 (σ = 15) at normal incidence.Colored sites represent nonzero amplitudes with blue and green sites having π phase difference.(b) ∆βy = 0, g = 0.005, the input keeps compactness with slight and symmetric oscillations for non-Hermitian lattices in the absence of y-gradient.(c) ∆βy = 0.12t, g = 0.005, the total energy stays constant over many periods and the energy distributes asymmetrically for non-Hermitian lattices in the presence of y-gradient.(d) ∆βy = 0.12t, g = 0.015, the energy of the wave-packet is amplified during oscillations.(e) ∆βy = 0.12t, g = 0, initial state evolves in a symmetric way for Hermitian lattices in the presence of y-gradient for direct comparison.The oscillation period is l = 2π/∆βx and the propagation length L = 5l is divided into 2000 steps.For all figures, the input peak intensity is normalized to 1.

FIG. 5 .
FIG. 5. LZB oscillations under a dispersive band excitation.(a) An initial input beam formed by combining the CLSs with a broad Gaussian distribution (red line) an = bn = cn = dn = en = e −n 2 /2σ 2 (σ = 15) at normal incidence.Colored sites represent nonzero amplitudes with equal phase structure.(b) Oscillations in a Hermitian system with g = 0, showing constant amplitude during propagation.(c, d) Oscillations in a non-Hermitian system with (c) g = 0.015 and (d) g = 0.03, showing secondary emissions and a decaying amplitude.(e) same as in (d) but the sign of x-gradient field is opposite, showing a growing amplitude.Other parameters are the same as in Fig. 4(c).