SCF and localized orbitals in ethylene: MBPT/CC results and comparison with one-million configuration Cl

Many-body perturbation theory (MBPT) and coupled-cluster (CC) calculations are performed on the ethylene molecule employing canonical SCF and simple bond-orbital localized orbitals (LO). Full fourth-order MBPT [i.e. SDTQ MBPT(4)], CC doubles (CCD) and CC singles and doubles (CCSD) energies are compared with the over one-million configuration ‘bench-mark” Cl calculation of Saxe et al. Though the SCF and LO reference determinant energies differ by 0.29706 hartree, the CCSD energy difference is only 1.7 mhartrees (mh). Our most extensive SCF orbital calculation, CCSD plus fourth-order triples, is found to be lower in energy than the CI result by 5.3 mh.     

Comparison of MBPT and coupled-cluster methods with full CI. Importance of triplet excitation and infinite summations

Results from full fourth-order perturbation theory [SDTQ MBPT(4)], and the coupled-cluster single- and double-excitation model (CCSD). are compared with recent full CI results for BH, HF, NH₃, and H₂O. For H₂O, studies include large symmetric displacements of the OH bonds, which offer a severe test for any single-reference MBPT/CC method. In every case. CCSD plus fourth-order triple-excitation terms provide agreement with the full CI to < 2 kcal/mole. SDTQ MBPT(4) has an error 10 kcal/mole for displaced H₂O.     

Explicit correlation and basis set superposition error: the structure and energy of carbon dioxide dimer

We have investigated the slipped parallel and t-shaped structures of carbon dioxide dimer [(CO(2))(2)] using both conventional and explicitly correlated coupled cluster methods, inclusive and exclusive of counterpoise (CP) correction. We have determined the geometry of both structures with conventional coupled cluster singles doubles and perturbative triples theory [CCSD(T)] and explicitly correlated cluster singles doubles and perturbative triples theory [CCSD(T)-F12b] at the complete basis set (CBS) limits using custom optimization routines. Consistent with previous investigations, we find that the slipped parallel structure corresponds to the global minimum and is 1.09 kJ mol(-1) lower in energy. For a given cardinal number, the optimized geometries and interaction energies of (CO(2))(2) obtained with the explicitly correlated CCSD(T)-F12b method are closer to the CBS limit than the corresponding conventional CCSD(T) results. Furthermore, the magnitude of basis set superposition error (BSSE) in the CCSD(T)-F12b optimized geometries and interaction energies is appreciably smaller than the magnitude of BSSE in the conventional CCSD(T) results. We decompose the CCSD(T) and CCSD(T)-F12b interaction energies into the constituent HF or HF CABS, CCSD or CCSD-F12b, and (T) contributions. We find that the complementary auxiliary basis set (CABS) singles correction and the F12b approximation significantly reduce the magnitude of BSSE at the HF and CCSD levels of theory, respectively. For a given cardinal number, we find that non-CP corrected, unscaled triples CCSD(T)-F12b/VXZ-F12 interaction energies are in overall best agreement with the CBS limit.     

Computational chemistry with transputers: a direct SCF program

By using transputers it is possible to build up networks of parallel processors with varying topology. Due to the architecture of the processors it is appropriate to use the MIMD (multiple instruction multiple data) concept of parallel computing. The most suitable programming language is OCCAM. We investigate the use of transputer networks in computational chemistry, starting with the direct SCF method. The most time consuming step, the calculation of the two electron integrals is executed parallelly. Each node in the network calculates whole batches of integrals. The main program is written in OCCAM. For some large-scale arithmetic processes running on a single node, however, we used FORTRAN subroutines out of standard ab-initio programs to reduce the programming effort. Test calculations show, that the integral calculation step can be parallelized very efficiently. We observe a speed-up of almost 8 using eight network processors. Even in consideration of the scalar part of the SCF iteration, the speed-up is not less than 7.1.     

Distributed data parallel coupled-cluster algorithm: Application to the 2-hydroxypyridine/2-pyridone tautomerism

The recently developed parallel coupled-cluster algorithm of Rendell, Lee, and Lindh [Chem. Phys. Lett., 194 , 84 (1992)] is extended to allow four-indexed quantities containing one or two indices in the virtual orbital space to be stored across the global memory of distributed-memory parallel processors. Quantities such as the double-excitation amplitudes can now be distributed over multiple nodes, with blocks of data retrieved from remote nodes by the use of interrupt handlers. As an application of the new code, we have investigated the potential energy surface of the 2-hydroxypyridine/2-pyridone tautomers. Using large basis sets, the structure of each tautomer and the transition state connecting the two minima has been determined at the SCF level. The relative energy difference and the activation energy were then redetermined using the MP2, CCSD, and CCSD(T) methods. All calculations have been performed on Intel distributed-memory supercomputers. The largest coupled-cluster calculations contained over 2 million double-excitation amplitudes. © John Wiley & Sons, Inc.     

Symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles method: improving upon CCSD(T) and CCSD(T)(Λ): preliminary application

Symmetric and asymmetric triple excitation corrections for the orbital-optimized coupled-cluster doubles (OO-CCD or simply "OD" for short) method are investigated. The conventional symmetric and asymmetric perturbative triples corrections [(T) and (T)(Λ)] are implemented, the latter one for the first time. Additionally, two new triples corrections, denoted as OD(Λ) and OD(Λ)(T), are introduced. We applied the new methods to potential energy surfaces of the BH, HF, C(2), N(2), and CH(4) molecules, and compare the errors in total energies, with respect to full configuration interaction, with those from the standard coupled-cluster singles and doubles (CCSD), with perturbative triples [CCSD(T)], and asymmetric triples correction (CCSD(T)(Λ)) methods. The CCSD(T) method fails badly at stretched geometries, the corresponding nonparallelity error is 7-281 kcal mol(-1), although it gives reliable results near equilibrium geometries. The new symmetric triples correction, CCSD(Λ), noticeably improves upon CCSD(T) (by 4-14 kcal mol(-1)) for BH, HF, and CH(4); however, its performance is worse than CCSD(T) (by 1.6-4.2 kcal mol(-1)) for C(2) and N(2). The asymmetric triples corrections, CCSD(T)(Λ) and CCSD(Λ)(T), perform remarkably better than CCSD(T) (by 5-18 kcal mol(-1)) for the BH, HF, and CH(4) molecules, while for C(2) and N(2) their results are similar to those of CCSD(T). Although the performance of CCSD and OD is similar, the situation is significantly different in the case of triples corrections, especially at stretched geometries. The OD(T) method improves upon CCSD(T) by 1-279 kcal mol(-1). The new symmetric triples correction, OD(Λ), enhances the OD(T) results (by 0.01-2.0 kcal mol(-1)) for BH, HF, and CH(4); however, its performance is worse than OD(T) (by 1.9-2.3 kcal mol(-1)) for C(2) and N(2). The asymmetric triples corrections, OD(T)(Λ) and OD(Λ)(T), perform better than OD(T) (by 2.0-6.2 kcal mol(-1)). The latter method is slightly better for the BH, HF, and CH(4) molecules. However, for C(2) and N(2) the new results are similar to those of OD(T). For the BH, HF, and CH(4) molecules, OD(Λ)(T) provides the best potential energy curves among the considered methods, while for C(2) and N(2) the OD(T) method prevails. Hence, for single-bond breaking the OD(Λ)(T) method appears to be superior, whereas for multiple-bond breaking the OD(T) method is better.     

The coupled cluster singles, doubles, and a hybrid treatment of connected triples based on the split virtual orbitals

We have proposed a simple strategy for splitting the virtual orbitals with a large basis set into two subgroups (active and inactive) by taking a smaller basis set as an auxiliary basis set. With the split virtual orbitals (SVOs), triple or higher excitations can be partitioned into active and inactive subgroups (according to the number of active virtual orbitals involved), which can be treated with different electron correlation methods. In this work, the coupled cluster (CC) singles, doubles, and a hybrid treatment of connected triples based on the SVO [denoted as SVO-CCSD(T)-h], has been implemented. The present approach has been applied to study the bond breaking potential energy surfaces in three molecules (HF, F(2), and N(2)), and the equilibrium properties in a number of open-shell diatomic molecules. For all systems under study, the SVO-CCSD(T)-h method based on the unrestricted Hartree-Fock (UHF) reference is an excellent approximation to the corresponding CCSDT (CC singles, doubles, and triples), and much better than the UHF-based CCSD(T) (CC singles, doubles, and perturbative triples). On the other hand, the SVO-CCSD(T)-h method based on the restricted HF (RHF) reference can also provide considerable improvement over the RHF-based CCSD(T).     

A gauge-independent zeroth-order regular approximation to the exact relativistic Hamiltonian--formulation and applications

A simple modification of the zeroth-order regular approximation (ZORA) in relativistic theory is suggested to suppress its erroneous gauge dependence to a high level of approximation. The method, coined gauge-independent ZORA (ZORA-GI), can be easily installed in any existing nonrelativistic quantum chemical package by programming simple one-electron matrix elements for the quasirelativistic Hamiltonian. Results of benchmark calculations obtained with ZORA-GI at the Hartree-Fock (HF) and second-order Moller-Plesset perturbation theory (MP2) level for dihalogens X(2) (X=F,Cl,Br,I,At) are in good agreement with the results of four-component relativistic calculations (HF level) and experimental data (MP2 level). ZORA-GI calculations based on MP2 or coupled-cluster theory with single and double perturbations and a perturbative inclusion of triple excitations [CCSD(T)] lead to accurate atomization energies and molecular geometries for the tetroxides of group VIII elements. With ZORA-GI/CCSD(T), an improved estimate for the atomization energy of hassium (Z=108) tetroxide is obtained.     

Is fifth-order MBPT enough?

Fifth- and higher-order MBPT results are reported for a series of examples, BH, Be₂, HF and H₂O, for which higher-order perturbation theory might be important. MBPT(5) differs from MBPT(4) by as much as 4.3 mh, and by constructing the size-extensive [2,1] Padé approximant, which is possible with E⁽⁵⁾, one can get exceptional agreement with the full CI results. Variational perturbation results are also reported.     

The coupled cluster approach with a hybrid treatment of connected triple excitations based on the restricted Hartree-Fock reference

A generalization of the coupled cluster (CC) singles, doubles, and a hybrid treatment of connected triples [denoted as CCSD(T)-h] [Shen et al., J. Chem. Phys. 132, 114115 (2010)] to the restricted Hartree-Fock (RHF) reference is presented. In this approach, active (or pseudoactive) RHF orbitals are constructed automatically by performing unitary transformations of canonical RHF orbitals so that they spatially mimic the natural orbitals of the unrestricted Hartree-Fock reference. The present RHF-based CCSD(T)-h approach has been applied to study the potential energy surfaces in several typical bond breaking processes and the singlet-triplet gaps in a diradical (HFH)(-1). For all systems under study, the overall performance of CCSD(T)-h is very close to that of the corresponding CCSD(T) (CC singles, doubles, and triples), and much better than that of CCSD(T) (CC singles, doubles, and perturbative triples).     

Ab initio study of Rg2I− (Rg = Ar,Kr, Xe)

The equilibrium geometries, vibrational frequencies, and dissociation energies of rare gas iodine clusters Rg₂I^?(Rg = Ar, Kr, Xe) were calculated at the Hartree–Fock (HF), second‐order Møller–Plesset (MP2), the coupled cluster method with single and double excitation and a noniterative correction for triple excitations method [CCSD(T)] levels. The title species have bent C_2v structure of about 60° angle. The electron correlation effects and relativistic effects on the geometry and stability were investigated at CCSD(T) level. Both effects stabilize title species. The calculated electron affinities are in good agreement with the experimental results available. The effect of high angular momentum functions (g and h) was studied. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Electron correlation effects on the angular momentum anisotropies of the dipole polarizabilities of the first-row stable atomic anions

The dipole polarizabilities of B⁻(³P) and O⁻(²P) are calculated using large basis sets and inclusion of high-order electron correlation effects. Special attention is given to the anisotropies derived from the different components of the orbital angular momentum |M_L| = 1, 0. A very large influence of the results on electron correlation is found. This influence is larger for the |M_L| = 1 component. However, whereas for B⁻ no low-order level gives acceptable results, for O⁻ due to cancellation effects, MBPT (4), and even MBPT (2), approximate the highest-level [CCSD(T)] result. Our best results are 463.1 and 400.3 au for the average and anisotropic polarizabilities of B⁻. For O⁻, these numbers are 37.7, and 22.2 au, respectively. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 459–463, 1997     

Near-degeneracy effects and efficiency of MBPT calculations of molecular properties. Polarizability curve of H2

The polarizability curve of H₂ is calculated by using the finite-field perturbation method. All self-consistency effects are accounted for at the HF level and many-body perturbation theory (MBPT) is used to evaluate the correlation contributions. Using a single HF determinant as a reference in MBPT calculations makes the near-degeneracy effects of essential importance on increasing the interatomic distance. Nevertheless, applying the MBPT scheme with appropriate fourth-order terms gives nearly exact values of both components of the polarizability tensor for interatomic distances up to ≈3.6 au.     

Reduced multireference coupled-cluster method: barrier heights for heavy atom transfer, nucleophilic substitution, association, and unimolecular reactions

The recently developed reduced multireference coupled-cluster method with singles and doubles (RMR CCSD) that is perturtatively corrected for triples [RMR CCSD(T)] is employed to compute the forward and reverse barrier heights for 19 non-hydrogen-transfer reactions. The method represents an extension of the conventional single-reference (SR) CCSD(T) method to multireference situations. The results are compared with a benchmark database, which is essentially based on the SR CCSD(T) results. With the exception of seven cases, the RMR CCSD(T) results are almost identical with those based on SR CCSD(T), implying the abatement of MR effects at the SD(T) level relative to the SD level. Using the differences between the RMR CCSD(T) and CCSD(T) barrier heights as a measure of MR effects, modified values for barrier heights of studied reactions are given.     

Interaction of dihydrogen with small and light molecules

Second-order M?ller-Plesset (MP2) calculations (using the approximate resolution of the identity, RI-MP2), explicitly correlated MP2 (MP2-R12) calculations, and coupled-cluster calculations including all single and double excitations with a perturbative estimate of triple excitations [CCSD(T)] are performed to study the interaction of molecular hydrogen with the small molecules HF, H2O, NH3, and LiOH. Different adsorption positions are studied. In the cases of H2O and NH3, the most favorable configuration places H2 in an end-on fashion on the O or N atom, respectively. In the cases of HF and LiOH, the H2 molecule takes a side-on position on the H atom of HF or the Li atom. With respect to MP2 calculations in a triple-zeta basis, both the enlargement of the basis set and the extension of the correlation treatment (CCSD(T) vs MP2) increase the interaction energy. The basis set limit CCSD(T) estimates of the interaction energy of H2 with the HF, H2O, NH3, and LiOH molecules amount to 4.40, 2.67, 3.02, and 10.74 kJ mol-1, respectively. The interaction energy for the simultaneous interaction of H2 with two LiOH molecules does not significantly exceed the value obtained for the interaction with a single LiOH molecule. Furthermore, the interaction energies (by MP2) of H2 with glycine, the glycine dimer, and imidazolium chloride amount to 2.78, 5.00, and 6.30 kJ mol-1, respectively.     

Quantum supercharger library: Hyper‐parallelism of the Hartree–Fock method

We present here a set of algorithms that completely rewrites the Hartree–Fock (HF) computations common to many legacy electronic structure packages (such as GAMESS‐US, GAMESS‐UK, and NWChem) into a massively parallel compute scheme that takes advantage of hardware accelerators such as Graphical Processing Units (GPUs). The HF compute algorithm is core to a library of routines that we name the Quantum Supercharger Library (QSL). We briefly evaluate the QSL's performance and report that it accelerates a HF 6‐31G Self‐Consistent Field (SCF) computation by up to 20 times for medium sized molecules (such as a buckyball) when compared with mature Central Processing Unit algorithms available in the legacy codes in regular use by researchers. It achieves this acceleration by massive parallelization of the one‐ and two‐electron integrals and optimization of the SCF and Direct Inversion in the Iterative Subspace routines through the use of GPU linear algebra libraries. © 2015 Wiley Periodicals, Inc.     

Peptide models. XXXIII. Extrapolation of low-level Hartree-Fock data of peptide conformation to large basis set SCF,MP2, DFT,and CCSD(T) results. The Ramachandran surface of alanine dipeptide computed at various levels of theory

At the dawn of the new millenium, new concepts are required for a more profound understanding of protein structures. Together with NMR and X-ray-based 3D-structure determinations in silico methods are now widely accepted. Homology-based modeling studies, molecular dynamics methods, and quantum mechanical approaches are more commonly used. Despite the steady and exponential increase in computational power, high level ab initio methods will not be in common use for studying the structure and dynamics of large peptides and proteins in the near future. We are presenting here a novel approach, in which low- and medium-level ab initio energy results are scaled, thus extrapolating to a higher level of information. This scaling is of special significance, because we observed previously on molecular properties such as energy, chemical shielding data, etc., determined at a higher theoretical level, do correlate better with experimental data, than those originating from lower theoretical treatments. The Ramachandran surface of an alanine dipeptide now determined at six different levels of theory [RHF and B3LYP 3-21G, 6-31+G(d) and 6-311++G(d,p)] serves as a suitable test. Minima, first-order critical points and partially optimized structures, determined at different levels of theory (SCF, DFT), were completed with high level energy calculations such as MP2, MP4D, and CCSD(T). For the first time three different CCSD(T) sets of energies were determined for all stable B3LYP/6-311++G(d,p) minima of an alanine dipeptide. From the simplest ab initio data (e.g., RHF/3-21G) to more complex results [CCSD(T)/6-311+G(d,p)//B3LYP/6-311++G(d,p)] all data sets were compared, analyzed in a comprehensive manner, and evaluated by means of statistics.     

Benchmark calculations on the adiabatic ionization potentials of M-NH(3) (M=Na,Al,Ga,In,Cu,Ag)

The ground states of the M-NH(3) (M=Na,Al,Ga,In,Cu,Ag) complexes and their cations have been studied with density functional theory and coupled cluster [CCSD(T)] methods. The adiabatic ionization potentials (AIPs) of these complexes are calculated, and these are compared to results from high-resolution zero-electron kinetic energy photoelectron spectroscopy. By extrapolating the CCSD(T) energies to the complete basis set (CBS) limit and including the core-valence, scalar relativistic, spin-orbit, and zero-point corrections, the CCSD(T) method is shown to be able to predict the AIPs of these complexes to better than 6 meV or 0.15 kcal/mol. 27 exchange-correlation functionals, including one in the local density approximation, 13 in the generalized gradient approximation (GGA), and 13 with hybrid GGAs, were benchmarked in the calculations of the AIPs. The B1B95, mPW1PW91, B98, B97-1, PBE1PBE, O3LYP, TPSSh, and HCTH93 functionals give an average error of 0.1 eV for all the complexes studied, with the B98 functional alone yielding a maximum error of 0.1 eV. In addition, the calculated metal-ammonia harmonic stretching frequencies with the CCSD(T) method are in excellent agreement with their experimental values, whereas the B3LYP method tends to underestimate these stretching frequencies. The metal-ammonia binding energies were also calculated at the CCSD(T)/CBS level, and are in excellent agreement with the available experimental values considering the error limits, except for Ag-NH(3) and Ag(+)-NH(3), where the calculations predict stronger bond energies than measured by about 4 kcal/mol, just outside the experimental error bars of +/-3 kcal/mol.     

13C carbonyl chemical shielding tensors: Comparing SCF,MBPT (2), and DFT predictions to experiment

In this work, we calculate the ¹³C nuclear magnetic resonance chemical shielding tensors for 18 carbonyl-containing compounds. The many-body perturbation theory (MBPT), self-consistent field (SCF), and density functional theory (DFT) formalisms were used with gauge including atomic orbitals (GIAO) to calculate the shielding tensors. Our data suggest that shielding tensors can be efficiently estimated by performing one MBPT(2) correlated calculation (e.g., at a reference geometry) and SCF-level calculations at other geometries and taking the SCF-to-correlated tensor element differences to be geometry independent. That is, the correlation contribution to the chemical shielding seems to be relatively constant over a considerable range of distortions. Treatment of correlation using DFT methods is shown to not be as systematically reliable as with MBPT(2). Data on 18 carbonyl compounds show that the single largest influence on the shielding tensor is the presence of nearby electron-withdrawing or electron-donating groups. Finally, although good agreement with powder or single-crystal experimental data is achieved for two or three tensor eigenvalues, systematic differences remain for one element; the origins of these differences are discussed. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 875–894, 1997     

Coupled-cluster singles,doubles and perturbative triples with density fitting approximation for massively parallel heterogeneous platforms

A high-performance implementation of the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] is developed in the Massively Parallel Quantum Chemistry program. Novel features include: (1) reduced memory requirements via a density-fitting (DF) CCSD implementation utilizing distributed lazy evaluation for tensors with more than two unoccupied indices and (2) the ability to utilize efficiently many-core nodes (Intel Xeon Phi) and heterogeneous nodes with multiple NVIDIA GPUs on each node. All data that are greater than quadratic in the system size are distributed among processes. Excellent strong scaling is observed on distributed-memory computers equipped with conventional CPUs, Intel Xeon Phi processors, and heterogeneous nodes with multiple NVIDIA GPUs Canonical CCSD(T) energies can be evaluated for systems containing 200 electrons and 1000 basis functions in a few days using a small size commodity cluster, with even larger computations possible on leadership-class computing resources.