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.

Direct ab initio calculation of the $^{4}$He nuclear electric dipole polarizability

Electric dipole polarizability from first principles calculations

Low-Energy Dipole Strength in He

Electric dipole polarizability of low-lying excited states in atomic nuclei

Low-energy dipole strength in 8He

Ab initio many-body calculations of the 3H(d,n)4He and 3He(d,p)4He fusion

Ab Initio Nuclear Many-Body Perturbation Calculations in the Hartree-Fock Basis

Approaching exact hyperpolarizabilities via sum-over-states Monte Carlo configuration interaction

Calculating the many-potential vacuum polarization density of the Dirac equation in the finite-basis approximation

Compton Scattering on 4He with Nuclear One- and Two-Body Densities

Magnetic polarizability of a charged pion from four-point functions in lattice QCD

Global trends of the electric dipole polarizability from shell-model calculations

Ab initio many-body calculations of nucleon-4He scattering with three-nucleon forces

QED of Bernabéu-Tarrach sum rule for electric polarizability and its implication for the Lamb shift

High-precision nonadiabatic calculations of dynamic polarizabilities and hyperpolarizabilities for the lowlying vibrational-rotational states of hydrogen molecular ions

Nuclear polarizability effects in $^3$He$^+$ hyperfine splitting

Accurate Calculation of Multipole Polarizabilities for One-Electron Atom in Debye Plasmas

Robust ab initio predictions for dimensionless ratios of E2 and radius observables. I. Electric quadrupole moments and deformation

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

High-precision Calculation of the Geometric Quantities of Two-Electron Atoms Based on the Hylleraas Configuration-Interaction Basis