Gauge-invariant frozen Gaussian approximation method for the Schrödinger equation with periodic potentials

Delgadillo, Ricardo; Lu, Jianfeng; Yang, Xu

doi:10.48550/arxiv.1509.05552

Cited by 3 publications

(5 citation statements)

References 15 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…While the FGA works for general strictly hyperbolic equations, we focus in this paper the case of semiclassical Schrödinger equation with periodic media (1.1). The computational algorithm and numerical results will be presented in a separate paper [5]. The rest of the paper is organized as follows.…”

Section: Introductionmentioning

confidence: 99%

“…While for the analysis, it suffices to assume smooth dependence of u n on ξ (which is possible as the n-th band is separated from the rest of the spectrum), this gauge freedom makes numerical computation nontrivial. We will further address this by designing a gauge-invariant algorithm in a companion paper [5] on the numerical algorithms. Differentiating (2.2) with respect to ξ produces (2.10)…”

Section: Introductionmentioning

confidence: 99%

“…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 [5].for Schwartz class functions, smooth functions, and compactly supported smooth functions, respectively. We will sometimes use subscripts to specify the dependence of a constant on the parameters, for instance, notations like C T to specify the dependence of a constant on a parameter T .…”

mentioning

confidence: 99%

mentioning

confidence: 99%

See 3 more Smart Citations

Frozen Gaussian approximation for high frequency wave propagation in periodic media

Delgadillo

Yang

2018

ASY

Self Cite

View full text Add to dashboard Cite

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 [5].1 2 RICARDO DELGADILLO, JIANFENG LU, AND XU YANG avoid the boundary effects. Therefore the total number of grid points is huge, which usually leads to unaffordable computational cost, especially in high (d > 1) dimensions.An alternative efficient approach is to solve (1.1) asymptotically by the Bloch decomposition and modified WKB methods [3,4,7], which lead to eikonal and transport equations in the semi-classical regime. An advantage of this method is that the computational cost is independent of ε. However, the eikonal equation can develop singularities which make the method break down at caustics. The Gaussian beam method (GBM) [31] was then introduced by Popov to overcome this drawback at caustics. The idea is to allow the phase function to be complex and choose the imaginary part properly so that the solution has a Gaussian profile; see [15-20, 28, 37, 38] for recent developments. Similar ideas can be also found in the Hagedorn wave packet method [8,9]. Unlike the geometric optics based method, the Gaussian beam method allows for accurate computation of wave function around caustics [6,33]. But the problem is that the constructed beam must stay near the geometric rays to maintain accuracy. This becomes a drawback when the solution spreads [23,28,32].The Herman-Kluk propagator [11,21,22] was proposed for Schrödinger equation without the oscillatory periodic background potential. The method was rigorously analyzed in [35,36] and further extended as the frozen Gaussian approximation (FGA) for general high frequency wave propagation in [23][24][25]. The FGA method uses Gaussian functions with fixed widths, instead of using those that might spread over time, to approximate the wave solution. Despite its superficial similarity with the Gaussian beam method, it is different at a fundamental level. FGA is based on phase plane analysis, while GBM is based on the asymptotic solution to a wave equation with Gaussian initial data. In FGA, the solution to the wave equation is approximated by a superposition of Gaussian functions living in phase space, and each function is not necessarily an asymptotic solution, while GBM uses Gaussian functions (called beams) in physical space, with each individual beam being an asymptotic solution to the wave equation. The main advantage of FGA over GBM is that the problem of beam spreading no longer exists.In this paper, we extend FGA for computation of hi...

show abstract

Section: Introductionmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

mentioning

confidence: 99%

mentioning

confidence: 99%

See 2 more Smart Citations

Frozen Gaussian approximation for high frequency wave propagation in periodic media

Delgadillo

Yang

2018

ASY

Self Cite

View full text Add to dashboard Cite

show abstract

“…Note that, since we study here WKB type solutions to the Schrödinger equation, the asymptotic solution we derive is valid only before caustics. If long time validity of the asymptotic solution is desired, one needs to consider instead for example the Gaussian beam methods [7,13], the Wigner functions [14,22], or the frozen Gaussian approximation for periodic media [5,6]. The derivations of Bloch dynamics with Berry phase corrections using these approaches are interesting future directions.…”

Section: Introductionmentioning

confidence: 99%