Communications in Mathematical Sciences

Volume 14 (2016)

Number 7

Sobolev and max norm error estimates for Gaussian beam superpositions

Pages: 2037 – 2072



Hailiang Liu (Department of Mathematics, Iowa State University, Ames, Ia., U.S.A.)

Olof Runborg (Department of Mathematics and Swedish e-Science Research Center (SeRC), Royalty Institute of Technology (KTH) Stockholm, Sweden)

Nicolay M. Tanushev (Z-Terra Inc., Houston Texas, U.S.A.)


This work is concerned with the accuracy of Gaussian beam superpositions, which are asymptotically valid high frequency solutions to linear hyperbolic partial differential equations and the Schrödinger equation. We derive Sobolev and max norms estimates for the difference between an exact solution and the corresponding Gaussian beam approximation, in terms of the short wavelength $\varepsilon$. The estimates are performed for the scalar wave equation and the Schrödinger equation. Our result demonstrates that a Gaussian beam superposition with kth order beams converges to the exact solution as $O(\varepsilon^{k / 2 - s})$ in order $s$ Sobolev norms. This result is valid in any number of spatial dimensions and it is unaffected by the presence of caustics in the solution. In max norm, we show that away from caustics the convergence rate is $O(\varepsilon^{\lceil k / 2 \rceil})$ and away from the essential support of the solution, the convergence is spectral in $\varepsilon$. However, in the neighborhood of a caustic point we are only able to show the slower, and dimensional dependent, rate $O(\varepsilon^{(k-n)/2)})$ in $n$ spatial dimensions.


high-frequency wave propagation, error estimates, Gaussian beams, Sobolev norm, max norm

2010 Mathematics Subject Classification

35A35, 35L05, 35L30, 41A60, 58J45

Published 14 September 2016