Linear scaling multireference singles and doubles configuration interaction

A linear scaling multireference singles and doubles configuration interaction (MRSDCI) method has been developed. By using localized bases to span the occupied and virtual subspace, local truncation schemes can be applied in tandem with integral screening to reduce the various bottlenecks in a MRSDCI calculation. Among these, the evaluation of electron repulsion integrals and their subsequent transformation, together with the diagonalization of the large CI Hamiltonian matrix, correspond to the most computationally intensive steps in a MRSDCI calculation. We show that linear scaling is possible within each step. The scaling of the method with system size is explored with a system of linear alkane chains and we proceed to demonstrate this method can produce smooth potential energy surfaces via calculating the dissociation of trans-6-dodecene (C(12)H(24)) along the central C[Double Bond]C bond.     

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.     

Solvent effects in electronically excited states using the continuum solvation model COSMO in combination with multireference configuration interaction with singles and doubles (MR-CISD)

An implementation of the COSMO continuum solvation model into the MCSCF and MR-CISD programs of the COLUMBUS program system is reported. Equilibrium solvation and non-equilibrium solvation models for the treatment of electronic excitations have been used. Solvatochromic effects have been computed for a representative set of n-* and -* states of formaldehyde, acrolein and pyrazine using several solvents ranging from some with apolar character to water. Agreement with experimental shifts is good within the limits of a continuum model.Contribution to the Jacopo Tomasi Honorary Issue     

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.     

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.     

An ab initio multireference doubles excitation configuration interaction study of low-lying electronic states of Cd2 using Slater-type orbitals

Potential-energy curves for the ground state and lower excited states of the Cd₂ dimer have been calculated. They are obtained using a multireference doubles excitation configuration interaction procedure and employing Slater basis sets, previously optimized at the self-consistent-field level for excited states of the Cd atom. The spectroscopic constants and excitation energies for the bound states of Cd₂ have been compared with experimental as well as other theoretical results. The ground state of Cd₂ is essentially repulsive and presents a shallow van der Waals minimum. The computed adiabatic electronic transitions are in good agreement with the experimental ones. Received: 16 September 1999 / Accepted: 3 February 2000 / Published online: 2 May 2000     

A multireference configuration interaction study with singles and doubles of some mesoionic rings: reaction and activation free energies for the ring-opening reaction

For the first time, multireference configuration interaction with singles and doubles (MR-CISD) calculations (including extensivity corrections, MR-CISD+Q) have been performed to study the ring-opening reactions of the following mesoionic rings and their 2, 3, and 4 methyl-substituted derivatives: ( a ) 1,3-oxazol-5-one, ( c ) 1,3-oxazol-5-thiolate, and ( d ) 1,3-thiazol-5-thiolate. The ring-opening reaction of the parent 1,3-thiazol-5-one mesoionic ring ( b ) has also been studied. The effect of methyl and S replacement on the reaction and activation free energies (ΔG and ΔG^†, respectively) is studied. For a , the effect of methyl replacement on C2 is almost negligible, while on N3 and C4, it is significantly larger, especially in the latter position. The open structure (ketene form) of a is considerably more stable, and the replacement of the exocyclic O by S stabilizes the ketene tautomer even more and increases ΔG^† considerably. On the other hand, replacement of the endocyclic O by S (yielding b ) completely prevents the formation of the open structure, while replacement of both endo and exo atoms strongly stabilizes the cyclic structure and leads to a remarkably high value of ΔG^†. In some cases, the relative polarity of the stationary points is used to estimate how ΔG and ΔG^† are expected to change as one goes from the gas phase to aprotic polar solvents. A linear correlation between the multiconfigurational character of the mesoionic rings and ΔG^† has been observed, and such correlation is used to estimate the ΔG^† values for the methyl-substituted c and d . The largest barrier (13.78 kcal/mol) has been obtained for the 4-methyl-1,3-thiazol-5-thiolate (4-m- d ).     

CASSCF and multireference CI with singles and doubles study of low-lying valence and Rydberg states of 2H-tetrazole

Complete active space self-consistent field (CASSCF) and multireference CI with singles and doubles (MR-CISD) calculations [including extensivity corrections, at MR-CISD+Q and multireference averaged quadratic coupled cluster (MR-AQCC) levels] have been performed to characterize the low-lying valence and the Rydberg states of 2H-tetrazole. The highest level results (MR-AQCC/d'-aug'-cc-pVDZ) indicate the following ordering of the valence singlet excited states: S(1) (n-pi*), 6.06 eV; S(2) (n-pi*), 6.55 eV; S(3) (pi-pi*), 6.55 eV. The MR-CISD+Q/d'-aug'-cc-pVDZ results indicate the same ordering, but at slight higher energies: 6.16, 6.68, and 6.69 eV, respectively. According to our MR-CISD+Q/d'-aug'-cc-pVDZ results, the next two states are Rydberg states, at 7.69 eV (pi-3s) and 7.89 eV (n-3s). The calculated energies of these two states, as well as their proximity, are consistent with the conclusion reached by Palmer and Beveridge (Chem Phys 1987, 111, 249) that the first band of the photoelectron spectrum of 2H-tetrazole is likely to be associated to the first two ionizations processes (of pi and N lone pair electrons), at energies close to 11.3 eV.     

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.     

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.     

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.     

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.     

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.     

One-electron properties and electrostatic interaction energies from the expectation value expression and wave function of singles and doubles coupled cluster theory

One-electron density matrices resulting from the explicitly connected commutator expansion of the expectation value were implemented at the singles and doubles coupled cluster (CCSD) level. In the proposed approach the one-electron density matrix is obtained at a little extra cost in comparison to the calculation of the CCSD correlation energy. Therefore, in terms of the computational time the new method is significantly less demanding than the conventional linear-response CCSD theory which requires additionally an expensive calculation of the left-hand solution of the CCSD equations. The quality of the new density matrices was investigated by computing a set of one-electron properties for a series of molecules of varying sizes and comparing the results with data obtained using the full configuration interaction method or higher level coupled cluster theory. It has been found that the results obtained using the new approach are of the same quality as those predicted by the linear-response CCSD method. The novel one-electron density matrices have also been applied to study the energy of the electrostatic interaction for a number of van der Waals complexes, including the benzene and azulene dimers.     

Iterative Cl general singles and doubles (ICIGSD) method for calculating the exact wave functions of the ground and excited states of molecules

The iterative configuration-interaction general singles and doubles (ICIGSD) method was applied to various closed- and open-shell electronic states of molecules within finite basis sets and was shown to give the exact results that are identical to the full CI ones. The structure of the ICIGSD is unique among the ICI formalisms, that is, the singularity problem intrinsic to atomic and molecular Hamiltonians can be avoided. The convergence of the ICIGSD method was fairly good regardless of the characters of the electronic states and the qualities of the basis sets; only several iterations were enough for obtaining microhartree accuracy. These favorable properties are attributed to the unique GSD structure. The present method was shown to be applicable to various spin states and to quasidegenerate states appearing in bond dissociation process. We have also applied the ICIGSD-CI method to calculate the excited states simultaneously. We have confirmed that the ICIGSD-CI method is accurate for calculating the excited states the symmetries of which are not only similar to but also different from that of the ground state.     

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.     

Polarizability and optical rotation calculated from the approximate coupled cluster singles and doubles CC2 linear response theory using Cholesky decompositions

A new implementation of the approximate coupled cluster singles and doubles CC2 linear response model using Cholesky decomposition of the two-electron integrals is presented. Significantly reducing storage demands and computational effort without sacrificing accuracy compared to the conventional model, the algorithm is well suited for large-scale applications. Extensive basis set convergence studies are presented for the static and frequency-dependent electric dipole polarizability of benzene and C60, and for the optical rotation of CNOFH2 and (-)-trans-cyclooctene (TCO). The origin-dependence of the optical rotation is calculated and shown to persist for CC2 even at basis set convergence.     

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