A Bloch band based level set method for computing the semiclassical limit of Schrödinger equations (Q1017593)

The authors propose a numerical method based on high-frequency asymptotics for the Schrödinger equation with a rapidly varying potential. The method starts with the computation of a number of Bloch eigenvalues. The density of each Bloch band is calculated from the zero level set in phase space of the derivative of the phase with respect to the slow spatial variable. The method breaks down in the vicinity of caustics.