Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEvaluating the GW Approximation with CCSD(T) for Charged Excitations Across the Oligoacenes
117
0
0.0
(
0
)
Added by
Tonatiuh Rangel Gordillo
Publication date
2016
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
Charged excitations of the oligoacene family of molecules, relevant for astrophysics and technological applications, are widely studied and therefore provide an excellent system for benchmarking theoretical methods. In this work, we evaluate the performance of many-body perturbation theory within the GW approximation relative to new high-quality CCSD(T) reference data for charged excitations of the acenes. We compare GW calculations with a number of hybrid density functional theory starting points and with eigenvalue self-consistency. Special focus is given to elucidating the trend of GW-predicted excitations with molecule length increasing from benzene to hexacene. We find that GW calculations with starting points based on an optimally tuned range-separated hybrid (OTRSH) density functional and eigenvalue self-consistency can yield quantitative ionization potentials for the acenes. However, for larger acenes, the predicted electron affinities can deviate considerably from reference values. Our work paves the way for predictive and cost-effective GW calculations of charged excitations of molecules and identifies certain limitations of current GW methods used in practice for larger molecules.
rate research
Read More
We describe a finite-field approach to compute density response functions, which allows for efficient $G_0W_0$ and $G_0W_0Gamma_0$ calculations beyond the random phase approximation. The method is easily applicable to density functional calculations performed with hybrid functionals. We present results for the electronic properties of molecules and solids and we discuss a general scheme to overcome slow convergence of quasiparticle energies obtained from $G_0W_0Gamma_0$ calculations, as a function of the basis set used to represent the dielectric matrix.
We discuss the analytic and diagrammatic structure of ionization potential (IP) and electron affinity (EA) equation-of-motion coupled-cluster (EOM-CC) theory, in order to put it on equal footing with the prevalent $GW$ approximation. The comparison is most straightforward for the time-ordered one-particle Greens function, and we show that the Greens function calculated by EOM-CC with single and double excitations (EOM-CCSD) includes fewer ring diagrams at higher order than does the $GW$ approximation, due to the formers unbalanced treatment of time-ordering. However, the EOM-CCSD Greens function contains a large number of vertex corrections, including ladder diagrams, mixed ring-ladder diagrams, and exchange diagrams. By including triple excitations, the EOM-CCSDT Greens function includes all diagrams contained in the $GW$ approximation, along with many high-order vertex corrections. In the same language, we discuss a number of common approximations to the EOM-CCSD equations, many of which can be classified as elimination of diagrams. Finally, we present numerical results by calculating the principal charged excitations energies of the molecules contained in the so-called $GW$100 test set [J. Chem. Theory Comput. 2015, 11, 5665-5687]. We argue that (in molecules) exchange is as important as screening, advocating for a Hartree-Fock reference and second-order exchange in the self-energy.
Using the simple (symmetric) Hubbard dimer, we analyze some important features of the $GW$ approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot $GW$ method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard $GW$ approximation and find, in particular, that i) some neutral excitation energies become complex when the electron-electron interaction $U$ increases, which can be traced back to the approximate nature of the $GW$ quasiparticle energies; ii) the BSE formalism yields accurate correlation energies over a wide range of $U$ when the trace (or plasmon) formula is employed; iii) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); iv) the trace formula has the correct behavior for weak (ie, small $U$) interaction, unlike the ACFDT expression.
We present a novel hybrid quantum/classical (QM/MM) approach to the calculation of charged excitations in molecular solids based on the many-body Greens function $GW$ formalism. Molecules described at the $GW$ level are embedded into the crystalline environment modeled with an accurate classical polarizable scheme. This allows the calculation of electron addition and removal energies in the bulk and at crystal surfaces where charged excitations are probed in photoelectron experiments. By considering the paradigmatic case of pentacene and perfluoropentacene crystals, we discuss the different contributions from intermolecular interactions to electronic energy levels, distinguishing between polarization, which is accounted for combining quantum and classical polarizabilities, and crystal field effects, that can impact energy levels by up to $pm0.6$ eV. After introducing band dispersion, we achieve quantitative agreement (within 0.2 eV) on the ionization potential and electron affinity measured at pentacene and perfluoropentacene crystal surfaces characterized by standing molecules.
Coupled cluster with single, double, and perturbative triple excitations [CCSD(T)] has been extensively employed as the reference method in benchmarking different quantum chemistry methods. In this work, we test the accuracy of CCSD(T) calculating ground state electric dipole moments at the extrapolated complete basis set (CBS) limit. The calculated dipole moments have been compared to an experimental dataset consisted of diatomic molecules with various kinds of bond natures and spin configurations. As a result, to reach a satisfactory agreement with experimental dipole moments, core-correlations should be included for some molecules. However, even when core-correlations are included, the predicted dipole moment deviates considerably from the experimental values for molecules involving transition metal atoms.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
