Gaussian Beam Methods for the Schrodinger Equation in the Semi-classical Regime

Shi Jin,Hao Wu,Xu Yang

2008-01-01

Abstract:The solution to the Schrodinger equation is highly oscillatory when the rescaled Planck constant " is small in the semiclassical regime. A direct numerical simulation requires the mesh size to be O("). The Gaussian beam method is an efficient way to solve the high frequency wave equations asymptotically, outperforming geometric optics method in that the Gaussian beam method is accurate even at caustics. In this paper, we solve the Schrodinger equation using both the Lagrangian and Eulerian formulations of the Gaussian beam methods. A new Eulerian Gaussian beam method is developed using the level set method based only on solving the (complex-valued) homogeneous Liouville equations. A major contribution here is that we are able to construct the Hessian matrices of the beams by using the level-set function's first derivatives. This greatly reduces the computational cost in computing the Hessian of the phase function in the Eulerian framework, yielding an Eulerian Gaussian beam method with compu- tational complexity comparable to that of the geometric optics but with a much better accuracy around caustics. We verify through several numerical experiments that our Gaussian beam solutions are good approximations to the Schrodinger solutions even at caustics. We also numerically study the optimal relation be- tween the number of beams and the rescaled Planck constant " in the Gaussian beam summation.