On the variable two-step IMEX BDF method for parabolic integro-differential equations with nonsmooth initial data arising in finance