Surface Hopping Dynamics Research Articles

Simulation of the non-adiabatic dynamics of an enone-Lewis acid complex in an explicit solvent.

Unlocking the full potential of Lewis acid catalysis for photochemical transformations requires a comprehensive understanding of the ultrafast dynamics of substrate-Lewis acid complexes. In a previous article [Peschel et al., Angew. Chem. Int. Ed., 2021, 60, 10155], time-resolved spectroscopy supported by static calculations revealed that the Lewis acid remains attached during the relaxation of the model complex cyclohexenone-BF3. In contrast to the experimental observation, surface-hopping dynamics in the gas phase predicted ultrafast heterolytic dissociation. We attributed the discrepancy to missing solvent interactions. Thus, in this work, we present an interface between the SHARC and FermiONs++ program packages, which enables us to investigate the ultrafast dynamics of cyclohexenone-BF3 in an explicit solvent environment. Our simulations demonstrate that the solvent prevents the dissociation of the complex, leading to an intriguing dissociation-reassociation mechanism. Comparing the dynamics with and without triplet states highlights their role in the relaxation process and shows that the Lewis acid inhibits intersystem crossing. These findings provide a clear picture of the relaxation process, which may aid in designing future Lewis acid catalysts for photochemical applications. They underscore that an explicit solvent model is required to describe relaxation processes in weakly bound states, as energy transfer to the solvent is crucial for the system to reach its minimum geometries.

Ultrafast photochemistry and electron diffraction for cyclobutanone in the S2 state: Surface hopping with time-dependent density functional theory.

We simulate the photodynamics of gas-phase cyclobutanone excited to the S2 state using fewest switches surface hopping (FSSH) dynamics powered by time-dependent density functional theory (TDDFT). We predict a total photoproduct yield of 8%, with a C3:C2 product ratio of 0 trajectories to 8 trajectories. One primary S2 → S1 conical intersection is identified involving the compression of an α-carbon-carbon-hydrogen bond angle. Excited state lifetimes computed with respect to electronic state populations were found to be 3.96ps (S2 → S1) and 498fs (S1 → S0). We also generate time-resolved difference pair distribution functions (ΔPDFs) from our TDDFT-FSSH dynamics results in order to generate direct comparisons with ultrafast electron diffraction experiment observables. Global and target analysis of time-resolved ΔPDFs produced a distinct set of lifetimes: (i) a 0.548ps decay and (ii) a 1.69ps decay, both resembling the S2 minimum, as well as (iii) a long decay that resembles the S1 minimum geometry and the fully separated C2 products. Finally, we contextualize our results by considering the impact of the most likely sources of significant errors.

A different photochromic mechanism of spirooxadiazine: Electronic structure calculations and nonadiabatic dynamics simulations

Previously, we conducted electronic structure calculations and surface hopping dynamics to investigate the photochromic mechanism of spirooxazine (SO). Herein, a very similar calculation strategy was employed to simulate the excited-state ring-opening, relaxation, and subsequent isomerization of spirooxadiazine (SOD). Interestingly, we discovered a significantly different photochromic behavior compared to SO. Firstly, the bright state of SOD is the S2(ππ*) state, while the S1 state is a weakly absorbing state with a nπ* character. This is in stark contrast to SO, of which the S1 state is a ππ* state. Secondly, we identified a new excited-state decay channel along the C2N3 dissociation within the benzoxadiazine ring, in addition to the closed-loop and C2O7 cleavage relaxation pathways. Thirdly, the out-of-plane hydrogen (HOOP) oscillations were no longer operative in the relaxation processes, as the CH group was replaced by the N atom in SOD. Finally, several nonadiabatic trajectories will lead to the ‘N2 release’ product in the ground state, which might result in deactivation of SOD and affect the light fatigue resistance.

How Do the Position and Number of Methyl Substituents Affect the Photochemical Process of Criegee Intermediate? Trajectory Surface-Hopping Dynamics of Four-Carbon CIs.

Electronic-structure calculations combined with nonadiabatic trajectory surface-hopping (TSH) dynamic simulations were carried out on two alkenyl-substituted Criegee intermediates (CIs), i.e., propenyl-substituted CI (PCI) and 1-methyl-propenyl substituted CI (MPCI), in order to investigate the influence of the position and number of substituents on the photochemical process of CI in S1 states. It is found that they play critical roles in the reactivity, dominant product channel, and mechanism of the CIs. More specifically, introducing a methyl group on either C1 (α-C) or C3 (γ-C) position of a vinyl-substituted CI (VCI) skeleton facilitates the rotation of the C1═O1 bond and leads to the formation of a three-membered dioxirane ring; meanwhile, it evidently enhances the reactively of the S1-state molecule. Meanwhile, methyl substitution on the vinyl moiety [i.e., C2 (β-C) and C3 (γ-C) positions] is beneficial for the rotation of the C2═C3 bond and thus facilitates the formation of the five-membered 1,2-dioxole ring, and the substitution on C2 site decreases the reactivity. The cosubstitution of C2 and C3 atoms by methyl groups well balances the features of VCI in the sense of high reactivity, consistently predominant channel, and possible dioxole side-product. The findings here not only deepen the knowledge on the photochemical processes of the CI but also inspire the rethinking of the "old" concept of substitution effect.

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

Machine learning is capable of effectively predicting the potential energies of molecules in the presence of high-quality data sets. Its application in the construction of ground- and excited-state potential energy surfaces is attractive to accelerate nonadiabatic molecular dynamics simulations of photochemical reactions. Because of the huge computational cost of excited-state electronic structure calculations, the construction of a high-quality data set becomes a bottleneck. In the present work, we first built two data sets. One was obtained from surface hopping dynamics simulations at the semiempirical OM2/MRCI level. Another was extracted from the dynamics trajectories at the CASSCF level, which was reported previously. The ground- and excited-state potential energy surfaces of ethylene-bridged azobenzene at the CASSCF computational level were constructed based on the former low-level data set. Although non-neural network machine learning methods can achieve good or modest performance during the training process, only neural network models provide reliable predictions on the latter external test data set. The BPNN and SchNet combined with the Δ-ML scheme and the force term in the loss functions are recommended for dynamics simulations. Then, we performed excited-state dynamics simulations of the photoisomerization of ethylene-bridged azobenzene on machine learning potential energy surfaces. Compared with the lifetimes of the first excited state (S1) estimated at different computational levels, our results on the E isomer are in good agreement with the high-level estimation. However, the overestimation of the Z isomer is unimproved. It suggests that smaller errors during the training process do not necessarily translate to more accurate predictions on high-level potential energies or better performance on nonadiabatic dynamics simulations, at least in the present case.

The development of the QM/MM interface and its application for the on-the-fly QM/MM nonadiabatic dynamics in JADE package: Theory, implementation, and applications.

Understanding the nonadiabatic dynamics of complex systems is a challenging task in computational photochemistry. Herein, we present an efficient and user-friendly quantum mechanics/molecular mechanics (QM/MM) interface to run on-the-fly nonadiabatic dynamics. Currently, this interface consists of an independent set of codes designed for general-purpose use. Herein, we demonstrate the ability and feasibility of the QM/MM interface by integrating it with our long-term developed JADE package. Tailored to handle nonadiabatic processes in various complex systems, especially condensed phases and protein environments, we delve into the theories, implementations, and applications of on-the-fly QM/MM nonadiabatic dynamics. The QM/MM approach is established within the framework of the additive QM/MM scheme, employing electrostatic embedding, link-atom inclusion, and charge-redistribution schemes to treat the QM/MM boundary. Trajectory surface-hopping dynamics are facilitated using the fewest switches algorithm, encompassing classical and quantum treatments for nuclear and electronic motions, respectively. Finally, we report simulations of nonadiabatic dynamics for two typical systems: azomethane in water and the retinal chromophore PSB3 in a protein environment. Our results not only illustrate the power of the QM/MM program but also reveal the important roles of environmental factors in nonadiabatic processes.

The photodissociation dynamics and ultrafast electron diffraction image of cyclobutanone from the surface hopping dynamics simulation

The comprehension of nonadiabatic dynamics in polyatomic systems relies heavily on the simultaneous advancements in theoretical and experimental domains. The gas-phase ultrafast electron diffraction (UED) technique has attracted significant attention as a unique tool for monitoring photochemical and photophysical processes at the all-atomic level with high temporal and spatial resolutions. In this work, we simulate the UED spectra of cyclobutanone using the trajectory surface hopping method at the extended multi-state complete active space second order perturbation theory (XMS-CASPT2) level and thereby predict the results of the upcoming UED experiments in the Stanford Linear Accelerator Laboratory. The simulated results demonstrate that a few pathways, including the C2 and C3 dissociation channels, as well as the ring opening channel, play important roles in the nonadiabatic reactions of cyclobutanone. We demonstrate that the simulated UED signal can be directly interpreted in terms of atomic motions, which provides a unique way of monitoring the evolution of the molecular structure in real time. Our work not only provides numerical data that help to determine the accuracy of the well-known surface hopping dynamics at the high XMS-CASPT2 electronic-structure level but also facilitates the understanding of the microscopic mechanisms of the photoinduced reactions in cyclobutanone.

MLatom Software Ecosystem for Surface Hopping Dynamics in Python with Quantum Mechanical and Machine Learning Methods.

We present an open-source MLatom@XACS software ecosystem for on-the-fly surface hopping nonadiabatic dynamics based on the Landau-Zener-Belyaev-Lebedev algorithm. The dynamics can be performed via Python API with a wide range of quantum mechanical (QM) and machine learning (ML) methods, including ab initio QM (CASSCF and ADC(2)), semiempirical QM methods (e.g., AM1, PM3, OMx, and ODMx), and many types of ML potentials (e.g., KREG, ANI, and MACE). Combinations of QM and ML methods can also be used. While the user can build their own combinations, we provide AIQM1, which is based on Δ-learning and can be used out-of-the-box. We showcase how AIQM1 reproduces the isomerization quantum yield of trans-azobenzene at a low cost. We provide example scripts that, in dozens of lines, enable the user to obtain the final population plots by simply providing the initial geometry of a molecule. Thus, those scripts perform geometry optimization, normal mode calculations, initial condition sampling, parallel trajectories propagation, population analysis, and final result plotting. Given the capabilities of MLatom to be used for training different ML models, this ecosystem can be seamlessly integrated into the protocols building ML models for nonadiabatic dynamics. In the future, a deeper and more efficient integration of MLatom with Newton-X will enable a vast range of functionalities for surface hopping dynamics, such as fewest-switches surface hopping, to facilitate similar workflows via the Python API.

Mechanisms of glycine formation in cold interstellar media: a theoretical study.

The possibility of the formation of glycine (Gly) from fundamental gas molecules in cold interstellar media was studied using quantum chemical methods, transition state theory and microcanonical molecular dynamics simulations with surface hopping dynamics (NVE-MDSH). This theoretical study emphasized five photochemical pathways in the lowest singlet-excited (S 1) state, thermochemical processes after non-radiative S 1→S 0 relaxations, and photo-to-thermal energy conversion in the NVE ensemble. The optimized reaction pathways suggested that to generate a reactive singlet dihydroxy carbene (HOCOH) intermediate, photochemical pathways involving the H2O…CO van der Waals and H2O-OC hydrogen bond precursors (Ch (1)_Step (1)) possess considerably lower energy barriers than the S 0 state pathways. The Gibbs free energy barriers (∆G ǂ ) calculated after the non-radiative S 1 →S 0 relaxations indicated higher spontaneous temperatures (T s) for the formation of the HOCOH intermediate (Ch (1)_Step (1)) than for Gly formation (Ch (1)_Step (2) and Ch (4)). Although the termolecular reaction in Ch (4) possesses a low energy barrier, and is thermodynamically favourable, the high exothermic S 1 →S 0 relaxation energy leads to the separation of the weakly associated H2O…CH2NH…CO complex into single molecules. The NVE-MDSH results also confirmed that the molecular processes after the S 1 →S 0 relaxations are thermally selective, and because the non-radiative S 1 →S 0 relaxation temperatures are exceedingly higher than T s, the formation of Gly on consecutive reaction pathways is non-synergistic with low yields and several side products. Based on the theoretical results, photo-to-thermal control strategies to promote desirable photochemical products are proposed. They could be used as guidelines for future theoretical and experimental research on photochemical reactions.

Prediction Challenge: Simulating Rydberg photoexcited cyclobutanone with surface hopping dynamics based on different electronic structure methods.

This research examines the nonadiabatic dynamics of cyclobutanone after excitation into the n → 3s Rydberg S2 state. It stems from our contribution to the Special Topic of the Journal of Chemical Physics to test the predictive capability of computational chemistry against unseen experimental data. Decoherence-corrected fewest-switches surface hopping was used to simulate nonadiabatic dynamics with full and approximated nonadiabatic couplings. Several simulation sets were computed with different electronic structure methods, including a multiconfigurational wavefunction [multiconfigurational self-consistent field (MCSCF)] specially built to describe dissociative channels, multireference semiempirical approach, time-dependent density functional theory, algebraic diagrammatic construction, and coupled cluster. MCSCF dynamics predicts a slow deactivation of the S2 state (10ps), followed by an ultrafast population transfer from S1 to S0 (<100fs). CO elimination (C3 channel) dominates over C2H4 formation (C2 channel). These findings radically differ from the other methods, which predicted S2 lifetimes 10-250 times shorter and C2 channel predominance. These results suggest that routine electronic structure methods may hold low predictive power for the outcome of nonadiabatic dynamics.

The photochemistry of Rydberg-excited cyclobutanone: Photoinduced processes and ground state dynamics.

Owing to ring strain, cyclic ketones exhibit complex excited state dynamics with multiple competing photochemical channels active on the ultrafast timescale. While the excited state dynamics of cyclobutanone after π* ← n excitation into the lowest-energy excited singlet (S1) state has been extensively studied, the dynamics following 3s ← n excitation into the higher-lying singlet Rydberg (S2) state are less well understood. Herein, we employ fully quantum multiconfigurational time-dependent Hartree (MCTDH) simulations using a model Hamiltonian as well as "on-the-fly" trajectory-based surface-hopping dynamics (TSHD) simulations to study the relaxation dynamics of cyclobutanone following 3s ← n excitation and to predict the ultrafast electron diffraction scattering signature of these relaxation dynamics. Our MCTDH and TSHD simulations indicate that relaxation from the initially-populated singlet Rydberg (S2) state occurs on the timescale of a few hundreds of femtoseconds to a picosecond, consistent with the symmetry-forbidden nature of the state-to-state transition involved. There is no obvious involvement of excited triplet states within the timeframe of our simulations (<2ps). After non-radiative relaxation to the electronic ground state (S0), vibrationally hot cyclobutanone has sufficient internal energy to form multiple fragmented products including C2H4 + CH2CO (C2; 20%) and C3H6 + CO (C3; 2.5%). We discuss the limitations of our MCTDH and TSHD simulations, how these may influence the excited state dynamics we observe, and-ultimately-the predictive power of the simulated experimental observable.

A simple one-electron expression for electron rotational factors.

Within the context of fewest-switch surface hopping (FSSH) dynamics, one often wishes to remove the angular component of the derivative coupling between states J and K. In a previous set of papers, Shu et al. [J. Phys. Chem. Lett. 11, 1135-1140 (2020)] posited one approach for such a removal based on direct projection, while we isolated a second approach by constructing and differentiating a rotationally invariant basis. Unfortunately, neither approach was able to demonstrate a one-electron operatorÔ whose matrix element JÔK was the angular component of the derivative coupling. Here, we show that a one-electron operator can, in fact, be constructed efficiently in a semi-local fashion. The present results yield physical insight into designing new surface hopping algorithms and are of immediate use for FSSH calculations.

Reconciling TD-DFT and CASPT2 electronic structure methods for describing the photophysics of DNA.

Time-dependent density functional theory (TD-DFT) and multiconfigurational second-order perturbation theory (CASPT2) are two of the most widely used methods to investigate photoinduced dynamics in DNA-based systems. These methods sometimes give diverse dynamics in physiological environments usually modeled by quantum mechanics/molecular mechanics (QM/MM) protocol. In this work, we demonstrate for the uridine test case that the underlying topology of the potential energy surfaces of electronic states involved in photoinduced relaxation is similar in both electronic structure methods. This is verified by analyzing surface-hopping dynamics performed at the QM/MM level on aqueous solvated uridine at TD-DFT and CASPT2 levels. By constraining the dynamics to remain on state we observe similar fluctuations in energy and relaxation lifetimes in surface-hopping dynamics in both TD-DFT and experimentally validated CASPT2 methods. This finding calls for a systematic comparison of the ES potential energy surfaces of DNA and RNA nucleosides at the single- and multi-reference levels of theory. The anomalous long excited state lifetime at the TD-DFT level is explained by trapping due to the tendency of TD-DFT in QM/MM schemes with electrostatic embedding to underestimate the energy of the state leading to a wrong energetic order. A study of the FC energetics suggests that improving the description of the surrounding environment through polarizable embedding or by the expansion of QM layer with hydrogen-bonded waters helps restore the correct state order at TD-DFT level. Thus by combining TDDFT with an accurate modeling of the environment, TD-DFT is positioned as the standout protocol to model photoinduced dynamics in DNA-based aggregates and multimers.

Recommendations for Velocity Adjustment in Surface Hopping.

This study investigates velocity adjustment directions after hopping in surface hopping dynamics. Using fulvene and a protonated Schiff base (PSB4) as case studies, we investigate the population decay and reaction yields of different sets of dynamics with the velocity adjusted in either the nonadiabatic coupling, gradient difference, or momentum directions. For the latter, in addition to the conventional algorithm, we investigated the performance of a reduced kinetic energy reservoir approach recently proposed. Our evaluation also considered velocity adjustment in the directions of approximate nonadiabatic coupling vectors. While results for fulvene are susceptible to the adjustment approach, PSB4 is not. We correlated this dependence to the topography near the conical intersections. When nonadiabatic coupling vectors are unavailable, the gradient difference direction is the best adjustment option. If the gradient difference is also unavailable, a semiempirical vector direction or the momentum direction with a reduced kinetic energy reservoir becomes an excellent option to prevent an artificial excess of back hoppings. The precise velocity adjustment direction is less crucial for describing the nonadiabatic dynamics than the kinetic energy reservoir's size.

Which Electronic Structure Method to Choose in Trajectory Surface Hopping Dynamics Simulations? Azomethane as a Case Study.

Nonadiabatic dynamics simulations have become a standard approach to explore photochemical reactions. Such simulations require underlying potential energy surfaces and couplings between them, calculated at a chosen level of theory, yet this aspect is rarely assessed. Here, in combination with the popular trajectory surface hopping dynamics method, we use a high-accuracy XMS-CASPT2 electronic structure level as a benchmark for assessing the performances of various post-Hartree-Fock methods (namely, CIS, ADC(2), CC2, and CASSCF) and exchange-correlation functionals (PBE, PBE0, and CAM-B3LYP) in a TD-DFT/TDA context, using the isomerization around a double bond as test case. Different relaxation pathways are identified, and the ability of the different methods to reproduce their relative importance and time scale is discussed. The results show that multireference electronic structure methods should be preferred, when studying nonadiabatic decay between excited and ground states. If not affordable, TD-DFT with TDA and hybrid functionals and ADC(2) are efficient alternatives but overestimate the nonradiative decay yield and thus may miss deexcitation pathways.

Photoluminescence mechanisms of BF2-formazanate dye sensitizers: a theoretical study.

Photodynamic therapy (PDT) is an alternative, minimally invasive treatment for human diseases such as cancer. PDT uses a photosensitizer to transfer photon energy directly to cellular 3O2 to generate 1O2 (Type II), the toxicity of which leads to cancer cell death. In this work, the photoluminescence mechanisms of a BF2-formazanate dye sensitizer (BF2-FORM) and its iodinated derivative (BF2-FORM-D) were studied using complementary theoretical approaches; the photoluminescence pathways in the S1 and T1 states were studied using density functional theory (DFT) and time-dependent (TD)-DFT methods, the kinetic and thermodynamic properties of the pathways using the transition state theory (TST), and the time evolution and dynamics of key processes using non-adiabatic microcanonical molecular dynamics simulations with surface-hopping dynamics (NVE-MDSH). Evaluation of the potential energy surfaces (PESs) in terms of the rotations of the phenyl rings suggested a pathway for the S1 → S0 transition for the perpendicular structure, whereas two pathways were anticipated for the T1 → S0 transition, namely, [T1 → S0]1 occurring immediately after the S1/T1 intersystem crossing (ISC) and [T1 → S0]2 occurring after the S1/T1 ISC and T1 equilibrium structure relaxation, with the T1 → S0 energy gap being comparable to the energy required for 3O2 → 1O2. The PESs also showed that because of the heavy-atom effect, BF2-FORM-D possessed a significantly smaller S1/T1 energy gap than BF2-FORM. The TST results revealed that at room temperature, BF2-FORM-D was thermodynamically more favorable than the parent molecule. Analysis of the NVE-MDSH results suggested that the librational motions of the phenyl rings play an important role in the internal conversion (IC) and ISC, and the S1/T1 ISC and T1 → S0 transitions could be enhanced by varying the irradiation wavelength and controlling the temperature. These findings can be used as guidelines to improve and/or design photosensitizers for PDT.

SpaiNN: equivariant message passing for excited-state nonadiabatic molecular dynamics.

Excited-state molecular dynamics simulations are crucial for understanding processes like photosynthesis, vision, and radiation damage. However, the computational complexity of quantum chemical calculations restricts their scope. Machine learning offers a solution by delivering high-accuracy properties at lower computational costs. We present SpaiNN, an open-source Python software for ML-driven surface hopping nonadiabatic molecular dynamics simulations. SpaiNN combines the invariant and equivariant neural network architectures of SchNetPack with SHARC for surface hopping dynamics. Its modular design allows users to implement and adapt modules easily. We compare rotationally-invariant and equivariant representations in fitting potential energy surfaces of multiple electronic states and properties arising from the interaction of two electronic states. Simulations of the methyleneimmonium cation and various alkenes demonstrate the superior performance of equivariant SpaiNN models, improving accuracy, generalization, and efficiency in both training and inference.

Photoinduced isomerization mechanism of isatin N&lt;sup&gt;2&lt;/sup&gt;-diphenylhydrazones molecular switch

Hydrazone molecular switches have significant application value in supramolecular chemistry. A new type of hydrazone molecular switch, named isatin N&lt;sup&gt;2&lt;/sup&gt;-diphenylhydrazone, has been synthesized. Owing to its cis-trans isomerization characteristics under visible light excitation, ease of synthesizing of derivatives, and sensitivity to external stimuli, it has important application value in the field of biochemistry. Because of its forward and backward visible light excitation characteristics, it is considered a class of compound that is very suitable for molecular switches, and it has a wide application value in fields such as biotechnology. In addition, the derivatives compound exhibits strong interactions with negative ions, which enhances its function as a molecular switch, making it a four-state molecular switch that can be achieved by a single molecule. However, the photo-induced isomerization mechanism of these new molecular switches is not yet clear, and whether there are novel phenomena in the isomerization process is also unknown. In this work, a semi empirical OM2/MRCI based trajectory surface hopping dynamics method is adopted to systematically study a photo induced isomerization mechanism based on the E-Z isomerization process of the isatin N&lt;sup&gt;2&lt;/sup&gt;-diphenylhydrazones molecular switch. Optimization configuration and the average lifetime of the first excited S&lt;sub&gt;1&lt;/sub&gt; state are obtained by using the semi-empirical OM2/MRCI method of molecular switch. It is found that the average lifetime of the S&lt;sub&gt;1&lt;/sub&gt; excited state of the E-configuration molecular switch is about 107 fs, and the quantum yield of E-Z isomerization of the molecular switch is 16.01%. By calculating the photo induced isomerization process of the molecular switch, two different isomerization mechanisms of the molecular switch are identified. In addition to the traditional molecular switch isomerization mechanism revolving around the C=N bond, a new isomerization mechanism, i.e. the face-to-face twisting of the molecular switch rotor part is elucidated. By calculating the time-resolved fluorescence radiation spectrum, it is predicted that there may be a very fast fluorescence quenching phenomenon occurring in about 75 fs in the isomerization process, slightly faster than the S&lt;sub&gt;1&lt;/sub&gt; average decay events (107 fs). The information about wavelength-resolved attenuation at different times is also calculated, which reflects the ultrafast fluorescence quenching process accompanied by fluorescence red shift, ranging from 2.1 × 10&lt;sup&gt;4&lt;/sup&gt; cm&lt;sup&gt;–1&lt;/sup&gt; to 3.4 × 10&lt;sup&gt;4&lt;/sup&gt; cm&lt;sup&gt;–1&lt;/sup&gt;. By comparing the calculated fluorescence spectra with the average lifetime of excited states, the existence of “dark states” is proposed, and possible explanations for the existence of “dark states” are provided, and those “dark states” may be related to lower quantum yields. The research results can provide theoretical guidance for the design and application of new molecular switches. The ease of synthesis and sensitivity to external stimuli of its derivatives make those compounds extremely valuable in molecular switching and light measurement applications.

Ab initio molecular dynamics study of intersystem crossing dynamics for MH2 (M = Si, Ge, Sn, Pb) on spin-pure and spin-mixed potential energy surfaces.

Recently, surface-hopping ab initio molecular dynamics (SH-AIMD) simulations have come to be used to discuss the mechanisms and dynamics of excited-state chemical reactions, including internal conversion and intersystem crossing. In dynamics simulations involving intersystem crossing, there are two potential energy surfaces (PESs) governing the motion of nuclei: PES in a spin-pure state and PES in a spin-mixed state. The former gives wrong results for molecular systems with large spin-orbit coupling (SOC), while the latter requires a potential gradient that includes a change in SOC at each point, making the computational cost very high. In this study, we systematically investigate the extent to which the magnitude of SOC affects the results of the spin-pure state-based dynamics simulations for the hydride MH2 (M = Si, Ge, Sn, Pb) by performing SH-AIMD simulations based on spin-pure and spin-mixed states. It is clearly shown that spin-mixed state PESs are indispensable for the dynamics simulation of intersystem crossing in systems containing elements Sn and Pb from the fifth period onward. Furthermore, in addition to the widely used Tully's fewest switches (TFS) algorithm, the Zhu-Nakamura (ZN) global switching algorithm, which is computationally less expensive, is applied to SH for comparison. The results from TFS- and ZN-SH-AIMD methods are in qualitative agreement, suggesting that the less expensive ZN-SH-AIMD can be successfully utilized to investigate the dynamics of photochemical reactions based on quantum chemical calculations.

Ultrafast ESIPT dynamics of 4-methyl-2,6-diformylphenol

We investigate the gas-phase photodynamics associated with the S1 (nπ*), S2 (ππ*) and S3 (nπ*) states of 4-methyl-2,6-diformylphenol using trajectory-based surface hopping dynamics simulations. A rapid internal conversion would populate S1 and S3 upon initiating the trajectories on “bright” S2. This nonadiabatic population behavior would lead to intramolecular proton-transfer tautomerization via three states, where the predominant tautomerization happens via S2. An average timescale of ∼ 44 fs is estimated for the tautomerization. The triplet formation can occur via energetically feasible S1-T3 (at the Franck-Condon geometry of normal tautomer) and S1-T2 (at the S1 equilibrium minimum of proton-transferred tautomer) pathways. Intersystem crossing and the forbidden nature of S1 (nπ*) makes the molecule exhibit weak fluorescence.