Non-Hermitian Floquet Phases with Even-Integer Topological Invariants in a Periodically Quenched Two-Leg Ladder

Periodically driven non-Hermitian systems could possess exotic nonequilibrium phases with unique topological, dynamical, and transport properties. In this work, we introduce an experimentally realizable two-leg ladder model subjecting to both time-periodic quenches and non-Hermitian effects, which belongs to an extended CII symmetry class. Due to the interplay between drivings and nonreciprocity, rich non-Hermitian Floquet topological phases emerge in the system, with each of them characterized by a pair of even-integer topological invariants (w0,wπ)∈2Z×2Z. Under the open boundary condition, these invariants further predict the number of zero- and π-quasienergy modes localized around the edges of the system. We finally construct a generalized version of the mean chiral displacement, which could be employed as a dynamical probe to the topological invariants of non-Hermitian Floquet phases in the CII symmetry class. Our work thus introduces a new type of non-Hermitian Floquet topological matter, and further reveals the richness of topology and dynamics in driven open systems.

In previous studies, non-Hermitian Floquet phases were explored mainly in two-band systems. According to the periodic table of topological matter, the CII symmetry class refers to the chiral symplectic class. A topological phase in this symmetry class is protected by the time-reversal, particle-hole, and chiral symmetries [64,65]. Moreover, the tight-binding Hamiltonian describing a system in the CII symmetry class possesses a minimum number of four bands. In this work, we uncover a new type of non-Hermitian Floquet topological matter in the extended CII symmetry class, which contains at least four quasienergy bands. Our system can be realized in a periodically quenched nonreciprocal two-leg ladder, as introduced in Section 2. Each topological phase of the system is characterized by a pair of even-integer winding numbers (w 0 , w π ) ∈ 2Z × 2Z. With the change of the nonreciprocal parameters of the model, we find rich non-Hermitian Floquet phases with large winding numbers, and various topological phase transitions induced by non-Hermitian effects, as presented in Section 3. In Section 4, we obtain multiple quartets of non-Hermitian Floquet edge modes in our system at zero and π quasienergies under the open boundary condition (OBC), and establish the correspondence between the number of these modes and bulk topological invariants (w 0 , w π ). In Section 5, we explore the dynamical aspects of our model by generalizing the mean chiral displacement (MCD) to non-Hermitian Floquet systems in the CII symmetry class, and demonstrate the relationship between the MCDs and topological winding numbers (w 0 , w π ). Finally, we summarize our findings and discuss potential future directions in Section 6.

Model and Symmetry
The model we are going to investigate can be viewed as a driven, non-Hermitian version of the Creutz ladder (CL) with spin-1/2 fermions and spin-orbit couplings (or spinless particles with four sublattice degrees of freedom). The CL model refers to a quasi-one-dimensional lattice formed by two coupled legs and subjected to a perpendicular magnetic flux [66]. It could possess symmetry-protected degenerate zero modes at its boundaries, and therefore belong to one of the earliest examples of a topological insulator [66]. In later studies, the CL model has been realized in photonic [67,68] and cold atom [69,70] systems, and utilized in the investigations of Aharonov-Bohm cages [71,72], topological pumping [73], localization [74,75], and many-body topological matter [76][77][78][79][80]. Recently, spin-1/2 extensions of the CL model have also been explored in several studies [81][82][83], leading to the discoveries of richer topological features. Furthermore, when time-periodic drivings are applied to the spin-1/2 CL, a series of Hermitian Floquet topological phases in the CII symmetry class were found [84]. Each of these phases is characterized by a pair of even-integer topological winding numbers, quantized dynamics of bulk states, together with degenerate quartets of zero and π Floquet edge modes under the OBC [84]. In this work, the construction of our system can be viewed as a non-Hermitian extension of the model studied in Ref. [84], and will be referred to as the non-Hermitian periodically quenched two-leg ladder (PQTLL).
The time-dependent Hamiltonian of the non-Hermitian PQTLL model takes the form: where j ∈ Z and T is the driving period. Within the first (second) half of each driving period, the system is described by the time independent Hamiltonian H (H ⊥ ). In the middle of the driving period, a sudden quench is applied to the system, so that its Hamiltonian is switched from H to H ⊥ . Throughout this work, we will seth = T = 1 as the convention of units. In the lattice representation, the Hamiltonian components H and H ⊥ are explicitly given by: Here n = 1, ..., N are the indices of unit cells, which are arranged along the horizontal (x) direction of the ladder. σ 0 and τ 0 are both 2 × 2 identity matrices. An illustration of the model is presented in Figure 1. Each unit cell of the ladder contains two spin and sublattice components, and σ x,y,z , τ x,y,z are Pauli matrices acting on the spin-1/2 and sublattice degrees of freedom, respectively. The system parameters J x and J y represent the intercell and intracell hopping amplitudes along the horizontal (x) and vertical (y) directions of the ladder. J d controls the coupling strength between different sublattices in adjacent unit cells, and V describes the amplitude of spin-orbit coupling among particles with opposite spins in the same sublattice of nearest-neighbor unit cells. In this work, we allow J y and J d to take complex values, i.e., J y = J r y + iJ i y and J d = J r d + iJ i d . Physically, the imaginary parts J i y and J i d could describe the nonreciprocity of hoppings along the vertical and diagonal directions of the ladder. The Floquet operator of the non-Hermitian PQTLL model, which corresponds to its evolution operator over a complete driving period (e.g., from t = j + 0 − to j + 1 + 0 − ), can be expressed as: where T is the time-ordering operator, which directly leads to the last equality as the system Hamiltonian in Equation (1) is time-independent within each half of the driving period. The quasienergy spectrum ε of the system can be obtained by solving the eigenvalue equation U|ψ = e −iε |ψ under a fixed boundary condition, where |ψ is the corresponding right eigenvector of U. With a ladder of N unit cells and under the period boundary condition (PBC), one can perform the Fourier transform |n = 1 √ N ∑ k e −ink |k to express U in momentum space as U = ∑ k |k U(k) k|. Here k ∈ [−π, π) is the quasimomentum, and Since J y and J d are in general complex parameters of the system, U(k) is not unitary. In terms of the biorthogonal eigenbasis of U(k), Equation (5) can be equivalently written as: where |ε η (k) ( ε η (k)|) is the right (left) eigenvector of U(k) with the quasienergy ε η (k) = ηε (k) ∈ C. = 1, 2 are the indices of the two quasienergy bands, whose real parts satisfy Re[ε (k)] ∈ (0, π]. The complex dispersion {ε η (k)} thus contains four Floquet bands, with two possible spectral gaps at quasienergies zero and π. A topological phase transition may happen when a gap closes at one of these quasienergies. This can be further captured by the vanishing of one of the two gap functions ∆ 0 and ∆ π , defined as: In the next section, these functions will be utilized to obtain the boundaries between different Floquet topological phases of the non-Hermitian PQTLL model. The topological invariants of the non-Hermitian Floquet phases in our system are determined by the symmetries of U(k). Following the usual strategy in the study of one-dimensional (1D) Floquet systems [85,86], we rewrite U(k) in a pair of symmetric time frames as: It is clear that U 1,2 (k) and U(k) are related by similarity transformations, which can be achieved by shifting the initial time of the driving forward or backward over half a period. The Floquet operators U 1,2 (k) thus share the same quasienergy dispersion with U(k), and they can be expressed in their corresponding biorthogonal basis and as: where α = 1, 2 denote the two time frames. Moreover, the effective Hamiltonians h 1,2 (k) in Equations (11) and (12) both possess the extended time-reversal symmetry T , the extended particle-hole symmetry C, and the sublattice (chiral) symmetry S, i.e., According to the symmetry classification of Floquet systems [87] and the periodic table of non-Hermitian topological phases [9,11], the non-Hermitian PQTLL model belongs to an extended CII symmetry class with even-integer topological invariants. In the meantime, the system also possesses the inversion According to Ref. [9], the coexistence of time-reversal and inversion symmetries allows a system to be immune to the non-Hermitian skin effect [88][89][90]. The topological characterization and bulk-boundary correspondence of our non-Hermitian PQTLL model can thus be treated in a standard manner, as will be presented in the following sections. Note in passing that under the combined action of P and T , we have . This is a different form of the P T -symmetry, which in general cannot guarantee the realness of the quasienergy spectrum of U α (k). Besides, even for a P T -symmetric system, the spectrum can only be real in the P T -invariant regime, and will in general be transformed from real to complex when the system parameters are varied across a P T -symmetry breaking transition [91,92]. Neverless, as will be shown in Section 4, the topological edge modes in our system always possess real quasienergies zero or π. These Floquet edge modes are thus immune to any possible P T -breaking transitions so long as the sublattice symmetry S is preserved.

Topological Invariants and Phase Diagrams
In this section, we introduce the topological invariants of our non-Hermitian PQTLL model, and construct its topological phase diagrams for typical situations.
Following the symmetry analysis in the last section and the topological characterizations of Hermitian Floquet phases [84,93], the Floquet operator U α (k) in the α's time frame possesses a topological winding number w α , which can be defined as: where α = 1, 2, k is the quasimomentum, S is the sublattice symmetry operator, and the trace is taken over all the internal degrees of freedom including spins and sublattices. Q α (k) is usually called the Q-matrix [94], which takes the form of a biorthogonal projector: Here = 1, 2 are the indices of the two Floquet quasienergy bands, whose real parts are positive. η = ± denote the signs of the real parts of quasienergies. The right (left) eigenvectors {|ε [94]. Therefore, Q α (k) can be viewed as a flattened effective Hamiltonian of U α (k), whose eigenphases with positive and negative real parts are set to zero and π, respectively. Note that since Q α is given by the difference between the projectors of two sets of bulk quasienergy bands (with η = ±), its formalism does not rely on the exact number of bands possessed by the system [9]. Indeed, through the formalism of the projector in Equation (18), only the net information of bands contained in the two quasienergy ranges Re(ε) ∈ (−π, 0) and Re(ε) ∈ (0, π) are taken into accout, which generalizes the winding number of two-band systems to multiple-band cases [9]. For a system with the sublattice symmetry S and the spectral gaps at ε = 0, π, the total number of bulk bands is even. Furthermore, the set of Pauli matrices σ x,y,z plus the 2 × 2 identity σ 0 could not satisfy all the symmetry requirements of the CII class (i.e., T T * = −1, CC * = −1 and S 2 = 1) simultaneously. Therefore, a 1D Floquet system in the CII symmetry class contains at least four quasienergy bands. The model we introduced satisfies this requirement for the minimal number of bulk bands.
With the help of (w 1 , w 2 ) in Equation (17) and the topological characterization of chiral symmetric Floquet systems [86], we can construct another pair of topological winding numbers (w 0 , w π ) as: According to Ref. [84], these invariants would only take even-integer values, and they provide a complete characterization for all 1D Hermitian Floquet topological phases in the CII symmetry class. Furthermore, the requirement of two invariants reveals the difference between Floquet and non-driven systems. Since the Floquet operator U possesses two quasienergy gaps at ε = 0 and π, there could be two types of degenerate edge modes at these quasienergies, whose numbers are thus counted by two distinct topological invariants.
In the following, we will demonstrate that the topological invariants (w 0 , w π ) in Equation (19) could also characterize the Floquet phases of the non-Hermitian PQTLL model, and they always take real and even-integer values for a gapped topological phase. Note in passing that for our system, the spectral gaps take the form of lines through zero and π quasienergies. For a non-driven system with the sublattice symmetry S, the non-Hermitian topological phases can be characterized by a winding number w ∈ 2Z [9], which count the number of zero-energy edge modes under the open boundary conditions. Our results extend this topological characterization to 2Z × 2Z, with the second even integer w π being related to the line-gap induced at the quasienergy π by the periodic driving fields. By evaluating (w 0 , w π ) numerically with Equations (17) and (18), we obtain the topological phase diagrams of the non-Hermitian PQTLL model for two typical cases, as presented in Figure 2 and 3. In Figure 2, we show the phase diagram of the system with respect to the real and imaginary parts of the vertical hopping amplitude J r y and J i y . The other system parameters are all chosen to be real and set as (J x , J d , V) = (0.5π, 4π, 0.1π). From the phase diagram, we see clearly that with the increase of the nonreciprocal parameter J i y , a series of topological phase transitions can be induced, with each of them being followed by the quantized change of w 0 or w π by an integer multiple of two. The resulting non-Hermitian Floquet topological phases could possess large and even-integer topological invariants due to the interplay between drivings and non-Hermitian effects. Moreover, phases carrying larger topological winding numbers can be realized when the diagonal hopping amplitude J d takes larger values. Therefore, the realization of non-Hermitian PQTLL model could also provide us with a convenient platform to explore non-Hermitian phases with large and even-integer topological numbers.
In Figure 3, we present the topological phase diagram of the non-Hermitian PQTLL model versus the imaginary parts of the vertical and diagonal hopping amplitudes J i y and J i d . The other system parameters are fixed at (J x , J r y , J r d , V) = (0.5π, 0.6π, 4π, 0.1π). From the phase diagram, we again observe rich non-Hermitian Floquet topological phases characterized by (w 0 , w π ) ∈ 2Z × 2Z, and multiple topological phase transitions induced by the change of the two non-Hermitian parameters. Furthermore, in certain regions of the phase diagram (e.g., around J i y = 6), we find phase transitions accompanied by the increase of topological winding numbers (w 0 , w π ) when the value of J i d rises.The emergence of such phases with stronger topological signatures in deeper non-Hermitian regimes (here at larger J i d ) is unique to Floquet non-Hermitian systems. In the meantime, it also suggests an approach to prepare topological phases with large winding numbers under the collaboration of drivings and nonreciprocity.  In the following two sections, we will present the edge states and bulk dynamics of the non-Hermitian PQTLL model, which would provide us with more transparent signatures about its topological properties.

Edge States and Bulk-Edge Correspondence
One of the key features for a gapped topological phase is the presence of degenerate edge states under the OBC [95]. In this section, we demonstrate the existence of Floquet topological edge modes at zero-and π-quasienergies in our non-Hermitian PQTLL model, and relate their numbers to the bulk topological winding numbers (w 0 , w π ) in Equation (19).
The Floquet quasienergy spectrum of our system under the OBC is obtained by solving the quasienergy eigenvalue equation U|ψ = e −iε |ψ , with the Floquet operator U given by Equation (4). With the quasienergy ε, we can define the gap functions under the OBC as: It is clear that ∆ 0 = 0 (∆ π = 0) only when the spectrum gap closes at the quasienergy 0 (π). (∆ 0 , ∆ π ) can thus be used to characterize the behaviors of the Floquet spectrum around the quasienergies ε = 0 and π. In Figure 4a,b, we present the gap functions ∆ 0 (red solid lines) and ∆ π (blue dashed lines) of the non-Hermitian PQTLL model versus the imaginary parts of the vertical and diagonal hopping amplitudes J i y and J i d for two typical sets of system parameters, respectively. In both panels, we see clearly that with the increase of the nonreciprocal hopping parameter J i y or J i d , the system undergoes a series of topological phase transitions. Each transition is accompanied by the closing and reopening of a line gap through the quasienergy zero or π, together with the increase or decrease of the number of Floquet zero or π edge modes by an integer multiple of four, as denoted in the figure. Intriguingly, by enhancing the strength of nonreciprocity, we observe transitions from topological phases with smaller winding numbers (w 0 , w π ) to larger ones with more edge modes in Figure 4b. The physical mechanism behind this interesting observation is again the interplay between drivings and non-Hermitian effects. Besides, it also indicates the possibility of preparing non-Hermitian Floquet topological states with the help of nonreciprocity.
Furthermore, comparing the number of quartets of the zero (π) edge modes n 0 (n π ) and the bulk winding number w 0 (w π ) in each regime of the non-Hermitian Floquet topological phase, we find the following bulk-edge correspondence relations: These relations hold so long as the symmetries that are protecting the non-Hermitian Floquet topological phases of the system are preserved. Experimentally, Equation (21) could also provide us with a window to look into the topological invariants of non-Hermitian Floquet systems in the CII symmetry class by imaging the edge modes. More generally, in the symmetric time frame α (= 1, 2), we can directly define a noncommutative winding number [84,94,96,97] under the OBC as: where S is again the sublattice symmetry operator, andn = ∑ N n=1 n|n n|σ 0 ⊗ τ 0 is the unit-cell position operator of the ladder. The total number of unit cells N contains two parts, i.e., N = N B + 2N E , where N B and N E are the number of unit cells in the bulk (n ∈ [N E + 1, N E + N B ]) and edge (n ∈ [1, N]) regions of the system, and the trace Tr B (·) is only taken over the bulk degrees of freedom. Different from the previous study [84], the Q-matrix for our non-Hermitian Floquet system in the α's time frame and under the OBC is expressed in the biorthogonal basis as: where n = 1, 2, ..., 2N, η = ±, and |ε η αn is the right Floquet eigenvector satisfying U α |ε η αn = e −iε η n |ε η αn . The left eigenvectors can be obtained by writing U α as U α = V α Λ α V −1 α , where Λ α is diagonal and {|ε η αn } correspond to the columns of (V −1 α ) † [94]. Similar to the Q-matrix under the PBC, Q α here can be viewed as an effective Hamiltonian of the Floquet operator U α , whose eigenvalues are set to 1 (−1) if the corresponding quasienergies of U α satisfying Re(ε η n ) > 0 [Re(ε η n ) < 0]. With the help ofw 1 andw 2 in Equation (23), we can construct another pair of topological invariants [84]: In a fixed time frame α, previous studies have shown thatw α = w α [84]. Therefore, we find the following bulk-edge correspondence relations for 1D non-Hermitian Floquet systems in the extended CII symmetry class, i.e., (|w 0 |, |w π |) = (|w 0 |, |w π |) = (2n 0 , 2n π ). (25) Since the second equality holds also under the OBC, it can be employed to investigate the effect of impurity on non-Hermitian Floquet topological phases, and characterize the non-Hermitian Floquet Anderson insulators that may appear due to the interplay between drivings, non-Hermiticity, and disorder. These topics are beyond the scope of the current work, and will be explored in the future. Despite edge states, the topological signatures of non-Hermitian Floquet phases can also be extracted from bulk dynamics, as will be discussed in the next section.   panels (a,b), respectively. The system parameters are (J x , J r y , J d , V) = (0.5π, 1.5π, 4π, 0.1π) for panel (a) and (J x , J y , J r d , V) = (0.5π, 0.6π + 6i, 4π, 0.1π) for panel (b). The number of quartets of zero and π edge modes (n 0 , n π ) are denoted explicitly near ∆ 0 = ∆ π = 0 in both panels, which are related to the winding numbers (w 0 , w π ) through the relations in Equation (21). The ticks along the horizontal axis denote the bulk gap closing points, whose numerical values are (y π 1 , y 0 1 , y 0 2 , y π 2 , y π 3 , y

Dynamical Probe to the Topological Phases
The mean chiral displacement (MCD) refers to the time-averaged chiral displacement Sn of a wavepacket in a lattice, where S is the sublattice symmetry operator andn is the position operator of the unit cell. The MCD was first introduced as a dynamical probe to the winding numbers of 1D topological insulators in the symmetry classes AIII and BDI [98], and later extended to Floquet systems [93,99,100], two-dimensional systems [101], many-body systems [102], systems in other symmetry classes [84], and recently also to non-Hermitian systems [50,51,53]. In the meantime, the MCD has also been measured experimentally in photonic [98,103] and cold atom [104,105] setups. In this section, we further generalize the MCD to non-Hermitian Floquet systems in the CII symmetry class, and employ it to dynamically characterize the topological phases found in the non-Hermitian PQTLL model.
For a non-Hermitian Floquet system with sublattice symmetry S, we define the MCD C α as the stroboscopic long-time average of the chiral displacement operator Sn in a given symmetric time frame α (= 1, 2), i.e., where m counts the number of driving periods, which has been set to 1 following our choice of units. |ψ(m) and |ψ(m) are the initial states evolved over m's driving periods in the right and left Hilbert spaces, respectively. The C α defined in this way is generally a complex number for finite M due to the implemented biorthogonal average. However, we will show that the imaginary part of C α tend to vanish in the long-time limit M → ∞.
Taking the Fourier transform from the position to momentum representation, and choosing the initial state to uniformly fill the non-Hermitian quasienergy band ( , η) ( = 1, 2, η = ±), we find the following form of MCD: Here |ε η α (k) ( ε η α (k)|) is the right (left) quasienergy eigenvector, and the corresponding Floquet operators can be expressed in the biorthogonal basis as: Note that in Equation (27), a normalization factor has been added to cancel the changing norm of the state during the nonunitary evolution. Inserting the identity in biorthogonal basis I = ∑ ,η |ε η α (k) ε η α (k)|, and using the orthonormality between left and right eigenvectors ε η α (k)|ε η α (k) = δ δ ηη , the denominator of c η α (k) in Equation (28) can be simplified as: where Im[ε η (k)] yields the imaginary part of the quasienergy ε η (k). Similarly, the numerator of c η α (k) can be expressed as: where we have also used the fact S|ε η α (k) ∝ |ε −η α (k) to arrive at the second term on the right hand side of Equation (32). Plugging Equations (31) and (32) into Equation (28), we find the integrand c η α (k) to be: The first term on the right hand side of Equation (33) will be related to the winding number of the system in the α's time frame. If Im[ε η (k)] > 0, the second term on the right hand side of Equation (33) will vanish in general after taking the sum over m and the limit M → ∞, as imposed in Equation (27). However, when Im[ε η (k)] < 0, the factor e i2ε η (k)m will contribute an exponentially growing factor to c η α (k) after the summation over m, making it diverge in the limit M → ∞.
To remove the divergence, we introduce another pair of Floquet propagators for the right and left initial states with Im[ε η (k)] < 0, which are given by: Comapring with Equation (29), it is clear thatÚ α (k) is just the inverse of Floquet operator U α (k), which can be obtained by simply reversing the signs of all the system parameters in our model. Note thatÚ α (k) and U α (k) correspond to the Floquet operators in the left and right Hilbert spaces, respectively. They have different expressions in the biorthogonal Floquet eigenbasis, and are therefore distinguished by different accents on their heads. With these considerations, we modify the definition of C η α in Equation (27) to: where c η α (k) is given by Equation (28), andč η α (k) takes the form: Following the steps in the derivations of Equations (31) and (32), we find the denominator and numerator ofč η α (k) to be: ε Plugging them into Equation (37), we further obtain: It is clear that under the condition Im[ε η (k)] < 0, the second term on the RHS of Equation (40) will in general vanish under the summation and long-time average lim M→∞ 1 M ∑ m , as imposed in Equation (36). Next, we extend the initial state of our system to an incoherent summation of all uniformly filled Floquet bands ( , η), which is equivalent to an "infinite-temperature" state at each quasimomentum k. In the lattice representation, such an initial state corresponds to the uniform population of all the four internal degrees of freedom (spins and sublattices) in the central unit cell of the ladder, which is relatively easy to prepare. For such an initial state, the MCD becomes C α = ∑ ,η C η α . With the help of Equations (33), (36) and (40), it can be written more compactly as: . It is not hard to see that the second term on the right hand side of Equation (41) will tend to vanish in the long-time limit M → ∞, so long as ε η (k) = ±π/2, ±π have solutions only at a discrete set of k-points in the first Brillouin zone, which is the case for our system. Finally, the relation between C α and the winding number w α in the α's time frame can be established as follows. Inserting the expression of projector Q α (k) in Equation (18) into the definition of w α in Equation (17), and taking the trace in the biorthogonal basis, we find: Using again the orthonormality between left and right eigenvectors and the fact S|ε η α (k) ∝ |ε −η α (k) , the expression for w α can be simplified to: Comparing Equation (43) with Equation (41), we find the relation between the long-time averaged MCD C α and winding number w α as: Furthermore, with the help of the relations between (w 1 , w 2 ) and the topological invariants (w 0 , w π ) in Equation (19), we arrive at the relations between the MCDs and the topological winding numbers of 1D non-Hermitian Floquet systems in the CII symmetry class, i.e., These relations establish a connection between the topology and dynamics of the non-Hermitian Floquet states in the system. They also provide us with a powerful way to probe the non-Hermitian Floquet topological phases in the CII symmetry class by measuring the MCDs experimentally in a pair of symmetric time frames.
For completeness, we demonstrate the relations in Equation (45) by numerically simulating the dynamics. The results for two typical cases are presented in Figure 5a,b. In both panels, the time average is taken over M = 20 driving periods, which is well within reach in current experiments. It is clear that the MCDs and topological winding numbers are consistent for all the non-Hermitian Floquet topological phases considered in the figure, and the small deviations are mainly originated from the finite-time effect. Furthermore, a quantized jump of the MCD is observed every time when the system passes through a topological phase transition point. Experimentally, the MCDs have been measured in both the cold atom [104,105] and photonic systems [98,103], in which non-Hermiticity and driving fields can also be implemented [1]. Furthermore, the MCDs may also be detected directly in momentum space with the help of a recently proposed setup certaining the nitrogen-vacancy-center in diamond [36]. Putting these together, we conclude that the MCD can indeed be employed as a dynamical probe to the topological phases and phase transitions in our non-Hermitian PQTLL model, and also in other 1D non-Hermitian Floquet systems in the CII symmetry class. (magenta circles), and C π = − C 1 −C 2 2 (black squares) versus the imaginary parts of vertical and diagonal hopping amplitudes J i y and J i d of the non-Hermitian PQTLL model in panels (a,b), respectively. The other system parameters are chosen as (J x , J r y , J d , V) = (0.5π, 1.5π, 4π, 0.1π) for panel (a) and (J x , J y , J r d , V) = (0.5π, 0.6π + 6i, 4π, 0.1π) for panel (b). The MCDs are averaged over M = 20 driving periods for the results in both panels.

Conclusions
In this work, we introduced a periodically quenched two-leg ladder model subjecting to nonreciprocal inter-leg hoppings. The system belongs to an extended CII symmetry class in the non-Hermitian periodic table [9,11], which is further characterized by a pair of even-integer topological winding numbers (w 0 , w π ) ∈ 2Z × 2Z due to the existence of time-periodic drivings. We established the topological phase diagrams of the model, and observed rich non-Hermitian Floquet topological phases with large winding numbers. In particular, Floquet phases carrying larger topological invariants can emerge in stronger non-Hermitian regimes thanks to the collaboration between drivings and non-Hermiticity. Under the open boundary condition, Floquet topological edge modes with zero and π quasienergies appear as fourfold degenerate quartets around the boundaries, whose exact numbers are determined by the bulk topological invariants (w 0 , w π ). Besides the bulk-edge correspondence, we introduced the generalized mean chiral displacement as another probe to the topological features of our system dynamically, and showed that the MCDs in the long-time limit yield the topological invariants of one-dimensional non-Hermitian Floquet systems in the CII symmetry class. Our work not only uncovers a new type of topological phase originated from the interplay between drivings and non-Hermitian effects, but also paves the way for the dynamical characterization of non-Hermitian Floquet topological matter. In future work, it would be interesting to extend our findings to other symmetry classes, higher spatial dimensions, and superconducting systems. Furthermore, intriguing non-Hermitian Floquet phases and phenomena are expected to appear under the effects of disorder and many-body interactions, which certainly deserve thorough explorations.