Showing papers on "Coupled cluster published in 2005"

Systematically convergent basis sets for transition metals. I. All-electron correlation consistent basis sets for the 3d elements Sc–Zn

Abstract: Sequences of basis sets that systematically converge towards the complete basis set (CBS) limit have been developed for the first-row transition metal elements Sc–Zn. Two families of basis sets, nonrelativistic and Douglas-Kroll-Hess (-DK) relativistic, are presented that range in quality from triple-ζ to quintuple-ζ. Separate sets are developed for the description of valence (3d4s) electron correlation (cc-pVnZ and cc-pVnZ-DK; n=T,Q, 5) and valence plus outer-core (3s3p3d4s) correlation (cc-pwCVnZ and cc-pwCVnZ-DK; n=T,Q, 5), as well as these sets augmented by additional diffuse functions for the description of negative ions and weak interactions (aug-cc-pVnZ and aug-cc-pVnZ-DK). Extensive benchmark calculations at the coupled cluster level of theory are presented for atomic excitation energies, ionization potentials, and electron affinities, as well as molecular calculations on selected hydrides (TiH, MnH, CuH) and other diatomics (TiF, Cu2). In addition to observing systematic convergence towards the C...

Theory and applications of computational chemistry : the first forty years

Abstract: Computing Technologies, Theories, and Algorithms. The Making of 40 Years and More of Theoretical and Computational. A Dynamical, Time-Dependent View of Molecular Theory. Computation of Non-covalent Binding. Electrodynamics in Computational Chemistry. Variational Transition State. Attempting to Simulate Large Molecular Systems. The Beginnings of Coupled Cluster Theory: An Eyewitness Account. Controlling Quantum Phenomena with Photonic Reagents. First-Principles Calculations of Anharmonic Vibrational Spectroscopy of Large Molecules. Finding Minima, Transition States, and Following Reaction Pathways on Ab Initio Potential Energy Surfaces. Progress in the Quantum Description of Vibrational Motion of Polyatomic Molecules. Toward Accurate Computations in Photobiology. The Nature of the Chemical Bond in the Light of an Energy Decomposition Analysis. Superoperator Many-body Theory of Molecular Currents: Non-equilibrium Green Functions in Real Time. Role of Computational Chemistry in the Theory of Unimolecular Reaction Rates. Molecular Dynamics: An Account of its Evolution. Equations of Motion (EOM) Methods for Computing Electron Affinities and Ionization Potentials. Multireference Coupled Cluster Method Based on the Brillouin-Wigner Perturbation Theory. Electronic Structure: The Momentum Perspective. Recent Advances in ab initio, DFT, and Relativistic Electronic Structure Theory . Semiempirical Quantum-Chemical Methods in Computational Chemistry. Size-consistent State-specific Multi-reference Methods: A Survey of Some Recent Developments. The Valence Bond Diagram Approach - A Paradigm for Chemical Reactivity. Development of Approximate Exchange-Correlation Functionals. Multiconfigurational Quantum Chemistry. Concepts of Perturbation, Orbital interaction, Orbital Mixing and Orbital Occupation. G2, G3 and Associated Quantum Chemical Models for Accurate Theoretical Thermochemistry. Factors that Affect Conductance at the Molecular Level. The CH O Hydrogen Bond. A Historical Account. Ab Initio and DFT Calculations on the Cope Rearrangement, a Reaction with a Chameleonic Transition State. High-Temperature Quantum Chemical Molecular Dynamics Simulations of Carbon Nanostructure Self-Assembly Processes. Computational Chemistry of Isomeric Fullerenes and Endofullerenes. On the importance of Many-Body Forces in Clusters and Condensed Phase. Clusters to Functional Molecules, Nanomaterials, and Molecular Devices: Theoretical Exploration. Monte Carlo Simulations of the Finite Temperature Properties of (H2O)6. Computational Quantum Chemistry on Polymer Chains: Aspects of the Last Half Century. Forty Years of Ab Initio Calculations on Intermolecular. Applied Density Functional Theory and the deMon Codes 1964 to 2004. SAC-CI Method Applied to Molecular Spectroscopy. Forty Years of Fenske-Hall Molecular Orbital Theory. Advances in Electronic Structure Theory: GAMESS a Decade Later. How and Why Coupled-Cluster Theory Became the Preeminent Method in Ab initio Quantum Chemistry.

Systematically convergent basis sets for transition metals. II. Pseudopotential-based correlation consistent basis sets for the group 11 (Cu, Ag, Au) and 12 (Zn, Cd, Hg) elements

Abstract: Sequences of basis sets that systematically converge towards the complete basis set (CBS) limit have been developed for the coinage metals (Cu, Ag, Au) and group 12 elements (Zn, Cd, Hg). These basis sets are based on recently published small-core relativistic pseudopotentials [Figgen D, Rauhut G, Dolg M, Stoll H (2005) Chem Phys 311:227] and range in size from double- through quintuple-ζ. Series of basis sets designed for valence-only and outer-core electron correlation are presented, as well as these sets augmented by additional diffuse functions for the accurate description of negative ions and weak interactions. Selected benchmark calculations at the coupled cluster level of theory are presented for both atomic and molecular properties. The latter include the calculation of both spectroscopic and thermochemical properties of the homonuclear dimers Cu2, Ag2, and Au2, as well as the van der Waals species Zn2, Cd2, and Hg2. The CBS limit results, including the effects of core-valence correlation and spin-orbit coupling, represent some of the most accurate carried out to date and result in new recommendations for the equilibrium bond lengths of the group 12 dimers. Comparisons are also made to a limited number of all-electron Douglas–Kroll–Hess (DKH) calculations (second and third order) carried out using new correlation consistent basis sets of triple-ζ quality.

Renormalized coupled-cluster methods exploiting left eigenstates of the similarity-transformed Hamiltonian.

Abstract: Completely renormalized (CR) coupled-cluster (CC) approaches, such as CR-CCSD(T), in which one corrects the standard CC singles and doubles (CCSD) energy for the effects of triply (T) and other higher-than-doubly excited clusters [K. Kowalski and P. Piecuch, J. Chem. Phys. 113, 18 (2000)], are reformulated in terms of the left eigenstates ⟨Φ∣L of the similarity-transformed Hamiltonian of CC theory. The resulting CR-CCSD(T)L or CR-CC(2,3) and other CR-CCL methods are derived from the new biorthogonal form of the method of moments of CC equations (MMCC) in which, in analogy to the original MMCC theory, one focuses on the noniterative corrections to standard CC energies that recover the exact, full configuration-interaction energies. One of the advantages of the biorthogonal MMCC theory, which will be further analyzed and extended to excited states in a separate paper, is a rigorous size extensivity of the basic ground-state CR-CCL approximations that result from it, which was slightly violated by the original CR-CCSD(T) and CR-CCSD(TQ) approaches. This includes the CR-CCSD(T)L or CR-CC(2,3) method discussed in this paper, in which one corrects the CCSD energy by the relatively inexpensive noniterative correction due to triples. Test calculations for bond breaking in HF, F2, and H2O indicate that the noniterative CR-CCSD(T)L or CR-CC(2,3) approximation is very competitive with the standard CCSD(T) theory for nondegenerate closed-shell states, while being practically as accurate as the full CC approach with singles, doubles, and triples in the bond-breaking region. Calculations of the activation enthalpy for the thermal isomerizations of cyclopropane involving the trimethylene biradical as a transition state show that the noniterative CR-CCSD(T)L approximation is capable of providing activation enthalpies which perfectly agree with experiment.

Approximate treatment of higher excitations in coupled-cluster theory.

Abstract: The possibilities for the approximate treatment of higher excitations in coupled-cluster (CC) theory are discussed. Potential routes for the generalization of corresponding approximations to lower-level CC methods are analyzed for higher excitations. A general string-based algorithm is presented for the evaluation of the special contractions appearing in the equations specific to those approximate CC models. It is demonstrated that several iterative and noniterative approximations to higher excitations can be efficiently implemented with the aid of our algorithm and that the coding effort is mostly reduced to the generation of the corresponding formulas. The performance of the proposed and implemented methods for total energies is assessed with special regard to quadruple and pentuple excitations. The applicability of our approach is illustrated by benchmark calculations for the butadiene molecule. Our results demonstrate that the proposed algorithm enables us to consider the effect of quadruple excitations for molecular systems consisting of up to 10–12 atoms.

Coupled-cluster methods including noniterative corrections for quadruple excitations

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

Ab initio quantum chemistry: Methodology and applications

Abstract: This Perspective provides an overview of state-of-the-art ab initio quantum chemical methodology and applications. The methods that are discussed include coupled cluster theory, localized second-order Moller–Plesset perturbation theory, multireference perturbation approaches, and density functional theory. The accuracy of each approach for key chemical properties is summarized, and the computational performance is analyzed, emphasizing significant advances in algorithms and implementation over the past decade. Incorporation of a condensed-phase environment by means of mixed quantum mechanical/molecular mechanics or self-consistent reaction field techniques, is presented. A wide range of illustrative applications, focusing on materials science and biology, are discussed briefly.

The extrapolation of one-electron basis sets in electronic structure calculations: how it should work and how it can be made to work.

Abstract: We consider the extrapolation of the one-electron basis to the basis set limit in the context of coupled cluster calculations. We produce extrapolation coefficients that produce much more accurate results than previous extrapolation forms. These are determined by fitting to accurate benchmark results. For coupled cluster singles doubles energies, we take our benchmark results from the work of Klopper that explicitly includes the interelectronic distance. For the perturbative triples energies, our benchmark results are obtained from large even-tempered basis set calculations.

A short-range gradient-corrected density functional in long-range coupled-cluster calculations for rare gas dimers

Abstract: A previously proposed scheme for coupling short-range (sr) density functionals with wavefunction-based long-range (lr) ab initio methods has been extended by (a) developing a new gradient-corrected sr functional of the Perdew–Burke–Ernzerhof (PBE) type and (b) introducing coupled-cluster (CC) approaches (CC with single and double excitations (CCSD), and with additional perturbative triples (CCSD(T))) at the ab initio side. The results show that mixing-in of lr-ab initio correlation helps to remove deficiencies of currently used density functionals for the treatment of van-der-Waals interactions. Compared to full ab initio calculations, the basis set dependence is weaker so that the accuracy of the mixed results surpasses that of the ab initio ones for basis sets of triple-zeta quality.

Coupled-cluster theory based upon the fragment molecular-orbital method.

Abstract: The fragment molecular-orbital (FMO) method was combined with the single-reference coupled-cluster (CC) theory. The developed method (FMO-CC) was applied at the CCSD and CCSD(T) levels of theory, for the cc-pVnZ family of basis sets (n=D,T,Q) to water clusters and glycine oligomers (up to 32 molecules/residues using as large basis sets as possible for the given system). The two- and three-body FMO-CC results are discussed at length, with emphasis on the basis-set dependence and three-body effects. Two- and three-body approximations based on interfragment distances were developed and the values appropriate for their accurate application carefully determined. The error in recovering the correlation energy was several millihartree for the two-body FMO-CC method and in the submillihartree range for the three-body FMO-CC method. In the largest calculations, we were able to perform the CCSD(T) calculations of (H2O)32 with the cc-pVQZ basis set (3680 basis functions) and (GLY)32 with the cc-VDZ basis set (712 co...

Where does the planar-to-nonplanar turnover occur in small gold clusters?

Abstract: Several levels of theory, including both Gaussian-based and plane wave density functional theory (DFT), second-order perturbation theory (MP2), and coupled cluster methods (CCSD(T)), are employed to study Au6 and Au8 clusters. All methods predict that the lowest energy isomer of Au6 is planar. For Au8, both DFT methods predict that the two lowest isomers are planar. In contrast, both MP2 and CCSD(T) predict the lowest Au8 isomers to be nonplanar.

An exponential multireference wave-function Ansatz.

Abstract: An exponential multireference wave-function Ansatz is formulated. In accordance with the state universal coupled-cluster Ansatz of Jeziorski and Monkhorst [Phys. Rev. A 24, 1668 (1981)] the approach uses a reference specific cluster operator. In order to achieve state selectiveness the excitation- and reference-related amplitude indexing of the state universal Ansatz is replaced by an indexing which is based on excited determinants. There is no reference determinant playing a particular role. The approach is size consistent, coincides with traditional single-reference coupled cluster if applied to a single-reference, and converges to full configuration interaction with an increasing cluster operator excitation level. Initial applications on BeH2, CH2, Li2, and nH2 are reported.

Thermodynamic Properties of Molecular Borane Amines and the [BH4-][NH4+] Salt for Chemical Hydrogen Storage Systems from ab Initio Electronic Structure Theory

Abstract: The heats of formation for the borane amines BH3NH3, BH2NH2, and HBNH, tetrahedral BH4-, and the BN molecule have been calculated by using ab initio molecular orbital theory. Coupled cluster calculations with single and double excitations and perturbative triples (CCSD(T)) were employed for the total valence electronic energies. Correlation consistent basis sets were used, up through the augmented quadruple-ζ, to extrapolate to the complete basis set limit. Core/valence, scalar relativistic, and spin−orbit corrections were included in an additive fashion to predict the atomization energies. Geometries were calculated at the CCSD(T) level up through at least aug-cc-pVTZ and frequencies were calculated at the CCSD(T)/aug-cc-pVDZ level. The heats of formation (in kcal/mol) at 0 K in the gas phase are ΔHf(BH3NH3) = −9.1, ΔHf(BH2NH2) = −15.9, ΔHf(BHNH) = 13.6, ΔHf(BN) = 146.4, and ΔHf(BH4-) = −11.6. The reported experimental value for ΔHf(BN) is clearly in error. The heat of formation of the salt [BH4-][NH4+](...

Coupled-cluster theory with simplified linear-r12 corrections: The CCSD(R12) model

Abstract: A simplified singles-and-doubles linear-r12 corrected coupled-cluster model, denoted CCSD(R12), is proposed and compared with the complete singles-and-doubles linear-r12 coupled-cluster method CCSD-R12. An orthonormal auxiliary basis set is used for the resolution-of-the-identity approximation to calculate three-electron integrals needed in the linear-r12 Ansatz. Basis-set convergence is investigated for a selected set of atoms and small molecules. In a large basis, the CCSD(R12) model provides an excellent approximation to the full linear-r12 energy contribution, whereas the magnitude of this contribution is significantly overestimated at the level of second-order perturbation theory.

How Well Can Density Functional Methods Describe Hydrogen Bonds to π Acceptors

Abstract: We employed four newly developed density functional theory (DFT) methods for the calculation of five π hydrogen bonding systems, namely, H2O−C6H6, NH3−C6H6, HCl−C6H6, H2O−indole, and H2O−methylindole. We report new coupled cluster calculations for HCl−C6H6 that support the experimental results of Gotch and Zwier. Using the best available theoretical and experimental results for all five systems, our calculations show that the recently proposed MPW1B95, MPWB1K, PW6B95, and PWB6K methods give accurate energetic and geometrical predictions for π hydrogen bonding interactions, for which B3LYP fails and PW91 is less accurate. We recommend the most recent DFT method, PWB6K, for investigating larger π hydrogen bonded systems, such as those that occur in molecular recognition, protein folding, and crystal packing.

Molecular equilibrium geometries based on coupled-cluster calculations including quadruple excitations

Abstract: Using analytic gradient techniques and an additivity scheme for the various electron correlation contributions, i.e. core-correlation, contribution due to full treatment of triple excitations and contributions due to quadruple excitations calculated with different basis sets, the accuracy of computed geometrical parameters are analysed in comparison with experiment. For a test set of 12 closed-shell and 5 open-shell molecules, it is found that inclusion of quadruple excitations is essential to reach agreement with experiment. The mean error of 0.002 pm and the standard deviation of 0.040 pm of the present CCSD(T)/cc-pV6Z + core(CCSD(T)/cc-pCVQZ) + T/cc-pVTZ + Q/cc-pVDZ results for the closed-shell systems underline the importance of quadruple excitations, in particular, as corresponding calculations without quadruple excitations exhibit significantly larger error. Quadruples contributions for multiply bonded systems as well as the challenging F2 molecule are as large as 0.1 to 0.3 pm, while for single bon...

Basis set and electron correlation effects on the polarizability and second hyperpolarizability of model open-shell π-conjugated systems

Abstract: The basis set and electron correlation effects on the static polarizability (alpha) and second hyperpolarizability (gamma) are investigated ab initio for two model open-shell pi-conjugated systems, the C(5)H(7) radical and the C(6)H(8) radical cation in their doublet state. Basis set investigations evidence that the linear and nonlinear responses of the radical cation necessitate the use of a less extended basis set than its neutral analog. Indeed, double-zeta-type basis sets supplemented by a set of d polarization functions but no diffuse functions already provide accurate (hyper)polarizabilities for C(6)H(8) whereas diffuse functions are compulsory for C(5)H(7), in particular, p diffuse functions. In addition to the 6-31G(*)+pd basis set, basis sets resulting from removing not necessary diffuse functions from the augmented correlation consistent polarized valence double zeta basis set have been shown to provide (hyper)polarizability values of similar quality as more extended basis sets such as augmented correlation consistent polarized valence triple zeta and doubly augmented correlation consistent polarized valence double zeta. Using the selected atomic basis sets, the (hyper)polarizabilities of these two model compounds are calculated at different levels of approximation in order to assess the impact of including electron correlation. As a function of the method of calculation antiparallel and parallel variations have been demonstrated for alpha and gamma of the two model compounds, respectively. For the polarizability, the unrestricted Hartree-Fock and unrestricted second-order Moller-Plesset methods bracket the reference value obtained at the unrestricted coupled cluster singles and doubles with a perturbative inclusion of the triples level whereas the projected unrestricted second-order Moller-Plesset results are in much closer agreement with the unrestricted coupled cluster singles and doubles with a perturbative inclusion of the triples values than the projected unrestricted Hartree-Fock results. Moreover, the differences between the restricted open-shell Hartree-Fock and restricted open-shell second-order Moller-Plesset methods are small. In what concerns the second hyperpolarizability, the unrestricted Hartree-Fock and unrestricted second-order Moller-Plesset values remain of similar quality while using spin-projected schemes fails for the charged system but performs nicely for the neutral one. The restricted open-shell schemes, and especially the restricted open-shell second-order Moller-Plesset method, provide for both compounds gamma values close to the results obtained at the unrestricted coupled cluster level including singles and doubles with a perturbative inclusion of the triples. Thus, to obtain well-converged alpha and gamma values at low-order electron correlation levels, the removal of spin contamination is a necessary but not a sufficient condition. Density-functional theory calculations of alpha and gamma have also been carried out using several exchange-correlation functionals. Those employing hybrid exchange-correlation functionals have been shown to reproduce fairly well the reference coupled cluster polarizability and second hyperpolarizability values. In addition, inclusion of Hartree-Fock exchange is of major importance for determining accurate polarizability whereas for the second hyperpolarizability the gradient corrections are large.

Coupled-cluster method tailored by configuration interaction.

Abstract: A method is presented which combines coupled cluster (CC) and configuration interaction (CI) to describe accurately potential-energy surfaces (PESs). We use the cluster amplitudes extracted from the complete active space CI calculation to manipulate nondynamic correlation to tailor a single reference CC theory (TCC). The dynamic correlation is then incorporated through the framework of the CC method. We illustrate the method by describing the PESs for HF, H2O, and N2 molecules which involve single, double, and triple bond-breaking processes. To the dissociation limit, this approach yields far more accurate PESs than those obtained from the conventional CC method and the additional computational cost is negligible compared with the CC calculation steps. We anticipate that TCC offers an effective and generally applicable approach for many problems.

On geometries of stacked and H-bonded nucleic acid base pairs determined at various DFT, MP2, and CCSD(T) levels up to the CCSD(T)/complete basis set limit level.

Abstract: The geometries and interaction energies of stacked and hydrogen-bonded uracil dimers and a stacked adeninecdots, three dots, centeredthymine pair were studied by means of high-level quantum chemical calculations. Specifically, standard as well as counterpoise-corrected optimizations were performed at second-order Moller-Plesset (MP2) and coupled cluster level of theory with single, double, and perturbative triple excitations [CCSD(T)] levels with various basis sets up to the complete basis set limit. The results can be summarized as follows: (i) standard geometry optimization with small basis set (e.g., 6-31G(*)) provides fairly reasonable intermolecular separation; (ii) geometry optimization with extended basis sets at the MP2 level underestimates the intermolecular distances compared to the reference CCSD(T) results, whereas the MP2/cc-pVTZ counterpoise-corrected optimization agrees well with the reference geometries and, therefore, is recommended as a next step for improving MP2/cc-pVTZ geometries; (iii) the stabilization energy of stacked nucleic acids base pairs depends considerably on the method used for geometry optimization, so the use of reliable geometries, such as counterpoise-corrected MP2/cc-pVTZ ones, is recommended; (iv) the density functional theory methods fail completely in locating the energy minima for stacked structures and when the geometries from MP2 calculations are used, the resulting stabilization energies are strongly underestimated; (v) the self-consistent charges-density functional tight binding method, with inclusion of the empirical dispersion energy, accurately reproduces interaction energies and geometries of dispersion-bonded (stacked) complexes; this method can thus be recommended for prescanning the potential energy surfaces of van der Waals complexes.

Extension of renormalized coupled-cluster methods including triple excitations to excited electronic states of open-shell molecules.

Abstract: The general-purpose open-shell implementation of the completely renormalized equation-of-motion coupled-cluster approach with singles, doubles, and noniterative triples [CR-EOMCCSD(T)] is reported. Benchmark calculations for the low-lying doublet and quartet states of the CH radical show that the CR-EOMCCSD(T) method is capable of providing a highly accurate description of ground and excited states of open-shell molecules. This includes states with strong double excitation character, for which the conventional EOMCCSD approach fails.

Frozen Natural Orbitals: Systematic Basis Set Truncation for Coupled-Cluster Theory

Abstract: The method of frozen natural orbital (FNO) basis set truncation for coupled-cluster theory is described. Numerical comparisons of the FNO potential energy surfaces of a group of small molecules at the CCSD(T) level in DZP, cc-pVTZ, cc-pVQZ bases show that truncation of up to 50% of the virtual space yields CC correlation energies that are accurate to 90 or 95% when added to the full MBPT(2) basis result. The FNO truncation method is also applied to dimethylnitramine (DMNA): both the equilibrium structure and dimer interactions, yielding results at the CCSD(T) level in both a DZP and cc-pVTZ basis set that agree with literature values. CCSD(T) calculations at two possible equilibrium structures of 1,3,5-trinitrohexahydro-1,3,5-triazine (RDX) in a truncated DZP basis are also reported.

Spin-conserving and spin-flipping equation-of-motion coupled-cluster method with triple excitations.

Abstract: We report the implementation of the spin-conserving and spin-flipping variants of the equation-of-motion (EOM) coupled-cluster (CC) model, which includes single and double excitations in the CC part and single, double, and triple excitations in the EOM part, ie, EOM-CC(2,3) [Hirata, Nooijen, Bartlett, Chem Phys Lett 326, 255 (2000)] for closed- and open-shell references Inclusion of triples significantly improves the accuracy of EOM-CCSD for excitation energies (EOM-EE-CCSD) and its spin-flip (SF) counterpart, EOM-SF-CCSD, especially when the reference wave function is strongly spin-contaminated A less computationally demanding active space variant with semi-internal triples has also been implemented The capabilities of full and active space EOM-CC(2,3) are demonstrated by applications to CO+ and CH radicals as well as to the methylene and trimethylenemethane diradicals and the dehydro-m-xylylene triradical

Molecular ionization energies and ground- and ionic-state properties using a non-Dyson electron propagator approach

Abstract: An earlier proposed propagator method for the treatment of molecular ionization is tested in first applications. The method referred to as the non-Dyson third-order algebraic-diagrammatic construction [nD-ADC(3)] approximation for the electron propagator represents a computationally promising alternative to the existing Dyson ADC(3) method. The advantage of the nD-ADC(3) scheme is that the (N+/-1)-electronic parts of the one-particle Green's function are decoupled from each other and the corresponding equations can be solved separately. For a test of the method the nD-ADC(3) results for the vertical ionization transitions in C(2)H(4), CO, CS, F(2), H(2)CO, H(2)O, HF, N(2), and Ne are compared with available experimental and theoretical data including results of full configuration interaction (FCI) and coupled cluster computations. The mean error of the nD-ADC(3) ionization energies relative to the experimental and FCI results is about 0.2 eV. The nD-ADC(3) method, scaling as n(5) with the number of orbitals, requires the solution of a relatively simple Hermitian eigenvalue problem. The method renders access to ground-state properties such as dipole moments. Moreover, also one-electron properties of (N+/-1) electron states can now be studied as a consequence of a specific intermediate-state representation (ISR) formulation of the nD-ADC approach. Corresponding second-order ISR equations are presented.

Coupled Cluster Calculation of the n → π* Electronic Transition of Acetone in Aqueous Solution

Abstract: The combined linear response coupled cluster/molecular mechanics (CC/MM) scheme including mutual polarization effects in the coupling Hamiltonian is applied together with supermolecular CC methods to the study of the gas-to-aqueous solution blue shift of the n → π* excitation energy in acetone. The aug-cc-pVDZ basis set is found to be adequate for the calculation of this excitation energy. In the condensed phase, the shift in the excitation energy is obtained by statistical averaging over 800 solute−solvent configurations extracted from a molecular dynamics simulation. We find the shift to be around 1100−1200 cm-1 depending on the specific model used to describe solvent polarization. The importance of including explicit polarization in both the molecular dynamics simulation as well as the CC/MM calculations is emphasized. Furthermore, the significant dependence of the excitation energy on the CO bond length of acetone is discussed.

The XΣg+1, BΔg1, and B′Σg+1 states of C2: A comparison of renormalized coupled-cluster and multireference methods with full configuration interaction benchmarks

Abstract: Unusual bonding and electronic near degeneracies make the lowest-lying singlet states of the C2 molecule particularly challenging for electronic structure theory. Here we compare two alternative approaches to modeling bond-breaking reactions and excited states: sophisticated multireference configuration interaction and multireference perturbation theory methods, and a more “black box,” single-reference approach, the completely renormalized coupled-cluster method. These approximate methods are assessed in light of their ability to reproduce the full configuration interaction potential energy curves for the XΣg+1, BΔg1, and B′Σg+1 states of C2, which are numerically exact solutions of the electronic Schrodinger equation within the space spanned by a 6-31G* basis set. Both the multireference methods and the completely renormalized coupled-cluster approach provide dramatic improvements over the standard single-reference methods. The multireference methods are nearly as reliable for this challenging test case ...

Extensive generalization of renormalized coupled-cluster methods.

Abstract: The recently developed completely renormalized (CR) coupled-cluster (CC) methods with singles, doubles, and noniterative triples or triples and quadruples [CR-CCSD(T) or CR-CCSD(TQ), respectively], which are based on the method of moments of CC equations (MMCC) [K. Kowalski and P. Piecuch, J. Chem. Phys. 113, 18 (2000)], eliminate the failures of the standard CCSD(T) and CCSD(TQ) methods at larger internuclear separations, but they are not rigorously size extensive. Although the departure from strict size extensivity of the CR-CCSD(T) and CR-CCSD(TQ) methods is small, it is important to examine the possibility of formulating the improved CR-CC methods, which are as effective in breaking chemical bonds as the existing CR-CCSD(T) and CR-CCSD(TQ) approaches, which are as easy to use as the CR-CCSD(T) and CR-CCSD(TQ) methods, and which can be made rigorously size extensive. This may be particularly useful for the applications of CR-CC methods and other MMCC approaches in calculations of potential energy surfaces of large many-electron systems and van der Waals molecules, where the additive separability of energies in the noninteracting limit is very important. In this paper, we propose different types of CR-CC approximations, termed the locally renormalized (LR) CCSD(T) and CCSD(TQ) methods, which become rigorously size extensive if the orbitals are localized on nointeracting fragments. The LR-CCSD(T) and LR-CCSD(TQ) methods rely on the form of the energy expression in terms of the generalized moments of CC equations, derived in this work, termed the numerator-denominator-connected MMCC expansion. The size extensivity and excellent performance of the LR-CCSD(T) and LR-CCSD(TQ) methods are illustrated numerically by showing the results for the dimers of stretched HF and LiH molecules and bond breaking in HF and H2O.

Active-space equation-of-motion coupled-cluster methods for excited states of radicals and other open-shell systems: EA-EOMCCSDt and IP-EOMCCSDt

Abstract: The previously developed active-space coupled-cluster (CC) and equation-of-motion (EOM) CC methods are extended to radicals and other open-shell systems by combining them with the electron attached (EA) and ionized (IP) EOMCC approaches. As illustrated by the calculations for the CH and OH radicals, the resulting EA-EOMCCSDt and IP-EOMCCSDt theories are capable of providing a highly accurate description of the electronic spectra of radical systems, including excited states displaying a manifestly multideterminantal nature, with the low costs that are not much greater that those characterizing the standard EOMCC singles and doubles method.

Ab initio calculations of structures and interaction energies of toluene dimers including CCSD(T) level electron correlation correction.

Abstract: The intermolecular interaction energy of the toluene dimer has been calculated with the ARS-F model (a model chemistry for the evaluation of intermolecular interaction energy between ARomatic Systems using Feller’s method), which was formerly called as the AIMI model III. The CCSD(T) (coupled cluster calculations with single and double substitutions with noniterative triple excitations) interaction energy at the basis set limit has been estimated from the second-order Moller–Plesset perturbation interaction energy at the basis set limit obtained by Feller’s method and the CCSD(T) correction term obtained using a medium-size basis set. The cross (C2) dimer has the largest (most negative) interaction energy (−4.08kcal∕mol). The antiparallel (C2h) and parallel (CS) dimers (−3.77 and −3.41kcal∕mol, respectively) are slightly less stable. The dispersion interaction is found to be the major source of attraction in the toluene dimer. The dispersion interaction mainly determines the relative stability of the stac...