PEXSI-: a Green's function embedding method for Kohn-Sham density functional theory

DOI10.4310/AMSA.2018.V3.N2.A3MaRDI QIDQ1797855zbMATH OpenFDO

Authors Xiantao Li, Lin Lin, Jianfeng Lu

Publication date 22 October 2018

Published in Annals of Mathematical Sciences and Applications (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1606.00515

zbMATH Keywords

grapheneDyson equationKohn-Sham density functional theoryGreen's function embedding method

Mathematics Subject Classification ID

Fundamental solutions, Green's function methods, etc. for boundary value problems involving PDEs (65N80) NLS equations (nonlinear Schrödinger equations) (35Q55) Quantum equilibrium statistical mechanics (general) (82B10) Statistical mechanics of nanostructures and nanoparticles (82D80)

Abstract: In this paper, we propose a new Green's function embedding method called PEXSI-$Sigma$ for describing complex systems within the Kohn-Sham density functional theory (KSDFT) framework, after revisiting the physics literature of Green's function embedding methods from a numerical linear algebra perspective. The PEXSI-$Sigma$ method approximates the density matrix using a set of nearly optimally chosen Green's functions evaluated at complex frequencies. For each Green's function, the complex boundary conditions are described by a self energy matrix $Sigma$ constructed from a physical reference Green's function, which can be computed relatively easily. In the linear regime, such treatment of the boundary condition can be numerically exact. The support of the $Sigma$ matrix is restricted to degrees of freedom near the boundary of computational domain, and can be interpreted as a frequency dependent surface potential. This makes it possible to perform KSDFT calculations with $mathcalO(N^2)$ computational complexity, where $N$ is the number of atoms within the computational domain. Green's function embedding methods are also naturally compatible with atomistic Green's function methods for relaxing the atomic configuration outside the computational domain. As a proof of concept, we demonstrate the accuracy of the PEXSI-$Sigma$ method for graphene with divacancy and dislocation dipole type of defects using the DFTB+ software package.

Recommendations

• Projection-based embedding theory for solving Kohn-Sham density functional theory
• Non-periodic finite-element formulation of Kohn-Sham density functional theory
• Density functional theory calculations using the finite element method
• Coarse-graining Kohn-Sham density functional theory
• Quasi-continuum orbital-free density-functional theory: a route to multi-million atom non-periodic DFT calculation

Cited in

(4)

• Thermodynamic limit of crystal defects with finite temperature tight binding
• Projection-based embedding theory for solving Kohn-Sham density functional theory
• QM/MM Methods for Crystalline Defects. Part 3: Machine-Learned MM Models
• Computing edge states without hard truncation

Describes a project that uses

Uses Software

• Quickstep
• SelInv
• DFTB+

This page was built for publication: PEXSI-\(\Sigma\): a Green's function embedding method for Kohn-Sham density functional theory

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1797855)