A PDE-based fast local level set method
Journal of Computational Physics
On time-splitting spectral approximations for the Schrödinger equation in the semiclassical regime
Journal of Computational Physics
Using K-branch entropy solutions for multivalued geometric optics computations
Journal of Computational Physics
Geometric optics in a phase-space-based level set and Eulerian framework
Journal of Computational Physics
Simplicial isosurfacing in arbitrary dimension and codimension
Journal of Computational Physics
Journal of Computational Physics
SIAM Journal on Scientific Computing
Journal of Scientific Computing
A Bloch band based level set method for computing the semiclassical limit of Schrödinger equations
Journal of Computational Physics
Computing multivalued physical observables for the semiclassical limit of the Schrödinger equation
Journal of Computational Physics
High Order Numerical Methods to Three Dimensional Delta Function Integrals in Level Set Methods
SIAM Journal on Scientific Computing
Hi-index | 31.45 |
The linear Schrodinger equation with periodic potentials is an important model in solid state physics. The most efficient direct simulation using a Bloch decomposition-based time-splitting spectral method [18] requires the mesh size to be O(@e) where @e is the scaled semiclassical parameter. In this paper, we generalize the Gaussian beam method introduced in Jin et al. [23] to solve this problem asymptotically. We combine the technique of Bloch decomposition and the Eulerian Gaussian beam method to arrive at an Eulerian computational method that requires mesh size of O(@e). The accuracy of this method is demonstrated via several numerical examples.