Reference CC3 Excitation Energies for Organic Chromophores: Benchmarking TD-DFT, BSE/, and Wave Function Methods.

To expand the QUEST database of highly accurate vertical transition energies, we consider a series of large organic chromogens ubiquitous in dye chemistry, such as anthraquinone, azobenzene, BODIPY, and naphthalimide. We compute, at the CC3 level of theory, the singlet and triplet vertical transition energies associated with the low-lying excited states. This leads to a collection of more than 120 new highly accurate excitation energies.

Polarizable Continuum Models and Green's Function Formalism: On the Dynamics of the Solvent Electrons.

The many-body formalism, for the calculation of ionization potentials or electronic affinities, relies on the frequency-dependent dielectric function built from the electronic degrees of freedom. Considering the case of water as a solvent treated within the polarizable continuum model, we explore the impact of restricting the full frequency-dependence of the solvent electronic dielectric response to a frequency-independent (ϵ) optical dielectric constant. For solutes presenting small to large highest-occupied to lowest-unoccupied molecular orbital energy gaps, we show that such a restriction induces errors no larger than a few percent on the energy level shifts from the gas to the solvated phase.

Static versus dynamically polarizable environments within the many-body GW formalism.

Continuum- or discrete-polarizable models for the study of optoelectronic processes in embedded subsystems rely mostly on the restriction of the surrounding electronic dielectric response to its low frequency limit. Such a description hinges on the assumption that the electrons in the surrounding medium react instantaneously to any excitation in the central subsystem, thus treating the environment in the adiabatic limit. Exploiting a recently developed embedded GW formalism with an environment described at the fully ab initio level, we assess the merits of the adiabatic limit with respect to an environment where the full dynamics of the dielectric response are considered.

Assessing the accuracy of TD-DFT excited-state geometries through optimal tuning with GW energy levels.

We study the accuracy of excited state (ES) geometries using optimally tuned LC-PBE functionals with tuning based on GW quasiparticle energies. We compare the results obtained with the PBE, PBE0, non-tuned, and tuned LC-PBE functionals with available high-level CC reference values as well as experimental data. First, we compare ES geometrical parameters obtained for three different types of systems: molecules composed of a few atoms, 4-(dimethylamino)benzonitrile (DMABN), and conjugated dyes.

Disentangling the multiorbital contributions of excitons by photoemission exciton tomography.

Excitons are realizations of a correlated many-particle wave function, specifically consisting of electrons and holes in an entangled state. Excitons occur widely in semiconductors and are dominant excitations in semiconducting organic and low-dimensional quantum materials. To efficiently harness the strong optical response and high tuneability of excitons in optoelectronics and in energy-transformation processes, access to the full wavefunction of the entangled state is critical, but has so far not been feasible.

Excess and excited-state dipole moments of real-life dyes: a comparison between wave-function, BSE/, and TD-DFT values.

In this work, we assess the accuracy of the Bethe-Salpeter equation (BSE) many-body Green's function formalism, adopting the eigenvalue-self-consistent ev exchange-correlation kernel, for the calculation of the excited-state () and excess dipole moments (Δ), the latter ones being the changes of dipole amplitude between the ground and excited states (ES), in organic dyes. We compare the results obtained with wave-function methods [ADC(2), CC2, and CCSD], time-dependent density functional theory (TD-DFT), and BSE/ev levels of theory. First, we compute the evolution of the dipole moments of the two lowest singlet excited states of 4-(dimethylamino)benzonitrile (DMABN) upon twisting of the amino group.

Many-body GW calculations with very large scale polarizable environments made affordable: A fully ab initio QM/QM approach.

We present a many-body GW formalism for quantum subsystems embedded in discrete polarizable environments containing up to several hundred thousand atoms described at a fully ab initio random phase approximation level. Our approach is based on a fragment approximation in the construction of the Green's function and independent-electron susceptibilities. Further, the environing fragments susceptibility matrices are reduced to a minimal but accurate representation preserving low order polarizability tensors through a constrained minimization scheme.

Lagrangian Z-vector approach to Bethe-Salpeter analytic gradients: Assessing approximations.

We present an implementation of excited-state analytic gradients within the Bethe-Salpeter equation formalism using an adapted Lagrangian Z-vector approach with a cost independent of the number of perturbations. We focus on excited-state electronic dipole moments associated with the derivatives of the excited-state energy with respect to an electric field. In this framework, we assess the accuracy of neglecting the screened Coulomb potential derivatives, a common approximation in the Bethe-Salpeter community, as well as the impact of replacing the GW quasiparticle energy gradients by their Kohn-Sham analogs.

Exploring Bethe-Salpeter Excited-State Dipoles: The Challenging Case of Increasingly Long Push-Pull Oligomers.

The change of molecular dipole moment induced by photon absorption is key to interpret the measured optical spectra. Except for compact molecules, time-dependent density functional theory (TD-DFT) remains the only theory allowing to quickly predict excited-state dipoles (μ), albeit with a strong dependency on the selected exchange-correlation functional. This Letter presents the first assessment of the performances of the many-body Green's function Bethe-Salpeter equation (BSE) formalism for the evaluation of the μ.

Excited state potential energy surfaces of -phenylpyrrole upon twisting: reference values and comparison between BSE/ and TD-DFT.

The puzzling case of the mixing between the charge transfer (CT) and local excited (LE) characters upon twisting of the geometry of -phenylpyrrole (-PP) is investigated considering the six low-lying singlet excited states (ES). The theoretical calculations of the potential energy surfaces (PES) have been performed for these states using a Coupled Cluster method accounting for the impact of the contributions from the triples, many-body Green's function and Bethe-Salpeter equation (BSE) formalisms, as well as Time-Dependent Density Functional Theory (TD-DFT) using various exchange-correlation functionals. Our findings confirm that the BSE formalism is more reliable than TD-DFT for close-lying ES with mixed CT/LE nature.

Modeling of excited state potential energy surfaces with the Bethe-Salpeter equation formalism: The 4-(dimethylamino)benzonitrile twist.

We present a benchmark study of excited state potential energy surfaces (PES) using the many-body Green's function GW and Bethe-Salpeter equation (BSE) formalisms, coupled cluster methods, as well as Time-Dependent Density Functional Theory (TD-DFT). More specifically, we investigate the evolution of the two lowest excited states of 4-(dimethylamino)benzonitrile (DMABN) upon the twisting of the amino group, a paradigmatic system for dual fluorescence and excited-state benchmarks. Our results demonstrate that the BSE/GW approach is able to reproduce the correct topology of excited state PES upon geometry changes in both gas and condensed phases.

Reference Energies for Intramolecular Charge-Transfer Excitations.

With the aim of completing our previous efforts devoted to local and Rydberg transitions in organic compounds, we provide a series of highly accurate vertical transition energies for intramolecular charge-transfer transitions occurring in (π-conjugated) molecular compounds. To this end, we apply a composite protocol consisting of linear-response CCSDT excitation energies determined with Dunning's double-ζ basis set corrected by CC3/CCSDT-3 energies obtained with the corresponding triple-ζ basis. Further basis set corrections (up to aug-cc-pVQZ) are obtained at the CCSD and CC2 levels.

Strongly Bound Excitons in Metal-Organic Framework MOF-5: A Many-Body Perturbation Theory Study.

During the past years, one of the most iconic metal-organic frameworks (MOFs), MOF-5, has been characterized as a semiconductor by theory and experiments. Here we employ the many-body perturbation theory in conjunction with the Bethe-Salpeter equation to compute the electronic structure and optical properties of this MOF. The calculations show that MOF-5 is a wide-band-gap insulator with a fundamental gap of ∼8 eV.

Cubic-Scaling All-Electron Calculations with a Separable Density-Fitting Space-Time Approach.

We present an implementation of the space-time approach that allows cubic-scaling all-electron calculations with standard Gaussian basis sets without exploiting any localization or sparsity considerations. The independent-electron susceptibility is constructed in a time representation over a nonuniform distribution of real-space locations {} optimized within a separable resolution-of-the-identity framework to reproduce standard Coulomb-fitting calculations with meV accuracy. The compactness of the obtained {} distribution leads to a crossover with the standard Coulomb-fitting scheme for system sizes below a few hundred electrons.

Dynamical correction to the Bethe-Salpeter equation beyond the plasmon-pole approximation.

The Bethe-Salpeter equation (BSE) formalism is a computationally affordable method for the calculation of accurate optical excitation energies in molecular systems. Similar to the ubiquitous adiabatic approximation of time-dependent density-functional theory, the static approximation, which substitutes a dynamical (i.e.

Orientation dependent molecular electrostatics drives efficient charge generation in homojunction organic solar cells.

Organic solar cells usually utilise a heterojunction between electron-donating (D) and electron-accepting (A) materials to split excitons into charges. However, the use of D-A blends intrinsically limits the photovoltage and introduces morphological instability. Here, we demonstrate that polycrystalline films of chemically identical molecules offer a promising alternative and show that photoexcitation of α-sexithiophene (α-6T) films results in efficient charge generation.

The Bethe-Salpeter Equation Formalism: From Physics to Chemistry.

The Bethe-Salpeter equation (BSE) formalism is steadily asserting itself as a new efficient and accurate tool in the ensemble of computational methods available to chemists in order to predict optical excitations in molecular systems. In particular, the combination of the so-called approximation, giving access to reliable ionization energies and electron affinities, and the BSE formalism, able to model UV/vis spectra, has shown to provide accurate singlet excitation energies with a typical error of 0.1-0.

Pros and Cons of the Bethe-Salpeter Formalism for Ground-State Energies.

The combination of the many-body Green's function approximation and the Bethe-Salpeter equation (BSE) formalism has shown to be a promising alternative to time-dependent density functional theory (TD-DFT) for computing vertical transition energies and oscillator strengths in molecular systems. The BSE formalism can also be employed to compute ground-state correlation energies thanks to the adiabatic-connection fluctuation-dissipation theorem (ACFDT). Here, we study the topology of the ground-state potential energy surfaces (PESs) of several diatomic molecules near their equilibrium bond length.

Accurate Prediction of the S Excitation Energy in Solvated Azobenzene Derivatives via Embedded Orbital-Tuned Bethe-Salpeter Calculations.

By employing the Bethe-Salpeter formalism coupled with a nonequilibrium embedding scheme, we demonstrate that the paradigmatic case of S band separation between and in azobenzene derivatives can be computed with excellent accuracy compared to experimental optical spectra. Besides embedding, we show that the choice of the Kohn-Sham exchange correlation functional for DFT is critical, despite the iterative convergence of quasiparticle energies. We address this by adopting an orbital-tuning approach via the global hybrid functional, PBEh, yielding an environment-consistent ionization potential.

Robust Analytic-Continuation Approach to Many-Body Calculations.

The analytic continuation of the self-energy from the imaginary to the real-energy axis is a central difficulty for approaches exploiting the favorable properties of response functions at imaginary frequencies. Within a scheme merging contour-deformation and analytic-continuation techniques, we show on the basis of extensive calculations for large molecular sets that it is preferable to perform an analytic continuation of the dynamically screened-Coulomb potential rather than the much more structured self-energy operator. The case of states lying far away from the gap, including core states, is addressed by generalizing the analytic-continuation scheme, accounting further for quasiparticle lifetimes.

Separable resolution-of-the-identity with all-electron Gaussian bases: Application to cubic-scaling RPA.

We explore a separable resolution-of-the-identity (RI) formalism built on quadratures over limited sets of real-space points designed for all-electron calculations. Our implementation preserves, in particular, the use of common atomic orbitals and their related auxiliary basis sets. The setup of the present density fitting scheme, i.