The nondiagonal matrix element
I no longer understand the numerics we did for circular trajectories. According to quasi-classics, nondiagonal elements should be 0. Perhaps there is an error. A check for that is to plot |\langle +|\sigma_z| -\rangle |
with respect to the radius and perhaps for different offsets.