No Arabic abstract
We propose a simple many-body based screening mixing strategy to considerably enhance the performance of the Bethe-Salpeter (BS) approach for prediction of excitation energies of molecular systems. This strategy enables us to nearly reproduce results of highly correlated equation of motion coupled cluster singles and doubles (EOM-CCSD) through optimal use of cancellation effects.
Inspired by Grimmes simplified Tamm-Dancoff density functional theory approach [S. Grimme, J. Chem. Phys. textbf{138}, 244104 (2013)], we describe a simplified approach to excited state calculations within the GW approximation to the self-energy and the Bethe-Salpeter equation (BSE), which we call sGW/sBSE. The primary simplification to the electron repulsion integrals yields the same structure as with tensor hypercontraction, such that our method has a storage requirement that grows quadratically with system size and computational timing that grows cubically with system size. The performance of sGW is tested on the ionization potential of the molecules in the GW100 test set, for which it differs from textit{ab intio} GW calculations by only 0.2 eV. The performance of sBSE (based on sGW input) is tested on the excitation energies of molecules in the Thiel set, for which it differs from textit{ab intio} GW/BSE calculations by about 0.5 eV. As examples of the systems that can be routinely studied with sGW/sBSE, we calculate the band gap and excitation energy of hydrogen-passivated silicon nanocrystals with up to 2650 electrons in 4678 spatial orbitals and the absorption spectra of two large organic dye molecules with hundreds of atoms.
We show that the origin of electronic transitions of molecular many-body systems can be revealed by a quantified natural transition orbitals (QNTO) analysis and the electronic excitations of the total system can be mapped onto a standard orbitals set of a reference system. We further illustrate QNTO on molecular systems by studying the origin of electronic transitions of DNA moiety, thymine and thymidine. This QNTO analysis also allows us to assess the performance of various functionals used in time-dependent density functional response theory.
We check the ab initio GW approximation and Bethe-Salpeter equation (BSE) many-body methodology against the exact solution benchmark of the hydrogen molecule H$_2$ ground state and excitation spectrum, and in comparison with the configuration interaction (CI) and time-dependent Hartree-Fock methods. The comparison is made on all the states we could unambiguously identify from the excitonic wave functions symmetry. At the equilibrium distance $R = 1.4 , a_0$, the GW+BSE energy levels are in good agreement with the exact results, with an accuracy of 0.1~0.2 eV. GW+BSE potential-energy curves are also in good agreement with the CI and the exact result up to $2.3 , a_0$. The solution no longer exists beyond $3.0 , a_0$ for triplets ($4.3 , a_0$ for singlets) due to instability of the ground state. We tried to improve the GW reference ground state by a renormalized random-phase approximation (r-RPA), but this did not solve the problem.
We present a Flexible Ansatz for N-body Configuration Interaction (FANCI) that includes any multideterminant wavefunction. This ansatz is a generalization of the Configuration Interaction (CI) wavefunction, where the coefficients are replaced by a specified function of certain parameters. By making an appropriate choice for this function, we can reproduce popular wavefunction structures like CI, Coupled-Cluster, Tensor Network States, and geminal-product wavefunctions. The universality of this framework suggests a programming structure that allows for the easy construction and optimization of arbitrary wavefunctions. Here, we will discuss the structures of the FANCI framework and its implications for wavefunction properties, particularly accuracy, cost, and size-consistency. We demonstrate the flexibility of this framework by reconstructing popular wavefunction ans{a}tze and modifying them to construct novel wavefunction forms. FANCI provides a powerful framework for exploring, developing, and testing new wavefunction forms.
We describe our efforts of the past few years to create a large set of more than 500 highly-accurate vertical excitation energies of various natures ($pi to pi^*$, $n to pi^*$, double excitation, Rydberg, singlet, doublet, triplet, etc) in small- and medium-sized molecules. These values have been obtained using an incremental strategy which consists in combining high-order coupled cluster and selected configuration interaction calculations using increasingly large diffuse basis sets in order to reach high accuracy. One of the key aspect of the so-called QUEST database of vertical excitations is that it does not rely on any experimental values, avoiding potential biases inherently linked to experiments and facilitating theoretical cross comparisons. Following this composite protocol, we have been able to produce theoretical best estimate (TBEs) with the aug-cc-pVTZ basis set for each of these transitions, as well as basis set corrected TBEs (i.e., near the complete basis set limit) for some of them. The TBEs/aug-cc-pVTZ have been employed to benchmark a large number of (lower-order) wave function methods such as CIS(D), ADC(2), CC2, STEOM-CCSD, CCSD, CCSDR(3), CCSDT-3, ADC(3), CC3, NEVPT2, and others (including spin-scaled variants). In order to gather the huge amount of data produced during the QUEST project, we have created a website [https://lcpq.github.io/QUESTDB_website] where one can easily test and compare the accuracy of a given method with respect to various variables such as the molecule size or its family, the nature of the excited states, the type of basis set, etc. We hope that the present review will provide a useful summary of our effort so far and foster new developments around excited-state methods.