Nonequilibrium H/D Isotope Effects from Trajectory-Based Nonadiabatic Dynamics

Computational Comparison of Chemical and Isotopic Approaches to Control the Photoisomerization Dynamics of Light-Driven Molecular Motors

Synthetic molecular motors driven by E/Z photoisomerization reactions are able to produce unidirectional rotary motion because of a structural asymmetry that makes one direction of rotation more probable than the other. In most such motors, this asymmetry is realized through the incorporation of a chemically asymmetric carbon atom. Here, we present molecular dynamics simulations based on multiconfigurational quantum chemistry to investigate whether the merits of this approach can be equaled by an alternative approach that instead exploits isotopic chirality. By first considering an N-methylpyrrolidine–cyclopentadiene motor design, it is shown that isotopically chiral variants of this design undergo faster photoisomerizations than a chemically chiral counterpart, while maintaining rotary photoisomerization quantum yields of similarly high magnitude. However, by subsequently considering a pyrrolinium–cyclopentene design, it is also found that the introduction of isotopic chirality does not provide any control of the directionality of the photoinduced rotations within this framework. Taken together, the results highlight both the potential usefulness of isotopic rather than chemical chirality for the design of light-driven molecular motors, and the need for further studies to establish the exact structural circumstances under which this asymmetry is best exploited.

Unidirectional Rotary Motion in Isotopically Chiral Molecular Motors: A Computational Analysis

Molecular dynamics simulations are performed to explore if isotopic chirality can induce unidirectional rotary motion in molecular motors operated through double-bond photoisomerizations. Using a high-quantum yield motor featuring a chemically asymmetric carbon atom as reference, it is found that isotopically chiral counterparts of this motor sustain such motion almost equally well. Overall, the study reveals a previously unexplored role for isotopic chirality in the design of rotary molecular motors.

Nonadiabatic Dynamics Simulations on Early-Time Photochemistry of Spirobenzopyran

Photoinduced ring-opening, decay, and isomerization of spirobenzopyran have been explored by the OM2/MRCI nonadiabatic dynamics simulations based on Tully's fewest-switches surface hopping scheme. The efficient S₁ to S₀ internal conversion as observed in experiments is attributed to the existence of two efficient excited-state decay pathways. The first one is related to the C-N dissociation, and the second one is done to the C-O dissociation. The C-O dissociation pathway is dominant, and more than 90% trajectories decay to the S₀ state via the C-O bond-fission related S₁/S₀ conical intersections. Near these regions in the S₀ state, trajectories can either return to spirobenzopyran or proceed to various intermediates including merocyanine via a series of bond rotations. Our nonadiabatic dynamics simulations also demonstrate that the hydrogen-out-of-plane (HOOP) motion is important for efficient and ultrafast excited-state deactivation. On the other hand, we have also found that the replacement of methyl groups by hydrogen atoms in spirobenzopyran can artificially introduce different intramolecular hydrogen transfers leading to hydrogen-transferred intermediates. This finding is important for the community and demonstrates that such a kind of structural truncation, sometimes, could be problematic, leading to incorrect photodynamics. Our present work provides valuable insights into the photodynamics of spirobenzopyran, which could be helpful for the design of spiropyran-based photochromic materials.

GFN2-xTB-An Accurate and Broadly Parametrized Self-Consistent Tight-Binding Quantum Chemical Method with Multipole Electrostatics and Density-Dependent Dispersion Contributions

An extended semiempirical tight-binding model is presented, which is primarily designed for the fast calculation of structures and noncovalent interaction energies for molecular systems with roughly 1000 atoms. The essential novelty in this so-called GFN2-xTB method is the inclusion of anisotropic second order density fluctuation effects via short-range damped interactions of cumulative atomic multipole moments. Without noticeable increase in the computational demands, this results in a less empirical and overall more physically sound method, which does not require any classical halogen or hydrogen bonding corrections and which relies solely on global and element-specific parameters (available up to radon, Z = 86). Moreover, the atomic partial charge dependent D4 London dispersion model is incorporated self-consistently, which can be naturally obtained in a tight-binding picture from second order density fluctuations. Fully analytical and numerically precise gradients (nuclear forces) are implemented. The accuracy of the method is benchmarked for a wide variety of systems and compared with other semiempirical methods. Along with excellent performance for the "target" properties, we also find lower errors for "off-target" properties such as barrier heights and molecular dipole moments. High computational efficiency along with the improved physics compared to its precursor GFN-xTB makes this method well-suited to explore the conformational space of molecular systems. Significant improvements are furthermore observed for various benchmark sets, which are prototypical for biomolecular systems in aqueous solution.

Ultrafast action chemistry in slow motion: atomistic description of the excitation and fluorescence processes in an archetypal fluorescent protein

We report quantum mechanical/molecular mechanical non-adiabatic molecular dynamics simulations on the electronically excited state of green fluorescent protein mutant S65T/H148D. We examine the driving force of the ultrafast (τ < 50 fs) excited-state proton transfer unleashed by absorption in the A band at 415 nm and propose an atomistic description of the two dynamical regimes experimentally observed [Stoner Ma et al., J. Am. Chem. Soc., 2008, 130, 1227]. These regimes are explained in terms of two sets of successive dynamical events: first the proton transfers quickly from the chromophore to the acceptor Asp148. Thereafter, on a slower time scale, there are geometrical changes in the cavity of the chromophore that involve the distance between the chromophore and Asp148, the planarity of the excited-state chromophore, and the distance between the chromophore and Tyr145. We find two different non-radiative relaxation channels that are operative for structures in the reactant region and that can explain the mismatch between the decay of the emission of A* and the rise of the emission of I*, as well as the temperature dependence of the non-radiative decay rate.

Electronic structure calculations and nonadiabatic dynamics simulations of excited-state relaxation of Pigment Yellow 101

Pigment Yellow 101 (PY101) is widely used as a typical pigment due to its excellent excited-state properties. However, the origin of its photostability is still elusive. In this work, we have systematically investigated the photodynamics of PY101 by performing combined electronic structure calculations and trajectory-based nonadiabatic dynamics simulations. On the basis of the results, we have found that upon photoexcitation to the S₁ state, PY101 undergoes an essentially barrierless excited-state intramolecular single proton transfer generating an S₁ keto species. In the keto region, there is an energetically accessible S₁/S₀ conical intersection that funnels the system to the S₀ state quickly. In the S₀ state, the keto species either goes back to its trans-enol species through a ground-state reverse hydrogen transfer or arrives at the cis-keto region. In addition, we have found an additional excited-state decay channel for the S₁ enol species, which is directly linked to an S₁/S₀ conical intersection located in the enol region. This mechanism has also been confirmed by our dynamics simulations, in which about 54% of the trajectories decay to the S₀ state via the enol S₁/S₀ conical intersection; while the remaining ones employ the keto S₁/S₀ conical intersection. The gained mechanistic information helps us understand the photostability of the PY101 chromophore and its variants with the same molecular scaffold.

Tuning ESIPT fluorophores into dual emitters

Using first-principle approaches, we show how ESIPT can be controlled by fine-tuning of substituents, hence leading to new dual emitters.

Dyes undergoing excited-state intramolecular proton transfer (ESIPT) are known to present large Stokes shifts as a result of the important geometrical reorganisation following photon absorption. When the ESIPT process is not quantitative, one can obtain dual emitters characterised by two distinct fluorescence bands, observed due to emissions from both the canonical and ESIPT isomers. However, dual emission generally requires to maintain a very specific balance, as the relative excited-state free energies of the two tautomers have to be within a narrow window to observe the phenomenon. Consequently, simple chemical intuition is insufficient to optimise dual emission. In the present contribution, we investigate, with the help of quantum-mechanical tools and more precisely, time-dependent density functional theory (TD-DFT) and algebraic diagrammatic construction (ADC), a wide panel of possible ESIPT/dual emitters with various substituents. The selected protocol is first shown to be very robust on a series of structures with known experimental behaviour, and next is applied to novel derivatives with various substituents located at different positions. This work encompasses the largest chemical library of potential ESIPT compounds studied to date. We pinpoint the most promising combinations for building dual emitters, highlight unexpected combination effects and rationalise the impact of the different auxochromes.

Excited-state intramolecular proton transfer to carbon atoms: nonadiabatic surface-hopping dynamics simulations

Excited-state intramolecular proton transfer (ESIPT) between two highly electronegative atoms, for example, oxygen and nitrogen, has been intensely studied experimentally and computationally, whereas there has been much less theoretical work on ESIPT to other atoms such as carbon. We have employed CASSCF, MS-CASPT2, RI-ADC(2), OM2/MRCI, DFT, and TDDFT methods to study the mechanistic photochemistry of 2-phenylphenol, for which such an ESIPT has been observed experimentally. According to static electronic structure calculations, irradiation of 2-phenylphenol populates the bright S1 state, which has a rather flat potential in the Franck-Condon region (with a shallow enol minimum at the CASSCF level) and may undergo an essentially barrierless ESIPT to the more stable S1 keto species. There are two S1/S0 conical intersections that mediate relaxation to the ground state, one in the enol region and one in the keto region, with the latter one substantially lower in energy. After S1 → S0 internal conversion, the transient keto species can return back to the S0 enol structure via reverse ground-state hydrogen transfer in a facile tautomerization. This mechanistic scenario is verified by OM2/MRCI-based fewest-switches surface-hopping simulations that provide detailed dynamic information. In these trajectories, ESIPT is complete within 118 fs; the corresponding S1 excited-state lifetime is computed to be 373 fs in vacuum. Most of the trajectories decay to the ground state via the S1/S0 conical intersection in the keto region (67%), and the remaining ones via the enol region (33%). The combination of static electronic structure computations and nonadiabatic dynamics simulations is expected to be generally useful for understanding the mechanistic photophysics and photochemistry of molecules with intramolecular hydrogen bonds.

Excited-State Proton-Transfer-Induced Trapping Enhances the Fluorescence Emission of a Locked GFP Chromophore

The chemical locking of the central single bond in core chromophores of green fluorescent proteins (GFPs) influences their excited-state behavior in a distinct manner. Experimentally, it significantly enhances the fluorescence quantum yield of GFP chromophores with an ortho-hydroxyl group, while it has almost no effect on the photophysics of GFP chromophores with a para-hydroxyl group. To unravel the underlying physical reasons for this different behavior, we report static electronic structure calculations and nonadiabatic dynamics simulations on excited-state intramolecular proton transfer, cis–trans isomerization, and excited-state deactivation in a locked ortho-substituted GFP model chromophore (o-LHBI). On the basis of our previous and present results, we find that the S₁ keto species is responsible for the fluorescence emission of the unlocked o-HBI and the locked o-LHBI species. Chemical locking does not change the parts of the S₁ and S₀ potential energy surfaces relevant to enol–keto tautomerization; hence, in both chromophores, there is an ultrafast excited-state intramolecular proton transfer that takes only 35 fs on average. However, the locking effectively hinders the S₁ keto species from approaching the keto S₁/S₀ conical intersections so that most of trajectories are trapped in the S₁ keto region for the entire 2 ps simulation time. Therefore, the fluorescence quantum yield of o-LHBI is enhanced compared with that of unlocked o-HBI, in which the S₁ excited-state decay is efficient and ultrafast. In the case of the para-substituted GFP model chromophores p-HBI and p-LHBI, chemical locking hardly affects their efficient excited-state deactivation via cis–trans isomerization; thus, the fluorescence quantum yields in these chromophores remain very low. The insights gained from the present work may help to guide the design of new GFP chromophores with improved fluorescence emission and brightness.

Photophysics of Auramine-O: electronic structure calculations and nonadiabatic dynamics simulations

Sequential vs. concerted S₁ relaxation pathways.

Excited-state proton transfer in 4-2′-hydroxyphneylpyridine: full-dimensional surface-hopping dynamics simulations

We have employed combined electronic structure calculations and “on-the-fly” fewest switches surface-hopping dynamics simulations to study the S₁ excited-state intramolecular proton transfer (ESIPT) and decay dynamics of 4-(2′-hydroxyphenyl)pyridine.

Generalized trajectory surface-hopping method for internal conversion and intersystem crossing

Trajectory-based fewest-switches surface-hopping (FSSH) dynamics simulations have become a popular and reliable theoretical tool to simulate nonadiabatic photophysical and photochemical processes. Most available FSSH methods model internal conversion. We present a generalized trajectory surface-hopping (GTSH) method for simulating both internal conversion and intersystem crossing processes on an equal footing. We consider hops between adiabatic eigenstates of the non-relativistic electronic Hamiltonian (pure spin states), which is appropriate for sufficiently small spin-orbit coupling. This choice allows us to make maximum use of existing electronic structure programs and to minimize the changes to available implementations of the traditional FSSH method. The GTSH method is formulated within the quantum mechanics (QM)/molecular mechanics framework, but can of course also be applied at the pure QM level. The algorithm implemented in the GTSH code is specified step by step. As an initial GTSH application, we report simulations of the nonadiabatic processes in the lowest four electronic states (S0, S1, T1, and T2) of acrolein both in vacuo and in acetonitrile solution, in which the acrolein molecule is treated at the ab initio complete-active-space self-consistent-field level. These dynamics simulations provide detailed mechanistic insight by identifying and characterizing two nonadiabatic routes to the lowest triplet state, namely, direct S1 → T1 hopping as major pathway and sequential S1 → T2 → T1 hopping as minor pathway, with the T2 state acting as a relay state. They illustrate the potential of the GTSH approach to explore photoinduced processes in complex systems, in which intersystem crossing plays an important role.

Assessment of approximate computational methods for conical intersections and branching plane vectors in organic molecules

Quantum-chemical computational methods are benchmarked for their ability to describe conical intersections in a series of organic molecules and models of biological chromophores. Reference results for the geometries, relative energies, and branching planes of conical intersections are obtained using ab initio multireference configuration interaction with single and double excitations (MRCISD). They are compared with the results from more approximate methods, namely, the state-interaction state-averaged restricted ensemble-referenced Kohn-Sham method, spin-flip time-dependent density functional theory, and a semiempirical MRCISD approach using an orthogonalization-corrected model. It is demonstrated that these approximate methods reproduce the ab initio reference data very well, with root-mean-square deviations in the optimized geometries of the order of 0.1 Å or less and with reasonable agreement in the computed relative energies. A detailed analysis of the branching plane vectors shows that all currently applied methods yield similar nuclear displacements for escaping the strong non-adiabatic coupling region near the conical intersections. Our comparisons support the use of the tested quantum-chemical methods for modeling the photochemistry of large organic and biological systems.

Photoswitching of salicylidene methylamine: a theoretical photodynamics study

Photoswitching of simple photochromic molecules attracts substantial attention because of its possible role in future photon-driven molecular electronics. Here we model the full photoswitching cycle of a minimal photochromic Schiff base–salicylidene methylamine (SMA). We perform semiempirical nonadiabatic on-the-fly photodynamics simulations at the OM2/MRCI level and thoroughly analyze the structural time evolution and switching efficiency of the system. We also identify and examine in detail the crucial steps in the SMA photochemistry ruled by excited-state intramolecular proton transfer. The results place the investigated model aromatic Schiff base among the promising candidates for novel photoswitching molecular materials. Our study also shows the potential of the semiempirical multireference photodynamics simulations as a tool for early stage molecular photodevice design.

Investigation of ESIPT in a panel of chromophores presenting N-H∙∙∙N intramolecular hydrogen bonds

Thermodynamic and kinetic aspects of excited state intramolecular proton transfer (ESIPT) are investigated in 11 chromophores harboring an intramolecular N-H∙∙∙N hydrogen bond [pyridyl pyrazole, pyridyl pyrrole, azaindole, pyridyl indole, pyrroloquinoline, and an analogue of the Blue Fluorescent Protein (BFP) chromophore] with the help of quantum mechanical calculations. For pyridyl pyrazoles, simulated spectra are used to help the interpretation of experimental ones and the effects of several substituents are investigated. Then it is shown that Time-Dependent Density Functional Theory fails to satisfactorily describe the energetic aspects of ESIPT for the BFP chromophore analogue. Equation-of-Motion Coupled Cluster theory is thus used to reach accurate insights for this challenging case.