The photoinduced isomerization mechanism of the 2-(1-(methylimino)methyl)-6-chlorophenol (SMAC): Nonadiabatic surface hopping dynamics simulations

Excited-state dynamics of 3-hydroxychromone in gas phase

In recent years, 3-hydroxychromone (3-HC) and its derivatives have attracted much interest for their applications as molecular photoswitches and fluorescent probes. A clear understanding of their excited-state dynamics is essential for their applications and further development of new functional 3-HC derivatives. However, the deactivation mechanism of the photoexcited 3-HC family is still puzzling as their spectral properties are sensitive to the surrounding medium and substituents. The excited-state relaxation channels of 3-HC have been a matter of intense debate. In the current work, we thoroughly investigated the excited-state decay process of the 3-HC system in the gas phase using high-level electronic structure calculations and on-the-fly excited-state dynamic simulations intending to provide insight into the intrinsic photochemical properties of the 3-HC system. A new deactivation mechanism is proposed in the gas phase, which is different from that in solvents. The excited-state intramolecular proton transfer (ESIPT) process that occurs in solutions is not preferred in the gas phase due to the existence of a sizable energy barrier (∼0.8 eV), and thus, no dual fluorescence is found. On the contrary, the non-radiative decay process is the dominant decay channel, which is driven by photoisomerization combined with ring-puckering and ring-opening processes. The results coincide with the observations of an experiment performed in a supersonic jet by Itoh (M. Itoh, Pure Appl. Chem., 1993, 65(8), 1629-1634). The current work indicates that the solution environment plays an important role in regulating the excited-state dynamic behaviour of the 3-HC system. This study thus provides theoretical guidance for the rational design and improvement of the photochemical properties of the 3-HC system and paves the way for further investigation into its photochemical properties in complex environments.

Revealing the excited-state dynamics of cytidine and the role of excited-state proton transfer process

The high photostability of DNAs and RNAs is inextricably related to the photochemical and photophysical properties of their building blocks, nucleobases and nucleosides, which can dissipate the absorbed UV light energy in a harmless manner. The deactivation mechanism of the nucleosides, especially the decay pathways of cytidine (Cyd), has been a matter of intense debate. In the current study, we employ high-level electronic structure calculations combined with excited state non-adiabatic dynamic simulations to provide a clear picture of the excited state deactivation of Cyd in both gas phase and aqueous solution. In both environments, a barrierless decay path driven by the ring-puckering motion and a relaxation channel with a small energy barrier driven by the elongation motion of CO bond are assigned to <200 fs and sub-picosecond decay time component, respectively. The presence of ribose group has a subtle effect on the dynamic behavior of Cyd in gas phase as the ribose-to-base hydrogen/proton transfer process is energetically inaccessible with a sizable energy barrier of about 1.4 eV. However, this energy barrier is significantly reduced in water, especially when an explicit water molecule is present. Therefore, we argue that the long-lived decay channel found in aqueous solution could be assigned to the Cyd-water intermolecular hydrogen/proton transfer process. The present study postulates a novel scenario toward deep understanding the intrinsic photostability of DNAs and RNAs and provides solid evidence to disclose the long history debate of cytidine excited-state decay mechanism, especially for the assignment of experimentally observed time components.

Position of the Benzene Ring Substituent Regulates the Excited-State Deactivation Process of the Benzyluracil Systems

A systematic theoretical study of the regulating effect of the substituent position on the photoinduced deactivation process of the benzyluracil systems has been performed based on the high-level static electronic structure calculations and on-the-fly full-dimensional excited-state dynamics simulations. Similarities and differences coexist for the two systems by comparative studies on the photoinduced deactivation process of the 5-benzyluracil (5-BU) and 6-benzyluracil (6-BU) systems. They both obey an S₂ → S₁ → S₀ two-step decay pattern, and the decay coordinates of the S₂ → S₁ and S₁ → S₀ processes are mainly driven by the elongation of the bridging bond and the out-of-plane ring deformation motion, respectively. However, the puckering motion occurring at the C2 atom in the uracil fragment dominates the decay pathway of the 5-BU system. On the contrary, the puckering motion at the C5 atom in the benzene fragment mainly drives the decay coordinate of the 6-BU system. Therefore, the substituent position could play significant roles in the deactivation process of the benzyluracil systems. Moreover, the S₁ → S₀ decay process of the 6-BU system consists of five pathways, possessing a more complex deactivation picture than the 5-BU system. The fitted time scale of the puckering motion is compatible with the experimentally observed lifetimes. This work provides a fundamental understanding of the photophysical and photochemical properties of the benzyluracil systems and can give rational suggestions to further design or regulate the bionic molecular systems.

Comparison of Spin-Flip TDDFT-Based Conical Intersection Approaches with XMS-CASPT2

Determining conical intersection geometries is of key importance to understanding the photochemical reactivity of molecules. While many small- to medium-sized molecules can be treated accurately using multireference approaches, larger molecules require a less computationally demanding approach. In this work, minimum energy crossing point conical intersection geometries for a series of molecules have been studied using spin-flip TDDFT (SF-TDDFT), within the Tamm-Dancoff Approximation, both with and without explicit calculation of nonadiabatic coupling terms, and compared with both XMS-CASPT2 and CASSCF calculated geometries. The less computationally demanding algorithms, which do not require explicit calculation of the nonadiabatic coupling terms, generally fare well with the XMS-CASPT2 reference structures, while the relative energetics are only reasonably replicated with the MECP structure as calculated with the BHHLYP functional and full nonadiabatic coupling terms. We also demonstrate that, occasionally, CASSCF structures deviate quantitatively from the XMS-CASPT2 structures, showing the importance of including dynamical correlation.

The effect of surfactants on hydrate particle agglomeration in liquid hydrocarbon continuous systems: a molecular dynamics simulation study

Anti-agglomerants (AAs), both natural and commercial, are currently being considered for gas hydrate risk management of petroleum pipelines in offshore operations. However, the molecular mechanisms of the interaction between the AAs and gas hydrate surfaces and the prevention of hydrate agglomeration remain critical and complex questions that need to be addressed to advance this technology. Here, we use molecular dynamics (MD) simulations to investigate the effect of model surfactant molecules (polynuclear aromatic carboxylic acids) on the agglomeration behaviour of gas hydrate particles and disruption of the capillary liquid bridge between hydrate particles. The results show that the anti-agglomeration pathway can be divided into two processes: the spontaneous adsorption effect of surfactant molecules onto the hydrate surface and the weakening effect of the intensity of the liquid bridge between attracted hydrate particles. The MD simulation results also indicate that the anti-agglomeration effectiveness of surfactants is determined by the intrinsic nature of their molecular functional groups. Additionally, we find that surfactant molecules can affect hydrate growth, which decreases hydrate particle size and correspondingly lower the risk of hydrate agglomeration. This study provides molecular-level insights into the anti-agglomeration mechanism of surfactant molecules, which can aid in the ultimate application of natural or commercial AAs with optimal anti-agglomeration properties.

Schematic of anti-agglomeration effect of surfactants promoting gas hydrate particle dispersion.