A new mixed self-consistent field procedure†

A new approach for the calculation of three-centre electronic repulsion integrals (ERIs) is developed, implemented and benchmarked in the framework of auxiliary density functional theory (ADFT). The so-called mixed self-consistent field (mixed SCF) divides the computationally costly ERIs in two sets: far-field and near-field. Far-field ERIs are calculated using the newly developed double asymptotic expansion as in the direct SCF scheme. Near-field ERIs are calculated only once prior to the SCF procedure and stored in memory, as in the conventional SCF scheme. Hence the name, mixed SCF. The implementation is particularly powerful when used in parallel architectures, since all RAM available are used for near-field ERI storage. In addition, the efficient distribution algorithm performs minimal intercommunication operations between processors, avoiding a potential bottleneck. One-, two- and three-dimensional systems are used for benchmarking, showing substantial time reduction in the ERI calculation for all of them. A Born–Oppenheimer molecular dynamics calculation for the Na+55 cluster is also shown in order to demonstrate the speed-up for small systems achievable with the mixed SCF.

[1]  G. G. Hall.  The molecular orbital theory of chemical valency VIII. A method of calculating ionization potentials , 1951, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[2]  W. Bocian,et al.  Structure of human insulin monomer in water/acetonitrile solution , 2008, Journal of biomolecular NMR.

[3]  Jürgen Gauss,et al.  An unconventional scf method for calculations on large molecules , 1986 .

[4]  Florian Janetzko,et al.  A MinMax self-consistent-field approach for auxiliary density functional theory. , 2009, The Journal of chemical physics.

[5]  S. H. Vosko,et al.  Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis , 1980 .

[6]  Evert Jan Baerends,et al.  Self-consistent molecular Hartree—Fock—Slater calculations I. The computational procedure , 1973 .

[7]  Benny G. Johnson,et al.  A simple yet powerful upper bound for Coulomb integrals , 1994 .

[8]  Asbjörn M. Burow,et al.  Resolution of identity approximation for the Coulomb term in molecular and periodic systems. , 2009, The Journal of chemical physics.

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

[10]  D. Salahub,et al.  QM/MM Calculations with deMon2k , 2015, Molecules.

[11]  A. Szabo,et al.  Modern quantum chemistry , 1982 .

[12]  S. F. Boys.  Electronic wave functions - I. A general method of calculation for the stationary states of any molecular system , 1950, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[13]  P. Calaminici,et al.  Density functional theory optimized basis sets for gradient corrected functionals: 3d transition metal systems. , 2007, The Journal of chemical physics.

[14]  Matthias Krack,et al.  AN ADAPTIVE NUMERICAL INTEGRATOR FOR MOLECULAR INTEGRALS , 1998 .

[15]  J. Connolly,et al.  On first‐row diatomic molecules and local density models , 1979 .

[16]  J. Mintmire,et al.  Fitting the Coulomb potential variationally in linear-combination-of-atomic-orbitals density-functional calculations , 1982 .

[17]  J. Mintmire,et al.  Local-density-functional methods in two-dimensionally periodic systems. Hydrogen and beryllium monolayers , 1982 .

[18]  Andreas M Köster,et al.  Efficient and reliable numerical integration of exchange-correlation energies and potentials. , 2004, The Journal of chemical physics.

[19]  D. Salahub,et al.  Contribution of high-energy conformations to NMR chemical shifts, a DFT-BOMD study. , 2013, Physical chemistry chemical physics : PCCP.

[20]  Kazuhiro Ishida,et al.  Efficient determination and characterization of transition states using ab-initio methods , 1977 .

[21]  A. Köster.  Efficient recursive computation of molecular integrals for density functional methods , 1996 .

[22]  Andreas M Köster,et al.  Calculation of exchange-correlation potentials with auxiliary function densities. , 2004, The Journal of chemical physics.

[23]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[24]  G. Diercksen,et al.  Methods in Computational Molecular Physics , 1983 .

[25]  M. Head‐Gordon,et al.  A multipole acceptability criterion for electronic structure theory , 1998 .

[26]  J. C. Slater,et al.  Self-Consistent-Field X α Cluster Method for Polyatomic Molecules and Solids , 1972 .

[27]  V. R. Saunders,et al.  A “Level–Shifting” method for converging closed shell Hartree–Fock wave functions , 1973 .

[28]  C. C. J. Roothaan,et al.  Self-Consistent Field Theory for Open Shells of Electronic Systems , 1960 .

[29]  Notker Rösch,et al.  Variational fitting methods for electronic structure calculations , 2010 .

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

[31]  A. Köster,et al.  Double asymptotic expansion of three-center electronic repulsion integrals. , 2013, The Journal of chemical physics.

[32]  E. Clementi.  Analytical Self‐Consistent Field Functions for Positive Ions. I. Isoelectronic Series with 2 to 10 Electrons , 1963 .

[33]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[34]  P. Calaminici,et al.  Heat capacities from Born–Oppenheimer molecular dynamics simulations: Al27+ and Al28+ , 2013 .

[35]  G. G. Stokes.  "J." , 1890, The New Yale Book of Quotations.

[36]  D. Salahub,et al.  Asymptotic Expansion for Electrostatic Embedding Integrals in QM/MM Calculations. , 2012, Journal of chemical theory and computation.

[37]  Clemens C. J. Roothaan,et al.  New Developments in Molecular Orbital Theory , 1951 .

[38]  B. M. Fulk.  MATH , 1992 .

[39]  J. Almlöf,et al.  Principles for a direct SCF approach to LICAO–MOab‐initio calculations , 1982 .

[40]  D. N. Laikov.  Fast evaluation of density functional exchange-correlation terms using the expansion of the electron density in auxiliary basis sets , 1997 .

[41]  Joseph E. Subotnik,et al.  Linear scaling density fitting. , 2006, The Journal of chemical physics.

[42]  Florian Janetzko,et al.  Parallelization of the deMon2k code , 2006, J. Comput. Chem..

[43]  L. Visscher,et al.  Accurate Coulomb Potentials for Periodic and Molecular Systems through Density Fitting. , 2014, Journal of chemical theory and computation.

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

[45]  A. Köster.  Hermite Gaussian auxiliary functions for the variational fitting of the Coulomb potential in density functional methods , 2003 .

[46]  Ronald H. Felton,et al.  A new computational approach to Slater’s SCF–Xα equation , 1975 .

[47]  Patrizia Calaminici,et al.  Robust and efficient density fitting. , 2009, The Journal of chemical physics.

[48]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[49]  Scarontefan Varga.  Long-range analysis of density fitting in extended systems† , 2008 .