Benchmarking two-photon absorption with CC3 quadratic response theory, and comparison with density-functional response theory

We present a detailed study of the effects of electron correlation on two-photon absorption calculated by coupled cluster quadratic response theory. The hierarchy of coupled cluster models CCS, CC2, CCSD, and CC3 has been used to investigate the effects of electron correlation on the two-photon absorption cross sections of formaldehyde (CH2O), diacetylene (C4H2), and water (H2O). In particular, the effects of triple excitations on two-photon transition cross sections are determined for the first time. In addition, we present a detailed comparison of the coupled cluster results with those obtained from Hartree-Fock and density-functional response theories. We have investigated the local-density approximation, the pure Becke-Lee-Yang-Parr (BLYP) functional, the hybrid Becke-3-parameter-Lee-Yang-Parr (B3LYP), and the Coulomb-attenuated B3LYP (CAM-B3LYP) functionals. Our results show that the CAM-B3LYP functional, when used in conjuction with a one-particle basis-set containing diffuse functions, has much promise; however, care must still be exercised for diffuse Rydberg-type states.     

Coupled cluster and density functional theory studies of the vibrational contribution to the optical rotation of (S)-propylene oxide

In a previous study (Chemical Physics Letters 2005, 401 , 385) we computed the optical rotatory dispersion of (S)-propylene oxide in gas phase and solution using the hierarchy of coupled cluster models CCS, CC2, CCSD, and CC3. Even for the highly correlated CC3 model combined with a flexible basis set, the theoretical gas-phase specific rotation at 355 nm was found to be negative in contrast to the experimental result. We argued that vibrational contributions could be crucial for obtaining a complete understanding of the experimental result. Here, we show that this indeed is the case by using coupled cluster models and density functional theory methods to calculate the vibrational contributions to the gas-phase specific rotation at 355, 589.3, and 633 nm. While density functional theory (B3LYP and SAOP functionals) overestimates the specific rotation at 355 nm by approximately 1 order of magnitude and yields an incorrect sign at 589.3 and 633 nm, the coupled cluster results are in excellent agreement with the experimentally measured optical rotations. We find that all vibrational modes contribute significantly to the optical rotation and that temperature effects must be taken into account.     

Calculated electronic transitions of the water ammonia complex

We have calculated vertical excitation energies and oscillator strengths of the low lying electronic transitions in H2O, NH3, and H2ONH3 using a hierarchy of coupled cluster response functions [coupled cluster singles (CCS), second order approximate coupled cluster singles and doubles (CC2), coupled cluster singles and doubles (CCSD), and third order approximate coupled cluster singles, doubles, and triples (CC3)] and correlation consistent basis functions (n-aug-cc-pVXZ, where n=s,d,t and X=D,T,Q). Our calculations indicate that significant changes in the absorption spectra of the photodissociative states of H2O and NH3 monomers occur upon complexation. In particular, we find that the electronic transitions originating from NH3 are blueshifted, whereas the electronic transitions originating from H2O are redshifted.     

State-of-the-art computations of dipole moments using analytic gradients of high-level density-fitted coupled-cluster methods with focal-point approximations

Using the analytic derivatives approach, dipole moments of high-level density-fitted coupled-cluster (CC) methods, such as coupled-cluster singles and doubles (CCSD), and coupled-cluster singles and doubles with perturbative triples [CCSD(T)], are presented. To obtain the high accuracy results, the computed dipole moments are extrapolated to the complete basis set (CBS) limits applying focal-point approximations. Dipole moments of the CC methods considered are compared with the experimental gas-phase values, as well as with the common DFT functionals, such as B3LYP, BP86, M06-2X, and BLYP. For all test sets considered, the CCSD(T) method provides substantial improvements over Hartree–Fock (HF), by 0.076–0.213 D, and its mean absolute errors are lower than 0.06 D. Furthermore, our results indicate that even though the performances of the common DFT functionals considered are significantly better than that of HF, their results are not comparable with the CC methods. Our results demonstrate that the CCSD(T)/CBS level of theory provides highly-accurate dipole moments, and its quality approaching the experimental results. © 2019 Wiley Periodicals, Inc.     

Energy landscapes of nucleophilic substitution reactions: a comparison of density functional theory and coupled cluster methods

We have carried out a detailed evaluation of the performance of all classes of density functional theory (DFT) for describing the potential energy surface (PES) of a wide range of nucleophilic substitution (SN2) reactions involving, amongst others, nucleophilic attack at carbon, nitrogen, silicon, and sulfur. In particular, we investigate the ability of the local density approximation (LDA), generalized gradient approximation (GGA), meta-GGA as well as hybrid DFT to reproduce high-level coupled cluster (CCSD(T)) benchmarks that are close to the basis set limit. The most accurate GGA, meta-GGA, and hybrid functionals yield mean absolute deviations of about 2 kcal/mol relative to the coupled cluster data, for reactant complexation, central barriers, overall barriers as well as reaction energies. For the three nonlocal DFT classes, the best functionals are found to be OPBE (GGA), OLAP3 (meta-GGA), and mPBE0KCIS (hybrid DFT). The popular B3LYP functional is not bad but performs significantly worse than the best GGA functionals. Furthermore, we have compared the geometries from several density functionals with the reference CCSD(T) data. The same GGA functionals that perform best for the energies (OPBE, OLYP), also perform best for the geometries with average absolute deviations in bond lengths of 0.06 A and 0.6 degrees, even better than the best meta-GGA and hybrid functionals. In view of the reduced computational effort of GGAs with respect to meta-GGAs and hybrid functionals, let alone coupled cluster, we recommend the use of accurate GGAs such as OPBE or OLYP for the study of SN2 reactions.     

Comparing ab initio density-functional and wave function theories: the impact of correlation on the electronic density and the role of the correlation potential

The framework of ab initio density-functional theory (DFT) has been introduced as a way to provide a seamless connection between the Kohn-Sham (KS) formulation of DFT and wave-function based ab initio approaches [R. J. Bartlett, I. Grabowski, S. Hirata, and S. Ivanov, J. Chem. Phys. 122, 034104 (2005)]. Recently, an analysis of the impact of dynamical correlation effects on the density of the neon atom was presented [K. Jankowski, K. Nowakowski, I. Grabowski, and J. Wasilewski, J. Chem. Phys. 130, 164102 (2009)], contrasting the behaviour for a variety of standard density functionals with that of ab initio approaches based on second-order M?ller-Plesset (MP2) and coupled cluster theories at the singles-doubles (CCSD) and singles-doubles perturbative triples [CCSD(T)] levels. In the present work, we consider ab initio density functionals based on second-order many-body perturbation theory and coupled cluster perturbation theory in a similar manner, for a range of small atomic and molecular systems. For comparison, we also consider results obtained from MP2, CCSD, and CCSD(T) calculations. In addition to this density based analysis, we determine the KS correlation potentials corresponding to these densities and compare them with those obtained for a range of ab initio density functionals via the optimized effective potential method. The correlation energies, densities, and potentials calculated using ab initio DFT display a similar systematic behaviour to those derived from electronic densities calculated using ab initio wave function theories. In contrast, typical explicit density functionals for the correlation energy, such as VWN5 and LYP, do not show behaviour consistent with this picture of dynamical correlation, although they may provide some degree of correction for already erroneous explicitly density-dependent exchange-only functionals. The results presented here using orbital dependent ab initio density functionals show that they provide a treatment of exchange and correlation contributions within the KS framework that is more consistent with traditional ab initio wave function based methods.     

Benchmark calculations on the electron detachment energies of MO3* and M2O6* (M = Cr, Mo, W)

Neutral and anionic molecules of the monomers and dimers of the group VIB transition metal oxides (MO3 and M2O6) were studied with density functional theory (DFT) and coupled cluster CCSD(T) theory. Franck-Condon simulations of the photoelectron spectra were carried out for the transition from the ground state of the anion to that of the neutral molecule. Molecular structures from the DFT and CCSD(T) methods are compared. Electron detachment energies reported in the literature were evaluated. The calculated adiabatic and vertical electron detachment energies (ADEs and VDEs) were compared with the experimental results. CCSD(T) gives results within 0.12 eV for the ADEs. CCSD(T) predicts VDEs that are in error by as much as 0.3 eV for M = Cr. DFT hybrid functionals were found to give poor results for the ADEs and VDEs for M = Cr due to the substantial amount of multireference character in the wavefunction, whereas the pure DFT functionals give superior results. For M = Mo and W, excellent agreement was found for both CCSD(T) and many DFT fucntionals. The BP86 functional yields the best overall results for the VDEs of all the metal oxide clusters considered. Heats of formation calculated at the CCSD(T) level extrapolated to the complete basis set limit are also in good agreement with available experimental data.     

Theoretical study of excited states of DNA base dimers and tetramers using optimally tuned range‐separated density functional theory

Excited states of various DNA base dimers and tetramers including Watson‐Crick H‐bonding and stacking interactions have been investigated by time‐dependent density functional theory using nonempirically tuned range‐separated exchange (RSE) functionals. Significant improvements are found in the prediction of excitation energies and oscillator strengths, with results comparable to those of high‐level coupled‐cluster (CC) models (RI‐CC2 and EOM‐CCSD(T)). The optimally‐tuned RSE functional significantly outperforms its non‐tuned (default) version and widely‐used B3LYP functional. Compared to those high‐level CC benchmarks, the large mean absolute deviations of conventional functionals can be attributed to their inappropriate amount of exact exchange and large delocalization errors which can be greatly eliminated by tuning approach. Furthermore, the impacts of H‐bonding and π‐stacking interactions in various DNA dimers and tetramers are analyzed through peak shift of simulated absorption spectra as well as corresponding change of absorption intensity. The result indicates the stacking interaction in DNA tetramers mainly contributes to the hypochromicity effect. The present work provides an efficient theoretical tool for accurate prediction of optical properties and excited states of nucleobase and other biological systems. © 2015 Wiley Periodicals, Inc.     

Dispersion coefficients for first hyperpolarizabilities using coupled cluster quadratic response theory

The frequency dependence of third-order properties can in the normal dispersion region be expanded in a Taylor series in the frequency arguments. The dispersion coefficients thus obtained provide an efficient way of expressing the dispersion of frequency-dependent properties and are transferable between different optical processes. We derive analytic expressions for the dispersion coefficients of third-order properties in coupled cluster quadratic response theory and report an implementation for the three coupled cluster models CCS, CC2, and CCSD. Calculations are performed for the first hyperpolarizability of the NH₃ molecule. The convergence of the dispersion expansion with the order of the coefficients is examined and we find good convergence up to about half the frequency at which the first pole in the hyperpolarizability occurs. Padé approximants improve the convergence dramatically and extend the application range of the dispersion expansion to frequencies close to the first pole. The sensitivity of the dispersion coefficients on the dynamic correlation treatment and on the choice of the one-electron basis set is investigated. The results demonstrate that, contrary to presumptions in the literature, the dispersion coefficients are sensitive to basis set effects and correlation treatment similar to the static hyperpolarizabilities. Received: 26 March 1998 / Accepted: 21 July 1998 / Published online: 19 October 1998     

A comparative assessment of the perturbative and renormalized coupled cluster theories with a noniterative treatment of triple excitations for thermochemical kinetics, including a study of basis set and core correlation effects

The CCSD, CCSD(T), and CR-CC(2,3) coupled cluster methods, combined with five triple-zeta basis sets, namely, MG3S, aug-cc-pVTZ, aug-cc-pV(T+d)Z, aug-cc-pCVTZ, and aug-cc-pCV(T+d)Z, are tested against the DBH24 database of diverse reaction barrier heights. The calculations confirm that the inclusion of connected triple excitations is essential to achieving high accuracy for thermochemical kinetics. They show that various noniterative ways of incorporating connected triple excitations in coupled cluster theory, including the CCSD(T) approach, the full CR-CC(2,3) method, and approximate variants of CR-CC(2,3) similar to the triples corrections of the CCSD(2) approaches, are all about equally accurate for describing the effects of connected triply excited clusters in studies of activation barriers. The effect of freezing core electrons on the results of the CCSD, CCSD(T), and CR-CC(2,3) calculations for barrier heights is also examined. It is demonstrated that to include core correlation most reliably, a basis set including functions that correlate the core and that can treat core-valence correlation is required. On the other hand, the frozen-core approximation using valence-optimized basis sets that lead to relatively small computational costs of CCSD(T) and CR-CC(2,3) calculations can achieve almost as high accuracy as the analogous fully correlated calculations.     

Static dipole polarizability and binding energy of sodium clusters Nan (n=1-10): A critical assessment of all-electron based post Hartree-Fock and density functional methods

A systematic all electron post Hartree-Fock as well as density functional theory (DFT) based calculations for the polarizability and binding energy of sodium metal clusters have been performed and an in-depth analysis of the discrepancy between the experimental and theoretical results is presented. A systematic investigation for the assessment of different DFT exchange-correlation functionals in predicting the polarizability values has also been reported. All the pure DFT functionals have been found to considerably underestimate the calculated polarizability values as compared to the MP2 results. DFT calculations using the full Hartree-Fock exchange along with one-parameter progressive correlation functional have, however, been shown to yield results in good agreement with the MP2 and experimental results. The possible sources of error present in the experimental measurements as well as in the different theoretical methods have also been analyzed. One of the most important conclusions of the present study is that the effect of electron correlation plays a significant role in determining the polarizability of the clusters and the MP2 method can be considered to be one of the most reliable methods for their prediction. It has also been noted that the polarizability value of the lower member clusters (Na2 and Na4) calculated by highly sophisticated methods such as, CCSD and CCSD(T) are found to be very close to the corresponding MP2 values. The polarizability and the binding energy of the clusters are found to be inversely related to each other and their correlation is rationalized by invoking the minimum polarizability principle. A good linear correlation between the polarizability and volume of the cluster has also been found to exist.     

Coupled cluster and DFT calculations of 14N nuclear quadrupole coupling constants

The dependence of ¹⁴N quadrupole coupling constants calculated using coupled cluster theory on the level of approximation is examined for a series of small molecules. For HCN, HNC, CH₃CN, and CH₃NC, we use the coupled cluster singles‐and‐doubles with a noniterative perturbative triples correction—CCSD(T)—approach, and we analyze the basis set dependence of the results. For aziridine, diazirine, and cyclopropyl cyanide, we use the CCSD(T) approach, but smaller basis sets, and for the largest studied molecules—quinuclidine and hexamine—we present CCSD results. The differences between computed and experimental values for the best basis sets used are ≈ 5% at the CCSD level and decrease noticeably at the CCSD(T) level. The ‐ N≡C bonds are an exception—in this case the quadrupole coupling constants are very small, hence the differences between theory and experiment become larger (up to 9%). We also consider the performance of density functional theory, comparing the results for different density functionals with the coupled cluster values of the same constants. Most of the functionals provide results systematically improved with respect to the Hartree–Fock values, with ¹⁴N coupling constants in ‐ N≡C bonds being again an exception. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Communication: convergence of anharmonic infrared intensities of hydrogen fluoride in traditional and explicitly correlated coupled cluster calculations

It is shown that the convergence of anharmonic infrared spectral intensities with respect to the basis set size is much enhanced in explicitly correlated calculations as compared to traditional configuration interaction type wave function expansion. Explicitly correlated coupled cluster (CC) calculations using Slater-type geminal correlation factor (CC-F12) yield well-converged dipole derivatives and vibrational intensities for hydrogen fluoride with basis set involving f functions on the heavy atom. Combination of CC-F12 with singles, doubles, and non-iterative triples (CCSD(T)-F12) with small corrections due to quadruple excitations, core-electron correlation, and relativistic effects yields vibrational line positions, dipole moments, and transition dipole matrix elements in good agreement with the best experimental values.     

Mono- and dinitrosyls on copper(I) site in a zeolite model: effects of static correlation

Multiconfigurational RASSCF/RASPT2 approach has been applied to investigate bonding of one and two nitric oxide (NO) molecules to a simple model of Cu(I) site in zeolite environment, Cu(I)[Al(OH)(4)]. Two binding modes were considered for the mononitrosyls and four alternative structures for the dinitrosyls (each one in either singlet or triplet state). Stabilities of the mono- and dinitrosyl complexes obtained from the multireference calculations were compared to the previously reported coupled cluster CCSD(T) results, as well as to DFT calculations performed here with various functionals, either hybrid or nonhybrid ones. RASSCF calculations provided also a qualitative insight into the electronic structure of the studied complexes, concerning mainly the interaction between the Cu and the NO ligand, and between the two NO fragments. Whereas the electronic structure of the mononitrosyls is dominated by a single configuration, the dinitrosyls have a considerably multireference character. Various effects of nondynamical correlation have been pointed out for these interesting species, trying to assess their impact on performance of the tested DFT methods.     

Ab initio calculations on halogen-bonded complexes and comparison with density functional methods

A systematic theoretical investigation on a series of dimeric complexes formed between some halocarbon molecules and electron donors has been carried out by employing both ab initio and density functional methods. Full geometry optimizations are performed at the Moller-Plesset second-order perturbation (MP2) level of theory with the Dunning's correlation-consistent basis set, aug-cc-pVDZ. Binding energies are extrapolated to the complete basis set (CBS) limit by means of two most commonly used extrapolation methods and the aug-cc-pVXZ (X = D, T, Q) basis sets series. The coupled cluster with single, double, and noniterative triple excitations [CCSD(T)] correction term, determined as a difference between CCSD(T) and MP2 binding energies, is estimated with the aug-cc-pVDZ basis set. In general, the inclusion of higher-order electron correlation effects leads to a repulsive correction with respect to those predicted at the MP2 level. The calculations described herein have shown that the CCSD(T) CBS limits yield binding energies with a range of -0.89 to -4.38 kcal/mol for the halogen-bonded complexes under study. The performance of several density functional theory (DFT) methods has been evaluated comparing the results with those obtained from MP2 and CCSD(T). It is shown that PBEKCIS, B97-1, and MPWLYP functionals provide accuracies close to the computationally very expensive ab initio methods.     

A coupled cluster approach with a hybrid treatment of connected triple excitations: implementation and applications for open-shell systems

An implementation of the coupled cluster (CC) singles, doubles, and a hybrid treatment of connected triples [denoted as CCSD(T)-h], based on the unrestricted Hartree-Fock (UHF) reference, is presented. Based on the spin-integrated formulation, we have developed a computer program to achieve the automatic derivation and implementation of the CCSD(T)-h approach. The CCSD(T)-h approach computationally scales as the seventh power of the system size, and is affordable for many medium-sized systems. The present approach has been applied to study the equilibrium geometries and harmonic vibrational frequencies in a number of open-shell diatomic molecules and bond breaking potential energy profiles in several open-shell molecules, including CH(3), NH(2), and SiH(2). For all systems under study, the overall performance of the UHF-based CCSD(T)-h approach is very close to that of the corresponding CCSDT (CC singles, doubles, and triples), and much better than that of the UHF-based CCSD(T) (CC singles, doubles, and perturbative triples).