Identifying the mobility edges in a one-dimensional incommensurate model with p-wave superfluid

The mobility edges which separate the localized energy eigenstates from the extended ones exist normally only in three dimensional systems. For one-dimensional systems with random on-site potentials, one never encounters mobility edges, where all the eigenstates are localized. However, there are two kinds of 1D systems such as correlated disordered models, and the systems of exponentially decaying hopping kinetics, features of mobility edges at some specific values become possible. We study in this paper the properties of the mobility edges in a one-dimensional p-wave superfluid on an incommensurate lattice with exponentially decaying hopping kinetics. Without the p-wave superluid, the system displays a single mobility edge, which separates the extended regime from the localized one at a certain energy. Without the exponentially decaying hopping term, the system displays a phase transition from a topological superconductor to an Anderson localization at a certain disorder strength, where no mobility edge exists. We are interested in the influence of the p-wave superfluid on the mobility edge. By solving the Bogoliubov-de Gennes equation, the eigenvalues and the eigenfunctions are obtained. In order to identify the extending or the localized properties of the eigenvectors, we define an inverse participation ratio IPR. For an extended state, IPRn1/L which goes to zero at a large L, and for a localized one, IPRn being constant. Therefore, the IPR can be taken as a criterion to distinguish the extended state from the localized one, while the mobility edge is defined as the boundary between two different states. We find that, with a p-wave superfluid, the system changes from a single mobility edge to a multiple one, and the number of mobility edges increases with the increased superfluid pairing order parameter. To further obtain the energy or the location of the mobility edge, we investigate the scaling behavior of wave functions by using a multifractal analysis, which is calculated through the scaling index . The minimum value of the index, with the values min= 1, 0min1, and min= 0, mean the extended, critical, and localized states, respectively. For the two consecutive states, the minima of the scaling index min when extrapolating to the large size limit between 0 and 1 signal the mobility edge. By exploring the corresponding Bogoliubov quasi-particle wave functions for the system under open boundary conditions together with the multifractal analysis for the system under periodic boundary conditions, we identify two mobility edges for the system of the p-wave superfluid pairing. Furthermore, we will investigate how the existence of the mobility edges influences the p-wave superfluid, and identify the phase diagram at the given parameters. We will in the future try to understand the relationship between the topological superfluid and the mobility edges.