-
Spherical-harmonic-analysis-based optimization of atomic weighting functions for multicenter numerical integration in molecules
Authors:
Dimitri N. Laikov
Abstract:
The well-known spatial integration schemes in molecular electronic structure theory, immune to cusps and point singularities of some kind at atomic positions, use a set of weighting functions to split the integrand into a sum of atom-centered parts, each dealt with in its own spherical coordinate system. Here, for a given set of integrands in the two-center case, a quality measure of the weighting…
▽ More
The well-known spatial integration schemes in molecular electronic structure theory, immune to cusps and point singularities of some kind at atomic positions, use a set of weighting functions to split the integrand into a sum of atom-centered parts, each dealt with in its own spherical coordinate system. Here, for a given set of integrands in the two-center case, a quality measure of the weighting functions is defined to compare, design, and optimize them, it is roughly proportional to the average number of angular quadrature points needed to reach a given integration accuracy. A study of Becke's fuzzy Voronoï cells has helped to improve their performance by a new modification. New spherically-symmetric unnormalized weighting functions are found in the form of a negative power times the negative exponential of the fourth power of the scaled distance to the atomic center, with the parameters related to the asymptotic decay of the integrand and the integration accuracy -- these are much simpler but no less efficient and naturally fit for linear-scaling calculations. Radial distribution of spherical quadrature orders is studied. A radial integration scheme of double exponential type is optimized. A symmetric analog of the pseudospectral approximation is used for the seminumerical evaluation of two-electron repulsion integrals. Taken together, this allows efficient calculation of all molecular integrals with well-controlled accuracy, as shown by tests on a set of molecules.
△ Less
Submitted 13 September, 2022; v1 submitted 10 June, 2020;
originally announced June 2020.
-
Optimization of atomic density-fitting basis functions for molecular two-electron integral approximations
Authors:
Dimitri N. Laikov
Abstract:
A general procedure for the optimization of atomic density-fitting basis functions is designed with the balance between accuracy and numerical stability in mind. Given one-electron wavefunctions and energies, weights are assigned to the product densities, modeling their contribution to the exchange and second-order correlation energy, and a simple weighted error measure is minimized. Generally-con…
▽ More
A general procedure for the optimization of atomic density-fitting basis functions is designed with the balance between accuracy and numerical stability in mind. Given one-electron wavefunctions and energies, weights are assigned to the product densities, modeling their contribution to the exchange and second-order correlation energy, and a simple weighted error measure is minimized. Generally-contracted Gaussian auxiliary basis sets are optimized to match the wavefunction basis sets [D. N. Laikov, Theor. Chem. Acc. 138, 40 (2019)] for all 102 elements in a scalar-relativistic approximation [D. N. Laikov, J. Chem. Phys. 150, 061102 (2019)].
△ Less
Submitted 8 September, 2020; v1 submitted 20 May, 2020;
originally announced May 2020.
-
Simple exchange hole models for long-range-corrected density functionals
Authors:
Dimitri N. Laikov
Abstract:
Density functionals with a range-separated treatment of the exchange energy are known to improve upon their semilocal forerunners and fixed-fraction hybrids. The conversion of a given semilocal functional into its short-range analog is not straightforward, however, and not even unique, because the latter has a higher information content that has to be recovered in some way. Simple models of the sp…
▽ More
Density functionals with a range-separated treatment of the exchange energy are known to improve upon their semilocal forerunners and fixed-fraction hybrids. The conversion of a given semilocal functional into its short-range analog is not straightforward, however, and not even unique, because the latter has a higher information content that has to be recovered in some way. Simple models of the spherically-averaged exchange hole as an interpolation between the uniform electron gas limit and a few-term Hermite function are developed here for use with generalized-gradient approximations, so that the energy density of the error-function-weighted Coulomb interaction is given by explicit closed-form expressions in terms of elementary and error functions. For comparison, some new non-oscillatory models in the spirit of earlier works are also built and studied, their energy densities match rather closely (within less than 5%) but do lack the exact uniform electron gas limit.
△ Less
Submitted 20 May, 2019;
originally announced May 2019.
-
Atomic effective potentials for starting molecular electronic structure calculations
Authors:
Dimitri N. Laikov,
Ksenia R. Briling
Abstract:
Atomic effective one-electron potentials in a compact analytic form in terms of a few Gaussian charge distributions are developed, for Hydrogen through Nobelium, for starting molecular electronic structure calculations by a simple diagonalization. For each element, all terms but one are optimized in an isolated-atom Hartree--Fock calculation, and the last one is parametrized on a set of molecules.…
▽ More
Atomic effective one-electron potentials in a compact analytic form in terms of a few Gaussian charge distributions are developed, for Hydrogen through Nobelium, for starting molecular electronic structure calculations by a simple diagonalization. For each element, all terms but one are optimized in an isolated-atom Hartree--Fock calculation, and the last one is parametrized on a set of molecules. This one-parameter-per-atom model gives a good starting guess for typical molecules and may be of interest even on its own.
△ Less
Submitted 6 January, 2020; v1 submitted 8 February, 2019;
originally announced February 2019.
-
Additive atomic approximation for relativistic effects: a two-component Hamiltonian for molecular electronic structure calculations
Authors:
Dimitri N. Laikov
Abstract:
An approximate relativistic two-component Hamiltonian for use in molecular electronic structure calculations is derived in the form of a sum of fixed atom-centered kinetic and spin-orbit operators added to the non-relativistic Hamiltonian. Starting from the well-known zeroth-order regular approximation, further steps are taken to get rid of its nonlinearity in the potential, ending up with a simpl…
▽ More
An approximate relativistic two-component Hamiltonian for use in molecular electronic structure calculations is derived in the form of a sum of fixed atom-centered kinetic and spin-orbit operators added to the non-relativistic Hamiltonian. Starting from the well-known zeroth-order regular approximation, further steps are taken to get rid of its nonlinearity in the potential, ending up with a simple formulation with easily computable integrals that can seamlessly work with any traditional electronic structure method. Molecular tests show a good accuracy of this approximation.
△ Less
Submitted 20 November, 2018;
originally announced November 2018.
-
Atomic basis functions for molecular electronic structure calculations
Authors:
Dimitri N. Laikov
Abstract:
Electronic structure methods for accurate calculation of molecular properties have a high cost that grows steeply with the problem size, therefore, it is helpful to have the underlying atomic basis functions that are less in number but of higher quality. Following our earlier work [Chem. Phys. Lett. 416, 116 (2005)] where general correlation-consistent basis sets are defined, for any atom, as solu…
▽ More
Electronic structure methods for accurate calculation of molecular properties have a high cost that grows steeply with the problem size, therefore, it is helpful to have the underlying atomic basis functions that are less in number but of higher quality. Following our earlier work [Chem. Phys. Lett. 416, 116 (2005)] where general correlation-consistent basis sets are defined, for any atom, as solutions of purely atomic functional minimization problems, and which are shown to work well for chemical bonding in molecules, we take a further step here and define a new kind of atomic polarization functionals, the minimization of which yields additional sets of diffuse functions that help to calculate better molecular electron affinities, polarizabilities, and intermolecular dispersion interactions. Analytical representations by generally-contracted Gaussian functions of up to microhartree numerical accuracy grades are developed for atoms Hydrogen through Nobelium within the four-component Dirac-Coulomb theory and its scalar-relativistic approximation, and also for Hydrogen through Krypton in the two-component nonrelativistic case. The convergence of correlation energy with the basis set size is studied, and complete-basis-set extrapolation formulas are developed.
△ Less
Submitted 16 August, 2018;
originally announced August 2018.
-
A new parametrizable model of molecular electronic structure
Authors:
Dimitri N. Laikov
Abstract:
A new electronic structure model is developed in which the ground state energy of a molecular system is given by a Hartree-Fock-like expression with parametrized one- and two-electron integrals over an extended (minimal + polarization) set of orthogonalized atom-centered basis functions, the variational equations being solved formally within the minimal basis but the effect of polarization functio…
▽ More
A new electronic structure model is developed in which the ground state energy of a molecular system is given by a Hartree-Fock-like expression with parametrized one- and two-electron integrals over an extended (minimal + polarization) set of orthogonalized atom-centered basis functions, the variational equations being solved formally within the minimal basis but the effect of polarization functions being included in the spirit of second-order perturbation theory. It is designed to yield good dipole polarizabilities and improved intermolecular potentials with dispersion terms. The molecular integrals include up to three-center one-electron and two-center two-electron terms, all in simple analytical forms. A method to extract the effective one-electron Hamiltonian of nonlocal-exchange Kohn-Sham theory from the coupled-cluster one-electron density matrix is designed and used to get its matrix representation in a molecule-intrinsic minimal basis as an input to the paramtrization procedure -- making a direct link to the correlated wavefunction theory. The model has been trained for 15 elements (H, Li--F, Na--Cl, 720 parameters) on a set of 5581 molecules (including ions, transition states, and weakly-bound complexes) whose first- and second-order properties were computed by the coupled-cluster theory as a reference, and a good agreement is seen. The model looks promising for the study of large molecular systems, it is believed to be an important step forward from the traditional semiempirical models towards higher accuracy at nearly as low a computational cost.
△ Less
Submitted 29 August, 2011; v1 submitted 23 June, 2011;
originally announced June 2011.