Accurate approximations of density functional theory for large systems with applications to defects in crystalline solids

This chapter presents controlled approximations of Kohn-Sham density functional theory (DFT) that enable very large scale simulations. The work is motivated by the study of defects in crystalline solids, though the ideas can be used in other applications. The key idea is to formulate DFT as a minimization problem over the density operator, and to cast spatial and spectral discretization as systematically convergent approximations. This enables efficient and adaptive algorithms that solve the equations of DFT with no additional modeling, and up to desired accuracy, for very large systems, with linear and sublinear scaling. Various approaches based on such approximations are presented, and their numerical performance demonstrated through selected examples. These examples also provide important insight about the mechanics and physics of defects in crystalline solids.