Theoretical S1 <- S0 absorption energies of the anionic forms of oxyluciferin by Variational Monte Carlo and Many Body Green's Function Theory

TL;DR

This study combines Variational Monte Carlo and Many Body Green's Function Theory to accurately compute absorption energies of oxyluciferin's anionic forms, revealing the importance of geometry and solvent effects on excitation energies.

Contribution

It introduces a novel approach combining VMC and MBGFT for precise absorption energy calculations of bioluminescent molecules, outperforming standard methods.

Findings

MBGFT results agree with experimental data for keto-1 form

Ground state geometry quality significantly affects excitation energies

Solvent embedding is crucial for accurate excited state calculations

Abstract

The structures of three negatively charged forms (anionic keto-1 and enol-1, dianonic enol-2) of oxyluciferin (OxyLuc), which are the most probable emitters responsible for the firefly bioluminescence, have been fully relaxed at the variational Monte Carlo (VMC) level. Absorption energies of the S1<-S0 vertical transition have been computed using different levels of theory, such as TDDFT, CC2 and many body Green's function Theory (MBGFT). The use of MBGFT, by means of the Bethe-Salpeter (BS) formalism, on VMC structures provides results in excellent agreement with the value (2.26(8) eV) obtained by action spectroscopy experiments for the keto-1 form (2.32 eV). To unravel the role of the quality of the optimized ground state geometry, BS excitation energies have also been computed on CASSCF geometries, inducing a non negligible blue shift (0.08 and 0.07 eV for keto-1 and enol-1 forms, respectively) with respect to the VMC ones. Structural effects have been analyzed in terms of over- or under-correlation along the conjugated bonds of OxyLuc by using different methods for the ground-state optimization. The relative stability of the S$_1$ state for the keto-1 and enol-1 forms depends on the method chosen for the excited state calculation, thus representing a fundamental caveat for any theoretical study on these systems. Finally, Kohn-Sham HOMO and LUMO orbitals of enol-2 are (nearly) bound only when the dianion is embedded into a solvent (water and toluene in the present work); excited state calculations are therefore meaningful only in presence of a dielectric medium which localizes the electronic density. The combination of VMC for the ground state geometry and BS formalism for the absorption spectra clearly outperforms standard TDDFT and quantum chemistry approaches.