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

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 multireference 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 the number of roots used for state averaging must be done carefully to avoid artifacts. 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.

Hydrogen-Bond Network Determines the Early Photoisomerization Processes of Cph1 and AnPixJ Phytochromes

Phytochrome proteins are light receptors that play a pivotal role in regulating the life cycles of plants and microorganisms. Intriguingly, while cyanobacterial phytochrome Cph1 and cyanobacteriochrome AnPixJ use the same phycocyanobilin (PCB) chromophore to absorb light, their excited‐state behavior is very different. We employ multiscale calculations to rationalize the different early photoisomerization mechanisms of PCB in Cph1 and AnPixJ. We found that their electronic S₁, T₁, and S₀ potential minima exhibit distinct geometric and electronic structures due to different hydrogen bond networks with the protein environment. These specific interactions influence the S₁ electronic structures along the photoisomerization paths, ultimately leading to internal conversion in Cph1 but intersystem crossing in AnPixJ. This explains why the excited‐state relaxation in AnPixJ is much slower (ca. 100 ns) than in Cph1 (ca. 30 ps). Further, we predict that efficient internal conversion in AnPixJ can be achieved upon protonating the carboxylic group that interacts with PCB.

Computer-Assisted Design of Environmentally Friendly and Light-Stable Fluorescent Dyes for Textile Applications

Five potentially environmentally friendly and light-stable hemicyanine dyes were designed based on integrated consideration of photo, environmental, and computational chemistry as well as textile applications. Two of them were synthesized and applied in dyeing polyacrylonitrile (PAN), cotton, and nylon fabrics, and demonstrated the desired properties speculated by the programs. The computer-assisted analytical processes includes estimation of the maximum absorption and emission wavelengths, aquatic environmental toxicity, affinity to fibers, and photo-stability. This procedure could effectively narrow down discovery of new potential dye structures, greatly reduce and prevent complex and expensive preparation processes, and significantly improve the development efficiency of novel environmentally friendly dyes.