Locally correlated equation-of-motion coupled cluster theory for the excited states of large molecules

We report an extension of the local correlation concept to electronically excited states via the equation-of-motion coupled cluster singles and doubles (EOM-CCSD) method. We apply the same orbital domain structure used successfully for ground-state CCSD by Werner and co-workers and find that the resulting localized excitation energies are in error generally by less than 0.2 eV relative to their canonical EOM-CCSD counterparts, provided the basis set is flexible and includes Rydberg-like functions. In addition, we account for weak-pair contributions efficiently using a correction to local-EOM-CCSD transition energies based on the perturbative (D) correction used with configuration interaction singles (CIS).     

Local CC2 response method for triplet states based on Laplace transform: excitation energies and first-order properties

A new multistate local CC2 response method for calculating excitation energies and first-order properties of excited triplet states in extended molecular systems is presented. The Laplace transform technique is employed to partition the left/right local CC2 eigenvalue problems as well as the linear equations determining the Lagrange multipliers needed for the properties. The doubles part in the equations can then be inverted on-the-fly and only effective equations for the singles part must be solved iteratively. The local approximation presented here is adaptive and state-specific. The density-fitting method is utilized to approximate the electron-repulsion integrals. The accuracy of the new method is tested by comparison to canonical reference values for a set of 12 test molecules and 62 excited triplet states. As an illustrative application example, the lowest four triplet states of 3-(5-(5-(4-(bis(4-(hexyloxy)phenyl)amino)phenyl)thiophene-2-yl)thiophene-2-yl)-2-cyanoacrylic acid, an organic sensitizer for solar-cell applications, are computed in the present work. No triplet charge-transfer states are detected among these states. This situation contrasts with the singlet states of this molecule, where the lowest singlet state has been recently found to correspond to an excited state with a pronounced charge-transfer character having a large transition strength.     

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

Vibrational excitation energies from vibrational coupled cluster response theory

Response theory in the context of vibrational coupled cluster (VCC) theory is introduced and used to obtain vibrational excitation energies. The relation to the vibrational configuration interaction (VCI) approach is described, and the increase in accuracy of VCC response energies relative to VCI energies is discussed theoretically in terms of a perturbational order expansion and demonstrated numerically. To illustrate the theory, a pilot implementation is used to obtain anharmonic vibrational frequencies for fundamental, first overtone and combination excitations of formaldehyde as well as for the fundamental transitions of ethylene.     

Study of approximate coupled cluster methods for first-order static properties

In this paper, we analyse the algebraic structure of the equations for calculating the first order static properties using several approximate versions of Coupled Cluster (CC) methods. In particular, the non-variational and the variational method using a CC wavefunction corresponding to an appropriately defined perturbed Hamiltonian as well as the simple expectation value expression using a CC stationary state are studied under different approximations. Two different models are proposed: (a) use of maximum overlap orbitals where the pertinent approximations are TT ₂, T ⁽¹⁾ T ₂ ⁽¹⁾, (b) use of Hartree-Fock orbitals and T(T ₁+T ₂), T ⁽¹⁾(T ₁ ⁽¹⁾ +T ₂ ⁽¹⁾ ) approximations. It is analytically shown that in both these models certain approximate versions of the methods under purview yield identical results for first order static properties.NCL Communication No. 3725     

Two-particle density matrix cumulant of coupled cluster theory

A two-particle density matrix obtained from the expectation value expression for the coupled cluster wave function is separated into the antisymmetrized product of the one-particle density matrices and the remaining cumulant part. It is demonstrated that the proposed formula for the coupled cluster cumulant is a valid cumulant expression, since it is a connected, and therefore size-extensive quantity. It is also shown that the density matrices from coupled cluster gradient theory, when used to define a cumulant, result in the expression containing disconnected terms. The proposed formulation of the coupled cluster cumulant makes it easy to develop size-extensive truncation schemes. As an example, explicit equations for the cumulant at the coupled cluster single, double and triple excitation level are presented.     

Localized optimized orbitals, coupled cluster theory, and chiroptical response properties

The impact of orbital localization on the efficiency and accuracy of the optimized-orbital coupled cluster model is examined for the prediction of chiroptical properties, in particular optical rotation. The specific rotations of several test cases-(P)-[4]triangulane, (S)-1-phenylethanol, and chiral conformers of 1-fluoropentane, heptane, and nonane-were computed using an approach in which localization is enforced throughout the orbital optimization and subsequent linear response computation. This method provides a robust local-correlation scheme for future production-level implementation. Although the cross-over point between the canonical and localized coupled cluster approach lies at larger molecules than for ground-state energies, the scheme presented should still provide reduced scaling sufficient to investigate much larger molecules than are presently accessible.     

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.     

Local CC2 electronic excitation energies for large molecules with density fitting

A new local method for the computation of electronic excitation energies of singlet states in extended molecular systems is presented. It is based on the CC2 model and local approximations to the wave functions. In the proposed method the singles excitations are treated nonlocally and local restrictions are imposed on doubles amplitudes only. The accuracy of the new method was tested by calculating several lowest excited states for 14 molecules and comparing them with canonical CC2 values. Deviations of the local excitation energies from the canonical reference values do not exceed 0.05 eV for all test molecules and all states in the lower energy range investigated in this work. The method uses the density-fitting approximation for all two-electron integrals, which considerably simplifies the computational complexity of the individual diagrams. A combination of the local approximations and the powerful density-fitting technique leads to a low-scaling method, capable to treat molecular systems comprised of 100 atoms and more in a basis of a polarized double zeta quality. A test calculation for a system consisting of 127 atoms and 370 active electrons without symmetry is presented to show the efficiency of the new method.     

Double-hybrid density functional theory for excited electronic states of molecules

Double-hybrid density functionals are based on a mixing of standard generalized gradient approximations (GGAs) for exchange and correlation with Hartree-Fock (HF) exchange and a perturbative second-order correlation part (PT2) that is obtained from the Kohn-Sham (GGA) orbitals and eigenvalues. This virtual orbital-dependent functional (dubbed B2PLYP) contains only two empirical parameters that describe the mixture of HF and GGA exchange (ax) and of the PT2 and GGA correlation (ac), respectively. Extensive testing has recently demonstrated the outstanding accuracy of this approach for various ground state problems in general chemistry applications. The method is extended here without any further empirical adjustments to electronically excited states in the framework of time-dependent density functional theory (TD-DFT) or the closely related Tamm-Dancoff approximation (TDA-DFT). In complete analogy to the ground state treatment, a scaled second-order perturbation correction to configuration interaction with singles (CIS(D)) wave functions developed some years ago by Head-Gordon et al. [Chem. Phys. Lett. 219, 21 (1994)] is computed on the basis of density functional data and added to the TD(A)-DFTGGA excitation energy. The method is implemented by applying the resolution of the identity approximation and the efficiency of the code is discussed. Extensive tests for a wide variety of molecules and excited states (of singlet, triplet, and doublet multiplicities) including electronic spectra are presented. In general, rather accurate excitation energies (deviations from reference data typically <0.2 eV) are obtained that are mostly better than those from standard functionals. Still, systematic errors are obtained for Rydberg (too low on average by about 0.3 eV) and charge-transfer transitions but due to the relatively large ax parameter (0.53), B2PLYP outperforms most other functionals in this respect. Compared to conventional HF-based CIS(D), the method is more robust in electronically complex situations due to the implicit account of static correlation effects by the GGA parts. The (D) correction often works in the right direction and compensates for the overestimation of the transition energy at the TD level due to the elevated fraction of HF exchange in the hybrid GGA part. Finally, the limitations of the method are discussed for challenging systems such as transition metal complexes, cyanine dyes, and multireference cases.     

Benchmarks for electronically excited states: CASPT2, CC2, CCSD, and CC3

A benchmark set of 28 medium-sized organic molecules is assembled that covers the most important classes of chromophores including polyenes and other unsaturated aliphatic compounds, aromatic hydrocarbons, heterocycles, carbonyl compounds, and nucleobases. Vertical excitation energies and one-electron properties are computed for the valence excited states of these molecules using both multiconfigurational second-order perturbation theory, CASPT2, and a hierarchy of coupled cluster methods, CC2, CCSD, and CC3. The calculations are done at identical geometries (MP26-31G*) and with the same basis set (TZVP). In most cases, the CC3 results are very close to the CASPT2 results, whereas there are larger deviations with CC2 and CCSD, especially in singlet excited states that are not dominated by single excitations. Statistical evaluations of the calculated vertical excitation energies for 223 states are presented and discussed in order to assess the relative merits of the applied methods. CC2 reproduces the CC3 reference data for the singlets better than CCSD. On the basis of the current computational results and an extensive survey of the literature, we propose best estimates for the energies of 104 singlet and 63 triplet excited states.     

Explicitly correlated second-order perturbation theory using density fitting and local approximations

Three major obstacles in electronic structure theory are the steep scalings of computer time with respect to system size and basis size and the slow convergence of correlation energies in orbital basis sets. Three solutions to these are, respectively, local methods, density fitting, and explicit correlation; in this work, we combine all three to produce a low-order scaling method that can achieve accurate MP2 energies for large systems. The errors introduced by the local approximations into the R12 treatment are analyzed for 16 chemical reactions involving 21 molecules. Weak pair approximations, as well as local resolution of the identity approximations, are tested for molecules with up to 49 atoms, over 100 correlated electrons, and over 1000 basis functions.     

Theoretical prediction of the electronic excited states and resonance Raman intensities in formamide from coupled cluster calculations

Electronic excitations and the resonance Raman spectrum of formamide were obtained from ab initio electron correlation calculations using the equation of motion coupled cluster (EOM-CCSD) method. Interpretation of the UV spectrum on the basis of calculated vertical excitation energies and oscillator strengths accounts for all experimental bands previously assigned. Our assignment, however, suggests an additional Rydberg band at about 7.4 eV which may be hidden under the main absorption. We also show that the Rydberg states appear pairwise, corresponding to n and π hole states, respectively. Using analytic derivative techniques, derivatives of the excited state energies with respect to normal coordinates of the ground state were calculated. Approximate resonance Raman intensities have been determined.     

Vibrationally excited states of DC5N: Millimeter-wave spectroscopy and coupled cluster calculations

The rotational spectrum of DC₅N has been investigated in the millimeter-wave region for 16 vibrationally excited states which approximately lie below 760 cm⁻¹, namely (v₆ v₇ v₈ v₉ v₁₀ v₁₁)=(000001), (000002), (000003), (000005), (000010), (000020), (000100) (001000), (010000), (100000), (000011), (000101), (001001), (010001), (001010), and (010010). Gas-phase copyrolysis of fully deuterated pyridine and phosphorus trichloride was used to produce the semi-stable DC₅N molecule. In addition to the usual l-type resonances, several vibrational interactions have been taken into account to fit properly the measured transition frequencies. The most important perturbations are caused by the cubic anharmonic interactions which mix the v₆ stretching state with the 2v₁₀ overtone and the v₈+v₁₁ and v₇+v₁₁ bending combination states. The analysis of the spectra was facilitated by theoretical predictions from CCSD(T) calculations with the cc-pVQZ basis.     

The accuracy of dipole moments from spin-component scaled CC2 in ground and electronically excited states

The accuracy of dipole moments calculated from wave function methods based on second-order perturbation theory is investigated in the ground and electronically excited states. Results from the approximate coupled-cluster singles-and-doubles model, CC2, M?ller-Plesset perturbation theory, MP2, and the algebraic diagrammatic construction through second-order, ADC(2), are discussed together with the spin-component scaled and the scaled opposite-spin variants of these methods. The computed dipole moments show a very good correlation with data from high-resolution spectroscopy. Compared to the unscaled methods, the spin-component scaling increases the accuracy of the results and improves the robustness of the calculations. An accuracy about 0.2 to 0.1 D in the ground state and about 0.3 to 0.2 D in the electronically excited states can be achieved with these approaches.     

Scaled opposite-spin CC2 for ground and excited states with fourth order scaling computational costs

An implementation of scaled opposite-spin CC2 (SOS-CC2) for ground and excited state energies is presented that requires only fourth order scaling computational costs. The SOS-CC2 method yields results with an accuracy comparable to the unscaled method. Furthermore the time-determining fifth order scaling steps in the algorithm can be replaced by only fourth order scaling computational costs using a "resolution of the identity" approximation for the electron repulsion integrals and a Laplace transformation of the orbital energy denominators. This leads to a significant reduction of computational costs especially for large systems. Timings for ground and excited state calculations are shown and the error of the Laplace transformation is investigated. An application to a chlorophyll molecule with 134 atoms results in a speed-up by a factor of five and demonstrates how the new implementation extends the applicability of the method. A SOS variant of the algebraic diagrammatic construction through second order ADC(2), which arises from a simplification of the SOS-CC2 model, is also presented. The SOS-ADC(2) model is a cost-efficient alternative in particular for future extensions to spectral intensities and excited state structure optimizations.     

Analytical Hessian of electronic excited states in time-dependent density functional theory with Tamm-Dancoff approximation

We present the analytical expression and computer implementation for the second-order energy derivatives of the electronic excited state with respect to the nuclear coordinates in the time-dependent density functional theory (TDDFT) with Gaussian atomic orbital basis sets. Here, the Tamm-Dancoff approximation to the full TDDFT is adopted, and therefore the formulation process of TDDFT excited-state Hessian is similar to that of configuration interaction singles (CIS) Hessian. However, due to the replacement of the Hartree-Fock exchange integrals in CIS with the exchange-correlation kernels in TDDFT, many quantitative changes in the derived equations are arisen. The replacement also causes additional technical difficulties associated with the calculation of a large number of multiple-order functional derivatives with respect to the density variables and the nuclear coordinates. Numerical tests on a set of test molecules are performed. The simulated excited-state vibrational frequencies by the analytical Hessian approach are compared with those computed by CIS and the finite-difference method. It is found that the analytical Hessian method is superior to the finite-difference method in terms of the computational accuracy and efficiency. The numerical differentiation can be difficult due to root flipping for excited states that are close in energy. TDDFT yields more exact excited-state vibrational frequencies than CIS, which usually overestimates the values.     

Addition by subtraction in coupled cluster theory. II. Equation-of-motion coupled cluster method for excited, ionized, and electron-attached states based on the nCC ground state wave function

New iterative double and triple excitation corrections to the equation-of-motion coupled cluster (EOM-CC) based upon the recently developed nCC methods [Bartlett and Musia?, J. Chem. Phys. 125, 204105-1 (2006)] are applied to excitation energies (EEs), ionization potentials (IPs), and electron affinities (EAs). The methods have been tested by the evaluation of the vertical EEs, IPs, and EAs for Ne, BH, CH(2), H(2)O, N(2), C(2), CH(+), CO, and C(2)H(4) compared to full configuration interaction, EOM-CCSD, EOM-CCSDT, and experimental data.