Quasiparticle interference patterns in bilayer graphene with trigonal warping

We calculate the form of quasiparticle interference patterns in bilayer graphene within a low-energy description, taking into account perturbatively the trigonal warping terms. We introduce four different types of impurities localized on the A and B sublattices of the first and the second layer, and we obtain closed-form analytical expressions both in real and Fourier spaces for the oscillatory corrections to the local density of states generated by the impurities. Finally, we compare our findings with recent experimental and semianalytical T-matrix results from Joucken et al. [Nano Lett. 21, 7100 (2021)], and we show that there is a very good agreement between our findings and the previous results, as well as with the experimental data.