Analytic gradients for density cumulant functional theory: The DCFT-06 model

Density cumulant functional theory (DCFT) is one of a number of nascent electron correlation methods that are derived from reduced density matrices and cumulants thereof, instead of the wavefunction. Deriving properties from the density cumulant naturally yields methods that are size extensive and size consistent. In this work, we derive expressions for the analytic gradient, with respect to an external perturbation, for the DCFT-06 variant of density cumulant functional theory. Despite the fact that the DCFT-06 energy functional is stationary with respect to the density cumulant, the analytic gradients of the energy require the solution of perturbation-independent equations for both orbital and cumulant response. These two sets of linear response equations are coupled in nature and are solved iteratively with the solution of orbital and cumulant response equations each macroiteration, exhibiting rapid convergence. The gradients are implemented and benchmarked against coupled cluster theory with single and double excitations (CCSD) and CCSD with perturbative triple excitations [CCSD(T)], as well as accurate empirically corrected experimental data, for a test set comprising 15 small molecules. For most of the test cases, results from DCFT-06 are closer to CCSD(T) and empirical data than those from CCSD. Although the total energy and analytic gradient have the same asymptotic scaling, the present experience shows that the computational cost of the gradient is significantly lower.     

Correction to DFT interaction energies by an empirical dispersion term valid for a range of intermolecular distances

The computation of intermolecular interaction energies via commonly used density functionals is hindered by their inaccurate inclusion of medium and long range dispersion interactions. Practical computation of inter- and intra-macrobiomolecule interaction energies, in particular, requires a fairly accurate yet not overly expensive methodology. It is also desirable to compute intermolecular energies not only at their equilibrium (lowest energy) configurations but also over a range of biophysically relevant distances. We present a method to compute intermolecular interaction energies by including an empirical correction for dispersion which is valid over a range of intermolecular distances. This is achieved by optimizing parameters that moderate the empirical correction by explicit comparison of density functional (B3LYP) energies with distance-dependent (DD) reference values obtained at the CCSD(T)/CBS limit. The resulting method, hereafter referred to as B3LYP-DD, yields interaction energies with an accuracy generally better than 1 kcal mol(-1) for different types of noncovalent complexes, over a range of intermolecular distances and interaction strengths, relative to the expensive CCSD(T)/CBS standard. For a training set of dispersion interacting complexes, B3LYP-DD interaction energies in combination with diffuse functions display absolute errors equal to or smaller than 0.68 kcal mol(-1). The empirical correction does not significantly increase the computational cost as compared to standard density functional calculations. Applications relevant to biomolecular energy and structure, such as prediction of DNA base-pair interactions, are also presented.     

One-electron properties and electrostatic interaction energies from the expectation value expression and wave function of singles and doubles coupled cluster theory

One-electron density matrices resulting from the explicitly connected commutator expansion of the expectation value were implemented at the singles and doubles coupled cluster (CCSD) level. In the proposed approach the one-electron density matrix is obtained at a little extra cost in comparison to the calculation of the CCSD correlation energy. Therefore, in terms of the computational time the new method is significantly less demanding than the conventional linear-response CCSD theory which requires additionally an expensive calculation of the left-hand solution of the CCSD equations. The quality of the new density matrices was investigated by computing a set of one-electron properties for a series of molecules of varying sizes and comparing the results with data obtained using the full configuration interaction method or higher level coupled cluster theory. It has been found that the results obtained using the new approach are of the same quality as those predicted by the linear-response CCSD method. The novel one-electron density matrices have also been applied to study the energy of the electrostatic interaction for a number of van der Waals complexes, including the benzene and azulene dimers.     

Toward a less costly but accurate calculation of the CCSD(T)/CBS noncovalent interaction energy

The popular method of calculating the noncovalent interaction energies at the coupled-cluster single-, double-, and perturbative triple-excitations [CCSD(T)] theory level in the complete basis set (CBS) limit was to add a CCSD(T) correction term to the CBS second-order Møller-Plesset perturbation theory (MP2). The CCSD(T) correction term is the difference between the CCSD(T) and MP2 interaction energies evaluated in a medium basis set. However, the CCSD(T) calculations with the medium basis sets are still very expensive for systems with more than 30 atoms. Comparatively, the domain-based local pair natural orbital coupled-cluster method [DLPNO-CCSD(T)] can be applied to large systems with over 1,000 atoms. Considering both the computational accuracy and efficiency, in this work, we propose a new scheme to calculate the CCSD(T)/CBS interaction energies. In this scheme, the MP2/CBS term keeps intact and the CCSD(T) correction term is replaced by a DLPNO-CCSD(T) correction term which is the difference between the DLPNO-CCSD(T) and DLPNO-MP2 interaction energies evaluated in a medium basis set. The interaction energies of the noncovalent systems in the S22, HSG, HBC6, NBC10, and S66 databases were recalculated employing this new scheme. The consistent and tight settings of the truncation parameters for DLPNO-CCSD(T) and DLPNO-MP2 in this noncanonical CCSD(T)/CBS calculations lead to the maximum absolute deviation and root-mean-square deviation from the canonical CCSD(T)/CBS interaction energies of less than or equal to 0.28 kcal/mol and 0.09 kcal/mol, respectively. The high accuracy and low cost of this new computational scheme make it an excellent candidate for the study of large noncovalent systems.     

An efficient algorithm for energy gradients and orbital optimization in valence bond theory

An efficient algorithm for energy gradients in valence bond theory with nonorthogonal orbitals is presented. A general Hartree-Fock-like expression for the Hamiltonian matrix element between valence bond (VB) determinants is derived by introducing a transition density matrix. Analytical expressions for the energy gradients with respect to the orbital coefficients are obtained explicitly, whose scaling for computational cost is m(4), where m is the number of basis functions, and is thus approximately the same as in HF method. Compared with other existing approaches, the present algorithm has lower scaling, and thus is much more efficient. Furthermore, the expression for the energy gradient with respect to the nuclear coordinates is also presented, and it provides an effective algorithm for the geometry optimization and the evaluation of various molecular properties in VB theory. Test applications show that our new algorithm runs faster than other methods.     

Variational formulation of perturbative explicitly-correlated coupled-cluster methods

We present a variational formulation of the recently-proposed CCSD(2)(R12) method [Valeev, Phys. Chem. Chem. Phys., 2008, 10, 106]. The centerpiece of this approach is the CCSD(2)(R12) Lagrangian obtained via L?wdin partitioning of the coupled-cluster singles and doubles (CCSD) Hamiltonian. Extremization of the Lagrangian yields the second-order basis set incompleteness correction for the CCSD energy. We also developed a simpler Hylleraas-type functional that only depends on one set of geminal amplitudes by applying screening approximations. This functional is used to develop a diagonal orbital-invariant version of the method in which the geminal amplitudes are fixed at the values determined by the first-order cusp conditions. Extension of the variational method to include perturbatively the effect of connected triples produces the method that approximates the complete basis-set limit of the standard CCSD plus perturbative triples [CCSD(T)] method. For a set of 20 small closed-shell molecules, the method recovered at least 94.5/97.3% of the CBS CCSD(T) correlation energy with the aug-cc-pVDZ/aug-cc-pVTZ orbital basis set. For 12 isogyric reactions involving these molecules, combining the aug-cc-pVTZ correlation energies with the aug-cc-pVQZ Hartree-Fock energies produces the electronic reaction energies with a mean absolute deviation of 1.4 kJ mol(-1) from the experimental values. The method has the same number of optimized parameters as the corresponding CCSD(T) model, does not require any modification of the coupled-cluster computer program, and only needs a small triple-zeta basis to match the precision of the considerably more expensive standard quintuple-zeta CCSD(T) computation.     

A truncated version of reduced multireference coupled-cluster method with singles and doubles and noniterative triples: application to F2 and Ni(CO)n (n=1, 2, and 4)

A perturbatively truncated version of the reduced multireference coupled-cluster method with singles and doubles and noniterative triples RMR CCSD(T) is described. In the standard RMR CCSD method, the effect of all triples and quadruples that are singles or doubles relative to references spanning a chosen multireference (MR) model space is accounted for via the external corrections based on the MR CISD wave function. In the full version of RMR CCSD(T), the remaining triples are then handled via perturbative corrections as in the standard, single-reference (SR) CCSD(T) method. By using a perturbative threshold in the selection of MR CISD configuration space, we arrive at the truncated version of RMR CCSD(T), in which the dimension of the MR CISD problem is significantly reduced, thus leaving more triples to be treated perturbatively. This significantly reduces the computational cost. We illustrate this approach on the F2 molecule, in which case the computational cost of the truncated version of RMR CCSD(T) is only about 10%-20% higher than that of the standard CCSD(T), while still eliminating the failure of CCSD(T) in the bond breaking region of geometries. To demonstrate the capabilities of the method, we have also used it to examine the structure and binding energy of transition metal complexes Ni(CO)n with n=1, 2, and 4. In particular, Ni(CO)2 is shown to be bent rather than linear, as implied by some earlier studies. The RMR CCSD(T) binding energy differs from the SR CCSD(T) one by 1-2 kcal/mol, while the energy barrier separating the linear and bent structures of Ni(CO)2 is smaller than 1 kcal/mol.     

Approximately size extensive local Multireference Singles and Doubles Configuration Interaction

Multi-reference Configuration Interaction (MRCI) is often used to predict the electronic structures and reaction energetics of small molecules with very high accuracy. Unfortunately, MRCI is inapplicable to large or even medium-sized molecules for two reasons: its computational cost scales poorly with molecule size and MRCI methods are not size extensive, leading to large energy errors. We have developed a new local (L) and approximately size extensive MRCI method that addresses both shortcomings. Truncating long-range electron correlation in a local orbital basis as well as efficient processing of two-electron integrals via Cholesky decomposition (CD) and integral screening reduce the computational cost to O(N(3)) with a small prefactor. A priori and a posteriori size extensivity corrections are both considered. The Multi-reference Averaged Coupled-Pair Functional (MRACPF) provides approximate size extensivity by modifying the energy functional. The very inexpensive Davidson-Silver and Pople a posteriori schemes also produce quite accurate corrections over a large range of molecular size. Employing CD-LMRACPF is slightly more expensive than using a Davidson-type correction, but the former gives superior results. Molecules with up to 50 heavy atoms can be treated with our CD-LMRACPF method thus far.     

Improved correlation energy extrapolation schemes based on local pair natural orbital methods

It is well-known that the basis set limit is difficult to reach in correlated post Hartree-Fock ab initio calculations. One possible route forward is to employ basis set extrapolation schemes. In order to avoid prohibitively expensive calculations, the highest level calculation (typically based on the "gold standard" coupled cluster theory with single, double, and perturbative triple excitations, CCSD(T)) is only performed with the smallest basis set, and the remaining basis set incompleteness is estimated at a lower level of theory, typically second-order M?ller-Plesset perturbation theory (MP2). In this work, we provide a comprehensive investigation of alternative schemes where the MP2 extrapolation is replaced by the coupled-electron pair approximation, version 1 (CEPA/1) or the local pair natural orbital version of this method (LPNO-CEPA/1). It is shown that the MP2 method achieves apparent accuracy only due to error cancellation. Systematically more accurate results at small additional computational cost are obtained if the MP2 step is replaced by LPNO-CEPA/1. The errors of LPNO-CEPA/1 relative to canonical CEPA/1 are negligible. Owing to the highly systematic nature of the deviations between canonical and LPNO methods, basis set extrapolation reduces the LPNO errors in the total energies by 1 order of magnitude (~0.2 kcal/mol) and errors in energy differences to essentially zero. Using the CCSD(T)/LPNO-CEPA/1-based extrapolation scheme, new reference values are proposed for the recently published S66 set of interaction energies. The deviations between the new values and the original interactions energies are mostly very small but reach values up to 0.3 kcal/mol.     

Quintuple-zeta quality coupled-cluster correlation energies with triple-zeta basis sets

The explicitly-correlated coupled-cluster method CCSD(T)(R12) is extended to include F12 geminal basis functions that decay exponentially with the interelectronic distance and reproduce the form of the average Coulomb hole more accurately than linear-r(12). Equations derived using the Ansatz 2 strong orthogonality projector are presented. The convergence of the correlation energy with orbital basis set for the new CCSD(T)(F12) method is studied and found to be rapid, 98% of the basis set limit correlation energy is typically recovered using triple-zeta orbital basis sets. The performance for reaction enthalpies is assessed via a test set of 15 reactions involving 23 molecules. The title statement is found to hold equally true for total and relative correlation energies.     

Energy screening for the incremental scheme: application to intermolecular interactions

A systematic screening procedure for small contributions in the incremental expansion of the correlation energy is presented. The performance of the proposed scheme is checked for the calculation of intermolecular interactions in realistic test systems as large as a guanine-cytosine base pair. It is found that the computational cost for the incremental expansion can be reduced considerably without significant loss of accuracy. Typically, the errors of the systems investigated here amount to <3.4, 0.22, and 0.06% for second-, third-, and fourth-order expansions, respectively. For almost all cases, the error in the total correlation energy can be kept below 1 kcal/mol with respect to the canonical CCSD result if the incremental series is truncated in a proper way.     

An efficient closed-shell singles and doubles coupled-cluster method

A reformulated set of equations for the closed-shell singles and doubles coupled-cluster (CCSD) method is presented. A computational cost of n_v⁴n₀²+7n_v³n₀³+1n_v²n₀⁴ for the n⁶ steps is obtained, where n_v is the number of virtual molecular orbitals included in the CCSD procedure, n₀ is the number of doubly occupied molecular orbitals and n=n₀+n_v. Test calculations for the cis and trans isomers of FNNF and planar and pyramidal CH₃⁻ are presented. Equilibrium structures determined with large Gaussian basis sets at the second-order Møller-Plesset (MP2) perturbation level of theory are reported and used for the other electron correlation methods. With the largest one-particle basis set (144 contracted Gaussian functions), the equilibrium geometries of cis- and trans-FNNF agree with experiment. Based on analyses of planar and pyramidal CH₃⁻ wavefunctions and the calculated inversion barrier, it is suggested that the molecular anion may not exist in a planar configuration but that autodetachment of an electron occurs before the transition state is reached. Comparisons of our new CCSD procedure demonstrate that coupled-cluster methods are not significantly more expensive than similar electron correlation techniques.     

Rapid and accurate assessment of GPCR–ligand interactions Using the fragment molecular orbital‐based density‐functional tight‐binding method

The reliable and precise evaluation of receptor–ligand interactions and pair‐interaction energy is an essential element of rational drug design. While quantum mechanical (QM) methods have been a promising means by which to achieve this, traditional QM is not applicable for large biological systems due to its high computational cost. Here, the fragment molecular orbital (FMO) method has been used to accelerate QM calculations, and by combining FMO with the density‐functional tight‐binding (DFTB) method we are able to decrease computational cost 1000 times, achieving results in seconds, instead of hours. We have applied FMO‐DFTB to three different GPCR–ligand systems. Our results correlate well with site directed mutagenesis data and findings presented in the published literature, demonstrating that FMO‐DFTB is a rapid and accurate means of GPCR–ligand interactions. © 2017 Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

Assessment of the orbital‐optimized coupled‐electron pair theory for thermochemistry and kinetics: Improving on CCSD and CEPA(1)

An assessment of the orbital‐optimized coupled‐electron pair theory [or simply “optimized CEPA(0),” OCEPA(0)] [Bozkaya and Sherrill, J. Chem. Phys. 2013, 139, 054104] for thermochemistry and kinetics is presented. The OCEPA(0) method is applied to closed‐ and open‐shell reaction energies, barrier heights, and radical stabilization energies (RSEs). The performance of OCEPA(0) is compared with those of the MP2, CEPA(0), OCEPA(0), CEPA(1), coupled‐cluster singles and doubles (CCSD), and CCSD(T) methods [at complete basis set limits employing cc‐pVTZ and cc‐pVQZ basis sets]. For the most of the test sets, the OCEPA(0) method performs better than CEPA(0), CEPA(1), and CCSD, and provides accurate results. Especially, for open‐shell reaction energies and barrier heights, the OCEPA(0)–CEPA(1) and OCEPA(0)–CCSD differences become obvious. Similarly, for barrier heights and RSEs, the OCEPA(0) method improves on CEPA(0) by 1.6 and 2.3 kcal mol^?1. Our results demonstrate that the CEPA(0) method dramatically fails when the reference wave function suffers from the spin‐contamination problem. Conversely, the OCEPA(0) method can annihilate spin‐contamination in the unrestricted‐Hartree–Fock initial guess orbitals and can yield stable solutions. For overall evaluation, we conclude that the OCEPA(0) method is quite helpful not only for problematic open‐shell systems and transition‐states but also for closed‐shell molecules. Hence, one may prefer OCEPA(0) over CEPA(0), CEPA(1), and CCSD as an method, where N is the number of basis functions, for thermochemistry and kinetics. As discussed previously, the cost of the OCEPA(0) method is as much as of CCSD and CEPA(1) for energy computations. However, for analytic gradient computations, the OCEPA(0) method is two times less expensive than CCSD and CEPA(1). Further, the stationary properties of the OCEPA(0) method making it promising for excited state properties via linear response theory. © 2014 Wiley Periodicals, Inc.     

Improving the accuracy of ab initio methods with summation approximants and singularity analysis

The accuracy of Møller–Plesset (MP) perturbation theory and coupled‐cluster (CC) theory can be significantly improved, at essentially no increase in computational cost, by using summation approximants that model the way in which these theories converge to the full configuration interaction limit. Approximants for MP4 and CCSD(T) are presented, their size scaling is analyzed, and the functional analysis of the MP energy, on which the MP4 approximant is based, is discussed. The MP approximants are shown to have a form that is appropriate for describing resonance energies. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2003     

The formulation and performance of a perturbative correction to the perfect quadruples model

A recently published alternative hierarchy of coupled-cluster approximations is reformulated as a perturbative correction. A single variant, a model for the total electronic energy based on the perfect quadruples model, is explored in detail. The computational scaling of the method developed is the same as canonical second order Mo?ller-Plesset perturbation theory (fifth order in the number of molecular orbitals), but its accuracy competes with the high-accuracy, high-cost standard CCSD(T), even when the latter is allowed to break spin-symmetry. The variation presented can be implemented without explicit calculation and storage of the most expensive energy contributions, thereby improving the range of systems which can be treated. The performance and scaling of the method are demonstrated with calculations on the water, fluorine, and oxirane molecules, and compared to the parent model.     

Quantum control of molecular motion including electronic polarization effects with a two-stage toolkit

A method for incorporating strong electric field polarization effects into optimal control calculations is presented. A Born-Oppenheimer-type separation, referred to as the electric-nuclear Born-Oppenheimer (ENBO) approximation, is introduced in which variations of both the nuclear geometry and the external electric field are assumed to be slow compared with the speed at which the electronic degrees of freedom respond to these changes. This assumption permits the generation of a potential energy surface that depends not only on the relative geometry of the nuclei but also on the electric field strength and on the orientation of the molecule with respect to the electric field. The range of validity of the ENBO approximation is discussed in the paper. A two-stage toolkit implementation is presented to incorporate the polarization effects and reduce the cost of the optimal control dynamics calculations. As an illustration of the method, it is applied to optimal control of vibrational excitation in a hydrogen molecule aligned along the field direction. Ab initio configuration interaction calculations with a large orbital basis set are used to compute the H-H interaction potential in the presence of the electric field. The significant computational cost reduction afforded by the toolkit implementation is demonstrated.     

Explicit correlation and intermolecular interactions: investigating carbon dioxide complexes with the CCSD(T)-F12 method

We have optimized the lowest energy structures and calculated interaction energies for the CO(2)-Ar, CO(2)-N(2), CO(2)-CO, CO(2)-H(2)O, and CO(2)-NH(3) dimers with the recently developed explicitly correlated coupled cluster singles doubles and perturbative triples [CCSD(T)]-F12 methods and the associated VXZ-F12 (where X = D,T,Q) basis sets. For a given cardinal number, we find that results obtained with the CCSD(T)-F12 methods are much closer to the CCSD(T) complete basis set limit than the conventional CCSD(T) results. The relatively modest increase in the computational cost between explicit and conventional CCSD(T) is more than compensated for by the impressive accuracy of the CCSD(T)-F12 method. We recommend use of the CCSD(T)-F12 methods in combination with the VXZ-F12 basis sets for the accurate determination of equilibrium geometries and interaction energies of weakly bound electron donor acceptor complexes.