Source author record

Eric Cancès

Eric Cancès 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

14works
10topics
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

14 published item(s)

preprint2022arXiv

Photoionization and core resonances from range-separated density-functional theory: General formalism and example of the beryllium atom

We explore the merits of linear-response range-separated time-dependent density-functional theory (TDDFT) for the calculation of photoionization spectra. We consider two variants of range-separated TDDFT, namely the timedependent range-separated hybrid (TDRSH) scheme which uses a global range-separation parameter and the timedependent locally range-separated hybrid (TDLRSH) which uses a local range-separation parameter, and compare with standard time-dependent local-density approximation (TDLDA) and time-dependent Hartree-Fock (TDHF). We show how to calculate photoionization spectra with these methods using the Sternheimer approach formulated in a nonorthogonal B-spline basis set with appropriate frequency-dependent boundary conditions. We illustrate these methods on the photoionization spectrum of the Be atom, focusing in particular on the core resonances. Both the TDRSH and TDLRSH photoionization spectra are found to constitute a large improvement over the TDLDA photoionization spectrum and a more modest improvement over the TDHF photoionization spectrum.

preprint2022arXiv

Photoionization and core resonances from range-separated time-dependent density-functional theory for open-shell states: Example of the lithium atom

We consider the calculations of photoionization spectra and core resonances of open-shell systems using range-separated time-dependent density-functional theory. Specifically, we use the time-dependent range-separated hybrid (TDRSH) scheme, combining a long-range Hartree-Fock (HF) exchange potential and kernel with a short-range potential and kernel from a local density-functional approximation, and the time-dependent locally range-separated hybrid (TDLRSH) scheme, which uses a local range-separation parameter. To efficiently perform the calculations, we formulate a spin-unrestricted linear-response Sternheimer approach in a non-orthogonal B-spline basis set and using appropriate frequency-dependent boundary conditions. We illustrate this approach on the Li atom, which suggests that TDRSH and TDLRSH are adequate simple methods for estimating single-electron photoionization spectra of open-shell systems.

preprint2021arXiv

Multi-center decomposition of molecular densities: a mathematical perspective

The aim of this paper is to analyze from a mathematical perspective some existing schemes to partition a molecular density into several atomic contributions, with a specific focus on Iterative Stockholder Atom (ISA) methods. We provide a unified mathematical framework to describe the latter family of methods and propose a new scheme, named L-ISA (for linear approximation of ISA). We prove several important mathematical properties of the ISA and L-ISA minimization problems and show that the so-called ISA algorithms can be viewed as alternating minimization schemes, which in turn enables us to obtain new convergence results for these numerical methods. Specific mathematical properties of the ISA decomposition for diatomic systems are also presented. We also review the basis-space oriented Distributed Multipole Analysis method, the mathematical formulation of which is also clarified. Different schemes are numerically compared on different molecules and we discuss the advantages and drawbacks of each approach.

preprint2020arXiv

A posteriori error estimation for the non-self-consistent Kohn-Sham equations

We address the problem of bounding rigorously the errors in the numerical solution of the Kohn-Sham equations due to (i) the finiteness of the basis set, (ii) the convergence thresholds in iterative procedures, (iii) the propagation of rounding errors in floating-point arithmetic. In this contribution, we compute fully-guaranteed bounds on the solution of the non-self-consistent equations in the pseudopotential approximation in a plane-wave basis set. We demonstrate our methodology by providing band structure diagrams of silicon annotated with error bars indicating the combined error.

preprint2020arXiv

Guaranteed a posteriori bounds for eigenvalues and eigenvectors: multiplicities and clusters

This paper presents a posteriori error estimates for conforming numerical approximations of eigenvalue clusters of second-order self-adjoint elliptic linear operators with compact resolvent. Given a cluster of eigenvalues, we estimate the error in the sum of the eigenvalues, as well as the error in the eigenvectors represented through the density matrix, i.e., the orthogonal projector on the associated eigenspace. This allows us to deal with degenerate (multiple) eigenvalues within the framework. All the bounds are valid under the only assumption that the cluster is separated from the surrounding smaller and larger eigenvalues; we show how this assumption can be numerically checked. Our bounds are guaranteed and converge with the same speed as the exact errors. They can be turned into fully computable bounds as soon as an estimate on the dual norm of the residual is available, which is presented in two particular cases: the Laplace eigenvalue problem discretized with conforming finite elements, and a Schr{ö}dinger operator with periodic boundary conditions of the form $--$Δ$ + V$ discretized with planewaves. For these two cases, numerical illustrations are provided on a set of test problems.

preprint2017arXiv

Variational projector augmented-wave method: theoretical analysis and preliminary numerical results

In Kohn-Sham electronic structure computations, wave functions have singularities at nuclear positions. Because of these singularities, plane-wave expansions give a poor approximation of the eigenfunctions. In conjunction with the use of pseudo-potentials, the PAW (projector augmented-wave) method circumvents this issue by replacing the original eigenvalue problem by a new one with the same eigenvalues but smoother eigenvectors. Here a slightly different method, called VPAW (variational PAW), is proposed and analyzed. This new method allows for a better convergence with respect to the number of plane-waves. Some numerical tests on an idealized case corroborate this efficiency.

preprint2013arXiv

Greedy algorithms for high-dimensional eigenvalue problems

In this article, we present two new greedy algorithms for the computation of the lowest eigenvalue (and an associated eigenvector) of a high-dimensional eigenvalue problem, and prove some convergence results for these algorithms and their orthogonalized versions. The performance of our algorithms is illustrated on numerical test cases (including the computation of the buckling modes of a microstructured plate), and compared with that of another greedy algorithm for eigenvalue problems introduced by Ammar and Chinesta.

preprint2012arXiv

Non-consistent approximations of self-adjoint eigenproblems: Application to the supercell method

In this article, we introduce a general theoretical framework to analyze non-consistent approximations of the discrete eigenmodes of a self-adjoint operator. We focus in particular on the discrete eigenvalues laying in spectral gaps. We first provide a priori error estimates on the eigenvalues and eigenvectors in the absence of spectral pollution. We then show that the supercell method for perturbed periodic Schrödinger operators falls into the scope of our study. We prove that this method is spectral pollution free, and we derive optimal convergence rates for the planewave discretization method, taking numerical integration errors into account. Some numerical illustrations are provided.

preprint2011arXiv

Periodic Schrödinger operators with local defects and spectral pollution

This article deals with the numerical calculation of eigenvalues of perturbed periodic Schrödinger operators located in spectral gaps. Such operators are encountered in the modeling of the electronic structure of crystals with local defects, and of photonic crystals. The usual finite element Galerkin approximation is known to give rise to spectral pollution. In this article, we give a precise description of the corresponding spurious states. We then prove that the supercell model does not produce spectral pollution. Lastly, we extend results by Lewin and Séré on some no-pollution criteria. In particular, we prove that using approximate spectral projectors enables one to eliminate spectral pollution in a given spectral gap of the reference periodic Schödinger operator.

preprint2010arXiv

Local defects are always neutral in the Thomas-Fermi-von Weiszäcker theory of crystals

The aim of this article is to propose a mathematical model describing the electronic structure of crystals with local defects in the framework of the Thomas-Fermi-von Weizsäcker (TFW) theory. The approach follows the same lines as that used in {\it E. Cancès, A. Deleurence and M. Lewin, Commun. Math. Phys., 281 (2008), pp. 129--177} for the reduced Hartree-Fock model, and is based on thermodynamic limit arguments. We prove in particular that it is not possible to model charged defects within the TFW theory of crystals. We finally derive some additional properties of the TFW ground state electronic density of a crystal with a local defect, in the special case when the host crystal is modelled by a homogeneous medium.

preprint2010arXiv

Numerical analysis of the planewave discretization of some orbital-free and Kohn-Sham models

We provide a priori error estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermi-von Weizsäcker (TFW) model and for the spectral discretization of the Kohn-Sham model, within the local density approximation (LDA). These models allow to compute approximations of the ground state energy and density of molecular systems in the condensed phase. The TFW model is stricly convex with respect to the electronic density, and allows for a comprehensive analysis. This is not the case for the Kohn-Sham LDA model, for which the uniqueness of the ground state electronic density is not guaranteed. Under a coercivity assumption on the second order optimality condition, we prove that for large enough energy cut-offs, the discretized Kohn-Sham LDA problem has a minimizer in the vicinity of any Kohn-Sham ground state, and that this minimizer is unique up to unitary transform. We then derive optimal a priori error estimates for the spectral discretization method.

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

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.

preprint2006arXiv

An efficient sampling algorithm for Variational Monte Carlo

We propose a new algorithm for sampling the $N$-body density $|Ψ({\bf R})|^2/\int_{\mathbb{R}^{3N}} |Ψ|^2$ in the Variational Monte Carlo (VMC) framework. This algorithm is based upon a modified Ricci-Ciccotti discretization of the Langevin dynamics in the phase space $({\bf R},{\bf P})$ improved by a Metropolis acceptation/rejection step. We show through some representative numerical examples (Lithium, Fluorine and Copper atoms, and phenol molecule), that this algorithm is superior to the standard sampling algorithm based on the biased random walk (importance sampling).