Floquet Topological Superfluid and Majorana Zero Modes in Two-Dimensional Periodically Driven Fermi Systems

We propose a simple approach to realize two-dimensional Floquet topological superfluid by periodically tuning the depth of square optical lattice potentials. We show that the periodic driving can induce topological phase transitions between trivial superfluid and Floquet topological superfluid. For this systems we verify the anomalous bulk-boundary correspondence, namely that the robust chiral Floquet edge states can appear even when the winding number of all the bulk Floquet bands is zero. We establish the existence of two Floquet Majorana zero modes separated in the quasienergy space, with ε0,π = 0,π/T at the topological defects.


Floquet Topological Superfluid and Majorana Zero Modes in Two-Dimensional Periodically Driven Fermi Systems
Xiaosen Yang 1 , Beibing Huang 2 & Zhengling Wang 1 We propose a simple approach to realize two-dimensional Floquet topological superfluid by periodically tuning the depth of square optical lattice potentials. We show that the periodic driving can induce topological phase transitions between trivial superfluid and Floquet topological superfluid. For this systems we verify the anomalous bulk-boundary correspondence, namely that the robust chiral Floquet edge states can appear even when the winding number of all the bulk Floquet bands is zero. We establish the existence of two Floquet Majorana zero modes separated in the quasienergy space, with ε 0,π = 0,π/T at the topological defects.
In this paper, we shall study two-dimensional Floquet Fermi systems with cold atoms. The two dimensionality is special in that non-Abelian statistics can be realized here. With the motivation of potential applications in non-Abelian statistics and topological quantum computation, we propose a simple approach to realize two-dimensional Floquet topological superfluids by periodically tuning the depth of optical lattice potentials. We demonstrate that Floquet topological superfluid phase arises when the driving frequency is lowered below the bandwidth, when robust chiral edge states span the gap at ε π . When the frequency if further lowered to below half bandwidth, robust chiral edge states span both the gaps around ε 0,π , with the Chern number of all the Floquet bands vanishing. We also show that the Floquet topological nontrivial phases are weak pairing phases and have rather clear Fermi surfaces, thus the topological phase transition is a transition from the strong pairing phase to the weak pairing phase. Additionally, we show that two flavor Floquet Majorana zero modes are localized at the topological defects in Floquet topological superfluid, moreover, the two flavors Floquet Majorana zero modes are separated at quasienergy space.

Results
Floquet Theorem. We consider a spin-polarized Fermi gas loaded on square optical lattice potentials The topological properties can be characterized by a winding number, which is given by the following well-known formulaˆˆ∫ The superfluid is topological nontrivial (W ≠ 0) for 0 < μ < 8J 0 and trivial (W = 0) for μ > 8J 0 and μ < 0. Now let us tune the depth of the lattice potential with period T, as a result, the nearest neighbor tunneling amplitude is periodically varying. For simplicity, we assume that the tunneling amplitude varies like J(t) = J 0 − J D cos(ωt). The pairing amplitude(Δ) changes much smaller than the tunneling amplitude(J) for intermediate and strong attraction. Additionally, the fermions cannot follow the shaking of the tunneling quickly in high frequency limit 45 . Therefore, we can assume the Δ is time independent. Consequentially, the above Hamiltonian becomes periodically time-dependent and satisfies H(t + T) = H(t) with period T = 2π/ω. The periodic time-dependent Hamiltonian can be rewritten aŝ † To study the properties of the driven system, we start from the Schrödinger equation 57 : According to the Floquet theorem, the wave function satisfies Ψ(k, t) = e −iε(k)t Φ(k, t) with Floquet states Φ(k,t) satisfying Φ(k, t) = Φ(k, T + t) and the Floquet equation Here the quasienergies ε(k) are defined modulo 2π/T, which is analogous to the lattice momentum in the Bloch band theory. Therefore, there exist two particle-hole conjugated quasienergies ε 0, π = 0, π/T in the quasienergy spectrum, in contrast to the static systems with only one particle-hole conjugated energy ε 0 . This feature implies nontrivial effects on the Floquet systems.
The Floquet states can be expanded as . The topological properties of the periodically driven superfluid are defined in terms of the time-independent Floquet Hamiltonian in Eq. (7).

Floquet Topological Superfluid.
In accordance with the general principle of bulk-boundary correspondence, the topological properties of a bulk system can be characterized by the robust edge states. To investigate the topological properties of the periodically driven systems, we consider the Floquet Hamiltonian in Eq. (7) in a strip geometry, in which the edges are along the x direction at L = 0, 60.
As a comparison, first we present numerical results without periodic driving, as shown in Fig. 1(a). As discussed above, without the driving the superfluid is topological trivial when the chemical potential μ > 8J 0 or μ < 0. Accordingly, there should be no robust edge states at the open boundaries, which is verified in our numerical calculations. On the other hand, when the system is driven with a frequency much greater than other energy scales, there exist a large gap at ε π , and the system is again topologically trivial. Suppose that we gradually decrease the driving frequency. The gap at ε π can be closed and reopened at Max(2 ) k c , wherein the nearest Floquet bands (up branch of the m-th Floquet band and down branch of the m + 1-th Floquet band) become inverted. When ω < 2W, robust chiral Floquet edge states span the gap near ε π , as shown in Fig. 1(b). The chiral Floquet edge states are different from those appearing in static topological superfluid in that edge states span the gap at ε π instead at ε 0 . A pair chiral Floquet edge states are localized at the two boundaries of the system, which can been found in Fig. 1(c). Therefore, changing driving frequency induces a topological phase transition from a trivial superfluid phase to Floquet topological superfluid phase.
To see the topological phase transition from a more transparent perspective, we can derive a time-independent effective Hamiltonian H eff (k) by evolution operator over a period T: − −  To derive the effective Hamiltonian near ε π for adjacent Floquet bands(such as between up branch of m Floquet band and down branch of m + 1 Floquet band). Following the method of ref 24 . we go to a rotating frame through an time-dependent unitary transformation Taking advantage of this effective Hamiltonian, we can use the winding number given in the previous section to analyze the topological properties of our system. It is found that, for the parameters of Fig. 1, the driving induces a topological phase transition from topological trivial superfluid to Floquet topological superfluid at ω = 2E (π,π) = 20, which is consistent with our previous statement.
Decreasing the driving frequency further, the gap at ε 0 can be closed and reopened again at The robust Floquet edge states will span both the gaps at ε 0 and ε π as shown in Fig. 1(d). These chiral edge states also localized at the boundaries of the system.
An interesting phenomenon is that both the ε 0,π gaps are spanned by a pair chiral Floquet edge states, as shown in Fig. 1(d). Due to the bulk-boundary correspondence, the winding number of a band is equal to the difference between the number of edge states at the gaps above and below the band, in other words, we have C εε′ = n edge (ε) − n edge (ε′). Therefore, the winding number of all the Floquet bands in Fig. 1(d) are zero. This is very different from the static cases, for which the robust chiral edge states only span the gap at ε 0 , and only when the total winding number of bands below zero energy is nonzero. This anomalous bulk-edge correspondence also emerge in Floquet topological insulator 23 . As another prominent difference, the Floquet topological superfluid can possess two types robust Floquet edge states within the gaps at ε 0 and ε π respectively 30 . This will induce two types Floquet Majorana zero modes at topological defects of Floquet topological superfluid, which we will discuss in the next sections. Certainly, we can also create Floquet topological superfluid by adding suitable periodic driving to a topological nontrivial superfluid for 0 < μ < 8J 0 . There also exist two flavors of robust chiral edge states at the boundaries, even though the winding number of all the Foquet bands are zero.
In above analysis the topological phase transition does not depend on the driving strength J 0 , which only relates the values of the inverse band gaps at ε 0,π .

Phase Transition between Strong and Weak Pairing phases.
For p x + ip y paired system 51 , the superfluid is topological trivial for strong pairing phase (μ < 0) and nontrivial for weak pairing phase (μ > 0). For the topological nontrivial phase, the pairing only changes the Fermi surface slightly. So, the Fermi surfaces is clear. Nevertheless, there is no Fermi surface in strong paired topological trivial phase. Now, we show that the above phase transition induced by periodic driving is transition from strong to weak pairing phase.
Firstly, we use the changing of the Fermi surface to show the transition between strong and weak pairing. In the driven superfluid system, we can get the density distribution from the Floquet equation (6). The Floquet state can be written as Ψ(k, t) = (u(k, t), υ(k, t)) T and φ m (k) = (u m (k), υ m (k)) T the density distribution is T m m 0 2 Figure 2 shows the dispersions and density distributions for the two cases of Fig. 1(a) and (b). The driving only changes the dispersion slightly as shown in Fig. 2(a),(b) but completely changes the density distribution as shown in Fig. 2(c),(d). For the topological trivial cases, the pairing is strong and there is no Fermi surface as shown in Fig. 2(c). In the presence of periodic driving, the Floquet topological superfluid has clear Fermi surfaces in the density distribution as shown in Fig. 2(d). Therefore, the periodic driving drives the strong paired topological trivial phase into weak paired Floquet topological phase. Figure 2 also shows the density of TRI momentum (π, π) has a universal jump at ω = 2E (π,π) . The density of other TRI momenta will also have jump. The density will jump from 0 to 1 at ω = 2E (0,π) for [(π, 0), (0, π)] and Lastly, we use the effective Hamiltonian to analyze the topological phase transition at ω = 2E (π,π) . Here, we only consider the properties of system near the phase transition point and let ω = 2E (π,π) + δ. So, the effective Hamiltonian can be expanded at k c = (π, π) and can be written as n A k A k ( ( ) , ( ) , sgn( )/2) 2 . Rewriting Δ′ = A(δ) and μ δ ε ′ = sgn( )/2 k c , the effective Hamiltonian has the same form with ref. 51 . Thus, the superfluid is strong pairing phase for ω > 2E (π,π) , weak pairing phase for ω < 2E (π,π) and the transition is at ω = 2E (π,π) .
Floquet Majorana Zero Modes. Fundamentally different from the static topological superfluid phase, Floquet topological superfluid phase has nontrivial robust Floquet edge states spanning the gap at ε π 30 . Now, we will show that this can be used to generate Floquet Majorana zero modes with finite quasienergy ε π at topological defects. To generate the Floquet Majorana zero modes, we add two 'π-Flux' to the driven systems by changing the sign of the links cut by the line between two separate sites ' A' and 'B' . Fermions will get a π phase by circling ' A' or 'B' sites. Figure 3(a) and (b) show the quasienergies, only around ε 0 and ε π respectively, of the topological defected system with the parameters as Fig. 1(d). For the particle-hole symmetry, there are two degenerate inner gap quasienergies ε 0 in Fig. 3(a) and ε π in Fig. 3(b). The wave functions of the two degenerate quasienergies also conjugate with each other. Figure 3(c) and (d) show the wave function with quasienergies ε 0 and ε π respectively. The wave functions Φ 0 and Φ π are localized at the same 'Flux' sites but separated in the quasienergy space at ε 0 and ε π respectively. These localized wave functions are the Floquet Majorana zero modes. Majorana zero modes are the states that they are their own conjugate. As for the particle-hole symmetry of the superfluid, the Majorana zero modes only can exist at ε 0 in static systems for only the state of ε 0 can be their own particle-hole conjugate. In the presence of periodic driving, the quasienergies of the driven systems are periodic. Therefore, the states of ε π = π/T ≡ −π/T also can be their own particle-hole conjugate. Thus, there exists two types Floquet Majorana zero modes with quasienergy ε 0,π in Floquet topological superfluid. Here, the Floquet Majorana zero modes can be thought as the localization of the chiral Floquet edge states at the topological defects.

Conclusion
In the present paper, we have proposed a simple scheme to realize the two-dimensional Floquet topological nontrivial superfluid by periodically tuning the depth of square optical lattice potentials. The periodic driving can induces a transition from strong pairing phase to weak pairing phase. The weak pairing phases are Floquet

Methods
The evolution operator can be written as To derive the effective Hamiltonian, we can expand the time-periodic Hamiltonian in a Fourier expansion: Then,    . For our driving, sin ,  Data Availability. All data generated or analysed during this study are included in this published article.