Matlab codes used to generate figures in the paper 'Cusps in the quench dynamics of a Bloch state', EPL 114, 60001 (2016).
https://iopscience.iop.org/article/10.1209/0295-5075/114/60001
https://iopscience.iop.org/article/10.1209/0295-5075/116/10008
We report some nonsmooth dynamics of a Bloch state in a one-dimensional tight binding model with the periodic boundary condition. After a sudden change of the potential of an arbitrary site, quantities like the survival probability of the particle in the initial Bloch state show cusps periodically, with the period being the Heisenberg time associated with the energy spectrum. This phenomenon is a nonperturbative counterpart of the nonsmooth dynamics observed previously (Zhang J. M. and Haque M., arXiv:1404.4280) in a periodically driven tight binding model. Underlying the cusps is an exactly solvable model, which consists of equally spaced levels extending from
Just run the matlab scripts and you will get all the figures in the paper.