Abstract:The use of flavins and flavoproteins in photocatalytic, sensing, and biotechnological applications has led to a growing interest in computationally modeling the excited state electronic structure and photophysics of flavin. However, there is limited consensus regarding which computational methods are appropriate for modeling flavin’s photophysics. We compare the energies of low-lying excited states of flavin computed with time-dependent density functional theory (TD-DFT), equation of motion coupled cluster (EOM-EE-CCSD), scaled opposite-spin configuration interaction (SOS-CIS(D)), multiconfiguration pair-density functional theory (MC-PDFT), and several multi-reference perturbation theory (MR-PT2) methods. In the first part, we focus on excitation energies of the first singlet excited state (S1) of five different redox and protonation states of flavin, with the goal of finding a suitable active space for MR-PT2 calculations. In the second part, we construct two sets of one-dimensional potential energy surfaces connecting the S0 and S1 equilibrium geometries (S0–S1 path) and the S1 (π,π*) and S2 (n,π*) equilibrium geometries (S1–S2 path). The first path therefore follows a Franck-Condon active mode of flavin while the second path maps crossings points between low-lying singlet and triplet states in flavin. We discuss the similarities and differences in the TD-DFT, EOM-EE-CCSD, SOS-CIS(D), MC-PDFT and MR-PT2 energy profiles along these paths. We find that (TD-)DFT methods are suitable for applications such as simulating the spectra of flavins but are inconsistent with several other methods when used for some geometry optimizations and when describing the energetics of dark (n,π*) states. MR-PT2 methods show promise for the simulation of flavin’s low-lying excited states, but the selection of orbitals for the active space and in the number of roots used for state averaging must be done carefully to avoid artefacts. Some properties, such as the intersystem crossing geometry and energy between the S1 (π,π*) and T2 (n,π*) states, may require additional benchmarking before they can be determined quantitatively.

Calculating Absorption and Fluorescence Spectra for Chromophores in Solution with Ensemble Franck-Condon Methods

Combining the Ensemble and Franck-Condon Approaches for Spectral Shapes of Molecules in Solution

Optical spectra in the condensed phase: Capturing anharmonic and vibronic features using dynamic and static approaches

Electronic Structure Methods for Simulating Flavin's Spectroscopy and Photophysics: Comparison of Multi-reference, TD-DFT, and Single-Reference Wave Function Methods

Electronic Structure Methods for Simulating Flavin’s Spectroscopy and Photophysics: Comparison of Multi-Reference, Single-Reference, and TD-DFT Methods

Multiscale Modeling of Electronic Spectra Including Nuclear Quantum Effects

Accelerated Computer-Aided Screening of Optical Materials: Investigating the Potential of Δ-SCF Methods to Predict Emission Maxima of Large Dye Molecules

Simulation of Condensed-Phase Spectroscopy with Near-Term Digital Quantum Computers.

Exploring the Influence of Approximations for Simulating Valence Excited X-ray Spectra

Explicit solvation modeling in the accurate TD-DFT prediction of absorption spectra for natural nucleobases and fluorescent nucleobase analogues

First-principles Modeling of the Absorption spectrum of Crystal Violet in Solution - The Importance of Environmentally-driven Symmetry Breaking

Influence of non-adiabatic effects on linear absorption spectra in the condensed phase: Methylene blue

Predicting fluorescence quantum yields for molecules in solution: A critical assessment of the harmonic approximation and the choice of the lineshape function

Nonlinear spectroscopy in the condensed phase: The role of Duschinsky rotations and third order cumulant contributions

Explicit Modelling of Spectral Bandshapes by a Mixed Quantum‐Classical Approach: Solvent Order and Temperature Effects in the Optical Spectra of Distryrylbenzene

Anharmonic Franck-Condon Simulation of the Absorption and Fluorescence Spectra for the Low-Lying S-1 and S-2 Excited States of Pyridine

Forster resonance energy transfer, absorption and emission spectra in multichromophoric systems: III. Exact stochastic path integral evaluation

Exploring the Mechanisms Behind Non Aromatic Fluorescence with Density Functional Tight Binding Method

Quantum-Classical Protocol for Efficient Characterization of Absorption Lineshape and Fluorescence Quenching upon Aggregation: The Case of Zinc Phthalocyanine Dyes

An analytical continuation approach for evaluating emission lineshapes of molecular aggregates and the adequacy of multichromophoric Förster theory