On the equivalence of ring-coupled cluster and adiabatic connection fluctuation-dissipation theorem random phase approximation correlation energy expressions

The correlation energy in the direct random phase approximation (dRPA) can be written, among other possibilities, either in terms of the interaction strength averaged correlation density matrix, or in terms of the coupled cluster doubles amplitudes obtained in the direct ring approximation (drCCD). Although the corresponding dRPA correlation density matrix on the one hand, and the drCCD amplitude matrix on the other hand, differ significantly, they yield identical energies. Similarly, the analogous RPA and rCCD correlation energies calculated from antisymmetrized two-electron integrals are identical to each other despite very different underlying working equations. In the present communication, a direct correspondence between amplitudes and densities is established and investigated with perturbation theory arguments. Our analysis also sheds some light on the properties of recently proposed RPA/rCCD variants which use antisymmetrized integrals in part of the equations and nonantisymmetrized integrals in others.     

Increasing the applicability of density functional theory. II. Correlation potentials from the random phase approximation and beyond

Density functional theory (DFT) results are mistrusted at times due to the presence of an unknown exchange correlation functional, with no practical way to guarantee convergence to the right answer. The use of a known exchange correlation functional based on wave-function theory helps to alleviate such mistrust. The exchange correlation functionals can be written exactly in terms of the density-density response function using the adiabatic-connection and fluctuation-dissipation framework. The random phase approximation (RPA) is the simplest approximation for the density-density response function. Since the correlation functional obtained from RPA is equivalent to the direct ring coupled cluster doubles (ring-CCD) correlation functional, meaning only Coulomb interactions are included, one can bracket RPA between many body perturbation theory (MBPT)-2 and CCD with the latter having all ring, ladder, and exchange contributions. Using an optimized effective potential strategy, we obtain correlation potentials corresponding to MBPT-2, RPA (ring-CCD), linear-CCD, and CCD. Using the suitable choice of the unperturbed Hamiltonian, Kohn-Sham self-consistent calculations are performed. The spatial behavior of the resulting potentials, total energies, and the HOMO eigenvalues are compared with the exact values for spherical atoms. Further, we demonstrate that the self-consistent eigenvalues obtained from these consistent potentials used in ab initio dft approximate all principal ionization potentials as demanded by ionization potential theorem.     

Third-order corrections to random-phase approximation correlation energies

Several random-phase approximation (RPA) correlation methods were compared in third order of perturbation theory. While all of the considered approaches are exact in second order of perturbation theory, it is found that their corresponding third-order correlation energy contributions strongly differ from the exact third-order correlation energy contribution due to missing interactions of the particle-particle-hole-hole type. Thus a simple correction method is derived which makes the different RPA methods also exact to third-order of perturbation theory. By studying the reaction energies of 16 chemical reactions for 21 small organic molecules and intermolecular interaction energies of 23 intermolecular complexes comprising weakly bound and hydrogen-bridged systems, it is found that the third-order correlation energy correction considerably improves the accuracy of RPA methods if compared to coupled-cluster singles doubles with perturbative triples as a reference.     

Linear-response theory for Mukherjee's multireference coupled-cluster method: Excitation energies

The recently presented linear-response function for Mukherjee's multireference coupled-cluster method (Mk-MRCC) [T.-C. Jagau and J. Gauss, J. Chem. Phys. 137, 044115 (2012)] is employed to determine vertical excitation energies within the singles and doubles approximation (Mk-MRCCSD-LR) for ozone as well as for o-benzyne, m-benzyne, and p-benzyne, which display increasing multireference character in their ground states. In order to assess the impact of a multireference ground-state wavefunction on excitation energies, we compare all our results to those obtained at the single-reference coupled-cluster level of theory within the singles and doubles as well as within the singles, doubles, and triples approximation. Special attention is paid to the artificial splitting of certain excited states which arises from the redundancy intrinsic to Mk-MRCC theory and hinders the straightforward application of the Mk-MRCC-LR method.     

Linear-response theory for Mukherjee's multireference coupled-cluster method: Static and dynamic polarizabilities

The formalism of response theory is applied to derive expressions for static and dynamic polarizabilities within the state-specific multireference coupled-cluster theory suggested by Mukherjee and co-workers (Mk-MRCC) [J. Chem. Phys. 110, 6171 (1998)]. We show that the redundancy problem inherent to Mk-MRCC theory gives rise to spurious poles in the Mk-MRCC response functions, which hampers the reliable calculation of dynamic polarizabilities. Furthermore, we demonstrate that in the case of a symmetry-breaking perturbation a working response theory is obtained only if certain internal excitations are included in the responses of the cluster amplitudes. Exemplary calculations within the singles and doubles approximation (Mk-MRCCSD) are carried out on aryne compounds to illustrate the impact of a multireference ansatz on the polarizability.     

Tensor factorizations of local second-order Møller-Plesset theory

Efficient electronic structure methods can be built around efficient tensor representations of the wavefunction. Here we first describe a general view of tensor factorization for the compact representation of electronic wavefunctions. Next, we use this language to construct a low-complexity representation of the doubles amplitudes in local second-order M?ller-Plesset perturbation theory. We introduce two approximations--the direct orbital-specific virtual approximation and the full orbital-specific virtual approximation. In these approximations, each occupied orbital is associated with a small set of correlating virtual orbitals. Conceptually, the representation lies between the projected atomic orbital representation in Pulay-Saeb? local correlation theories and pair natural orbital correlation theories. We have tested the orbital-specific virtual approximations on a variety of systems and properties including total energies, reaction energies, and potential energy curves. Compared to the Pulay-Saeb? ansatz, we find that these approximations exhibit favorable accuracy and computational times while yielding smooth potential energy curves.     

Coupled-cluster methods including noniterative corrections for quadruple excitations

A new method is presented for treating the effects of quadruple excitations in coupled-cluster theory. In the approach, quadruple excitation contributions are computed from a formula based on a non-Hermitian perturbation theory analogous to that used previously to justify the usual noniterative triples correction used in the coupled cluster singles and doubles method with a perturbative treatment of the triple excitations (CCSD(T)). The method discussed in this paper plays a parallel role in improving energies obtained with the full coupled-cluster singles, doubles, and triples method (CCSDT) by adding a perturbative treatment of the quadruple excitations (CCSDT(Q)). The method is tested for an extensive set of examples, and is shown to provide total energies that compare favorably with those obtained with the full singles, doubles, triples, and quadruples (CCSDTQ) method.     

Explicitly correlated coupled-cluster theory using cusp conditions. I. Perturbation analysis of coupled-cluster singles and doubles (CCSD-F12)

Geminal functions based on Slater-type correlation factors and fixed expansion coefficients, determined by cusp conditions, have in recent years been forwarded as an efficient and numerically stable method for introducing explicit electron correlation into coupled-cluster theory. In this work, we analyze the equations of explicitly correlated coupled-cluster singles and doubles (CCSD-F12) theory and introduce an ordering scheme based on perturbation theory which can be used to characterize and understand the various approximations found in the literature. Numerical results for a test set of 29 molecules support our analysis and give additional insight. In particular, our results help rationalize the success of the CCSD(F12) approximation which is based on a very systematic cancellation of the neglected, otherwise individually large third-order geminal-geminal coupling terms. Further approximations to CCSD(F12) can be introduced without sacrificing the accuracy if the entire set of third-order coupling terms between the conventional doubles cluster amplitudes and the geminal doubles amplitudes is retained, leading to the recently proposed CCSD[F12] and CCSD(F12(?)) models, which have negligible overhead compared to conventional CCSD calculations. Particularly, the importance of the ring-term type contribution is pointed out which may be used to improve on other existing approximations such as CCSD-F12b. For small basis sets, it might be advantageous to keep certain higher-order terms leading to CCSD-F12(?), which, for the case of the SP ansatz, merely involves a noniterative correction to CCSD(F12(?)).     

Are insertion compounds of CH2CHF and the rare gases stable? A computational study

Ab initio calculations, using second-order M?ller-Plesset perturbation theory with a 6-311+ +G(2d,2p) basis set, predict the stability of two novel compounds of monofluoroethene, CH2CHF, with the rare-gas atoms Ar and Kr. The dissociation energies to the lowest-energy fragmentation products, CH2CHF + Rg (Rg = Ar,Kr), were computed to be -528 and -449 kJ mol(-1), respectively, at the coupled-cluster singles, doubles, and triples level of theory. Possible transition states (at second-order M?ller-Plesset theory) via a C-Rg-F bending mode for these fragmentation reactions were also located with barrier heights of about 76 and 106 kJ mol(-1), for the Ar- and Kr-containing species, respectively. However, the Ar-containing species may not exist at all as it is less stable than the fragments CH2CH + F + Ar at the higher level of theory and may possibly dissociate via this route.     

Coupled-cluster calculation of dispersion contributions to interaction energies and polarizabilities

 The induced dipole dispersion-type contributions to the two-body and nonadditive three-body energies and electric dipole polarizabilities are studied for long-range interactions involving the He, Ne, Ar and Kr atoms and the H₂ and N₂ molecules. The coupled-cluster singles and doubles model and large basis sets are used. Comparison of the energy contributions with data derived from experiment shows in most cases the deviations to be less than 1%; therefore, it may be expected that the calculated polarizability increments are accurately determined and can be used to estimate the accuracy of approximate methods. Received: 20 March 2001 / Accepted: 5 April 2001 / Published online: 27 June 2001     

Extrapolating to the one-electron basis-set limit in electronic structure calculations

A simple, yet reliable, scheme based on treating uniformly singlet-pair and triplet-pair interactions is suggested to extrapolate atomic and molecular electron correlation energies calculated at two basis-set levels of ab initio theory to the infinite one-electron basis-set limit. The novel dual-level method is first tested on extrapolating the full correlation in single-reference coupled-cluster singles and doubles energies for the closed-shell systems CH2((1)A1), H2O, HF, N2, CO, Ne, and F2 with correlation-consistent basis sets of the type cc-pVXZ (X=D,T,Q,5,6) reported by Klopper [Mol. Phys. 6, 481 (2001)] against his own benchmark calculations with large uncontracted basis sets obtained from explicit correlated singles and doubles coupled-cluster theory. Comparisons are also reported for the same data set but using both single-reference Moller-Plesset and coupled-cluster doubles methods. The results show a similar, often better, accordance with the target results than Klopper's extrapolations where singlet-pair and triplet-pair energies are extrapolated separately using the popular X(-3) and X(-5) dual-level laws, respectively. Applications to the extrapolation of the dynamical correlation in multireference configuration interaction calculations carried out anew for He, H2, HeH+, He2 ++, H3+(1 (1)A'), H3+(1 (3)A'), BH, CH, NH, OH, FH, B2, C2, N2, O2, F2, BO, CO, NO, BN, CN, SH, H2O, and NH3 with standard augmented correlation-consistent basis sets of the type aug-cc-pVXZ (X=D,T,Q,5,6) are also reported. Despite lacking accurate theoretical or experimental data for comparison in the case of most diatomic systems, the new method also shows in this case a good performance when judged from the results obtained with the traditional schemes which extrapolate using the two largest affordable basis sets. For the Hartree-Fock and complete-active space self-consistent field energies, a simple pragmatic extrapolation rule is examined whose results are shown to compare well with the ones obtained from the best reported schemes.     

Accurate all-electron correlation energies for the closed-shell atoms from Ar to Rn and their relationship to the corresponding MP2 correlation energies

All-electron correlation energies E(c) are not very well-known for atoms with more than 18 electrons. Hence, coupled-cluster calculations in carefully designed basis sets are combined with fully converged second-order M?ller-Plesset perturbation theory (MP2) computations to obtain fairly accurate, nonrelativistic E(c) values for the 12 closed-shell atoms from Ar to Rn. These energies will be useful for the evaluation and parameterization of density functionals. The results show that MP2 overestimates ∣E(c)∣ for heavy atoms. Spin-component scaling of the MP2 correlation energy is used to provide a simple explanation for this overestimation.     

Improved supermolecular second order Møller-Plesset intermolecular interaction energies using time-dependent density functional response theory

The supermolecular second order Moller-Plesset (MP2) intermolecular interaction energy is corrected by employing time-dependent density functional (TDDFT) response theory. This is done by replacing the uncoupled second order dispersion contribution contained in the supermolecular MP2 energy with the coupled dispersion energy obtained from the TDDFT approach. Preliminary results for the rare gas dimers He2, Ne2, and Ar2 and a few structures of the (HF)2 and (H2O)2 dimers show that the conventional MP2 interaction energies are considerably improved by this procedure if compared to coupled cluster singles doubles with perturbative triples [CCSD(T)] interaction energies. However, the quality of the interaction energies obtained in this way strongly depends on the exchange-correlation potential employed in the monomer calculations: It is shown that an exact exchange-only potential surprisingly often performs better than an asymptotically corrected hybrid exchange-correlation potential. Therefore the method proposed in this work is similar to the method by Cybulski and Lytle [J. Chem. Phys., 127, 141102 (2007)] which corrects the supermolecular MP2 energies with a scaled dispersion energy from time-dependent Hartree-Fock. The results in this work are also compared to the combination of density functional theory and intermolecular perturbation theory.     

Nuclear spin-spin coupling constants in a coupled-cluster polarization propagator approach

A recently derived coupled-cluster doubles polarization propagator approximation (CCDPPA) for excitation energies and transition moments is extended to the calculation of nuclear spin-spin coupling constants. As expected, in cases where the perturbation is small (e.g. C₂H₂) the difference between results obtained using the present scheme and results obtained using the perturbative second-order polarization propagator approach (SOPPA) is small (1-5%). However, for CH⁺ the difference is found to be more than 40%.     

A local correlation model that yields intrinsically smooth potential-energy surfaces

We demonstrate an algorithm for computing local coupled-cluster doubles (LCCD) energies that form rigorously smooth potential-energy surfaces and which should be fast enough for application to large systems in the future. Like previous LCCD algorithms, our method solves iteratively for only a limited number of correlation amplitudes, treating the remaining amplitudes with second-order perturbation theory. However, by employing bump functions, our method smoothes the transition from iteratively solved amplitude to perturbation-treated amplitude, invoking the implicit function theorem to prove that our LCCD energy is an infinitely differentiable function of nuclear coordinates. We make no explicit amplitude domains nor do we rely on the existence of atom-centered, redundant orbitals in order to get smooth potential-energy curves. In fact, our algorithm employs only localized orthonormal occupied and virtual orbitals. Our approach should be applicable to many other electron correlation methods.     

Assessing the performance of density functional theory for the electronic structure of metal-salens: the 3d(0)-metals

A series of metal-salen complexes of the 3d(0) metals Sc(III), Ti(IV), V(V), Cr(VI), and Mn(VII) have been explored using high-level electronic structure methods including coupled-cluster theory with singles, doubles, and perturbative triples as well as complete active-space third-order perturbation theory. The performance of three common density functional theory approaches has been assessed for both the geometries and the relative energies of the low-lying electronic states. The nondynamical correlation effects are demonstrated to be extremely large in all of the systems examined. Although density functional theory provides reasonable results for some of the systems, the overall agreement is quite poor. This said, the density functional theory approaches are shown to outperform the single-reference perturbation theory and coupled-cluster theory approaches for cases of strong nondynamical correlation.     

Identification of the matrix shift: a fingerprint for neutral neon complex?

Ar-NiCO and Ne-NiCO have been predicted as novel neutral noble gas charge-transfer complexes, with binding energies of 7.70 and 2.16 kcal/mol, respectively, by the highly correlated coupled-cluster singles and doubles including a perturbational estimate of triple excitations calculations. The calculated shifts in the Ni-C-O bending frequency are 48 and 36 cm(-1) for Ar-NiCO and Ne-NiCO, while the corresponding experimental matrix shifts are 46 and 36 cm(-1), respectively. The anharmonicity effects for these frequencies are verified to be very small. The interaction between a noble gas atom and NiCO is discussed through natural population analyses and the electron density difference map. We further examined the noble gas matrix effects on the geometrical structure and vibrational frequencies of NiCO by performing density functional theory calculations for the Ng31-NiCO (Ng = Ar, Ne, He) system. The present results will inspire the further experimental investigation on the complexes of noble gas and transition metal compounds generated in the matrix isolation experiments.     

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.     

Electronic structure and vertical excitation spectrum of methylene amidogen CH2N

Ab initio electronic structure calculations are reported for five electronic states of the methylene amidogen radical. Structure parameters for the ground electronic state are predicted by RHF and D -MBPT (4) calculations. Vertical excitation energies were determined using four different theoretical chemical models: complete active space (CAS ) MCSCF , CAS /MCSCF plus singles and doubles Cl, fourth-order many-body perturbation theory SDQ -MBPT (4), and coupled-cluster theory.