Chemical bonding in large systems using projected population analysis from real-space density functional theory calculations

Abstract

We present an efficient and scalable computational approach for conducting projected population analysis from real-space finite-element (FE) based Kohn-Sham density functional theory calculations (DFT-FE). This work provides an important direction towards extracting chemical bonding information from large-scale DFT calculations on materials systems involving thousands of atoms while accommodating periodic, semi-periodic or fully non-periodic boundary conditions. Towards this, we derive the relevant mathematical expressions and develop efficient numerical implementation procedures that are scalable on multi-node CPU architectures to compute the projected overlap and Hamilton populations. The population analysis is accomplished by projecting either the self-consistently converged FE discretized Kohn-Sham orbitals, or the FE discretized Hamiltonian onto a subspace spanned by a localized atom-centred basis set. The proposed methods are implemented in a unified framework within DFT-FE code where the ground-state DFT calculations and the population analysis are performed on the same FE grid. We further benchmark the accuracy and performance of this approach on representative material systems involving periodic and non-periodic DFT calculations with LOBSTER, a widely used projected population analysis code. Finally, we discuss a case study demonstrating the advantages of our scalable approach to extract the quantitative chemical bonding information of hydrogen chemisorbed in large silicon nanoparticles alloyed with carbon, a candidate material for hydrogen storage.

arXiv comment24 Figures, 6 Tables, 57 pages with references and supplementary information

References in corpus (11)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Optimized norm-conserving Vanderbilt pseudopotentials
• Maximally localized Wannier functions: Theory and applications
• The PseudoDojo: Training and grading a 85 element optimized norm-conserving pseudopotential table
• Projection of plane-wave calculations into atomic orbitals
• DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations
• 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: Isolated clusters
• DFT-FE 1.0: A massively parallel hybrid CPU-GPU density functional theory code using finite-element discretization
• Theory and applications of generalized Pipek--Mezey Wannier functions
• Population Analysis with Wannier Orbitals