Combined influence of the QM methods, active space size, Franck–Condon approximation, Herzberg–Teller effect and Duschinsky effect on vibrationally resolved electronic spectra: insights from firefly oxyluciferin

Theoretical simulations of vibrationally resolved electronic spectra not only allow direct comparison with experimental results but also enable simulations under conditions that are experimentally challenging, such as in the gas phase. For instance, firefly oxyluciferin (OL) exists in six neutral and ionic forms in equilibrium within a physiologically relevant pH range in aqueous solutions. Investigating the absorption and fluorescence spectra of OL helps in elucidating the mechanism of firefly bioluminescence and reveals its photophysical properties, thereby enabling control of the emission color. Recent experiments have investigated its absorption and fluorescence spectra in the gas phase, revealing that oxyluciferin anion (OL<sup>-</sup>) can exist in both phenolate-enol and phenolate-keto forms. In this study, we present a simulation of the vibrationally resolved absorption and fluorescence spectra of four OL<sup>-</sup> derivatives, namely, phenolate-enol, phenolate-keto, meo-OL<sup>-</sup> (enol-locked) and dm-OL<sup>-</sup> (keto-locked). Our spectra simulation employs several electronic structure methods, including time-dependent density functional theory (TDDFT), complete active space second-order perturbation theory (CASPT2), and density matrix renormalization group (DMRG) CASPT2. In addition, we studied the influence of quantum mechanics (QM) methods, active space size, the Herzberg-Teller effect, and the Duschinsky effect on the spectra simulation. The simulated spectra not only agree well with experimental results but also detail the experimental findings, revealing key vibronic transitions of OL<sup>-</sup> derivatives. The combination of CASPT2 and TDDFT calculated that the 0-0 transition energy can be directly used as a reference for the band origin when it is not reported experimentally. Our findings can also be extended to simulate the absorption and fluorescence spectra of other semirigid organic molecules.

Authors

Institutions

Publication Details

Journal
Physical Chemistry Chemical Physics
Published
2025-12-02
DOI
https://doi.org/10.1039/d5cp04222d
Primary Topic
bioluminescence and chemiluminescence research
Type
article
Field-Weighted Citation Impact
0.00
Controls
|||
ALL TIME
JAN
FEB
MAR
APR
MAY
JUN
JUL
AUG
SEP
OCT
article

Combined influence of the QM methods, active space size, Franck–Condon approximation, Herzberg–Teller effect and Duschinsky effect on vibrationally resolved electronic spectra: insights from firefly oxyluciferin

Peng Cheng, Ya‐Jun Liu, Haibo Ma, S.-Y. Pi
Physical Chemistry Chemical Physics
bioluminescence and chemiluminescence research
article

Combined influence of the QM methods, active space size, Franck–Condon approximation, Herzberg–Teller effect and Duschinsky effect on vibrationally resolved electronic spectra: insights from firefly oxyluciferin

Peng Cheng, Ya‐Jun Liu, Haibo Ma, S.-Y. Pi
article en

Abstract

Theoretical simulations of vibrationally resolved electronic spectra not only allow direct comparison with experimental results but also enable simulations under conditions that are experimentally challenging, such as in the gas phase. For instance, firefly oxyluciferin (OL) exists in six neutral and ionic forms in equilibrium within a physiologically relevant pH range in aqueous solutions. Investigating the absorption and fluorescence spectra of OL helps in elucidating the mechanism of firefly bioluminescence and reveals its photophysical properties, thereby enabling control of the emission color. Recent experiments have investigated its absorption and fluorescence spectra in the gas phase, revealing that oxyluciferin anion (OL<sup>-</sup>) can exist in both phenolate-enol and phenolate-keto forms. In this study, we present a simulation of the vibrationally resolved absorption and fluorescence spectra of four OL<sup>-</sup> derivatives, namely, phenolate-enol, phenolate-keto, meo-OL<sup>-</sup> (enol-locked) and dm-OL<sup>-</sup> (keto-locked). Our spectra simulation employs several electronic structure methods, including time-dependent density functional theory (TDDFT), complete active space second-order perturbation theory (CASPT2), and density matrix renormalization group (DMRG) CASPT2. In addition, we studied the influence of quantum mechanics (QM) methods, active space size, the Herzberg-Teller effect, and the Duschinsky effect on the spectra simulation. The simulated spectra not only agree well with experimental results but also detail the experimental findings, revealing key vibronic transitions of OL<sup>-</sup> derivatives. The combination of CASPT2 and TDDFT calculated that the 0-0 transition energy can be directly used as a reference for the band origin when it is not reported experimentally. Our findings can also be extended to simulate the absorption and fluorescence spectra of other semirigid organic molecules.

Physical Chemistry Chemical PhysicsVol. 28(2)
Shandong University (CN), Qingdao Binhai University (CN), Beijing Normal University (CN)
Openalex Percentile: Top 63%
bioluminescence and chemiluminescence research
AI Navigator

Ask Laika to Summarize, Analyze, and Connect papers live on the map.

Summarize Papers & Methodologies

Extract key findings, datasets, and comparative methods across publications.

Benchmark Rankings & Visual Analytics

Rank top research institutions, authors, funders, topics, and journals by Field-Weighted Citation Impact (FWCI) and paper volume with instant charts.

Connect Distant Disciplines

Bridge topological clusters on the map to find hidden collaborative intersections.