We study the regularity of stationary and time-dependent solutions to strong
perturbations of the free Schr\"odinger equation on two-dimensional flat tori.
This is achieved by performing a second microlocalization related to the size
of the perturbation and by analysing concentration and nonconcentration
properties at this new scale. In particular, we show that sufficiently accurate
quasimodes can only concentrate on the set of critical points of the average of
the potential along geodesics