Refine
Has Fulltext
- no (4)
Document Type
- Article (4)
Language
- English (4)
Is part of the Bibliography
- yes (4)
Institute
We propose an entirely nonempirical and computationally efficient scheme to calculate highly reliable vibrationally resolved photoelectron spectra for molecules from first principles. To this end, we combine nonempirically tuned long-range corrected hybrid functionals with non-self-consistent many-body perturbation theory in the G(0)W(0) approximation and a Franck-Condon multimode analysis based on DFT-calculated frequencies. The vibrational analysis allows for a direct comparison of the GW-calculated spectra to gas-phase ultraviolet photoelectron measurements of neutral and anionic molecules, respectively. Direct comparison of the calculated peak maxima with experiment yields mean absolute errors below 0.1 eV for ionization potentials, electron affinities, and fundamental gaps, clearly outperforming commonly used G(0)W(0) approaches at similar numerical costs.
The performance of non-empirically tuned long-range corrected hybrid functionals for the prediction of vertical ionization potentials (IPs) and electron affinities (EAs) is assessed for a set of 24 organic acceptor molecules. Basis set extrapolated coupled cluster singles, doubles, and perturbative triples [CCSD(T)] calculations serve as a reference for this study. Compared to standard exchange-correlation functionals, tuned long-range corrected hybrid functionals produce highly reliable results for vertical IPs and EAs, yielding mean absolute errors on par with computationally more demanding GW calculations. In particular, it is demonstrated that long-range corrected hybrid functionals serve as ideal starting points for non-self-consistent GW calculations.
The performance of different GW methods is assessed for a set of 24 organic acceptors. Errors are evaluated with respect to coupled cluster singles, doubles, and perturbative triples [CCSD(T)] reference data for the vertical ionization potentials (IPs) and electron affinities (EAs), extrapolated to the complete basis set limit. Additional comparisons are made to experimental data, where available. We consider fully self-consistent GW (scGW), partial self-consistency in the Green’s function (scGW0), non-self-consistent G0W0 based on several mean-field starting points, and a “beyond GW” second-order screened exchange (SOSEX) correction to G0W0. We also describe the implementation of the self-consistent Coulomb hole with screened exchange method (COHSEX), which serves as one of the mean-field starting points. The best performers overall are G0W0+SOSEX and G0W0 based on an IP-tuned long-range corrected hybrid functional with the former being more accurate for EAs and the latter for IPs. Both provide a balanced treatment of localized vs delocalized states and valence spectra in good agreement with photoemission spectroscopy (PES) experiments.
Vibrationally resolved lowest-energy bands of the photoelectron spectra (PES) of adamantane, diamantane, and urotropine were simulated by a time-dependent correlation function approach within the harmonic approximation. Geometries and normal modes for neutral and cationic molecules were obtained from B3LYP hybrid density functional theory (DFT). It is shown that the simulated spectra reproduce the experimentally observed vibrational finestructure (or its absence) quite well. Origins of the finestructure are discussed and related to recurrences of autocorrelation functions and dominant vibrations. Remaining quantitative and qualitative errors of the DFT-derived PES spectra refer to (i) an overall redshift by ∼0.5 eV and (ii) the absence of satellites in the high-energy region of the spectra. The former error is shown to be due to the neglect of many-body corrections to ordinary Kohn-Sham methods, while the latter has been argued to be due to electron-nuclear couplings beyond the Born-Oppenheimer approximation [Gali et al., Nat. Commun. 7, 11327 (2016)].