Source author record

Mathieu Lewin

Mathieu Lewin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

43works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

43 published item(s)

preprint2022arXiv

Coulomb and Riesz gases: The known and the unknown

We review what is known, unknown and expected about the mathematical properties of Coulomb and Riesz gases. Those describe infinite configurations of points in $\mathbb{R}^d$ interacting with the Riesz potential $\pm |x|^{-s}$ (resp. $-\log|x|$ for $s=0$). Our presentation follows the standard point of view of statistical mechanics, but we also mention how these systems arise in other important situations (e.g. in random matrix theory). The main question addressed in the article is how to properly define the associated infinite point process and characterize it using some (renormalized) equilibrium equation. This is largely open in the long range case $s<d$. For the convenience of the reader we give the detail of what is known in the short range case $s>d$. In the last part we discuss phase transitions and mention what is expected on physical grounds.

preprint2020arXiv

The double-power nonlinear Schrödinger equation and its generalizations: uniqueness, non-degeneracy and applications

In this paper we first prove a general result about the uniqueness and non-degeneracy of positive radial solutions to equations of the form $Δu+g(u)=0$. Our result applies in particular to the double power non-linearity where $g(u)=u^q-u^p-μu$ for $p>q>1$ and $μ>0$, which we discuss with more details. In this case, the non-degeneracy of the unique solution $u_μ$ allows us to derive its behavior in the two limits $μ\to0$ and $μ\toμ_*$ where $μ_*$ is the threshold of existence. This gives the uniqueness of energy minimizers at fixed mass in certain regimes. We also make a conjecture about the variations of the $L^2$ mass of $u_μ$ in terms of $μ$, which we illustrate with numerical simulations. If valid, this conjecture would imply the uniqueness of energy minimizers in all cases and also give some important information about the orbital stability of $u_μ$.

preprint2016arXiv

Bose Gases at Positive Temperature and Non-Linear Gibbs Measures

We summarize recent results on positive temperature equilibrium states of large bosonic systems. The emphasis will be on the connection between bosonic grand-canonical thermal states and the (semi-) classical Gibbs measures on one-body quantum states built using the corresponding mean-field energy functionals. An illustrative comparison with the case of "distinguishable" particles (boltzons) is provided.

preprint2016arXiv

Coulomb potentials and Taylor expansions in Time-Dependent Density Functional Theory

We investigate when Taylor expansions can be used to prove the Runge-Gross Theorem, which is at the foundation of Time-Dependent Density Functional Theory (TDDFT). We start with a general analysis of the conditions for the Runge-Gross argument, especially the time-differentiability of the density. The latter should be questioned in the presence of singular (e.g. Coulomb) potentials. Then, we show that a singular potential in a one-body operator considerably decreases the class of time-dependent external potentials to which the original argument can be applied. A two-body singularity has an even stronger impact and an external potential is essentially incompatible with it. For the Coulomb interaction and all reasonable initial many-body states, the Taylor expansion only exists to a finite order, except for constant external potentials. Therefore, high-order Taylor expansions are not the right tool to study atoms and molecules in TDDFT.

preprint2016arXiv

Semi-classical Dirac vacuum polarisation in a scalar field

We study vacuum polarisation effects of a Dirac field coupled to an external scalar field and derive a semi-classical expansion of the regu-larised vacuum energy. The leading order of this expansion is given by a classical formula due to Chin, Lee-Wick and Walecka, for which our result provides the first rigorous proof. We then discuss applications to the non-relativistic large-coupling limit of an interacting system, and to the stability of homogeneous systems.

preprint2015arXiv

A many-body RAGE theorem

We prove a generalized version of the RAGE theorem for N-body quantum systems. The result states that only bound states of systems with $0\leq n\leq N$ particles persist in the long time average. The limit is formulated by means of an appropriate weak topology for many-body systems, which was introduced by the second author in a previous work, and is based on reduced density matrices. This topology is connected to the weak-* topology of states on the algebras of canonical commutation or anti-commutation relations, and we give a formulation of our main result in this setting.

preprint2015arXiv

A note on 2D focusing many-boson systems

We consider a 2D quantum system of $N$ bosons in a trapping potential $|x|^s$, interacting via a pair potential of the form $N^{2β-1} w(N^βx)$. We show that for all $0 \textless{} β\textless{} (s+1)/(s+2)$, the leading order behavior of ground states of the many-body system is described in the large $N$ limit by the corresponding cubic nonlinear Schr{ö}dinger energy functional. Our result covers the focusing case ($w \leq 0$) where even the stability of the many-body system is not obvious. This answers an open question mentioned by X.~Chen and J.~Holmer for harmonic traps ($s=2$). Together with the BBGKY hierarchy approach used by these authors, our result implies the convergence of the many-body quantum dynamics to the focusing NLS equation with harmonic trap for all $0 \textless{} β\textless{} 3/4$.

preprint2015arXiv

Derivation of nonlinear Gibbs measures from many-body quantum mechanics

We prove that nonlinear Gibbs measures can be obtained from the corresponding many-body, grand-canonical, quantum Gibbs states, in a mean-field limit where the temperature T diverges and the interaction behaves as 1/T. We proceed by characterizing the interacting Gibbs state as minimizing a functional counting the free-energy relatively to the non-interacting case. We then perform an infinite-dimensional analogue of phase-space semiclassical analysis, using fine properties of the quantum relative entropy, the link between quantum de Finetti measures and upper/lower symbols in a coherent state basis, as well as Berezin-Lieb type inequalities. Our results cover the measure built on the defocusing nonlinear Schr{ö}dinger functional on a finite interval, as well as smoother interactions in dimensions d\textgreater{}1.

preprint2015arXiv

Fluctuations around Hartree states in the mean-field regime

We consider the dynamics of a large system of N interacting bosons in the mean-field regime where the interaction is of order 1/N. We prove that the fluctuations around the nonlinear Hartree state are generated by an effective quadratic Hamiltonian in Fock space, which is derived from Bogoliubov's approximation. We use a direct method in the N-particle space, which is different from the one based on coherent states in Fock space.

preprint2015arXiv

Improved Lieb-Oxford exchange-correlation inequality with gradient correction

We prove a Lieb-Oxford-type inequality on the indirect part of the Coulomb energy of a general many-particle quantum state, with a lower constant than the original statement but involving an additional gradient correction. The result is similar to a recent inequality of Benguria, Bley and Loss, except that the correction term is purely local, which is more usual in density functional theory. In an appendix, we discuss the connection between the indirect energy and the classical Jellium energy for constant densities. We show that they differ by an explicit shift due to the long range of the Coulomb potential.

preprint2015arXiv

The Crystallization Conjecture: A Review

In this article we describe the crystallization conjecture. It states that, in appropriate physical conditions, interacting particles always place themselves into periodic configurations, breaking thereby the natural translation-invariance of the system. This famous problem is still largely open. Mathematically, it amounts to studying the minima of a real-valued function defined on $\mathbb{R}^{3N}$ where $N$ is the number of particles, which tends to infinity. We review the existing literature and mention several related open problems, of which many have not been thoroughly studied.

preprint2015arXiv

The mean-field approximation and the non-linear Schrödinger functional for trapped Bose gases

We study the ground state of a trapped Bose gas, starting from the full many-body Schr{ö}dinger Hamiltonian, and derive the nonlinear Schr{ö}dinger energy functional in the limit of large particle number, when the interaction potential converges slowly to a Dirac delta function. Our method is based on quantitative estimates on the discrepancy between the full many-body energy and its mean-field approximation using Hartree states. These are proved using finite dimensional localization and a quantitative version of the quantum de Finetti theorem. Our approach covers the case of attractive interactions in the regime of stability. In particular, our main new result is a derivation of the 2D attractive nonlinear Schr{ö}dinger ground state.

preprint2014arXiv

Bogoliubov spectrum of interacting Bose gases

We study the large-N limit of a system of N bosons interacting with a potential of intensity 1/N. When the ground state energy is to the first order given by Hartree's theory, we study the next order, predicted by Bogoliubov's theory. We show the convergence of the lower eigenvalues and eigenfunctions towards that of the Bogoliubov Hamiltonian (up to a convenient unitary transform). We also prove the convergence of the free energy when the system is sufficiently trapped. Our results are valid in an abstract setting, our main assumptions being that the Hartree ground state is unique and non-degenerate, and that there is complete Bose-Einstein condensation on this state. Using our method we then treat two applications: atoms with ''bosonic'' electrons on one hand, and trapped 2D and 3D Coulomb gases on the other hand.

preprint2014arXiv

Remarks on the quantum de Finetti theorem for bosonic systems

The quantum de Finetti theorem asserts that the k-body density matrices of a N-body bosonic state approach a convex combination of Hartree states (pure tensor powers) when N is large and k fixed. In this note we review a construction due to Christandl, Mitchison, König and Renner valid for finite dimensional Hilbert spaces, which gives a quantitative version of the theorem. We first propose a variant of their proof that leads to a slightly improved estimate. Next we provide an alternative proof of an explicit formula due to Chiribella, which gives the density matrices of the constructed state as a function of those of the original state.

preprint2014arXiv

The Hartree equation for infinitely many particles. I. Well-posedness theory

We show local and global well-posedness results for the Hartree equation $$i\partial_tγ=[-Δ+w*ρ_γ,γ],$$ where $γ$ is a bounded self-adjoint operator on $L^2(\R^d)$, $ρ_γ(x)=γ(x,x)$ and $w$ is a smooth short-range interaction potential. The initial datum $γ(0)$ is assumed to be a perturbation of a translation-invariant state $γ_f=f(-Δ)$ which describes a quantum system with an infinite number of particles, such as the Fermi sea at zero temperature, or the Fermi-Dirac and Bose-Einstein gases at positive temperature. Global well-posedness follows from the conservation of the relative (free) energy of the state $γ(t)$, counted relatively to the stationary state $γ_f$. We indeed use a general notion of relative entropy, which allows to treat a wide class of stationary states $f(-Δ)$. Our results are based on a Lieb-Thirring inequality at positive density and on a recent Strichartz inequality for orthonormal functions, which are both due to Frank, Lieb, Seiringer and the first author of this article.

preprint2014arXiv

Uniqueness and non-degeneracy for a nuclear nonlinear Schrödinger equation

We prove the uniqueness and non-degeneracy of positive solutions to a cubic nonlinear Schrödinger (NLS) type equation that describes nucleons. The main difficulty stems from the fact that the mass depends on the solution itself. As an application, we construct solutions to the $σ$--$ω$ model, which consists of one Dirac equation coupled to two Klein-Gordon equations (one focusing and one defocusing).

preprint2013arXiv

Derivation of Hartree's theory for generic mean-field Bose systems

In this paper we provide a novel strategy to prove the validity of Hartree's theory for the ground state energy of bosonic quantum systems in the mean-field regime. For the well-known case of trapped Bose gases, this can be shown using the strong quantum de Finetti theorem, which gives the structure of infinite hierarchies of k-particles density matrices. Here we deal with the case where some particles are allowed to escape to infinity, leading to a lack of compactness. Our approach is based on two ingredients: (1) a weak version of the quantum de Finetti theorem, and (2) geometric techniques for many-body systems. Our strategy does not rely on any special property of the interaction between the particles. In particular, our results cover those of Benguria-Lieb and Lieb-Yau for, respectively, bosonic atoms and boson stars.

preprint2013arXiv

Strichartz inequality for orthonormal functions

We prove a Strichartz inequality for a system of orthonormal functions, with an optimal behavior of the constant in the limit of a large number of functions. The estimate generalizes the usual Strichartz inequality, in the same fashion as the Lieb-Thirring inequality generalizes the Sobolev inequality. As an application, we consider the Schrödinger equation in a time-dependent potential and we show the existence of the wave operator in Schatten spaces.

preprint2013arXiv

The Hartree equation for infinitely many particles. II. Dispersion and scattering in 2D

We consider the nonlinear Hartree equation for an interacting gas containing infinitely many particles and we investigate the large-time stability of the stationary states of the form $f(-Δ)$, describing an homogeneous Fermi gas. Under suitable assumptions on the interaction potential and on the momentum distribution $f$, we prove that the stationary state is asymptotically stable in dimension 2. More precisely, for any initial datum which is a small perturbation of $f(-Δ)$ in a Schatten space, the system weakly converges to the stationary state for large times.

preprint2012arXiv

A Numerical Perspective on Hartree-Fock-Bogoliubov Theory

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.

preprint2012arXiv

Construction of the Pauli-Villars-regulated Dirac vacuum in electromagnetic fields

Using the Pauli-Villars regularization and arguments from convex analysis, we construct solutions to the classical time-independent Maxwell equations in Dirac's vacuum, in the presence of small external electromagnetic sources. The vacuum is not an empty space, but rather a quantum fluctuating medium which behaves as a nonlinear polarizable material. Its behavior is described by a Dirac equation involving infinitely many particles. The quantum corrections to the usual Maxwell equations are nonlinear and nonlocal. Even if photons are described by a purely classical electromagnetic field, the resulting vacuum polarization coincides to first order with that of full Quantum Electrodynamics.

preprint2012arXiv

Dynamical Ionization Bounds for Atoms

We study the long-time behavior of the 3-dimensional repulsive nonlinear Hartree equation with an external attractive Coulomb potential $-Z/|x|$, which is a nonlinear model for the quantum dynamics of an atom. We show that, after a sufficiently long time, the average number of electrons in any finite ball is always smaller than 4Z (respectively 2Z in the radial case). This is a time-dependent generalization of a celebrated result by E.H. Lieb on the maximum negative ionization of atoms in the stationary case. Our proof involves a novel positive commutator argument (based on the cubic weight $|x|^3$) and our findings are reminiscent of the RAGE theorem. In addition, we prove a similar universal bound on the local kinetic energy. In particular, our main result means that, in a weak sense, any solution is attracted to a bounded set in the energy space, whatever the size of the initial datum. Moreover, we extend our main result to Hartree--Fock theory and to the linear many-body Schrödinger equation for atoms.

preprint2012arXiv

Ground state properties of graphene in Hartree-Fock theory

We study the Hartree-Fock approximation of graphene in infinite volume, with instantaneous Coulomb interactions. First we construct its translation-invariant ground state and we recover the well-known fact that, due to the exchange term, the effective Fermi velocity is logarithmically divergent at zero momentum. In a second step we prove the existence of a ground state in the presence of local defects and we discuss some properties of the linear response to an external electric field. All our results are non perturbative.

preprint2012arXiv

On the binding of polarons in a mean-field quantum crystal

We consider a multi-polaron model obtained by coupling the many-body Schrödinger equation for N interacting electrons with the energy functional of a mean-field crystal with a localized defect, obtaining a highly non linear many-body problem. The physical picture is that the electrons constitute a charge defect in an otherwise perfect periodic crystal. A remarkable feature of such a system is the possibility to form a bound state of electrons via their interaction with the polarizable background. We prove first that a single polaron always binds, i.e. the energy functional has a minimizer for N=1. Then we discuss the case of multi-polarons containing two electrons or more. We show that their existence is guaranteed when certain quantized binding inequalities of HVZ type are satisfied.

preprint2011arXiv

Comment on `Solutions to quasi-relativistic multi-configurative Hartree-Fock equations in quantum chemistry', by C. Argaez & M. Melgaard

In a recent paper published in Nonlinear Analysis: Theory, Methods & Applications, C. Argaez and M. Melgaard studied excited states for pseudo-relativistic multi-configuration methods. Their paper follows a previous work of mine in the non-relativistic case (Arch. Rat. Mech. Anal., 171, 2004). The main results of the paper of C. Argaez and M. Melgaard are correct, but the proofs are both wrong and incomplete.

preprint2011arXiv

Derivation of Pekar's Polarons from a Microscopic Model of Quantum Crystals

A polaron is an electron interacting with a polar crystal, which is able to form a bound state by using the distortions of the crystal induced by its own density of charge. In this paper we derive Pekar's famous continuous model for polarons (in which the crystal is replaced by a simple effective Coulomb self-attraction) by studying the macroscopic limit of the reduced Hartree-Fock theory of the crystal. The macroscopic density of the polaron converges to that of Pekar's nonlinear model, with a possibly anisotropic dielectric matrix. The polaron also exhibits fast microscopic oscillations which contribute to the energy at the same order, but whose characteristic length is small compared to the scale of the polaron. These oscillations are described by a simple periodic eigenvalue equation. Our approach also covers multi-polarons composed of several electrons, repelling each other by Coulomb forces.

preprint2011arXiv

Generalised Weyl theorems and spectral pollution in the Galerkin method

We consider a general framework for investigating spectral pollution in the Galerkin method. We show how this phenomenon is characterised via the existence of particular Weyl sequences which are singular in a suitable sense. For a semi-bounded selfadjoint operator A we identify relative compactness conditions on a selfadjoint perturbation B ensuring that the limiting set of spectral pollution of A and B coincide. Our results show that, under perturbation, this limiting set behaves in a similar fashion as the essential spectrum.

preprint2011arXiv

On Singularity formation for the L^2-critical Boson star equation

We prove a general, non-perturbative result about finite-time blowup solutions for the $L^2$-critical boson star equation $i\partial_t u = \sqrt{-Δ+m^2} \, u - (|x|^{-1} \ast |u|^2) u$ in 3 space dimensions. Under the sole assumption that the solution blows up in $H^{1/2}$ at finite time, we show that $u(t)$ has a unique weak limit in $L^2$ and that $|u(t)|^2$ has a unique weak limit in the sense of measures. Moreover, we prove that the limiting measure exhibits minimal mass concentration. A central ingredient used in the proof is a "finite speed of propagation" property, which puts a strong rigidity on the blowup behavior of $u$. As the second main result, we prove that any radial finite-time blowup solution $u$ converges strongly in $L^2$ away from the origin. For radial solutions, this result establishes a large data blowup conjecture for the $L^2$-critical boson star equation, similar to a conjecture which was originally formulated by F. Merle and P. Raphael for the $L^2$-critical nonlinear Schrödinger equation in [CMP 253 (2005), 675-704]. We also discuss some extensions of our results to other $L^2$-critical theories of gravitational collapse, in particular to critical Hartree-type equations.

preprint2010arXiv

Geometric methods for nonlinear many-body quantum systems

Geometric techniques have played an important role in the seventies, for the study of the spectrum of many-body Schrödinger operators. In this paper we provide a formalism which also allows to study nonlinear systems. We start by defining a weak topology on many-body states, which appropriately describes the physical behavior of the system in the case of lack of compactness, that is when some particles are lost at infinity. We provide several important properties of this topology and use them to provide a simple proof of the famous HVZ theorem in the repulsive case. In a second step we recall the method of geometric localization in Fock space as proposed by Dereziński and Gérard, and we relate this tool to our weak topology. We then provide several applications. We start by studying the so-called finite-rank approximation which consists in imposing that the many-body wavefunction can be expanded using finitely many one-body functions. We thereby emphasize geometric properties of Hartree-Fock states and prove nonlinear versions of the HVZ theorem, in the spirit of works of Friesecke. In the last section we study translation-invariant many-body systems comprising a nonlinear term, which effectively describes the interactions with a second system. As an example, we prove the existence of the multi-polaron in the Pekar-Tomasevich approximation, for certain values of the coupling constant.

preprint2010arXiv

Renormalization and asymptotic expansion of Dirac's polarized vacuum

We perform rigorously the charge renormalization of the so-called reduced Bogoliubov-Dirac-Fock (rBDF) model. This nonlinear theory, based on the Dirac operator, describes atoms and molecules while taking into account vacuum polarization effects. We consider the total physical density including both the external density of a nucleus and the self-consistent polarization of the Dirac sea, but no `real' electron. We show that it admits an asymptotic expansion to any order in powers of the physical coupling constant $\alphaph$, provided that the ultraviolet cut-off behaves as $Λ\sim e^{3π(1-Z_3)/2\alphaph}\gg1$. The renormalization parameter $0

preprint2010arXiv

Symmetry Breaking of Relativistic Multiconfiguration Methods in the Nonrelativistic Limit

The multiconfiguration Dirac-Fock method allows to calculate the state of relativistic electrons in atoms or molecules. This method has been known for a long time to provide certain wrong predictions in the nonrelativistic limit. We study in full mathematical details the nonlinear model obtained in the nonrelativistic limit for Be-like atoms. We show that the method with sp+pd configurations in the J=1 sector leads to a symmetry breaking phenomenon in the sense that the ground state is never an eigenvector of L^2 or S^2. We thereby complement and clarify some previous studies.

preprint2010arXiv

The Microscopic Origin of the Macroscopic Dielectric Permittivity of Crystals: A Mathematical Viewpoint

The purpose of this paper is to provide a mathematical analysis of the Adler-Wiser formula relating the macroscopic relative permittivity tensor to the microscopic structure of the crystal at the atomic level. The technical level of the presentation is kept at its minimum to emphasize the mathematical structure of the results. We also briefly review some models describing the electronic structure of finite systems, focusing on density operator based formulations, as well as the Hartree model for perfect crystals or crystals with a defect.

preprint2009arXiv

On Blowup for time-dependent generalized Hartree-Fock equations

We prove finite-time blowup for spherically symmetric and negative energy solutions of Hartree-Fock and Hartree-Fock-Bogoliubov type equations, which describe the evolution of attractive fermionic systems (e. g. white dwarfs). Our main results are twofold: First, we extend the recent blowup result of [Hainzl and Schlein, Comm. Math. Phys. \textbf{287} (2009), 705--714] to Hartree-Fock equations with infinite rank solutions and a general class of Newtonian type interactions. Second, we show the existence of finite-time blowup for spherically symmetric solutions of a Hartree-Fock-Bogoliubov model, where an angular momentum cutoff is introduced. We also explain the key difficulties encountered in the full Hartree-Fock-Bogoliubov theory.

preprint2009arXiv

Strongly correlated phases in rapidly rotating Bose gases

We consider a system of trapped spinless bosons interacting with a repulsive potential and subject to rotation. In the limit of rapid rotation and small scattering length, we rigorously show that the ground state energy converges to that of a simplified model Hamiltonian with contact interaction projected onto the Lowest Landau Level. This effective Hamiltonian models the bosonic analogue of the Fractional Quantum Hall Effect (FQHE). For a fixed number of particles, we also prove convergence of states; in particular, in a certain regime we show convergence towards the bosonic Laughlin wavefunction. This is the first rigorous justification of the effective FQHE Hamiltonian for rapidly rotating Bose gases. We review previous results on this effective Hamiltonian and outline open problems.

preprint2009arXiv

The Dielectric Permittivity of Crystals in the reduced Hartree-Fock approximation

In a recent article (Cancès, Deleurence and Lewin, Commun. Math. Phys., 281 (2008), pp. 129-177), we have rigorously derived, by means of bulk limit arguments, a new variational model to describe the electronic ground state of insulating or semiconducting crystals in the presence of local defects. In this so-called reduced Hartree-Fock model, the ground state electronic density matrix is decomposed as $γ= γ^0_{\rm per} + Q_{ν,ε_{\rm F}}$, where $γ^0_{\rm per}$ is the ground state density matrix of the host crystal and $Q_{ν,ε_{\rm F}}$ the modification of the electronic density matrix generated by a modification $ν$ of the nuclear charge of the host crystal, the Fermi level $ε_{\rm F}$ being kept fixed. The purpose of the present article is twofold. First, we study more in details the mathematical properties of the density matrix $Q_{ν,ε_{\rm F}}$ (which is known to be a self-adjoint Hilbert-Schmidt operator on $L^2(\R^3)$). We show in particular that if $\int_{\RR^3} ν\neq 0$, $Q_{ν,ε_{\rm F}}$ is not trace-class. Moreover, the associated density of charge is not in $L^1(\R^3)$ if the crystal exhibits anisotropic dielectric properties. These results are obtained by analyzing, for a small defect $ν$, the linear and nonlinear terms of the resolvent expansion of $Q_{ν,ε_{\rm F}}$. Second, we show that, after an appropriate rescaling, the potential generated by the microscopic total charge (nuclear plus electronic contributions) of the crystal in the presence of the defect, converges to a homogenized electrostatic potential solution to a Poisson equation involving the macroscopic dielectric permittivity of the crystal. This provides an alternative (and rigorous) derivation of the Adler-Wiser formula.

preprint2008arXiv

Spectral Pollution and How to Avoid It (With Applications to Dirac and Periodic Schrödinger Operators)

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.