A theoretical model chemistry designed to achieve high accuracy for enthalpies of formation of atoms and small molecules is described. This approach is entirely independent of experimental data and contains no empirical scaling factors, and includes a treatment of electron correlation up to the full coupled-cluster singles, doubles, triples and quadruples approach. Energies are further augmented by anharmonic zero-point vibrational energies, a scalar relativistic correction, first-order spin-orbit coupling, and the diagonal Born-Oppenheimer correction. The accuracy of the approach is assessed by several means. Enthalpies of formation (at 0 K) calculated for a test suite of 31 atoms and molecules via direct calculation of the corresponding elemental formation reactions are within 1 kJ mol(-1) to experiment in all cases. Given the quite different bonding environments in the product and reactant sides of these reactions, the results strongly indicate that even greater accuracy may be expected in reactions that preserve (either exactly or approximately) the number and types of chemical bonds.

Schemes for the analytic calculation of the diagonal Born-Oppenheimer correction (DBOC) are formulated and implemented for use with general single-reference configuration-interaction and coupled-cluster wave function models. Calculations are reported to demonstrate the convergence of the DBOC with respect to electron-correlation treatment and basis set as well as to investigate the size-consistency error in configuration-interaction calculations of the DBOC. The importance of electron-correlation contributions to the DBOC is illustrated in the computation of the corresponding corrections for the reaction energy and activation barrier of the F + H2 --> FH + H reaction as well as of the atomization energy for trans-butadiene.

The numerous existing publications on benchmarking quantum chemistry methods for excited states rarely include Charge Transfer (CT) states, although many interesting phenomena in, e.g., biochemistry and material physics involve the transfer of electrons between fragments of the system. Therefore, it is timely to test the accuracy of quantum chemical methods for CT states, as well. In this study we first propose a new benchmark set consisting of dimers having low-energy CT states. On this set, the vertical excitation energy has been calculated with Coupled Cluster methods including triple excitations (CC3, CCSDT-3, CCSD(T)(a)*), as well as with methods including full or approximate doubles (CCSD, STEOM-CCSD, CC2, ADC(2), EOM-CCSD(2)). The results show that the popular CC2 and ADC(2) methods are much less accurate for CT states than for valence states. On the other hand, EOM-CCSD seems to have similar systematic overestimation of the excitation energies for both types of states. Among the triples methods the novel EOM-CCSD(T)(a)* method including noniterative triple excitations is found to stand out with its consistently good performance for all types of states, delivering essentially EOM-CCSDT quality results.

We present a comprehensive statistical analysis on the accuracy of various excited state Coupled Cluster methods, accentuating the effect of diffuse basis sets on vertical excitation energies of valence and Rydberg-type states. Many popular approximate doubles and triples methods are benchmarked with basis sets up to aug-cc-pVTZ, with high level EOM-CCSDT results used as reference. The results reveal a serious deficiency of CC2 linear response and CIS(D) techniques in the description of Rydberg states, a feature not shown by the EOM-CCSD(2) and EOM-CCSD variants. The CC3 theory proves to be an accurate choice among the iterative approximate triples methods, while the novel perturbation-based CCSD(T)(a)* variant turns out to be the best way to include the effect of triple excitations in a noniterative way.

Theory and implementation for evaluation of the nonadiabatic coupling vector between excited electronic states described by equation-of-motion excitation energy coupled-cluster singles and doubles (EOMEE-CCSD) method is presented. Problems arising from the non-Hermitian nature of the theory are discussed in detail. The performance of the new approach is demonstrated by the nice agreement of the nonadiabatic coupling curves for LiH obtained at the EOMEE-CCSD and MR-CISD levels. Using the tools developed we also present a computational procedure to evaluate the interstate coupling constants used in vibronic coupling theories. As an application of this part of the implementation we present simulation of the electronic absorption spectrum of the pyrazine molecule within the linear vibronic coupling model.

4-(N,N-Dimethylamino)benzonitrile (DMABN) is the archetypal system for dual fluorescence. Several past studies, both experimental and theoretical, have examined the mechanism of its relaxation in the gas phase following photoexcitation to the S2 state, without converging to a single description. In this contribution, we report first-principles simulations of the early events involved in this process performed using the nonadiabatic trajectory surface hopping (TSH) approach in combination with the ADC(2) electronic structure method. ADC(2) is verified to reproduce the ground- and excited-state structures of DMABN in reasonably close agreement with previous theoretical benchmarks. The TSH simulations predict that internal conversion from the S2 state to the S1 takes place as early as 8.5 fs, on average, after the initial photoexcitation, and with no significant torsion of the dimethylamino group relative to the aromatic ring. As evidenced by supporting EOM-CCSD calculations, the population transfer from S2 to S1 can be attributed to the skeletal deformation modes of the aromatic ring and the stretching of the ring-dimethylamino nitrogen bond. The non- or slightly twisted locally excited structure is the predominant product of the internal conversion, and the twisted intramolecular charge transfer structure is formed through equilibration with the locally excited structure with no change of adiabatic state. These findings point toward a new interpretation of data from previous time-resolved experiments.

Conical intersections between electronic states are of great importance for the understanding of radiationless ultrafast relaxation processes. In particular, accidental degeneracies of hypersurfaces, i.e., between states of the same symmetry, become increasingly relevant for larger molecular systems. Coupled-cluster theory, including both single and multireference based schemes, offers a size-extensive description of the electronic wave function, but it sacrifices the Hermitian character of the theory. In this contribution, we examine the consequences of anti-Hermitian contributions to the coupling matrix element between near-degenerate states such as linear dependent eigenvectors and complex eigenvalues. Numerical examples are given for conical intersections between two excited states calculated at the equation-of-motion coupled-cluster level which indeed show the predicted artifacts. A simple method is suggested which allows physically meaningful potential energy surfaces to be extracted from the otherwise ill-behaved results. This provides a perspective for obtaining potential energy surfaces near conical intersections at the coupled-cluster level.

Describing electronically excited states of molecules accurately poses a challenging problem for theoretical methods. Popular second order techniques like Linear Response CC2 (CC2-LR), Partitioned Equation-of-Motion MBPT(2) (P-EOM-MBPT(2)), or Equation-of-Motion CCSD(2) (EOM-CCSD(2)) often produce results that are controversial and are ill-balanced with their accuracy on valence and Rydberg type states. In this study, we connect the theory of these methods and, to investigate the origin of their different behavior, establish a series of intermediate variants. The accuracy of these on excitation energies of singlet valence and Rydberg electronic states is benchmarked on a large sample against high-accuracy Linear Response CC3 references. The results reveal the role of individual terms of the second order similarity transformed Hamiltonian, and the reason for the bad performance of CC2-LR in the description of Rydberg states. We also clarify the importance of the T̂ transformation employed in the CC2 procedure, which is found to be very small for vertical excitation energies.