Long-range asymptotic expansion of the diagonal Born–Oppenheimer correction