Higher-order adaptive finite-element methods for orbital-free density functional theory

Abstract

In the present work, we investigate the computational efficiency afforded by higher-order finite-element discretization of the saddle-point formulation of orbital-free density functional theory. We first investigate the robustness of viable solution schemes by analyzing the solvability conditions of the discrete problem. We find that a staggered solution procedure where the potential fields are computed consistently for every trial electron-density is a robust solution procedure for higher-order finite-element discretizations. We next study the numerical convergence rates for various orders of finite-element approximations on benchmark problems. We obtain close to optimal convergence rates in our studies, although orbital-free density-functional theory is nonlinear in nature and some benchmark problems have Coulomb singular potential fields. We finally investigate the computational efficiency of various higher-order finite-element discretizations by measuring the CPU time for the solution of discrete equations on benchmark problems that include large Aluminum clusters. In these studies, we use mesh coarse-graining rates that are derived from error estimates and an a priori knowledge of the asymptotic solution of the far-field electronic fields. Our studies reveal a significant 100-1000 fold computational savings afforded by the use of higher-order finite-element discretization, alongside providing the desired chemical accuracy. We consider this study as a step towards developing a robust and computationally efficient discretization of electronic structure calculations using the finite-element basis.

43 Pages

References in corpus (2)

• Adaptive local basis set for Kohn-Sham density functional theory in a discontinuous Galerkin framework I: Total energy calculation
• All-electron density functional theory and time-dependent density functional theory with high-order finite elements

Cited by in corpus (15)

• Higher-order adaptive finite-element methods for Kohn-Sham density functional theory
• SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Extended systems
• SPARC: Accurate and efficient finite-difference formulation and parallel implementation of Density Functional Theory: Isolated clusters
• Large-scale all-electron density functional theory calculations using an enriched finite element basis
• Augmented Lagrangian formulation of Orbital-Free Density Functional Theory
• A subquadratic-scaling subspace projection method for large-scale Kohn-Sham density functional theory calculations using spectral finite-element discretization
• Electronic-structure study of an edge dislocation in Aluminum and the role of macroscopic deformations on its energetics
• Configurational forces in electronic structure calculations using Kohn-Sham density functional theory
• Highly Accurate Local Pseudopotentials of Li, Na, and Mg for Orbital Free Density Functional Theory
• Electronic structure study of screw dislocation core energetics in Aluminum and core energetics informed forces in a dislocation aggregate
• Real-time time-dependent density functional theory using higher order finite element methods
• Orbital-free DFT study of the energetics of vacancy clustering and prismatic dislocation loop nucleation in aluminum
• Ab initio framework for systems with helical symmetry: theory, numerical implementation and applications to torsional deformations in nanostructures
• Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements
• Density functional theory method for twisted geometries with application to torsional deformations in group-IV nanotubes