Ricardo Delgadillo

Ricardo Delgadillo, Jianfeng Lu, Xu Yang

We develop a gauge-invariant frozen Gaussian approximation (GIFGA) method for the linear Schrödinger equation (LSE) with periodic potentials in the semiclassical regime. The method generalizes the Herman-Kluk propagator for LSE to the case with periodic media. It provides an efficient computational tool based on asymptotic analysis on phase space and Bloch waves to capture the high-frequency oscillations of the solution. Compared to geometric optics and Gaussian beam methods, GIFGA works in both scenarios of caustics and beam spreading. Moreover, it is invariant with respect to the gauge choice of the Bloch eigenfunctions, and thus avoids the numerical difficulty of computing gauge-dependent Berry phase. We numerically test the method by several one-dimensional examples, in particular, the first order convergence is validated, which agrees with our companion analysis paper [Delgadillo, Lu and Yang, arXiv:1504.08051].

Ricardo Delgadillo, Jianfeng Lu, Xu Yang

Propagation of high-frequency wave in periodic media is a challenging problem due to the existence of multiscale characterized by short wavelength, small lattice constant and large physical domain size. Conventional computational methods lead to extremely expensive costs, especially in high dimensions. In this paper, based on Bloch decomposition and asymptotic analysis in the phase space, we derive the frozen Gaussian approximation for high-frequency wave propagation in periodic media and establish its converge to the true solution. The formulation leads to efficient numerical algorithms, which are presented in a companion paper [Delgadillo, Lu and Yang, arXiv:1509.05552].