Mobility edge as one of the most important concepts in a disordered system in which there exists an energy dependent conductor-to-insulator transition has aroused great interest. Unlike an arbitrarily small disorder inducing the Anderson localization in one-dimensional random potential, the well-known Aubry-André model presents a metal-to-insulator transition without mobility edges. Some generalized Aubry-André models are proposed whose the mobility edges in compactly analytic forms are found. However, the existence of the many-body mobility edges in thermodynamic limit for an interacting disordered system is still an open question due to the dimension of the Hilbert space beyond the numerical capacity. In this paper, we demonstrate the existence of the mobility edges of bosonic pairs trapped in one dimensional quasi-periodical lattices subjected to strongly interactions. We believe that our theory will provide a new insight into the studying of the many-body mobility edges.

Two strongly interacting bosons are trapped in an incommensurate model, which is described as \hat H = - J\sum\limits_j \left( \hat c_j^\dagger \hat c_j + 1 + \rmh\rm.c\rm. \right) + 2\lambda \sum\limits_j \dfrac\cos \left( 2\textπ\alpha j \right)1 - b\cos \left( 2\textπ\alpha j \right) \hat n_j + \dfracU2\sum\limits_j \hat n_j\left( \hat n_j - 1 \right) , where there exists no interaction, the system displays mobility edges at b\varepsilon = 2(J - \lambda ), which separates the extended regime from the localized one and b = 0 is the standard Aubry-André model. By applying the perturbation method to the third order in a strong interaction case, we can induce an effective Hamiltonian for bosonic pairs. In the small b case, the bosonic pairs present the mobility edges in a simple closed expression form b\left( \dfracE^2U - E - \dfrac4E \right) = - 4\left(\dfrac1E + \lambda \right), which is the central result of the paper. In order to identify our results numerically, we define a normalized participation ratio (NPR) \eta (E) to discriminate between the extended properties of the many-body eigenvectors and the localized ones. In the thermodynamic limit, the NPR tends to 0 for a localized state, while it is finite for an extended state. The numerical calculations finely coincide with the analytic results for b = 0 and small b cases. Especially, for the b = 0 case, the mobility edges of the bosonic pairs are described as \lambda = - 1/E. The extended regime and the one with the mobility edges will vanish with the interaction U increasing to infinity. We also study the scaling of the NPR with system size in both extended and localized regimes. For the extended state the NPR \eta (E) \propto 1/L tends to a finite value with the increase of L and L \to \infty , while for the localized case, \eta (E) \propto (1/L)^2 tends to zero when L \to \infty . The b \to 1 limit is also considered. As the modulated potential approaches to a singularity when b \to 1, the analytic expression does not fit very well. However, the numerical results indicate that the mobility edges of bosonic pairs still exist. We will try to consider the detection of the mobility edges of the bosonic pairs in the future.