Striving to define very accurate vertical transition energies, we perform both high-level coupled cluster (CC) calculations (up to CCSDTQP) and selected configuration interaction (sCI) calculations (up to several millions of determinants) for 18 small compounds (water, hydrogen sulfide, ammonia, hydrogen chloride, dinitrogen, carbon monoxide, acetylene, ethylene, formaldehyde, methanimine, thioformaldehyde, acetaldehyde, cyclopropene, diazomethane, formamide, ketene, nitrosomethane, and the smallest streptocyanine). By systematically increasing the order of the CC expansion, the number of determinants in the CI expansion as well as the size of the one-electron basis set, we have been able to reach near full CI (FCI) quality transition energies. These calculations are carried out on CC3/aug-cc-pVTZ geometries, using a series of increasingly large atomic basis sets systematically including diffuse functions. In this way, we define a list of 110 transition energies for states of various characters (valence, ...

A Mountaineering Strategy to Excited States: Highly Accurate Reference Energies and Benchmarks.

Development of exponentially scaling methods has seen great progress in tackling larger systems than previously thought possible. One such technique, full configuration interaction quantum Monte Carlo, is a useful algorithm that allows exact diagonalization through stochastically sampling determinants. The method derives its utility from the information in the matrix elements of the Hamiltonian, along with a stochastic projected wave function, to find the important parts of Hilbert space. However, the stochastic representation of the wave function is not required to search Hilbert space efficiently, and here we describe a highly efficient deterministic method that can achieve chemical accuracy for a wide range of systems, including the difficult Cr2 molecule. We demonstrate for systems like Cr2 that such calculations can be performed in just a few cpu hours which makes it one of the most efficient and accurate methods that can attain chemical accuracy for strongly correlated systems. In addition our method also allows efficient calculation of excited state energies, which we illustrate with benchmark results for the excited states of C2.

/pdf/a-deterministic-alternative-to-the-full-configuration-3xzfnquuyp.pdf

A deterministic alternative to the full configuration interaction quantum Monte Carlo method

Development of exponentially scaling methods has seen great progress in tackling larger systems than previously thought possible. One such technique, full configuration interaction quantum Monte Carlo, is a useful algorithm that allows exact diagonalization through stochastically sampling determinants. The method derives its utility from the information in the matrix elements of the Hamiltonian, along with a stochastic projected wave function, to find the important parts of Hilbert space. However, the stochastic representation of the wave function is not required to search Hilbert space efficiently, and here we describe a highly efficient deterministic method to achieve chemical accuracy for a wide range of systems, including the difficult Cr$_{2}$ dimer. In addition our method also allows efficient calculation of excited state energies, for which we illustrate with benchmark results for the excited states of C$_{2}$.

It is shown both theoretically and numerically that the minimal multireference configuration interaction (CI) approach [Liu, W.; Hoffmann, M. R. Theor. Chem. Acc. 2014, 133, 1481] converges quickly and monotonically from above to full CI by updating the primary, external, and secondary states that describe the respective static, dynamic, and again static components of correlation iteratively, even when starting with a rather poor description of a strongly correlated system. In short, the iterative CI (iCI) is a very effective means toward highly correlated wave functions and, ultimately, full CI.

iCI: Iterative CI toward full CI.

Excited states exhibiting double-excitation character are notoriously difficult to model using conventional single-reference methods, such as adiabatic time-dependent density functional theory (TD-DFT) or equation-of-motion coupled cluster (EOM-CC). In addition, these states are typical experimentally "dark", making their detection in photoabsorption spectra very challenging. Nonetheless, they play a key role in the faithful description of many physical, chemical, and biological processes. In the present work, we provide accurate reference excitation energies for transitions involving a substantial amount of double excitation using a series of increasingly large diffuse-containing atomic basis sets. Our set gathers 20 vertical transitions from 14 small- and medium-size molecules (acrolein, benzene, beryllium atom, butadiene, carbon dimer and trimer, ethylene, formaldehyde, glyoxal, hexatriene, nitrosomethane, nitroxyl, pyrazine, and tetrazine). Depending on the size of the molecule, selected configuration interaction (sCI) and/or multiconfigurational (CASSCF, CASPT2, (X)MS-CASPT2, and NEVPT2) calculations are performed in order to obtain reliable estimates of the vertical transition energies. In addition, coupled cluster approaches including at least contributions from iterative triples (such as CC3, CCSDT, CCSDTQ, and CCSDTQP) are assessed. Our results clearly evidence that the error in CC methods is intimately related to the amount of double-excitation character of the transition. For "pure" double excitations (i.e., for transitions which do not mix with single excitations), the error in CC3 can easily reach 1 eV, while it goes down to a few tenths of an electronvolt for more common transitions (such as in trans-butadiene) involving a significant amount of singles. As expected, CC approaches including quadruples yield highly accurate results for any type of transition. The quality of the excitation energies obtained with multiconfigurational methods is harder to predict. We have found that the overall accuracy of these methods is highly dependent on both the system and the selected active space. The inclusion of the σ and σ* orbitals in the active space, even for transitions involving mostly π and π* orbitals, is mandatory in order to reach high accuracy. A theoretical best estimate (TBE) is reported for each transition. We believe that these reference data will be valuable for future methodological developments aiming at accurately describing double excitations.

Reference Energies for Double Excitations.

Defining accurate and compact trial wavefunctions leading to small statistical and fixed-node errors in quantum Monte Carlo (QMC) calculations is still a challenging problem. Here we propose to mak...

/pdf/using-perturbatively-selected-configuration-interaction-in-ong7rwq1vg.pdf

CaffarelMichel

Papers

Using perturbatively selected configuration interaction in quantum Monte Carlo calculations