Hybrid coupled cluster methods: combining active space coupled cluster methods with coupled cluster singles, doubles, and perturbative triples

Based on the coupled-cluster singles, doubles, and a hybrid treatment of triples (CCSD(T)-h) method developed by us [J. Shen, E. Xu, Z. Kou, and S. Li, J. Chem. Phys. 132, 114115 (2010); and ibid. 133, 234106 (2010); and ibid. 134, 044134 (2011)], we developed and implemented a new hybrid coupled cluster (CC) method, named CCSD(T)q-h, by combining CC singles and doubles, and active triples and quadruples (CCSDtq) with CCSD(T) to deal with the electronic structures of molecules with significant multireference character. These two hybrid CC methods can be solved with non-canonical and canonical MOs. With canonical MOs, the CCSD(T)-like equations in these two methods can be solved directly without iteration so that the storage of all triple excitation amplitudes can be avoided. A practical procedure to divide canonical MOs into active and inactive subsets is proposed. Numerical calculations demonstrated that CCSD(T)-h with canonical MOs can well reproduce the corresponding results obtained with non-canonical MOs. For three atom exchange reactions, we found that CCSD(T)-h can offer a significant improvement over the popular CCSD(T) method in describing the reaction barriers. For the bond-breaking processes in F(2) and H(2)O, our calculations demonstrated that CCSD(T)q-h is a good approximation to CCSDTQ over the entire bond dissociation processes.     

Reduced multireference coupled cluster method with singles and doubles: Perturbative corrections for triples

The reduced multireference coupled-cluster method with singles and doubles (RMR CCSD) that employs multireference configuration interaction wave function as an external source for a small subset of approximate connected triples and quadruples, is perturbatively corrected for the remaining triples along the same lines as in the standard CCSD(T) method. The performance of the resulting RMR CCSD(T) method is tested on four molecular systems, namely, the HF and F(2) molecules, the NO radical, and the F(2) (+) cation, representing distinct types of molecular structure, using up to and including a cc-pVQZ basis set. The results are compared with those obtained with the standard CCSD(T), UCCSD(T), CCSD(2), and CR CCSD(T) methods, wherever applicable or available. An emphasis is made on the quality of the computed potentials in a broad range of internuclear separations and on the computed equilibrium spectroscopic properties, in particular, harmonic frequencies omega(e). It is shown that RMR CCSD(T) outperforms other triply corrected methods and is widely applicable.     

The orbital-specific-virtual local coupled cluster singles and doubles method

We extend the orbital-specific-virtual tensor factorization, introduced for local M?ller-Plesset perturbation theory in Ref. [J. Yang, Y. Kurashige, F. R. Manby and G. K. L. Chan, J. Chem. Phys. 134, 044123 (2011)], to local coupled cluster singles and doubles theory (OSV-LCCSD). The method is implemented by modifying an efficient projected-atomic-orbital local coupled cluster program (PAO-LCCSD) described recently, [H.-J. Werner and M. Schu?tz, J. Chem. Phys. 135, 144116 (2011)]. By comparison of both methods we find that the compact representation of the amplitudes in the OSV approach affords various advantages, including smaller computational time requirements (for comparable accuracy), as well as a more systematic control of the error through a single energy threshold. Overall, the OSV-LCCSD approach together with an MP2 correction yields small domain errors in practical calculations. The applicability of the OSV-LCCSD is demonstrated for molecules with up to 73 atoms and realistic basis sets (up to 2334 basis functions).     

Multireference Brillouin-Wigner coupled cluster method with singles, doubles, and triples: efficient implementation and comparison with approximate approaches

We have developed an efficient implementation of the multireference Brillouin-Wigner coupled cluster method with full iterative treatment of connected singles, doubles, and triples. Its computational costs are too high for applications to larger molecules; however, it can be used as a useful tool for benchmarking approximate methods. Performance of the method has been tested on the ground and low-lying excited states of the oxygen molecule and on the singlet-triplet gap in methylene; the results are in good agreement with experimental data.     

Gauge-origin-independent coupled cluster singles and doubles calculation of magnetic circular dichroism of azabenzenes and phosphabenzene using London orbitals

A computational study of the Faraday B term of magnetic circular dichroism at the coupled cluster singles and doubles level is presented for pyridine, pyrazine, pyrimidine, and phosphabenzene. Gauge-origin independence is obtained by expressing the B term as a total derivative of the one-photon dipole transition strength and using London orbitals. The high quality of the coupled cluster singles and doubles (CCSD) B terms makes these useful for the assignment of experimental spectra. Previous assignments of the experimental spectra based on the qualitative perimeter model are confirmed by the CCSD results for the three azines, while a reassignment is proposed for phosphabenzene. For non-overlapping bands, the B terms calculated at the equilibrium geometries are in good agreement with the experimental values. For overlapping bands, large deviations occur. Attributing a line width to the calculated equilibrium B terms leads to a large cancellation of positive and negative contributions. This cancellation may result in a large displacement of the band center maximum, leading to a large uncertainty in the assignments of "vertical experimental excitation energies" (pyridine). Bands may also completely vanish due to such cancellation (phosphabenzene). Explicit consideration of the cancellation yields simulated theoretical spectra that are in good agreement with experiment once the theoretical spectra are parallel displaced. A major contribution for this parallel displacement is the shift in the excitation energies due to correlation beyond CCSD, as seen when comparing vertical CCSD and CC3 equilibrium-geometry excitation energies.     

Explicitly time-dependent coupled cluster singles doubles calculations of laser-driven many-electron dynamics

We report explicitly time-dependent coupled cluster singles doubles (TD-CCSD) calculations, which simulate the laser-driven correlated many-electron dynamics in molecular systems. Small molecules, i.e., HF, H(2)O, NH(3), and CH(4), are treated mostly with polarized valence double zeta basis sets. We determine the coupled cluster ground states by imaginary time propagation for these molecules. Excited state energies are obtained from the Fourier transform of the time-dependent dipole moment after an ultrashort, broadband laser excitation. The time-dependent expectation values are calculated from the complex cluster amplitudes using the corresponding configuration interaction singles doubles wave functions. Also resonant laser excitations of these excited states are simulated, in order to explore the limits for the numerical stability of our current TD-CCSD implementation, which uses time-independent molecular orbitals to form excited configurations.     

Dispersion energy from density-fitted density susceptibilities of singles and doubles coupled cluster theory

A new method of calculation of the second-order dispersion energy is proposed. It is based on the Longuet-Higgins formula [Faraday Discuss. Chem. Soc. 40, 7 (1965)], which describes the dispersion interaction in terms of frequency-dependent density susceptibilities of monomers. In this study, the density susceptibilities are obtained from the coupled cluster theory at the singles and doubles level. Density fitting is applied in order to reduce the computational effort for the evaluation of density susceptibilities. It is shown that density fitting improves the scaling of the computational resources with molecular size by one order of magnitude without affecting the accuracy of the resulting dispersion energy. Numerical results are presented for several van der Waals molecules to illustrate the performance of the new approach.     

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).     

Extension of linear-scaling divide-and-conquer-based correlation method to coupled cluster theory with singles and doubles excitations

This paper describes the extension of the linear-scaling divide-and-conquer (DC)-based correlation method to the coupled cluster with singles and doubles excitations (CCSD) theory. In this DC-CCSD method, the CCSD equations are solved for all subsystems including their buffer regions with the use of the subsystem orbitals, which are obtained by the DC-Hartree-Fock method. Then, the correlation energy of the total system is evaluated by summing up the subsystem contributions other than the buffer regions by the energy density analysis technique. Numerical applications demonstrate that the present DC-CCSD gives highly accurate results with drastically less computational costs with regard to the required computer memory, scratch-disk capacity, and calculation time.     

Correction to constrained coupled cluster doubles models based on the second coupled cluster central moment

We develop a correction for the coupled cluster version of the perfect pairing (PP) model. The correction is based on finding modified values of the PP amplitudes such that the second coupled cluster central moment defined in the space of all valence single and double substitutions vanishes and, subject to this constraint, minimizing the deviation between the modified and unmodified PP amplitudes with respect to a chosen metric. We discuss how this correction can be generalized to other constrained doubles models, such as local correlation and active-space models. While the correction is not strictly size consistent and retains some of the deficiencies of the PP model, numerical results indicate that much of the missing active-space coupled cluster singles and doubles correlation energy is recovered.     

A truncation hierarchy of coupled cluster models of strongly correlated systems based on perfect-pairing references: the singles+doubles models

Paired, active-space treatments of static correlation are augmented with additional amplitudes to produce a hierarchy of parsimonious and efficient cluster truncations that approximate the total energy. The number of parameters introduced in these models grow with system size in a tractable way: two powers larger than the static correlation model it is built upon: for instance cubic for the models built on perfect pairing, fourth order for a perfect quadruples (PQ) reference, and fifth order for the models built on perfect hextuples. These methods are called singles+doubles (SD) corrections to perfect pairing, PQ, perfect hextuples, and two variants are explored. An implementation of the SD methods is compared to benchmark results for F(2) and H(2)O dissociation problems, the H(4) and H(8) model systems, and the insertion of beryllium into hydrogen. In the cases examined even the quartic number of parameters associated with PQSD is able to provide results which meaningfully improve on coupled-cluster singles doubles (CCSD) (which also has quartic amplitudes) and compete with existing multi-reference alternatives.     

Linear response coupled cluster singles and doubles approach with modified spectral resolution of the similarity transformed Hamiltonian

This paper discusses practical scheme for correcting the linear response coupled cluster with singles and doubles (CCSD) equations by shifting the poles corresponding to the equation-of-motion CCSD excitation energies by adding noniterative corrections due to triples. A simple criterion is derived for the excited states to be corrected in the spectral resolution of similarity transformed Hamiltonian on the CCSD level. Benchmark calculations were performed to compare the accuracies of static and dynamic polarizabilities obtained in this way with the CC3 and CCSDT counterparts.     

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.     

Calculation of X‐ray scattering intensities by means of the coupled cluster singles and doubles model

Total X‐ray scattering intensity σ_ee(q) is very sensitive to electron correlation effects. In this study σ_ee(q) of N₂, CO, and N₂O have been computed by the coupled cluster singles and doubles (CCSD) method and compared with configuration interaction singles and doubles (CISD) calculations as well as experimental observations. σ_ee(q) curves by CCSD calculations are rather close to those by CISD, but although small, there still exist some discrepancies between calculated and observed values. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1315–1320, 2001     

Tailored coupled cluster singles and doubles method applied to calculations on molecular structure and harmonic vibrational frequencies of ozone

To assess the separation of dynamic and nondynamic correlations and orbital choice, we calculate the molecular structure and harmonic vibrational frequencies of ozone with the recently developed tailored coupled cluster singles and doubles method (TCCSD). We employ the Hartree-Fock and complete active space (CAS) self-consistent field (SCF) orbitals to perform TCCSD calculations. When using the Hartree-Fock orbitals, it is difficult to reproduce the experimental vibrational frequency of the asymmetric stretching mode. On the other hand, the TCCSD based on the CASSCF orbitals in a correlation consistent polarized valence triple zeta basis yields excellent results with the two symmetric vibrations differing from the experimental harmonic values by 2 cm(-1) and the asymmetric vibration differing by 9 cm(-1).     

Second- and third-order triples and quadruples corrections to coupled-cluster singles and doubles in the ground and excited states

Second- and third-order perturbation corrections to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) incorporating excited configurations in the space of triples [EOM-CCSD(2)T and (3)T] or in the space of triples and quadruples [EOM-CCSD(2)TQ] have been implemented. Their ground-state counterparts--third-order corrections to coupled-cluster singles and doubles (CCSD) in the space of triples [CCSD(3)T] or in the space of triples and quadruples [CCSD(3)TQ]--have also been implemented and assessed. It has been shown that a straightforward application of the Rayleigh-Schrodinger perturbation theory leads to perturbation corrections to total energies of excited states that lack the correct size dependence. Approximations have been introduced to the perturbation corrections to arrive at EOM-CCSD(2)T, (3)T, and (2)TQ that provide size-intensive excitation energies at a noniterative O(n(7)), O(n(8)), and O(n(9)) cost (n is the number of orbitals) and CCSD(3)T and (3)TQ size-extensive total energies at a noniterative O(n(8)) and O(n(10)) cost. All the implementations are parallel executable, applicable to open and closed shells, and take into account spin and real Abelian point-group symmetries. For excited states, they form a systematically more accurate series, CCSD1 eV) and the ground-state wave function has single-determinant character. In other cases, however, the corrections tend to overestimate the triples and quadruples effects, the origin of which is discussed. For ground states, the third-order corrections lead to a rather small improvement over the highly effective second-order corrections [CCSD(2)T and (2)TQ], which is a manifestation of the staircase convergence of perturbation series.     

Calculation of two-photon absorption strengths with the approximate coupled cluster singles and doubles model CC2 using the resolution-of-identity approximation

An implementation of two-photon absorption matrix elements using the approximate second-order coupled-cluster singles and doubles model CC2 is presented. In this implementation we use the resolution-of-the-identity approximation for the two-electron repulsion integrals to reduce the computational cost. To avoid storage of large arrays we introduce in addition a numerical Laplace transformation of orbital energy denominators for the response of the doubles amplitudes. The error due to the numerical Laplace transformation is found to be negligible. Using this new implementation, we performed a series of benchmark calculations on substituted benzene and azobenzene derivatives to get reference values for TD-DFT results. We show that results obtained with the Coulomb-attenuated B3LYP functional are in reasonable agreement with the coupled-cluster results, whereas other density functionals which do not have a long-range correction give considerably less accurate results. Applications to the AF240 dye molecule and a weakly bound molecular tweezer complex demonstrate that this new RI-CC2 implementation allows for the first time to compute two-photon absorption cross sections with a correlated wave function method for molecules with more than 70 atoms and to apply this method for benchmarking TD-DFT calculations on molecules which are of particular relevance for experimental studies of two-photon absorption.     

Application of quadratic CI with singles,doubles, and triples (QCISDT): An attractive alternative to CCSDT

The size-extensive quadratic CI method with single(s), double(D ), and triple (T ) excitations, QCISDT , is compared with QCISD, QCISD (T ), CCSDT -n, and CCSDT . It is shown that QCISDT results are more accurate than are either QCISD or QCISD (T ) results. In particular, QCISDT turns out to be more stable than are QCISD and QCISD (T ) in cases with considerable multireference character. QCISDT and CCSDT results are of similar accuracy with slight advantages for the former method. Since QCISDT is much easier to implement on a computer than is CCSDT , it is an attractive alternative to CCSDT . © 1996 John Wiley & Sons, Inc.     

Fock space multireference coupled cluster method with full inclusion of connected triples for excitation energies

We report the initial Fock space multireference coupled cluster method with the full inclusion of single, double, and triple excitations (FS-CCSDT) for the (1,1) sector. We present pilot applications for calculating excitation energies for the N(2) molecule and the Ne atom. The performance of the current model, along with the FS-CCSD one, has been studied in comparison with the equation-of-motion coupled-cluster and the similarity transformed methods.