Treecode-accelerated Green Iteration for Kohn-Sham Density Functional Theory
arXiv:2003.01833 · doi:10.1016/j.jcp.2020.110101
Abstract
We present a real-space computational method called treecode-accelerated Green Iteration (TAGI) for all-electron Kohn-Sham Density Functional Theory. TAGI is based on a reformulation of the Kohn-Sham equations in which the eigenvalue problem in differential form is converted into a fixed-point problem in integral form by convolution with the modified Helmholtz Green's function. In each self-consistent field (SCF) iteration, the fixed-points are computed by Green Iteration, where the discrete convolution sums are efficiently evaluated by a GPU-accelerated barycentric Lagrange treecode. Other techniques used in TAGI include adaptive mesh refinement, Fejér quadrature, singularity subtraction, gradient-free eigenvalue update, and Anderson mixing to accelerate convergence of the SCF and Green Iterations. Ground state energy computations of several atoms (Li, Be, O) and small molecules (H, CO, CH) demonstrate TAGI's ability to efficiently achieve chemical accuracy.
References in corpus (14)
- Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
- Advanced capabilities for materials modelling with Quantum ESPRESSO
- Optimized norm-conserving Vanderbilt pseudopotentials
- Perspective on density functional theory
- Libxc: a library of exchange and correlation functionals for density functional theory
- DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations
- Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration
- Higher-order adaptive finite-element methods for Kohn-Sham density functional theory
- Real-space local polynomial basis for solid-state electronic-structure calculations: A finite-element approach
- All-electron density functional theory and time-dependent density functional theory with high-order finite elements
- Large-scale all-electron density functional theory calculations using an enriched finite element basis
- A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization
- A kernel-independent treecode based on barycentric Lagrange interpolation
- Grid-based electronic structure calculations: the tensor decomposition approach