Abstract:The density functional theory has become a popular tool for electronic structure calculations in various fields include chemistry, material science and biology in recent years. The thesis is devoted to the mathematical analysis and algorithm developments for density functional theory. The dissertation consists of three parts. In the first part, we introduce the density functional theory, in particular the Kohn-Sham density functional theory. We begin with the mathematical formulation and structure of density functional theory. By exploiting the invariance properties of the model, we reformulate density functional theory as a variational problem on subspaces, the minimizer is the occupied space in quantum mechanics terminology. We will then focus on the difference between metallic and insulating systems, in particular, the locality principle for insulators: The exponential decay of the Green's function, the density matrix and the Wannier functions. The second part is some mathematical results related with the Kohn-Sham density functional theory. To gain better understanding of the Kohn-Sham map arising from the theory, we analyze the linearized Kohn-Sham map in detail, in particular its action on different function spaces. The rest of the second part is devoted to the study of the electronic structure of elastically deformed crystals. For a reduced model, we validate the Cauchy-Born rule in the continuum limit under stability conditions. The notion of Wannier functions is generalized to elastically deformed systems, and we also establish the Cauchy-Born rule for Wannier functions, and hence for the electronic structure. We discuss about numerical issues related with density functional theory in the third part. We first discuss the self-consistent iteration and show different behaviors in convergence for insulating and metallic systems. Then we focus on the developments of more efficient algorithms for solving the model: Linear scaling algorithms for insulating systems based on the locality principle, and also efficient algorithms for metallic systems based on better representation of the Fermi operator and fast algorithm to extract the diagonal of an inverse matrix. We then discuss the sublinear scaling algorithms based on the knowledge of elastically deformed system. This is an example of multiscale methods, and as such, we present a thorough discussion about the treatment of boundary between the fine-scale and the coarse-scale regions and how to correctly formulate the problem to avoid artifacts on the boundary.

Linear Scaling DFT for defects in metals

Linear scaling DFT calculations for large Tungsten systems using an optimized local basis

A Sub-Linear Scaling Algorithm for Computing the Electronic Structure of Materials

Large scale ab-initio simulations of dislocations

Spectral quadrature for the first principles study of crystal defects: Application to magnesium

Linear scaling Krylov subspace method for large scale {\it ab initio} electronic structure calculations of metals

Unveiling Defect Physics in Gapped Metals: A Theoretical Investigation into Defect Formation and Electronic Structure Interplay

Evolutionary computing and machine learning for the discovering of low-energy defect configurations

Importance of elastic finite-size effects: neutral defects in ionic compounds

Energetics of point defects in aluminum via orbital-free density functional theory

The decoupled DFT-$\frac{1}{2}$ method for defect excitation energies

Defect Energy Levels in Density Functional Calculations: Alignment and Band Gap Problem

Aperiodic defects in periodic solids

Excited state properties of point defects in semiconductors and insulators investigated with time-dependent density functional theory

Analysis of nearly planar defects using the Thomas--Fermi--von Weiszacker model

Optical Properties of Defects in Solids via Quantum Embedding with Good Active Space Orbitals

Defect energetics of cubic hafnia from quantum Monte Carlo simulations

Multiscale Analysis of Density Functional Theory

Optimizing density-functional simulations for two-dimensional metals

Finite element approximation of the fields of bulk and interfacial line defects

Multiscale Quantum-Defect Theory and Its Application to Atomic Spectrum