Ergo: An open-source program for linear-scaling electronic structure calculations

Abstract Ergo is a C++ program for all-electron Hartree–Fock and Kohn–Sham density functional theory electronic structure calculations using Gaussian basis sets. The program uses algorithms for which the computational cost increases linearly with system size for all parts of the calculation, including computation of the Fock/Kohn–Sham matrix and density matrix construction. Both spin-restricted and unrestricted calculations are supported, and both pure and hybrid density functionals. The program also supports linear-scaling computation of highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) vectors. This paper briefly describes how the code is organized and provides examples of how it can be used.

[1]  M. Ratner.  Molecular electronic-structure theory , 2000 .

[2]  Frank Neese,et al.  Sparse maps--A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory. , 2016, The Journal of chemical physics.

[3]  Emanuel H. Rubensson,et al.  Locality-aware parallel block-sparse matrix-matrix multiplication using the Chunks and Tasks programming model , 2015, Parallel Comput..

[4]  Elias Rudberg,et al.  Difficulties in applying pure Kohn–Sham density functional theory electronic structure methods to protein molecules , 2011, Journal of physics. Condensed matter : an Institute of Physics journal.

[5]  A. Akimov,et al.  Large-Scale Computations in Chemistry: A Bird's Eye View of a Vibrant Field. , 2015, Chemical reviews.

[6]  Jun Li,et al.  Basis Set Exchange: A Community Database for Computational Sciences , 2007, J. Chem. Inf. Model..

[7]  P. Pulay.  Improved SCF convergence acceleration , 1982 .

[8]  Alexey V. Akimov,et al.  Libra: An open‐Source “methodology discovery” library for quantum and classical dynamics simulations , 2016, J. Comput. Chem..

[9]  Emanuel H. Rubensson,et al.  Hartree-Fock calculations with linearly scaling memory usage. , 2008, The Journal of chemical physics.

[10]  Paweł Sałek,et al.  Nonlocal exchange interaction removes half-metallicity in graphene nanoribbons. , 2007, Nano letters.

[11]  Trygve Helgaker,et al.  Hartree–Fock and Kohn–Sham atomic-orbital based time-dependent response theory , 2000 .

[12]  Emanuel H. Rubensson,et al.  Bringing about matrix sparsity in linear‐scaling electronic structure calculations , 2011, J. Comput. Chem..

[13]  Emanuel H. Rubensson,et al.  On-the-Fly Computation of Frontal Orbitals in Density Matrix Expansions. , 2017, Journal of chemical theory and computation.

[14]  Emanuel H. Rubensson,et al.  Parameterless Stopping Criteria for Recursive Density Matrix Expansions. , 2015, Journal of chemical theory and computation.

[15]  Emanuel H. Rubensson,et al.  Chunks and Tasks: A programming model for parallelization of dynamic algorithms , 2012, Parallel Comput..

[16]  Abdul-Rahman Allouche,et al.  Gabedit—A graphical user interface for computational chemistry softwares , 2011, J. Comput. Chem..

[17]  Dmitri A Romanov,et al.  A time-dependent Hartree-Fock approach for studying the electronic optical response of molecules in intense fields. , 2005, Physical chemistry chemical physics : PCCP.

[18]  Paweł Sałek,et al.  Efficient implementation of the fast multipole method. , 2006, The Journal of chemical physics.

[19]  D. Bowler,et al.  O(N) methods in electronic structure calculations. , 2011, Reports on progress in physics. Physical Society.

[20]  David Robinson,et al.  A QM/MM study of the nature of the entatic state in plastocyanin , 2016, J. Comput. Chem..

[21]  Paweł Sałek,et al.  A linear scaling study of solvent-solute interaction energy of drug molecules in aqua solution. , 2007, The journal of physical chemistry. B.

[22]  Emanuel H. Rubensson,et al.  Density matrix purification with rigorous error control. , 2008, The Journal of chemical physics.

[23]  Nicholas D M Hine,et al.  Electrostatic considerations affecting the calculated HOMO–LUMO gap in protein molecules , 2013, Journal of physics. Condensed matter : an Institute of Physics journal.

[24]  Emanuel H. Rubensson,et al.  A hierarchic sparse matrix data structure for large‐scale Hartree‐Fock/Kohn‐Sham calculations , 2007, J. Comput. Chem..

[25]  Emanuel H. Rubensson,et al.  Automatic Selection of Integral Thresholds by Extrapolation in Coulomb and Exchange Matrix Constructions. , 2009, Journal of chemical theory and computation.

[26]  Emanuel H. Rubensson,et al.  Truncation of small matrix elements based on the Euclidean norm for blocked data structures , 2009, J. Comput. Chem..

[27]  Emanuel H. Rubensson,et al.  Canonical density matrix perturbation theory. , 2015, Physical review. E, Statistical, nonlinear, and soft matter physics.

[28]  Emanuel H. Rubensson,et al.  Kohn-Sham Density Functional Theory Electronic Structure Calculations with Linearly Scaling Computational Time and Memory Usage. , 2011, Journal of chemical theory and computation.

[29]  Emanuel H. Rubensson,et al.  Rotations of occupied invariant subspaces in self-consistent field calculations , 2008 .

[30]  Thierry Deutsch,et al.  Challenges in large scale quantum mechanical calculations , 2016, 1609.00252.

[31]  Emanuel H. Rubensson,et al.  Nonmonotonic Recursive Polynomial Expansions for Linear Scaling Calculation of the Density Matrix. , 2010, Journal of chemical theory and computation.