Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Topological characterization and stability of Floquet Majorana modes in Rashba nanowires

Debashish Mondal, Arnob Kumar Ghosh, Tanay Nag, and Arijit Saha
Phys. Rev. B 107, 035427 – Published 24 January 2023
PDFHTMLExport Citation
Abstract
Authors
Article Text
  • INTRODUCTION
  • MODEL
  • FLOQUET MAJORANA MODES AND THEIR…
  • STABILITY OF FLOQUET MEMs IN THE…
  • GENERATION OF FLOQUET MEMs USING…
  • POSSIBLE EXPERIMENTAL CONNECTION
  • SUMMARY AND CONCLUSIONS
  • ACKNOWLEDGMENTS
    • References

    Abstract

    We theoretically investigate a practically realizable Floquet topological superconductor model based on a one-dimensional Rashba nanowire and proximity-induced s-wave superconductivity in the presence of a Zeeman field. The driven system hosts regular 0-Majorana end modes and anomalous π-Majorana end modes (MEMs). By tuning the chemical potential and the frequency of the drive, we illustrate the generation of multiple MEMs in our theoretical setup. We utilize the chiral symmetry operator to topologically characterize these MEMs via a dynamical winding number constructed out of the periodized evolution operator. Interestingly, the robustness of the 0- and π-MEMs is established in the presence of on-site time-independent random disorder potential. We employ the twisted boundary condition to define the dynamical topological invariant for this translational-symmetry broken system. The interplay between the Floquet driving and the weak disorder can stabilize the MEMs, giving rise to a quantized value of the dynamical winding number for a finite range of drive parameters. This observation might be experimentally helpful in scrutinizing the topological nature of the Floquet MEMs. We showcase another driving protocol, namely, a periodic kick in the chemical potential, to study the generation of Floquet MEMs in our setup. Our work paves a realistic way to engineer multiple MEMs in a driven system.

    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Figure
    • Received 28 July 2022
    • Revised 6 January 2023
    • Accepted 10 January 2023

    DOI:https://doi.org/10.1103/PhysRevB.107.035427

    ©2023 American Physical Society

    Physics Subject Headings (PhySH)

    1. Research Areas
    Majorana fermionsTopological superconductors
    1. Physical Systems
    Floquet systemsNanowires
    Condensed Matter, Materials & Applied Physics

    Authors & Affiliations

    Debashish Mondal1,2,*, Arnob Kumar Ghosh1,2,†, Tanay Nag3,‡, and Arijit Saha1,2,§

    • 1Institute of Physics, Sachivalaya Marg, Bhubaneswar 751005, India
    • 2Homi Bhabha National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India
    • 3Department of Physics and Astronomy, Uppsala University, Box 516, 75120 Uppsala, Sweden

    • *debashish.m@iopb.res.in
    • †arnob@iopb.res.in
    • ‡tanay.nag@physics.uu.se
    • §arijit@iopb.res.in

    Article Text (Subscription Required)

    Click to Expand

    References (Subscription Required)

    Click to Expand
    Issue

    Vol. 107, Iss. 3 — 15 January 2023

    Reuse & Permissions
    Access Options
    Author publication services for translation and copyediting assistance advertisement

    Authorization Required

    Log In
    Other Options
    ×

    Images

    • Figure 1
      Figure 1

      (a) Schematic representation of our setup is shown here. One-dimensional nanowire (NW) (blue) with strong spin-orbit coupling is placed on top of an s-wave superconducting slab (green). A magnetic field Bx is applied along the x direction and a gate voltage Vg is applied across the cross section of the NW to control the chemical potential. The Majorana end modes (MEMs) [red] appear at the two ends of the NW. (b) Schematic demonstration of our periodic three-step drive protocol [see Eq. (5)] is presented.

      Reuse & Permissions
    • Figure 2
      Figure 2

      (a) Energy eigenvalue spectra for the static Hamiltonian [Eq. (1)], employing the open boundary condition, are shown as a function of the magnetic field Bx. The corresponding winding number W is depicted as a function of Bx in panel (b) for a fixed value of c0=1.0. (c) We demonstrate W in the c0−Bx plane to illustrate the topological phase diagram. Other model parameters are chosen as (t,u,Δ)=(1.0,0.5,1.0). MZMs are obtained for |Bx1|≤Bx≤|Bx2| (see text for details).

      Reuse & Permissions
    • Figure 3
      Figure 3

      Bulk quasienergy spectra for the three-step drive protocol [Eq. (5)] are shown for Ω=1.5 and Ω=2.5 in panels (a) and (b), respectively. In the insets I1 and I2 of panel (a) [inset of panel (b)], we depict the bulk gap around quasienergy E(k)=π and E(k)=0, respectively [quasienergy E(k)=π]. In panels (c) and (d), we illustrate the quasienergy spectra Em as a function of the state index m for our system obeying OBC, corresponding to panels (a) and (b), respectively, while in the insets, we portray the zoomed-in quasienergies for better clarity. The appearance of 0- and π-Majorana end modes (MEMs) are evident from these figures. We show the normalized site-resolved probability (|Ψ|2) corresponding to the 0- and π-MEMs in panels (e) and (f), respectively. Here, we consider N=4000 lattice sites to obtain sharp Majorana localization at individual ends of the NW. All the other model parameters are chosen as (c0,c1,Δ,Bx,t,u)=(1.0,0.2,1.0,1.0,1.0,0.5).

      Reuse & Permissions
    • Figure 4
      Figure 4

      We demonstrate the dynamical winding numbers W0 and Wπ (quasienergy gap corresponding to 0-energy, G0 and π-energy, Gπ) in the c0−Ω plane in panels (a) and (b) [panels (c) and (d)], respectively. Here the color bar represents the winding number (energy gap) for panels (a) and (b) [panels (c) and (d)]. All the other parameters are chosen to be the same as those mentioned in Fig. 3.

      Reuse & Permissions
    • Figure 5
      Figure 5

      In panels (a) and (b), we depict the dynamical winding numbers W0 and Wπ, respectively, in the presence of small (w=0.1) on-site disorder as a function of the driving frequency Ω. We repeat panels (a) and (b) in panels (c) and (d) and panels (e) and (f) for moderate (w=0.5) and strong (w=0.8) disorder strength, respectively. Here, green (dashed), red (solid), and blue (dotted-dashed) lines represent the winding number (Wεc) for the corresponding clean system, the disorder-averaged winding number (Wεd) for the disordered system, and the variance (σ2) of the winding number of the disordered system from the mean value Wεd, respectively, for quasienergy ε. We choose c0=0.2, while all the other model parameter values remain the same as those mentioned in Fig. 3.

      Reuse & Permissions
    • Figure 6
      Figure 6

      Panels (a) and (b) represent the bulk quasienergy spectra for the periodic kick drive protocol [Eq. (17)] for Ω=3.0 and Ω=5.0, respectively. The insets I1 and I2 of panel (a) [inset of panel (b)] exhibit(s) the bulk gap around quasienergies E(k)=π and E(k)=0, respectively [quasienergy E(k)=π]. In panels (c) and (d), we illustrate the quasienergy spectra Em as a function of the state index m while system obeys the OBC, corresponding to panels (a) and (b), respectively. In the insets of these panels, we depict the zoomed-in quasienergies for better clarity. Here, we consider N=4000 lattice sites. All the other model parameters are chosen as (c0,c1,Δ,Bx,t,u)=(1.0,0.2,1.0,1.0,1.0,0.5).

      Reuse & Permissions
    • Figure 7
      Figure 7

      We demonstrate the dynamical winding numbers W0 and Wπ in the c0−Ω plane in panels (a) and (b), respectively, for periodic-kick drive protocol. All the model parameter values remain the same as those mentioned in Fig. 6.

      Reuse & Permissions
    ×

    Sign up to receive regular email alerts from Physical Review B

    Log In

    Cancel
    ×

    Search

    Article Lookup

    Paste a citation or DOI
    Enter a citation
    ×