Search

In this article, we derive fully guaranteed error bounds for the energy of convex nonlinear mean-field models. These results apply in particular to Kohn-Sham equations with convex density functionals, which includes the reduced Hartree-Fock (rHF) model, as well as the Kohn-Sham model with exact exchange-density functional (which is unfortunately not explicit and therefore not usable in practice). We then decompose the obtained bounds into two parts, one depending on the chosen discretization and one depending on the number of iterations performed in the self-consistent algorithm used to solve the nonlinear eigenvalue problem, paving the way for adaptive refinement strategies. The accuracy of the bounds is demonstrated on a series of test cases, including a Silicon crystal and an Hydrogen Fluoride molecule simulated with the rHF model and discretized with planewaves. We also show that, although not anymore guaranteed, the error bounds remain very accurate for a Silicon crystal simulated with the Kohn-Sham model using nonconvex exchangecorrelation functionals of practical interest.

We provide a mathematical analysis of the Dynamical Mean-Field Theory, a celebrated representative of a class of approximations in quantum mechanics known as embedding methods. We start by a pedagogical and self-contained mathematical formulation of the Dynamical Mean-Field Theory equations for the finite Hubbard model. After recalling the definition and properties of one-body time-ordered Green's functions and self-energies, and the mathematical structure of the Hubbard and Anderson impurity models, we describe a specific impurity solver, namely the Iterated Perturbation Theory solver, which can be conveniently formulated using Matsubara's Green's functions. Within this framework, we prove under certain assumptions that the Dynamical Mean-Field Theory equations admit a solution for any set of physical parameters. Moreover, we establish some properties of the solution(s)., Comment: 40 pages, 5 figures

In this study, we analyze various Iterative Stockholder Analysis (ISA) methods for molecular density partitioning, focusing on the numerical performance of the recently proposed Linear approximation of Iterative Stockholder Analysis model (LISA) [J. Chem. Phys. 156, 164107 (2022)]. We first provide a systematic derivation of various iterative solvers to find the unique LISA solution. In a subsequent systematic numerical study, we evaluate their performance on 48 organic and inorganic, neutral and charged molecules and also compare LISA to two other well-known ISA variants: the Gaussian Iterative Stockholder Analysis (GISA) and Minimum Basis Iterative Stockholder analysis (MBIS). The study reveals that LISA-family methods can offer a numerically more efficient approach with better accuracy compared to the two comparative methods. Moreover, the well-known issue with the MBIS method, where atomic charges obtained for negatively charged molecules are anomalously negative, is not observed in LISA-family methods. Despite the fact that LISA occasionally exhibits elevated entropy as a consequence of the absence of more diffuse basis functions, this issue can be readily mitigated by incorporating additional or integrating supplementary basis functions within the LISA framework. This research provides the foundation for future studies on the efficiency and chemical accuracy of molecular density partitioning schemes.

We explore Riemannian optimization methods for Restricted-Open-shell Hartree-Fock (ROHF) and Complete Active Space Self-Consistent Field (CASSCF) methods. After showing that ROHF and CASSCF can be reformulated as optimization problems on so-called flag manifolds, we review Riemannian optimization basics and their application to these specific problems. We compare these methods to traditional ones and find robust convergence properties without fine-tuning of numerical parameters. Our study suggests Riemannian optimization as a valuable addition to orbital optimization for ROHF and CASSCF, warranting further investigation., Comment: 21 pages

This paper investigates the mathematical properties of independent-electron models for twisted bilayer graphene by examining the density-of-states of corresponding single-particle Hamiltonians using tools from semiclassical analysis. This study focuses on a specific atomic-scale Hamiltonian $H_{d,\theta}$ constructed from Density-Functional Theory, and a family of moir\'e-scale Hamiltonians $H_{d,K,\theta}^{\rm eff}$ containing the Bistritzer-MacDonald model. The parameter $d$ represents the interlayer distance, and $\theta$ the twist angle. It is shown that the density-of-states of $H_{d,\theta}$ and $H_{d,K,\theta}^{\rm eff}$ admit asymptotic expansions in the twist angle parameter $\epsilon:=\sin(\theta/2)$. The proof relies on a twisted version of the Weyl calculus and a trace formula for an exotic class of pseudodifferential operators suitable for the study of twisted 2D materials. We also show that the density-of-states of $H_{d,\theta}$ admits an asymptotic expansion in $\eta:=\tan(\theta/2)$ and comment on the differences between the expansions in $\epsilon$ and $\eta$., Comment: 58 pages, 3 figures

This article provides the first mathematical analysis of the Density Matrix Embedding Theory (DMET) method. We prove that, under certain assumptions, (i) the exact ground-state density matrix is a fixed-point of the DMET map for non-interacting systems, (ii) there exists a unique physical solution in the weakly-interacting regime, and (iii) DMET is exact at first order in the coupling parameter. We provide numerical simulations to support our results and comment on the physical meaning of the assumptions under which they hold true. We show that the violation of these assumptions may yield multiple solutions of the DMET equations. We moreover introduce and discuss a specific N-representability problem inherent to DMET., Comment: 40 pages, 11 figures

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., Comment: arXiv admin note: substantial text overlap with arXiv:2203.05951

Response calculations in density functional theory aim at computing the change in ground-state density induced by an external perturbation. At finite temperature these are usually performed by computing variations of orbitals, which involve the iterative solution of potentially badly-conditioned linear systems, the Sternheimer equations. Since many sets of variations of orbitals yield the same variation of density matrix this involves a choice of gauge. Taking a numerical analysis point of view we present the various gauge choices proposed in the literature in a common framework and study their stability. Beyond existing methods we propose a new approach, based on a Schur complement using extra orbitals from the self-consistent-field calculations, to improve the stability and efficiency of the iterative solution of Sternheimer equations. We show the success of this strategy on nontrivial examples of practical interest, such as Heusler transition metal alloy compounds, where savings of around 40% in the number of required cost-determining Hamiltonian applications have been achieved.

In solid state physics, electronic properties of crystalline materials are often inferred from the spectrum of periodic Schr\"odinger operators. As a consequence of Bloch's theorem, the numerical computation of electronic quantities of interest involves computing derivatives or integrals over the Brillouin zone of so-called energy bands, which are piecewise smooth, Lipschitz continuous periodic functions obtained by solving a parametrized elliptic eigenvalue problem on a Hilbert space of periodic functions. Classical discretization strategies for resolving these eigenvalue problems produce approximate energy bands that are either non-periodic or discontinuous, both of which cause difficulty when computing numerical derivatives or employing numerical quadrature. In this article, we study an alternative discretization strategy based on an ad hoc operator modification approach. While specific instances of this approach have been proposed in the physics literature, we introduce here a systematic formulation of this operator modification approach. We derive a priori error estimates for the resulting energy bands and we show that these bands are periodic and can be made arbitrarily smooth (away from band crossings) by adjusting suitable parameters in the operator modification approach. Numerical experiments involving a toy model in 1D, graphene in 2D, and silicon in 3D validate our theoretical results and showcase the efficiency of the operator modification approach., Comment: 40 pages, 12 figures

In this article, we propose general criteria to construct optimal atomic centered basis sets in quantum chemistry. We focus in particular on two criteria, one based on the ground-state one-body density matrix of the system and the other based on the ground-state energy. The performance of these two criteria are then numerically tested and compared on a parametrized eigenvalue problem, which corresponds to a one-dimensional toy version of the ground-state dissociation of a diatomic molecule.

We provide a formal derivation of a reduced model for twisted bilayer graphene (TBG) from Density Functional Theory. Our derivation is based on a variational approximation of the TBG Kohn-Sham Hamiltonian and asymptotic limit techniques. In contrast with other approaches, it does not require the introduction of an intermediate tight-binding model. The so-obtained model is similar to that of the Bistritzer-MacDonald (BM) model but contains additional terms. Its parameters can be easily computed from Kohn-Sham calculations on single-layer graphene and untwisted bilayer graphene with different stackings. It allows one in particular to estimate the parameters $w_{\rm AA}$ and $w_{\rm AB}$ of the BM model from first-principles. The resulting numerical values, namely $w_{\rm AA}= w_{\rm AB} \simeq 126$ meV for the experimental interlayer mean distance are in good agreement with the empirical values $w_{\rm AA}= w_{\rm AB}=110$ meV obtained by fitting to experimental data. We also show that if the BM parameters are set to $w_{\rm AA}= w_{\rm AB} \simeq 126$ meV, the BM model is an accurate approximation of our reduced model.

This paper is concerned with the numerical analysis of linear and nonlinear Schr{\"o}dinger equations with analytic potentials. While the regularity of the potential (and the source term when there is one) automatically conveys to the solution in the linear cases, this is no longer true in general in the nonlinear case. We also study the rate of convergence of the planewave (Fourier) discretization method for computing numerical approximations of the solution.

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.

We consider the homogenization at second-order in $\varepsilon$ of $\mathbb{L}$-periodic Schr\"odinger operators with rapidly oscillating potentials of the form $H^\varepsilon =-\Delta + \varepsilon^{-1} v(x,\varepsilon^{-1}x ) + W(x)$ on $L^2(\mathbb{R}^d)$, where $\mathbb{L}$ is a Bravais lattice of $\mathbb{R}^d$, $v$ is $(\mathbb{L} \times \mathbb{L})$-periodic, $W$ is $\mathbb{L}$-periodic, and $\varepsilon \in \mathbb{N}^{-1}$. We treat both the linear equation with fixed right-hand side and the eigenvalue problem, as well as the case of physical observables such as the integrated density of states. We illustrate numerically that these corrections to the homogenized solution can significantly improve the first-order ones, even when $\varepsilon$ is not small.

We propose accurate computable error bounds for quantities of interest in plane-wave electronic structure calculations, in particular ground-state density matrices and energies, and interatomic forces. These bounds are based on an estimation of the error in terms of the residual of the solved equations, which is then efficiently approximated with computable terms. After providing coarse bounds based on an analysis of the inverse Jacobian, we improve on these bounds by solving a linear problem in a small dimension that involves a Schur complement. We numerically show how accurate these bounds are on a few representative materials, namely silicon, gallium arsenide and titanium dioxide.

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.

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{\"o}dinger operator with periodic boundary conditions of the form $--$\Delta$ + V$ discretized with planewaves. For these two cases, numerical illustrations are provided on a set of test problems.

We extend a method (E. Canc\`es and L.R. Scott, SIAM J. Math. Anal., 50, 2018, 381--410) to compute more terms in the asymptotic expansion of the van der Waals attraction between two hydrogen atoms. These terms are obtained by solving a set of modified Slater-Kirkwood partial differential equations. The accuracy of the method is demonstrated by numerical simulations and comparison with other methods from the literature. It is also shown that the scattering states of the hydrogen atom, that are the states associated with the continuous spectrum of the Hamiltonian, have a major contribution to the C$_6$ coefficient of the van der Waals expansion.

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.

This article is concerned with the numerical solution of subspace optimization problems, consisting of minimizing a smooth functional over the set of orthogonal projectors of fixed rank. Such problems are encountered in particular in electronic structure calculation (Hartree-Fock and Kohn-Sham Density Functional Theory -DFT- models). We compare from a numerical analysis perspective two simple representatives, the damped self-consistent field (SCF) iterations and the gradient descent algorithm, of the two classes of methods competing in the field: SCF and direct minimization methods. We derive asymptotic rates of convergence for these algorithms and analyze their dependence on the spectral gap and other properties of the problem. Our theoretical results are complemented by numerical simulations on a variety of examples, from toy models with tunable parameters to realistic Kohn-Sham computations. We also provide an example of chaotic behavior of the simple SCF iterations for a nonquadratic functional.

We consider independent electrons in a periodic crystal in their ground state, and turn on a uniform electric field at some prescribed time. We rigorously define the current per unit volume and study its properties using both linear response and adiabatic theory. Our results provide a unified framework for various phenomena such as the quantification of Hall conductivity of insulators with broken time-reversibility, the ballistic regime of electrons in metals, Bloch oscillations in the long-time response of metals, and the static conductivity of graphene. We identify explicitly the regime in which each holds.

We study by means of Monte-Carlo numerical simulations the resistance of two-dimensional random percolating networks of stick, widthless nanowires. We use the multi-nodal representation (MNR) to model a nanowire network as a graph. We derive numerically from this model the expression of the total resistance as a function of all meaningful parameters, geometrical and physical, over a wide range of variation for each. We justify our choice of non-dimensional variables applying Buckingham $\pi-$theorem. The effective resistance of 2D random percolating networks of nanowires is found to write as $R^{eq}(\rho,R_c,R_{m,w})=A\left(N,\frac{L}{l^{*}}\right) \rho l^* + B\left(N,\frac{L}{l^{*}}\right) R_c+C\left(N,\frac{L}{l^*} \right) R_{m,w}$ where $N$, $\frac{L}{l^{*}}$ are the geometrical parameters (number of wires, aspect ratio of electrode separation over wire length) and $\rho$, $R_c$, $R_{m,w}$ are the physical parameters (nanowire linear resistance per unit length, nanowire/nanowire contact resistance, metallic electrode/nanowire contact resistance). The dependence of the resistance on the geometry of the network, one the one hand, and on the physical parameters (values of the resistances), on the other hand, is thus clearly separated thanks to this expression, much simpler than the previously reported analytical expressions., Comment: 16 pages, 4 figures, regular article

This contribution is the numerically oriented companion article of the work [E. Canc\`es, V. Ehrlacher, F. Legoll, B. Stamm and S. Xiang, arxiv preprint 1807.05131]. We focus here on the numerical resolution of the embedded corrector problem introduced in [E. Canc\`es, V. Ehrlacher, F. Legoll and B. Stamm, CRAS 2015; E. Canc\`es, V. Ehrlacher, F. Legoll, B. Stamm and S. Xiang, arxiv preprint 1807.05131] in the context of homogenization of diffusion equations. Our approach consists in considering a corrector-type problem, posed on the whole space, but with a diffusion matrix which is constant outside some bounded domain. In [E. Canc\`es, V. Ehrlacher, F. Legoll, B. Stamm and S. Xiang, arxiv preprint 1807.05131], we have shown how to define three approximate homogenized diffusion coefficients on the basis of the embedded corrector problems. We have also proved that these approximations all converge to the exact homogenized coefficients when the size of the bounded domain increases. We show here that, under the assumption that the diffusion matrix is piecewise constant, the corrector problem to solve can be recast as an integral equation. In case of spherical inclusions with isotropic materials, we explain how to efficiently discretize this integral equation using spherical harmonics, and how to use the fast multipole method (FMM) to compute the resulting matrix-vector products at a cost which scales only linearly with respect to the number of inclusions. Numerical tests illustrate the performance of our approach in various settings.

Studying the electronic structure of defects in materials is an important subject in condensed matter physics. From a mathematical point of view, nonlinear mean-field models of localized defects in insulators are well understood. We present here a mean-field model to study a particular instance of extended defects in metals. These extended defects typically correspond to taking out a slab of finite width in the three-dimensional homogeneous electron gas. We work in the framework of the reduced Hartree-Fock model with either Yukawa or Coulomb interactions. Using techniques developed in~[Frank2011, Frank2013] to study local perturbations of the free-electron gas, we show that our model admits minimizers, and that Yukawa ground state energies and density matrices converge to ground state Coulomb energies and density matrices as the Yukawa parameter tends to zero. We moreover present numerical simulations where we observe Friedel oscillations in the total electronic density.

This article is the first part of a two-fold study, the objective of which is the theoretical analysis and numerical investigation of new approximate corrector problems in the context of stochastic homogenization. We present here three new alternatives for the approximation of the homogenized matrix for diffusion problems with highly-oscillatory coefficients. These different approximations all rely on the use of an embedded corrector problem (that we previously introduced in [Canc\`es, Ehrlacher, Legoll and Stamm, C. R. Acad. Sci. Paris, 2015]), where a finite-size domain made of the highly oscillatory material is embedded in a homogeneous infinite medium whose diffusion coefficients have to be appropriately determined. The motivation for considering such embedded corrector problems is made clear in the companion article [Canc\`es, Ehrlacher, Legoll, Stamm and Xiang, in preparation], where a very efficient algorithm is presented for the resolution of such problems for particular heterogeneous materials. In the present article, we prove that the three different approximations we introduce converge to the homogenized matrix of the medium when the size of the embedded domain goes to infinity.

As a consequence of Bloch's theorem, the numerical computation of the fermionic ground state density matrices and energies of periodic Schrodinger operators involves integrals over the Brillouin zone. These integrals are difficult to compute numerically in metals due to discontinuities in the integrand. We perform an error analysis of several widely-used quadrature rules and smearing methods for Brillouin zone integration. We precisely identify the assumptions implicit in these methods and rigorously prove error bounds. Numerical results for two-dimensional periodic systems are also provided. Our results shed light on the properties of these numerical schemes, and provide guidance as to the appropriate choice of numerical parameters.

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., Comment: 41 pages, 6 figures

In this article, we consider the extended Kohn-Sham model for atoms subjected to cylindrically-symmetric external potentials. The variational approximation of the model and the construction of appropriate discretization spaces are detailed together with the algorithm to solve the discretized Kohn-Sham equations used in our code. Using this code, we compute the occupied and unoccupied energy levels of all the atoms of the first four rows of the periodic table for the reduced Hartree-Fock (rHF) and the extended Kohn-Sham X$\alpha$ models. These results allow us to test numerically the assumptions on the negative spectra of atomic rHF and Kohn-Sham Hamiltonians used in our previous theoretical works on density functional perturbation theory and pseudopotentials. Interestingly, we observe accidental degeneracies between s and d shells or between p and d shells at the Fermi level of some atoms. We also consider the case of an atom subjected to a uniform electric-field. For various magnitudes of the electric field, we compute the response of the density of the carbon atom confined in a large ball with Dirichlet boundary conditions, and we check that, in the limit of small electric fields, the results agree with the ones obtained with first-order density functional perturbation theory.

It is often claimed that error cancellation plays an essential role in quantum chemistry and first-principle simulation for condensed matter physics and materials science. Indeed, while the energy of a large, or even medium-size, molecular system cannot be estimated numerically within chemical accuracy (typically 1 kcal/mol or 1 mHa), it is considered that the energy difference between two configurations of the same system can be computed in practice within the desired accuracy. The purpose of this paper is to provide a quantitative study of discretization error cancellation. The latter is the error component due to the fact that the model used in the calculation (e.g. Kohn-Sham LDA) must be discretized in a finite basis set to be solved by a computer. We first report comprehensive numerical simulations performed with Abinit on two simple chemical systems, the hydrogen molecule on the one hand, and a system consisting of two oxygen atoms and four hydrogen atoms on the other hand. We observe that errors on energy differences are indeed significantly smaller than errors on energies, but that these two quantities asymptotically converge at the same rate when the energy cut-off goes to infinity. We then analyze a simple one-dimensional periodic Schr\"odinger equation with Dirac potentials, for which analytic solutions are available. This allows us to explain the discretization error cancellation phenomenon on this test case with quantitative mathematical arguments.

We give an exact formulation for the transport coefficients of incommensurate two-dimensional atomic multilayer systems in the tight-binding approximation. This formulation is based upon the C* algebra framework introduced by Bellissard and collaborators to study aperiodic solids (disordered crystals, quasicrystals, and amorphous materials), notably in the presence of magnetic fields (quantum Hall effect). We also present numerical approximations and test our methods on a one-dimensional incommensurate bilayer system.

We propose an algorithm to determine Maximally Localized Wannier Functions (MLWFs). This algorithm, based on recent theoretical developments, does not require any physical input such as initial guesses for the Wannier functions, unlike popular schemes based on the projection method. We discuss how the projection method can fail on fine grids when the initial guesses are too far from MLWFs. We demonstrate that our algorithm is able to find localized Wannier functions through tests on two-dimensional systems, simplified models of semiconductors, and realistic DFT systems by interfacing with the Wannier90 code. We also test our algorithm on the Haldane and Kane-Mele models to examine how it fails in the presence of topological obstructions.

Molecular dynamics (MD) simulations involving reactive potentials can be used to model material failure. The empirical potentials which are used in such simulations are able to adapt to the atomic environment, at the expense of a significantly higher computational cost than non-reactive potentials. However, during a simulation of failure, the reactive ability is needed only in some limited parts of the system, where bonds break or form and the atomic environment changes. Therefore, simpler non-reactive potentials can be used in the remainder of the system, provided that such potentials reproduce correctly the behavior of the reactive potentials in this region, and that seamless coupling is ensured at the interface between the reactive and non-reactive regions. In this article, we propose a methodology to combine a reactive potential with a non-reactive approximation thereof, made of a set of harmonic pair and angle interactions and whose parameters are adjusted to predict the same energy, geometry and Hessian in the ground state of the potential. We present a methodology to construct the non-reactive approximation of the reactive potential, and a way to couple these two potentials. We also propose a criterion for on-the-fly substitution of the reactive potential by its non-reactive approximation during a simulation. We illustrate the correctness of this hybrid technique for the case of MD simulation of failure in two-dimensional graphene originally modeled with REBO potential.

This paper analyses the GW method for finite electronic systems. In a first step, we provide a mathematical framework for the usual one-body operators that appear naturally in many-body perturbation theory. We then discuss the GW equations which construct an approximation of the one-body Green's function, and give a rigorous mathematical formulation of these equations. Finally, we study the well-posedness of the GW0 equations, proving the existence of a unique solution to these equations in a perturbative regime.

We consider a diffusion equation with highly oscillatory coefficients that admits a homogenized limit. As an alternative to standard corrector problems, we introduce here an embedded corrector problem, written as a diffusion equation in the whole space in which the diffusion matrix is uniform outside some ball of radius $R$. Using that problem, we next introduce three approximations of the homogenized coefficients. These approximations, which are variants of the standard approximations obtained using truncated (supercell) corrector problems, are shown to converge when $R \to \infty$. We also discuss efficient numerical methods to solve the embedded corrector problem.

This article is concerned with the mathematical analysis of the perturbation method for extended Kohn-Sham models, in which fractional occupation numbers are allowed. All our results are established in the framework of the reduced Hartree-Fock (rHF) model, but our approach can be used to study other kinds of extended Kohn-Sham models, under some assumptions on the mathematical structure of the exchange- correlation functional. The classical results of Density Functional Perturbation Theory in the non-degenerate case (that is when the Fermi level is not a degenerate eigenvalue of the mean-field Hamiltonian) are formalized, and a proof of Wigner's (2n + 1) rule is provided. We then focus on the situation when the Fermi level is a degenerate eigenvalue of the rHF Hamiltonian, which had not been considered so far.

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., Comment: 33 pages, 5 figures

In this article, we present a family of numerical approaches to solve high-dimensional linear non-symmetric problems. The principle of these methods is to approximate a function which depends on a large number of variates by a sum of tensor product functions, each term of which is iteratively computed via a greedy algorithm. There exists a good theoretical framework for these methods in the case of (linear and nonlinear) symmetric elliptic problems. However, the convergence results are not valid any more as soon as the problems considered are not symmetric. We present here a review of the main algorithms proposed in the literature to circumvent this difficulty, together with some new approaches. The theoretical convergence results and the practical implementation of these algorithms are discussed. Their behaviors are illustrated through some numerical examples., Comment: 57 pages, 9 figures

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\"odinger 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., Comment: 29 pages, 5 figures

In this article, we set up a functional setting for mean-field electronic structure models of Hartree-Fock or Kohn-Sham types for disordered crystals. The electrons are quantum particles and the nuclei are classical point-like articles whose positions and charges are random. We prove the existence of a minimizer of the energy per unit volume and the uniqueness of the ground state density of such disordered crystals, for the reduced Hartree-Fock model (rHF). We consider both (short-range) Yukawa and (long-range) Coulomb interactions. In the former case, we prove in addition that the rHF ground state density matrix satisfies a self-consistent equation, and that our model for disordered crystals is the thermodynamic limit of the supercell model.

This article deals with the numerical calculation of eigenvalues of perturbed periodic Schr\"odinger 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\'er\'e 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\"odinger operator., Comment: 23 pages, 5 figures

This works extends the recent study on the dielectric permittivity of crystals within the Hartree model [E. Cances and M. Lewin, Arch. Rational Mech. Anal., 197 (2010) 139--177] to the time-dependent setting. In particular, we prove the existence and uniqueness of the nonlinear Hartree dynamics (also called the random phase approximation in the physics literature), in a suitable functional space allowing to describe a local defect embedded in a perfect crystal. We also give a rigorous mathematical definition of the microscopic frequency-dependent polarization matrix, and derive the macroscopic Maxwell-Gauss equation for insulating and semiconducting crystals, from a first order approximation of the nonlinear Hartree model, by means of homogenization arguments.

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., Comment: Proceedings of the Workshop "Numerical Analysis of Multiscale Computations" at Banff International Research Station, December 2009

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\"acker (TFW) theory. The approach follows the same lines as that used in {\it E. Canc\`es, 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., Comment: 34 pages

In this article, we present a greedy algorithm based on a tensor product decomposition, whose aim is to compute the global minimum of a strongly convex energy functional. We prove the convergence of our method provided that the gradient of the energy is Lipschitz on bounded sets. The main interest of this method is that it can be used for high-dimensional nonlinear convex problems. We illustrate this method on a prototypical example for uncertainty propagation on the obstacle problem., Comment: 36 pages, 9 figures, accepted in Mathematical Models and Methods for Applied Sciences

We provide a priori error estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermi-von Weizs\"{a}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., Comment: 50 pages

We provide {\it a priori} error estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermi-von Weizs\"acker (TFW) model and 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 (Part I). 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 in Part II 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 {\it a priori} error estimates for both the spectral and the pseudospectral discretization methods.