0000000000133087
AUTHOR
Poul Jørgensen
Communication: The performance of non-iterative coupled cluster quadruples models
We compare the numerical performance of various non-iterative coupled cluster (CC) quadruples models. The results collectively show how approaches that attempt to correct the CC singles and doubles energy for the combined effect of triple and quadruple excitations all fail at recovering the correlation energy of the full CC singles, doubles, triples, and quadruples (CCSDTQ) model to within sufficient accuracy. Such a level of accuracy is only achieved by models that make corrections to the full CC singles, doubles, and triples (CCSDT) energy for the isolated effect of quadruple excitations of which the CCSDT(Q–3) and CCSDT(Q–4) models of the Lagrangian-based CCSDT(Q–n) perturbation series a…
On the convergence of perturbative coupled cluster triples expansions:Error cancellations in the CCSD(T) model and the importance of amplitude relaxation
Recently, we proposed a novel Lagrangian-based perturbation series-the CCSD(T-n) series-which systematically corrects the coupled cluster singles and doubles (CCSD) energy in orders of the Møller-Plesset fluctuation potential for effects due to triple excitations. In the present study, we report numerical results for the CCSD(T-n) series up through fourth order which show the predicted convergence trend throughout the series towards the energy of its target, the coupled cluster singles, doubles, and triples (CCSDT) model. Since effects due to the relaxation of the CCSD singles and doubles amplitudes enter the CCSD(T-n) series at fourth order (the CCSD(T-4) model), we are able to separate th…
The prediction of molecular equilibrium structures by the standard electronic wave functions
A systematic investigation has been carried out of the accuracy of molecular equilibrium structures of 19 small closed-shell molecules containing first-row atoms as predicted by the following standard electronic ab initio models: Hartree–Fock (HF) theory, Mo/ller–Plesset theory to second, third, and fourth orders (MP2, MP3, and MP4), coupled-cluster singles and doubles (CCSD) theory; CCSD theory with perturbational triples corrections [CCSD(T)], and the configuration-interaction singles and doubles (CISD) model. For all models, calculations were carried out using the correlation-consistent polarized valence double-zeta (cc-pVDZ) basis, the correlation-consistent polarized valence triple-zet…
A Lagrangian framework for deriving triples and quadruples corrections to the CCSD energy.
Using the coupled cluster Lagrangian technique, we have determined perturbative corrections to the coupled cluster singles and doubles (CCSD) energy that converge towards the coupled cluster singles, doubles, and triples (CCSDT) and coupled cluster singles, doubles, triples, and quadruples (CCSDTQ) energies, considering the CCSD state as the unperturbed reference state and the fluctua- tion potential as the perturbation. Since the Lagrangian technique is utilized, the energy corrections satisfy Wigner’s 2n + 1 rule for the cluster amplitudes and the 2n + 2 rule for the Lagrange multi- pliers. The energy corrections define the CCSD perturbation series, CCSD(T–n) and CCSD(TQ–n), which are ter…
Ground and excited state polarizabilities and dipole transition properties of benzene from coupled cluster response theory
Abstract The electronic properties and transition properties have been investigated for the lowest singlet electronic states of benzene using coupled cluster response theory. The polarizabilities have been calculated for the ground state and the 1 1 B 2u , 1 1 B 1u 1 1 E 1u and 2 1 E 1u excited states. The dipole allowed transitions out of these states have also been calculated and discussed in the context of the calculated polarizabilities. Oscillator strengths and the second electronic moments of the charge distributions have been used to characterize and identify qualitative features of the individual states. The performance of coupled cluster singles (CCS), the recently proposed CC2 mod…
CCSDT calculations of molecular equilibrium geometries
Abstract CCSDT equilibrium geometries of CO, CH2, F2, HF, H2O and N2 have been calculated using the correlation-consistent cc-pVXZ basis sets. Similar calculations have been performed for SCF, CCSD and CCSD(T). In general, bond lengths decrease when improving the basis set and increase when improving the N-electron treatment. CCSD(T) provides an excellent approximation to CCSDT for bond lengths as the largest difference between CCSDT and CCSD(T) is 0.06 pm. At the CCSDT/cc-pVQZ level, basis set deficiencies, neglect of higher-order excitations, and incomplete treatment of core-correlation all give rise to errors of a few tenths of a pm, but to a large extent, these errors cancel. The CCSDT/…
Comparison of full-configuration interaction and coupled-cluster harmonic and fundamental frequencies for BH and HF
The harmonic and fundamental frequencies are calculated for the potential-energy curves of BH and HF using the full-configuration interaction model and two hierarchies of coupled-cluster wavefunction models. The anharmonic contributions are also obtained using second-order vibrational perturbation theory. A consistent and systematic improvement is seen for both the harmonic and anharmonic contributions when increasing the level of the correlation treatment. The changes are largest for the harmonic contributions. This is also the case when including valence or diffuse functions in the basis set. Second-order perturbation theory gives a good approximation to the anharmonic contribution and in…
Triple excitation effects in coupled cluster calculations of Verdet constants
Abstract The CC3 approach has been employed to calculate the Verdet constants of N 2 ,C 2 H 2 , and CH 4 . For N 2 and C 2 H 2 , relatively large triples contributions are found which need to be included in order to reach close agreement with the experimental constants.
Equation-of-motion coupled cluster perturbation theory revisited
The equation-of-motion coupled cluster (EOM-CC) framework has been used for deriving a novel series of perturbative corrections to the coupled cluster singles and doubles energy that formally con- verges towards the full configuration interaction energy limit. The series is based on a Møller-Plesset partitioning of the Hamiltonian and thus size extensive at any order in the perturbation, thereby rem- edying the major deficiency inherent to previous perturbation series based on the EOM-CC ansatz. © 2014 AIP Publishing LLC. [http://dx.doi.org/10.1063/1.4873138]
Assessment of the accuracy of coupled cluster perturbation theory for open-shell systems. I. Triples expansions
The accuracy at which total energies of open-shell atoms and organic radicals may be calculated is assessed for selected coupled cluster perturbative triples expansions, all of which augment the coupled cluster singles and doubles (CCSD) energy by a non-iterative correction for the effect of triple excitations. Namely, the second- through sixth-order models of the recently proposed CCSD(T-n) triples series [J. Chem. Phys. 140, 064108 (2014)] are compared to the acclaimed CCSD(T) model for both unrestricted as well as restricted open-shell Hartree-Fock (UHF/ROHF) reference determinants. By comparing UHF- and ROHF-based statistical results for a test set of 18 modest-sized open-shell species …
The Dalton quantum chemistry program system
Dalton is a powerful general-purpose program system for the study of molecular electronic structure at the Hartree-Fock, Kohn-Sham, multiconfigurational self-consistent-field, MOller-Plesset, confi ...
The CC3 model : An iterative coupled cluster approach including connected triples
An alternative derivation of many-body perturbation theory (MBPT) has been given, where a coupled cluster parametrization is used for the wave function and the method of undetermined Lagrange multipliers is applied to set up a variational coupled cluster energy expression. In this variational formulation, the nth-order amplitudes determine the energy to order 2n+1 and the nth-order multipliers determine the energy to order 2n+2. We have developed an iterative approximate coupled cluster singles, doubles, and triples model CC3, where the triples amplitudes are correct through second order and the singles amplitudes are treated without approximations due to the unique role of singles as appro…
Convergence of coupled cluster perturbation theory.
The convergence of a recently proposed coupled cluster (CC) family of perturbation series [Eriksen, J. J. et al., J. Chem. Phys. 140, 064108 (2014)], in which the energetic difference between two CC models - a low-level parent and a high-level target model - is expanded in orders of the M{\o}ller-Plesset (MP) fluctuation potential, is investigated for four prototypical closed-shell systems (Ne, singlet methylene, distorted HF, and the fluoride anion) in standard and augmented basis sets. In these investigations, energy corrections of the various series have been calculated to high orders and their convergence radii determined by probing for possible front- and back-door intruder states, the…
The Cotton-Mouton effect of Neon and Argon: a benchmark study using highly correlated coupled cluster wave functions
The Cotton-Mouton effect (magnetic field induced linear birefringence) has been studied for neon and argon using state-of-the-art coupled cluster techniques. The coupled cluster singles, doubles and triples (CCSDT) approach has been used to obtain static benchmark results and the CC3 model with an approximate treatment of triple excitations to obtain frequency-dependent results. In the case of neon the effect of excitations beyond triples has also been estimated via coupled cluster calculations including quadruple excitations (CCSDTQ), pentuple excitations (CCSDTQP), etc. up to the full configuration-interaction level. The results obtained for the anisotropy of the hypermagnetizability Delt…
Assessment of the accuracy of coupled cluster perturbation theory for open-shell systems. II. Quadruples expansions
We extend our assessment of the potential of perturbative coupled cluster (CC) expansions for a test set of open-shell atoms and organic radicals to the description of quadruple excitations. Namely, the second- through sixth-order models of the recently proposed CCSDT(Q-n) quadruples series [J. Chem. Phys. 140, 064108 (2014)] are compared to the prominent CCSDT(Q) and lambda-CCSDT(Q) models. From a comparison of the models in terms of their recovery of total CC singles, doubles, triples, and quadruples (CCSDTQ) energies, we find that the performance of the CCSDT(Q-n) models is independent of the reference used (unrestricted or restricted (open-shell) Hartree-Fock), in contrast to the CCSDT(…
Coupled-cluster singles, doubles and triples (CCSDT) calculations of atomization energies
Atomization energies have been calculated for CO, H2O, F-2, HF, N-2 and CH2 (the (1)A(1) state) using the coupled-duster singles, doubles and triples (CCSDT) model as well as the coupled-cluster singles and doubles model with a perturbative correction for triples [CCSD(T)]. TheCCSD(T) model provides an excellent approximation to the CCSDT model; at the cc-pV5Z basis set level, the CCSDT valence triples contribution is underestimated by 9.1% (0.8 kJ/mol) for CH, and overestimated for the remaining molecules by as little as 4.3%(1.3 kJ/mol) for F-2,and as much as 8.4% (3.0 kJ/mol) for N-2. At the CCSDT level, the agreement with experiment is not improved, suggesting that some cancellation of …
Large-scale calculations of excitation energies in coupled cluster theory: The singlet excited states of benzene
Algorithms for calculating singlet excitation energies in the coupled cluster singles and doubles (CCSD) model are discussed and an implementation of an atomic-integral direct algorithm is presented. Each excitation energy is calculated at a cost comparable to that of the CCSD ground-state energy. Singlet excitation energies are calculated for benzene using up to 432 basis functions. Basis-set effects of the order of 0.2 eV are observed when the basis is increased from augmented polarized valence double-zeta (aug-cc-pVDZ) to augmented polarized valence triple-zeta (aug-cc-pVTZ) quality. The correlation problem is examined by performing calculations in the hierarchy of coupled cluster models…
The electronic spectrum of pyrrole
The electronic spectrum of pyrrole has been investigated by performing calculations using a hierarchy of coupled-cluster models consisting of CCS, CC2, CCSD, and CC3. Basis-set effects have been investigated by carrying out calculations using correlation-consistent basis sets augmented with functions especially designed for the description of Rydberg states. Oscillator strengths, excited state dipole moments, and second moments of the electronic charge distributions have been used to characterize the electronic transitions and final states. Structures and vibrational frequencies have been calculated for a few selected states, and the importance of distinguishing between vertical and adiabat…
Restricted and complete-active-space multiconfiguration linear response calculations of the polarizability of formamide and urea
Abstract Using the polarized basis sets of Sadlej, we have carried out multiconfiguration linear response (MCLR) calculations of static and dynamic polarizabilities of water, carbon dioxide, formamide and urea. It is found that the polarized basis sets give a good description of the polarizabilities. The uncorrelates (self-consistent field) polarizabilities are in general 10% or more lower than the experimental values. The correlation as introduced in the complete-active-space (CAS) and restricted-active-space (RAS) MCLR calculations recovers the major part of this deviation.
Accurate Nonlinear Optical Properties for Small Molecules
During the last decade it became possible to calculate by quantum chemical ab initio methods not only static but also frequency-dependent properties with high accuracy. Today, the most important tools for such calculations are coupled cluster response methods in combination with systematic hierarchies of correlation consistent basis sets. Coupled cluster response methods combine a computationally efficient treatment of electron correlation with a qualitatively correct pole structure and frequency dispersion of the response functions. Both are improved systematically within a hierarchy of coupled cluster models. The present contribution reviews recent advances in the highly accurate calculat…
The accuracy of molecular dipole moments in standard electronic structure calculations
Abstract A systematic investigation has been carried out of the accuracy of calculated molecular equilibrium dipole moments of 11 polar closed-shell molecules, using the HF, MP2, CCSD and CCSD(T) models and correlation-consistent basis sets. Augmented basis sets are important for improving the basis-set convergence, but the quality of the results depends more on the correlation treatment than on the cardinal number of the basis set. Augmented triple-zeta basis sets are sufficient for most calculations. The mean absolute error of the HF calculations is 0.16 D, which is reduced at the MP2 and CCSD levels to 0.048 and 0.025 D, respectively. The CCSD(T) errors are small – typically
Molecular equilibrium structures from experimental rotational constants and calculated vibration–rotation interaction constants
A detailed study is carried out of the accuracy of molecular equilibrium geometries obtained from least-squares fits involving experimental rotational constants B(0) and sums of ab initio vibration-rotation interaction constants alpha(r)(B). The vibration-rotation interaction constants have been calculated for 18 single-configuration dominated molecules containing hydrogen and first-row atoms at various standard levels of ab initio theory. Comparisons with the experimental data and tests for the internal consistency of the calculations show that the equilibrium structures generated using Hartree-Fock vibration-rotation interaction constants have an accuracy similar to that obtained by a dir…
The equilibrium structure of trans-glyoxal from experimental rotational constants and calculated vibration–rotation interaction constants
A total of six high-resolution FT-IR spectra for trans-glyoxal-d2, trans-glyoxal-d1 and trans-glyoxal-13C2 were recorded with a resolution ranging from 0.003 to 0.004 cm−1. By means of a simultaneous ground state combination difference analysis for each of these isotopologues using the Watson Hamiltonian in A-reduction and Ir-representation the ground state rotational constants are obtained. An empirical equilibrium structure is determined for trans-glyoxal using these experimental ground state rotational constants and vibration–rotation interaction constants calculated at the CCSD(T)/cc-pVTZ level of theory. The least-squares fit yields the following structural parameters for trans-glyoxal…