Avoiding Negligible Shell Pairs and Quartets in Electronic Structure Calculations

We define a significant shell pair in an electronic structure calculation as one that generates at least one two-electron integral larger than a preset threshold. We define a significant shell quartet similarly. We then explore several methods for identifying nonsignificant pairs and quartets so that they can be avoided and computational efficiency improved. We find that the widely used Cauchy-Schwarz bound identifies most nonsignificant quartets but that the Hölder bound is slightly more powerful for identifying nonsignificant pairs.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Q-MP2-OS: Møller-Plesset Correlation Energy by Quadrature

We present a quadrature-based algorithm for computing the opposite-spin component of the MP2 correlation energy which scales quadratically with basis set size and is well-suited to large-scale parallelization. The key ideas, which are rooted in the earlier work of Hirata and co-workers, are to abandon all two-electron integrals, recast the energy as a seven-dimensional integral, approximate that integral by quadrature, and employ a cutoff strategy to minimize the number of intermediate quantities. We discuss our implementation in detail and show that it parallelizes almost perfectly on 840 cores for cyclosporine (a molecule with roughly 200 atoms), exhibits [Formula: see text] scaling for a sequence of polyglycines, and is principally limited by the accuracy of its quadrature.

Efficient Method for Calculating Effective Core Potential Integrals

Effective core potential (ECP) integrals are among the most difficult one-electron integrals to calculate due to the projection operators. The radial part of these operators may include r⁰, r^-1, and r^-2 terms. For the r⁰ terms, we exploit a simple analytic expression for the fundamental projected integral to derive new recurrence relations and upper bounds for ECP integrals. For the r^-1 and r^-2 terms, we present a reconstruction method that replaces these terms by a sum of r⁰ terms and show that the resulting errors are chemically insignificant for a range of molecular properties. The new algorithm is available in Q-Chem 5.0 and is significantly faster than the ECP implementations in Q-Chem 4.4, GAMESS (US) and Dalton 2016.

Accurate Electron Densities at Nuclei Using Small Ramp-Gaussian Basis Sets

Electron densities at nuclei are difficult to calculate accurately with all-Gaussian basis sets because they lack an electron-nuclear cusp. The newly developed mixed ramp-Gaussian basis sets, such as R-31G, possess electron-nuclear cusps due to the presence of ramp functions in the basis. The R-31G basis set is a general-purpose mixed ramp-Gaussian basis set modeled on the 6-31G basis set. The prediction of electron densities at nuclei using R-31G basis sets for Li-F outperforms Dunning, Pople, and Jensen general purpose all-Gaussian basis sets of triple-ζ quality or lower and the cc-pVQZ basis set. It is of similar quality to the specialized pcJ-0 basis set which was developed with partial decontraction of core functions and extra high exponent s-Gaussians to predict electron density at the nucleus. These results show significant advantages in the properties of mixed ramp-Gaussian basis sets compared to all-Gaussian basis sets.

Mixed Ramp-Gaussian Basis Sets

We discuss molecular orbital basis sets that contain both Gaussian and polynomial (ramp) functions. We show that, by modeling ramp-Gaussian products as sums of ramps, all of the required one- and two-electron integrals can be computed quickly and accurately. To illustrate our approach, we construct R-31+G, a mixed ramp-Gaussian basis in which the core basis functions of the 6-31+G basis are replaced by ramps. By performing self-consistent Hartree-Fock calculations, we show that the thermochemical predictions of R-31+G and 6-31+G are similar but the former has the potential to be significantly faster.

Interfacing Q-Chem and CHARMM to perform QM/MM reaction path calculations

A hybrid quantum mechanical/molecular mechanical (QM/MM) potential energy function with Hartree-Fock, density functional theory (DFT), and post-HF (RIMP2, MP2, CCSD) capability has been implemented in the CHARMM and Q-Chem software packages. In addition, we have modified CHARMM and Q-Chem to take advantage of the newly introduced replica path and the nudged elastic band methods, which are powerful techniques for studying reaction pathways in a highly parallel (i.e., parallel/parallel) fashion, with each pathway point being distributed to a different node of a large cluster. To test our implementation, a series of systems were studied and comparisons were made to both full QM calculations and previous QM/MM studies and experiments. For instance, the differences between HF, DFT, MP2, and CCSD QM/MM calculations of H2O...H2O, H2O...Na+, and H2O...Cl- complexes have been explored. Furthermore, the recently implemented polarizable Drude water model was used to make comparisons to the popular TIP3P and TIP4P water models for doing QM/MM calculations. We have also computed the energetic profile of the chorismate mutase catalyzed Claisen rearrangement at various QM/MM levels of theory and have compared the results with previous studies. Our best estimate for the activation energy is 8.20 kcal/mol and for the reaction energy is -23.1 kcal/mol, both calculated at the MP2/6-31+G(d)//MP2/6-31+G(d)/C22 level of theory.

BFW: A Density Functional for Transition Metal Clusters

Ionization potentials (IPs) or electron affinities (EAs) for transition metal clusters are an important property that can be used to identify and differentiate between clusters. Accurate calculation of these values is therefore vital. Previous attempts using a variety of DFT models have correctly predicted trends, but have relied on the use of scaling factors to compare to experimental IPs. In this paper, we introduce a new density functional (BFW) that is explicitly designed to yield accurate, absolute IPs for transition metal clusters. This paper presents the numerical results for a selection of transition metal clusters and their carbides, nitrides, and oxides for which experimental IPs are known. When tested on transition metal clusters, the BFW functional is found to be significantly more accurate than B3LYP and B3PW91.

Advances in methods and algorithms in a modern quantum chemistry program package

Advances in theory and algorithms for electronic structure calculations must be incorporated into program packages to enable them to become routinely used by the broader chemical community. This work reviews advances made over the past five years or so that constitute the major improvements contained in a new release of the Q-Chem quantum chemistry package, together with illustrative timings and applications. Specific developments discussed include fast methods for density functional theory calculations, linear scaling evaluation of energies, NMR chemical shifts and electric properties, fast auxiliary basis function methods for correlated energies and gradients, equation-of-motion coupled cluster methods for ground and excited states, geminal wavefunctions, embedding methods and techniques for exploring potential energy surfaces.

Decay behavior of least-squares coefficients in auxiliary basis expansions

Many quantum chemical methods, both wave function and density based, rely on an expansion of elements of the electron density in an auxiliary basis. However, little is known about the analytical behavior of the expansion coefficients and, in particular, about their rate of decay with distance. We discuss an exactly solvable model system and characterize the expansion coefficients for various fitting metrics and various dimensionalities of the auxiliary basis. In the case of Coulomb fitting, we find that the decay rate depends critically on the effective dimensionality D of the auxiliary basis, varying from O(r(-1)) to O(r(-3)) to O(e(-zetar)) for D = 1, 2, or 3.

Extracting atoms from molecular electron densities via integral equations

The observation that a molecular electron density is close to the superposition of its constituent atoms leads naturally to the idea of modeling a density by a sum of nuclear-centered, spherically symmetric functions. The functions that are optimal in a least-squares sense are known as Stewart atoms. Previous attempts to construct Stewart atoms by expanding them in an auxiliary basis have been thwarted by slow convergence with respect to the size of the auxiliary basis used. We present a method for constructing Stewart atoms via convolution integrals which bypasses the need for an auxiliary basis, and is able to produce highly accurate approximations to Stewart atoms.