Small Basis Set Density Functional Theory Method for Cost-efficient, Large-scale Condensed Matter Simulations

Johannes T. Margraf,Elisabeth Keller,Jack Morgenstein,Karsten Reuter
DOI: https://doi.org/10.26434/chemrxiv-2024-x59gl
2024-06-20
Abstract:We present an efficient first-principles based method geared towards reliably predicting the structures of solid materials across the periodic table. To this end, we use a density functional theory (DFT) baseline with a compact, near-minimal \emph{min+s} basis set, yielding low computational costs and memory demands. Since the use of such small basis set leads to systematic errors in chemical bond lengths, we develop a linear pairwise correction (LPC), available for elements $Z$ = 1-86 (excluding the lanthanide series), parameterized for use with the PBE exchange-correlation functional. We demonstrate the reliability of this corrected approach for equilibrium volumes across the periodic table and the transferability to differently coordinated environments and multi-elemental crystals. We examine relative energies, forces and stresses in geometry optimizations and MD simulations.
Chemistry
What problem does this paper attempt to address?
This paper mainly discusses the problem of how to simulate condensed matter materials on a large scale in a low-cost and efficient manner in materials science. The researchers proposed a first-principles method based on density functional theory (DFT) using a compact approximate minimal basis set (min+s) to reduce computational cost and memory requirements. However, the small basis set can lead to systematic errors in chemical bond lengths. To address this, they developed a linear pairwise correction (LPC) applicable to elements with atomic numbers 1-86 (excluding the lanthanide series) and parameterized it for the PBE exchange-correlation functional. The paper demonstrates the reliability of this method in predicting various solid material structures in the periodic table and proves the transferability of the correction method in different coordination environments and multi-element crystals. The performance of relative energy, forces, and stresses is evaluated through geometric optimization and molecular dynamics simulations. The research reveals that although the small basis set reduces computational cost, it results in underbinding in structural predictions, where lattice constants are overestimated. To correct this issue, they introduce LPC, a linear correction based on element pairs, which can shorten excessively long chemical bonds due to incomplete basis sets. After parameterization of LPC, it effectively reduces errors in equilibrium volumes in the system and improves the accuracy of material structure predictions. In addition, the paper discusses the optimization of computational cost and points out that the cubic complexity of solving the Kohn-Sham equations is the bottleneck when dealing with large simulation units. Despite the proposed iterative solvers and density matrix solvers with linear or subcubic complexity, direct solvers are still preferable for small to medium-sized systems and materials with small band gaps, such as metals and semiconductors. Overall, this paper aims to provide an economical and reliable method for predicting solid material structures by combining small basis set DFT and linear pairwise correction. It is particularly suitable for preliminary screening in materials discovery and generating training data for machine learning potentials.