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

From MaRDI portal
Publication:2446907

DOI10.1016/J.JCP.2012.04.036zbMATH Open1284.65171arXiv1110.1280OpenAlexW1963908375MaRDI QIDQ2446907FDOQ2446907


Authors: P. Motamarri, Mrinal Iyer, V. Gavini, Jaroslaw Knap Edit this on Wikidata


Publication date: 23 April 2014

Published in: Journal of Computational Physics (Search for Journal in Brave)

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.


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




Recommendations




Cites Work


Cited In (16)

Uses Software





This page was built for publication: Higher-order adaptive finite-element methods for orbital-free density functional theory

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