Understanding the Excited-State Relaxation Mechanisms of Xanthophyll Lutein by Multi-configurational Electronic Structure Calculations

Abstract

The contradictory behaviors in light harvesting and non-photochemical quenching make xanthophyll lutein the most attractive functional molecule in photosynthesis. Despite several theoretical simulations on the spectral properties and excited-state dynamics, the atomic-level photophysical mechanisms need to be further studied and established, especially for an accurate description of geometric and electronic structures of conical intersections for the lowest several electronic states of lutein. In the present work, semiempirical OM2/MRCI and multi-configurational restricted active space self-consistent field methods were performed to optimize the minima and conical intersections in and between the 1Ag − , 2Ag − , 1Bu + , and 1Bu − states. Meanwhile, the relative energies were refined by MS-CASPT2(10,8)/6-31G*, which can reproduce correct electronic state properties as those in the spectroscopic experiments. Based on the above calculation results, we proposed a possible excitedstate relaxation mechanism for lutein from its initially populated 1Bu + state. Once excited to the optically bright 1Bu + state, the system will propagate along the key reaction coordinate, i.e., the stretching vibration of the conjugated carbon chain. During this period of time, the 1Bu − state will participate in and forms a resonance state between the 1Bu − and 1Bu + states. Later, the system will rapidly hop to the 2Ag − state via the 1Bu + /2Ag − conical intersection. Finally, the lutein molecule will survive in the 2Ag − state for a relatively long time before it internally converts to the ground state directly or via a twisted S 1 /S 0 conical intersection. Notably, though the photophysical picture may be very different in solvents and proteins, the current theoretical study proposed a promising calculation protocol and also provided many valuable mechanistic insights for lutein and similar carotenoids.

Construction of Highly Accurate Machine Learning Potential Energy Surfaces for Excited-State Dynamics Simulations Based on Low-Level Data Sets

“…Many successful applications, such as the study on the photochromic mechanism of bridged diarylethene, the prediction of electron ionization mass spectra, and the simulation of ultrafast intersystem crossing processes, have been reported previously. − However, the impact of electronic structure methods on the on-the-fly nonadiabatic dynamics may be strongly dependent on the simulated systems. On the one hand, Shen and co-workers recently reported that the excited-state relaxation mechanisms of xanthophyll lutein can be well captured using OM2/MRCI . They also observed that the nonadiabatic dynamics results of cis -azobenzene photoisomerization at the OM2/MRCI level are comparable with experimental values even though the most stable isomer of azobenzene in the ground state (S 0 ) is misidentified. , On the other hand, sensitivity analysis of nonadiabatic dynamics to electronic structure methods reveals that the importance of high-quality PESs is much more significant than expected before, − which queries the reliability of OM2/MRCI calculations on photochemical systems.…”