Atomic Electronic Structure Calculations with Hermite Interpolating Polynomials

We have recently described the implementation of atomic electronic structure calculations within the finite element method with numerical radial basis functions of the form χμ(r) = r–1Bμ(r), where high-order Lagrange interpolating polynomials (LIPs) were used as the shape functions Bμ(r). In this work, we discuss how χμ(r) can be evaluated in a stable manner at small r and also revisit the choice of the shape functions Bμ(r). Three kinds of shape functions are considered: in addition to the continuous LIPs, we consider the analytical implementation of first-order Hermite interpolating polynomials (HIPs) that are continuous, as well as numerical implementations of n-th order ( continuous) HIPs that are expressed in terms of an underlying high-order LIP basis. Furnished with the new implementation, we demonstrate that the first-order HIPs are reliable even with large numbers of nodes and that they also work with nonuniform element grids, affording even better results in atomic electronic structure calculations than LIPs with the same total number of basis functions. We demonstrate that discontinuities can be observed in the spin-σ local kinetic energy τσ in small LIP basis sets, while HIP basis sets do not suffer from such issues; however, either set can be used to reach the complete basis set limit with smooth τσ. Moreover, we discuss the implications of HIPs on calculations with meta-GGA functionals with a number of recent meta-GGA functionals, and we find most Minnesota functionals to be ill-behaved. We also examine the potential usefulness of the explicit control over the derivative in HIPs for forming numerical atomic orbital basis sets, but we find that confining potentials are still likely a better option.

[1]  S. Lehtola.  Meta-GGA Density Functional Calculations on Atoms with Spherically Symmetric Densities in the Finite Element Formalism , 2023, Journal of chemical theory and computation.

[2]  F. Gygi.  All-Electron Plane-Wave Electronic Structure Calculations , 2023, Journal of chemical theory and computation.

[3]  S. Lehtola,et al.  How good are recent density functionals for ground and excited states of one-electron systems? , 2022, The Journal of chemical physics.

[4]  S. Lehtola,et al.  Many recent density functionals are numerically ill-behaved. , 2022, The Journal of chemical physics.

[5]  S. Kümmel,et al.  First steps towards achieving both ultranonlocality and a reliable description of electronic binding in a meta-generalized gradient approximation , 2022, Physical Review Research.

[6]  M. Côté,et al.  Cubic spline solver for generalized density functional treatments of atoms and generation of atomic datasets for use with exchange-correlation functionals including meta-GGA , 2022, Physical Review B.

[7]  Susi Lehtola,et al.  Free and open source software for computational chemistry education , 2021, WIREs Computational Molecular Science.

[8]  S. Lehtola,et al.  Meta-Local Density Functionals: A New Rung on Jacob’s Ladder , 2020, Journal of chemical theory and computation.

[9]  J. Perdew,et al.  Correction to "Accurate and Numerically Efficient r2SCAN Meta-Generalized Gradient Approximation". , 2020, The journal of physical chemistry letters.

[10]  J. Perdew,et al.  Accurate and Numerically Efficient r 2 SCAN Meta-Generalized Gradient Approximation , 2020 .

[11]  J. Perdew,et al.  Accurate and numerically efficient r$^2$SCAN meta-generalized gradient approximation , 2020, 2008.03374.

[12]  E. Engel,et al.  Efficient implementation of the superposition of atomic potentials initial guess for electronic structure calculations in Gaussian basis sets. , 2020, The Journal of chemical physics.

[13]  D. Truhlar,et al.  M06-SX screened-exchange density functional for chemistry and solid-state physics , 2020, Proceedings of the National Academy of Sciences of the United States of America.

[14]  S. Lehtola.  Polarized Gaussian basis sets from one-electron ions. , 2020, The Journal of chemical physics.

[15]  S. Lehtola.  Accurate reproduction of strongly repulsive interatomic potentials , 2019, Physical Review A.

[16]  S. Lehtola,et al.  An Overview of Self-Consistent Field Calculations Within Finite Basis Sets † , 2019, Molecules.

[17]  S. Lehtola.  Fully numerical calculations on atoms with fractional occupations and range-separated exchange functionals , 2019, Physical Review A.

[18]  S. Kümmel,et al.  Ultranonlocality and accurate band gaps from a meta-generalized gradient approximation , 2019, Physical Review Research.

[19]  S. Lehtola.  A review on non‐relativistic, fully numerical electronic structure calculations on atoms and diatomic molecules , 2019, International Journal of Quantum Chemistry.

[20]  Xiao He,et al.  Revised M11 Exchange-Correlation Functional for Electronic Excitation Energies and Ground-State Properties. , 2019, The journal of physical chemistry. A.

[21]  S. Lehtola,et al.  Fully numerical electronic structure calculations on diatomic molecules in weak to strong magnetic fields , 2018, Molecular Physics.

[22]  Susi Lehtola,et al.  Fully numerical Hartree‐Fock and density functional calculations. II. Diatomic molecules , 2018, International Journal of Quantum Chemistry.

[23]  S. Lehtola.  Fully numerical Hartree‐Fock and density functional calculations. I. Atoms , 2018, International Journal of Quantum Chemistry.

[24]  Susi Lehtola,et al.  Assessment of Initial Guesses for Self-Consistent Field Calculations. Superposition of Atomic Potentials: Simple yet Efficient , 2018, Journal of chemical theory and computation.

[25]  L. Goerigk,et al.  The Nonlocal Kernel in van der Waals Density Functionals as an Additive Correction: An Extensive Analysis with Special Emphasis on the B97M-V and ωB97M-V Approaches. , 2018, Journal of chemical theory and computation.

[26]  D. Truhlar,et al.  Revised M06 density functional for main-group and transition-metal chemistry , 2018, Proceedings of the National Academy of Sciences.

[27]  Micael J. T. Oliveira,et al.  Recent developments in libxc - A comprehensive library of functionals for density functional theory , 2018, SoftwareX.

[28]  Haoyu S. Yu,et al.  Revised M06-L functional for improved accuracy on chemical reaction barrier heights, noncovalent interactions, and solid-state physics , 2017, Proceedings of the National Academy of Sciences.

[29]  Xiao He,et al.  Correction: MN15: A Kohn–Sham global-hybrid exchange–correlation density functional with broad accuracy for multi-reference and single-reference systems and noncovalent interactions , 2016, Chemical science.

[30]  M. Head‐Gordon,et al.  ωB97M-V: A combinatorially optimized, range-separated hybrid, meta-GGA density functional with VV10 nonlocal correlation. , 2016, The Journal of chemical physics.

[31]  Xiao He,et al.  MN15-L: A New Local Exchange-Correlation Functional for Kohn-Sham Density Functional Theory with Broad Accuracy for Atoms, Molecules, and Solids. , 2016, Journal of chemical theory and computation.

[32]  M. Head‐Gordon,et al.  Mapping the genome of meta-generalized gradient approximation density functionals: the search for B97M-V. , 2015, The Journal of chemical physics.

[33]  F. Marsiglio,et al.  The importance of basis states: an example using the Hydrogen basis , 2015, 1505.03510.

[34]  L. Kronik,et al.  A self-interaction-free local hybrid functional: accurate binding energies vis-à-vis accurate ionization potentials from Kohn-Sham eigenvalues. , 2014, The Journal of chemical physics.

[35]  Narbe Mardirossian,et al.  ωB97X-V: a 10-parameter, range-separated hybrid, generalized gradient approximation density functional with nonlocal correlation, designed by a survival-of-the-fittest strategy. , 2014, Physical chemistry chemical physics : PCCP.

[36]  Narbe Mardirossian,et al.  x B 97 XV : A 10-parameter , range-separated hybrid , generalized gradient approximation density functional with nonlocal correlation , designed by a survival-ofthe-fittest strategy , 2014 .

[37]  M. Head‐Gordon,et al.  Characterizing and Understanding the Remarkably Slow Basis Set Convergence of Several Minnesota Density Functionals for Intermolecular Interaction Energies. , 2013, Journal of chemical theory and computation.

[38]  Fabiano Corsetti,et al.  Optimal finite-range atomic basis sets for liquid water and ice , 2013, Journal of physics. Condensed matter : an Institute of Physics journal.

[39]  Frank Jensen,et al.  Atomic orbital basis sets , 2013 .

[40]  J. Grant Hill,et al.  Gaussian basis sets for molecular applications , 2013 .

[41]  D. Truhlar,et al.  Screened-exchange density functionals with broad accuracy for chemistry and solid-state physics. , 2012, Physical chemistry chemical physics : PCCP.

[42]  Donald G Truhlar,et al.  An improved and broadly accurate local approximation to the exchange-correlation density functional: the MN12-L functional for electronic structure calculations in chemistry and physics. , 2012, Physical chemistry chemical physics : PCCP.

[43]  G. Rothenberg,et al.  Transferable basis sets of numerical atomic orbitals , 2012 .

[44]  Donald G. Truhlar,et al.  M11-L: A Local Density Functional That Provides Improved Accuracy for Electronic Structure Calculations in Chemistry and Physics , 2012 .

[45]  Donald G. Truhlar,et al.  Improving the Accuracy of Hybrid Meta-GGA Density Functionals by Range Separation , 2011 .

[46]  Zhenyu Li,et al.  Linear scaling electronic structure calculations with numerical atomic basis set , 2010 .

[47]  Troy Van Voorhis,et al.  Nonlocal van der Waals density functional: the simpler the better. , 2010, The Journal of chemical physics.

[48]  S. Saito.  Hartree–Fock–Roothaan energies and expectation values for the neutral atoms He to Uuo: The B-spline expansion method , 2009 .

[49]  Matthias Scheffler,et al.  Ab initio molecular simulations with numeric atom-centered orbitals , 2009, Comput. Phys. Commun..

[50]  Fernando Nogueira,et al.  Generating relativistic pseudo-potentials with explicit incorporation of semi-core states using APE, the Atomic Pseudo-potentials Engine , 2008, Comput. Phys. Commun..

[51]  D. Truhlar,et al.  The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals , 2008 .

[52]  Xinghua Wang,et al.  On the Hermite interpolation , 2007 .

[53]  D. Truhlar,et al.  A new local density functional for main-group thermochemistry, transition metal bonding, thermochemical kinetics, and noncovalent interactions. , 2006, The Journal of chemical physics.

[54]  Donald G Truhlar,et al.  Design of Density Functionals by Combining the Method of Constraint Satisfaction with Parametrization for Thermochemistry, Thermochemical Kinetics, and Noncovalent Interactions. , 2006, Journal of chemical theory and computation.

[55]  Yan Zhao,et al.  Exchange-correlation functional with broad accuracy for metallic and nonmetallic compounds, kinetics, and noncovalent interactions. , 2005, The Journal of chemical physics.

[56]  T. Ozaki,et al.  Variationally optimized basis orbitals for biological molecules. , 2004, The Journal of chemical physics.

[57]  Taisuke Ozaki,et al.  Numerical atomic basis orbitals from H to Kr , 2004 .

[58]  G. Scuseria,et al.  Climbing the density functional ladder: nonempirical meta-generalized gradient approximation designed for molecules and solids. , 2003, Physical review letters.

[59]  Taisuke Ozaki,et al.  Variationally optimized atomic orbitals for large-scale electronic structures , 2003 .

[60]  L. Ram-Mohan.  Finite Element and Boundary Element Applications in Quantum Mechanics , 2002 .

[61]  J. Soler,et al.  Systematic generation of finite-range atomic basis sets for linear-scaling calculations , 2002, cond-mat/0207548.

[62]  Frank Jensen,et al.  Polarization consistent basis sets. II. Estimating the Kohn-Sham basis set limit , 2002 .

[63]  D. Sánchez-Portal,et al.  Numerical atomic orbitals for linear-scaling calculations , 2001, cond-mat/0104170.

[64]  H. Tatewaki,et al.  Chemically reliable uncontracted Gaussian-type basis sets for atoms H to Lr , 2000 .

[65]  Hideaki Fujitani,et al.  Transferable atomic-type orbital basis sets for solids , 2000 .

[66]  Andrew P. Horsfield,et al.  Efficient AB Initio Tight Binding , 1997 .

[67]  K. Burke,et al.  Generalized Gradient Approximation Made Simple [Phys. Rev. Lett. 77, 3865 (1996)] , 1997 .

[68]  Daniel Sánchez-Portal,et al.  Density‐functional method for very large systems with LCAO basis sets , 1997 .

[69]  Burke,et al.  Generalized Gradient Approximation Made Simple. , 1996, Physical review letters.

[70]  Seifert,et al.  Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon. , 1995, Physical review. B, Condensed matter.

[71]  Emilio San-Fabián,et al.  Automatic numerical integration techniques for polyatomic molecules , 1994 .

[72]  Ishikawa,et al.  Universal Gaussian basis set for accurate ab initio /P relat ivistic Dirac-Fock calculations. , 1993, Physical review. A, Atomic, molecular, and optical physics.

[73]  Wang,et al.  Accurate and simple analytic representation of the electron-gas correlation energy. , 1992, Physical review. B, Condensed matter.

[74]  O. Sankey,et al.  Ab initio multicenter tight-binding model for molecular-dynamics simulations and other applications in covalent systems. , 1989, Physical review. B, Condensed matter.

[75]  David Feller,et al.  Basis Set Selection for Molecular Calculations , 1986 .

[76]  R. Cookson,et al.  An explicit formula for the generalized Hermitian interpolation problem , 1985 .

[77]  A. Ruffa.  Continuum Wave Functions in the Calculation of Sums Involving Off-Diagonal Matrix Elements , 1973 .

[78]  W. Kohn,et al.  Self-Consistent Equations Including Exchange and Correlation Effects , 1965 .

[79]  A. Spitzbart.  A Generalization of Hermite's Interpolation Formula , 1960 .

[80]  P. Löwdin,et al.  Role of the Continuum in Superposition of Configurations , 1955 .

[81]  P. Dirac.  Note on Exchange Phenomena in the Thomas Atom , 1930, Mathematical Proceedings of the Cambridge Philosophical Society.

[82]  F. Bloch,et al.  Bemerkung zur Elektronentheorie des Ferromagnetismus und der elektrischen Leitfähigkeit , 1929 .