Review of the finite difference Hartree-Fock method for atoms and diatomic molecules, and its implementation in the x2dhf program

Abstract

We present an extensive review of the two-dimensional finite difference Hartree--Fock (FD HF) method, and present its implementation in the newest version of X2DHF, the FD HF program for atoms and diatomic molecules. The program was originally published in Comput. Phys. Commun. in 1996, and was last revised in 2013. X2DHF can be used to obtain HF limit values of total energies and multipole moments for a wide range of diatomic molecules and their ions, using either point nuclei or a finite nuclear model. Polarizabilities (${\alpha }_{zz}$) and hyperpolarizabilities (${\beta }_{zzz}$, ${\gamma }_{zzzz}$, $A_{z,zz}$, $B_{zz,zz}$) can also be computed by the program with the finite-field method. X2DHF has been extensively used in the literature to assess the accuracy of existing atomic basis sets and to help in developing new ones. As a new feature since the last revision, the program can now also perform Kohn-Sham density functional calculations with local and generalized gradient exchange-correlation functionals with the Libxc library of density functionals, enabling new types of studies. Furthermore, the initialization of calculations has been greatly simplified. As before, X2DHF can also perform one-particle calculations with (smooth) Coulomb, Green-Sellin-Zachor and Krammers-Henneberger potentials, while calculations with a superposition of atomic potentials have been added as a new feature. The program is easy to install from the GitHub repository and build via CMake using the x2dhfctl script that facilitates creating its single- and multiple-threaded versions, as well as building in Libxc support. Calculations can be carried out with X2DHF in double- or quadruple-precision arithmetic.

119 pages, 16 figures

References in corpus (19)

• Piecewise Linearity of Approximate Density Functionals Revisited: Implications for Frontier Orbital Energies
• High-harmonic generation from arbitrarily oriented diatomic molecules including nuclear motion and field-free alignment
• An overview of self-consistent field calculations within finite basis sets
• A self-interaction-free local hybrid functional: Accurate binding energies vis-à-vis accurate ionization potentials from Kohn-Sham eigenvalues
• Many recent density functionals are numerically ill-behaved
• Fully numerical calculations on atoms with fractional occupations. Range-separated exchange functionals
• Meta-GGA density functional calculations on atoms with spherically symmetric densities in the finite element formalism
• A call to arms: making the case for more reusable libraries
• Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets
• Atomic electronic structure calculations with Hermite interpolating polynomials
• Polarized Gaussian basis sets from one-electron ions
• On the accurate reproduction of strongly repulsive interatomic potentials
• Reproducibility of density functional approximations: how new functionals should be reported
• Short range repulsive interatomic interactions in energetic processes in solids
• Insight on Gaussian basis set truncation errors in weak to intermediate magnetic fields with an approximate Hamiltonian
• Efficient Algorithm for Two-Center Coulomb and Exchange Integrals of Electronic Prolate Spheroidal Orbitals
• On the accuracy of a recent regularized nuclear potential
• Nuclear cusps and singularities in the non-additive kinetic potential bi-functional from analytical inversion
• Importance profiles. Visualization of atomic basis set requirements