Many-body Green's function GW and Bethe-Salpeter study of the optical excitations in a paradigmatic model dipeptide

A Guide to Molecular Properties from the Bethe-Salpeter Equation

The Bethe-Salpeter equation (BSE) combined with the Green's function GW method has been successfully transformed into a robust computational tool to describe light-matter interactions and excitation spectra for molecules, solids, and materials from first principles. Due to its ability to accurately describe charge transfer and Rydberg excitations, GW-BSE is already an established and cost-efficient alternative to time-dependent density functional theory. This raises the question whether the GW-BSE approach can become a more general framework for molecular properties beyond excitation energies. In this Mini-Review, we recapitulate recent endeavors along this point in terms of both theoretical and practical developments for quantum chemistry, physical chemistry, and related fields. In doing so, we provide guidelines for current applications to chemical challenges in collaboration with experimentalists as well as to future developments to extended the GW-BSE toolkit.

Excited state properties from the Bethe-Salpeter equation: State-to-state transitions and spin-orbit coupling

The formalism to calculate excited state properties from the GW-Bethe-Salpeter equation (BSE) method is introduced, providing convenient access to excited state absorption, excited state circular dichroism, and excited state optical rotation in the framework of the GW-BSE method. This is achieved using the second-order transition density, which can be obtained by solving a set of auxiliary equations similar to time-dependent density functional theory (TD-DFT). The proposed formulation therefore leads to no increase in the formal computational complexity when compared to the corresponding ground state properties. We further outline the calculation of fully relaxed spin-orbit coupling matrix elements within the GW-BSE method, allowing us to include perturbative corrections for spin-orbit coupling in aforementioned properties. These corrections are also extended to TD-DFT. Excited state absorption and perturbative spin-orbit coupling corrections within GW-BSE are evaluated for a selected set of molecular systems, yielding promising results.

Beyond Electrons: Correlation and Self-Energy in Multicomponent Density Functional Theory

Post-Kohn-Sham methods are used to evaluate the ground-state correlation energy and the orbital self-energy of systems consisting of multiple flavors of different fermions. Starting from multicomponent density functional theory, suitable ways to arrive at the corresponding multicomponent random-phase approximation and the multicomponent Green's functionG W ${GW}$ approximation, including relativistic effects, are outlined. Given the importance of both of this methods in the development of modern Kohn-Sham density functional approximations, this work will provide a foundation to design advanced multicomponent density functional approximations. Additionally, theG W ${GW}$ quasiparticle energies are needed to study light-matter interactions with the Bethe-Salpeter equation.

NWChem: Recent and Ongoing Developments

This paper summarizes developments in the NWChem computational chemistry suite since the last major release (NWChem 7.0.0). Specifically, we focus on functionality, along with input blocks, that is accessible in the current stable release (NWChem 7.2.0) and in the "master" development branch, interfaces to quantum computing simulators, interfaces to external libraries, the NWChem github repository, and containerization of NWChem executable images. Some ongoing developments that will be available in the near future are also discussed.

TDDFT versus GW/BSE Methods for Prediction of Light Absorption and Emission in a TADF Emitter

Design concepts for organic light emitting diode (OLED) emitters, which exhibit thermally activated delayed fluorescence (TADF) and thereby achieve quantum yields exceeding 25%, depend on singlet-triplet splitting energies of order kT to allow reverse intersystem crossing at ambient temperatures. Simulation methods for these systems must be able to treat relatively large organic molecules, as well as predict their excited state energies, transition energies, singlet-triplet splittings, and absorption and emission cross sections with reasonable accuracy, in order to prove useful in the design process. Here we compare predictions of TDDFT with M06-2X and ωB97X-D exchange-correlation functionals and a G_oW_o@HF/BSE method for these quantities in the well-studied DPTZ-DBTO2 TADF emitter molecule. Geometry optimization is performed for ground state (GS) and lowest donor-acceptor charge transfer (CT) state for each functional. Optical absorption and emission cross sections and energies are calculated at these geometries. Relaxation energies are on the order of 0.5 eV, and the importance of obtaining excited state equilibrium geometries in predicting delayed fluorescence is demonstrated. There are clear trends in predictions of G_oW_o@HF/BSE, and TDDFT/ωB97X-D and M06-2X methods in which the former method favors local exciton (LE) states while the latter favors DA CT states and ωB97X-D makes intermediate predictions. G_oW_o@HF/BSE suffers from triplet instability for LE states but not CT states relevant for TADF. Shifts in HOMO and LUMO levels on adding a conductor-like polarizable continuum model dielectric background are used to estimate changes in excitation energies on going from the gas phase to a solvated molecule.

Combining Renormalized Singles GW Methods with the Bethe-Salpeter Equation for Accurate Neutral Excitation Energies

We apply the renormalized singles (RS) Green's function in the Bethe-Salpeter equation (BSE)/GW approach to predict accurate neutral excitation energies of molecular systems. The BSE calculations are performed on top of the G_RSW_RS method, which uses the RS Green's function also for the computation of the screened Coulomb interaction W. We show that the BSE/G_RSW_RS approach significantly outperforms BSE/G₀W₀ for predicting excitation energies of valence, Rydberg, and charge-transfer (CT) excitations by benchmarking the Truhlar-Gagliardi set, Stein CT set, and an atomic Rydberg test set. For the Truhlar-Gagliardi test set, BSE/G_RSW_RS provides comparable accuracy to time-dependent density functional theory (TDDFT) and is slightly better than BSE starting from eigenvalue self-consistent GW (evGW). For the Stein CT test set, BSE/G_RSW_RS significantly outperforms BSE/G₀W₀ and TDDFT with the accuracy comparable to BSE/evGW. We also show that BSE/G_RSW_RS predicts Rydberg excitation energies of atomic systems well. Besides the excellent accuracy, BSE/G_RSW_RS largely eliminates the dependence on the choice of the density functional approximation. This work demonstrates that the BSE/G_RSW_RS approach is accurate and efficient for predicting excitation energies for a broad range of systems, which expands the applicability of the BSE/GW approach.

Quasiparticle Self-Consistent GW-Bethe-Salpeter Equation Calculations for Large Chromophoric Systems

The GW-Bethe–Salpeter equation (BSE) method is promising for calculating the low-lying excitonic states of molecular systems. However, so far it has only been applied to rather small molecules and in the commonly implemented diagonal approximations to the electronic self-energy, it depends on a mean-field starting point. We describe here an implementation of the self-consistent and starting-point-independent quasiparticle self-consistent (qsGW)-BSE approach, which is suitable for calculations on large molecules. We herein show that eigenvalue-only self-consistency can lead to an unfaithful description of some excitonic states for chlorophyll dimers while the qsGW-BSE vertical excitation energies (VEEs) are in excellent agreement with spectroscopic experiments for chlorophyll monomers and dimers measured in the gas phase. Furthermore, VEEs from time-dependent density functional theory calculations tend to disagree with experimental values and using different range-separated hybrid (RSH) kernels does change the VEEs by up to 0.5 eV. We use the new qsGW-BSE implementation to calculate the lowest excitation energies of the six chromophores of the photosystem II (PSII) reaction center (RC) with nearly 2000 correlated electrons. Using more than 11,000 (6000) basis functions, the calculation could be completed in less than 5 (2) days on a single modern compute node. In agreement with previous TD-DFT calculations using RSH kernels on models that also do not include environmental effects, our qsGW-BSE calculations only yield states with local characters in the low-energy spectrum of the hexameric complex. Earlier works with RSH kernels have demonstrated that the protein environment facilitates the experimentally observed interchromophoric charge transfer. Therefore, future research will need to combine correlation effects beyond TD-DFT with an explicit treatment of environmental electrostatics.

Many-body effects in the X-ray absorption spectra of liquid water

SignificanceIn X-ray absorption spectroscopy, an electron-hole excitation probes the local atomic environment. The interpretation of the spectra requires challenging theoretical calculations, particularly in a system like liquid water, where quantum many-body effects and molecular disorder play an important role. Recent advances in theory and simulation make possible new calculations that are in good agreement with experiment, without recourse to commonly adopted approximations. Based on these calculations, the three features observed in the experimental spectra are unambiguously attributed to excitonic effects with different characteristic correlation lengths, which are distinctively affected by perturbations of the underlying H-bond structure induced by temperature changes and/or by isotopic substitution. The emerging picture of the water structure is fully consistent with the conventional tetrahedral model.

Combining localized orbital scaling correction and Bethe-Salpeter equation for accurate excitation energies

We applied localized orbital scaling correction (LOSC) in Bethe-Salpeter equation (BSE) to predict accurate excitation energies for molecules. LOSC systematically eliminates the delocalization error in the density functional approximation and is capable of approximating quasiparticle (QP) energies with accuracy similar to or better than GW Green's function approach and with much less computational cost. The QP energies from LOSC, instead of commonly used G0W0 and evGW, are directly used in BSE. We show that the BSE/LOSC approach greatly outperforms the commonly used BSE/G0W0 approach for predicting excitations with different characters. For the calculations of Truhlar-Gagliardi test set containing valence, charge transfer, and Rydberg excitations, BSE/LOSC with the Tamm-Dancoff approximation provides a comparable accuracy to time-dependent density functional theory (TDDFT) and BSE/evGW. For the calculations of Stein CT test set and Rydberg excitations of atoms, BSE/LOSC considerably outperforms both BSE/G0W0 and TDDFT approaches with a reduced starting point dependence. BSE/LOSC is, thus, a promising and efficient approach to calculate excitation energies for molecular systems.

Assessing the Performances of CASPT2 and NEVPT2 for Vertical Excitation Energies

Methods able to simultaneously account for both static and dynamic electron correlations have often been employed, not only to model photochemical events but also to provide reference values for vertical transition energies, hence allowing benchmarking of lower-order models. In this category, both the complete-active-space second-order perturbation theory (CASPT2) and the N-electron valence state second-order perturbation theory (NEVPT2) are certainly popular, the latter presenting the advantage of not requiring the application of the empirical ionization-potential-electron-affinity (IPEA) and level shifts. However, the actual accuracy of these multiconfigurational approaches is not settled yet. In this context, to assess the performances of these approaches, the present work relies on highly accurate (±0.03 eV) aug-cc-pVTZ vertical transition energies for 284 excited states of diverse character (174 singlet, 110 triplet, 206 valence, 78 Rydberg, 78 n → π*, 119 π → π*, and 9 double excitations) determined in 35 small- to medium-sized organic molecules containing from three to six non-hydrogen atoms. The CASPT2 calculations are performed with and without IPEA shift and compared to the partially contracted (PC) and strongly contracted (SC) variants of NEVPT2. We find that both CASPT2 with IPEA shift and PC-NEVPT2 provide fairly reliable vertical transition energy estimates, with slight overestimations and mean absolute errors of 0.11 and 0.13 eV, respectively. These values are found to be rather uniform for the various subgroups of transitions. The present work completes our previous benchmarks focused on single-reference wave function methods ( J. Chem. Theory Comput. 2018, 14, 4360; J. Chem. Theory Comput. 2020, 16, 1711), hence allowing for a fair comparison between various families of electronic structure methods. In particular, we show that ADC(2), CCSD, and CASPT2 deliver similar accuracies for excited states with a dominant single-excitation character.

NMR Coupling Constants Based on the Bethe-Salpeter Equation in the GW Approximation

We present the first steps to extend the Green's function GW method and the Bethe-Salpeter equation (BSE) to molecular response properties such as nuclear magnetic resonance (NMR) indirect spin-spin coupling constants. We discuss both a nonrelativistic one-component and a quasi-relativistic two-component formalism. The latter describes scalar-relativistic and spin-orbit effects and allows us to study heavy-element systems with reasonable accuracy. Efficiency is maintained by the application of the resolution of the identity approximation throughout. The performance is demonstrated using conventional central processing units (CPUs) and modern graphics processing units (GPUs) for molecules involving several thousand basis functions. Our results show that a large amount of Hartree-Fock exchange is vital to provide a sufficient Kohn-Sham starting point to compute the GW quasi-particle energies. As the GW-BSE approach is generally less accurate for triplet excitations or related properties such as the Fermi-contact interaction, the admixture of the Kohn-Sham correlation kernel through the contracted BSE (cBSE) method improves the results for NMR coupling constants. This leads to remarkable results when combined with the eigenvalue-only self-consistent variant (evGW) and Becke's half and half functional (BH&HLYP) or the CAM-QTP family. The developed methodology is used to calculate the Karplus curve of tin molecules, illustrating its applicability to extended chemically relevant molecules. Here, the GW-cBSE method improves upon the chosen BH&HLYP Kohn-Sham starting points.

Low-Order Scaling G0W0 by Pair Atomic Density Fitting

We derive a low-scaling G0W0 algorithm for molecules using pair atomic density fitting (PADF) and an imaginary time representation of the Green's function and describe its implementation in the Slater type orbital (STO)-based Amsterdam density functional (ADF) electronic structure code. We demonstrate the scalability of our algorithm on a series of water clusters with up to 432 atoms and 7776 basis functions and observe asymptotic quadratic scaling with realistic threshold qualities controlling distance effects and basis sets of triple-ζ (TZ) plus double polarization quality. Also owing to a very small prefactor, a G0W0 calculation for the largest of these clusters takes only 240 CPU hours with these settings. We assess the accuracy of our algorithm for HOMO and LUMO energies in the GW100 database. With errors of 0.24 eV for HOMO energies on the quadruple-ζ level, our implementation is less accurate than canonical all-electron implementations using the larger def2-QZVP GTO-type basis set. Apart from basis set errors, this is related to the well-known shortcomings of the GW space-time method using analytical continuation techniques as well as to numerical issues of the PADF approach of accurately representing diffuse atomic orbital (AO) products. We speculate that these difficulties might be overcome by using optimized auxiliary fit sets with more diffuse functions of higher angular momenta. Despite these shortcomings, for subsets of medium and large molecules from the GW5000 database, the error of our approach using basis sets of TZ and augmented double-ζ (DZ) quality is decreasing with system size. On the augmented DZ level, we reproduce canonical, complete basis set limit extrapolated reference values with an accuracy of 80 meV on average for a set of 20 large organic molecules. We anticipate our algorithm, in its current form, to be very useful in the study of single-particle properties of large organic systems such as chromophores and acceptor molecules.

Understanding real-time time-dependent density-functional theory simulations of ultrafast laser-induced dynamics in organic molecules

Real-time time-dependent density functional theory, in conjunction with the Ehrenfest molecular dynamics scheme, is becoming a popular methodology to investigate ultrafast phenomena on the nanoscale. Thanks to recent developments, it is also possible to explicitly include in the simulations a time-dependent laser pulse, thereby accessing the transient excitation regime. However, the complexity entailed in these calculations calls for in-depth analysis of the accessible and yet approximate (either "dressed" or "bare") quantities in order to evaluate their ability to provide us with a realistic picture of the simulated processes. In this work, we analyze the ultrafast dynamics of three small molecules (ethylene, benzene, and thiophene) excited by a resonant laser pulse in the framework of the adiabatic local-density approximation. The electronic response to the laser perturbation in terms of induced dipole moment and excited-state population is compared to the results given by an exactly solvable two-level model. In this way, we can interpret the charge-carrier dynamics in terms of simple estimators, such as the number of excited electrons. From the computed transient absorption spectra, we unravel the appearance of nonlinear effects such as excited-state absorption and vibronic coupling. In this way, we observe that the laser excitation affects the vibrational spectrum by enhancing the anharmonicities therein, while the coherent vibrational motion contributes to stabilizing the electronic excitation already within a few tens of femtoseconds.

Ab initio modelling of local interfaces in doped organic semiconductors

Doping in organic semiconductors remains a debated issue from both an experimental and ab initio perspective. Due to the complexity of these systems, which exhibit a low degree of crystallinity and high level of disorder, modelling doped organic semiconductors from first-principles calculations is not a trivial task, as their electronic and optical properties are sensitive to the choice of initial geometries. A crucial aspect to take into account, in view of rationalizing the electronic structure of these materials through ab initio calculations, is the role of local donor/acceptor interfaces. We address this problem in the framework of state-of-the-art density-functional theory and many-body perturbation theory, investigating the structural, electronic, and optical properties of quaterthiophene and sexithiophene oligomers doped by 2,3,5,6-tetrafluoro-7,7,8,8-tetracyano-quinodimethane (F4TCNQ). We consider different model structures ranging from isolated dimers and trimers, to periodic stacks. Our results demonstrate that the choice of the initial geometry critically impacts the resulting electronic structure and the degree of charge transfer in the materials, depending on the amount and on the nature of the local interfaces between donor and acceptor species. The optical spectra appear less sensitive to these parameters at least from a first glance, although a quantitative analysis of the excitations reveals that their Frenkel or charge-transfer character is affected by the characteristics of the donor/acceptor interfaces as well as by the donor length. Our findings represent an important step forward towards an insightful first-principles description of the microscopic properties of doped organic semiconductors complementary to experiments.

Many-body perturbation theory calculations using the yambo code

yambo is an open source project aimed at studying excited state properties of condensed matter systems from first principles using many-body methods. As input, yambo requires ground state electronic structure data as computed by density functional theory codes such as Quantum ESPRESSO and Abinit. yambo's capabilities include the calculation of linear response quantities (both independent-particle and including electron-hole interactions), quasi-particle corrections based on the GW formalism, optical absorption, and other spectroscopic quantities. Here we describe recent developments ranging from the inclusion of important but oft-neglected physical effects such as electron-phonon interactions to the implementation of a real-time propagation scheme for simulating linear and non-linear optical properties. Improvements to numerical algorithms and the user interface are outlined. Particular emphasis is given to the new and efficient parallel structure that makes it possible to exploit modern high performance computing architectures. Finally, we demonstrate the possibility to automate workflows by interfacing with the yambopy and AiiDA software tools.

Combining Wave Function Methods with Density Functional Theory for Excited States

We review state-of-the-art electronic structure methods based both on wave function theory (WFT) and density functional theory (DFT). Strengths and limitations of both the wave function and density functional based approaches are discussed, and modern attempts to combine these two methods are presented. The challenges in modeling excited-state chemistry using both single-reference and multireference methods are described. Topics covered include background, combining density functional theory with single-configuration wave function theory, generalized Kohn-Sham (KS) theory, global hybrids, range-separated hybrids, local hybrids, using KS orbitals in many-body theory (including calculations of the self-energy and the GW approximation), Bethe-Salpeter equation, algorithms to accelerate GW calculations, combining DFT with multiconfigurational WFT, orbital-dependent correlation functionals based on multiconfigurational WFT, building multiconfigurational wave functions from KS configurations, adding correlation functionals to multiconfiguration self-consistent-field (MCSCF) energies, combining DFT with configuration-interaction singles by means of time-dependent DFT, using range separation to combine DFT with MCSCF, embedding multiconfigurational WFT in DFT, and multiconfiguration pair-density functional theory.

Assessing GW Approaches for Predicting Core Level Binding Energies

Here we present a systematic study on the performance of different GW approaches: G₀W₀, G₀W₀ with linearized quasiparticle equation (lin-G₀W₀), and quasiparticle self-consistent GW (qsGW), in predicting core level binding energies (CLBEs) on a series of representative molecules comparing to Kohn-Sham (KS) orbital energy-based results. KS orbital energies obtained using the PBE functional are 20-30 eV lower in energy than experimental values obtained from X-ray photoemission spectroscopy (XPS), showing that any Koopmans-like interpretation of KS core level orbitals fails dramatically. Results from qsGW lead to CLBEs that are closer to experimental values from XPS, yet too large. For the qsGW method, the mean absolute error is about 2 eV, an order of magnitude better than plain KS PBE orbital energies and quite close to predictions from ΔSCF calculations with the same functional, which are accurate within ∼1 eV. Smaller errors of ∼0.6 eV are found for qsGW CLBE shifts, again similar to those obtained using ΔSCF PBE. The computationally more affordable G₀W₀ approximation leads to results less accurate than qsGW, with an error of ∼9 eV for CLBEs and ∼0.9 eV for their shifts. Interestingly, starting G₀W₀ from PBE0 reduces this error to ∼4 eV with a slight improvement on the shifts as well (∼0.4 eV). The validity of the G₀W₀ results is however questionable since only linearized quasiparticle equation results can be obtained. The present results pave the way to estimate CLBEs in periodic systems where ΔSCF calculations are not straightforward although further improvement is clearly needed.

The Bethe–Salpeter equation in chemistry: relations with TD-DFT, applications and challenges

We review the Bethe–Salpeter formalism and analyze its performances for the calculation of the excited state properties of molecular systems.

Accurate description of the electronic structure of organic semiconductors by GW methods

Electronic properties associated with charged excitations, such as the ionization potential (IP), the electron affinity (EA), and the energy level alignment at interfaces, are critical parameters for the performance of organic electronic devices. To computationally design organic semiconductors and functional interfaces with tailored properties for target applications it is necessary to accurately predict these properties from first principles. Many-body perturbation theory is often used for this purpose within the GW approximation, where G is the one particle Green's function and W is the dynamically screened Coulomb interaction. Here, the formalism of GW methods at different levels of self-consistency is briefly introduced and some recent applications to organic semiconductors and interfaces are reviewed.

Modeling the Photochrome-TiO2 Interface with Bethe-Salpeter and Time-Dependent Density Functional Theory Methods

Hybrid organic-inorganic semiconductor systems have important applications in both molecular electronics and photoresponsive materials. The characterizations of the interface and of the electronic excited-states of these hybrid systems remain a challenge for state-of-the-art computational methods, as the systems of interest are large. In the present investigation, we present for the first time a many-body Green's function Bethe-Salpeter investigation of a series of photochromic molecules adsorbed onto TiO₂ nanoclusters. On the basis of these studies, the performance of time-dependent density functional theory (TD-DFT) calculations is assessed. In addition, the photochromic properties of different hybrid systems are also evaluated. This work shows that qualitatively different conclusions can be reached with TD-DFT relying on various exchange-correlation functionals for such organic-inorganic interfaces and paves the way to more accurate simulation of many hybrid materials.

Combining the GW formalism with the polarizable continuum model: A state-specific non-equilibrium approach

We have implemented the polarizable continuum model within the framework of the many-body Green's function GW formalism for the calculation of electron addition and removal energies in solution. The present formalism includes both ground-state and non-equilibrium polarization effects. In addition, the polarization energies are state-specific, allowing to obtain the bath-induced renormalisation energy of all occupied and virtual energy levels. Our implementation is validated by comparisons with ΔSCF calculations performed at both the density functional theory and coupled-cluster single and double levels for solvated nucleobases. The present study opens the way to GW and Bethe-Salpeter calculations in disordered condensed phases of interest in organic optoelectronics, wet chemistry, and biology.

Benchmark of Bethe-Salpeter for Triplet Excited-States

We have evaluated the accuracy of the Bethe-Salpeter singlet-triplet transition energies as well as singlet-triplet and triplet-triplet splittings for 20 organic molecules, using as reference the CC3 values determined by Thiel and co-workers with both the TZVP and aug-cc-pVTZ atomic basis sets. Our excitation energies are obtained on the basis of GW quasiparticle energy levels that are self-consistently converged with respect to the starting DFT eigenvalues. In its current form, BSE/GW is often unable to provide a balanced description of both singlet and triplet excited-states. While the singlet-singlet and triplet-triplet energy separations are obtained accurately, triplets are located too close in energy from the ground-state, by typically -0.55 eV when using standard functionals to generate the starting eigenstates. Applying the Tamm-Dancoff approximation upshifts the BSE triplet energies and allows reducing this error to ca. -0.40 eV, while using M06-HF eigenstates allows a further increase and hence a reduction of the error for triplet states, but at the cost of larger errors for the singlet excited-states. At this stage, the most accurate TD-DFT estimates therefore remain competitive for computing singlet-triplet transition energies. Indeed, with M06-2X, irrespective of the application or not of the Tamm-Dancoff approximation and of the selected atomic basis set, the deviations obtained with TD-DFT are rather small.

Theoretical description of protein field effects on electronic excitations of biological chromophores

Photoinitiated phenomena play a crucial role in many living organisms. Plants, algae, and bacteria absorb sunlight to perform photosynthesis, and convert water and carbon dioxide into molecular oxygen and carbohydrates, thus forming the basis for life on Earth. The vision of vertebrates is accomplished in the eye by a protein called rhodopsin, which upon photon absorption performs an ultrafast isomerisation of the retinal chromophore, triggering the signal cascade. Many other biological functions start with the photoexcitation of a protein-embedded pigment, followed by complex processes comprising, for example, electron or excitation energy transfer in photosynthetic complexes. The optical properties of chromophores in living systems are strongly dependent on the interaction with the surrounding environment (nearby protein residues, membrane, water), and the complexity of such interplay is, in most cases, at the origin of the functional diversity of the photoactive proteins. The specific interactions with the environment often lead to a significant shift of the chromophore excitation energies, compared with their absorption in solution or gas phase. The investigation of the optical response of chromophores is generally not straightforward, from both experimental and theoretical standpoints; this is due to the difficulty in understanding diverse behaviours and effects, occurring at different scales, with a single technique. In particular, the role played by ab initio calculations in assisting and guiding experiments, as well as in understanding the physics of photoactive proteins, is fundamental. At the same time, owing to the large size of the systems, more approximate strategies which take into account the environmental effects on the absorption spectra are also of paramount importance. Here we review the recent advances in the first-principle description of electronic and optical properties of biological chromophores embedded in a protein environment. We show their applications on paradigmatic systems, such as the light-harvesting complexes, rhodopsin and green fluorescent protein, emphasising the theoretical frameworks which are of common use in solid state physics, and emerging as promising tools for biomolecular systems.

A systematic benchmark of the ab initio Bethe-Salpeter equation approach for low-lying optical excitations of small organic molecules

The predictive power of the ab initio Bethe-Salpeter equation (BSE) approach, rigorously based on many-body Green's function theory but incorporating information from density functional theory, has already been demonstrated for the optical gaps and spectra of solid-state systems. Interest in photoactive hybrid organic/inorganic systems has recently increased and so has the use of the BSE for computing neutral excitations of organic molecules. However, no systematic benchmarks of the BSE for neutral electronic excitations of organic molecules exist. Here, we study the performance of the BSE for the 28 small molecules in Thiel's widely used time-dependent density functional theory benchmark set [Schreiber et al., J. Chem. Phys. 128, 134110 (2008)]. We observe that the BSE produces results that depend critically on the mean-field starting point employed in the perturbative approach. We find that this starting point dependence is mainly introduced through the quasiparticle energies obtained at the intermediate GW step and that with a judicious choice of starting mean-field, singlet excitation energies obtained from BSE are in excellent quantitative agreement with higher-level wavefunction methods. The quality of the triplet excitations is slightly less satisfactory.

Benchmarking the Bethe-Salpeter Formalism on a Standard Organic Molecular Set

We perform benchmark calculations of the Bethe-Salpeter vertical excitation energies for the set of 28 molecules constituting the well-known Thiel’s set, complemented by a series of small molecules representative of the dye chemistry field. We show that Bethe-Salpeter calculations based on a molecular orbital energy spectrum obtained with non-self-consistent G0W0 calculations starting from semilocal DFT functionals dramatically underestimate the transition energies. Starting from the popular PBE0 hybrid functional significantly improves the results even though this leads to an average -0.59 eV redshift compared to reference calculations for Thiel’s set. It is shown, however, that a simple self-consistent scheme at the GW level, with an update of the quasiparticle energies, not only leads to a much better agreement with reference values, but also significantly reduces the impact of the starting DFT functional. On average, the Bethe-Salpeter scheme based on self-consistent GW calculations comes close to the best time-dependent DFT calculations with the PBE0 functional with a 0.98 correlation coefficient and a 0.18 (0.25) eV mean absolute deviation compared to TD-PBE0 (theoretical best estimates) with a tendency to be red-shifted. We also observe that TD-DFT and the standard adiabatic Bethe-Salpeter implementation may differ significantly for states implying a large multiple excitation character.