Evidence for a vibrational phase-dependent isotope effect on the photochemistry of vision

Comparative quantum-classical dynamics of natural and synthetic molecular rotors show how vibrational synchronization modulates the photoisomerization quantum efficiency

We use quantum-classical trajectories to investigate the origin of the different photoisomerization quantum efficiency observed in the dim-light visual pigment Rhodopsin and in the light-driven biomimetic molecular rotor para-methoxy N-methyl indanylidene-pyrrolinium (MeO-NAIP) in methanol. Our results reveal that effective light-energy conversion requires, in general, an auxiliary molecular vibration (called promoter) that does not correspond to the rotary motion but synchronizes with it at specific times. They also reveal that Nature has designed Rhodopsin to exploit two mechanisms working in a vibrationally coherent regime. The first uses a wag promoter to ensure that ca. 75% of the absorbed photons lead to unidirectional rotations. The second mechanism ensures that the same process is fast enough to avoid directional randomization. It is found that MeO-NAIP in methanol is incapable of exploiting the above mechanisms resulting into a 50% quantum efficiency loss. However, when the solvent is removed, MeO-NAIP rotation is predicted to synchronize with a ring-inversion promoter leading to a 30% increase in quantum efficiency and, therefore, biomimetic behavior.

Picosecond quantum-classical dynamics reveals that the coexistence of light-induced microbial and animal chromophore rotary motion modulates the isomerization quantum yield of heliorhodopsin

Rhodopsins are light-responsive proteins forming two vast and evolutionary distinct superfamilies whose functions are invariably triggered by the photoisomerization of a single retinal chromophore. In 2018 a third widespread superfamily of rhodopsins called heliorhodopsins was discovered using functional metagenomics. Heliorhodopsins, with their markedly different structural features with respect to the animal and microbial superfamilies, offer an opportunity to study how evolution has manipulated the chromophore photoisomerization to achieve adaptation. One question is related to the mechanism of such a reaction and how it differs from that of animal and microbial rhodopsins. To address this question, we use hundreds of quantum-classical trajectories to simulate the spectroscopically documented picosecond light-induced dynamics of a heliorhodopsin from the archaea thermoplasmatales archaeon (TaHeR). We show that, consistently with the observations, the trajectories reveal two excited state decay channels. However, inconsistently with previous hypotheses, only one channel is associated with the -C13C14- rotation of microbial rhodopsins while the second channel is characterized by the -C11C12- rotation typical of animal rhodopsins. The fact that such -C11C12- rotation is aborted upon decay and ground state relaxation, explains why illumination of TaHeR only produces the 13-cis isomer with a low quantum efficiency. We argue that the documented lack of regioselectivity in double-bond excited state twisting motion is the result of an "adaptation" that could be completely lost via specific residue substitutions modulating the steric hindrance experienced along the isomerization motion.

Hydrogen Isotope Effect Endows a Breakthrough in Photoluminescent Covalent Organic Frameworks

Luminescent covalent organic frameworks (LCOFs) have emerged as indispensable candidates in various applications due to their greater tunable emitting properties and structural robustness compared to small molecule emitters. An unsolved issue in this area is developing highly luminescent LCOFs of which the nonradiative quenching pathways were suppressed as much as possible. Here, a robust aminal-linked COF (DD-COF) possessing perdeuterated light-emitting monomers was designed and synthesized. The solid-state photoluminescence quantum yield of the DD-COF reaches 81%, significantly outcompeting all state-of-the-art LCOFs reported so far. The exceptional luminescent efficiency is attributed to the inhibition of different pathways of nonradiative decay, especially from bond vibrations where only substitution by a heavier isotope with a lower zero-point vibration frequency works. Furthermore, the prepared deuterated COF not only boosts higher photostability under UV irradiation but also enables superior fluorescence sensing performance for iodine detection compared to nondeuterated COF.

Assessment of the Electron Correlation Treatment on the Quantum-Classical Dynamics of Retinal Protonated Schiff Base Models: XMS-CASPT2, RMS-CASPT2, and REKS Methods

We compare the performance of three different multiconfigurational wave function-based electronic structure methods and two implementations of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method. The study is characterized by three features: (i) it uses a small set of quantum-classical trajectories rather than potential energy surface mapping, (ii) it focuses, exclusively, on the photoisomerization of retinal protonated Schiff base models, and (iii) it probes the effect of both methyl substitution and the increase in length of the conjugate π-system. For each tested method, the corresponding analytical gradients are used to drive the quantum-classical (Tully's FSSH method) trajectory propagation, including the recent multistate XMS-CASPT2 and RMS-CASPT2 gradients. It is shown that while CASSCF, XMS-CASPT2, and RMS-CASPT2 yield consistent photoisomerization dynamics descriptions, REKS produces, in some of these systems, qualitatively different behavior that is attributed to a flatter and topographically different excited state potential energy surface. The origin of this behavior can be traced back to the effect of the employed density functional approximation. The above studies are further expanded by benchmarking, at the CASSCF and REKS levels, the electronic structure methods using a QM/MM model of the visual pigment rhodopsin.

Addressing high excitation conditions in time-resolved X-ray diffraction experiments and issues of biological relevance

One of the most important fundamental questions connecting chemistry to biology is how chemistry scales in complexity up to biological systems where there are innumerable possible pathways and competing processes. With the development of ultrabright electron and x-ray sources, it has been possible to literally light up atomic motions to directly observe the reduction in dimensionality in the barrier crossing region to a few key reaction modes. How do these chemical processes further couple to the surrounding protein or macromolecular assembly to drive biological functions? Optical methods to trigger photoactive biological processes are needed to probe this issue on the relevant timescales. However, the excitation conditions have been in the highly nonlinear regime, which questions the biological relevance of the observed structural dynamics.

Deciphering Methylation Effects on S2(ππ*) Internal Conversion in the Simplest Linear α,β-Unsaturated Carbonyl

Chemical substituents can influence photodynamics by altering the location of critical points and the topography of the potential energy surfaces (electronic effect) and by selectively modifying the inertia of specific nuclear modes (inertial effects). Using nonadiabatic dynamics simulations, we investigate the impact of methylation on S₂(ππ*) internal conversion in acrolein, the simplest linear α,β-unsaturated carbonyl. Consistent with time constants reported in a previous time-resolved photoelectron spectroscopy study, S₂ → S₁ deactivation occurs on an ultrafast time scale (∼50 fs). However, our simulations do not corroborate the sequential decay model used to fit the experiment. Instead, upon reaching the S₁ state, the wavepacket bifurcates: a portion undergoes ballistic S₁ → S₀ deactivation (∼90 fs) mediated by fast bond-length alternation motion, while the remaining decays on the picosecond time scale. Our analysis reveals that methyl substitution, generally assumed to mainly exert inertial influence, is also manifested in important electronic effects due to its weak electron-donating ability. While methylation at the β C atom gives rise to effects principally of an inertial nature, such as retarding the twisting motion of the terminal -CHCH₃ group and increasing its coupling with pyramidalization, methylation at the α or carbonyl C atom modifies the potential energy surfaces in a way that also contributes to altering the late S₁-decay behavior. Specifically, our results suggest that the observed slowing of the picosecond component upon α-methylation is a consequence of a tighter surface and reduced amplitude along the central pyramidalization, effectively restricting the access to the S₁/S₀-intersection seam. Our work offers new insight into the S₂(ππ*) internal conversion mechanisms in acrolein and its methylated derivatives and highlights site-selective methylation as a tuning knob to manipulate photochemical reactions.

Counterion at an atypical position: Investigating the mechanism of photoisomerization in jellyfish rhodopsin

The position of the counterion in animal rhodopsins plays a crucial role in maintaining visible light sensitivity and facilitating the photoisomerization of their retinal chromophore. The counterion displacement is thought to be closely related to the evolution of rhodopsins, with different positions found in invertebrates and vertebrates. Interestingly, box jellyfish rhodopsin (JelRh) acquired the counterion in transmembrane 2 (TM2) independently. This is a unique feature, as in most animal rhodopsins, the counterion is found in a different location. In this study, we used Fourier Transform Infrared spectroscopy to examine the structural changes that occur in the early photointermediate state of JelRh. We aimed to determine whether the photochemistry of JelRh is similar to that of other animal rhodopsins by comparing its spectra to those of vertebrate bovine rhodopsin (BovRh) and invertebrate squid rhodopsin (SquRh). We observed that the N-D stretching band of the retinal Schiff base was similar to that of BovRh, indicating the interaction between the Schiff base and the counterion is similar in both rhodopsins, despite their different counterion positions. Furthermore, we found that the chemical structure of the retinal in JelRh is similar to that in BovRh, including the changes in the hydrogen-out-of-plane band that indicates a retinal distortion. Overall, the protein conformational changes induced by the photoisomerization of JelRh yielded spectra that resemble an intermediate between BovRh and SquRh, suggesting a unique spectral property of JelRh, and making it the only animal rhodopsin with a counterion in TM2 and an ability to activate G_s protein.

New Gradient Correction Scheme for Electronically Nonadiabatic Dynamics Involving Multiple Spin States

It has been recommended that the best representation to use for trajectory surface hopping (TSH) calculations is the fully adiabatic basis in which the Hamiltonian is diagonal. Simulations of intersystem crossing processes with conventional TSH methods require an explicit computation of nonadiabatic coupling vectors (NACs) in the molecular-Coulomb-Hamiltonian (MCH) basis, also called the spin-orbit-free basis, in order to compute the gradient in the fully adiabatic basis (also called the diagonal representation). This explicit requirement destroys some of the advantages of the overlap-based algorithms and curvature-driven algorithms that can be used for the most efficient TSH calculations. Therefore, although these algorithms allow one to perform NAC-free simulations for internal conversion processes, one still requires NACs for intersystem crossing. Here, we show that how the NAC requirement is circumvented by a new computation scheme called the time-derivative-matrix scheme.

A Nitrogen Out-of-Plane (NOOP) Mechanism for Imine-Based Light-Driven Molecular Motors

Light-driven molecular motors have generated considerable interest due to their potential applications in material and biological systems. Recently, Greb and Lehn reported a new class of molecular motors, chiral N-alkyl imines, which undergo unidirectional rotation induced by light and heat. The mechanism of unidirectional motion in molecular motors containing a C═N group has been assumed to consist of photoinduced torsion about the double bond. In this work, we present a computational study of the photoisomerization dynamics of a chiral N-alkyl imine motor. We find that the location and energetics of minimal energy conical intersections (MECIs) alone are insufficient to understand the mechanism of the motor. Furthermore, a key part of the mechanism consists of out-of-plane distortions of the N atom (followed by isomerization about the double bond). Dynamic effects and out-of-plane distortions are critical to understand the observed (rather low) quantum yield for photoisomerization. Our results provide hints as to how the photoisomerization quantum yield might be increased, improving the efficiency of this class of molecular motors.

Ultrafast structural changes direct the first molecular events of vision

Vision is initiated by the rhodopsin family of light-sensitive G protein-coupled receptors (GPCRs)1. A photon is absorbed by the 11-cis retinal chromophore of rhodopsin, which isomerizes within 200 femtoseconds to the all-trans conformation2, thereby initiating the cellular signal transduction processes that ultimately lead to vision. However, the intramolecular mechanism by which the photoactivated retinal induces the activation events inside rhodopsin remains experimentally unclear. Here we use ultrafast time-resolved crystallography at room temperature3 to determine how an isomerized twisted all-trans retinal stores the photon energy that is required to initiate the protein conformational changes associated with the formation of the G protein-binding signalling state. The distorted retinal at a 1-ps time delay after photoactivation has pulled away from half of its numerous interactions with its binding pocket, and the excess of the photon energy is released through an anisotropic protein breathing motion in the direction of the extracellular space. Notably, the very early structural motions in the protein side chains of rhodopsin appear in regions that are involved in later stages of the conserved class A GPCR activation mechanism. Our study sheds light on the earliest stages of vision in vertebrates and points to fundamental aspects of the molecular mechanisms of agonist-mediated GPCR activation.

From a one-mode to a multi-mode understanding of conical intersection mediated ultrafast organic photochemical reactions

This review discusses how ultrafast organic photochemical reactions are controlled by conical intersections, highlighting that decay to the ground-state at multiple points of the intersection space results in their multi-mode character.

Direct structural observation of ultrafast photoisomerization dynamics in sinapate esters

Sinapate esters have been extensively studied for their potential application in 'nature-inspired' photoprotection. There is general consensus that the relaxation mechanism of sinapate esters following photoexcitation with ultraviolet radiation is mediated by geometric isomerization. This has been largely inferred through indirect studies involving transient electronic absorption spectroscopy in conjunction with steady-state spectroscopies. However, to-date, there is no direct experimental evidence tracking the formation of the photoisomer in real-time. Using transient vibrational absorption spectroscopy, we report on the direct structural changes that occur upon photoexcitation, resulting in the photoisomer formation. Our mechanistic analysis predicts that, from the photoprepared ππ* state, internal conversion takes place through a conical intersection (CI) near the geometry of the initial isomer. Our calculations suggest that different CI topographies at relevant points on the seam of intersection may influence the isomerization yield. Altogether, we provide compelling evidence suggesting that a sinapate ester's geometric isomerization can be a more complex dynamical process than originally thought.

On the fluorescence enhancement of arch neuronal optogenetic reporters

The lack of a theory capable of connecting the amino acid sequence of a light-absorbing protein with its fluorescence brightness is hampering the development of tools for understanding neuronal communications. Here we demonstrate that a theory can be established by constructing quantum chemical models of a set of Archaerhodopsin reporters in their electronically excited state. We found that the experimentally observed increase in fluorescence quantum yield is proportional to the computed decrease in energy difference between the fluorescent state and a nearby photoisomerization channel leading to an exotic diradical of the protein chromophore. This finding will ultimately support the development of technologies for searching novel fluorescent rhodopsin variants and unveil electrostatic changes that make light emission brighter and brighter.

Evidence of exciton-libron coupling in chirally adsorbed single molecules

Interplay between motion of nuclei and excitations has an important role in molecular photophysics of natural and artificial structures. Here we provide a detailed analysis of coupling between quantized librational modes (librons) and charged excited states (trions) on single phthalocyanine dyes adsorbed on a surface. By means of tip-induced electroluminescence performed with a scanning probe microscope, we identify libronic signatures in spectra of chirally adsorbed phthalocyanines and find that these signatures are absent from spectra of symmetrically adsorbed species. We create a model of the libronic coupling based on the Franck-Condon principle to simulate the spectral features. Experimentally measured librational spectra match very well the theoretically calculated librational eigenenergies and peak intensities (Franck-Condon factors). Moreover, the comparison reveals an unexpected depopulation channel for the zero libron of the excited state that can be effectively controlled by tuning the size of the nanocavity. Our results showcase the possibility of characterizing the dynamics of molecules by their low-energy molecular modes using µeV-resolved tip-enhanced spectroscopy.

Vibronic coupling in molecules plays an essential role in photophysics. Here, the authors observe optical fingerprints of the coupling between librational states and charged excited states in a single phthalocyanine molecule chirally absorbed on a surface.

Few-Femtosecond Isotope Effect in Polyatomic Molecules Ionized by Extreme Ultraviolet Attosecond Pulse Trains

Following ionization by an extreme ultraviolet (XUV) attosecond pulse train, a polyatomic molecule can be promoted to more-than-one excited states of the residual ion. The ensuing relaxation dynamics is often facilitated by several reaction coordinates, making them difficult to disentangle by the usual spectroscopic means. Here, we show that in atto-chemistry isotope labeling can be an efficient tool for unraveling the relaxation pathways in highly excited photoionized molecules. Employing an XUV pump pulse and a near-infrared probe pulse, we found the nuclear as well as coupled electron-nuclear dynamics in ethylene to be almost 40% faster compared to that of its deuterated counterpart. The findings, which are supported by advanced nonadiabatic dynamics calculations, led to the identification of the relevant nuclear coordinates controlling the relaxation. Our experiment highlights the relevance of ultrashort XUV pulses to capture the isotopic effect in few-femtosecond molecular photodynamics.

Rhodopsins: An Excitingly Versatile Protein Species for Research, Development and Creative Engineering

The first member and eponym of the rhodopsin family was identified in the 1930s as the visual pigment of the rod photoreceptor cell in the animal retina. It was found to be a membrane protein, owing its photosensitivity to the presence of a covalently bound chromophoric group. This group, derived from vitamin A, was appropriately dubbed retinal. In the 1970s a microbial counterpart of this species was discovered in an archaeon, being a membrane protein also harbouring retinal as a chromophore, and named bacteriorhodopsin. Since their discovery a photogenic panorama unfolded, where up to date new members and subspecies with a variety of light-driven functionality have been added to this family. The animal branch, meanwhile categorized as type-2 rhodopsins, turned out to form a large subclass in the superfamily of G protein-coupled receptors and are essential to multiple elements of light-dependent animal sensory physiology. The microbial branch, the type-1 rhodopsins, largely function as light-driven ion pumps or channels, but also contain sensory-active and enzyme-sustaining subspecies. In this review we will follow the development of this exciting membrane protein panorama in a representative number of highlights and will present a prospect of their extraordinary future potential.

Steady State Photoisomerization Quantum Yield of Model Rhodopsin: Insights from Wavepacket Dynamics?

We simulate the nonequilibrium steady state cis-trans photoisomerization of retinal chromophore in rhodopsin on the basis of a two-state, two-mode model coupled to a thermal environment. By analyzing the systematic trends within an inhomogeneously broadened ensemble of systems, we find that the steady state reaction quantum yield (QY) correlates strongly with the excess energy above the crossing point of the system, in agreement with the prediction of the short-time dynamical wavepacket picture. However, the nontrivial dependence of the QY on the system-environment interaction indicates that a pure dynamical picture is insufficient and that environment-induced partial internal energy redistribution takes place before the reaction concludes. These results imply that a proper treatment of the photoisomerization reaction, particularly its high QY, must account for the redistribution and dissipation of energy beyond the dynamical wavepacket motion that is typically employed in the literature and that is appropriate only in the transient regime.

Generalized Formulation of the Density Functional Tight Binding-Based Restricted Ensemble Kohn-Sham Method with Onsite Correction to Long-Range Correction

We present a generalized formulation for the combination of the density functional tight binding (DFTB) approach and the state-interaction state-average spin-restricted ensemble-referenced Kohn-Sham (SI-SA-REKS or SSR) method by considering onsite correction (OC) as well as the long-range corrected (LC) functional. The OC contribution provides more accurate energies and analytic gradients for individual microstates, while the multireference character of the SSR provides the correct description for conical intersections. We benchmark the LC-OC-DFTB/SSR method against various DFTB calculation methods for excitation energies and conical intersection structures with π/π* or n/π* characters. Furthermore, we perform excited-state molecular dynamics simulations with a molecular rotary motor with variations of LC-OC-DFTB/SSR approaches. We show that the OC contribution to the LC functional is crucial to obtain the correct geometry of conical intersections.

Independent Trajectory Mixed Quantum-Classical Approaches Based on the Exact Factorization

Mixed quantum-classical dynamics based on the exact factorization exploits the "derived" electron-nuclear correlation (ENC) term aiming for the description of quantum coherences. The ENC contains interactions between the phase of electronic states and nuclear quantum momenta which depend on the spatial shape of the nuclear density.The original surface hopping based on the exact factorization (SHXF) [\textit{J. Phys. Chem. Lett.} \textbf{2018}, \textit{9}, 1097] exploits frozen Gaussian functions to construct the nuclear density in the ENC term while the phase of electronic states is approximated as a fictitious nuclear momentum change.However, in reality, the width of nuclear wave packets varies in time depending on the shape of potential energy surfaces.In this work, we present a modified SHXF approach and a newly-developed Ehrenfest dynamics based on the exact factorization (EhXF) with time-dependent Gaussian functions and phases by enforcing total energy conservation.We perform numerical tests for various one-dimensional two-state model Hamiltonians.Overall, the time-dependent width of Gaussian functions and the energy conserving phase show a reliable decoherence compared to the original frozen Gaussian-based SHXF and the exact quantum mechanical calculation.Especially, the energy conserving phase is crucial for EhXF to reproduce the correct quantum dynamics.

Quantum-classical simulations of rhodopsin reveal excited-state population splitting and its effects on quantum efficiency

The activation of rhodopsin, the light-sensitive G-protein coupled receptor responsible for dim-light vision in vertebrates, is driven by an ultrafast excited-state double-bond isomerization with a quantum efficiency of almost 70%. The origin of such light sensitivity is not understood and a key question is whether in-phase nuclear motion controls the quantum efficiency value. Here, we use hundreds of quantum-classical trajectories to show that, 15 femtoseconds after light absorption, a degeneracy between the reactive excited state and a neighboring state causes the splitting of the rhodopsin population into subpopulations. These subpopulations propagate with different velocities and lead to distinct contributions to the quantum efficiency. We also show that such splitting is modulated by protein electrostatics, thus linking amino-acid sequence variations to quantum efficiency modulation. Finally, we discuss how such a linkage that in principle could be exploited to achieve higher quantum efficiencies, would simultaneously increase the receptor thermal noise leading to a trade-off that may have played a role in rhodopsin evolution.

Coherent Two-Dimensional and Broadband Electronic Spectroscopies

Over the past few decades, coherent broadband spectroscopy has been widely used to improve our understanding of ultrafast processes (e.g., photoinduced electron transfer, proton transfer, and proton-coupled electron transfer reactions) at femtosecond resolution. The advances in femtosecond laser technology along with the development of nonlinear multidimensional spectroscopy enabled further insights into ultrafast energy transfer and carrier relaxation processes in complex biological and material systems. New discoveries and interpretations have led to improved design principles for optimizing the photophysical properties of various artificial systems. In this review, we first provide a detailed theoretical framework of both coherent broadband and two-dimensional electronic spectroscopy (2DES). We then discuss a selection of experimental approaches and considerations of 2DES along with best practices for data processing and analysis. Finally, we review several examples where coherent broadband and 2DES were employed to reveal mechanisms of photoinitiated ultrafast processes in molecular, biological, and material systems. We end the review with a brief perspective on the future of the experimental techniques themselves and their potential to answer an even greater range of scientific questions.

Analytical Nonadiabatic Couplings and Gradients within the State-Averaged Orbital-Optimized Variational Quantum Eigensolver

We introduce several technical and analytical extensions to our recent state-averaged orbital-optimized variational quantum eigensolver (SA-OO-VQE) algorithm (see Yalouz et al. Quantum Sci. Technol. 2021, 6, 024004). Motivated by the limitations of current quantum computers, the first extension consists of an efficient state-resolution procedure to find the SA-OO-VQE eigenstates, and not just the subspace spanned by them, while remaining in the equi-ensemble framework. This approach avoids expensive intermediate resolutions of the eigenstates by postponing this problem to the very end of the full algorithm. The second extension allows for the estimation of analytical gradients and nonadiabatic couplings, which are crucial in many practical situations ranging from the search of conical intersections to the simulation of quantum dynamics, in, for example, photoisomerization reactions. The accuracy of our new implementations is demonstrated on the formaldimine molecule CH₂NH (a minimal Schiff base model relevant for the study of photoisomerization in larger biomolecules), for which we also perform a geometry optimization to locate a conical intersection between the ground and first-excited electronic states of the molecule.

Light Driven Ultrafast Bioinspired Molecular Motors: Steering and Accelerating Photoisomerization Dynamics of Retinal

Photoisomerization of retinal protonated Schiff base in microbial and animal rhodopsins are strikingly ultrafast and highly specific. Both protein environments provide conditions for fine-tuning the photochemistry of their chromophores. Here, by combining time-resolved action absorption spectroscopy and high-level electronic structure theory, we show that similar control can be gained in a synthetically engineered retinal chromophore. By locking the dimethylated retinal Schiff base at the C11═C12 double bond in its trans configuration (L-RSB), the excited-state decay is rendered from a slow picosecond to an ultrafast subpicosecond regime in the gas phase. Steric hindrance and pretwisting of L-RSB are found to be important for a significant reduction in the excited-state energy barriers, where isomerization of the locked chromophore proceeds along C9═C10 rather than the preferred C11═C12 isomerization path. Remarkably, the accelerated excited-state dynamics also becomes steered. We show that L-RSB is capable of unidirectional 360° rotation from all-trans to 9-cis and from 9-cis to all-trans in only two distinct steps induced by consecutive absorption of two 600 nm photons. This opens a way for the rational design of red-light-driven ultrafast molecular rotary motors based on locked retinal chromophores.

Internal conversion of the anionic GFP chromophore: in and out of the I-twisted S1/S0 conical intersection seam

The functional diversity of the green fluorescent protein (GFP) family is intimately connected to the interplay between competing photo-induced transformations of the chromophore motif, anionic p-hydroxybenzylidene-2,3-dimethylimidazolinone (HBDI^-). Its ability to undergo Z/E-isomerization is of particular importance for super-resolution microscopy and emerging opportunities in optogenetics. Yet, key dynamical features of the underlying internal conversion process in the native HBDI^- chromophore remain largely elusive. We investigate the intrinsic excited-state behavior of isolated HBDI^- to resolve competing decay pathways and map out the factors governing efficiency and the stereochemical outcome of photoisomerization. Based on non-adiabatic dynamics simulations, we demonstrate that non-selective progress along the two bridge-torsional (i.e., phenolate, P, or imidazolinone, I) pathways accounts for the three decay constants reported experimentally, leading to competing ultrafast relaxation primarily along the I-twisted pathway and S₁ trapping along the P-torsion. The majority of the population (∼70%) is transferred to S₀ in the vicinity of two approximately enantiomeric minima on the I-twisted intersection seam (MECI-Is). Despite their sloped, reactant-biased topographies (suggesting low photoproduct yields), we find that decay through these intersections leads to products with a surprisingly high quantum yield of ∼30%. This demonstrates that E-isomer generation results at least in part from direct isomerization on the excited state. A photoisomerization committor analysis reveals a difference in intrinsic photoreactivity of the two MECI-Is and that the observed photoisomerization is the combined result of two effects: early, non-statistical dynamics around the less reactive intersection followed by later, near-statistical behavior around the more reactive MECI-I. Our work offers new insight into internal conversion of HBDI^- that both establishes the intrinsic properties of the chromophore and enlightens principles for the design of chromophore derivatives and protein variants with improved photoswitching properties.

Torsionally broken symmetry assists infrared excitation of biomimetic charge-coupled nuclear motions in the electronic ground state

The concerted interplay between reactive nuclear and electronic motions in molecules actuates chemistry. Here, we demonstrate that out-of-plane torsional deformation and vibrational excitation of stretching motions in the electronic ground...

Nonadiabatic Dynamics with Coupled Trajectories

In this paper, we discuss coupled-trajectory schemes for molecular-dynamics simulations of excited-state processes. New coupled-trajectory strategies to capture decoherence effects, revival of coherence and nonadiabatic interferences in long-time dynamics are proposed, and compared to independent-trajectory schemes. The working framework is provided by the exact factorization of the electron-nuclear wave function, and it exploits ideas emanating from various surface-hopping schemes. The new coupled-trajectory algorithms are tested on a one-dimensional two-state system using different model parameters which allow one to induce different dynamics. The benchmark is provided by the numerically exact solution of the time-dependent Schrödinger equation.

Deuterated Covalent Organic Frameworks with Significantly Enhanced Luminescence

Luminescent covalent organic frameworks (COFs) find promising applications in chemical sensing, photocatalysis, and optoelectronic devices, however, the majority of COFs are non or weakly emissive owing to the aggregation-caused quenching (ACQ) or the molecular thermal motion-based energy dissipation. Here, we report a previously unperceived approach to improve luminescence performance of COFs by introducing isotope effect, which is achieved through substitution of hydrogen from high-frequency oscillators X-H (X=O, N, C) by heavier isotope deuterium. Combining the "bottom-up" and in situ deuteration methods generates the first deuterated COF, which exhibits an impressively 19-fold enhancement in quantum yield over that of the non-deuterated counterpart. These results are interpreted by theoretical calculations as the consequence of slower C/N-D and OD⋅⋅⋅O vibrations that impede the nonradiative deactivation process. The proposed strategy is proved applicable to many other types of emissive COFs.

Insights into the Early-Time Excited-State Dynamics of Structurally Inhomogeneous Rhodopsin KR2

The light-driven sodium-pump rhodopsin KR2 exhibits ultrafast photoisomerization dynamics of its all-trans protonated Schiff-base retinal (PSBR). However, the excited-state decay of KR2 also shows slow picosecond time constants, which are attributed to nonreactive states. The mechanism that produces long-lived states is unclear. Here, by using molecular dynamics simulations and large-scale XMCQDPT2-based QM/MM modeling, we explore the origin of reactive and nonreactive states in KR2. By calculating the S₀-S₁ vibronic band shapes, we gain insight into the early-time excited-state dynamics of PSBR and show that the protein environment can significantly alter vibrational modes that are active upon photoexcitation, thus facilitating photoisomerization from all-trans to 13-cis PSBR. Importantly, we reveal structural heterogeneity of the retinal-binding pocket of KR2, characterized by three distinct conformations, and conclude that the formation of a strong hydrogen bond between the retinal Schiff base and its counterion is essential for the ultrafast reaction dynamics.

PyUNIxMD: A Python-based excited state molecular dynamics package

Theoretical/computational description of excited state molecular dynamics is nowadays a crucial tool for understanding light-matter interactions in many materials. Here we present an open-source Python-based nonadiabatic molecular dynamics program package, namely PyUNIxMD, to deal with mixed quantum-classical dynamics for correlated electron-nuclear propagation. The PyUNIxMD provides many interfaces for quantum chemical calculation methods with commercial and noncommercial ab initio and semiempirical quantum chemistry programs. In addition, the PyUNIxMD offers many nonadiabatic molecular dynamics algorithms such as fewest-switch surface hopping and its derivatives as well as decoherence-induced surface hopping based on the exact factorization (DISH-XF) and coupled-trajectory mixed quantum-classical dynamics (CTMQC) for general purposes. Detailed structures and flows of PyUNIxMD are explained for the further implementations by developers. We perform a nonadiabatic molecular dynamics simulation for a molecular motor system as a simple demonstration.

Effect of charge-transfer enhancement on the efficiency and rotary mechanism of an oxindole-based molecular motor

Harvesting energy and converting it into mechanical motion forms the basis for both natural and artificial molecular motors. Overcrowded alkene-based light-driven rotary motors are powered through sequential photochemical and thermal steps. The thermal helix inversion steps are well characterised and can be manipulated through adjustment of the chemical structure, however, the insights into the photochemical isomerisation steps still remain elusive. Here we report a novel oxindole-based molecular motor featuring pronounced electronic push-pull character and a four-fold increase of the photoisomerization quantum yield in comparison to previous motors of its class. A multidisciplinary approach including synthesis, steady-state and transient absorption spectroscopies, and electronic structure modelling was implemented to elucidate the excited state dynamics and rotary mechanism. We conclude that the charge-transfer character of the excited state diminishes the degree of pyramidalisation at the alkene bond during isomerisation, such that the rotational properties of this oxindole-based motor stand in between the precessional motion of fluorene-based molecular motors and the axial motion of biomimetic photoswitches.

Extreme Parametric Sensitivity in the Steady-State Photoisomerization of Two-Dimensional Model Rhodopsin

We computationally studied the photoisomerization reaction of the retinal chromophore in rhodopsin using a two-state two-mode model coupled to thermal baths. Reaction quantum yields at the steady state (10 ps and beyond) were found to be considerably different than their transient values, suggesting a weak correlation between transient and steady-state dynamics in these systems. Significantly, the steady-state quantum yield was highly sensitive to minute changes in system parameters, while transient dynamics was nearly unaffected. Correlation of such sensitivity with standard level spacing statistics of the nonadiabatic vibronic system suggests a possible origin in quantum chaos. The significance of this observation of quantum yield parametric sensitivity in biological models of vision has profound conceptual and fundamental implications.

Frontiers in Multiscale Modeling of Photoreceptor Proteins

This perspective article highlights the challenges in the theoretical description of photoreceptor proteins using multiscale modeling, as discussed at the CECAM workshop in Tel Aviv, Israel. The participants have identified grand challenges and discussed the development of new tools to address them. Recent progress in understanding representative proteins such as green fluorescent protein, photoactive yellow protein, phytochrome, and rhodopsin is presented, along with methodological developments.

Intermolecular conical intersections in molecular aggregates

Conical intersections (CoIns) of multidimensional potential energy surfaces are ubiquitous in nature and control pathways and yields of many photo-initiated intramolecular processes. Such topologies can be potentially involved in the energy transport in aggregated molecules or polymers but are yet to be uncovered. Here, using ultrafast two-dimensional electronic spectroscopy (2DES), we reveal the existence of intermolecular CoIns in molecular aggregates relevant for photovoltaics. Ultrafast, sub-10-fs 2DES tracks the coherent motion of a vibrational wave packet on an optically bright state and its abrupt transition into a dark state via a CoIn after only 40 fs. Non-adiabatic dynamics simulations identify an intermolecular CoIn as the source of these unusual dynamics. Our results indicate that intermolecular CoIns may effectively steer energy pathways in functional nanostructures for optoelectronics.

Quantum and Quantum-Classical Studies of the Photoisomerization of a Retinal Chromophore Model

We report an in-depth analysis of the photo-induced isomerization of the 2-cis-penta-2,4-dieniminium cation: a minimal model of the 11-cis retinal protonated Schiff base chromophore of the dim-light photoreceptor rhodopsin. Based on recently developed three-dimensional potentials parametrized on ab initio multi-state multi-configurational second-order perturbation theory data, we perform quantum-dynamical studies. In addition, simulations based on various quantum-classical methods, among which Tully surface hopping and the coupled-trajectory approach derived from the exact factorization, allow us to validate their performance against vibronic wavepacket propagation and, therefore, a purely quantum treatment. Quantum-dynamics results uncover qualitative differences with respect to the two-dimensional Hahn-Stock potentials, widely used as model potentials for the isomerization of the same chromophore, due to the increased dimensionality and three-mode correlation. Quantum-classical simulations show, instead, that three-dimensional model potentials are capable of capturing a number of features revealed by atomistic simulations and experimental observations. In particular, a recently reported vibrational phase relationship between double-bond torsion and hydrogen-out-of-plane modes critical for rhodopsin isomerization efficiency is correctly reproduced.

A systematic model study quantifying how conical intersection topography modulates photochemical reactions

Despite their important role in photochemistry and expected presence in most polyatomic molecules, conical intersections have been thoroughly characterized in a comparatively small number of systems. Conical intersections can confer molecular photoreactivity or photostability, often with remarkable efficacy, due to their unique structure: at a conical intersection, the adiabatic potential energy surfaces of two or more electronic states are degenerate, enabling ultrafast decay from an excited state without radiative emission, known as nonadiabatic transfer. Furthermore, the precise conical intersection topography determines fundamental properties of photochemical processes, including excited-state decay rate, efficacy, and molecular products that are formed. However, these relationships have yet to be defined comprehensively. In this article, we use an adaptable computational model to investigate a variety of conical intersection topographies, simulate resulting nonadiabatic dynamics, and calculate key photochemical observables. We varied the vibrational mode frequencies to modify conical intersection topography systematically in four primary classes of conical intersections and quantified the resulting rate, total yield, and product yield of nonadiabatic decay. The results reveal that higher vibrational mode frequencies reduce nonadiabatic transfer, but increase the transfer rate and resulting photoproduct formation. These trends can inform progress toward experimental control of photochemical reactions or tuning of molecules' photochemical properties based on conical intersections and their topography.

Molecular Photochemistry: Recent Developments in Theory

Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.

Uncovering the Functional Role of Coherent Phonons during the Photoinduced Phase Transition in a Molecular Crystal

The atomic motions that make up phonons and molecular vibrations in molecular crystals influence their photophysical and electronic properties, including polaron formation, carrier mobility, and phase transitions. Discriminating between spectator and driving motions is a significant challenge hindering optimization. Unlocking this information and developing fine-tuned controls over actively participating phonon modes would not only lead to a stronger understanding of photochemistry but also provide a significant new tool in controlling solid state chemistry. We present a strategy using rationally designed double pulses to unveil the unique function of specific excited state phonon modes. Using ultrafast spectroscopy, we identified 50 and 90 cm^-1 phonons involved in modulating the photoinduced spin-Peierls melting of potassium tetracyanoquinodimethane crystals. We show that the 50 cm^-1 phonon specifically corresponds to the coherent nuclear wavepacket involved in the charge transfer component of the overall spin-Peierls phase melting process, while the 90 cm^-1 phonon facilitates the phase transition component.

Modern quantum chemistry with [Open]Molcas

MOLCAS/OpenMolcas is an ab initio electronic structure program providing a large set of computational methods from Hartree-Fock and density functional theory to various implementations of multiconfigurational theory. This article provides a comprehensive overview of the main features of the code, specifically reviewing the use of the code in previously reported chemical applications as well as more recent applications including the calculation of magnetic properties from optimized density matrix renormalization group wave functions.

Combining Graphics Processing Units, Simplified Time-Dependent Density Functional Theory, and Finite-Difference Couplings to Accelerate Nonadiabatic Molecular Dynamics

Starting from our recently published implementation of nonadiabatic molecular dynamics (NAMD) on graphics processing units (GPUs), we explore further approaches to accelerate ab initio NAMD calculations at the time-dependent density functional theory (TDDFT) level of theory. We employ (1) the simplified TDDFT schemes of Grimme et al. and (2) the Hammes-Schiffer-Tully approach to obtain nonadiabatic couplings from finite-difference calculations. The resulting scheme delivers an accurate physical picture while virtually eliminating the two computationally most demanding steps of the algorithm. Combined with our GPU-based integral routines for SCF, TDDFT, and TDDFT derivative calculations, NAMD simulations of systems of a few hundreds of atoms at a reasonable time scale become accessible on a single compute node. To demonstrate this and to present a first, illustrative example, we perform TDDFT/MM-NAMD simulations of the rhodopsin protein.

Effect of Oriented External Electric Fields on the Photo and Thermal Isomerization of Azobenzene

Azobenzene is a prototype molecule with potential applications in molecular switches, solar thermal batteries, sensors, photoresponsive membranes, molecular electronics, data storage, and nonlinear optics. Photo and thermal isomerization pathways exhibit different charge-transfer character and dipole moments, implying that the use of electric fields can be used to modulate the reactivity of azobenzene. This article examines the differential effect of orientated electric fields on the rotation and inversion thermal and photoisomerization pathways of azobenzene to explore the feasibility of using electric fields in the design of azobenzene-based molecular devices. Our findings demonstrate that the application of orientated electric fields modifies the accessibility of the S₀/S₁ seam of electronic degeneracy, as well as changes the energetically favored relaxation pathway in the branching space to yield different photoproducts. In addition, we observed strong-field dipole-inversion effects that cause a topographical change in the response of the potential energy surface to the applied field and can result in geometric minima that do not exist under field-free conditions. On the S₀ surface, transition barriers can be modified on the order of ±10 kcal mol^-1, enabling control of thermal isomerization rates.

Mapping Structural Dynamics of Proteins with Femtosecond Stimulated Raman Spectroscopy

The structure-function relationships of biomolecules have captured the interest and imagination of the scientific community and general public since the field of structural biology emerged to enable the molecular understanding of life processes. Proteins that play numerous functional roles in cellular processes have remained in the forefront of research, inspiring new characterization techniques. In this review, we present key theoretical concepts and recent experimental strategies using femtosecond stimulated Raman spectroscopy (FSRS) to map the structural dynamics of proteins, highlighting the flexible chromophores on ultrafast timescales. In particular, wavelength-tunable FSRS exploits dynamic resonance conditions to track transient-species-dependent vibrational motions, enabling rational design to alter functions. Various ways of capturing excited-state chromophore structural snapshots in the time and/or frequency domains are discussed. Continuous development of experimental methodologies, synergistic correlation with theoretical modeling, and the expansion to other nonequilibrium, photoswitchable, and controllable protein systems will greatly advance the chemical, physical, and biological sciences.

Non-adiabatic Excited-State Molecular Dynamics: Theory and Applications for Modeling Photophysics in Extended Molecular Materials

Optically active molecular materials, such as organic conjugated polymers and biological systems, are characterized by strong coupling between electronic and vibrational degrees of freedom. Typically, simulations must go beyond the Born-Oppenheimer approximation to account for non-adiabatic coupling between excited states. Indeed, non-adiabatic dynamics is commonly associated with exciton dynamics and photophysics involving charge and energy transfer, as well as exciton dissociation and charge recombination. Understanding the photoinduced dynamics in such materials is vital to providing an accurate description of exciton formation, evolution, and decay. This interdisciplinary field has matured significantly over the past decades. Formulation of new theoretical frameworks, development of more efficient and accurate computational algorithms, and evolution of high-performance computer hardware has extended these simulations to very large molecular systems with hundreds of atoms, including numerous studies of organic semiconductors and biomolecules. In this Review, we will describe recent theoretical advances including treatment of electronic decoherence in surface-hopping methods, the role of solvent effects, trivial unavoided crossings, analysis of data based on transition densities, and efficient computational implementations of these numerical methods. We also emphasize newly developed semiclassical approaches, based on the Gaussian approximation, which retain phase and width information to account for significant decoherence and interference effects while maintaining the high efficiency of surface-hopping approaches. The above developments have been employed to successfully describe photophysics in a variety of molecular materials.

Observe while it happens: catching photoactive proteins in the act with non-adiabatic molecular dynamics simulations

Organisms use photo-receptors to react to light. The first step is usually the absorption of a photon by a prosthetic group embedded inside the photo-receptor, often a conjugated chromophore. The electronic changes in the chromophore induced by photo-absorption can trigger a cascade of structural or chemical transformations that culminate into a response to light. Understanding how these proteins have evolved to mediate their activation process has remained challenging because the required time and spacial resolutions are notoriously difficult to achieve experimentally. Therefore, mechanistic insights into photoreceptor activation have been predominantly obtained with computer simulations. Here we briefly outline the challenges associated with such computations and review the progress made in this field.

Vibrational coherence and quantum yield of retinal-chromophore-inspired molecular switches

UV-Vis transient absorption (TA) spectroscopy is used to carry out a systematic investigation of the ultrafast CC double photoisomerization dynamics and quantum yield of each isomer of a set of six chromophores based on the same retinal-inspired, indanylidene pyrrolinium (IP) molecular framework.

Photoinduced Electron Transfer in Donor-Acceptor Complexes: Isotope Effect and Dynamic Symmetry Breaking

Electron-nuclear (vibronic) coupling has emerged as an important factor in determining the efficiency of energy transfer and charge separation in natural and artificial photosynthetic systems. Here we investigate the photoinduced charge-transfer process in a hydrogen-bonded donor-acceptor molecular complex. By using real-time quantum-classical simulations based on time-dependent Kohn-Sham equations, we follow in detail the relaxation from the Franck-Condon point to the region of strong nonadiabatic coupling where electron transfer occurs. We elucidate how the charge transfer is coupled to specific vibrational modes and how it is affected by isotope substitution. The importance of resonance in nuclear and electron dynamics and the role of dynamic symmetry breaking are emphasized. Using the dipole moment as a descriptive parameter, exchange of angular momentum between nuclear and electronic subsystems in an electron-nuclear resonant process is inferred. The performed simulations support a nonadiabatic conversion via adiabatic passage process that was recently put forward. These results are relevant in deriving rational design principles for solar-to-fuel conversion devices.