Do Photochemical Ring-Openings Occur in the Spectroscopic State? 1B2 Pathways for the Cyclohexadiene/Hexatriene Photochemical Interconversion

Revisiting fulgide photochromism: Mechanistic decoding and electron transport from computational exploration

The photochromic behavior of the fulgide molecule relies on ring-closure and ring-opening processes involving conical intersections during excited state transformation between isomers. The precise location and topography of these conical intersections significantly shape the decay process and fluorescence phenomena inherent to the molecule. This work combines electronic structure theory calculations using the density functional theory and wavefunction methods, as well as surface hopping simulation to analyze the photochemical behavior of an experimentally synthesized fulgide molecule, (E)-p-methylacetophenylisopropylidenesuccinic anhydride (1E). Our study reveals the conical intersection between the first excited state (S1) and the ground electronic state (S0), which emerges beyond the S1 minimum of 1E to the ring-closing side. The distinctive topography of this conical intersection appears to be sloped. These findings suggest a reduced quantum yield for the formation of the closed isomer, indicating a higher likelihood of reformation of the open isomer(s). The surface hopping simulation further supports this observation, revealing a mere ∼8% quantum yield for the formation of the closed isomer. In addition, the photoisomerization of the fulgide molecule initiates a cascade of conduction switching and holds great potential for applications in molecular electronics. Delving into the realm of molecular electronics, we have further examined the electron transport properties, disclosing the higher conductivity of the closed isomer.

Light-Induced Dyotropic Rearrangement of Diarylethenes: Scope, Mechanism, and Prospects

Irreversible two-photon photorearrangement of 1,2-diarylethenes is a unique process providing access to complex 2a1 ,5a-dihydro-5,6-dithiaacenaphthylene (DDA) heterocyclic core. This reaction was serendipitously discovered during studies on photoswitchable diarylethenes and was initially considered as a highly undesired process. However, in recent years, it has been recognized as an efficient photochemical reaction, interesting by itself and as a promising synthetic method for the synthesis of challenging molecules. Herein, we discuss the state-of-the-art in studies on this notable process.

Photochemical Ring-Opening Reaction of 1,3-Cyclohexadiene: Identifying the True Reactive State

The photochemically induced ring-opening isomerization reaction of 1,3-cyclohexadiene to 1,3,5-hexatriene is a textbook example of a pericyclic reaction and has been amply investigated with advanced spectroscopic techniques. The main open question has been the identification of the single reactive state which drives the process. The generally accepted description of the isomerization pathway starts with a valence excitation to the lowest lying bright state, followed by a passage through a conical intersection to the lowest lying doubly excited state, and finally a branching between either the return to the ground state of the cyclic molecule or the actual ring-opening reaction leading to the open-chain isomer. Here, in a joint experimental and computational effort, we demonstrate that the evolution of the excitation-deexcitation process is much more complex than that usually described. In particular, we show that an initially high-lying electronic state smoothly decreasing in energy along the reaction path plays a key role in the ring-opening reaction.

Photochemically Induced 1,3-Butadiene Ring-Closure from the Topological Analysis of the Electron Localization Function Viewpoint

The electronic rearrangement featuring the photochemically-induced 1,3-cis-butadiene is discussed within a bonding evolution theory (BET) perspective based on the topological analysis of the electron localization function and Thom's catastrophe theory. The process involves the vertical singlet-singlet excitation S₀ →S₂ , and the subsequent deactivation implying the S₂ /S₁ and S₁ /S₀ conical intersection regions. BET results reveal that the new CC bond is finally formed on the S₀ surface, as also recently found in the photochemical addition of two ethylenes [Phys. Chem. Chem. Phys. 23, 20598, (2021)].

Chiral photochemistry of achiral molecules

Chirality is a molecular property governed by the topography of the potential energy surface (PES). Thermally achiral molecules interconvert rapidly when the interconversion barrier between the two enantiomers is comparable to or lower than the thermal energy, in contrast to thermally stable chiral configurations. In principle, a change in the PES topography on the excited electronic state may diminish interconversion, leading to electronically prochiral molecules that can be converted from achiral to chiral by electronic excitation. Here we report that this is the case for two prototypical examples – cis-stilbene and cis-stiff stilbene. Both systems exhibit unidirectional photoisomerization for each enantiomer as a result of their electronic prochirality. We simulate an experiment to demonstrate this effect in cis-stilbene based on its interaction with circularly polarized light. Our results highlight the drastic change in chiral behavior upon electronic excitation, opening up the possibility for asymmetric photochemistry from an effectively nonchiral starting point.

The authors report non-adiabatic first principles molecular dynamics to show how an achiral molecule can be converted to a chiral one upon photoexcitation. These results demonstrate the possibility of asymmetric photochemistry starting from achiral reactants.

On Electron Pair Rearrangements in Photochemical Reactions: 1,3-Cyclohexadiene Ring Opening

1,3-Cyclohexadiene ring opening has been studied within the bonding evolution theory (BET) framework. We have focused on describing for the first time the electron pair rearrangements leading to the cis-1,3,5-hexatriene (HT) product from CHD. The nature of bonding in this process begins with the weakening of the double bonds in the Franck-Condon region. Along the 1¹B surface, the C-C sigma bond weakens. Meanwhile, its density redistributes toward the whole CHD ring, mainly over double bonds. Breaking of this bond occurs on the 2¹A surface due to the symmetrical splitting of pair density from this region. This density redistributes toward the reaction center once the pericyclic minimum is reached. The formation of the double bonds that characterize HT occurs gradually in the ground state. However, near the 2¹A/1¹A intersection, these bonds are partially established.

Structural or population dynamics: what is revealed by the time-resolved photoelectron spectroscopy of 1,3-cyclohexadiene? A study with an ensemble density functional theory method

Time-resolved photoelectron spectra during the photochemical ring-opening reaction of 1,3-cyclohexadiene (CHD) are modeled by an ensemble density functional theory (eDFT) method. The computational methodology employed in this work is capable of correctly describing the multi-reference effects arising in the ground and excited electronic states of molecules, which is important for the correct description of the ring-opening reaction of CHD. The geometries of molecular species along the non-adiabatic molecular dynamics (NAMD) trajectories reported in a previous study of the CHD photochemical ring-opening were used in this work to calculate the ionization energies and the respective Dyson orbitals for all possible ionization channels. The obtained theoretical time-resolved spectra display decay characteristics in a reasonable agreement with the experimental observations; i.e., the decay (and rise) of the most mechanistically significant signals occurs on the timescale of 100-150 fs. This is very different from the excited state population decay characteristics (τ_S1 = 234 ± 8 fs) obtained in the previous NAMD study. The difference between the population decay and the decay of the photoelectron signal intensity is traced back to the geometric transformation that the molecule undergoes during the photoreaction. This demonstrates the importance of including the geometric information in interpretation of the experimental observations.

Excited state dynamics of cis,cis-1,3-cyclooctadiene: Non-adiabatic trajectory surface hopping

We have performed trajectory surface hopping dynamics for cis,cis-1,3-cyclooctadiene to investigate the photochemical pathways involved after being excited to the S₁ state. Our calculations reveal ultrafast decay to the ground state, facilitated by conical intersections involving distortions around the double bonds. The main distortions are localized on one double bond, involving twisting and pyramidalization of one of the carbons of that double bond (similar to ethylene), while a limited number of trajectories decay via delocalized (non-local) twisting of both double bonds. The interplay between local and non-local distortions is important in our understanding of photoisomerization in conjugated systems. The calculations show that a broad range of the conical intersection seam space is accessed during the non-adiabatic events. Several products formed on the ground state have also been observed.

Ultrafast ring closing of a diarylethene-based photoswitchable nucleoside

Deoxyuridine nucleosides embodied into diarylethenes form an especial class of photoswitchable compounds that are designed to stack and pair with DNA bases. The molecular geometry can be switched between "open" and "closed" isomers by a pericyclic reaction that affects the stability of the surrounding double helix. This potentially enables light-induced control of DNA hybridization at microscopic resolution. Despite its importance for the optimization of DNA photoswitches, the ultrafast photoisomerization mechanism of these diarylethenes is still not well understood. In this work, femtosecond transient absorption spectroscopy is applied to study the ring closing reaction upon UV excitation with 45 fs pulses. Excited-state absorption decays rapidly and gives rise to the UV-Vis difference spectrum of the "closed" form within ≈15 ps. Time constants of 0.09, 0.49 and 6.6 ps characterize the multimodal dynamics, where a swift recurrence in the signal anisotropy indicates transient population of the intermediate 21A-like state.

Ultrafast Z → E photoisomerisation of structurally modified furylfulgides

Femtosecond transient absorption spectroscopy of Z-fulgides with selected structural motifs revealed fast and direct excited-state dynamics independent of chemical modifications.

Optimal control theory--closing the gap between theory and experiment

Optimal control theory and optimal control experiments are state-of-the-art tools to control quantum systems. Both methods have been demonstrated successfully for numerous applications in molecular physics, chemistry and biology. Modulated light pulses could be realized, driving these various control processes. Next to the control efficiency, a key issue is the understanding of the control mechanism. An obvious way is to seek support from theory. However, the underlying search strategies in theory and experiment towards the optimal laser field differ. While the optimal control theory operates in the time domain, optimal control experiments optimize the laser fields in the frequency domain. This also implies that both search procedures experience a different bias and follow different pathways on the search landscape. In this perspective we review our recent developments in optimal control theory and their applications. Especially, we focus on approaches, which close the gap between theory and experiment. To this extent we followed two ways. One uses sophisticated optimization algorithms, which enhance the capabilities of optimal control experiments. The other is to extend and modify the optimal control theory formalism in order to mimic the experimental conditions.

Unravelling the details of vitamin D photosynthesis by non-adiabatic molecular dynamics simulations

We investigate the photodynamics of vitamin D derivatives by a fully analytical implementation of the linear response time-dependent density functional theory surface hopping method (LR-TDDFT-SH). Our study elucidates the dynamics of the processes involved in vitamin D formation at the molecular level and with femtosecond resolution. We explain the major experimental findings and provide new insights that cannot directly be obtained from experiments: firstly, we investigate the dynamics of the photoinduced ring-opening of provitamin D (Pro) and cyclohexadiene (CHD) and the subsequent rotational isomerization. In agreement with recent experiments and CC2 calculations, only the bright S(1) state is involved in the ring-opening reaction. Our calculations confirm the experimentally reported 5 : 1 ratio between the excited state lifetimes of Pro and CHD. The longer lifetimes of Pro are attributed to steric constraints of the steroid skeleton and to temperature effects, both emerging directly from our simulations. For CHD and Pro, we present an explanation of the biexponential decay recently reported by Sension and coworkers [Tang et al., J. Phys. Chem., 2011, 134, 104503]: our calculations suggest that the fast and slow components arise from a reactive and an unreactive reaction pathway, respectively. Secondly, we assess the wavelength dependent photochemistry of previtamin D (Pre). Using replica exchange molecular dynamics we sample the Pre conformers present at thermal equilibrium. Based on this ensemble we explain the conformation dependent absorption and the essential features of Pre photochemistry. Consistent with the experiments, we find ring-closure to occur mostly after excitation of the cZc conformers and at lower energies, whereas Z/E isomerization of the central double bond preferably occurs after excitation at higher energies. For the isomerization we provide the first theoretical evidence of the proposed hula-twist mechanism. Our results show that LR-TDDFT-SH is a highly valuable tool for studying the photochemistry of moderately large systems, even though challenges remain in the vicinity of conical intersections.

The Ultrafast Pathway of Photon-Induced Electrocyclic Ring-Opening Reactions: The Case of 1,3-Cyclohexadiene

The photochemically induced electrocyclic ring-opening reaction of 1,3-cyclohexadiene to 1,3,5-hexatriene serves as a prototype for many important reactions in chemistry and in biological systems. Based on experimental and computational studies, a detailed picture of the reaction has now emerged: Excitation to the Franck-Condon region places the molecule on a steeply repulsive part of the 1B potential energy surface, which propels the molecule in exactly the conrotatory direction that conforms to the Woodward-Hoffmann rules of orbital symmetry. Bypassing a cusp in a symmetry-breaking direction, the wave packet enters the 2A state within 55 fs. It continues to move directly toward the 2A/1A conical intersection, where it crosses in approximately 80 fs to the ground state. This article summarizes the published experimental and theoretical work to describe the current understanding of the reaction while pointing to important questions that remain to be addressed.

Quantitative spatial distribution of sirolimus and polymers in drug-eluting stents using confocal Raman microscopy

Raman spectroscopy was used to differentiate each component found in the CYPHER Sirolimus-eluting Coronary Stent. The unique spectral features identified for each component were then used to develop three separate calibration curves to describe the solid phase distribution found on drug-polymer coated stents. The calibration curves were obtained by analyzing confocal Raman spectral depth profiles from a set of 16 unique formulations of drug-polymer coatings sprayed onto stents and planar substrates. The sirolimus model was linear from 0 to 100 wt % of drug. The individual polymer calibration curves for poly(ethylene-co-vinyl acetate) [PEVA] and poly(n-butyl methacrylate) [PBMA] were also linear from 0 to 100 wt %. The calibration curves were tested on three independent drug-polymer coated stents. The sirolimus calibration predicted the drug content within 1 wt % of the laboratory assay value. The polymer calibrations predicted the content within 7 wt % of the formulation solution content. Attenuated total reflectance Fourier transform infrared (ATR-FTIR) spectra from five formulations confirmed a linear response to changes in sirolimus and polymer content.