A refined cluster-in-molecule local correlation approach for predicting the relative energies of large systems

A refined cluster-in-molecule (CIM) method for local correlation calculations of large molecules is presented. In the present work, two new strategies are introduced to further improve the CIM approach: (1) Some medium-range electron correlation energies, which are neglected in the previous CIM approach, are taken into account. (2) A much simpler procedure using only a distance threshold is used to construct various clusters. To cover the medium-range correlation effect as much as possible, some two-atom-centered clusters are built, in addition to one-atom-centered clusters. Our test calculations at the second order perturbation theory (MP2) level show that the refined CIM method can recover about 99.9% of the conventional MP2 correlation energy using an appropriate distance threshold. The accuracy of the present CIM method is capable of providing reliable relative energies of medium-sized systems such as polyalanines with 10 residues, and water molecules with 50 water molecules. For polyalanines with up to 30 residues, we have demonstrated that the computational cost of the CIM-MP2 calculation increases linearly with the molecular size, but the required memory and disc-space do not need to increase for large systems. The improved CIM method has been used to compute the relative energy of ice-like (H(2)O)(96) clusters (with 2400 basis functions) and to predict the dimerization energy of a double-helical foldamer (with 2330 basis functions). The present CIM method is expected to be a practical local correlation method for describing the relative energies of large systems.     

Fully optimized implementation of the cluster-in-molecule local correlation approach for electron correlation calculations of large systems

A fully optimized implementation of the cluster-in-molecule (CIM) local correlation method for faster and more accurate electron correlation calculations of large systems is reported. The speedup comes from the new procedure of constructing virtual localized molecular orbitals of clusters. In the new procedure, Boughton–Pulay projection method is employed instead of a much more expensive Boys localization procedure. In addition, basis set superposition error correction for binding energy calculations and parallelized electron correlation calculations of clusters are now implemented. Benchmark calculations and illustrative applications at the Møller–Plesset perturbation theory, coupled cluster singles and doubles (CCSD), and CCSD with perturbative triples correction levels show that this newly optimized CIM approach is a reliable theoretical tool for electron correlation calculations of various large chemical systems. © 2018 Wiley Periodicals, Inc.     

An efficient implementation of the "cluster-in-molecule" approach for local electron correlation calculations

An efficient implementation of the "cluster-in-molecule" (CIM) approach is presented for performing local electron correlation calculations in a basis of orthogonal occupied and virtual localized molecular orbitals (LMOs). The main idea of this approach is that significant excitation amplitudes can be approximately obtained by solving the coupled cluster (or Moller-Plesset perturbation theory) equations of a series of "clusters," each of which contains a subset of occupied and virtual LMOs. In the present implementation, we have proposed a simple approach for constructing virtual LMOs of clusters, and new ways of constructing clusters and extracting the correlation contributions from calculations on clusters, which are more efficient than those suggested in the original work. More importantly, linear scaling of computational time of the CIM approach is achieved by evaluating the transformed two-electron integrals over LMOs using simple truncation techniques in limited operations (independent of the molecular size). With typical thresholds, for a variety of molecules our test calculations demonstrate that more than 99% of the conventional MP2 or coupled cluster with doubles correlation energies can be recovered in the present CIM approach.     

Large scale polarizability calculations using the approximate coupled cluster model CC2 and MP2 combined with the resolution-of-the-identity approximation

We present an implementation of static and frequency-dependent polarizabilities for the approximate coupled cluster singles and doubles model CC2 and static polarizabilities for second-order Mo?ller-Plesset perturbation theory. Both are combined with the resolution-of-the-identity approximation for electron repulsion integrals to achieve unprecedented low operation counts, input-output, and disc space demands. To avoid the storage of double excitation amplitudes during the calculation of derivatives of density matrices, we employ in addition a numerical Laplace transformation for orbital energy denominators. It is shown that the error introduced by this approximation is negligible already with a small number of sampling points. Thereby an implementation of second-order one-particle properties is realized, which avoids completely the storage of quantities scaling with the fourth power of the system size. The implementation is tested on a set of organic molecules including large fused aromatic ring systems and the C(60) fullerene. It is demonstrated that exploiting symmetry and shared memory parallelization, second-order properties for such systems can be evaluated at the CC2 and MP2 level within a few hours of calculation time. As large scale applications, we present results for the 7-, 9-, and 11-ring helicenes.     

Comparison of the Hartree-Fock,Møller-Plesset,and Hartree-Fock-Slater method with respect to electrostatic properties of small molecules

Summary The results of various quantum chemical calculations, the Hartree-Fock (HF) method, the Møller-Plesset perturbation theory (MP2), and the Hartree-Fock-Slater (HFS) method are compared. Atomic charges, dipole moments, topological properties of the electron density distribution and polarizabilities, and first hyperpolarizabilities are calculated. Atomic charges obtained with the HFS method are found to be very close to those calculated with the MP2 method, from which we conclude that the HFS method describes to some extent electron correlation effects. Performing an MP2 calculation after an HF calculation improves the molecular dipole moments considerably, yielding values close to the experimental ones. HFS calculations are computationally less demanding than MP2 and yield comparable results for the electron density distributions, dipole moments and polarizabilities.     

An overview of coupled cluster theory and its applications in physics

Summary What has since become known as the normal coupled cluster method (NCCM) was invented about thirty years ago to calculate ground-state energies of closed-shell atomic nuclei. Coupled cluster (CC) techniques have since been developed to calculate excited states, energies of open-shell systems, density matrices and hence other properties, sum rules, and the sub-sum-rules that follow from imbedding linear response theory within the NCCM. Further extensions deal both with systems at nonzero temperature and with general dynamical behaviour. More recently, a new version of CC theory, the so-called extended coupled cluster method (ECCM) has been introduced. It has the potential to describe such global phenomena as phase transitions, spontaneous symmetry breaking, states of topological excitation, and nonequilibrium behaviour. CC techniques are now widely recognized as providing one of the most universally applicable, most powerful, and most accurate of all microscopicab initio methods in quantum many-body theory. The number of successful applications within physics is now impressively large. In most such cases the numerical results are either the best or among the best available. A typical case is the electron gas, where the CC results for the correlation energy agree over the entire metallic density range to within less than 1 millihartree (or <1%) with the essentially exact Green's function Monte Carlo results. The role of CC theory within modern quantum many-body theory is first surveyed, by a comparison with other techniques. Its full range of applications in physics is then reviewed. These include problems in nuclear physics, both for finite nuclei and infinite nuclear matter; the electron gas; various integrable and nonintegrable models; various relativistic quantum field theories; and quantum spin chain and lattice models. Particular applications of the ECCM include the quantum hydrodynamics of a zero-temperature, strongly-interacting condensed Bose fluid; a charged impurity in a polarizable medium (e.g., positron annihilation in metals); and various anharmonic oscillator and spin systems.     

Static dipole polarizability and binding energy of sodium clusters Nan (n=1-10): A critical assessment of all-electron based post Hartree-Fock and density functional methods

A systematic all electron post Hartree-Fock as well as density functional theory (DFT) based calculations for the polarizability and binding energy of sodium metal clusters have been performed and an in-depth analysis of the discrepancy between the experimental and theoretical results is presented. A systematic investigation for the assessment of different DFT exchange-correlation functionals in predicting the polarizability values has also been reported. All the pure DFT functionals have been found to considerably underestimate the calculated polarizability values as compared to the MP2 results. DFT calculations using the full Hartree-Fock exchange along with one-parameter progressive correlation functional have, however, been shown to yield results in good agreement with the MP2 and experimental results. The possible sources of error present in the experimental measurements as well as in the different theoretical methods have also been analyzed. One of the most important conclusions of the present study is that the effect of electron correlation plays a significant role in determining the polarizability of the clusters and the MP2 method can be considered to be one of the most reliable methods for their prediction. It has also been noted that the polarizability value of the lower member clusters (Na2 and Na4) calculated by highly sophisticated methods such as, CCSD and CCSD(T) are found to be very close to the corresponding MP2 values. The polarizability and the binding energy of the clusters are found to be inversely related to each other and their correlation is rationalized by invoking the minimum polarizability principle. A good linear correlation between the polarizability and volume of the cluster has also been found to exist.     

Comparison of ab initio and DFT electronic structure methods for peptides containing an aromatic ring: effect of dispersion and BSSE

We establish that routine B3LYP and MP2 methods give qualitatively wrong conformations for flexible organic systems containing pi systems and that recently developed methods can overcome the known inadequacies of these methods. This is illustrated for a molecule (a conformer of the Tyr-Gly dipeptide) for which B3LYP/6-31+G(d) and MP2/6-31+G(d) geometry optimizations yield strikingly different structures [Mol. Phys. 2006, 104, 559-570]: MP2 predicts a folded "closed-book" conformer with the glycine residue located above the tyrosine ring, whereas B3LYP predicts a more open conformation. By employing different levels of theory, including the local electron correlation methods LMP2 (local MP2) and LCCSD(T0) (local coupled cluster with single, double, and noniterative local triple excitations) and large basis sets (aug-cc-pVnZ, n=D, T, Q), it is shown that the folded MP2 minimum is an artifact caused by large intramolecular BSSE (basis set superposition error) effects in the MP2/6-31+G(d) calculations. The B3LYP functional gives the correct minimum, but the potential energy apparently rises too steeply when the glycine and tyrosine residues approach each other, presumably due to missing dispersion effects in the B3LYP calculations. The PWB6K and M05-2X functionals, designed to give good results for weak interactions, remedy this to some extent. The reduced BSSE in the LMP2 calculations leads to faster convergence with increasing basis set quality, and accurate results can be obtained with smaller basis sets as compared to canonical MP2. We propose LMP2 as a suitable method to study interactions with pi-electron clouds.     

Accuracy and limitations of second-order many-body perturbation theory for predicting vertical detachment energies of solvated-electron clusters

Vertical electron detachment energies (VDEs) are calculated for a variety of (H(2)O)(n)(-) and (HF)(n)(-) isomers, using different electronic structure methodologies but focusing in particular on a comparison between second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster theory with noniterative triples, CCSD(T). For the surface-bound electrons that characterize small (H(2)O)(n)(-) clusters (n< or = 7), the correlation energy associated with the unpaired electron grows linearly as a function of the VDE but is unrelated to the number of monomers, n. In every example considered here, including strongly-bound "cavity" isomers of (H(2)O)(24)(-), the correlation energy associated with the unpaired electron is significantly smaller than that associated with typical valence electrons. As a result, the error in the MP2 detachment energy, as a fraction of the CCSD(T) value, approaches a limit of about -7% for (H(2)O)(n)(-) clusters with VDEs larger than about 0.4 eV. CCSD(T) detachment energies are bounded from below by MP2 values and from above by VDEs calculated using second-order many-body perturbation theory with molecular orbitals obtained from density functional theory. For a variety of both strongly- and weakly-bound isomers of (H(2)O)(20)(-) and (H(2)O)(24)(-), including both surface states and cavity states, these bounds afford typical error bars of +/-0.1 eV. We have found only one case where the Hartree-Fock and density functional orbitals differ qualitatively; in this case the aforementioned bounds lie 0.4 eV apart, and second-order perturbation theory may not be reliable.     

Accurate calculation of the heats of formation for large main group compounds with spin-component scaled MP2 methods

Three MP2-type electron correlation treatments and standard density functional theory (DFT) approaches are used to predict the heats of formation for a wide variety of different molecules. The SCF and MP2 calculations are performed efficiently using the resolution-of-the-identity (RI) approximation such that large basis set (i.e., polarized valence quadruple-zeta quality) treatments become routinely possible for systems with 50-100 atoms. An atom equivalent scheme that corrects the calculated atomic energies is applied to extract the "real" accuracy of the methods for chemically relevant problems. It is found that the spin-component-scaled MP2 method (SCS-MP2, J. Chem. Phys, 2003, 118, 9095) performs best and provides chemical accuracy (MAD of 1.18 kcal/mol) for a G2/97 test set of molecules. The computationally more economical SOS-MP2 variant, which retains only the opposite-spin part of the correlation energy, is slightly less accurate (MAD of 1.36 kcal/mol) than SCS-MP2. Both spin-component-scaled MP2 treatments perform significantly better than standard MP2 (MAD of 1.77 kcal/mol) and DFT-B3LYP (MAD of 2.12 kcal/mol). These conclusions are supported by results obtained for a second test set of complex systems containing 70 molecules, including charged, strained, polyhalogenated, hypervalent, and large unsaturated species (e.g. C60). For this set, DFT-B3LYP performs badly (MAD of 8.6 kcal/mol) with many errors >10-20 kcal/mol while the spin-component-scaled MP2 methods are still very accurate (MAD of 2.8 and 3.7 kcal/mol, respectively). DFT-B3LYP shows an obvious tendency to underestimate molecular stability as the system size increases. Out of six density functionals tested, the hybrid functional PBE0 performs best. All in all, the SCS-MP2 method, together with large AO basis sets, clearly outperforms current DFT approaches and seems to be the most accurate quantum chemical model that routinely can predict the thermodynamic properties of large main group compounds.     

On the relation between basis set convergence and electron correlation: a critical test for modern ab initio quantum chemistry on a ��mindless�� data set

The correlation of the only two error sources in the solution of the electronic Schrödinger equation is addressed: the basis set convergence (incompleteness) error (BSIE) and the electron correlation effect. The electron correlation effect and basis set incompleteness error are found to be correlated for all of the molecules in Grimme??s ??mindless?? data set (MB08-165). One can use an extrapolation to the HF or MP2 complete basis set (CBS) limit to see with which type of quantum chemical problem (??simple?? and ??hard??) the researcher is dealing. The origin of the slow convergence of the partial wave expansion can be the Kato cusp condition for electron?Celectron coalescence. Such an extrapolation is possible for many large molecular systems and would give the researcher an idea about the expected electron correlation level that would lead to the desired theoretical accuracy. In other words, it is possible to use not only the CBS energy value itself but the speed with which it is reached to get extra information about the molecular system under study.     

<Emphasis Type="Italic">Ab initio</Emphasis> study of small coinage metal telluride clusters Au<Subscript>n</Subscript>Te<Subscript>m</Subscript> (<Emphasis Type="Italic">n,m</Emphasis> = 1, 2)

The geometries of the most stable isomers of gold telluride systems AuTe, Au₂Te, and AuTe₂ are determined using the MP2 method. The aspect of gold—telluride interaction, the electron correlation, and relativistic effects on geometry and stability are investigated at the MP2 and CCSD(T) theoretical levels. The results show that the electron correlation and relativistic effects are responsible not only for gold—gold attraction but also for additional gold—telluride interaction. The gold—telluride interaction is strong enough to modify the known pattern of bare gold clusters. Both effects are essential for determining the geometry and relative stability of this type of systems.     

Structures and properties of stable Al4, Al4+, and Al4- comparatively studied by ab initio theories

Previous high-level theoretical calculations of aluminum clusters mostly relied on density functional theory (DFT) or theories less sophisticated than it. Here, we point out that a second-order M?ller-Plesset perturbation (MP2) method is more appropriate and is the minimum level of theory in the predictions of property such as geometries, electronic structures, and IR and Raman spectra of Al4, Al4+, and Al4- clusters. The theoretical electron affinities and ionization potentials predicted with MP2 geometries are closer to experimental ones than those by DFT. The p-electron characters and single valence state of aluminum atoms and IR and Raman spectra of the aluminum clusters were also reliably predicted by MP2 and could be based on for further experimental and theoretical studies.     

Towards benchmark second-order correlation energies for large atoms: Zn2+ revisited

To provide very accurate reference results for the second-order M?ller-Plesset (MP2) energy and its various components for Zn(2+), which plays for 3d-electron systems a similar role as Ne for smaller atoms and molecules, we have performed extensive calculation by two completely different implementations of the MP2 method: the finite element method (FEM) and the variation-perturbation (VP) method. The FEM and VP calculations yield partial wave contributions up to l(max)=45 and 12, respectively. Detailed comparison of all FEM and VP energy components for l(max)=12 has disclosed an extraordinary similarity, which justifies using the present results as benchmarks. The present correlation energies are compared with other works. The dependability of an earlier version of FEM, already applied to very large closed-shell atoms, is confirmed. It has been found that for larger atoms the accuracy of the analytical Hartree-Fock results has an impact on the accuracy of the MP2 energies greater than for smaller atoms. Fields of applications of the present results in studies of various electron correlation effects in 3d-electron atoms and molecules are indicated.     

Divide-and-conquer local correlation approach to the correlation energy of large molecules

A divide-and-conquer local correlation approach for correlation energy calculations on large molecules is proposed for any post-Hartree-Fock correlation method. The main idea of this approach is to decompose a large system into various fragments capped by their local environments. The total correlation energy of the whole system can be approximately obtained as the summation of correlation energies from all capped fragments, from which correlation energies from all adjacent caps are removed. This approach computationally achieves linear scaling even for medium-sized systems. Our test calculations for a wide range of molecules using the 6-31G or 6-31G( * *) basis set demonstrate that this simple approach recovers more than 99.0% of the conventional second-order Moller-Plesset perturbation theory and coupled cluster with single and double excitations correlation energies.     

Rapid topography mapping of scalar fields: Large molecular clusters

An efficient and rapid algorithm for topography mapping of scalar fields, molecular electron density (MED) and molecular electrostatic potential (MESP) is presented. The highlight of the work is the use of fast function evaluation by Deformed-atoms-in-molecules (DAM) method. The DAM method provides very rapid as well as sufficiently accurate function and gradient evaluation. For mapping the topography of large systems, the molecular tailoring approach (MTA) is invoked. This new code is tested out for mapping the MED and MESP critical points (CP's) of small systems. It is further applied to large molecular clusters viz. (H(2)O)(25), (C(6)H(6))(8) and also to a unit cell of valine crystal at MP2∕6-31+G(d) level of theory. The completeness of the topography is checked by extensive search as well as applying the Poincare?-Hopf relation. The results obtained show that the DAM method in combination with MTA provides a rapid and efficient route for mapping the topography of large molecular systems.     

Visualizing internal stabilization in weakly bound systems using atomic energies: hydrogen bonding in small water clusters

Atomic energies are used to visualize the local stabilizing and destabilizing energy changes in water clusters. Small clusters, (H(2)O)(n), from n = 2-5, at MP2/aug-cc-pVTZ geometries are evaluated using energies defined by the quantum theory of atoms in molecules (QTAIM). The atomic energies reproduce MP2 total energies to within 0.005 kcal mol(-1). Oxygen atoms are stabilized for all systems and hydrogen atoms are destabilized. The increased stability of the water clusters due to hydrogen bond cooperativity is demonstrated at an atomic level. Variations in atomic energies within the clusters are correlated to the geometry of the waters and reveal variations in the hydrogen bond strengths. The method of visualization of the energy changes applied here is especially suited for application to large biomolecules.     

Strategies for electron correlation calculations on large molecular systems

The two main problems in applying common methods for electron correlation to large systems are (1) the steep dependency of the computational task with the size of the system and (2) the large and often prohibitive number of two-electron integrals that must be stored and retrieved during the calculation. The direct local correlation method eliminates both these problems. This method is a combination of the local correlation method that eliminates the steep computational dependency with the size of the system and extension of the direct SCF approach to electron correlation methods. The latter method eliminates external storage of the electron repulsion integrals. In this review, the direct local correlation approach and its computer implementation will be described in detail, including computational examples and comparisons to similar methods developed in other groups. The feasibility of future applications to extended systems is discussed.     

Analysis of fourth-order Møller–Plesset limit energies: the importance of three-electron correlation effects

Fourth-order M?ller–Plesset (MP4) correlation energies are computed for 28 atoms and simple molecules employing Dunning's correlation-consistent polarized-valence m-zeta basis sets for m=2, 3, 4, and 5. Extrapolation formulas are used to predict MP4 energies for infinitely large basis sets. It is shown that both total and partial MP4 correlation energies can be extrapolated to limit values and that the sum of extrapolated partial MP4 energies equals the extrapolated total MP4 correlation energy within calculational accuracy. Therefore, partial MP4 correlation energies can be presented in the form of an MP4 spectrum reflecting the relative importance of different correlation effects. Typical trends in calculated correlation effects for a given class of electron systems are independent of the basis set used. As first found by Cremer and He [(1996) J Phys Chem 100:6173], one can use MP4 spectra to distinguish between electron systems with well-separated electron pairs and systems for which electrons cluster in a confined region of atomic or molecular space. MP4 spectra for increasing size of the basis set reveal that smaller basis set calculations underestimate the importance of three-electron correlation effects for both classes by overestimating the importance of pair correlation effects. The minimum size of a basis set required for reliable MP4 calculations is given by a valence triple-zeta polarized basis, which even in the case of anions performs better than a valence double-zeta basis augmented by diffuse functions. Received: 14 June 2000 / Accepted: 16 June 2000 / Published online: 24 October 2000