Nonadiabatic dynamics: The SHARC approach

Reverse Intersystem Crossing Dynamics in Vibronically Modulated Inverted Singlet-Triplet Gap System: A Wigner Phase Space Study

We inspect the origin of the inverted singlet-triplet gap (INVEST) and slow change in the reverse intersystem crossing (rISC) rate with temperature, as recently observed. A Wigner phase space study reveals that, though INVEST is found at equilibrium geometry, variation in the exchange interaction and the doubles-excitation for other geometries in the harmonic region leads to non-INVEST behavior. This highlights the importance of nuclear degrees of freedom for the INVEST phenomenon, and in this case, geometric puckering of the studied molecule determines INVEST and the associated rISC dynamics.

Untangling the catalytic importance of the Se oxidation state in organoselenium-mediated oxygen-transfer reactions: the conversion of aniline to nitrobenzene

Seleninic acids and their precursors are well-known oxygen-transfer agents that can catalyze several oxidations with H2O2 as the final oxidant. Until very recently, the Se(iv) "peroxyseleninic" acid species has been considered the only plausible catalytic oxidant. Conversely, in 2020, the involvement of Se(vi) "peroxyselenonic" acid has been proposed for the selenium mediated epoxidation of alkenes. In this work, we theoretically probe different mechanisms of H2O2 activation and of Se(iv) to Se(vi) interconversion. In addition, we investigate through a combined theoretical (DFT) and experimental approach the mechanistic steps leading to the oxidation of aniline to nitrobenzene, when Se(iv) seleninic acid or Se(vi) selenonic acids are used as catalysts and H2O2 as the oxidant. This process encompasses three subsequent organoselenium mediated oxidations by H2O2. These results provide a mechanistic explanation of the advantages and disadvantages of both oxidation states (iv and vi) in the different stages of catalytic oxygen-transfer reactions: hydrogen peroxide activation and actual substrate oxidation. While the Se(vi) "peroxyselenonic" acid is found to be a better oxidant, the privileged role of "peroxyseleninic" acid as the main active species is assessed and its origin is identified in the lower catalyst-distortion that seleninic acid undergoes when activating H2O2. Conversely, the higher catalyst-distortion that characterizes the reaction of selenonic acid with H2O2 supports an inactivating role of Se(iv) to Se(vi) interconversion.

Implementation of Nonadiabatic Molecular Dynamics for Intersystem Crossing Based on a Time-Dependent Density-Functional Tight-Binding Method

Intersystem crossing (ISC) and internal conversion (IC) are types of nonadiabatic transitions that play important roles in a wide range of fields, including photochemistry, photophysics, and photobiology. The nonadiabatic molecular dynamics (NA-MD) method is a powerful tool for computational simulations of dynamic phenomena involving nonadiabatic transitions. In this study, we implemented the NA-MD method, which treats ISC and IC on an equal footing, where the electronic structure is treated at the level of the time-dependent (TD) density-functional tight-binding (DFTB) method, a low-cost semiempirical analog of TD density functional theory (DFT). In particular, the spin-orbit coupling calculation algorithm was implemented in the TD-DFTB framework, and the results showed trends similar to those obtained using TD-DFT. In addition, the NA-MD method successfully reproduced ultrafast ISC of 2-nitronaphthalene.

Ultrafast Photoinduced Dynamics in 1,3-Cyclohexadiene: A Comparison of Trajectory Surface Hopping Schemes†

Photoinduced nonadiabatic processes play a crucial role in a wide range of disciplines, from fundamental steps in biology to modern applications in advanced materials science. A theoretical understanding of these processes is highly desirable, and trajectory surface hopping (TSH) has proven to be a well-suited framework for a wide range of systems. In this work, we present a comprehensive comparison between two TSH algorithms, the conventional Tully's fewest switches surface hopping (FSSH) scheme and the Landau-Zener surface hopping (LZSH), to study the photoinduced ring-opening of 1,3-cyclohexadiene (CHD) to 1,3,5-hexatriene at the spin-flip time-dependent density functional theory (SF-TDDFT) level of theory. Additionally, we compare our results with a literature study at the extended multistate complete active space second-order perturbation theory method (XMS-CASPT2) level of theory. Our results show that the average population and lifetimes estimated with LZSH using SF-TDDFT are closer to the literature (using multireference methods) than those estimated with FSSH using SF-TDDFT. The latter speaks in favor of applying LZSH in combination with the SF-TDDFT method to study larger and more complex systems such as molecular photoswitches where the CHD molecule acts as a backbone. In addition, we present an implementation of Tully's FSSH algorithm as an extension to the PySurf software package.

Stepwise Excited-State Intramolecular Double Proton Transfer in 1,8-Dihydroxynaphthalene-2,7-dicarbaldehyde

We theoretically investigate the salient features of stepwise excited-state intramolecular double proton transfer (ESIDPT) in 1,8-dihydroxynaphthalene-2,7-dicarbaldehyde (DHDA). Surface trajectory simulations using the TD-B3LYP/6-31G(d) level of theory reveal that the proton transfer primarily happens via S1, wherein about ∼42% of trajectories (40 out of 95) show the single proton transfer alone and another ∼32% (30 out of 95) show double proton transfer. The average time scale for the single proton transfer originating from those ∼42% trajectories is ∼147 fs. In the case of double proton transfer, the average time for the first step, i.e., single proton transfer, is ∼54 fs, and the subsequent step, i.e., double proton transfer, completes in ∼151 fs. All three tautomers, i.e., normal, single, and double proton-transferred tautomers, possess a stable minimum in their first singlet excited state. This state has a ππ* character in the former two tautomers, resulting in a dual fluorescence emission phenomenon upon photoexcitation of DHDA.

Perspective on Theoretical and Experimental Advances in Atmospheric Photochemistry

Research that explores the chemistry of Earth's atmosphere is central to the current understanding of global challenges such as climate change, stratospheric ozone depletion, and poor air quality in urban areas. This research is a synergistic combination of three established domains: earth observation, for example, using satellites, and in situ field measurements; computer modeling of the atmosphere and its chemistry; and laboratory measurements of the properties and reactivity of gas-phase molecules and aerosol particles. The complexity of the interconnected chemical and photochemical reactions which determine the composition of the atmosphere challenges the capacity of laboratory studies to provide the spectroscopic, photochemical, and kinetic data required for computer models. Here, we consider whether predictions from computational chemistry using modern electronic structure theory and nonadiabatic dynamics simulations are becoming sufficiently accurate to supplement quantitative laboratory data for wavelength-dependent absorption cross-sections, photochemical quantum yields, and reaction rate coefficients. Drawing on presentations and discussions from the CECAM workshop on Theoretical and Experimental Advances in Atmospheric Photochemistry held in March 2024, we describe key concepts in the theory of photochemistry, survey the state-of-the-art in computational photochemistry methods, and compare their capabilities with modern experimental laboratory techniques. From such considerations, we offer a perspective on the scope of computational (photo)chemistry methods based on rigorous electronic structure theory to become a fourth core domain of research in atmospheric chemistry.

Excitonic Configuration Interaction: Going Beyond the Frenkel Exciton Model

We present the excitonic configuration interaction (ECI) method─a fragment-based analogue of the CI method for electronic structure calculations of multichromophoric systems. It can also be viewed as a generalization of the exciton approach, with the following properties: (i) It constructs the effective Hamiltonian exclusively from monomer calculations. (ii) It employs the strong orthogonality assumption and is exact within McWeeny's group function theory, thus requiring only one-electron density matrices of the monomer states. (iii) It is agnostic of the monomer electronic structure method, allowing us to use/combine different methods. (iv) It includes an embedding point charge scheme (called excitonic Hartree-Fock, EHF) to improve the accuracy of the monomer states, but such that the effective full-system Hamiltonian is not explicitly dependent on the embedding. (v) It is systematically improvable, by expanding the set of monomer states and by including configurations where two or more monomers are excited (defining the ECIS, ECISD, etc., methods). The performance of ECI is assessed by computing the absorption spectrum of two exemplary multichromophoric systems, using CIS as the monomer electronic structure method. The accuracy of ECI significantly depends on the chosen embedding charges and the ECI expansion. The most accurate assessed combinations─ECIS or ECISD with EHF embedding─yield spectra that agree qualitatively and quantitatively with full-system direct calculations, with deviations of the excitation energies below 0.1 eV. We also show that ECISD based on CIS monomer calculations can predict states where two monomers are excited simultaneously (e.g., triplet-triplet double-local excitations) that are inaccessible in a full-system CIS calculation.

A size-consistent multi-state mapping approach to surface hopping

We develop a multi-state generalization of the recently proposed mapping approach to surface hopping (MASH) for the simulation of electronically nonadiabatic dynamics. This new approach extends the original MASH method to be able to treat systems with more than two electronic states. It differs from previous approaches in that it is size consistent and rigorously recovers the original two-state MASH in the appropriate limits. We demonstrate the accuracy of the method by applying it to a series of model systems for which exact benchmark results are available, and we find that the method is well suited to the simulation of photochemical relaxation processes.

Nonadiabatic Dynamics with Exact Factorization: Implementation and Assessment

In this work, we report our implementation of several independent-trajectory mixed-quantum-classical (ITMQC) nonadiabatic dynamics methods based on exact factorization (XF) in the Libra package for nonadiabatic and excited-state dynamics. Namely, the exact factorization surface hopping (SHXF), mixed quantum-classical dynamics (MQCXF), and mean-field (MFXF) are introduced. Performance of these methods is compared to that of several traditional surface hopping schemes, such as the fewest-switches surface hopping (FSSH), branching-corrected surface hopping (BCSH), and the simplified decay of mixing (SDM), as well as conventional Ehrenfest (mean-field, MF) method. Based on a comprehensive set of 1D model Hamiltonians, we find the ranking SHXF ≈ MQCXF > BCSH > SDM > FSSH ≫ MF, with the BCSH sometimes outperforming the XF methods in terms of describing coherences. Although the MFXF method can yield reasonable populations and coherences for some cases, it does not conserve the total energy and is therefore not recommended. We also find that the branching correction for auxiliary trajectories is important for the XF methods to yield accurate populations and coherences. However, the branching correction can worsen the quality of the energy conservation in the MQCXF. Finally, we find that using the time-dependent Gaussian width approximation used in the XF methods for computing decoherence correction can improve the quality of energy conservation in the MQCXF dynamics. The parameter-free scheme of Subotnik for computing the Gaussian widths is found to deliver the best performance in situations where such widths are not known a priori.

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.

Resolving Photoinduced Femtosecond Three-Dimensional Solute-Solvent Dynamics through Surface Hopping Simulations

Photoinduced dynamics in solution is governed by mutual solute-solvent interactions, which give rise to phenomena like solvatochromism, the Stokes shift, dual fluorescence, or charge transfer. Understanding these phenomena requires simulating the solute's photoinduced dynamics and simultaneously resolving the three-dimensional solvent distribution dynamics. If using trajectory surface hopping (TSH) to this aim, thousands of trajectories are required to adequately sample the time-dependent three-dimensional solvent distribution functions, and thus resolve the solvent dynamics with sub-Ångstrom and femtosecond accuracy and sufficiently low noise levels. Unfortunately, simulating thousands of trajectories with TSH in the framework of hybrid quantum mechanical/molecular mechanical (QM/MM) can be prohibitively expensive when employing ab initio electronic structure methods. To tackle this challenge, we recently introduced a computationally efficient approach that combines efficient linear vibronic coupling models with molecular mechanics (LVC/MM) via electrostatic embedding [Polonius et al., JCTC 2023, 19, 7171-7186]. This method provides solvent-embedded, nonadiabatically coupled potential energy surfaces while scaling similarly to MM force fields. Here, we employ TSH with LVC/MM to unravel the photoinduced dynamics of two small thiocarbonyl compounds solvated in water. We describe how to estimate the number of trajectories required to produce nearly noise-free three-dimensional solvent distribution functions and present an analysis based on approximately 10,000 trajectories propagated for 3 ps. In the electronic ground state, both molecules exhibit in-plane hydrogen bonds to the sulfur atom. Shortly after excitation, these bonds are broken and reform perpendicular to the molecular plane on timescales that differ by an order of magnitude due to steric effects. We also show that the solvent relaxation dynamics is coupled to the electronic dynamics, including intersystem crossing. These findings are relevant to advance the understanding of the coupled solute-solvent dynamics of solvated photoexcited molecules, e.g., biologically relevant thio-nucleobases.

Generalized Semiclassical Ehrenfest Method: A Route to Wave Function-Free Photochemistry and Nonadiabatic Dynamics with Only Potential Energies and Gradients

We reconsider recent methods by which direct dynamics calculations of electronically nonadiabatic processes can be carried out while requiring only adiabatic potential energies and their gradients. We show that these methods can be understood in terms of a new generalization of the well-known semiclassical Ehrenfest method. This is convenient because it eliminates the need to evaluate electronic wave functions and their matrix elements along the mixed quantum-classical trajectories. The new approximations and procedures enabling this advance are the curvature-driven approximation to the time-derivative coupling, the generalized semiclassical Ehrenfest method, and a new gradient correction scheme called the time-derivative matrix (TDM) scheme. When spin-orbit coupling is present, one can carry out dynamics calculations in the fully adiabatic basis using potential energies and gradients calculated without spin-orbit coupling plus the spin-orbit coupling matrix elements. Even when spin-orbit coupling is neglected, the method is useful because it allows calculations by electronic structure methods for which nonadiabatic coupling vectors are unavailable. In order to place the new considerations in context, the article starts out with a review of background material on trajectory surface hopping, the semiclassical Ehrenfest scheme, and methods for incorporating decoherence. We consider both internal conversion and intersystem crossing. We also review several examples from our group of successful applications of the curvature-driven approximation.

Quantum Quality with Classical Cost: Ab Initio Nonadiabatic Dynamics Simulations Using the Mapping Approach to Surface Hopping

Nonadiabatic dynamics methods are an essential tool for investigating photochemical processes. In the context of employing first-principles electronic structure techniques, such simulations can be carried out in a practical manner using semiclassical trajectory-based methods or wave packet approaches. While all approaches applicable to first-principles simulations are necessarily approximate, it is commonly thought that wave packet approaches offer inherent advantages over their semiclassical counterparts in terms of accuracy and that this trait simply comes at a higher computational cost. Here we demonstrate that the mapping approach to surface hopping (MASH), a recently introduced trajectory-based nonadiabatic dynamics method, can be efficiently applied in tandem with ab initio electronic structure. Our results even suggest that MASH may provide more accurate results than on-the-fly wave packet techniques, all at a much lower computational cost.

Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone

Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α-CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.

Modeling the temperature dependence of the fluorescence properties of Indole in aqueous solution

In a recent paper, we proposed a scheme to describe the relaxation mechanism of the excited Indole in aqueous solution, involving the fluctuations among the diabatic electronic states 1Lb, 1La and 1πσ∗. Such a theoretical and computational model reproduced accurately the available experimental data at room temperature. Following these results, in the present work, we model the complex temperature dependence of the fluorescence properties of Indole in aqueous solution, with results further validating the proposed relaxation scheme. This scheme is able to explain the temperature effects on the fluorescence behavior indicating the water fluctuations as the main cause of (i) the stabilization of the dark state (1πσ∗) and (ii) the increase in temperature of the kinetics of the irreversible transition towards such a state.

Nuclear Quantum Effects on Nonadiabatic Dynamics of a Green Fluorescent Protein Chromophore Analogue: Ring-Polymer Surface-Hopping Simulation

Herein, we have used the "on-the-fly" ring-polymer surface-hopping simulation method with the centroid approximation (RPSH-CA), in combination with the multireference OM2/MRCI electronic structure calculations to study the photoinduced dynamics of a green fluorescent protein (GFP) chromophore analogue in the gas phase, i.e., o-HBI, at 50, 100, and 300 K with 1, 5, 10, and 15 beads (3600 1 ps trajectories). The electronic structure calculations identified five new minimum-energy conical intersection (MECI) structures, which, together with the previous one, play crucial roles in the excited-state decay dynamics of o-HBI. It is also found that the excited-state intramolecular proton transfer (ESIPT) occurs in an ultrafast manner and is completed within 20 fs in all the simulation conditions because there is no barrier associated with this ESIPT process in the S1 state. However, the other excited-state dynamical results are strongly related to the number of beads. At 50 and 100 K, the nuclear quantum effects (NQEs) are very important; therefore, the excited-state dynamical results change significantly with the bead number. For example, the S1 decay time deduced from time-dependent state populations becomes longer as the bead number increases. Nevertheless, an essentially convergent trend is observed when the bead number is close to 10. In contrast, at 300 K, the NQEs become weaker and the above dynamical results converge very quickly even with 1 bead. Most importantly, the NQEs seriously affect the excited-state decay mechanism of o-HBI. At 50 and 100 K, most trajectories decay to the S0 state via perpendicular keto MECIs, whereas, at 300 K, only twisted keto MECIs are responsible for the excited-state decay. The present work not only comprehensively explores the temperature-dependent photoinduced dynamics of o-HBI, but also demonstrates the importance and necessity of NQEs in nonadiabatic dynamics simulations, especially at relatively low temperatures.

Machine-Learned Kohn-Sham Hamiltonian Mapping for Nonadiabatic Molecular Dynamics

In this work, we report a simple, efficient, and scalable machine-learning (ML) approach for mapping non-self-consistent Kohn-Sham Hamiltonians constructed with one kind of density functional to the nearly self-consistent Hamiltonians constructed with another kind of density functional. This approach is designed as a fast surrogate Hamiltonian calculator for use in long nonadiabatic dynamics simulations of large atomistic systems. In this approach, the input and output features are Hamiltonian matrices computed from different levels of theory. We demonstrate that the developed ML-based Hamiltonian mapping method (1) speeds up the calculations by several orders of magnitude, (2) is conceptually simpler than alternative ML approaches, (3) is applicable to different systems and sizes and can be used for mapping Hamiltonians constructed with arbitrary density functionals, (4) requires a modest training data, learns fast, and generates molecular orbitals and their energies with the accuracy nearly matching that of conventional calculations, and (5) when applied to nonadiabatic dynamics simulation of excitation energy relaxation in large systems yields the corresponding time scales within the margin of error of the conventional calculations. Using this approach, we explore the excitation energy relaxation in C60 fullerene and Si75H64 quantum dot structures and derive qualitative and quantitative insights into dynamics in these systems.

Formulation of transition dipole gradients for non-adiabatic dynamics with polaritonic states

A general formulation of the strong coupling between photons confined in a cavity and molecular electronic states is developed for the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham method. The light-matter interaction is included in the Jaynes-Cummings model, which requires the derivation and implementation of the analytical derivatives of the transition dipole moments between the molecular electronic states. The developed formalism is tested in the simulations of the nonadiabatic dynamics in the polaritonic states resulting from the strong coupling between the cavity photon mode and the ground and excited states of the penta-2,4-dieniminium cation, also known as PSB3. Comparison with the field-free simulations of the excited-state decay dynamics in PSB3 reveals that the light-matter coupling can considerably alter the decay dynamics by increasing the excited state lifetime and hindering photochemically induced torsion about the C=C double bonds of PSB3. The necessity of obtaining analytical transition dipole gradients for the accurate propagation of the dynamics is underlined.

A Portrait of the Chromophore as a Young System-Quantum-Derived Force Field Unraveling Solvent Reorganization upon Optical Excitation of Cyclocurcumin Derivatives

The study of fast non-equilibrium solvent relaxation in organic chromophores is still challenging for molecular modeling and simulation approaches, and is often overlooked, even in the case of non-adiabatic dynamics simulations. Yet, especially in the case of photoswitches, the interaction with the environment can strongly modulate the photophysical outcomes. To unravel such a delicate interplay, in the present contribution we resorted to a mixed quantum-classical approach, based on quantum mechanically derived force fields. The main task is to rationalize the solvent reorganization pathways in chromophores derived from cyclocurcumin, which are suitable for light-activated chemotherapy to destabilize cellular lipid membranes. The accurate and reliable decryption delivered by the quantum-derived force fields points to important differences in the solvent's reorganization, in terms of both structure and time scale evolution.

Bifurcation of Excited-State Population Leads to Anti-Kasha Luminescence in a Disulfide-Decorated Organometallic Rhenium Photosensitizer

We report a rhenium diimine photosensitizer equipped with a peripheral disulfide unit on one of the bipyridine ligands, [Re(CO)3(bpy)(S-Sbpy4,4)]+ (1+, bpy = 2,2'-bipyridine, S-Sbpy4,4 = [1,2]dithiino[3,4-c:6,5-c']dipyridine), showing anti-Kasha luminescence. Steady-state and ultrafast time-resolved spectroscopies complemented by nonadiabatic dynamics simulations are used to disclose its excited-state dynamics. The calculations show that after intersystem crossing the complex evolves to two different triplet minima: a (S-Sbpy4,4)-ligand-centered excited state (3LC) lying at lower energy and a metal-to-(bpy)-ligand charge transfer (3MLCT) state at higher energy, with relative yields of 90% and 10%, respectively. The 3LC state involves local excitation of the disulfide group into the antibonding σ* orbital, leading to significant elongation of the S-S bond. Intriguingly, it is the higher-lying 3MLCT state, which is assigned to display luminescence with a lifetime of 270 ns: a signature of anti-Kasha behavior. This assignment is consistent with an energy barrier ≥ 0.6 eV or negligible electronic coupling, preventing reaction toward the 3LC state after the population is trapped in the 3MLCT state. This study represents a striking example on how elusive excited-state dynamics of transition-metal photosensitizers can be deciphered by synergistic experiments and state-of-the-art calculations. Disulfide functionalization lays the foundation of a new design strategy toward harnessing excess energy in a system for possible bimolecular electron or energy transfer reactivity.

Prediction challenge: First principles simulation of the ultrafast electron diffraction spectrum of cyclobutanone

Computer simulation has long been an essential partner of ultrafast experiments, allowing the assignment of microscopic mechanistic detail to low-dimensional spectroscopic data. However, the ability of theory to make a priori predictions of ultrafast experimental results is relatively untested. Herein, as a part of a community challenge, we attempt to predict the signal of an upcoming ultrafast photochemical experiment using state-of-the-art theory in the context of preexisting experimental data. Specifically, we employ ab initio Ehrenfest with collapse to a block mixed quantum-classical simulations to describe the real-time evolution of the electrons and nuclei of cyclobutanone following excitation to the 3s Rydberg state. The gas-phase ultrafast electron diffraction (GUED) signal is simulated for direct comparison to an upcoming experiment at the Stanford Linear Accelerator Laboratory. Following initial ring-opening, dissociation via two distinct channels is observed: the C3 dissociation channel, producing cyclopropane and CO, and the C2 channel, producing CH2CO and C2H4. Direct calculations of the GUED signal indicate how the ring-opened intermediate, the C2 products, and the C3 products can be discriminated in the GUED signal. We also report an a priori analysis of anticipated errors in our predictions: without knowledge of the experimental result, which features of the spectrum do we feel confident we have predicted correctly, and which might we have wrong?

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.

Suppression of Dexter transfer by covalent encapsulation for efficient matrix-free narrowband deep blue hyperfluorescent OLEDs

Hyperfluorescence shows great promise for the next generation of commercially feasible blue organic light-emitting diodes, for which eliminating the Dexter transfer to terminal emitter triplet states is key to efficiency and stability. Current devices rely on high-gap matrices to prevent Dexter transfer, which unfortunately leads to overly complex devices from a fabrication standpoint. Here we introduce a molecular design where ultranarrowband blue emitters are covalently encapsulated by insulating alkylene straps. Organic light-emitting diodes with simple emissive layers consisting of pristine thermally activated delayed fluorescence hosts doped with encapsulated terminal emitters exhibit negligible external quantum efficiency drops compared with non-doped devices, enabling a maximum external quantum efficiency of 21.5%. To explain the high efficiency in the absence of high-gap matrices, we turn to transient absorption spectroscopy. It is directly observed that Dexter transfer from a pristine thermally activated delayed fluorescence sensitizer host can be substantially reduced by an encapsulated terminal emitter, opening the door to highly efficient 'matrix-free' blue hyperfluorescence.

Ultrafast electronic relaxation pathways of the molecular photoswitch quadricyclane

The light-induced ultrafast switching between molecular isomers norbornadiene and quadricyclane can reversibly store and release a substantial amount of chemical energy. Prior work observed signatures of ultrafast molecular dynamics in both isomers upon ultraviolet excitation but could not follow the electronic relaxation all the way back to the ground state experimentally. Here we study the electronic relaxation of quadricyclane after exciting in the ultraviolet (201 nanometres) using time-resolved gas-phase extreme ultraviolet photoelectron spectroscopy combined with non-adiabatic molecular dynamics simulations. We identify two competing pathways by which electronically excited quadricyclane molecules relax to the electronic ground state. The fast pathway (<100 femtoseconds) is distinguished by effective coupling to valence electronic states, while the slow pathway involves initial motions across Rydberg states and takes several hundred femtoseconds. Both pathways facilitate interconversion between the two isomers, albeit on different timescales, and we predict that the branching ratio of norbornadiene/quadricyclane products immediately after returning to the electronic ground state is approximately 3:2.

Universal Measure for the Impact of Adiabaticity on Quantum Transitions

Adiabaticity is crucial for our understanding of complex quantum dynamics and thus for advancing fundamental physics and technology, but its impact cannot yet be quantified in complex but common cases where dynamics is only partially adiabatic, several eigenstates are simultaneously populated and transitions between noneigenstates are of key interest. We construct a universally applicable measure that can quantify the adiabaticity of quantum transitions in an arbitrary basis. Our measure distinguishes transitions that occur due to the adiabatic change of populated system eigenstates from transitions that occur due to beating between several eigenstates and can handle nonadiabatic events. While all quantum dynamics fall within the scope of the measure, we demonstrate its usage and utility through two important material science problems-energy and charge transfer-where adiabaticity could be effected by nuclear motion and its quantification will aid not only in unraveling mechanisms but also in system design, for example, of light harvesting systems.

Impact of large A-site cations on electron-vibrational interactions in 2D halide perovskites: Ab initio quantum dynamics

Using ab initio nonadiabatic molecular dynamics, we study the effect of large A-site cations on nonradiative electron-hole recombination in two-dimensional Ruddlesden-Popper perovskites HA2APb2I7, HA = n-hexylammonium, A = methylammonium (MA), or guanidinium (GA). The steric hindrance created by large GA cations distorts and stiffens the inorganic Pb-I lattice, reduces thermal structural fluctuations, and maintains the delocalization of electrons and holes at ambient and elevated temperatures. The delocalized charges interact more strongly in the GA system than in the MA system, and the charge recombination is accelerated. In contrast, replacement of only some MA cations with GA enhances disorder and increases charge lifetime, as seen in three-dimensional perovskites. This study highlights the key influence of structural fluctuations and disorder on the properties of charge carriers in metal halide perovskites, providing guidance for tuning materials' optoelectronic performance.

Dissociation Time, Quantum Yield, and Dynamic Reaction Pathways in the Thermolysis of trans-3,4-Dimethyl-1,2-dioxetane

The thermolysis of trans-3,4-dimethyl-1,2-dioxetane is studied by trajectory surface hopping. The significant difference between long and short dissociation times is rationalized by frustrated dissociations and the time spent in triplet states. If the C-C bond breaks through an excited state channel, then the trajectory passes over a ridge of the potential energy surface of that state. The calculated triplet quantum yields match the experimental results. The dissociation half-times and quantum yields follow the same ascending order as per the product states, justifying the conjecture that the longer dissociation time leads to a higher quantum yield, proposed in the context of the methylation effect. The populations of the molecular Coulomb Hamiltonian and diagonal states reach equilibrium, but the triplet populations with different Sz components fluctuate indefinitely. Certain initial velocities, leading the trajectories to given product states, can be identified as the most characteristic features for sorting trajectories according to their product states.

Monitoring the Evolution of Relative Product Populations at Early Times during a Photochemical Reaction

Identifying multiple rival reaction products and transient species formed during ultrafast photochemical reactions and determining their time-evolving relative populations are key steps toward understanding and predicting photochemical outcomes. Yet, most contemporary ultrafast studies struggle with clearly identifying and quantifying competing molecular structures/species among the emerging reaction products. Here, we show that mega-electronvolt ultrafast electron diffraction in combination with ab initio molecular dynamics calculations offer a powerful route to determining time-resolved populations of the various isomeric products formed after UV (266 nm) excitation of the five-membered heterocyclic molecule 2(5H)-thiophenone. This strategy provides experimental validation of the predicted high (∼50%) yield of an episulfide isomer containing a strained three-membered ring within ∼1 ps of photoexcitation and highlights the rapidity of interconversion between the rival highly vibrationally excited photoproducts in their ground electronic state.

Simulating Non-Adiabatic Dynamics of Photoexcited Phenyl Azide: Investigating Electronic and Structural Relaxation en Route to the Formation of Phenyl Nitrene

Excited state molecular dynamics simulations of the photoexcited phenyl azide have been performed. The semi-classical surface hopping approximation has enabled an unconstrained analysis of the electronic and nuclear degrees of freedom which contribute to the molecular dissociation of phenyl azide into phenyl nitrene and molecular nitrogen. The significance of the second singlet excited state in leading the photodissociation has been established through electronic structure calculations, based on multi-configurational schemes, and state population dynamics. The investigations on the structural dynamics have revealed the N-N bond separation to be accompanied by synchronous changes in the azide N-N-N bond angle. The 100 fs simulation results in a nitrene fragment that is electronically excited in the singlet manifold.

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.

Linear and angular momentum conservation in surface hopping methods

We demonstrate that, for systems with spin-orbit coupling and an odd number of electrons, the standard fewest switches surface hopping algorithm does not conserve the total linear or angular momentum. This lack of conservation arises not so much from the hopping direction (which is easily adjusted) but more generally from propagating adiabatic dynamics along surfaces that are not time reversible. We show that one solution to this problem is to run along eigenvalues of phase-space electronic Hamiltonians H(R, P) (i.e., electronic Hamiltonians that depend on both nuclear position and momentum) with an electronic-nuclear coupling Γ · P [see Eq. (25)], and we delineate the conditions that must be satisfied by the operator Γ. The present results should be extremely useful as far as developing new semiclassical approaches that can treat systems where the nuclear, electronic orbital, and electronic spin degrees of freedom altogether are all coupled together, hopefully including systems displaying the chiral-induced spin selectivity effect.

Spin-Vibronic Dynamics in Open-Shell Systems beyond the Spin Hamiltonian Formalism

Vibronic coupling has a dramatic influence over a large number of molecular processes, ranging from photochemistry to spin relaxation and electronic transport. The simulation of vibronic coupling with multireference wave function methods has been largely applied to organic compounds, and only early efforts are available for open-shell systems such as transition metal and lanthanide complexes. In this work, we derive a numerical strategy to differentiate the molecular electronic Hamiltonian in the context of multireference ab initio methods and inclusive of spin-orbit coupling effects. We then provide a formulation of open quantum system dynamics able to predict the time evolution of the electron density matrix under the influence of a Markovian phonon bath up to fourth-order perturbation theory. We apply our method to Co(II) and Dy(III) molecular complexes exhibiting long spin relaxation times and successfully validate our strategy against the use of an effective spin Hamiltonian. Our study sheds light on the nature of vibronic coupling, the importance of electronic excited states in spin relaxation, and the need for high-level computational chemistry to quantify it.

A Proof-of-Principle Design for Through-Space Transmission of Unidirectional Rotary Motion by Molecular Photogears

The construction of molecular photogears that can achieve through-space transmission of the unidirectional double-bond rotary motion of light-driven molecular motors onto a remote single-bond axis is a formidable challenge in the field of artificial molecular machines. Here, we present a proof-of-principle design of such photogears that is based on the possibility of using stereogenic substituents to control both the relative stabilities of two helical forms of the photogear and the double-bond photoisomerization reaction that connects them. The potential of the design was verified by quantum-chemical modeling through which photogearing was found to be a favorable process compared to free-standing single-bond rotation ("slippage"). Overall, our study unveils a surprisingly simple approach to realizing unidirectional photogearing.

Energy-Conserving and Thermally Corrected Neglect of Back-Reaction Approximation Method for Nonadiabatic Molecular Dynamics

In this work, the energy-conserving and thermally corrected neglect of the back-reaction approximation approach for nonadiabatic molecular dynamics in extended atomistic systems is developed. The new approach introduces three key corrections to the original method: (1) it enforces the total energy conservation, (2) it introduces an explicit coupling of the system to its environment, and (3) it introduces a renormalization of nonadiabatic couplings to account for a difference between the instantaneous nuclear kinetic energy and the kinetic energy of guiding trajectories. In the new approach, an auxiliary kinetic energy variable is introduced as an independent dynamical variable. The new approach produces nonzero equilibrium populations, whereas the original neglect of the back-reaction approximation method does not. It yields population relaxation time scales that are favorably comparable to the reference values, and it introduces an explicit and controllable way of dissipating energy into a bath without an assumption of the bath being at equilibrium.

Controlling Charge Carrier Dynamics in Porphyrin Nanorings by Optically Active Templates

Understanding the dynamics of photogenerated charge carriers is essential for enhancing the performance of solar and optoelectronic devices. Using atomistic quantum dynamics simulations, we demonstrate that a short π-conjugated optically active template can be used to control hot carrier relaxation, charge carrier separation, and carrier recombination in light-harvesting porphyrin nanorings. Relaxation of hot holes is slowed by 60% with an optically active template compared to that with an analogous optically inactive template. Both systems exhibit subpicosecond electron transfer from the photoactive core to the templates. Notably, charge recombination is suppressed 6-fold by the optically active template. The atomistic time-domain simulations rationalize these effects by the extent of electron and hole localization, modification of the density of states, participation of distinct vibrational motions, and changes in quantum coherence. Extension of the hot carrier lifetime and reduction of charge carrier recombination, without hampering charge separation, demonstrate a strategy for enhancing efficiencies of energy materials with optically active templates.

Vibrational Coherences of the Photoinduced Mixed-Valent Creutz-Taube Ion Revealed by Excited State Dynamics

A recent study of photoinduced mixed-valency in the one-electron reduced form (μ-pz)[RuII(NH3)5]24+ of the Creutz-Taube ion used transient absorption spectroscopy with vis-NIR broadband detection to uncover a mixed-valent excited state with a typical intervalence charge transfer band and a nanosecond lifetime [Pieslinger et al. Angew. Chem., Int. Ed. 2022, 61, e202211747]. Herein, we use excited state dynamics simulations with implicit solvation to elucidate the electronic and vibrational evolution in the first 10 ps after the optical excitation. A manifold of excited states with weak interaction between the metal centers is populated already at time zero due to the breakdown of the Condon approximation and dominates the population of electronic states at short time scales (<0.5 ps). A long-lived vibrational wave packet mostly confined to oscillations of the metal center-bridge distances is observed. The oscillations are traced to the electronic structure properties of states with weak metal-metal coupling. The long-lived mixed-valent excited state of the Creutz-Taube ion analogue is formed vibrationally cold and has a more compact geometry. While experimentally, intersystem crossing and vibrational relaxation were deduced to be completed within 1 ps, our analysis indicates that both processes might persist at longer times.

What Controls the Quality of Photodynamical Simulations? Electronic Structure Versus Nonadiabatic Algorithm

The field of nonadiabatic dynamics has matured over the last decade with a range of algorithms and electronic structure methods available at the moment. While the community currently focuses more on developing and benchmarking new nonadiabatic dynamics algorithms, the underlying electronic structure controls the outcome of nonadiabatic simulations. Yet, the electronic-structure sensitivity analysis is typically neglected. In this work, we present a sensitivity analysis of the nonadiabatic dynamics of cyclopropanone to electronic structure methods and nonadiabatic dynamics algorithms. In particular, we compare wave function-based CASSCF, FOMO-CASCI, MS- and XMS-CASPT2, density-functional REKS, and semiempirical MRCI-OM3 electronic structure methods with the Landau-Zener surface hopping, fewest switches surface hopping, and ab initio multiple spawning with informed stochastic selection algorithms. The results clearly demonstrate that the electronic structure choice significantly influences the accuracy of nonadiabatic dynamics for cyclopropanone even when the potential energy surfaces exhibit qualitative and quantitative similarities. Thus, selecting the electronic structure solely on the basis of the mapping of potential energy surfaces can be misleading. Conversely, we observe no discernible differences in the performance of the nonadiabatic dynamics algorithms across the various methods. Based on the above results, we discuss the present-day practice in computational photodynamics.

Excited-State Dynamics Simulations of a Light-Driven Molecular Motor in Solution

Molecular motors, where light can be transformed into motion, are promising in the design of nanomechanical devices. For applications, however, finding relationships between molecular motion and the environment is important. Here, we report the study of excited-state dynamics of an overcrowded alkene in solution using a hybrid quantum mechanics/molecular mechanics (QM/MM) approach combined with excited-state molecular dynamics simulations. Using QM/MM surface-hopping trajectories, we calculated time-resolved emission and transient absorption spectra. These show the rise of a short-lived Franck-Condon state, followed by the formation of a dark state in the first 150 fs before the molecular motor relaxes to the ground state in about 1 ps. From the analysis of radial distribution functions, we infer that the orientation of the solvent with respect to the molecular motor in the electronic excited state is similar to that in the ground state during the photoisomerization.

Twisting Aromaticity and Photoinduced Dynamics in Hexapole Helicenes

Curved aromatic molecules are attractive electronic materials, where an additional internal strain uniquely modifies their structure, aromaticity, dynamics, and optical properties. Helicenes are examples of such twisted conjugated systems. Herein, we analyze the photoinduced dynamics in different stereoisomers of a hexapole helicene by using nonadiabatic excited-state molecular dynamics simulations. We explore how changes in symmetry and structural distortion modulate the intramolecular energy redistribution. We find that distinct helical assembly leads to different rigid distorted structures that in turn impact the nonradiative energy relaxation and ultimately formation of the self-trapped exciton. Subsequently, the value of the twisting angles relative to the central triphenylene core structure controls the global molecular aromaticity and electronic localization during the internal conversion process. Our work sheds light on how the future synthesis of novel curved aromatic compounds can be directed to attain specific desired electronic properties through the modulation of their twisted aromaticity.

Substituent-Induced Hyperconjugation: Origin of the Structural Effects on the Efficiency of Photochemical Ring Opening

Photochemical ring-opening reactions are among the most extensively employed chemical reactions in the field of chemistry. Owing to their significance, molecular-level studies of these reactions have been widely conducted. One of the major considerations in investigating the ring-opening dynamics of complex molecules on the molecular scale is the differences in dynamics between different conformers because the number of conformers arising from a specific substrate rapidly increases with the complexity of the substrate. However, to date, studies dealing with this problem have been limited to specific individual cases. That is, a rule applicable to arbitrary conformers to estimate and explain the effects of the molecular structure, such as substituents and conformations, on photochemical ring opening has not been established. Herein, we propose the concept of substituent-induced electron density leakage via hyperconjugation as a candidate for this general rule. Based on our hypothesis, we present an indicator that can predict the efficiency of the photochemical ring-opening reactions of various conformers. The relative error between the ring-opening efficiency as obtained from the indicator and that obtained from the nonadiabatic simulations was less than 25% in 56 of the 66 conformers arising from 1,3-cyclohexadiene and 12 distinct analogues. This approach offers the possibility of accurately and quickly predicting the photochemical ring-opening efficiency of arbitrary molecules in arbitrary conformations.

Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry

The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.

A Perspective on Sustainable Computational Chemistry Software Development and Integration

The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.

LVC/MM: A Hybrid Linear Vibronic Coupling/Molecular Mechanics Model with Distributed Multipole-Based Electrostatic Embedding for Highly Efficient Surface Hopping Dynamics in Solution

We present a theoretical framework for a hybrid linear vibronic coupling model electrostatically embedded into a molecular mechanics environment, termed the linear vibronic coupling/molecular mechanics (LVC/MM) method, for the surface hopping including arbitrary coupling (SHARC) molecular dynamics package. Electrostatic embedding is realized through the computation of interactions between environment point charges and distributed multipole expansions (DMEs, up to quadrupoles) that represent each electronic state and transition densities in the diabatic basis. The DME parameters are obtained through a restrained electrostatic potential (RESP) fit, which we extended to yield higher-order multipoles. We also implemented in SHARC a scheme for achieving roto-translational invariance of LVC models as well as a general quantum mechanics/molecular mechanics (QM/MM) interface, an OpenMM interface, and restraining potentials for simulating liquid droplets. Using thioformaldehyde in water as a test case, we demonstrate that LVC/MM can accurately reproduce the solvation structure and energetics of rigid solutes, with errors on the order of 1-2 kcal/mol compared to a BP86/MM reference. The implementation in SHARC is shown to be very efficient, enabling the simulation of trajectories on the nanosecond time scale in a matter of days.

The impact of the chalcogen-substitution element and initial spectroscopic state on excited-state relaxation pathways in nucleobase photosensitizers: a combination of static and dynamic studies

The substitution of oxygen with chalcogen in carbonyl group(s) of canonical nucleobases gives an impressive triplet generation, enabling their promising applications in medicine and other emerging techniques. The excited-state relaxation S2(ππ*) → S1(nπ*) → T1(ππ*) has been considered the preferred path for triplet generation in these nucleobase derivatives. Here, we demonstrate enhanced quantum efficiency of direct intersystem crossing from S2 to triplet manifold upon substitution with heavier chalcogen elements. The excited-state relaxation dynamics of sulfur/selenium substituted guanines in a vacuum is investigated using a combination of static quantum chemical calculations and on-the-fly excited-state molecular dynamics simulations. We find that in sulfur-substitution the S2 state predominantly decays to the S1 state, while upon selenium-substitution the S2 state deactivation leads to simultaneous population of the S1 and T2,3 states in the same time scale and multi-state quasi-degeneracy region S2/S1/T2,3. Interestingly, the ultrafast deactivation of the spectroscopic S3 state of both studied molecules to the S1 state occurs through a successive S3 → S2 → S1 path involving a multi-state quasi-degeneracy S3/S2/S1. The populated S1 and T2 states will cross the lowest triplet state, and the S1 → T intersystem crossing happens in a multi-state quasi-degeneracy region S1/T2,3/T1 and is accelerated by selenium-substitution. The present study reveals the influence of both the chalcogen substitution element and initial spectroscopic state on the excited-state relaxation mechanism of nucleobase photosensitizers and also highlights the important role of multi-state quasi-degeneracy in mediating the complex relaxation process. These theoretical results provide additional insights into the intrinsic photophysics of nucleobase-based photosensitizers and are helpful for designing novel photo-sensitizers for real applications.

Nature of Hops, Coordinates, and Detailed Balance for Nonadiabatic Simulations in the Condensed Phase

Photoinduced processes play a crucial role in a multitude of important molecular phenomena. Accurately modeling these processes in an environment other than a vacuum requires a detailed description of the electronic states involved as well as how energy flows are coupled to the surroundings. Nonadiabatic effects must also be included in order to describe the exchange of energy between electronic and nuclear degrees of freedom correctly. In this work, we revisit the ring-opening reaction 1,3-cylohexadiene (CHD) in a solvent environment. Using our newly developed Interface for Non-Adiabatic Quantum mechanics/molecular mechanics in Solvent (INAQS) we trace the evolution of the reaction via hybrid quantum mechanics/molecular mechanics (QM/MM) surface hopping with a focus on the solvent's participation in the nonadiabatic relaxation process and the long-time approach to equilibrium. We explicitly include the MM solvent contribution to the nonadiabatic coupling vector─enabling an accurate approach to equilibrium at long times─and find that in highly multidimensional systems gradients can have little or nothing to do with the nonadiabatic couplings.

Nonadiabatic Coupling in Trajectory Surface Hopping: Accurate Time Derivative Couplings by the Curvature-Driven Approximation

Trajectory surface hopping (TSH) is a widely used mixed quantum-classical dynamics method that is used to simulate molecular dynamics with multiple electronic states. In TSH, time-derivative coupling is employed to propagate the electronic coefficients and in that way to determine when the electronic state on which the nuclear trajectory is propagated switches. In this work, we discuss nonadiabatic TSH dynamics algorithms employing the curvature-driven approximation and overlap-based time derivative couplings, and we report test calculations on six photochemical reactions where we compare the results to one another and to calculations employing analytic nonadiabatic coupling vectors. We correct previous published results thanks to a bug found in the software. We also provide additional, more detailed studies of the time-derivative couplings. Our results show good agreement between curvature-driven algorithms and overlap-based algorithms.

HORTENSIA, a program package for the simulation of nonadiabatic autoionization dynamics in molecules

We present a program package for the simulation of ultrafast vibration-induced autoionization dynamics in molecular anions in the manifold of the adiabatic anionic states and the discretized ionization continuum. This program, called HORTENSIA (Hopping Real-time Trajectories for Electron-ejection by Nonadiabatic Self-Ionization in Anions), is based on the nonadiabatic surface-hopping methodology, wherein nuclei are propagated as an ensemble along classical trajectories in the quantum-mechanical potential created by the electronic density of the molecular system. The electronic Schrödinger equation is numerically integrated along the trajectory, providing the time evolution of electronic state coefficients, from which switching probabilities into discrete electronic states are determined. In the case of a discretized continuum state, this hopping event is interpreted as the ejection on an electron. The derived diabatic and nonadiabatic couplings in the time-dependent electronic Schrödinger equation are calculated from anionic and neutral wavefunctions obtained from quantum-chemical calculations with commercially available program packages interfaced with our program. Based on this methodology, we demonstrate the simulation of autoionization electron kinetic energy spectra that are both time- and angle-resolved. In addition, the program yields data that can be interpreted easily with respect to geometric characteristics, such as bonding distances and angles, which facilitate the detection of molecular configurations important for the autoionization process. Furthermore, several useful extensions are included, namely, tools for the generation of initial conditions and input files as well as for the evaluation of output files, all of this both through console commands and a graphical user interface.

Ultrafast Photocontrolled Rotation in a Molecular Motor Investigated by Machine Learning-Based Nonadiabatic Dynamics Simulations

The thermal helix inversion (THI) of the overcrowded alkene-based molecular motors determines the speed of the unidirectional rotation due to the high reaction barrier in the ground state, in comparison with the ultrafast photoreaction process. Recently, a phosphine-based motor has achieved all-photochemical rotation experimentally, promising to be controlled without a thermal step. However, the mechanism of this photochemical reaction has not yet been fully revealed. The comprehensive computational studies on photoisomerization still resort to nonadiabatic molecular dynamics (NAMD) simulations based on electronic structure calculations, which remains a high computational cost for large systems such as molecular motors. Machine learning (ML) has become an accelerating tool in NAMD simulations recently, where excited-state potential energy surfaces (PESs) are constructed analytically with high accuracy, providing an efficient approach for simulations in photochemistry. Herein the reaction pathway is explored by a spin-flip time-dependent density functional theory (SF-TDDFT) approach in combination with ML-based NAMD simulations. According to our computational simulations, we notice that one of the key factors of fulfilling all-photochemical rotation in the phosphine-based motor is that the excitation energies of four isomers are similar. Additionally, a shortcut photoinduced transformation between unstable isomers replaces the THI step, which shares the conical intersection (CI) with photoisomerization. In this study, we provide a practical approach to speed up the NAMD simulations in photochemical reactions for a large system that could be extended to other complex systems.

Efficient Modeling of Quantum Dynamics of Charge Carriers in Materials Using Short Nonequilibrium Molecular Dynamics

Nonadiabatic molecular dynamics provides essential insights into excited-state processes, but it is computationally intense and simplifications are needed. The classical path approximation provides critical savings. Still, long heating and equilibration steps are required. We demonstrate that practical results can be obtained with short, partially equilibrated ab initio trajectories. Once the system's structure is adequate and essential fluctuations are sampled, the nonadiabatic Hamiltonian can be constructed. Local structures require only 1-2 ps trajectories, as demonstrated with point defects in metal halide perovskites. Short trajectories represent anharmonic motions common in defective structures, an essential improvement over the harmonic approximation around the optimized geometry. Glassy systems, such as grain boundaries, require different simulation protocols, e.g., involving machine learning force fields. 10-fold shorter trajectories generate 10-20% time scale errors, which are acceptable, given experimental uncertainties and other approximations. The practical NAMD protocol enables fast screening of excited-state dynamics for rapid exploration of new materials.

Significance of Nonadiabatic Effects on Efficient Triplet Generation in Lumazines

The photophysics of lumazines leading to triplet formation and the effect of thionation are explored in the presence of near-degenerate electronic states. Wave packet simulations are performed on model potential energy surfaces to understand the nonadiabatic population transfer among close-lying excited states. Ultrafast population transfer among singlets opens up new intersystem crossing channels from the higher states. An increased spin-orbit coupling strength originating from thionation enhances intersystem crossing and populates the higher triplets first. The rapid internal conversion in the triplet manifold eventually brings the molecules to their respective low-lying long-lived triplet state.