Skip to Main content Skip to Navigation
Journal articles

Post-processing of the planewave approximation of Schrödinger equations. Part I: linear operators

Abstract : In this article, we prove a priori error estimates for the perturbation-based post-processing of the plane-wave approximation of Schrödinger equations introduced and tested numerically in previous works [6, 7]. We consider here a Schrödinger operator $H = − 1/2 ∆ + V$ on $L^2 (Ω$), where $Ω$ is a cubic box with periodic boundary conditions, and where V is a multiplicative operator by a regular-enough function V. The quantities of interest are, on the one hand, the ground-state energy defined as the sum of the lowest N eigenvalues of H , and, on the other hand, the ground-state density matrix, that is the spectral projector on the vector space spanned by the associated eigenvectors. Such a problem is central in first-principle molecular simulation, since it corresponds to the so-called linear subproblem in Kohn-Sham density functional theory (DFT). Interpreting the exact eigenpairs of H as perturbations of the numerical eigenpairs obtained by a variational approximation in a plane-wave (i.e. Fourier) basis, we compute first-order corrections for the eigenfunctions, which are turned into corrections on the ground-state density matrix. This allows us to increase the accuracy by a factor proportional to the inverse of the kinetic energy cutoff $E_c^{-1}$ of both the ground-state energy and the ground-state density matrix in Hilbert--Schmidt norm at a low computational extra-cost. Indeed, the computation of the corrections only requires the computation of the residual of the solution in a larger planewave basis and two Fast Fourier Transforms per eigenvalue. Schrödinger operator, perturbation method, planewave approximation, eigenvalue problem, post-processing.
Document type :
Journal articles
Complete list of metadata

Cited literature [23 references]  Display  Hide  Download
Contributor : Geneviève Dusson Connect in order to contact the contributor
Submitted on : Tuesday, March 17, 2020 - 12:53:44 PM
Last modification on : Monday, October 11, 2021 - 10:04:10 AM
Long-term archiving on: : Thursday, June 18, 2020 - 2:34:14 PM


Files produced by the author(s)



Eric Cancès, Geneviève Dusson, Yvon Maday, Benjamin Stamm, Martin Vohralík. Post-processing of the planewave approximation of Schrödinger equations. Part I: linear operators. IMA Journal of Numerical Analysis, Oxford University Press (OUP), 2020, ⟨10.1093/imanum/draa044⟩. ⟨hal-01908039v4⟩



Record views


Files downloads