Configurational forces in electronic structure calculations using Kohn-Sham density functional theory

Abstract

We derive the expressions for configurational forces in Kohn-Sham density functional theory, which correspond to the generalized variational force computed as the derivative of the Kohn-Sham energy functional with respect to the position of a material point $\text{x}$. These configurational forces that result from the inner variations of the Kohn-Sham energy functional provide a unified framework to compute atomic forces as well as stress tensor for geometry optimization. Importantly, owing to the variational nature of the formulation, these configurational forces inherently account for the Pulay corrections. The formulation presented in this work treats both pseudopotential and all-electron calculations in single framework, and employs a local variational real-space formulation of Kohn-Sham DFT expressed in terms of the non-orthogonal wavefunctions that is amenable to reduced-order scaling techniques. We demonstrate the accuracy and performance of the proposed configurational force approach on benchmark all-electron and pseudopotential calculations conducted using higher-order finite-element discretization. To this end, we examine the rates of convergence of the finite-element discretization in the computed forces and stresses for various materials systems, and, further, verify the accuracy from finite-differencing the energy. Wherever applicable, we also compare the forces and stresses with those obtained from Kohn-Sham DFT calculations employing plane-wave basis (pseudopotential calculations) and Gaussian basis (all-electron calculations). Finally, we verify the accuracy of the forces on large materials systems involving a metallic aluminum nanocluster containing 666 atoms and an alkane chain containing 902 atoms, where the Kohn-Sham electronic ground state is computed using a reduced-order scaling subspace projection technique (P. Motamarri and V. Gavini, Phys. Rev. B 90, 115127).

32 pages and 20 figures

References in corpus (5)

• Recent progress with large-scale ab initio calculations: the CONQUEST code
• Parallel Self-Consistent-Field Calculations via Chebyshev-Filtered Subspace Acceleration
• 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
• Spectrum-splitting approach for Fermi-operator expansion in all-electron Kohn-Sham DFT calculations

Cited by in corpus (10)

• DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations
• DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization
• Real-time time-dependent density functional theory using higher order finite element methods
• On the calculation of the stress tensor in real-space Kohn-Sham Density Functional Theory
• Fast and robust all-electron density functional theory calculations in solids using orthogonalized enriched finite elements
• Ab initio framework for systems with helical symmetry: theory, numerical implementation and applications to torsional deformations in nanostructures
• Real-space formulation of the stress tensor for $\mathcal{O}(N)$ density functional theory: application to high temperature calculations
• Quasicrystal bulk and surface energies from density functional theory
• Ionic forces and stress tensor in all-electron DFT calculations using enriched finite element basis
• Finite-element methods for noncollinear magnetism and spin-orbit coupling in real-space pseudopotential density functional theory