Multiconfigurational Quantum Chemistry Determinations of Absorption Cross Sections (σ) in the Gas Phase and Molar Extinction Coefficients (ε) in Aqueous Solution and Air-Water Interface.

Theoretical determinations of absorption cross sections (σ) in the gas phase and molar extinction coefficients (ε) in condensed phases (water solution, interfaces or surfaces, protein or nucleic acids embeddings, etc.) are of interest when rates of photochemical processes, J = ∫ ϕ(λ) σ(λ) I(λ) dλ, are needed, where ϕ(λ) and I(λ) are the quantum yield of the process and the irradiance of the light source, respectively, as functions of the wavelength λ. Efficient computational strategies based on single-reference quantum-chemistry methods have been developed enabling determinations of line shapes or, in some cases, achieving rovibrational resolution. Developments are however lacking for strongly correlated problems, with many excited states, high-order excitations, and/or near degeneracies between states of the same and different spin multiplicities. In this work, we define and compare the performance of distinct computational strategies using multiconfigurational quantum chemistry, nuclear sampling of the chromophore (by means of molecular dynamics, ab initio molecular dynamics, or Wigner sampling), and conformational and statistical sampling of the environment (by means of molecular dynamics). A new mathematical approach revisiting previous absolute orientation algorithms is also developed to improve alignments of geometries. These approaches are benchmarked through the nπ* band of acrolein not only in the gas phase and water solution but also in a gas-phase/water interface, a common situation for instance in atmospheric chemistry. Subsequently, the best strategy is used to compute the absorption band for the adduct formed upon addition of an OH radical to the C6 position of uracil and compared with the available experimental data. Overall, quantum Wigner sampling of the chromophore with molecular dynamics sampling of the environment with CASPT2 electronic-structure determinations arise as a powerful methodology to predict meaningful σ(λ) and ε(λ) band line shapes with accurate absolute intensities.

Light emission colour modulation study of oxyluciferin synthetic analogues via QM and QM/MM approaches.

Firefly oxyluciferin is the chemical product of bioluminescence responsible for light emission. Experiments have already shown that different analogues of natural oxyluciferin, exhibit different emission colours. In particular, the structure of natural oxyluciferin has been modified by atom or group substitutions. However, a rationalization of the origin of the bioluminescence emission colour modulation of these analogues has still not been reported. For these reasons, the aim of this study is to explain the influence of structural modifications within the natural oxyluciferin on the colour modulation of bioluminescence. To do this, natural firefly oxyluciferin and three synthetic analogues whose experimental bioluminescence spectra are red- and blue-shifted compared to the natural one were studied. The absorption and emission transition energies have been calculated at the Time-dependent density functional theory (TD-DFT) level using both quantum mechanics (QM) and quantum mechanics/molecular mechanics (QM/MM) methods. Moreover, the solvent (water using the PCM model) and the protein surrounding effect have also been considered. The predicted emission spectra are in quite good agreement with the available experimental spectra, validating the methodology followed in this study. In particular, it was demonstrated that using the QM/MM approach, and considering explicitly the protein environment, the experimental bioluminescence spectra can be reproduced. Furthermore, this study shows that the substitution within the oxyluciferin structure causes a change of its electronic distribution and energies of the HOMO and LUMO orbitals involved in the vertical transitions, leading to different light emission colours. This work will promote future studies focused on luciferin mutations guided by the prediction of their bioluminescence emission spectra.

The role of solvation models on the computed absorption and emission spectra: the case of fireflies oxyluciferin.

Surrounding effects are crucial to successfully simulate the absorption and emission spectra of molecular systems. In this work we test different solvation models to compute transition energies and to simulate the spectra of oxyluciferin responsible for the light emission in fireflies and its derivatives. We demonstrate that, within the PCM model, the IBSF formalism is suitable for computing the transition energies of the oxyluciferin chemical forms characterized by a charge transfer character. On the other hand, the LR approach could be used for the chemical forms where an almost negligible charge transfer takes place. Moreover, we demonstrate that explicit solvation models, applied by QM/MM calculations, are needed to accurately reproduce the experimental shape of the spectra. Finally, the vibrationally resolved spectra using a solvation model (implicit or microsolvation) is computed. Some noticeable differences arise when considering the implicit solvation with respect to gas phase vibrational spectra, while small changes were found when explicit water molecules within a microsolvated model are considered.

Simulation and Analysis of the Spectroscopic Properties of Oxyluciferin and Its Analogues in Water.

Firefly bioluminescence is a quite efficient process largely used for numerous applications. However, some fundamental photochemical properties of the light emitter are still to be analyzed. Indeed, the light emitter, oxyluciferin, can be in six different forms due to interexchange reactions. In this work, we present the simulation of the absorption and emission spectra of the possible natural oxyluciferin forms in water and some of their analogues considering both the solvent/oxyluciferin interactions and the dynamical effects by using MD simulations and QM/MM methods. On the one hand, the absorption band shapes have been rationalized by analyzing the electronic nature of the transitions involved. On the other hand, the simulated and experimental emission spectra have been compared. In this case, an ultrafast excited state proton transfer (ESPT) occurs in oxyluciferin and its analogues, which impairs the detection of the emission from the protonated state by steady-state fluorescence spectroscopy. Transient absorption spectroscopy was used to evidence this ultrafast ESPT and rationalize the comparison between simulated and experimental steady-state emission spectra. Finally, this work shows the suitability of the studied oxyluciferin analogues to mimic the corresponding natural forms in water solution, as an elegant way to block the desired interexchange reactions allowing the study of each oxyluciferin form separately.