Communications in Mathematical Sciences

Volume 5 (2007)

Number 4

A sub-linear scaling algorithm for computing the electronic structure of materials

Pages: 999 – 1026

DOI: https://dx.doi.org/10.4310/CMS.2007.v5.n4.a14

Authors

Weinan E

Carlos J. García-Cervera

Jianfeng Lu

Abstract

We introduce a class of sub-linear scaling algorithms for analyzing the electronic structure of crystalline solids with isolated defects. We divide the localized orbitals of the electrons into two sets: one set associated with the atoms in the region where the deformation of the material is smooth (smooth region), and the other set associated with the atoms around the defects (non- smooth region). The orbitals associated with atoms in the smooth region can be approximated accurately using asymptotic analysis. The results can then be used in the original formulation to find the orbitals in the non-smooth region. For orbital-free density functional theory, one can simply partition the electron density into a sum of the density in the smooth region and a density in the non-smooth region. This kind of partition is not used for Kohn-Sham density functional theory and one uses instead the partition of the set of orbitals. As a byproduct, we develop the necessary real space formulations and we present a formulation of the electronic structure problem for a subsystem, when the electronic structure for the remaining part is known.

Keywords

sub-linear scaling algorithms, asymptotics, DFT-continuum approximation, density functional theory

2010 Mathematics Subject Classification

34E05, 35Q40, 74Q05

Published 1 January 2007