Theoretical study of cyclohexadiene/hexatriene photochemical interconversion using spin-Flip time-Dependent density functional theory

Deciphering the Internal Conversion Processes Involved in the Photochemical Ring-Opening of 1,3-Cyclohexadiene by Symmetric sp2-Carbon Substitutions

Chemical substituents hold the potential to markedly influence the photochemical behavior in molecular systems and assist in gaining a comprehensive understanding of nonadiabatic phenomena. In this study, we have conducted a comparative analysis of the influence of chemical substituents on the photochemical ring-opening of 1,3-cyclohexadiene (CHD), considering four systems: CHD, 2,3-dimethylcyclohexadiene (CHD-Me2-1), 1,4-dimethylcyclohexadiene (CHD-Me2-2), and 1,2,3,4-tetramethylcyclohexadiene (CHD-Me4), using electronic structure theory calculations and nonadiabatic molecular dynamics simulations. Employing extended multistate complete active space second-order perturbation (XMS-CASPT2) theory, we optimized reactants, S1 states, conical intersections (CIns), and products, revealing structural and energetic variations consistent with prior research. Nonadiabatic molecular dynamics simulation was used to gain insights into photochemical dynamics at state-averaged complete active space self-consistent field (SA-CASSCF) theory. CHD-Me4 exhibited reduced carbon-carbon single bond rupture rates, responsible for ring-opening, due to substituent proximity. Further, CHD-Me2-2 and CHD-Me4 displayed prolonged excited-state relaxation times, highlighting notable substituents' impact. Analysis of kinetic energy profiles of specific carbon atoms also revealed restrained atomic displacements, particularly in CHD-Me2-2 and CHD-Me4. These findings advance our understanding of how substituents modulate photochemical reactions in cyclohexadiene derivatives, guiding new molecular design and future research.

Ab Initio Investigation of Intramolecular Charge Transfer States in DMABN by Calculation of Excited State X-ray Absorption Spectra

Dual fluorescence in 4-(dimethylamino)benzonitrile (DMABN) and its derivatives in polar solvents has been studied extensively for the past several decades. An intramolecular charge transfer (ICT) minimum on the excited state potential energy surface, in addition to the localized low-energy (LE) minimum, has been proposed as a mechanism for this dual fluorescence, with large geometric relaxation and molecular orbital reorganization a key feature of the ICT pathway. Herein, we have used both equation-of-motion coupled-cluster with single and double excitations (EOM-CCSD) and time-dependent density functional (TDDFT) methods to investigate the landscape of excited state potential energy surfaces across a number of geometric conformations proposed as ICT structures. In order to correlate these geometries and valence excited states in terms of potential experimental observables, we have calculated the nitrogen K-edge ground and excited state absorption spectra for each of the predicted "signpost" structures and identified several key spectral features that could be used to interpret a future time-resolved X-ray absorption experiment.

Photochemistry of Thymine in Solution and DNA Revealed by an Electrostatic Embedding QM/MM Combined with Mixed-Reference Spin-Flip TDDFT

The photochemistry of nucleobases, important for their role as building blocks of DNA, is largely affected by the electrostatic environment in which they are soaked. For example, despite the numerous studies of thymine in solution and DNA, there is still a debate on the photochemical deactivation pathways after UV absorption. Many theoretical models are oversimplified due to the lack of computationally accurate and efficient electronic structure methodologies that capture excited state electron correlation effects when nucleobases are embedded in large electrostatic media. Here, we combine mixed-reference spin-flip time-dependent density functional theory (MRSF-TDDFT) with electrostatic embedding QM/MM using electrostatic potential fittingfitted (ESPF) atomic charges, as a strategy to accurately and efficiently describe the electronic structure of chromophores polarized by an electrostatic medium. In particular, we develop analytic expressions for the energy and gradient of MRSF/MM based on the ESPF coupling using atom-centered grids and total charge conservation. We apply this methodology to the study of solvation effects on thymine photochemistry in water and thymine dimers in DNA. In the former, the combination of trajectory surface hopping (TSH) nonadiabatic molecular dynamics (NAMD) with MRSF/MM remarkably revealed accelerated deactivation decay pathways, which is consistent with the experimental decay time of ∼400 fs. The enhanced hopping rate can be explained by the preferential stabilization of corresponding conical interactions due to their increased dipole moments. Structurally, it is a consequence of characteristic methyl puckered geometries near the conical intersection region. For the thymine dimer in B-DNA, we found new photochemical pathways through conical intersections that could explain the formation of cyclobutadiene dimers and 6-4 photoproducts.

Analytical gradients and derivative couplings for the TDDFT-1D method

We derive and implement analytic gradients and derivative couplings for time-dependent density functional theory plus one double (TDDFT-1D) which is a semiempirical configuration interaction method whereby the Hamiltonian is diagonalized in a basis of all singly excited configurations and one doubly excited configuration as constructed from a set of reference Kohn-Sham orbitals. We validate the implementation by comparing against finite difference values. Furthermore, we show that our implementation can locate both optimized geometries and minimum-energy crossing points along conical seams of S₁/S₀ surfaces for a set of test cases.

Nonadiabatic Dynamics of 1,3-Cyclohexadiene by Curvature-Driven Coherent Switching with Decay of Mixing

The photoinduced ring-opening reaction of 1,3-cyclohexadiene to produce 1,3,5-hexatriene is a classic electrocyclic reaction and is also a prototype for many reactions of biological and synthetic importance. Here, we simulate the ultrafast nonadiabatic dynamics of the reaction in the manifold of the three lowest valence electronic states by using extended multistate complete-active-space second-order perturbation theory (XMS-CASPT2) combined with the curvature-driven coherent switching with decay of mixing (κCSDM) dynamical method. We obtain an excited-state lifetime of 79 fs, and a product quantum yield of 40% from the 500 trajectories initiated in the S₁ excited state. The obtained lifetime and quantum yield values are very close to previously reported experimental and computed values, showing the capability of performing a reasonable nonadiabatic ring-opening dynamics with the κCSDM method that does not require nonadiabatic coupling vectors, time derivatives, or diabatization. In addition, we study the ring-opening reaction by initiating the trajectories in the dark state S₂. We also optimize the S₀/S₁ and S₁/S₂ minimum-energy conical intersections (MECIs) by XMS-CASPT2; for S₁/S₂, we optimized both an inner and an outer local-minimum-energy conical intersections (LMECIs). We provide the potential energy profile along the ring-opening coordinate by joining selected critical points via linear synchronous transit paths. We find the inner S₁/S₂ LMECI to be more crucial than the outer one.

Providing theoretical insight into the role of symmetry in the photoisomerization mechanism of a non-symmetric dithienylethene photoswitch

Dithienylethene (DTE) molecular photoswitches have shown to be excellent candidates in the design of efficient optoelectronic devices, due to their high photoisomerization quantum yield (QY), for which symmetry is suggested to play a crucial role. Here, we present a theoretical study on the photochemistry of a non-symmetric dithienylethene photoswitch, with a special emphasis on the effect of asymmetric substitution on the photocyclization and photoreversion mechanisms. We used the Spin-Flip Time Dependent Density Functional Theory (SF-TDDFT) method to locate and characterize the main structures (conical intersections and minima) of the ground state and the first two excited states, S₁ and S₂, along the ring-opening/closure reaction coordinate of the photocyclization and photoreversion processes, and to identify the important coordinates governing the radiationless decay pathways. Our results suggest that while the main features that characterize the photoisomerization of symmetric DTEs are also present for the photoisomerization of the non-symmetric DTE, the lower energy barrier on S₁ along the cycloreversion reaction speaks in favor of a more efficient and therefore a higher cycloreversion QY for the non-symmetric DTEs, making them a better candidate for molecular optoelectronic devices than their symmetric counterparts.

On the inclusion of one double within CIS and TDDFT

We present an improved approach for generating a set of optimized frontier orbitals (HOMO and LUMO) that minimizes the energy of one double configuration. We further benchmark the effect of including such a double within a rigorous configuration interaction singles or a parameterized semi-empirical time-dependent density functional theory Hamiltonian for a set of test cases. Although we cannot quite achieve quantitative accuracy, the algorithm is quite robust and routinely delivers an enormous qualitative improvement to standard single-reference electronic structure calculations.

Nonadiabatic dynamics with spin-flip vs linear-response time-dependent density functional theory: A case study for the protonated Schiff base C5H6NH2. J Chem Phys 2021;155:124111. [PMID: 34598550 DOI: 10.1063/5.0062757]

Nonadiabatic trajectory surface hopping simulations are reported for trans-C₅H₆NH₂ ⁺, a model of the rhodopsin chromophore, using the augmented fewest-switches algorithm. Electronic structure calculations were performed using time-dependent density functional theory (TDDFT) in both its conventional linear-response (LR) and its spin-flip (SF) formulations. In the SF-TDDFT case, spin contamination in the low-lying singlet states is removed by projecting out the lowest triplet component during iterative solution of the TDDFT eigenvalue problem. The results show that SF-TDDFT qualitatively describes the photoisomerization of trans-C₅H₆NH₂ ⁺, with favorable comparison to previous studies using multireference electronic structure methods. In contrast, conventional LR-TDDFT affords qualitatively different photodynamics due to an incorrect excited-state potential surface near the Franck-Condon region. In addition, the photochemistry (involving pre-twisting of the central double bond) appears to be different for SF- and LR-TDDFT, which may be a consequence of different conical intersection topographies afforded by these two methods. The present results contrast with previous surface-hopping studies suggesting that the LR-TDDFT method's incorrect topology around S₁/S₀ conical intersections is immaterial to the photodynamics.

Theoretical investigation of a novel xylene-based light-driven unidirectional molecular motor

In this study, the working mechanism of the first light-driven rotary molecular motors used to control an eight-base-pair DNA hairpin has been investigated. In particular, this linker was reported to have promising photophysical properties under physiological conditions, which motivated our work at the quantum mechanical level. Cis-trans isomerization is triggered by photon absorption at wavelengths ranging 300 nm-400 nm, promoting the rotor to the first excited state, and it is mediated by an energy-accessible conical intersection from which the ground state is reached back. The interconversion between the resulting unstable isomer and its stable form occurs at physiological conditions in the ground state and is thermally activated. Here, we compare three theoretical frameworks, generally used in the quantum description of medium-size chemical systems: Linear-Response Time-Dependent Density Functional Theory (LR-TDDFT), Spin-Flip TDDFT (SF-TDDFT), and multistate complete active space second-order perturbation theory on state-averaged complete active space self consistent field wavefunctions (MS-CASPT2//SA-CASSCF). In particular, we show the importance of resorting to a multireference approach to study the rotational cycle of light-driven molecular motors due to the occurrence of geometries described by several configurations. We also assess the accuracy and computational cost of the SF-TDDFT method when compared to MS-CASPT2 and LR-TDDFT.

Analytic gradients and derivative couplings for configuration interaction with all single excitations and one double excitation-En route to nonadiabatic dynamics

We present analytic gradients and derivative couplings for the simplest possible multireference configuration interaction method, CIS-1D, an electronic structure Ansatz that includes all single excitations and one lone double excitation on top of a Hartree-Fock reference state. We show that the resulting equations are numerically stable and require the evaluation of a similar number of integrals as compared to standard CIS theory; one can easily differentiate the required frontier orbitals (h and l) with minimal cost. The resulting algorithm has been implemented within the Q-Chem electronic structure package and should be immediately useful for understanding photochemistry with S₀-S₁ crossings.