### Principal investigators

Prof. Dr. Guido Schneider | (7/2015 - 6/2019) | |

JProf. Dr. Katharina Schratz | (7/2015 - 6/2019) |

### Project summary

The aim of this project is to analyze the Klein-Gordon-Zakharov system in high plasma frequency and (simultaneous) subsonic regimes. Due to the highly oscillatory nature of the problem classical numerical schemes break down, in particular in certain singular limits. Therefore, we need to develop new analytic techniques and efficient numerical schemes to overcome severe step size restrictions and huge computational costs.

## Numerical challenges of highly oscillatory problems

If a solution of the underlying equation is highly oscillatory, it becomes very challenging for numerical methods to resolve the oscillations and to yield a good numerical approximation. In order to compute a numerical approximation of a highly oscillatory function (see Figure 1 below) we have to discretize the interval with a finite number of points. At these grid points we compute approximations of the function values.

Afterwards we construct our numerical solution by interpolating between the approximated values. The figure shows an example where the approximation fails completely even though we evaluate the function values exactly at the discrete (finite) points. We can improve the approximation by introducing a finer grid, which, however, requires much more computational and memory costs and therefore will not yield an efficient numerical scheme. This makes the efficient approximation of highly oscillatory differential equations an ongoing challenge in numerical analysis.

### The Klein-Gordon-Zakharov system

The starting point of this project was the Klein Gordon Zakhavov (KGZ) system \[ \begin{equation*} c^{-2}\partial_{tt}^2 z = \Delta z - c^2 z \pm nz, \\ \alpha^{-2}\partial_{tt}^2 n = \Delta n + \Delta \vert z\vert^2\quad\,\,\, \end{equation*} \] in certain high frequency regimes, i.e., for large \(c\) and/or large \(\alpha\), where the solution becomes highly oscillatory in time (see the Figure 2 below).

We are interested in a robust analytic and numerical description of the KGZ system for such values of \(c\) and \(\alpha\). Such a description can be obtained via the limit systems for \(c\to\infty\) and/or \(\alpha\to\infty\). This analysis is particularly of great interest from a numerical point of view: Correctly resolving this highly oscillatory behavior in the exact solution in these regimes is numerically very delicate.

Severe time step restrictions need to be imposed and this requires a great deal of computational costs. An asymptotic ansatz allows to construct efficient time integrators as the highly oscillatory parts in the exact solution could be filtered out explicitly. Thus, the numerical task could be reduced to the time integration of the corresponding non-oscillatory limit systems, which can be carried out very efficiently without any additional time step restriction.

The construction of such robust, efficient time integrators involves a careful analysis of the numerical time integrators for the corresponding limit systems and rigorous error estimates between true solutions of the original system and the approximations constructed via the limit system.

The KGZ system turned out to be a well chosen prototype problem since the tools developed for the KGZ system can be used to establish similar results for a wider class of problems. Moreover, extending the original question, uniformly accurate methods for KGZ like systems have been developed, i.e., the numerical costs for such schemes are independent of the \(c\) and \(\alpha\) and not only valid for \(c\) large and/or \(\alpha\) large.

- Analytic results that the various limit systems make correct predictions about the dynamics for \(c\) large and/or \(\alpha\) large can be found in the

**References**:

[1] Markus Daub, Guido Schneider, and Katharina Schratz. From the Klein-Gordon-Zakharov system to the Klein-Gordon equation.*Math. Methods Appl. Sci.*, 39(18):5371--5380, 2016.

[2] Simon Baumstark, Guido Schneider, Katharina Schratz, and Dominik Zimmermann. Effective slow dynamics models for a class of dispersive systems. Preprint 2018.

In the first paper it turned out that limit systems in the limit \(c\) fixed and \(\alpha\to\infty\) depends whether the KGZ system is posed on \(\mathbb{R}^3\) or \(\mathbb{T}^3\). The abstract approximation theorem in the second paper applies to a number of semilinear systems, such as the Dirac-Klein-Gordon system, the Klein-Gordon-Zakharov system, and a mean field polaron model. It extracts the common features of scattered results in the literature, but also gains an approximation result for the Dirac-Klein-Gordon system which has not been documented in the literature before. The abstract approximation theorem is sharp in the sense that there exists a quasilinear system of the same structure, namely the Zakharov system with a 'wrong' sign in the nonlinearity, where the regular limit system, namely the NLS equation, makes wrong predictions. - We propose asymptotic consistent exponential-type integrators for the Klein-Gordon-Schrödinger system. This novel class of integrators allows us to solve the system from slowly varying relativistic up to challenging highly oscillatory non-relativistic regimes without any step size restriction. In particular, our first- and second-order exponential-type integrators are asymptotically consistent in the sense of asymptotically converging to the corresponding decoupled free Schröodinger limit system.

**Reference**:

[3] Simon Baumstark, Georgia Kokkala, and Katharina Schratz. Asymptotic consistent exponential-type integrators for Klein-Gordon-Schrödinger systems from relativistic to non-relativistic regimes.*ETNA, Electron. Trans. Numer. Anal.*, 48:63--80, 2018. - We introduce efficient and robust exponential-type integrators for Klein-Gordon equations which resolve the solution in the relativistic regime as well as in the highly-oscillatory nonrelativistic regime without any step-size restriction under the same regularity assumptions on the initial data required for the integration of the corresponding nonlinear Schrödinger (NLS) limit system. In contrast to previous works we do not employ any asymptotic/multiscale expansion of the solution. This allows us to derive uniform convergent schemes under far weaker regularity assumptions on the exact solution. In addition, the newly derived first- and second-order exponential-type integrators converge to the classical Lie, respectively, Strang splitting in the nonlinear Schrödinger limit.

**Reference**:

[4] Simon Baumstark, Erwan Faou, and Katharina Schratz. Uniformly accurate exponential-type integrators for Klein-Gordon equations with asymptotic convergence to the classical NLS splitting.*Math. Comput.*, 87(311):1227--1254, 2018. - We present a novel class of oscillatory integrators for the KGZ system which are uniformly accurate with respect to the plasma frequency \(c\). Convergence holds from the slowly-varying low-plasma up to the highly-oscillatory high-plasma frequency regimes without any step size restriction and, especially, uniformly in \(c\). The introduced scheme is moreover asymptotic consistent and approximates the solutions of the corresponding Zakharov limit system in the high-plasma frequency limit \(c\to\infty\). We establish rigorous error estimates for the the introduced oscillatory integrator and numerically underline its uniform convergence property. The derivation of uniformly accurate methods for the KGZ system in the subsonic limit regime (\(\alpha\to\infty\)) and also in the simultaneous limit regimes (\((c,\alpha)\to\infty\)) will be the subject of future research.

**Reference**:

[5] Simon Baumstark and Katharina Schratz. Uniformly accurate oscillatory integrators for the Klein-Gordon-Zakharov system from low- to high-plasma frequency regimes. https://arxiv.org/abs/1803.09126.

## Publications

- Splitting methods for nonlinear Dirac equations with Thirring type interaction in the nonrelativistic limit regime,

P. Krämer, K. Schratz, X. Zhao,

*J. Computat. Appl. Math.*, pp. 112494, online first, September 2019 – BibTeX - Uniformly accurate exponential-type integrators for Klein-Gordon equations with asymptotic convergence to the classical NLS splitting,

S. Baumstark, E. Faou, K. Schratz,

*Math. Comput.*, vol. 87, no. 311, pp. 1227-1254, May 2018 – BibTeX - Asymptotic consistent exponential-type integrators for Klein-Gordon-Schrödinger systems from relativistic to non-relativistic regimes,

S. Baumstark, G. Kokkala, K. Schratz,

*Electron. Trans. Numer. Anal.*, vol. 48, pp. 63-80, March 2018 – BibTeX - Efficient time integration of the Maxwell–Klein–Gordon equation in the non-relativistic limit regime,

P. Krämer, K. Schratz,

*J. Computat. Appl. Math.*, vol. 316, pp. 247-259, May 2017 – BibTeX - From the Klein–Gordon–Zakharov system to the Klein–Gordon equation,

M. Daub, G. Schneider, K. Schratz,

*Math. Methods Appl. Sci.*, vol. 23, issue 26, pp. 33044–33064, December 2016 – BibTeX

## Preprints

- Asymptotic preserving trigonometric integrators for the quantum Zakharov system,

S. Baumstark, K. Schratz, January 2019 – BibTeX - A model for the periodic water wave problem and its long wave amplitude equations,

R. Bauer, P. Cummings, G. Schneider, November 2018 – BibTeX - The KdV approximation for a system with unstable resonances,

G. Schneider, November 2018 – BibTeX - Effective slow dynamics models for a class of dispersive systems,

S. Baumstark, G. Schneider. K. Schratz, D. Zimmermann, November 2018 – BibTeX - Uniformly accurate oscillatory integrators for the Klein-Gordon-Zakharov system from low- to high-plasma frequency regimes,

S. Baumstark, K. Schratz, November 2018 – BibTeX - Splitting methods for nonlinear Dirac equations with thirring type interaction in the nonrelativistic limit regime,

P. Krämer, K. Schratz, X. Zhao, November 2018 – BibTeX

## Theses

- Uniformly accurate methods for Klein-Gordon type equations,

Simon Baumstark, Dissertation, July 2018, BibTeX