scispace - formally typeset
Search or ask a question
Book ChapterDOI

An Introduction to Coupled Cluster Theory for Computational Chemists

05 Jan 2007-pp 33-136
About: The article was published on 2007-01-05. It has received 566 citations till now. The article focuses on the topics: Coupled cluster.
Citations
More filters
Journal ArticleDOI
TL;DR: 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.
Abstract: 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.

2,527 citations

Journal ArticleDOI
TL;DR: The Psi4 program is a new approach to modern quantum chemistry, encompassing Hartree–Fock and density‐functional theory to configuration interaction and coupled cluster and offers flexible user input built on the Python scripting language that enables both new and experienced users to make full use of the program's capabilities.
Abstract: The Psi4 program is a new approach to modern quantum chemistry, encompassing Hartree–Fock and density-functional theory to configuration interaction and coupled cluster. The program is written entirely in C++ and relies on a new infrastructure that has been designed to permit high-efficiency computations of both standard and emerging electronic structure methods on conventional and high-performance parallel computer architectures. Psi4 offers flexible user input built on the Python scripting language that enables both new and experienced users to make full use of the program's capabilities, and even to implement new functionality with moderate effort. To maximize its impact and usefulness, Psi4 is available through an open-source license to the entire scientific community. © 2011 John Wiley & Sons, Ltd.

902 citations


Cites background or methods from "An Introduction to Coupled Cluster ..."

  • ...The levels of SAPT termed SAPT2, SAPT2+, SAPT2+(3), and SAPT2+3 are all fully functional.37 The largest computations performed with higher-order SAPT are SAPT2+(3)/aug-cc-pVTZ computations on the adenine–thymine complex....

    [...]

  • ...D evelopers of modern quantum chemistry programs face three formidable challenges: (1) they must provide efficient code for established electronic structures methods, such as Hartree–Fock (HF), density functional theory, and advanced coupled cluster (CC) models, for ‘conventional’ computing hardware—ranging from multicore laptop and desktop computers to moderate-scale parallel clusters containing thousands of cores; (2) they must be prepared to implement new methods quickly and at ‘production level’; and (3) they must be ready to adapt their programs to newly emerging computing hardware, such as the cutting-edge massively parallel petaand exascale systems, or the nascent graphical processing unit (GPU) machines that represent an entirely new archetype for scientific computing....

    [...]

  • ...37 The largest computations performed with higher-order SAPT are SAPT2+(3)/aug-cc-pVTZ computations on the adenine–thymine complex....

    [...]

  • ...The levels of SAPT termed SAPT2, SAPT2+, SAPT2+(3), and SAPT2+3 are all fully functional....

    [...]

Journal ArticleDOI
TL;DR: This Review collects and emphasizes the importance of inter- and intramolecular dispersion for molecules consisting mostly of first row atoms and forces us to reconsider the authors' perception of steric hindrance and stereoelectronic effects.
Abstract: London dispersion, which constitutes the attractive part of the famous van der Waals potential, has long been underappreciated in molecular chemistry as an important element of structural stability, and thus affects chemical reactivity and catalysis. This negligence is due to the common notion that dispersion is weak, which is only true for one pair of interacting atoms. For increasingly larger structures, the overall dispersion contribution grows rapidly and can amount to tens of kcal mol(-1) . This Review collects and emphasizes the importance of inter- and intramolecular dispersion for molecules consisting mostly of first row atoms. The synergy of experiment and theory has now reached a stage where dispersion effects can be examined in fine detail. This forces us to reconsider our perception of steric hindrance and stereoelectronic effects. The quantitation of dispersion energy donors will improve our ability to design sophisticated molecular structures and much better catalysts.

614 citations

Journal ArticleDOI
TL;DR: It is shown that for medium sized molecules the total wall clock time required to complete the LPNO-CCSD calculations is only two to four times that of the preceding self-consistent field (SCF) and these methods are highly suitable for large-scale computational chemistry applications.
Abstract: A production level implementation of the closed-shell local quadratic configuration interaction and coupled cluster methods with single and double excitations (QCISD and CCSD) based on the concept of pair natural orbitals [local pair natural orbital LPNO-QCISD and LPNO-CCSD) is reported, evaluated, and discussed. This work is an extension of the earlier developed LPNO coupled-electron pair approximation (LNPO-CEPA) method [F. Neese et al., Chem. Phys. 130, 114108 (2009)] and makes extended use of the resolution of the identity (RI) or density fitting (DF) approximation. Two variants of each method are compared. The less accurate approximations (LPNO2-QCISD/LPNO2-CCSD) still recover 98.7%-99.3% of the correlation energy in the given basis and have modest disk space requirements. The more accurate variants (LPNO1-QCISD/LPNO1-CCSD) typically recover 99.75%-99.95% of the correlation energy in the given basis but require the Coulomb and exchange operators with up to two-external indices to be stored on disk. Both variants have comparable computational efficiency. The convergence of the results with respect to the natural orbital truncation parameter (T(CutPNO)) has been studied. Extended numerical tests have been performed on absolute and relative correlation energies as function of basis set size and T(CutPNO) as well as on reaction energies, isomerization energies, and weak intermolecular interactions. The results indicate that the errors of the LPNO methods compared to the canonical QCISD and CCSD methods are below 1 kcal/mol with our default thresholds. Finally, some calculations on larger molecules are reported (ranging from 40-86 atoms) and it is shown that for medium sized molecules the total wall clock time required to complete the LPNO-CCSD calculations is only two to four times that of the preceding self-consistent field (SCF). Thus these methods are highly suitable for large-scale computational chemistry applications. Since there are only three thresholds involved that have been given conservative default values, the methods can be confidentially used in a "black-box" fashion in the same way as their canonical counterparts.

467 citations

Journal ArticleDOI
TL;DR: A general-purpose neural network potential is trained that approaches CCSD(T)/CBS accuracy on benchmarks for reaction thermochemistry, isomerization, and drug-like molecular torsions.
Abstract: Computational modeling of chemical and biological systems at atomic resolution is a crucial tool in the chemist’s toolset. The use of computer simulations requires a balance between cost and accuracy: quantum-mechanical methods provide high accuracy but are computationally expensive and scale poorly to large systems, while classical force fields are cheap and scalable, but lack transferability to new systems. Machine learning can be used to achieve the best of both approaches. Here we train a general-purpose neural network potential (ANI-1ccx) that approaches CCSD(T)/CBS accuracy on benchmarks for reaction thermochemistry, isomerization, and drug-like molecular torsions. This is achieved by training a network to DFT data then using transfer learning techniques to retrain on a dataset of gold standard QM calculations (CCSD(T)/CBS) that optimally spans chemical space. The resulting potential is broadly applicable to materials science, biology, and chemistry, and billions of times faster than CCSD(T)/CBS calculations. Computational modelling of chemical systems requires a balance between accuracy and computational cost. Here the authors use transfer learning to develop a general purpose neural network potential that approaches quantum-chemical accuracy for reaction thermochemistry, isomerization, and drug-like molecular torsions.

400 citations

References
More filters
Journal ArticleDOI
TL;DR: In this paper, a detailed study of correlation effects in the oxygen atom was conducted, and it was shown that primitive basis sets of primitive Gaussian functions effectively and efficiently describe correlation effects.
Abstract: In the past, basis sets for use in correlated molecular calculations have largely been taken from single configuration calculations. Recently, Almlof, Taylor, and co‐workers have found that basis sets of natural orbitals derived from correlated atomic calculations (ANOs) provide an excellent description of molecular correlation effects. We report here a careful study of correlation effects in the oxygen atom, establishing that compact sets of primitive Gaussian functions effectively and efficiently describe correlation effects i f the exponents of the functions are optimized in atomic correlated calculations, although the primitive (s p) functions for describing correlation effects can be taken from atomic Hartree–Fock calculations i f the appropriate primitive set is used. Test calculations on oxygen‐containing molecules indicate that these primitive basis sets describe molecular correlation effects as well as the ANO sets of Almlof and Taylor. Guided by the calculations on oxygen, basis sets for use in correlated atomic and molecular calculations were developed for all of the first row atoms from boron through neon and for hydrogen. As in the oxygen atom calculations, it was found that the incremental energy lowerings due to the addition of correlating functions fall into distinct groups. This leads to the concept of c o r r e l a t i o n c o n s i s t e n t b a s i s s e t s, i.e., sets which include all functions in a given group as well as all functions in any higher groups. Correlation consistent sets are given for all of the atoms considered. The most accurate sets determined in this way, [5s4p3d2f1g], consistently yield 99% of the correlation energy obtained with the corresponding ANO sets, even though the latter contains 50% more primitive functions and twice as many primitive polarization functions. It is estimated that this set yields 94%–97% of the total (HF+1+2) correlation energy for the atoms neon through boron.

26,705 citations

Journal ArticleDOI
TL;DR: In this article, a perturbation theory for treating a system of n electrons in which the Hartree-Fock solution appears as the zero-order approximation was developed, and it was shown by this development that the first order correction for the energy and the charge density of the system is zero.
Abstract: A perturbation theory is developed for treating a system of n electrons in which the Hartree-Fock solution appears as the zero-order approximation. It is shown by this development that the first order correction for the energy and the charge density of the system is zero. The expression for the second-order correction for the energy greatly simplifies because of the special property of the zero-order solution. It is pointed out that the development of the higher approximation involves only calculations based on a definite one-body problem.

12,067 citations

Journal ArticleDOI
TL;DR: In this paper, a new augmented version of coupled-cluster theory, denoted as CCSD(T), is proposed to remedy some of the deficiencies of previous augmented coupledcluster models.

7,249 citations

Journal ArticleDOI
TL;DR: The coupled cluster singles and doubles model (CCSD) as discussed by the authors is derived algebraically, presenting the full set of equations for a general reference function explicitly in spin-orbital form, and the computational implementation of the CCSD model, which involves cubic and quartic terms, is discussed and results are compared with full CI calculations for H2O and BeH2.
Abstract: The coupled‐cluster singles and doubles model (CCSD) is derived algebraically, presenting the full set of equations for a general reference function explicitly in spin–orbital form. The computational implementation of the CCSD model, which involves cubic and quartic terms, is discussed and results are reported and compared with full CI calculations for H2O and BeH2. We demonstrate that the CCSD exponential ansatz sums higher‐order correlation effects efficiently even for BeH2, near its transition state geometry where quasidegeneracy efforts are quite large, recovering 98% of the full CI correlation energy. For H2O, CCSD plus the fourth‐order triple excitation correction agrees with the full CI energy to 0.5 kcal/mol. Comparisons with low‐order models provide estimates of the effect of the higher‐order terms T1T2, T21T2, T31, and T41 on the correlation energy.

5,603 citations

Journal ArticleDOI
TL;DR: In this article, a general procedure for calculation of the electron correlation energy, starting from a single Hartree-Fock determinant, is introduced, and the relation of this method to coupled-cluster (CCSD) theory is discussed.
Abstract: A general procedure is introduced for calculation of the electron correlation energy, starting from a single Hartree–Fock determinant. The normal equations of (linear) configuration interaction theory are modified by introducing new terms which are quadratic in the configuration coefficients and which ensure size consistency in the resulting total energy. When used in the truncated configuration space of single and double substitutions, the method, termed QCISD, leads to a tractable set of quadratic equations. The relation of this method to coupled‐cluster (CCSD) theory is discussed. A simplified method of adding corrections for triple substitutions is outlined, leading to a method termed QCISD(T). Both of these new procedures are tested (and compared with other procedures) by application to some small systems for which full configuration interaction results are available.

4,194 citations