Mathieu Lewin

Mathieu Lewin, Eric Séré

This paper, devoted to the study of spectral pollution, contains both abstract results and applications to some self-adjoint operators with a gap in their essential spectrum occuring in Quantum Mechanics. First we consider Galerkin basis which respect the decomposition of the ambient Hilbert space into a direct sum $H=PH\oplus(1-P)H$, given by a fixed orthogonal projector $P$, and we localize the polluted spectrum exactly. This is followed by applications to periodic Schrödinger operators (pollution is absent in a Wannier-type basis), and to Dirac operator (several natural decompositions are considered). In the second part, we add the constraint that within the Galerkin basis there is a certain relation between vectors in $PH$ and vectors in $(1-P)H$. Abstract results are proved and applied to several practical methods like the famous "kinetic balance" of relativistic Quantum Mechanics.

Mathieu Lewin, Séverine Paul

The method of choice for describing attractive quantum systems is Hartree-Fock-Bogoliubov (HFB) theory. This is a nonlinear model which allows for the description of pairing effects, the main explanation for the superconductivity of certain materials at very low temperature. This paper is the first study of Hartree-Fock-Bogoliubov theory from the point of view of numerical analysis. We start by discussing its proper discretization and then analyze the convergence of the simple fixed point (Roothaan) algorithm. Following works by Cancès, Le Bris and Levitt for electrons in atoms and molecules, we show that this algorithm either converges to a solution of the equation, or oscillates between two states, none of them being a solution to the HFB equations. We also adapt the Optimal Damping Algorithm of Cancès and Le Bris to the HFB setting and we analyze it. The last part of the paper is devoted to numerical experiments. We consider a purely gravitational system and numerically discover that pairing always occurs. We then examine a simplified model for nucleons, with an effective interaction similar to what is often used in nuclear physics. In both cases we discuss the importance of using a damping algorithm.