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

Landau-Zener transitions with spin flipping in chiral structures

We study Landau-Zener transitions with the spin flipping of electron states in donor-chiral bridge-acceptor structures due to electron spin coupling with the strain tensor induced by the anisotropic phonon field of the chiral bridge molecule, where the transition rates near conical intersections between two potential energy surfaces for the right-handed and left-handed bridge structures are calculated. We find that the flipping rates sensitively depend on the handedness of the bridge molecules and could be significantly enhanced with the increasing of phonon wave vectors, displaying the crucial rules of the electron-phonon coupling for the spin selectivity effect in chiral structures. These results not only shed light on the microscopic mechanisms of chirality-induced spin selectivity but also enlighten the role of spin selectivity in electron transfer in chemical and biological systems.

Tuning the Photoisomerization Mechanism of Oxindole Switches with Electron-Donating Substituents

This study investigates the photoreaction mechanism of a hydroxy-substituted oxindole photoswitch using femtosecond transient absorption, fluorescence up-conversion, and computational chemistry. Deprotonation of the hydroxyl group enhances the push-pull character in the molecule, allowing tuning of the photoisomerization mechanism from a precessional to an axial motion. The neutral form of the switch exhibits longer excited-state lifetimes, while the anionic form decays rapidly within 200 fs. Computational models show that deprotonation increases the charge transfer and accessibility to conical intersections. This work highlights how varying the electron-donating strength of a substituent in a push-pull photoswitch tunes the photoreaction mechanism in designing photoswitches.

Rational Construction of a Pyrimido[2,1-b] Benzothiazole-Based Photoswitchable Smart AIE Material: A Theoretical Insight into Fluorescence Switching Mechanism of the Chimeric Dyad

Currently, AIEgen-photochromics conjugates and intrinsic photochromic AIEgens are the two major molecular design strategies for photoresponsive AIE materials. However, these two strategies still have their own limitations. In our previous research, we discovered pyrimido[2,1-b][1,3]benzothiazole (PBT) as a novel core structure of AIEgens. We herein explored a chimeric strategy to rationally integrate photochromic bisthienylethene (BTE) as the orthogonal head group into a new PBT molecule (PBTE). Compared to the conventional design strategies, the compact chimeric design of PBTE not only well maintains the AIE and photochromic properties of the parent fragments but also leads to excellent AIE-photoswitching capability in both films and single crystals. Theoretical calculations revealed that MO energy level arrangement of PBTE reorganized upon photoisomerization and the fast vibrational relaxation (VR) along with internal conversion (IC) from S5 to S1 state may serve as a competitive channel for the fluorescence quenching of PBTE-c. Comparison of the kinetics of the nonradiative decay with those of the excited-state energy transfer (EET) processes clearly showed that the ultrafast intramolecular Förster resonance energy transfer (FRET) is the dominant cause of fluorescence-off state. The applications of PBTE in erasable optical memory material and multi-dimensional anti-counterfeiting have also been demonstrated.

Spectroscopy and Excited-State Dynamics of Methyl Ferulate in Molecular Beams

The spectroscopic and dynamic properties of methyl ferulate─a naturally occurring ultraviolet-protecting filter─and microsolvated methyl ferulate have been studied under molecular beam conditions using resonance-enhanced multiphoton ionization spectroscopy in combination with quantum chemical calculations. We demonstrate and rationalize how the phenyl substitution pattern affects the state ordering of the lower excited singlet state manifold and what the underlying reason is for the conformation-dependent Franck-Condon (FC) activity in the UV-excitation spectra. Studies on microsolvated methyl ferulate reveal potential coordination sites and the influence of such coordination on the spectroscopic properties. Our quantum chemical studies also enable us to obtain a quantitative understanding of the dominant excited-state decay routes of the photoexcited ππ* state involving a ∼3 ns intersystem crossing pathway to the triplet manifold─which is much slower than found for coumarates─and a relatively fast intersystem crossing back to the ground state (∼30 ns). We show that a common T1/S0 crossing can very well explain the observation that T1 lifetimes are quasi-independent of the phenyl substitution pattern.

Vibronic Conical Intersection Trajectory Signatures in Wave Packet Coherences

Conical intersections are ubiquitous in the energy landscape of chemical systems, drive photochemical reactivity, and are extremely challenging to observe spectroscopically. Using two-dimensional electronic spectroscopy, we observe the nonadiabatic dynamics in Wurster's Blue after excitation to the lowest two vibronic excited states. The excited populations relax ballistically through a conical intersection in 55 fs to the electronic ground state potential energy surface as the molecule undergoes an intramolecular electron transfer. While the kinetics are identical on both vibronic energy surfaces, we observe different patterns of coherent oscillations after traversing the conical intersection indicating distinct nonadiabatic relaxation pathways through the conical energetic funnel. These coherences are not created directly by the excitation pulses but are the result of the dynamical trajectories projecting differently on the conical intersection vibrational space. Our spectroscopic data offers a fresh perspective into the complex conical intersection topology and dynamics that emphasizes the critical involvement of the intersection space in dictating the dynamics.

Complementing Adiabatic and Nonadiabatic Methods To Understand Internal Conversion Dynamics in Porphyrin Derivatives

We analyze the internal conversion dynamics within the Qy and Qx excited states of both bare and functionalized porphyrins, which are known to exhibit significantly different time constants experimentally. Through the integration of two complementary approaches, static calculation of per-mode reorganization energies and nonadiabatic molecular dynamics, we achieve a comprehensive understanding of the factors determining the different behavior of the two molecules. We identify the key normal and essential modes responsible for the population transfer between excited states and discuss the efficacy of different statistical and nonstatistical analyses in providing a full physics-based description of the phenomenon.

Study of Photoselectivity in Linear Conjugated Chromophores Using the XMS-CASPT2 Method

Photoisomerization, the structural alteration of molecules upon absorption of light, is crucial for the function of biological chromophores such as retinal in opsins, proteins vital for vision and other light-sensitive processes. The intrinsic selectivity of this isomerization process (i.e., which double bond in the chromophore is isomerized) is governed by both the inherent properties of the chromophore and its surrounding environment. In this study, we employ the extended multistate complete active space second-order perturbation theory (XMS-CASPT2) method to investigate photoisomerization selectivity in linear conjugated chromophores, focusing on two simple molecular models resembling retinal. By analyzing electronic energies, intramolecular charge separation, and conical intersection topographies in the gas phase, we show that the photoproduct formed by rotation around the double bond near the Schiff base is energetically favored. The topographic differences at the conical intersections leading to different photoproducts reveal differences in photodynamics. In multiphoton excitation, the primary photoproduct typically reverts to the initial configuration rather than rotating around a different double bond. Our study offers new insights into the photodynamics of photoisomerizing double bonds in π-conjugated chromophores. We anticipate that our findings will provide valuable perspectives for advancing the understanding of biological chromophores and for designing efficient photochemical switches with applications in molecular electronics and phototherapy.

Nondirect-Product Local Diabatic Representation with Smolyak Sparse Grids

Modeling nonadiabatic conical intersection dynamics is critical for understanding a wide range of photophysical, photochemical, and biological phenomena. Here we develop a nonadiabatic conical intersection wave packet dynamic method in the local diabatic representation using Smolyak sparse grids. Employing sparse grids avoids the direct-product grids in configuration space and alleviates the exponential scaling of computation costs with the molecular size. Numerical demonstrations are first performed for a two-dimensional vibronic model of pyrazine, where the results using sparse grids are in excellent agreement with those using direct-product grids, with sparse grids being much faster. Moreover, we demonstrate that for a four-dimensional pyrazine model, where direct-product grids are computationally infeasible, sparse grids can provide almost exact results. The sparse grid local diabatic representation method is further applied to a realistic model system of phenol photodissociation with much more complex potential energy surfaces; the results using sparse grids still agree very well with the direct-product grids. Finally, by combining with electronic structure calculations, we apply our method to the Shin-Metiu model without quasi-diabatization. The sparse grid and direct-product grid results are in good agreement, with the sparse grid computational cost being half of the full grid.

Reversible Photochromic Reactions of Bacteriorhodopsin from Halobacterium salinarum at Femto- and Picosecond Times

The operation of bacteriorhodopsin (BR) from the archaeon Halobacterium salinarum is based on the photochromic reaction of isomerization of the chromophore group (the retinal protonated Schiff base, RPSB) from the all-trans to the 13-cis form. The ultrafast dynamics of the reverse 13-cis → all-trans photoreaction was studied using femtosecond transient absorption spectroscopy in comparison with the forward photoreaction. The forward photoreaction was initiated by photoexcitation of BR by pulse I (540 nm). The reverse photoreaction was initiated by photoexcitation of the product K590 at an early stage of its formation (5 ps) by pulse II (660 nm). The conversion of the excited K590 to the ground state proceeds at times of 0.19, 1.1, and 16 ps with the relative contributions of ~20/60/20, respectively. All these decay channels lead to the formation of the initial state of BR as a product with a quantum yield of ~1. This state is preceded by vibrationally excited intermediates, the relaxation of which occurs in the 16 ps time range. Likely, the heterogeneity of the excited state of K590 is determined by the heterogeneity of its chromophore center. The forward photoreaction includes two components-0.52 and 3.5 ps, with the relative contributions of 91/9, respectively. The reverse photoreaction initiated from K590 proceeds more efficiently in the conical intersection (CI) region but on the whole at a lower rate compared to the forward photoreaction, due to significant heterogeneity of the potential energy surface.

Real-time capture of nuclear motions influencing photoinduced electron transfer

Although vibronic coupling phenomena have been recognized in the excite state dynamics of transition metal complexes, its impact on photoinduced electron transfer (PET) remains largely unexplored. This study investigates coherent wavepacket (CWP) dynamics during PET processes in a covalently linked electron donor-acceptor complex featuring a cyclometalated Pt(ii) dimer as the donor and naphthalene diimide (NDI) as the acceptors. Upon photoexciting the Pt(ii) dimer electron donor, ultrafast broadband transient absorption spectroscopy revealed direct modulation of NDI radical anion formation through certain CWP motions and correlated temporal evolutions of the amplitudes for these CWPs with the NDI radical anion formation. These results provide clear evidence that the CWP motions are the vibronic coherences coupled to the PET reaction coordinates. Normal mode analysis identified that the CWP motions originate from vibrational modes associated with the dihedral angles and bond lengths between the planes of the cyclometalating ligand and the NDI, the key modes altering their π-interaction, consequently influencing PET dynamics. The findings highlight the pivotal role of vibrations in shaping the favorable trajectories for the efficient PET processes.

Structure, Dynamics, and Reactivity of Encapsulated Molecules in Restricted Spaces: Arylazoisoxazoles within an Octa Acid Capsule

In this study, a well-defined organic capsule assembled from two octa acid (OA) molecules acting as host and select arylazoisoxazoles (AAIO) acting as guests were employed to demonstrate that confined molecules have restricted freedom that translates into reaction selectivity in both ground and excited states. The behavior of these AAIO guests in confined capsules was found to be different from that found in both crystals, where there is very little freedom, and in isotropic solvents, where there is complete freedom. Through one-dimensional (1D) and two-dimensional (2D) 1H NMR spectroscopic experiments, we have established a relationship between structure, dynamics and reactivity of molecules confined in an OA capsule. Introduction of CF3 and CH3 substitution at the 4-position of the aryl group of AAIO reveals that in addition to space confinement, weak interactions between the guest and the OA capsule control the dynamics and reactivity of guest molecules. 1H NMR studies revealed that there is a temperature-dependence to guest molecules tumbling (180° rotation along the capsular short axis) within an OA capsule. While 1H NMR points to the occurrence of tumbling motion, MD simulations and simulation of the temperature-dependent NMR signals provide an insight into the mechanism of tumbling within OA capsules. Thermal and photochemical isomerization of AAIO were found to occur within an OA capsule just as in organic solvents. The observed selectivity noted during thermal and photo induced isomerization of OA encapsulated AAIOs can be qualitatively understood in terms of the well-known concepts due to Bell-Evans-Polanyi (BEP principle), Hammond and Zimmerman.

In-silico predicted mouse melanopsins with blue spectral shifts deliver efficient subcellular signaling

Melanopsin is a photopigment belonging to the G Protein-Coupled Receptor (GPCR) family expressed in a subset of intrinsically photosensitive retinal ganglion cells (ipRGCs) and responsible for a variety of processes. The bistability and, thus, the possibility to function under low retinal availability would make melanopsin a powerful optogenetic tool. Here, we aim to utilize mouse melanopsin to trigger macrophage migration by its subcellular optical activation with localized blue light, while simultaneously imaging the migration with red light. To reduce melanopsin's red light sensitivity, we employ a combination of in silico structure prediction and automated quantum mechanics/molecular mechanics modeling to predict minimally invasive mutations to shift its absorption spectrum towards the shorter wavelength region of the visible spectrum without compromising the signaling efficiency. The results demonstrate that it is possible to achieve melanopsin mutants that resist red light-induced activation but are activated by blue light and display properties indicating preserved bistability. Using the A333T mutant, we show that the blue light-induced subcellular melanopsin activation triggers localized PIP3 generation and macrophage migration, which we imaged using red light, demonstrating the optogenetic utility of minimally engineered melanopsins.

Tracking the conical intersection dynamics for the photoinduced Jahn-Teller switch of a Mn(iii) complex

Octahedral Mn(iii) ions predominantly exhibit an axially elongated Jahn-Teller (JT) distortion, which is responsible for their large uniaxial magnetic anisotropy. As a result, they are often used in the synthesis of single-molecule magnets (SMMs). Modulation of the JT distortion using femtosecond laser pulses could offer a route to controlling the magnetisation direction, and therefore is promising for the development of data storage devices that work on ultrafast timescales. Photoinduced switching of the distortion from an axially elongated to an axially compressed structure has been demonstrated for various Mn(iii) complexes. However, the dynamics around the region of the conical intersection for the photoinduced JT switch remains unclear. Here, ultrafast transient absorption spectra were recorded for solutions of tris(2,2,6,6-tetramethyl-3,5-heptanedionato)manganese(iii) (Mn(dpm)3) in ethanol to further explore the dynamics of the photoinduced JT switch. We observe the generation of a vibrational wavepacket on the excited state surface, which has a frequency of approximately 155 cm-1 and encompasses a JT-active vibrational mode. This coherent motion is maintained after passage through the conical intersection back to the ground state, which launches wavepackets along the ground state potential energy surface (PES) with frequencies of approximately 180 and 110 cm-1 that we assign to the elongated and compressed state, respectively. Inspection of the relative phases of the frequencies reveals phase shifts that are consistent with a one-mode reaction coordinate, and passes through the conical intersection at 1/4 and 3/4 of the excited state vibrational period. Our results provide direct insights into the non-adiabatic dynamics of Mn(iii) complexes, which can be used to guide the synthesis of optically controlled SMMs.

Excited State Dynamics in Unidirectional Photochemical Molecular Motors

Unidirectional photochemically driven molecular motors (PMMs) convert the energy of absorbed light into continuous rotational motion. As such they are key components in the design of molecular machines. The prototypical and most widely employed class of PMMs is the overcrowded alkenes, where rotational motion is driven by successive photoisomerization and thermal helix inversion steps. The efficiency of such PMMs depends upon the speed of rotation, determined by the rate of ground state thermal helix inversion, and the quantum yield of photoisomerization, which is dependent on the excited state energy landscape. The former has been optimized by synthetic modification across three generations of overcrowded alkene PMMs. These improvements have often been at the expense of photoisomerization yield, where there remains room for improvement. In this perspective we review the application of ultrafast spectroscopy to characterize the excited state dynamics in PMMs. These measurements lead to a general mechanism for all generations of PMMs, involving subpicosecond decay of a Franck-Condon excited state to populate a dark excited state which decays within picoseconds via conical intersections with the electronic ground state. The model is discussed in the context of excited state dynamics calculations. Studies of PMM photochemical dynamics as a function of solvent suggest exploitation of intramolecular charge transfer and solvent polarity as a route to controlling photoisomerization yield. A test of these ideas for a first generation motor reveals a high degree of solvent control over isomerization yield. These results suggest a pathway to fine control over the performance of future PMMs.

Exploring Exact-Factorization-Based Trajectories for Low-Energy Dynamics near a Conical Intersection

We study low-energy dynamics generated by a two-dimensional two-state Jahn-Teller Hamiltonian in the vicinity of a conical intersection using quantum wave packet and trajectory dynamics. Recently, these dynamics were studied by comparing the adiabatic representation and the exact factorization, with the purpose to highlight the different nature of topological-phase and geometric-phase effects arising in the two theoretical representations of the same problem. Here, we employ the exact factorization to understand how to accurately model low-energy dynamics in the vicinity of a conical intersection using an approximate description of the nuclear motion that uses trajectories. We find that since nonadiabatic effects are weak but non-negligible, the trajectory-based description that invokes the classical approximation struggles to capture the correct 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.

On the description of conical intersections between excited electronic states with LR-TDDFT and ADC(2)

Conical intersections constitute the conceptual bedrock of our working understanding of ultrafast, nonadiabatic processes within photochemistry (and photophysics). Accurate calculation of potential energy surfaces within the vicinity of conical intersections, however, still poses a serious challenge to many popular electronic structure methods. Multiple works have reported on the deficiency of methods like linear-response time-dependent density functional theory within the adiabatic approximation (AA LR-TDDFT) or algebraic diagrammatic construction to second-order [ADC(2)]-approaches often used in excited-state molecular dynamics simulations-to describe conical intersections between the ground and excited electronic states. In the present study, we focus our attention on conical intersections between excited electronic states and probe the ability of AA LR-TDDFT and ADC(2) to describe their topology and topography, using protonated formaldimine and pyrazine as two exemplar molecules. We also take the opportunity to revisit the performance of these methods in describing conical intersections involving the ground electronic state in protonated formaldimine-highlighting in particular how the intersection ring exhibited by AA LR-TDDFT can be perceived either as a (near-to-linear) seam of intersection or two interpenetrating cones, depending on the magnitude of molecular distortions within the branching space.

Application of molecular dynamic simulations in modeling the excited state behavior of confined molecules

Relative to isotropic organic solvent medium, the structure and conformation of a reactant molecule in an organized and confining medium are often different. In addition, because of the rigidity of the immediate environment, the reacting molecule have a little freedom to undergo large changes even upon gaining energy or modifications in the electronic structure. These alterations give rise to differences in the photochemistry of a molecular and supramolecular species. In this study, one such example is presented. α-Alkyl dibenzylketones upon excitation in isotropic solvents give products via Norrish type I and type II reactions that are independent of the chain length of the alkyl substituent. On the other hand, when these molecules are enclosed within an organic capsule of volume ~ 550 Å3, they give products that are strikingly dependent on the length of the α-alkyl substitution. These previously reported experimental observations are rationalized based on the structures generated by molecular modeling (docking and molecular dynamics (MD) simulations). It is shown that MD simulations that are utilized extensively in biologically important macromolecules can also be useful to understand the excited state behavior of reactive molecules that are part of supramolecular assemblies. These simulations can provide structural information of the reactant molecule and the surroundings complementing that with the one obtained from 1 and 2D NMR experiments. MD simulated structures of seven α-alkyl dibenzylketones encapsulated within the octa acid capsule provide a clear understanding of their unique behavior in this restricted medium. Because of the rigidity of the medium, these structures although generated in the ground state can rationalize the photochemical behavior of the molecules in the excited state.

Understanding the Impact of Symmetrical Substitution on the Photodynamics of Sinapate Esters Using Gas-Phase Ultrafast Spectroscopy

Two model biomimetic systems, ethyl sinapate (ES) and its symmetrical analogue, diethyl 2-(4-hydroxy-3,5-dimethoxybenzylidene)malonate (or diethyl sinapate, DES), are stripped to their core fundamentals through gas-phase spectroscopy to understand the underlying photophysics of photothermal materials. Following photoexcitation to the optically bright S1(ππ*) state, DES is found to repopulate the electronic ground state over 3 orders of magnitude quicker than its nonsymmetrical counterpart, ES. Our XMS-CASPT2 calculations shed light on the experimental results, revealing crucial differences in the potential energy surfaces and conical intersection topography between ES and DES. From this work, a peaked conical intersection, seen for DES, shows vital importance for the nonradiative ground-state recovery of photothermal materials. This fundamental comparative study highlights the potential impact that symmetrical substitution can have on the photodynamics of sinapate esters, providing a blueprint for future advancement in photothermal technology.

A visible-light-driven molecular motor based on barbituric acid

We present a class of visible-light-driven molecular motors based on barbituric acid. Due to a serendipitous reactivity we observed during their synthesis, these motors possess a tertiary stereogenic centre on the upper half, characterised by a hydroxy group. Using a combination of femto- and nanosecond transient absorption spectroscopy, molecular dynamics simulations and low-temperature 1H NMR experiments we found that these motors operate similarly to push-pull second-generation overcrowded alkene-based molecular motors. Interestingly, the hydroxy group at the stereocentre enables a hydrogen bond with the carbonyl groups of the barbituric acid lower half, which drives a sub-picosecond excited-state isomerisation, as observed spectroscopically. Computational simulations predict an excited state "lasso" mechanism where the intramolecular hydrogen bond pulls the molecule towards the formation of the metastable state, with a high predicted quantum yield of isomerisation (68%) in gas phase.

Serial Femtosecond Crystallography Reveals that Photoactivation in a Fluorescent Protein Proceeds via the Hula Twist Mechanism

Chromophore cis/trans photoisomerization is a fundamental process in chemistry and in the activation of many photosensitive proteins. A major task is understanding the effect of the protein environment on the efficiency and direction of this reaction compared to what is observed in the gas and solution phases. In this study, we set out to visualize the hula twist (HT) mechanism in a fluorescent protein, which is hypothesized to be the preferred mechanism in a spatially constrained binding pocket. We use a chlorine substituent to break the twofold symmetry of the embedded phenolic group of the chromophore and unambiguously identify the HT primary photoproduct. Through serial femtosecond crystallography, we then track the photoreaction from femtoseconds to the microsecond regime. We observe signals for the photoisomerization of the chromophore as early as 300 fs, obtaining the first experimental structural evidence of the HT mechanism in a protein on its femtosecond-to-picosecond timescale. We are then able to follow how chromophore isomerization and twisting lead to secondary structure rearrangements of the protein β-barrel across the time window of our measurements.

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.

Smooth Things Come in Threes: A Diabatic Surrogate Model for Conical Intersection Optimization

The optimization of conical intersection structures is complicated by the nondifferentiability of the adiabatic potential energy surfaces. In this work, we build a pseudodiabatic surrogate model, based on Gaussian process regression, formed by three smooth and differentiable surfaces that can adequately reproduce the adiabatic surfaces. Using this model with the restricted variance optimization method results in a notable decrease of the overall computational effort required to obtain minimum energy crossing points.