Applied/PDE Seminar: Geneviève Dusson (Université Pierre et Marie Curie (Paris VI)) - 'A Posteriori Error Estimates for Nonlinear Schrödinger Equations'

Event Date: 

Friday, September 23, 2016 - 1:30pm to 2:30pm

Event Location: 

  • South Hall 4607

Joint work with E. Cancès, Y. Maday, B. Stamm, and M. Vohralík.

The determination of electronic structures often requires to solve nonlinear eigenvalue problems. These problems are computationally very costly, and therefore approximations have to be resorted to, among which the chosen model, the chosen discretization, and the chosen (possibly iterative) algorithm. The a priori analysis for the discretization of such problems is quite recent (see e.g. [1]). It proves the convergence and the optimality of the method used. For a given approximation, the a posteriori analysis, once it is performed, provides a guaranteed upper bound on the total error. It may also enable to separate error components stemming from the different sources of approximation and control each of them. This makes possible to iteratively fit the approximation parameters leading to small errors at low computational cost. In this talk, we shall present an a posteriori analysis for a one-dimensional Gross-Pitaevskii type equation in a periodic setting using a planewave discretization [3]. We will provide a computable bound of the error between the exact and approximate solutions. This bound will then be decomposed into two components, each of them depending mainly on one approximation parameter, which are here the number of degrees of freedom and the number of iterations in this iterative algorithm used to solve the problem numerically. Numerical simulations will be presented to illustrate the possibility of balancing the different error components. We shall also present a post-processing on the eigenvalues and eigenfunctions, using a linear perturbation theory based on residual computation [2]. This theoretically and numerically reduces the error significantly both for eigenfunctions and for the energy. Finally, possible extensions of this theory to the 3D Kohn–Sham model will be discussed. References [1] E. Cancès, R. Chakir and Y. Maday, Numerical analysis of the planewave discretization of orbital-free and Kohn-Sham models, M2AN (highlight article) 46 (2012) 341-388. [2] E. Cancès, G. Dusson, Y. Maday, B. Stamm, and M. Vohralík, A perturbation-methodbased post-processing for the planewave discretization of Kohn–Sham models, Journal of Computational Physics 307 (2016) 446-459. [3] G. Dusson, and Y. Maday, A posteriori analysis of a nonlinear GrossPitaevskii-type eigenvalue problem, IMA Journal of Numerical Analysis (2016), drw001.