No Arabic abstract
We study within the many-body Greens function GW and Bethe-Salpeter approaches the neutral singlet excitations of the zinctetraphenylporphyrin and C70 fullerene donor-acceptor complex. The lowest transition is a charge-transfer excitation between the donor and the acceptor with an energy in excellent agreement with recent constrained density functional theory calculations. Beyond the lowest charge-transfer state, of which the energy can be determined with simple electrostatic models that we validate, the Bethe-Salpeter approach provides the full excitation spectrum. We evidence the existence of hot electron-hole states which are resonant in energy with the lowest donor intramolecular excitation and show an hybrid intramolecular and charge-transfer character, favouring the transition towards charge separation. These findings support the recently proposed scenario for charge separation at donor-acceptor interfaces through delocalized hot charge-transfer states.
We study using the Bethe-Salpeter formalism the excitation energies of the zincbacteriochlorinbacteriochlorin dyad, a paradigmatic photosynthetic complex. In great contrast with standard timedependent density functional theory calculations with (semi)local kernels, charge transfer excitations are correctly located above the intramolecular Q-bands transitions found to be in excellent agreement with experiment. Further, the asymptotic Coulomb behavior towards the true quasiparticle gap for charge transfer excitations at long distance is correctly reproduced, showing that the present scheme allows to study with the same accuracy intramolecular and charge transfer excitations at various spatial range and screening environment without any adjustable parameter.
We study within the many-body Greens function GW and Bethe-Salpeter formalisms the excitation energies of a paradigmatic model dipeptide, focusing on the four lowest-lying local and charge-transfer excitations. Our GW calculations are performed at the self-consistent level, updating first the quasiparticle energies, and further the single-particle wavefunctions within the static Coulomb-hole plus screened-exchange approximation to the GW self-energy operator. Important level crossings, as compared to the starting Kohn-Sham LDA spectrum, are identified. Our final Bethe-Salpeter singlet excitation energies are found to agree, within 0.07 eV, with CASPT2 reference data, except for one charge-transfer state where the discrepancy can be as large as 0.5 eV. Our results agree best with LC-BLYP and CAM-B3LYP calculations with enhanced long-range exchange, with a 0.1 eV mean absolute error. This has been achieved employing a parameter-free formalism applicable to metallic or insulating extended or finite systems.
Helium atom is the simplest many-body electronic system provided by nature. The exact solution to the Schrodinger equation is known for helium ground and excited states, and represents a workbench for any many-body methodology. Here, we check the ab initio many-body GW approximation and Bethe-Salpeter equation (BSE) against the exact solution for helium. Starting from Hartree-Fock, we show that GW and BSE yield impressively accurate results on excitation energies and oscillator strength, systematically improving time-dependent Hartree-Fock. These findings suggest that the accuracy of BSE and GW approximations is not significantly limited by self-interaction and self-screening problems even in this few electron limit. We further discuss our results in comparison to those obtained by time-dependent density-functional theory.
We introduce a new computational method to study porphyrin-like transition metal complexes, bridging density functional theory and exact many-body techniques, such as the density matrix renormalization group (DMRG). We first derive a multi-orbital Anderson impurity Hamiltonian starting from first principles considerations that qualitatively reproduce GGA+U results when ignoring inter-orbital Coulomb repulsion $U$ and Hund exchange $J$. An exact canonical transformation is used to reduce the dimensionality of the problem and make it amenable to DMRG calculations, including all many-body terms (both intra, and inter-orbital), which are treated in a numerically exact way. We apply this technique to FeN$_4$ centers in graphene and show that the inclusion of these terms has dramatic effects: as the iron orbitals become single occupied due to the Coulomb repulsion, the inter-orbital interaction further reduces the occupation yielding a non-monotonic behavior of the magnetic moment as a function of the interactions, with maximum polarization only in a small window at intermediate values of the parameters. Furthermore, $U$ changes the relative position of the peaks in the density of states, particularly on the iron $d_{z^2}$ orbital, which is expected to greatly affect the binding of ligands.
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.