“…Given the complexity of the single particle problem, and the critical behavior produced by the many body interaction on free particles, it is not surprising that quantitative predictions for the interacting Aubry-André model are a difficult task, and have been pursued mostly by numerical methods in the fermionic [14], [15], [16], [17] and bosonic case [18], [19], [20]; non perturbative effects could however be missed due to size limitations. On the analytical side, in [21] it was studied a system of interacting fermions subject to a weak quasi-periodic potential with exponentially fast decaying harmonics, and insulating behavior was rigorously proved for chemical potentials in the gaps, provided that the corresponding harmonic is non vanishing; for Aubry-André potential, this says that the largest gap persists.…”