
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

39232007
70995
10.1038/s41598-024-70995-2
Article
Quasiperiodic disorder induced critical phases in a periodically driven dimerized p-wave Kitaev chain
Roy Koustav koustav.roy@iitg.ac.in

1
Roy Shilpi 12
Basu Saurabh 1
1 https://ror.org/0022nd079 grid.417972.e 0000 0001 1887 8311 Department of Physics, Indian Institute of Technology Guwahati, Guwahati, 781039 Assam India
2 https://ror.org/01tgyzw49 grid.4280.e 0000 0001 2180 6431 Department of Physics, National University of Singapore, 117542 Singapore, Singapore
4 9 2024
4 9 2024
2024
14 206031 3 2024
22 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
The intricate relationship between topology and disorder in non-equilibrium quantum systems presents a captivating avenue for exploring localization phenomenon. Here, we look for a suitable platform that enables an in-depth investigation of the topic. To this end, we delve into the nuanced analysis of the topological and localization characteristics exhibited by a one-dimensional dimerized Kitaev chain under periodic driving and perform detailed analyses of the Floquet Majorana modes. Such a non-equilibrium scenario is made further interesting by including a spatially varying quasiperiodic potential with a temporally modulated amplitude. Apriori, the motivation is to explore an interplay between dimerization and a quasiperiodic disorder in a topological setting which is also known to demonstrate unique (re-entrant) localization properties. While the topological properties of the driven system confirm the presence of zero and π Majorana modes, the phase diagram obtained by constructing a pair of topological invariants (Z×Z), also referred to as the real space winding numbers, at different driving frequencies reveal intriguing features that are distinct from the static scenario. In particular, at either low or intermediate frequency regimes, the phase diagram concerning the zero mode involves two distinct phase transitions, one from a topologically trivial to a non-trivial phase, and another from a topological phase to an Anderson localized phase. On the other hand, the study of the Majorana π mode unveils the emergence of a unique topological phase, characterized by complete localization of both the bulk and the edge modes, which may be called as the Floquet topological Anderson phase. Moreover, different frequency regimes showcase distinct localization features which can be examined via the localization toolbox, namely, the inverse and the normalized participation ratios. Specifically, the low and high-frequency regimes demonstrate the existence of completely extended and localized phases, respectively. While at intermediate frequencies, we observe the critical (multifractal) phase of the model which is further investigated via a finite-size scaling analysis of the fractal dimension. Finally, to add depth into our study, we have performed a mean level spacing analyses and computed the Hausdorff dimension which yields specific characteristics inherent to the critical phase, offering profound insights into its underlying properties.

Subject terms

Phase transitions and critical phenomena
Superconducting properties and materials
Topological matter
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Anderson localization (AL) is a fundamental phenomenon involving the complete vanishing of transport properties of systems due to the presence of random disorder1,2. Consequently, all the single-particle eigenstates of the non-interacting system suffer a transition from a completely extended (metallic) to a localized (insulating) phase. While this transition depends upon the dimensionality of the system3, nevertheless, this captivating topic gained paramount interest in a broad range of physical systems, such as matter waves, light waves, optical lattices, photonic lattices, etc 4–8. On the other hand, quantum systems with incommensurate potential, such as a quasiperiodic (QP) potential, which lies between periodic and random, can exhibit localization transitions in a one-dimensional system9. Recent developments in experimental accessibility to control over quantum systems and engineering new models in need provide a new era of opportunities for both the experimentalists and theoreticians working in the field. As a consequence, quantum systems with QP potential have been studied in a vast range of experimental setups, including optical10–12 and photonic lattices13,14, optical cavities15,16, and moiré lattices17, etc. To understand the localization transition in QP systems, the Aubry-André (AA) model18,19 is the most studied one. Later, several general models of the AA model were introduced, comprising of many interesting non-trivial results12,20–25. Moreover, in specific generalized AA models, the transition from an extended to a localized phase is often linked through a critical region, which is characterized as an intermediate or a critical phase26. Numerous studies have shown that an energy-dependent transition that is a mobility edge may be present in this phase27–29. The emergence of a mobility edge in a one-dimensional system has sparked substantial curiosity owing to its experimental realizations30,31.

On a parallel front, topological quantum matter, such as topological insulators (TIs) and topological superconductors (TSCs), has been receiving enormous attention due to its possible applications in topological quantum computation and spintronics devices32. The scientific community believes the fundamental building blocks of TSCs are the Majorana zero modes (MZMs), which could be a potential candidate for qubits33,34. A simple toy model of the TSCs is provided by the Kitaev chain model35. It is a one-dimensional spinless p-wave superconductor chain with MZMs localized at the edges. Various experimental proposals have been projected so far in order to obtain 1D TSCs36–41. Moreover, several research have been conducted using this model to study the role of QP potential on driving the system towards complete localization42–46. In a similar way, a lot of progress has gained attention from the theoretical aspects. Among them, we are particularly interested in a variant of the generalized Kitaev models. Specifically, a dimerized Kitaev chain47–49, which is a hybrid of a Su-Schriffer-Heeger (SSH) model50 and a Kitaev chain, owing to its exciting features. There are rich topological and localization properties (protected by the particle-hole symmetry) that make the study of this model intriguing51. Furthermore, it is crucial to remember that, as some of us have demonstrated that the interplay between dimerization and the quasiperiodic (QP) potential has been found to have a major impact on a re-entrant localization transition20. This led us to consider such a particular choice of the Hamiltonian which involves both the hopping and superconducting pairing being modulated via a dimerization parameter δ. Other versions of this model that also include dimerization in its spin dual representation may be found in the literature, where conventional ferromagnetic and anti-ferromagnetic phases have been studied using several local and non-local order parameters52–54.Figure 1 Schematic illustration of dimerized Kitaev chain. The red and green circles denote two sublattices, A and B respectively, within a unit cell (dotted line). The thick and thin chain indicates the dimerization of the model. For example, the intracell (strong) and intercell (weak) hopping strengths are defined as t(1+δ) and t(1-δ), respectively. Whereas the intracell (strong) and intercell (weak) superconducting pairing strengths are defined as Δ(1+δ) and Δ(1-δ), respectively.

Further, Quantum systems driven periodically far from equilibrium are known to exhibit new phases that are otherwise inaccessible in a static setup. These kinds of periodically driven systems can be understood by means of Floquet formalism55–57. By using an external periodic drive, one can pave the way for creating topologically non-trivial materials with substantial tunability, even from those that are topologically trivial in the undriven case. Additionally, the energy bands of the driven systems can be backfolded to a Floquet Brillouin zone (FBZ) at the boundary of which, new types of edge modes, namely the so-called π modes, appear58–61. Numerous studies using ultracold atoms in driven optical lattices, acoustic, and photonic devices have successfully employed the concept of Floquet engineering62–68. Among others, silicon-on-insulator-like materials with lattices of tightly linked octagonal resonators have been used in the experimental realization of Floquet topological insulators (FTIs) on nano-photonics platform69. Additionally, photo-induced band gaps can be considered to study the temporal periodicity of systems, where these band gaps can be resolved using a method called the time and angle-resolved photoemission spectroscopy (t-ARPES)70. In recent years, there has been a remarkable surge in interest revealed to the study of non-equilibrium dynamics of closed systems. These include generation of higher winding or the Chern numbers in 1D, quasi-1D, and 2D systems71–79, emergence of time crystalline phases along with period doubling oscillations80,81, Floquet analysis of higher-order topological insulators82,83, Floquet topological characterization of quantum chaos model84,85, etc. In this context, the periodically driven Kitaev chain has been explored in several studies86,87. In addition to that, several studies have been reported on the periodically driven AA model88–91.

Very recently, the interplay of QP potential and topology has been well explored using this model51, where different phases are characterized in terms of topologically trivial and non-trivial regimes. In addition to this, the presence of an entire region comprising of multifractal states makes the model more promising in the context of non-ergodic physics92–94. Moreover, in the presence of random disorder, the model can show behavior similar to the Anderson model49. Further, a recent study95 has shown that the Majorana modes created using spatially QP driving are more robust to decoherence than those created using a spatially uniform one. Deriving motivations from the above inputs, an interesting proposal arose here to bridge the two important aspects, namely the topological properties and the localization-delocalization transition for a time-periodic system, which has been considered by us. Thus, we want to ask very specific questions, such as, how the periodically modulated QP potential affects the above-mentioned properties of the model. Our primary goal here is to compare and contrast the topological and localization properties of a periodically kicked dimerized Kitaev chain corresponding to different frequency regimes. In particular, we observe intriguing nontrivial behavior in the driving setting, which is not present in the static scenario. While the high-frequency regime captures the properties of the static counterpart, the low-frequency regime demands great attention to study deeply.

The layout of the subsequent discussions is as follows. In section “The Hamiltonian and the Floquet formalism”, we describe the static version of the model to recapitulate its properties, including dimerization in both hopping and the superconducting pairing term, and benchmark against the results for the driven case. Afterwards, we shall introduce the Floquet tool to construct an effective time-independent Hamiltonian. In section “Results”, we shall discuss our results on both the topological and the localization properties of the system. Specifically, we have included detailed discussions on the fractal and Hausdorff dimensions and the mean level spacing to distinguish between the localized, critical and extended states. At the end, we summarize and conclude our findings in section “Conclusions”.Figure 2 A static phase diagram is shown using the real-space winding number as a function of dimerization strength (δ) and the onsite QP potential strength (λ). The blue region represents a trivial phase with ν=0, whereas the red region represents a topological phase with ν=1.

Figure 3 The real-space winding number corresponding to the Majorana-zero modes as a function of the driving period (T) and the onsite QP potential strength (λ) is shown. Here, the dimerization strength is δ=0.5. The blue region represents a trivial phase with ν=0, whereas the red region represents a topological phase with ν=1.

The Hamiltonian and the Floquet formalism

We consider a one-dimensional tight-binding model describing the dimerized Kitaev chain of spinless particles with p-wave superconductivity in the presence of onsite QP potential, illustrated in Fig. 1. The corresponding Hamiltonian is given by,1 H=-t∑j=1N[(1+δ)c^j,B†c^j,A+(1-δ)c^j+1,A†c^j,B+H.c.]-Δ∑j=1N[(1+δ)c^j,B†c^j,A†+(1-δ)c^j+1,A†c^j,B†+H.c.]-∑j=1Nμ[c^A,j†c^A,j+c^B,j†c^B,j].

Here, A and B represent sublattice indices. The number of unit cells denoted as N corresponds to the unit cell index j (j=1,2,3,…N). Thus, the length of the chain is given by L=2N. The creation (annihilation) operator to create (annihilate) a fermion at the sublattice site (j, A) and (j, B) is given by c^j,A† (c^j,A) and c^j,B† (c^j,B), respectively. Further, t is the nearest-neighbor hopping integral, and Δ denotes the nearest-neighbor superconducting pairing term, which is taken to be real without any loss of generality. It is assumed that the strengths of the hopping integral (t) and the p-wave superconducting pairing term (Δ) alternate between strong (inside the unit cell) and weak (between the unit cells) bonds. Consequently, a dimerization tuning parameter δ is introduced in the model to distinguish between them. Hence, the intra (inter) cell hopping integral and the superconducting pairing term are represented by t(1+δ) (t(1-δ)) and Δ(1+δ) (Δ(1-δ)), respectively, as shown in Fig. 1. We enforce a restriction on δ to ensure the hopping terms to assume only positive values, namely, |δ|<1. Further, in addition to the onsite term, we modulate the chemical potential at sublattice A (μA) and sublattice B (μB) quasiperiodically, given as, 2a μA=λAcos[2πβ(2j-1)+ϕ],

2b μB=λBcos[2πβ(2j)+ϕ].

The periodicity of the QP potential is described by 1/β, where β is an irrational number usually chosen to be the golden ratio, namely, β=5-12. The phase term of the potential is represented as ϕ, which is set to be zero. The potential strengths at the two sublattice sites are denoted by λA and λB. In this paper, we shall mainly focus on the staggered case, namely, λA=-λB=λ. The staggered case is particularly interesting owing to the possibility of an interplay between the dimerization parameter and strength of the QP potential, etc., on the localization properties20. By considering the limiting cases, the dimerized Kitaev chain Hamiltonian is reduced to the Kitaev chain corresponding to δ=0 and to the SSH chain for μ=0 and Δ=0. Throughout this paper, we have fixed t as a unit of energy, that is, t=1, and the length of the chain is set as L=1220. Additionally, the superconducting pairing strength is set as |Δ|<t, that is, Δ=0.5.Figure 4 The real-space winding number corresponding to the zero and π modes as a function of dimerization strength (δ) and the onsite QP potential strength (λ) are shown in (a) and (b), respectively. Here, the driving frequency is ω=2.5. The system size taken for the calculation is L=1220. The blue region represents a trivial phase with ν0,π=0, whereas the red region represents a topological phase with ν0,π=1.

Figure 5 The Floquet quasi-energy spectrum for a particular frequency, say, ω=2.5 is shown as a function of λ corresponding to (a) δ=0.25 and (b) δ=0.5.

Now, to elucidate the effects of periodic driving, we introduce temporal modulations to the QP potentials, given as,3 μA,B→μA,B∑m=-∞∞δ(t-mT).

Here, the driving protocol (Eq. 3) indicates applying the onsite QP potential at the sublattice sites A and B periodically at times t=mT where m is an integer that counts the number of kicks.

In general, to tackle any periodically driven system, one can adopt the Floquet formalism approach. Which provides a tool to construct an effective time-independent Floquet Hamiltonian whose stroboscopic dynamics are equivalent to the static system. Note that, a detailed analysis of the static system in the presence of QP potential has been studied by some of us in51. According to the Floquet theorem, the dynamical evolution of a periodically kicked system is obtained via a time-ordered product of the Floquet evolution operators. As a result, it can be written as a product of two terms, that is,4 U^(T,0)=T∫0Te-iH(t)dt=U^′U^″

where,5 U^′=e-i∑j=1NμAc^A,j†c^A,j+μBc^B,j†c^B,j

and6 U^″=e-iH^0T.

Here H0 is the static part of the Hamiltonian. We can now numerically diagonalize U^(T,0) to obtain its eigenvectors ψm and eigenvalues e-iEm using,7 U^(T,0)ψm=e-iH^effTψm=e-iEmTψm,

where Heff is the Floquet effective Hamiltonian. Em denotes the quasi-energies which lie within the first FBZ. In the following section, we shall present the exact numerical results obtained by using the Floquet effective Hamiltonian.

Results

Our main aim here is to study the topological and localization properties of the driven system induced by the interplay of the dimerization strength (δ) and the periodic driving amplitude (λ). Hence, based on our numerical analysis, we shall investigate the behavior of the system at different frequencies. In particular, we show different phase transitions and properties of the edge modes and bulk states that are discernible at high- and low-frequency regimes. The system size taken for the numerical calculation is L=1220.

Topological properties

In this section, we shall start our discussion on the topological properties based on the study of zero-energy edge modes, that is, MZMs that emerge in our system. Before focusing on our periodically kicked setting, a brief discussion on the static scenario in the presence of the QP potential is useful and presented in the following.

The existence of p-wave superconductivity in our system implies the particle-hole symmetry is inherently present in the Hamiltonian. Consequently, the quasi-particle energy spectrum is symmetric in nature with respect to the Fermi level (EF=0), even in the presence of an onsite QP potential. Thus, for each particle-like solution with energy +E, there will be a hole-like solution with energy -E. Only the zero-energy states (E=0) are self-conjugated with each other. In this scenario, the topologically non-trivial phase is distinct from the topologically trivial phase by the presence of gapless zero-energy modes. Additionally, this distinction can also be captured via the bulk properties of the Hamiltonian. This yields a topological invariant. Thus, a protected chiral, although with broken translational symmetry (due to the onsite QP potential), hints towards the calculation of the real-space winding number as the topological invariant in our system. In fact, the momentum space formula for the winding number turns out to be useful in this regard. Hence, the real-space winding number (ν) can then be written by drawing an analogy in the momentum space as,51,96,8 ν=1L′Tr(Γ^Q^[Q^,X^]),

where Q^=∑jNjj-j~j~ is obtained by solving for a generic chiral symmetric Hamiltonian Hj=Ejj corresponding to j~=Γ^-1j, where j is the eigenstate index. Γ^ and X^ are operators corresponding to the chiral symmetry and position, respectively. Note that, in the momentum space, the chiral symmetry operators are defined as C^=σ^0⊗σ^z and C^=σ^x⊗σ^0 for the cases, corresponding to μ=0 and μ≠0, respectively. Hence, the real space representation of the chiral symmetry operator (Γ^) can be written as the tensor product of C^ with the corresponding Nth order identity matrix, IN. Finally, Tr denotes the trace over the lattice sites corresponding to half of the length of the chain, namely, L′=L/2, where half the number of sites are considered from the middle of the chain to eliminate edge effects.

The analysis based on the real-space winding number, with the inclusion of a constant chemical potential (μ) in the present scenario results in a distinction between the sublattice and chiral symmetries47. Further, the phase diagram depicting the topological regime via real-space winding numbers in the δ-λ plane is shown in Fig. 2. Subsequently, the onset of the Majorana zero-modes emerges at λ=0, which is used as a benchmark for making a comparison between the static and the driven scenarios. A disorder-free dimerized Kitaev chain shows a phase transition at δ=0.9 corresponding to μ=1.548. On the other hand, introducing the onsite QP potential in the no dimerization limit, that is, δ=0, implies a phase transition from a topologically non-trivial to a trivial at λ<2. While in the strong dimerization limit (δ=1), the system hosts a trivial phase. As the dimerization strength is increased further from δ=0, it helps in exhibiting the phase transition from the topological to the trivial at the lower value of λ with respect to the δ=0 condition.Figure 6 Majorana Zero modes (MZMs) for a particular frequency, say ω=2.5 are shown as a function of λ in (a) and (c) corresponding to δ=0.25 and δ=0.5, respectively. The real-space winding numbers for zero-energy modes are shown in (b) and (d), corresponding to δ=0.25 and δ=0.5, respectively.

Expectedly, owing to the unharmed chiral symmetry and periodic drive, the system can host Majorana zero modes. In addition to that, a Majorana π mode that has no static counterpart can appear at higher values of the driving strength (λ). A study of the topological properties in this system is pursued as per the periodic table of FTIs97. Accordingly, it says that each non-trivial phase of the system can be characterized by a pair of winding numbers (ν0,νπ) ∈ (Z×Z). The classification of the two non-commutative winding numbers relies on the mechanism of building a pair of symmetric time frames73,98,99. Consequently, the Floquet evolution operator acquires a form,9 U^=F^G^,

where, F^ and G^ are related by the chiral symmetry operator (C^) as,10 C^F^C^=G^-1.

It is also easy to verify that if a symmetric time frame corresponding to a Floquet evolution operator, U^1=F^G^ exists, then there must exist another symmetric time frame corresponding to the Floquet operator, U^2=G^F^. Now the Floquet evolution operator in one symmetric time frame from t=T/2 to t=3T/2 reads as,11 U^1=e-iλV^/2e-iH^0Te-iλV^/2,

where,12 λV^=∑j=1NμAc^A,j†c^A,j+μBc^B,j†c^B,j.

Figure 7 The probability distribution of eigenstates as a function of site index (j) are shown corresponding to (a) λ=0.65, (b) λ=0.7, (c) λ=1.9, and (d) λ=2 for a fixed frequency, say, ω=2.5 and dimerization, say, δ=0.25.

Similarly, using Eqs. (9) and (10), the Floquet evolution in the second time frame takes the form,13 U^2=e-iH^0T/2e-iλV^e-iH^0T/2.

Given that U^1 and U^2 are chiral symmetric partners and hence share the same quasi-energy spectrum with that of U^(T,0). Thus, a suitable combination of their winding numbers should be able to provide the topological invariants of the system via the following equations,14 ν0=ν′+ν″2;νπ=ν′-ν″2.

Here, ν′ and ν″ are the winding numbers for the two effective Hamiltonians corresponding to the two symmetric time frames, U^1 and U^2 respectively. Meanwhile, we add some pedagogical details by depicting the procedure for obtaining the bulk invariants in a uniformly driven system as outlined in the supplementary material.

Now, we start our analysis by enumerating the scenario corresponding to different frequency regimes. Moreover, we want to obtain an upper limit for the frequency which will demarcate between low and high frequency regimes. This can be achieved by plotting the winding number corresponding to the Majorana zero mode (ν0) as a function of the time period T and the driving amplitude λ for an arbitrary value of the dimerization strength δ=0.5 in Fig. 3. The phase diagram demonstrates a smooth boundary separating the topological non-trivial phase from a trivial phase in a high-frequency (small period) limit. The observation should emulate the static scenario. On the other hand, in the low-frequency (large period) regime, an intriguing, however highly complex, situation emerges distinct from that of the static counterpart. This non-triviality can be understood by expanding the Hamiltonian using the Baker-Campbell-Hausdorff (BCH) formula given as,15 ln(eXeY)=X+Y+12[X,Y]+112[X-Y,[X,Y]]-124[Y,[X,[X,Y]]]+⋯

Such that the effective Hamiltonian assumes the form,16 H^eff=H^0+λV^T+λT2H^0,V^+λT12H^0,H^0,V^-λ212V^,H^0,V^+⋯

In the high frequency limit (T≪1) and small driving strength (λ), Heff can be truncated upto the first order as,17 H^eff=H^0+λV^T.

Such an expansion tells us that the Hamiltonian corresponding to a small time period (high frequency) regime is equivalent to the static Hamiltonian plus a renormalized potential (denoted by the second term in Eq. 17) that increases linearly with the frequency. The inclusion of the second term along with the static Hamiltonian demarcates the different topological phases via a straight line-like boundary in the λ-T plane, as shown in Fig. 3. However, in the limit of low-frequency (large time period), one can not ignore the effects of the additional nested commutators in Eq. (15) that become more important with increasing power of T. Clearly, in the low-frequency regime, the drive can induce longer-range interactions. As a result, one can get topologically protected zero-energy modes even where the QP potential strength is high. Furthermore, multiple phase transitions induced by the disorder can be observed in the limit of low frequency. For further study in the low-frequency regime and in order to get insightful results, we fix the value of frequency corresponding to this regime. For example, we take ω=2.5.

We show the phase diagrams via the real space winding numbers (ν0 and νπ) in the δ-λ plane corresponding to the above-chosen driving frequency, namely, ω=2.5 in Fig. 4a and b. This value really denotes a representative point in the low-frequency regime, and it hosts a critical phase (See Fig. 9 and associated discussions below). We shall explain both the phase diagrams and associated analysis in the following.

The Majorana zero-energy phase diagram (Fig. 4a) of the driven system is seen to host an interesting topological behavior and can be perceived by comparing it with the static scenario (Fig. 2). Specifically, it shows a trivial phase corresponding to weak values for the dimerization strength (δ) and driving amplitude (λ), which initially was topologically non-trivial in its static counterpart (Fig. 2).Figure 8 The real-space winding number for π-energies is shown as a function of λ in (a). (b) The probability distributions as a function of the site index are shown corresponding to the π modes and the bulk state. The dimerization strength for all the cases is considered as δ=0.5 and the frequency is fixed at ω=2.5.

Here, we choose two representative dimerization values, namely, δ=0.25 and 0.5 to further discuss the behavior of the system. In the first case (δ=0.25), the result demarcates a phase transition from a topologically trivial phase to a topologically non-trivial phase at some critical driving amplitude, (shown in Fig. 4). Later, upon increasing the driving strength, the system undergoes another phase transition through a second critical point, (shown in Fig. 4) where a transition from a topological non-trivial phase to a trivial Anderson phase occurs. Hence, we may conclude that phase transitions induced by period driving are possible, resulting in the emergence of non-trivial phases, and yields a feature that is unseen in the static scenario (Fig. 2). Further, in the second case (δ=0.5), the topological non-trivial phases expands in the large dimerization regime, resulting in a single-phase transition. Therefore, we observe only a single transition from topologically non-trivial to trivial phases at larger dimerization strengths and perceived only at strong disorder.

In order to acquire a concrete understanding of the topological features observed in the phase diagram, we have shown the Floquet quasi-energy spectrum as a function of λ corresponding to two different values of the dimerization strength, say, δ=0.25 and 0.5 in Fig. 5a and b, respectively. The choice of δ=0.25 derives motivation from the left panel of Fig. 4, where an initial trivial phase is driven into a topological phase beyond a certain value of λ, which eventually gets destroyed at larger values of λ. Thus, a topologically non-trivial phase corresponding to the Majorana zero modes appears in the spectrum, implying phase transitions occurring at two values of λ, that are, namely, λ1≃0.65 and λ2≃2.00. These values are shown in Fig. 4a by the intersection of the white dashed line with the red region. On the other hand, corresponding to δ=0.5, we observe the onset of the topological phase from the left edge of Fig. 4, that is, λ=0, implying the emergence of MZMs. The topological phase exists upto λ≃2.2, beyond which the MZMs hybridize with the bulk. Subsequently, the number of MZMs and the winding number (ν0) as a function of λ confirm the validity of the bulk-edge correspondence, which is shown in Fig. 6a and b corresponding to δ=0.25.Figure 9 The average value of the IPR and the NPR are plotted as a function of λ for (a) ω=0.5, (b) ω=2.5, (c) ω=3.5, and (d) ω=12.0.

In addition to that, we show the same for the other value of δ, namely, δ=0.5 in Fig. 6c and d. Note that all these observations validate the phase diagram presented in Fig. 4a.

Next, we have plotted the probability distribution of the zero-energy Floquet eigenstates as a function of the site index around the first transition point, λ1 corresponding to λ=0.6, and 0.7 in Fig. 7a and b, respectively for δ=0.25. A phase transition supported by the extended and localized nature of the zero-energy edge modes precisely indicates the existence of the first transition around the critical point λ1=0.65. On the other hand, we have plotted the same corresponding to other values of λ that lie on either side of the second transition point; namely, we have taken λ=1.9 and 2.0 in Fig. 7c and d, respectively. The complete localization of the states at the edges validates the signature of the MZMs, implying a topological non-trivial phase being present at λ=1.9. With an increase in the driving strength, the probability distribution at λ=2.0 shows a critical behavior at the transition point, thereby indicating a transition to an Anderson phase.

Further, we explore the non-trivial topological features that emerge corresponding to the Majorana π-mode in the system. Expectedly, it is clear from Fig. 4b that the Majorana π-modes appear at a strong driving amplitude λ with respect to the MZMs corresponding to the same values of the parameters as considered for the discussion of the MZMs. Here, we shall study the driving induced features by focusing on a specific value of δ, such as δ=0.5, where we observe a phase transition occurring at λ≃3 (See Fig. 4b). For further clarification of the situation, we plot the winding number corresponding to π-modes in Fig. 8a. The data clearly indicate a phase transition from the topologically trivial to a topologically non-trivial phase around λ≃3. Moreover, we observe the probability distribution of the π-modes and one of the bulk states as a function of the site index in Fig. 8b. Most interestingly, we find the π energies to be localized at the edges along with a completely localized bulk. We have checked that all the bulk states are localized in nature. Therefore, a non-trivial phase comprising of the π energy edge modes in addition to a complete localized bulk state, which may be called the Floquet topological Anderson phase, is observed in this system.Figure 10 A static phase diagram is shown in panel (a) using the parameter η as a function of dimerization strength (δ) and the onsite QP potential strength (λ). Whereas, panel (b) depicts the η phase diagram corresponding to a driven scenario (ω=2.5). Extended, intermediate, and localized phases are shown via ‘EP’, ‘IP’ and ‘LP’ respectively.

Localization properties

In this section, we shall explore the localization properties of bulk states in the periodically kicked setting. To do that, we employ two diagnostic tools, namely, the inverse participation ratio (IPR) and the normalized participation ratio (NPR), to distinguish between the extended, critical, and localized phases. Similar to the non-driven case, the effective Floquet Hamiltonian can be solved using the Bogoliubov transformation. This can be done by defining a quasiparticle operator in terms of the superposition of the single-particle creation (c†) and annihilation (c) operators via,18 Φn†=∑j=1,α=A,BN[uj,α(n)c^j,α†+vj,α(n)c^j,α],

where α and n denote the sublattice and the energy band indices, respectively. While un and vn denote the particle and hole coefficients. Hence, we can define the IPR and NPR corresponding to the nth Floquet eigenstate using ujn and vjn as,2019 IPR(n)=∑j=1,α=A,BN[|uj,α(n)|4+|vj,α(n)|4]

and,20 NPR(n)=[L∑j=1,α=A,BN[|uj,α(n)|4+|vj,α(n)|4]]-1.

In the thermodynamic limit, the IPR(NPR) values scale with the system size as L-1(L0) corresponding to an extended state. On the other hand, it varies as L0(L-1) for a localized state. Moreover, one can calculate the average values of the IPR and NPR averaged over all the Floquet eigenstates, given by,21 ⟨IPR⟩=1L∑n=1LIPR(n)and⟨NPR⟩=1L∑n=1LNPR(n).

In Fig. 9 we have shown the variation of ⟨IPR⟩ and ⟨NPR⟩ as a function of λ corresponding to different frequencies, such as ω=0.5,2.5,3.5, and 12.0. We found all the eigenstates to be extended in nature corresponding to the low-frequency limit, that is, for ω=0.5 in Fig. 9a. Also, all the states are localized for any value of λ at a large frequency, namely, ω=12, shown in Fig. 9d. The situation can be understood by the BCH expanded the Hamiltonian given in Eq. (16). Now, in the limit of high frequency, the renormalized potential dominates over the rest of the parameters in the static counterpart. As a result, we find that at relatively high (low) driving frequencies, all the states become localized (extended). Whereas, at some intermediate frequencies, such as ω=2.5 and 3.5, we find the emergence of all three phases, namely extended, critical, or the multifractal, and localized phases with increasing value of λ. Also, once the system is in the localized phase, it remains localized and is independent of the choice of λ. However, the multifractal nature of the states can not be uniquely determined via IPR or NPR alone. Additionally, we are interested in the global properties of the model. To this end, we shall introduce another quantity η, which is defined as20,26,22 η=log10⟨IPR⟩×⟨NPR⟩.

If either of the phases, such as an extended or a localized phase, is present in the system, the condition η≤-log10(L) arises from the fact that the value survives owing to the finiteness of either of the average values of IPR or NPR. Consequently, an intermediate phase which can be thought of as an admixture of localized and extended states can emerge when both ⟨IPR⟩ and ⟨NPR⟩ are finite. This leads to the condition, η≥-log10(L). Therefore, in our case, where we have considered L=1220, η≥-3.08 would denote an intermediate phase.Figure 11 The fractal dimension D2 is plotted as function of n/L corresponding to (a) λ=0.5, (b) λ=1.5 and (c) λ=3.5. The dimerization strength is taken as δ=0.5. The system sizes chosen for this calculation are L=754,1220, and 1974.

In Fig. 10a and b, we have shown the phase diagram using η value as a function of δ and λ corresponding to the non-driven and the driven cases with the same frequency ω=2.5 for the sake of completeness. The phase diagram contains overall information on the localized, intermediate, and extended phases that appear in the model. By comparing the driven system with the undriven one, it is clear that the effective potential strength decreases in the latter situation. The plot indicates a broader region for the intermediate phase (see Fig. 10b).

Furthermore, it is important to note that while η gives a global picture differentiating distinct phases, it can not uniquely separate distinct states. As a result, a finite-size analysis is required (shown in Fig. 11a–c) the quantity, namely, the fractal dimension, D2, which helps in identifying different states accurately and is defined as23,23 D2=-limL→∞log(IPR)log(L).

It has a value of 1(0) corresponding to the extended(localized) states, in the thermodynamic limit. The critical (multifractal) states have a value between 0 and 1. In Fig. 11a, we plot D2 as a function of the ratio of eigenstate index (n) and system size (L) corresponding to λ=0.5. The D2 values move towards a value 1 as the system size is increased, thereby signaling the onset of a completely extended bulk state. On the other hand, in Fig. 11b the presence of the multifractal states, in addition to both the localized and the delocalized states in the intermediate phase (corresponding to λ=1.5) confirm the emergence of multiple mobility edges. For instance, in addition to the mobility edge that separates critical and localized states, there are also mobility edges that distinguish between extended and localized states, as well as those that separate extended and critical states. Recent study100 has also proven that all three phases namely, extended, critical, and localized phases can co-exist in the same parameter space. In the end, at a higher value of the driving amplitude, namely for λ=3.5, we observe a complete localization to occur. In all of these three cases, we have chosen δ=0.5 and have considered half of the quasi-energy spectrum owing to the presence of the chiral symmetry.

Further, to discern the localization-delocalization transition in the energy spectrum, we chose to look at another numerically accessible quantity that clearly highlights the difference between these distinct phases and hence should widen our analysis of the results. Specifically, we focus on the spectral statistics of the consecutive energy levels89,101,102, which can be captured by the ratio between the adjacent energy gaps, namely, rn which is defined as,24 rn=min{sn,sn-1}max{sn,sn-1}.

Here, sn=En+1-En, En being the nth eigenvalue of the Floquet quasi-energy spectrum. In the localized phase, states that are nearby in energy do not interact or show level repulsion, and as a result nearby energy levels are Poisson distributed with average of rn over all the eigenvalues have ⟨rn⟩∼0.39102. Whereas, in the extended phase ⟨rn⟩ vanishes.

In Fig. 12, we have plotted ⟨rn⟩ as a function of disorder strength, λ for different frequency regimes. In the high-frequency limit (T<1), there is an abrupt transition to a completely localized phase, indicating that in the high-frequency limit, the system has no extended states. Such a localization behavior at high frequencies validates the effect of the renormalized potential obtained via BCH expansion in Eq. (16). While, in the low-frequency limit (T=2.5), such abrupt transition ceases to exist due to the presence of fractal states. For example in the range 0.5<λ<1.5, the rise of ⟨rn⟩ is more gradual and hence is an indicative of the presence of critical phase. Moreover, for λ≤0.5, ⟨rn⟩≈0, indicating the existence of a completely extended phase at low values of disorder.

We notice that ⟨rn⟩≈0 occurs due to the presence of double degeneracy of the eigenvalues in the extended region. To elucidate this, we introduce the concept of even-odd and odd-even level spacings of the quasi-energy spectrum, which we defined by103,25 sne-o=E2n-E2n-1;sno-e=E2n+1-E2n.

Figure 12 The figure shows the variation of ⟨rn⟩ as a function of λ for different frequency regimes that are, T=0.5 (red dashed line), and T=2.5 (green dashed line). the vertical black dotted lines mark the localization transition point. Other parameters are chosen as, δ=0.5 and L=1220.

Figure 13 The even-odd and odd-even level spacings as a function of eigenstate index (n/L) are shown in (a): λ=0.5 (extended phase), (b): λ=1.5 (intermediate phase), and (c): λ=3.5 (localized phase) respectively. Other parameters are chosen as, ω=2.5, δ=0.5 and L=1220.

In Fig. 13, we have plotted the even-odd and odd-even level spacings as a function of site index for a few values of the disorder, namely, λ=0.5,1.5,3.5, corresponding to extended, critical, and the localized phases, respectively. In the extended regime (Fig. 13a), the spectral features exhibit double degeneracy, resulting in a noticeable gap between sne-o and sno-e. In contrast, within the localized regime (Fig. 13c), both sne-o and sno-e show very similar behavior, causing the spectral gap to disappear. Finally, in the intermediate/critical region (Fig. 13b), distributions of sne-o and sno-e display significant fluctuations throughout the entire energy spectrum. Moreover, a small gap can be observed between sne-o and sno-e at the lower end of the spectrum, indicating that the critical phase can be thought of as an admixture of extended, and localized states.

Further, to unveil the distinct characteristics of the states in the critical regime, we utilize the concept of Hausdorff dimension21,104. Employing a direct box-counting method for this analysis, we examine the power-law behavior of the total number of boxes (total number of filled states), NB for a given box length, lB defined via,26 NB∝lB-DH

where DH denotes the Hausdorff dimension corresponding to the energy spectrum. In our case, we compute DH in two different critical phases corresponding to λ=1.2 and λ=1.8 respectively. The characteristics of non-trivial fractal can be captured by 0<DH<1. In Fig. 14, we have plotted NB as a function of box length lB. We find a linear scaling in the log-log scale, consistent with Eq. (26), and the Hausdorff dimensions are obtained as, DH=0.71 and DH=0.80, corresponding to λ=1.2 and λ=1.8 respectively. Both the values are notably less than the geometrical dimension, d=1. Additionally, we observe a significant difference between the Hausdorff dimension determined in our study and that obtained in previous works at the critical point of the AA model, that is DH∼0.5105. This distinction can be attributed to the nature of the eigenstates of the AA model, which is characterized by a sharp transition. Further, we get a variety of values for DH (we have shown two here) which implies different admixtures of the localized and the extended states present in the critical regime. Also, that leads to higher values of Hausdorff dimension compared to that found in the context of AA model. Therefore, we infer that the dimerized Kitaev chain exhibits a denser spectrum compared to that of the AA model.Figure 14 The figure shows the variation of NB as function of lB in the log-log scale corresponding to λ=1.2 (blue dots) and λ=1.8 (red dots). The slopes of these plots give the Hausdorff dimensions, which are obtained as DH=0.71 and 0.80 for λ=1.2 and 1.8, respectively. Other parameters are chosen as, ω=2.5, δ=0.5 and L=1220.

Summarizing, the interplay among dimerization (δ), QP potential (λ), and the drive (T) generates distinct phases within our system. To gain a comprehensive understanding of the origin of these phases and their transitions, it is crucial to simultaneously analyze the role of each parameter in the concurrent emergence of topological and localization properties in the system. To begin with, in the static (undriven) case, the system may host topological properties, accompanied by a delocalization-localization transition induced via a strong QP potential. Further, as the value of dimerization increases from zero (δ=0) to an extreme dimer limit (δ=1), the chain breaks up into isolated dimers. Thus, the following physical picture emerges: since the unit cells are “not talking to each other”, the QP disorder becomes irrelevant. Hence, as the dimerization becomes large, the localization phenomena (with no topological character) set up with ease, even at lower values of the on-site QP strength.

On the other hand, the inclusion of a periodic drive takes the system one step further by notably expediting the localization phenomena as depicted in Fig. 9. Moreover, the periodic drive renormalizes the QP strength (Eq. 17), which significantly influences the delocalization-localization transition. For instance, all the states are delocalized corresponding to a large time period, T (low-frequency regime) of the drive, while a complete localization of states occurs corresponding to a small period (high-frequency regime). The findings of the low-frequency regime imply the absence of topological characteristics, which may be attributed to the existence of a fully delocalized phase. As the frequency is increased, the potential gets renormalized, which induces some of the states to emerge as localized states. This yields an intermediate phase (with some localized and others delocalized) that holds the key to the delocalization-localization transition. Consequently, this phenomenon could offer prospects of inducing a transition from a trivial to a topological phase as a function of drive frequency. Finally, at very high frequencies, all the states become localized with no topological character. Hence, it is evident that in the driven scenario, the drive frequency, ω, plays a role akin to dimerization, δ, in the static system. It is also intriguing to observe that in the intermediate frequency range, a critical/intermediate phase emerges, characterized by the complete localization of low-energy states alongside the delocalization of high-energy states (Fig. 11). This indicates the emergence of an energy-dependent boundary between them, referred to as the mobility edge. As λ increases, the higher-energy states start getting localized. Additionally, since the π energy modes appear at the edges of the Floquet Brillouin zone, it explains why we observe π energy modes only when the disorder, λ, is very high.

Moreover, to attain a thorough comprehension of the diverse phases emerging within our system and the impact of different components of the model on these phases, we have provided a table (Table 1) that distinctly delineates the various topological and localization characteristics in small and large λ,T, and δ regimes.

Further, a comparison between our model and several others that exist in the literature will add insights on the roles played by different parameters that derive our model Hamiltonian. To be specific, our system comprises of the following key elements, namely, periodic driving (T), dimerization (δ), on-site quasiperiodic potential (λ), and the p-wave superconducting pairing term (Δ). By carefully selecting these parameters, the Hamiltonian can be effectively mapped onto several established models, yielding significant and noteworthy physically realizable outcomes. The associated models under specific conditions and their relevant references are listed below.The pure Aubry-André (AA) model (λA=λB=λ, Δ=0, and δ=0)19.

The SSH model (λ=0, and Δ=0)50.

Disordered SSH model (λA=-λB=λ, and Δ=0)20.

The Kitaev chain (λ=0, and δ=0)35.

The dimerized Kitaev chain (λ=0)47.

We commence our exploration with the pure 1D Aubry-André (AA) model renowned for its manifestation of a delocalization-localization transition at the critical transition point, λ=2t19. However, upon introducing p-wave superconducting pairing (Δ) to the AA model, the scenario is significantly altered. Notably, the critical transition shifts to a new value, λ=2(t+Δ), implying that achieving the localization transition necessitates a higher value of λ, thereby delaying the transition43. On the other hand, when we only have a dimerization (δ) parameter, the Hamiltonian reduces to the celebrated Su-Schrieffer-Heeger (SSH) model, known for its topological-trivial phase transition occurring at a critical dimerization strength. Further, upon introducing quasiperiodic on-site disorder with the initial condition being λA=-λB=λ yields a recently discovered phenomenon, namely, the re-entrant localization transition20. Now, if we take a re-look at the model under consideration, δ=0 implies a tight binding chain while, δ=1 breaks the lattice into isolated dimers resulting in a completely localized phase. Thus, despite the existence of non-zero pairing (Δ), the localization process, which initially required a higher λ in the absence of dimerization in the model43, is now expedited with the inclusion of δ, and even manifests at weak values of λ. Finally, the incorporation of time-periodic modulation manifests that the driving frequency plays a role similar to the dimerization in the static case.Table 1 Table depicting various topological and localization phases corresponding to static (left) and dynamic (right) version of the model.

Static system	
	Strong λ	Weak λ	
Strong δ	Trivial and

localized phase

	Trivial and

localized phase

	
Weak δ	Trivial and

localized phase

	Topological and

extended phase

	
Driven system	
	Strong λ	Weak λ	
High T	Trivial and

extended phase

	Trivial and

extended phase

	
Low T	Trivial and

localized phase

	Topological and

localized phase

	
Note that strong(weak) δ refers to δ∼1(0). As δ increases, the span of zero-energy modes decreases with respect to λ. In the extreme dimer limit (δ=1), all zero-energy modes vanish, and all bulk states become localized. Further, in a static system, the critical phase is only present within the range 0<δ<1, where the system undergoes various delocalization-localization transitions. Whereas, for a driven system the critical phase only emerges within the intermediate frequency regime. Additionally, it is noteworthy that in the driven system, the frequency, ω plays a role akin to that of dimerization, δ in the static system.

Conclusions

In this paper, we consider a one-dimensional dimerized Kitaev chain under a periodic drive of quasiperiodically modulated onsite chemical potential. We have analyzed the topological and localization properties due to the interplay of the periodic drive and the dimerization term present in the system. Based on the behavior of Majorana zero and π modes, we have found that low driving frequency induces intriguing observations as compared to the static (undriven) counterpart. Corresponding to the Majorana zero mode, driving induces a phase transition from a trivial to a topologically non-trivial phase, followed by another transition from topological to an Anderson localized phase, which is found for a specific range of the dimerization strength. Most interestingly, a phase, namely the Floquet topological Anderson phase, that consists of the localized π-modes at the edges associated with completely localized bulk states is found at large values of the driving strength. Further, the localization properties of the bulk states are analyzed. The observation indicates a fully extended phase corresponding to a low-frequency range, while a complete localization of all the states is established at the high-frequency regime. However, localized, critical, and extended phases co-exist at an intermediate frequency region. Further, we have established our results via a finite-size scaling analysis of the fractal dimension. Also, we have looked at the mean level spacing corresponding to low and high-frequency regimes, where distinct features are noted for the localized and the critical phases. Finally, the study of Hausdorff dimension elucidates distinct attributes unique to the critical phase, providing deeper insights into its properties.

Our insights that a static system in presence of QP potential simulates rich delocalization-localization phenomenon under non-equilibrium settings allow us to suggest a route for the experimental validation of our results. Interestingly, owing to the intrinsic simplicity of circuit design, significant advancements have been made in topolectrics, affirming the presence of Majorana modes106–108. Further, a dimerized Kitaev chain can be constructed by considering two arrays of inductors within an LC network such that each sublattice is simulated by an inductor109. The top array can then be contemplated as particles and the bottom array as antiparticles. Next, the inclusion of periodic drive as outlined in Ref.110, entails coupling each harmonic of the drive (also termed as Floquet replicas), to the primary circuit via capacitive connections, that needs to match the strength of the drive. Finally, to analyze the circuit’s localization behavior, one may opt for the simplest experimentally measurable parameter, the two-port impedance factor as illustrated in Ref.108.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-70995-2.

Author contributions

All the authors have contributed equally to the manuscript.

Data availability

All data and codes that support the findings of this study can be obtained from the corresponding author upon request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Anderson PW Absence of diffusion in certain random lattices Phys. Rev. 1958 109 1492 10.1103/PhysRev.109.1492
Anderson, P. W. Absence of diffusion in certain random lattices. Phys. Rev. 109, 1492. 10.1103/PhysRev.109.1492 (1958).10.1103/PhysRev.109.1492
2. Mott N The mobility edge since 1967 J. Phys. C: Solid State Phys. 1987 20 3075 10.1088/0022-3719/20/21/008
Mott, N. The mobility edge since 1967. J. Phys. C: Solid State Phys. 20, 3075. 10.1088/0022-3719/20/21/008 (1987).10.1088/0022-3719/20/21/008
3. Abrahams E Anderson PW Licciardello DC Ramakrishnan TV Scaling theory of localization: Absence of quantum diffusion in two dimensions Phys. Rev. Lett. 1979 42 673 10.1103/PhysRevLett.42.673
Abrahams, E., Anderson, P. W., Licciardello, D. C. & Ramakrishnan, T. V. Scaling theory of localization: Absence of quantum diffusion in two dimensions. Phys. Rev. Lett. 42, 673. 10.1103/PhysRevLett.42.673 (1979).10.1103/PhysRevLett.42.673
4. Billy J Direct observation of Anderson localization of matter waves in a controlled disorder Nature 2008 453 891 894 10.1038/nature07000 18548065
Billy, J. et al. Direct observation of Anderson localization of matter waves in a controlled disorder. Nature 453, 891–894. 10.1038/nature07000 (2008).18548065 10.1038/nature07000
5. Wiersma DS Bartolini P Lagendijk A Righini R Localization of light in a disordered medium Nature 1997 390 671 673 10.1038/37757
Wiersma, D. S., Bartolini, P., Lagendijk, A. & Righini, R. Localization of light in a disordered medium. Nature 390, 671–673. 10.1038/37757 (1997).10.1038/37757
6. Lahini Y Avidan A Pozzi F Sorel M Morandotti R Christodoulides DN Silberberg Y Anderson localization and nonlinearity in one-dimensional disordered photonic lattices Phys. Rev. Lett. 2008 100 013906 10.1103/PhysRevLett.100.013906 18232768
Lahini, Y. et al. Anderson localization and nonlinearity in one-dimensional disordered photonic lattices. Phys. Rev. Lett. 100, 013906. 10.1103/PhysRevLett.100.013906 (2008).18232768 10.1103/PhysRevLett.100.013906
7. Drese K Holthaus M Exploring a metal-insulator transition with ultracold atoms in standing light waves Phys. Rev. Lett. 1997 78 2932 10.1103/PhysRevLett.78.2932
Drese, K. & Holthaus, M. Exploring a metal-insulator transition with ultracold atoms in standing light waves. Phys. Rev. Lett. 78, 2932. 10.1103/PhysRevLett.78.2932 (1997).10.1103/PhysRevLett.78.2932
8. Schwartz T Bartal G Fishman S Segev M Transport and Anderson localization in disordered two-dimensional photonic lattices Nature 2007 446 52 55 10.1038/nature05623 17330037
Schwartz, T., Bartal, G., Fishman, S. & Segev, M. Transport and Anderson localization in disordered two-dimensional photonic lattices. Nature 446, 52–55. 10.1038/nature05623 (2007).17330037 10.1038/nature05623
9. Sokoloff J Unusual band structure, wave functions and electrical conductance in crystals with incommensurate periodic potentials Phys. Rep. 1985 126 189 10.1016/0370-1573(85)90088-2
Sokoloff, J. Unusual band structure, wave functions and electrical conductance in crystals with incommensurate periodic potentials. Phys. Rep. 126, 189. 10.1016/0370-1573(85)90088-2 (1985).10.1016/0370-1573(85)90088-2
10. Roati G D’Errico C Fallani L Fattori M Fort C Zaccanti M Modugno G Modugno M Inguscio M Anderson localization of a non-interacting Bose-Einstein condensate Nat. (Lond.) 2008 453 895 10.1038/nature07071 18548066
Roati, G. et al. Anderson localization of a non-interacting Bose-Einstein condensate. Nat. (Lond.) 453, 895. 10.1038/nature07071 (2008).18548066 10.1038/nature07071
11. Lüschen HP Scherg S Kohlert T Schreiber M Bordia P Li X Sarma SD Bloch I Single-particle mobility edge in a one-dimensional quasiperiodic optical lattice Phys. Rev. Lett. 2018 120 160404 10.1103/PhysRevLett.120.160404 29756950
Lüschen, H. P. et al. Single-particle mobility edge in a one-dimensional quasiperiodic optical lattice. Phys. Rev. Lett. 120, 160404. 10.1103/PhysRevLett.120.160404 (2018).29756950 10.1103/PhysRevLett.120.160404
12. Sanchez-Palencia L Santos L Bose-Einstein condensates in optical quasicrystal lattices Phys. Rev. A 2005 72 053607 10.1103/PhysRevA.72.053607
Sanchez-Palencia, L. & Santos, L. Bose-Einstein condensates in optical quasicrystal lattices. Phys. Rev. A 72, 053607. 10.1103/PhysRevA.72.053607 (2005).10.1103/PhysRevA.72.053607
13. Lahini Y Pugatch R Pozzi F Sorel M Morandotti R Davidson N Silberberg Y Observation of a localization transition in quasiperiodic photonic lattices Phys. Rev. Lett. 2009 103 013901 10.1103/PhysRevLett.103.013901 19659147
Lahini, Y. et al. Observation of a localization transition in quasiperiodic photonic lattices. Phys. Rev. Lett. 103, 013901. 10.1103/PhysRevLett.103.013901 (2009).19659147 10.1103/PhysRevLett.103.013901
14. Verbin M Zilberberg O Lahini Y Kraus YE Silberberg Y Topological pumping over a photonic Fibonacci quasicrystal Phys. Rev. B 2015 91 064201 10.1103/PhysRevB.91.064201
Verbin, M., Zilberberg, O., Lahini, Y., Kraus, Y. E. & Silberberg, Y. Topological pumping over a photonic Fibonacci quasicrystal. Phys. Rev. B 91, 064201. 10.1103/PhysRevB.91.064201 (2015).10.1103/PhysRevB.91.064201
15. Tanese D Gurevich E Baboux F Jacqmin T Lemaître A Galopin E Sagnes I Amo A Bloch J Akkermans E Fractal energy spectrum of a polariton gas in a Fibonacci quasiperiodic potential Phys. Rev. Lett. 2014 112 146404 10.1103/PhysRevLett.112.146404 24765996
Tanese, D. et al. Fractal energy spectrum of a polariton gas in a Fibonacci quasiperiodic potential. Phys. Rev. Lett. 112, 146404. 10.1103/PhysRevLett.112.146404 (2014).24765996 10.1103/PhysRevLett.112.146404
16. Goblot V Štrkalj A Pernet N Lado JL Dorow C Lemaître A Gratiet L Harouri A Sagnes I Ravets S Emergence of criticality through a cascade of delocalization transitions in quasiperiodic chains Nat. Phys. 2020 16 832 10.1038/s41567-020-0908-7
Goblot, V. et al. Emergence of criticality through a cascade of delocalization transitions in quasiperiodic chains. Nat. Phys. 16, 832. 10.1038/s41567-020-0908-7 (2020).10.1038/s41567-020-0908-7
17. Balents L Dean CR Efetov DK Young AF Superconductivity and strong correlations in moiré flat bands Nat. Phys. 2020 16 725 10.1038/s41567-020-0906-9
Balents, L., Dean, C. R., Efetov, D. K. & Young, A. F. Superconductivity and strong correlations in moiré flat bands. Nat. Phys. 16, 725. 10.1038/s41567-020-0906-9 (2020).10.1038/s41567-020-0906-9
18. Harper PG Single band motion of conduction electrons in a uniform magnetic field Phys. Soc. A 1955 68 874 10.1088/0370-1298/68/10/304
Harper, P. G. Single band motion of conduction electrons in a uniform magnetic field. Phys. Soc. A 68, 874. 10.1088/0370-1298/68/10/304 (1955).10.1088/0370-1298/68/10/304
19. Aubry S André G Analyticity breaking and Anderson localization in incommensurate lattices Ann. Israel Phys. Soc. 1980 3 133
Aubry, S. & André, G. Analyticity breaking and Anderson localization in incommensurate lattices. Ann. Israel Phys. Soc. 3, 133 (1980).
20. Roy S Mishra T Tanatar B Basu S Reentrant localization transition in a quasiperiodic chain Phys. Rev. Lett. 2021 126 106803 10.1103/PhysRevLett.126.106803 33784113
Roy, S., Mishra, T., Tanatar, B. & Basu, S. Reentrant localization transition in a quasiperiodic chain. Phys. Rev. Lett. 126, 106803. 10.1103/PhysRevLett.126.106803 (2021).33784113 10.1103/PhysRevLett.126.106803
21. Roy S Chattopadhyay S Mishra T Basu S Critical analysis of the reentrant localization transition in a one-dimensional dimerized quasiperiodic lattice Phys. Rev. B 2022 105 214203 10.1103/PhysRevB.105.214203
Roy, S., Chattopadhyay, S., Mishra, T. & Basu, S. Critical analysis of the reentrant localization transition in a one-dimensional dimerized quasiperiodic lattice. Phys. Rev. B 105, 214203. 10.1103/PhysRevB.105.214203 (2022).10.1103/PhysRevB.105.214203
22. Pradhan A Giri MK Mondal S Mishra T Emergence of multiple localization transitions in a one-dimensional quasiperiodic lattice Phys. Rev. B 2022 105 L220201 10.1103/PhysRevB.105.L220201
Pradhan, A., Giri, M. K., Mondal, S. & Mishra, T. Emergence of multiple localization transitions in a one-dimensional quasiperiodic lattice. Phys. Rev. B 105, L220201. 10.1103/PhysRevB.105.L220201 (2022).10.1103/PhysRevB.105.L220201
23. Borgnia DS Vishwanath A Slager R-J Rational approximations of quasiperiodicity via projected Green’s functions Phys. Rev. B 2022 106 054204 10.1103/PhysRevB.106.054204
Borgnia, D. S., Vishwanath, A. & Slager, R.-J. Rational approximations of quasiperiodicity via projected Green’s functions. Phys. Rev. B 106, 054204. 10.1103/PhysRevB.106.054204 (2022).10.1103/PhysRevB.106.054204
24. Szabó A Schneider U Non-power-law universality in one-dimensional quasicrystals Phys. Rev. B 2018 98 134201 10.1103/PhysRevB.98.134201
Szabó, A. & Schneider, U. Non-power-law universality in one-dimensional quasicrystals. Phys. Rev. B 98, 134201. 10.1103/PhysRevB.98.134201 (2018).10.1103/PhysRevB.98.134201
25. Purkayastha A Dhar A Kulkarni M Nonequilibrium phase diagram of a one-dimensional quasiperiodic system with a single-particle mobility edge Phys. Rev. B 2017 96 180204(R) 10.1103/PhysRevB.96.180204
Purkayastha, A., Dhar, A. & Kulkarni, M. Nonequilibrium phase diagram of a one-dimensional quasiperiodic system with a single-particle mobility edge. Phys. Rev. B 96, 180204(R). 10.1103/PhysRevB.96.180204 (2017).10.1103/PhysRevB.96.180204
26. Li X Sarma S Das, Mobility edge and intermediate phase in one-dimensional incommensurate lattice potentials Phys. Rev. B 2020 101 064203 10.1103/PhysRevB.101.064203
Li, X. & Sarma, S. Das, Mobility edge and intermediate phase in one-dimensional incommensurate lattice potentials. Phys. Rev. B 101, 064203. 10.1103/PhysRevB.101.064203 (2020).10.1103/PhysRevB.101.064203
27. Liu F Ghosh S Chong YD Localization and adiabatic pumping in a generalized Aubry-André-Harper model Phys. Rev. B 2015 91 014108 10.1103/PhysRevB.91.014108
Liu, F., Ghosh, S. & Chong, Y. D. Localization and adiabatic pumping in a generalized Aubry-André-Harper model. Phys. Rev. B 91, 014108. 10.1103/PhysRevB.91.014108 (2015).10.1103/PhysRevB.91.014108
28. Wang Y Xia X Zhang L Yao H Chen S You J Zhou Q Liu X-J One-dimensional quasiperiodic Mosaic lattice with exact mobility edges Phys. Rev. Lett. 2020 125 196604 10.1103/PhysRevLett.125.196604 33216579
Wang, Y. et al. One-dimensional quasiperiodic Mosaic lattice with exact mobility edges. Phys. Rev. Lett. 125, 196604. 10.1103/PhysRevLett.125.196604 (2020).33216579 10.1103/PhysRevLett.125.196604
29. Biddle J Sarma SD Predicted mobility edges in one-dimensional incommensurate optical lattices: An exactly solvable model of Anderson localization Phys. Rev. Lett. 2010 104 070601 10.1103/PhysRevLett.104.070601 20366867
Biddle, J. & Sarma, S. D. Predicted mobility edges in one-dimensional incommensurate optical lattices: An exactly solvable model of Anderson localization. Phys. Rev. Lett. 104, 070601. 10.1103/PhysRevLett.104.070601 (2010).20366867 10.1103/PhysRevLett.104.070601
30. An FA Padavic K Meier EJ Hegde S Ganeshan S Pixley JH Vishveshwara S Gadway B Interactions and mobility edges: Observing the generalized Aubry-André model Phys. Rev. Lett. 2021 126 040603 10.1103/PhysRevLett.126.040603 33576679
An, F. A. et al. Interactions and mobility edges: Observing the generalized Aubry-André model. Phys. Rev. Lett. 126, 040603. 10.1103/PhysRevLett.126.040603 (2021).33576679 10.1103/PhysRevLett.126.040603
31. An FA Meier EJ Gadway B Engineering a flux-dependent mobility edge in disordered Zigzag chains Phys. Rev. X 2018 8 031045 10.1103/PhysRevX.8.031045
An, F. A., Meier, E. J. & Gadway, B. Engineering a flux-dependent mobility edge in disordered Zigzag chains. Phys. Rev. X 8, 031045. 10.1103/PhysRevX.8.031045 (2018).10.1103/PhysRevX.8.031045
32. Hasan MZ Kane CL Colloquium: Topological insulators Rev. Mod. Phys. 2010 82 3045 10.1103/RevModPhys.82.3045
Hasan, M. Z. & Kane, C. L. Colloquium: Topological insulators. Rev. Mod. Phys. 82, 3045. 10.1103/RevModPhys.82.3045 (2010).10.1103/RevModPhys.82.3045
33. Sarma SD Freedman M Nayak C Topologically protected Qubits from a possible Non-Abelian fractional quantum Hall state Phys. Rev. Lett. 2005 94 166802 10.1103/PhysRevLett.94.166802 15904258
Sarma, S. D., Freedman, M. & Nayak, C. Topologically protected Qubits from a possible Non-Abelian fractional quantum Hall state. Phys. Rev. Lett. 94, 166802. 10.1103/PhysRevLett.94.166802 (2005).15904258 10.1103/PhysRevLett.94.166802
34. Nayak C Simon SH Stern A Freedman M Sarma SD Non-Abelian anyons and topological quantum computation Rev. Mod. Phys. 2008 80 1083 10.1103/RevModPhys.80.1083
Nayak, C., Simon, S. H., Stern, A., Freedman, M. & Sarma, S. D. Non-Abelian anyons and topological quantum computation. Rev. Mod. Phys. 80, 1083. 10.1103/RevModPhys.80.1083 (2008).10.1103/RevModPhys.80.1083
35. Kitaev AY Unpaired Majorana fermions in quantum wires Phys. Usp. 2001 44 131 10.1070/1063-7869/44/10S/S29
Kitaev, A. Y. Unpaired Majorana fermions in quantum wires. Phys. Usp. 44, 131. 10.1070/1063-7869/44/10S/S29 (2001).10.1070/1063-7869/44/10S/S29
36. Lutchyn RM Sau JD Sarma SD Majorana Fermions and a topological phase transition in semiconductor-superconductor heterostructures Phys. Rev. Lett. 2010 105 077001 10.1103/PhysRevLett.105.077001 20868069
Lutchyn, R. M., Sau, J. D. & Sarma, S. D. Majorana Fermions and a topological phase transition in semiconductor-superconductor heterostructures. Phys. Rev. Lett. 105, 077001. 10.1103/PhysRevLett.105.077001 (2010).20868069 10.1103/PhysRevLett.105.077001
37. Oreg Y Refael G Oppen FV Helical liquids and Majorana bound states in quantum wires Phys. Rev. Lett. 2010 105 177002 10.1103/PhysRevLett.105.177002 21231073
Oreg, Y., Refael, G. & Oppen, F. V. Helical liquids and Majorana bound states in quantum wires. Phys. Rev. Lett. 105, 177002. 10.1103/PhysRevLett.105.177002 (2010).21231073 10.1103/PhysRevLett.105.177002
38. Alicea J Majorana fermions in a tunable semiconductor device Phys. Rev. B 2010 81 125318 10.1103/PhysRevB.81.125318
Alicea, J. Majorana fermions in a tunable semiconductor device. Phys. Rev. B 81, 125318. 10.1103/PhysRevB.81.125318 (2010).10.1103/PhysRevB.81.125318
39. Rokhinson LP Liu X Furdyna JK The fractional a.c. Josephson effect in a semiconductor-superconductor nanowire as a signature of Majorana particles Nat. Phys. 2012 8 795 10.1038/nphys2429
Rokhinson, L. P., Liu, X. & Furdyna, J. K. The fractional a.c. Josephson effect in a semiconductor-superconductor nanowire as a signature of Majorana particles. Nat. Phys. 8, 795. 10.1038/nphys2429 (2012).10.1038/nphys2429
40. Das A Ronen Y Most Y Oreg Y Heiblum M Shtrikman H Zero-bias peaks and splitting in an Al-InAs nanowire topological superconductor as a signature of Majorana fermions Nat. Phys. 2012 8 887 10.1038/nphys2479
Das, A. et al. Zero-bias peaks and splitting in an Al-InAs nanowire topological superconductor as a signature of Majorana fermions. Nat. Phys. 8, 887. 10.1038/nphys2479 (2012).10.1038/nphys2479
41. Perge SN Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor Science 2014 346 602 10.1126/science.1259327 25278507
Perge, S. N. et al. Observation of Majorana fermions in ferromagnetic atomic chains on a superconductor. Science 346, 602. 10.1126/science.1259327 (2014).25278507 10.1126/science.1259327
42. Lang L-J Chen S Majorana fermions in density-modulated p-wave superconducting wires Phys. Rev. B 2012 86 205135 10.1103/PhysRevB.86.205135
Lang, L.-J. & Chen, S. Majorana fermions in density-modulated -wave superconducting wires. Phys. Rev. B 86, 205135. 10.1103/PhysRevB.86.205135 (2012).10.1103/PhysRevB.86.205135
43. Cai X Lang L-J Chen S Wang Y Topological superconductor to Anderson localization transition in one-dimensional incommensurate lattices Phys. Rev. Lett. 2013 110 176403 10.1103/PhysRevLett.110.176403 23679750
Cai, X., Lang, L.-J., Chen, S. & Wang, Y. Topological superconductor to Anderson localization transition in one-dimensional incommensurate lattices. Phys. Rev. Lett. 110, 176403. 10.1103/PhysRevLett.110.176403 (2013).23679750 10.1103/PhysRevLett.110.176403
44. Cai X Quantum phase transitions and phase diagram for a one-dimensional p-wave superconductor with an incommensurate potential J. Phys.: Condens. Matter 2014 26 155701 10.1088/0953-8984/26/15/155701 24675766
Cai, X. Quantum phase transitions and phase diagram for a one-dimensional p-wave superconductor with an incommensurate potential. J. Phys.: Condens. Matter 26, 155701. 10.1088/0953-8984/26/15/155701 (2014).24675766 10.1088/0953-8984/26/15/155701
45. Zeng Q-B Chen S Lü R Generalized Aubry-André-Harper model with p-wave superconducting pairing Phys. Rev. B 2016 94 125408 10.1103/PhysRevB.94.125408
Zeng, Q.-B., Chen, S. & Lü, R. Generalized Aubry-André-Harper model with -wave superconducting pairing. Phys. Rev. B 94, 125408. 10.1103/PhysRevB.94.125408 (2016).10.1103/PhysRevB.94.125408
46. Katsura H Schuricht D Takahashi M Exact ground states and topological order in interacting Kitaev/Majorana chains Phys. Rev. B 2015 92 115137 10.1103/PhysRevB.92.115137
Katsura, H., Schuricht, D. & Takahashi, M. Exact ground states and topological order in interacting Kitaev/Majorana chains. Phys. Rev. B 92, 115137. 10.1103/PhysRevB.92.115137 (2015).10.1103/PhysRevB.92.115137
47. Wang Y Miao J-J Jin H-K Chen S Characterization of topological phases of dimerized Kitaev chain via edge correlation functions Phys. Rev. B 2017 96 205428 10.1103/PhysRevB.96.205428
Wang, Y., Miao, J.-J., Jin, H.-K. & Chen, S. Characterization of topological phases of dimerized Kitaev chain via edge correlation functions. Phys. Rev. B 96, 205428. 10.1103/PhysRevB.96.205428 (2017).10.1103/PhysRevB.96.205428
48. Wakatsuki R Ezawa M Tanaka Y Nagaosa N Fermion fractionalization to Majorana fermions in a dimerized Kitaev superconductor Phys. Rev. B 2014 90 014505 10.1103/PhysRevB.90.014505
Wakatsuki, R., Ezawa, M., Tanaka, Y. & Nagaosa, N. Fermion fractionalization to Majorana fermions in a dimerized Kitaev superconductor. Phys. Rev. B 90, 014505. 10.1103/PhysRevB.90.014505 (2014).10.1103/PhysRevB.90.014505
49. Hua C-B Chen R Xu D-H Zhou B Disorder-induced Majorana zero modes in a dimerized Kitaev superconductor chain Phys. Rev. B 2019 100 205302 10.1103/PhysRevB.100.205302
Hua, C.-B., Chen, R., Xu, D.-H. & Zhou, B. Disorder-induced Majorana zero modes in a dimerized Kitaev superconductor chain. Phys. Rev. B 100, 205302. 10.1103/PhysRevB.100.205302 (2019).10.1103/PhysRevB.100.205302
50. Su WP Schrieffer JR Heeger AJ Solitons in Polyacetylene Phys. Rev. Lett. 1979 42 1698 10.1103/PhysRevLett.42.1698
Su, W. P., Schrieffer, J. R. & Heeger, A. J. Solitons in Polyacetylene. Phys. Rev. Lett. 42, 1698. 10.1103/PhysRevLett.42.1698 (1979).10.1103/PhysRevLett.42.1698
51. Roy S Nabi SN Basu S Critical and topological phases of dimerized Kitaev chain in presence of quasiperiodic potential Phys. Rev. B 2023 107 014202 10.1103/PhysRevB.107.014202
Roy, S., Nabi, S. N. & Basu, S. Critical and topological phases of dimerized Kitaev chain in presence of quasiperiodic potential. Phys. Rev. B 107, 014202. 10.1103/PhysRevB.107.014202 (2023).10.1103/PhysRevB.107.014202
52. Perk JHH Capel HW Zuilhof MJ Siskens TJ On a soluble model of an antiferromagnetic chain with alternating interactions and magnetic moments Phys. A 1975 81 319 10.1016/0378-4371(75)90052-7
Perk, J. H. H., Capel, H. W., Zuilhof, M. J. & Siskens, T. J. On a soluble model of an antiferromagnetic chain with alternating interactions and magnetic moments. Phys. A 81, 319. 10.1016/0378-4371(75)90052-7 (1975).10.1016/0378-4371(75)90052-7
53. de Lima JP Anisotropic X-Y model on the inhomogeneous periodic chain Phys. Rev. B 2007 75 214406 10.1103/PhysRevB.75.214406
de Lima, J. P. et al. Anisotropic model on the inhomogeneous periodic chain. Phys. Rev. B 75, 214406. 10.1103/PhysRevB.75.214406 (2007).10.1103/PhysRevB.75.214406
54. Chitov GY Pandey T Timonin PN String and conventional order parameters in the solvable modulated quantum chain Phys. Rev. B 2019 100 104428 10.1103/PhysRevB.100.104428
Chitov, G. Y., Pandey, T. & Timonin, P. N. String and conventional order parameters in the solvable modulated quantum chain. Phys. Rev. B 100, 104428. 10.1103/PhysRevB.100.104428 (2019).10.1103/PhysRevB.100.104428
55. Grifoni M Hänggi P Driven quantum tunneling Phys. Rep. 1998 304 229 10.1016/S0370-1573(98)00022-2
Grifoni, M. & Hänggi, P. Driven quantum tunneling. Phys. Rep. 304, 229. 10.1016/S0370-1573(98)00022-2 (1998).10.1016/S0370-1573(98)00022-2
56. Goldman N Dalibard J Periodically driven quantum systems: Effective Hamiltonians and engineered gauge fields Phys. Rev. X 2014 4 031027 10.1103/PhysRevX.4.031027
Goldman, N. & Dalibard, J. Periodically driven quantum systems: Effective Hamiltonians and engineered gauge fields. Phys. Rev. X 4, 031027. 10.1103/PhysRevX.4.031027 (2014).10.1103/PhysRevX.4.031027
57. Restrepo S Cerrillo J Bastidas VM Angelakis DG Brandes T Driven open quantum systems and Floquet stroboscopic dynamics Phys. Rev. Lett. 2016 117 250401 10.1103/PhysRevLett.117.250401 28036226
Restrepo, S., Cerrillo, J., Bastidas, V. M., Angelakis, D. G. & Brandes, T. Driven open quantum systems and Floquet stroboscopic dynamics. Phys. Rev. Lett. 117, 250401. 10.1103/PhysRevLett.117.250401 (2016).28036226 10.1103/PhysRevLett.117.250401
58. Cayssol J Dóra B Simon F Moessner R Floquet topological insulators Phys. Status Solidi 2013 7 101 10.1002/pssr.201206451
Cayssol, J., Dóra, B., Simon, F. & Moessner, R. Floquet topological insulators. Phys. Status Solidi 7, 101. 10.1002/pssr.201206451 (2013).10.1002/pssr.201206451
59. Rudner MS Lindner NH Berg E Levin M Anomalous edge states and the bulk-edge correspondence for periodically driven two-dimensional systems Phys. Rev. X 2013 3 031005 10.1103/PhysRevX.3.031005
Rudner, M. S., Lindner, N. H., Berg, E. & Levin, M. Anomalous edge states and the bulk-edge correspondence for periodically driven two-dimensional systems. Phys. Rev. X 3, 031005. 10.1103/PhysRevX.3.031005 (2013).10.1103/PhysRevX.3.031005
60. Gómez-León A Platero G Floquet-Bloch theory and topology in periodically driven lattices Phys. Rev. Lett. 2013 110 200403 10.1103/PhysRevLett.110.200403 25167386
Gómez-León, A. & Platero, G. Floquet-Bloch theory and topology in periodically driven lattices. Phys. Rev. Lett. 110, 200403. 10.1103/PhysRevLett.110.200403 (2013).25167386 10.1103/PhysRevLett.110.200403
61. Rudner MS Lindner NH Band structure engineering and non-equilibrium dynamics in Floquet topological insulators Nat. Rev. Phys. 2020 2 229 10.1038/s42254-020-0170-z
Rudner, M. S. & Lindner, N. H. Band structure engineering and non-equilibrium dynamics in Floquet topological insulators. Nat. Rev. Phys. 2, 229. 10.1038/s42254-020-0170-z (2020).10.1038/s42254-020-0170-z
62. Rechtsman M Photonic Floquet topological insulators Nature 2013 496 196 10.1038/nature12066 23579677
Rechtsman, M. et al. Photonic Floquet topological insulators. Nature 496, 196. 10.1038/nature12066 (2013).23579677 10.1038/nature12066
63. Wintersperger K Realization of an anomalous Floquet topological system with ultracold atoms Nat. Phys. 2020 16 1058 10.1038/s41567-020-0949-y
Wintersperger, K. et al. Realization of an anomalous Floquet topological system with ultracold atoms. Nat. Phys. 16, 1058. 10.1038/s41567-020-0949-y (2020).10.1038/s41567-020-0949-y
64. Maczewsky LJ Zeuner JM Nolte S Szameit A Observation of photonic anomalous Floquet topological insulators Nat. Commun. 2017 8 13756 10.1038/ncomms13756 28051080
Maczewsky, L. J., Zeuner, J. M., Nolte, S. & Szameit, A. Observation of photonic anomalous Floquet topological insulators. Nat. Commun. 8, 13756. 10.1038/ncomms13756 (2017).28051080 10.1038/ncomms13756
65. Mukherjee S Experimental observation of anomalous topological edge modes in a slowly driven photonic lattice Nat. Commun. 2017 8 13918 10.1038/ncomms13918 28051060
Mukherjee, S. et al. Experimental observation of anomalous topological edge modes in a slowly driven photonic lattice. Nat. Commun. 8, 13918. 10.1038/ncomms13918 (2017).28051060 10.1038/ncomms13918
66. Jotzu G Experimental realization of the topological Haldane model with ultracold fermions Nature 2014 515 237 10.1038/nature13915 25391960
Jotzu, G. et al. Experimental realization of the topological Haldane model with ultracold fermions. Nature 515, 237. 10.1038/nature13915 (2014).25391960 10.1038/nature13915
67. Aidelsburger M Measuring the Chern number of Hofstadter bands with ultracold bosonic atoms Nat. Phys. 2015 11 162 10.1038/nphys3171
Aidelsburger, M. et al. Measuring the Chern number of Hofstadter bands with ultracold bosonic atoms. Nat. Phys. 11, 162. 10.1038/nphys3171 (2015).10.1038/nphys3171
68. Eckardt A Colloquium: Atomic quantum gases in periodically driven optical lattices Rev. Mod. Phys. 2017 89 011004 10.1103/RevModPhys.89.011004
Eckardt, A. Colloquium: Atomic quantum gases in periodically driven optical lattices. Rev. Mod. Phys. 89, 011004. 10.1103/RevModPhys.89.011004 (2017).10.1103/RevModPhys.89.011004
69. Afzal S Zimmerling TJ Ren Y Perron D Van V Realization of anomalous Floquet insulators in strongly coupled nanophotonic lattices Phys Rev. Lett. 2020 124 253601 10.1103/PhysRevLett.124.253601
Afzal, S., Zimmerling, T. J., Ren, Y., Perron, D. & Van, V. Realization of anomalous Floquet insulators in strongly coupled nanophotonic lattices Phys. Rev. Lett. 124, 253601. 10.1103/PhysRevLett.124.253601 (2020).10.1103/PhysRevLett.124.253601
70. Wang Y Steinberg H Jarillo-Herrero P Gedik N Observation of Floquet-Bloch states on the surface of a topological insulator Science 2013 342 453 10.1126/science.1239834 24159040
Wang, Y., Steinberg, H., Jarillo-Herrero, P. & Gedik, N. Observation of Floquet-Bloch states on the surface of a topological insulator. Science 342, 453. 10.1126/science.1239834 (2013).24159040 10.1126/science.1239834
71. Agrawal A Bandyopadhyay JN Floquet topological phases with high Chern numbers in a periodically driven extended Su-Schrieffer-Heeger model J. Phys. Condens. Matter 2022 34 305401 10.1088/1361-648X/ac6eac
Agrawal, A. & Bandyopadhyay, J. N. Floquet topological phases with high Chern numbers in a periodically driven extended Su-Schrieffer-Heeger model. J. Phys. Condens. Matter 34, 305401. 10.1088/1361-648X/ac6eac (2022).10.1088/1361-648X/ac6eac
72. Li C-F Luan L-N Wang L-C Topological properties of an extend Su-Schrieffer-Heeger model under periodic Kickings Int. J. Theor. Phys. 2020 59 2852 10.1007/s10773-020-04545-7
Li, C.-F., Luan, L.-N. & Wang, L.-C. Topological properties of an extend Su-Schrieffer-Heeger model under periodic Kickings. Int. J. Theor. Phys. 59, 2852. 10.1007/s10773-020-04545-7 (2020).10.1007/s10773-020-04545-7
73. Roy K Basu S Topological properties of a periodically driven Creutz ladder Phys. Rev. B 2023 108 045415 10.1103/PhysRevB.108.045415
Roy, K. & Basu, S. Topological properties of a periodically driven Creutz ladder. Phys. Rev. B 108, 045415. 10.1103/PhysRevB.108.045415 (2023).10.1103/PhysRevB.108.045415
74. Jangjan M Foa Torres LEF Hosseini M-V Floquet topological phase transitions in a periodically quenched dimer Phys. Rev. B 2022 106 224306 10.1103/PhysRevB.106.224306
Jangjan, M., Foa Torres, L. E. F. & Hosseini, M.-V. Floquet topological phase transitions in a periodically quenched dimer. Phys. Rev. B 106, 224306. 10.1103/PhysRevB.106.224306 (2022).10.1103/PhysRevB.106.224306
75. Yang K Xu S Zhou L Zhao Z Xie T Ding Z Ma W Gong J Shi F Du J Observation of Floquet topological phases with large Chern numbers Phys. Rev. B 2022 106 184106 10.1103/PhysRevB.106.184106
Yang, K. et al. Observation of Floquet topological phases with large Chern numbers. Phys. Rev. B 106, 184106. 10.1103/PhysRevB.106.184106 (2022).10.1103/PhysRevB.106.184106
76. Xiong T-S Gong J An J-H Towards large-Chern-number topological phases by periodic quenching Phys. Rev. B 2016 93 184306 10.1103/PhysRevB.93.184306
Xiong, T.-S., Gong, J. & An, J.-H. Towards large-Chern-number topological phases by periodic quenching. Phys. Rev. B 93, 184306. 10.1103/PhysRevB.93.184306 (2016).10.1103/PhysRevB.93.184306
77. Molignini P Edge mode manipulation through commensurate multifrequency driving Phys. Rev. B 2020 102 235143 10.1103/PhysRevB.102.235143
Molignini, P. Edge mode manipulation through commensurate multifrequency driving. Phys. Rev. B 102, 235143. 10.1103/PhysRevB.102.235143 (2020).10.1103/PhysRevB.102.235143
78. Molignini P Chen W Chitra R Universal quantum criticality in static and Floquet-Majorana chains Phys. Rev. B 2018 98 125129 10.1103/PhysRevB.98.125129
Molignini, P., Chen, W. & Chitra, R. Universal quantum criticality in static and Floquet-Majorana chains. Phys. Rev. B 98, 125129. 10.1103/PhysRevB.98.125129 (2018).10.1103/PhysRevB.98.125129
79. Roy, K. & Basu, S. Different driving protocols and multiple Majorana modes in a Rashba coupled superconducting nanowire (2023). arXiv:2405.11307v1. 10.48550/arXiv.2405.11307.
80. Pan Y Wang B Time-crystalline phases and period-doubling oscillations in one-dimensional Floquet topological insulators Phys. Rev. Res. 2020 2 043239 10.1103/PhysRevResearch.2.043239
Pan, Y. & Wang, B. Time-crystalline phases and period-doubling oscillations in one-dimensional Floquet topological insulators. Phys. Rev. Res. 2, 043239. 10.1103/PhysRevResearch.2.043239 (2020).10.1103/PhysRevResearch.2.043239
81. Wang B Quan J Han J Shen X Wu H Pan Y Observation of photonic topological Floquet Time crystals Laser Photon. Rev. 2022 16 2100469 10.1002/lpor.202100469
Wang, B. et al. Observation of photonic topological Floquet Time crystals. Laser Photon. Rev. 16, 2100469. 10.1002/lpor.202100469 (2022).10.1002/lpor.202100469
82. Mondal S Sen D Dutta A Disconnected entanglement entropy as a marker of edge modes in a periodically driven Kitaev chain J. Phys. Condens. Matter 2022 35 085601 10.1088/1361-648X/aca7f7
Mondal, S., Sen, D. & Dutta, A. Disconnected entanglement entropy as a marker of edge modes in a periodically driven Kitaev chain. J. Phys. Condens. Matter 35, 085601. 10.1088/1361-648X/aca7f7 (2022).10.1088/1361-648X/aca7f7
83. Seshadri R Dutta A Sen D Generating a second-order topological insulator with multiple corner states by periodic driving Phys. Rev. B 2019 100 115403 10.1103/PhysRevB.100.115403
Seshadri, R., Dutta, A. & Sen, D. Generating a second-order topological insulator with multiple corner states by periodic driving. Phys. Rev. B 100, 115403. 10.1103/PhysRevB.100.115403 (2019).10.1103/PhysRevB.100.115403
84. Zhou L Gong J Floquet topological phases in a spin-1/2 double kicked rotor Phys. Rev. A 2018 97 063603 10.1103/PhysRevA.97.063603
Zhou, L. & Gong, J. Floquet topological phases in a spin-1/2 double kicked rotor. Phys. Rev. A 97, 063603. 10.1103/PhysRevA.97.063603 (2018).10.1103/PhysRevA.97.063603
85. Ho DYH Gong J Topological effects in chiral symmetric driven systems Phys. Rev. B 2014 90 195419 10.1103/PhysRevB.90.195419
Ho, D. Y. H. & Gong, J. Topological effects in chiral symmetric driven systems. Phys. Rev. B 90, 195419. 10.1103/PhysRevB.90.195419 (2014).10.1103/PhysRevB.90.195419
86. Yates DJ Mitra A Entanglement properties of the time-periodic Kitaev chain Phys. Rev. B 2017 96 115108 10.1103/PhysRevB.96.115108
Yates, D. J. & Mitra, A. Entanglement properties of the time-periodic Kitaev chain. Phys. Rev. B 96, 115108. 10.1103/PhysRevB.96.115108 (2017).10.1103/PhysRevB.96.115108
87. Mondal S Disconnected entanglement entropy as a marker of edge modes in a periodically driven Kitaev chain J. Phys.: Condens. Matter 2023 35 085601 10.1088/1361-648X/aca7f7
Mondal, S. et al. Disconnected entanglement entropy as a marker of edge modes in a periodically driven Kitaev chain. J. Phys.: Condens. Matter 35, 085601. 10.1088/1361-648X/aca7f7 (2023).10.1088/1361-648X/aca7f7
88. Qin P Yin C Chen S Dynamical Anderson transition in one-dimensional periodically kicked incommensurate lattices Phys. Rev. B 2014 90 054303 10.1103/PhysRevB.90.054303
Qin, P., Yin, C. & Chen, S. Dynamical Anderson transition in one-dimensional periodically kicked incommensurate lattices. Phys. Rev. B 90, 054303. 10.1103/PhysRevB.90.054303 (2014).10.1103/PhysRevB.90.054303
89. Zhang Y Zhou B Hu H Chen S Localization, multifractality, and many-body localization in periodically kicked quasiperiodic lattices Phys. Rev. B 2022 106 054312 10.1103/PhysRevB.106.054312
Zhang, Y., Zhou, B., Hu, H. & Chen, S. Localization, multifractality, and many-body localization in periodically kicked quasiperiodic lattices. Phys. Rev. B 106, 054312. 10.1103/PhysRevB.106.054312 (2022).10.1103/PhysRevB.106.054312
90. Sarkar M Ghosh R Sen A Sengupta K Mobility edge and multifractality in a periodically driven Aubry-André model Phys. Rev. B 2021 103 184309 10.1103/PhysRevB.103.184309
Sarkar, M., Ghosh, R., Sen, A. & Sengupta, K. Mobility edge and multifractality in a periodically driven Aubry-André model. Phys. Rev. B 103, 184309. 10.1103/PhysRevB.103.184309 (2021).10.1103/PhysRevB.103.184309
91. Aditya S Sengupta K Sen D Periodically driven model with quasiperiodic potential and staggered hopping amplitudes: Engineering of mobility gaps and multifractal states Phys. Rev. B 2023 107 035402 10.1103/PhysRevB.107.035402
Aditya, S., Sengupta, K. & Sen, D. Periodically driven model with quasiperiodic potential and staggered hopping amplitudes: Engineering of mobility gaps and multifractal states. Phys. Rev. B 107, 035402. 10.1103/PhysRevB.107.035402 (2023).10.1103/PhysRevB.107.035402
92. Evers F Mirlin AD Anderson transitions Rev. Mod. Phys. 2008 80 1355 10.1103/RevModPhys.80.1355
Evers, F. & Mirlin, A. D. Anderson transitions. Rev. Mod. Phys. 80, 1355. 10.1103/RevModPhys.80.1355 (2008).10.1103/RevModPhys.80.1355
93. Roy, S., Basu, S. & Khyamovich, I.-M. Ergodicity-breaking phase diagram and fractal dimensions in long-range models with generically correlated disorder (2023). arXiv:2307.03085v1. 10.48550/arXiv.2307.03085.
94. Kravtsov VE A random matrix model with localization and ergodic transitions New J. Phys. 2015 17 122002 10.1088/1367-2630/17/12/122002
Kravtsov, V. E. et al. A random matrix model with localization and ergodic transitions. New J. Phys. 17, 122002. 10.1088/1367-2630/17/12/122002 (2015).10.1088/1367-2630/17/12/122002
95 ČadeŽ T Mondaini R Sacramento PD Edge and bulk localization of Floquet topological superconductors Phys. Rev. B 2019 99 014301 10.1103/PhysRevB.99.014301
ČadeŽ, T., Mondaini, R. & Sacramento, P. D. Edge and bulk localization of Floquet topological superconductors. Phys. Rev. B 99, 014301. 10.1103/PhysRevB.99.014301 (2019).10.1103/PhysRevB.99.014301
96. Mondragon-Shem I Hughes TL Song J Prodan E Topological criticality in the chiral-symmetric AIII class at strong disorder Phys. Rev. Lett. 2014 113 046802 10.1103/PhysRevLett.113.046802 25105642
Mondragon-Shem, I., Hughes, T. L., Song, J. & Prodan, E. Topological criticality in the chiral-symmetric AIII class at strong disorder. Phys. Rev. Lett. 113, 046802. 10.1103/PhysRevLett.113.046802 (2014).25105642 10.1103/PhysRevLett.113.046802
97. Roy R Harper F Periodic table for Floquet topological insulators Phys. Rev. B 2017 96 155118 10.1103/PhysRevB.96.155118
Roy, R. & Harper, F. Periodic table for Floquet topological insulators. Phys. Rev. B 96, 155118. 10.1103/PhysRevB.96.155118 (2017).10.1103/PhysRevB.96.155118
98. Asbóth JK Tarasinski B Delplace P Chiral symmetry and bulk-boundary correspondence in periodically driven one-dimensional systems Phys. Rev. B 2014 90 125143 10.1103/PhysRevB.90.125143
Asbóth, J. K., Tarasinski, B. & Delplace, P. Chiral symmetry and bulk-boundary correspondence in periodically driven one-dimensional systems. Phys. Rev. B 90, 125143. 10.1103/PhysRevB.90.125143 (2014).10.1103/PhysRevB.90.125143
99. Asbóth JK Obuse H Bulk-boundary correspondence for chiral symmetric quantum walks Phys. Rev. B 2013 88 121406 10.1103/PhysRevB.88.121406
Asbóth, J. K. & Obuse, H. Bulk-boundary correspondence for chiral symmetric quantum walks. Phys. Rev. B 88, 121406. 10.1103/PhysRevB.88.121406 (2013).10.1103/PhysRevB.88.121406
100. Wang Y Zhang L Sun W Ting-Fung PJ Liu X-J Quantum phase with coexisting localized, extended, and critical zones Phys. Rev. B 2022 106 L140203 10.1103/PhysRevB.106.L140203
Wang, Y., Zhang, L., Sun, W., Ting-Fung, P. J. & Liu, X.-J. Quantum phase with coexisting localized, extended, and critical zones. Phys. Rev. B 106, L140203. 10.1103/PhysRevB.106.L140203 (2022).10.1103/PhysRevB.106.L140203
101. Oganesyan V Huse DA Localization of interacting fermions at high temperature Phys. Rev. B 2007 75 155111 10.1103/PhysRevB.75.155111
Oganesyan, V. & Huse, D. A. Localization of interacting fermions at high temperature. Phys. Rev. B 75, 155111. 10.1103/PhysRevB.75.155111 (2007).10.1103/PhysRevB.75.155111
102. D’Alessio L Rigol M Long-time behavior of isolated periodically driven interacting lattice systems Phys. Rev. X 2014 4 041048 10.1103/PhysRevX.4.041048
D’Alessio, L. & Rigol, M. Long-time behavior of isolated periodically driven interacting lattice systems. Phys. Rev. X 4, 041048. 10.1103/PhysRevX.4.041048 (2014).10.1103/PhysRevX.4.041048
103. Deng X Ray S Sinha S Shlyapnikov GV Santos L One-dimensional quasicrystals with power-law Hopping Phys. Rev. Lett. 2019 123 025301 10.1103/PhysRevLett.123.025301 31386526
Deng, X., Ray, S., Sinha, S., Shlyapnikov, G. V. & Santos, L. One-dimensional quasicrystals with power-law Hopping. Phys. Rev. Lett. 123, 025301. 10.1103/PhysRevLett.123.025301 (2019).31386526 10.1103/PhysRevLett.123.025301
104. Yao H Khoudli H Bresque L Sanchez-Palencia L Critical behavior and fractality in shallow one-dimensional quasiperiodic potentials Phys. Rev. Lett. 2019 123 070405 10.1103/PhysRevLett.123.070405 31491103
Yao, H., Khoudli, H., Bresque, L. & Sanchez-Palencia, L. Critical behavior and fractality in shallow one-dimensional quasiperiodic potentials. Phys. Rev. Lett. 123, 070405. 10.1103/PhysRevLett.123.070405 (2019).31491103 10.1103/PhysRevLett.123.070405
105. Tang C Kohmoto M Global scaling properties of the spectrum for a quasiperiodic schrödinger equation Phys. Rev. B 1986 34 2041(R) 10.1103/PhysRevB.34.2041
Tang, C. & Kohmoto, M. Global scaling properties of the spectrum for a quasiperiodic schrödinger equation. Phys. Rev. B 34, 2041(R). 10.1103/PhysRevB.34.2041 (1986).10.1103/PhysRevB.34.2041
106. Helbig T Hofmann T Lee CH Thomale R Imhof S Molenkamp LW Kiessling T Band structure engineering and reconstruction in electric circuit networks Phys. Rev. B 2019 99 161114 10.1103/PhysRevB.99.161114
Helbig, T. et al. Band structure engineering and reconstruction in electric circuit networks. Phys. Rev. B 99, 161114. 10.1103/PhysRevB.99.161114 (2019).10.1103/PhysRevB.99.161114
107. Ezawa M Braiding of Majorana-like corner states in electric circuits and its non-Hermitian generalization Phys. Rev. B 2019 100 045407 10.1103/PhysRevB.100.045407
Ezawa, M. Braiding of Majorana-like corner states in electric circuits and its non-Hermitian generalization. Phys. Rev. B 100, 045407. 10.1103/PhysRevB.100.045407 (2019).10.1103/PhysRevB.100.045407
108. Halder D Thomale R Basu S Circuit realization of a two-orbital non-Hermitian tight-binding chain Phys. Rev. B 2024 109 115407 10.1103/PhysRevB.109.115407
Halder, D., Thomale, R. & Basu, S. Circuit realization of a two-orbital non-Hermitian tight-binding chain. Phys. Rev. B 109, 115407. 10.1103/PhysRevB.109.115407 (2024).10.1103/PhysRevB.109.115407
109. Yao J Majorana-like end states in one-dimensional dimerized Kitaev topolectrical circuit New J. Phys. 2022 24 043032 10.1088/1367-2630/ac6230
Yao, J. et al. Majorana-like end states in one-dimensional dimerized Kitaev topolectrical circuit. New J. Phys. 24, 043032. 10.1088/1367-2630/ac6230 (2022).10.1088/1367-2630/ac6230
110. Dabiri SS Cheraghchi H Electric circuit simulation of Floquet topological insulators in Fourier space J. Appl. Phys. 2023 134 084303 10.1063/5.0150118
Dabiri, S. S. & Cheraghchi, H. Electric circuit simulation of Floquet topological insulators in Fourier space. J. Appl. Phys. 134, 084303. 10.1063/5.0150118 (2023).10.1063/5.0150118
