Abstract:We present the implementation of a variational finite element solver in the HelFEM program for benchmark calculations on diatomic systems. A basis set of the form $\chi_{nlm}(\mu,\nu,\phi)=B_{n}(\mu)Y_{l}^{m}(\nu,\phi)$ is used, where $(\mu,\nu,\phi)$ are transformed prolate spheroidal coordinates, $B_{n}(\mu)$ are finite element shape functions, and $Y_{l}^{m}$ are spherical harmonics. The basis set allows for an arbitrary level of accuracy in calculations on diatomic molecules, which can be performed at present with either nonrelativistic Hartree--Fock (HF) or density functional (DF) theory. Hundreds of DFs at the local spin-density approximation (LDA), generalized gradient approximation (GGA) and the meta-GGA level can be used through an interface with the Libxc library; meta-GGA and hybrid DFs aren't available in other fully numerical diatomic program packages. Finite electric fields are also supported in HelFEM, enabling access to electric properties. We introduce a powerful tool for adaptively choosing the basis set by using the core Hamiltonian as a proxy for its completeness. HelFEM and the novel basis set procedure are demonstrated by reproducing the restricted open-shell HF limit energies of 68 diatomic molecules from the $1^{\text{st}}$ to the $4^{\text{th}}$ period with excellent agreement with literature values, despite requiring \emph{orders of magnitude} fewer parameters for the wave function. Then, the electric properties of the BH and N2 molecules under finite field are studied, again yielding excellent agreement with previous HF limit values for energies, dipole moments, and dipole polarizabilities, again with much more compact wave functions than what were needed in the literature references. Finally, HF, LDA, GGA, and meta-GGA calculations of the atomization energy of N2 are performed, demonstrating the superb accuracy of the present approach.

Accurate Hartree-Fock energy of extended systems using large Gaussian basis sets

Fully numerical Hartree-Fock and density functional calculations. II. Diatomic molecules

Accurate Hellmann-Feynman forces from density functional calculations with augmented Gaussian basis sets

Exact Exchange with Range-separated Algorithm for Thermodynamic Limit of Periodic Hartree-Fock Theory

Adaptive Local Basis Set for Kohn-Sham Density Functional Theory in a Discontinuous Galerkin Framework II: Force, Vibration, and Molecular Dynamics Calculations

Efficient Hartree-Fock Exchange Algorithm with Coulomb Range Separation and Long-Range Density Fitting

Accurate Localized Resolution Of Identity Approach For Linear-Scaling Hybrid Density Functionals And For Many-Body Perturbation Theory

Efficient Calculation of the Lattice Thermal Conductivity by Atomistic Simulations with Ab Initio Accuracy

Efficient periodic resolution-of-the-identity Hartree–Fock exchange method with k -point sampling and Gaussian basis sets

Resolution-of-identity approach to Hartree–Fock, hybrid density functionals, RPA, MP2 and GW with numeric atom-centered orbital basis functions

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

Extended Hartree-Fock method based on pair density functional theory

Accurate Many-Body Electronic Structure Near the Basis Set Limit: Application to the Chromium Dimer

Energy Density Functionals From the Strong-Coupling Limit Applied to the Anions of the He Isoelectronic Series

Geometrical correction for the inter- and intramolecular basis set superposition error in periodic density functional theory calculations

Accurate Prediction of Heat of Formation by Combining Hartree-Fock/Density Functional Theory Calculation with Linear Regression Correction Approach

Building a New Platform for Significantly Improving Performance of Hartree-Fock and CCSD(T) Correlation Energy Based on Two-Point Complete Basis Set Extrapolation Schemes

Density fitting in periodic systems: application to TDHF in diamond and oxides

Ab Initio Calculation of Coupling-Constant Averaged Exchange-Correlation Holes for Spherically Symmetric Atoms

Fully selfconsistent GW calculations for molecules