Multiconfigurational photodynamics simulations reveal the mechanism of photodecarbonylations of cyclopropenones in explicit aqueous environments.

Gas-evolving photochemical reactions use light and mild conditions to access strained organic compounds irreversibly. Cyclopropenones are a class of light-responsive molecules used in bioorthogonal photoclick reactions; their excited-state decarbonylation reaction mechanisms are misunderstood due to their ultrafast (<100 femtosecond) lifetimes. We have combined multiconfigurational quantum mechanical (QM) calculations and non-adiabatic molecular dynamics (NAMD) simulations to uncover the excited-state mechanism of cyclopropenone and a photoprotected cyclooctyne-(COT)-precursor in gaseous and explicit aqueous environments. We explore the role of H-bonding with fully quantum mechanical explicitly solvated NAMD simulations for the decarbonylation reaction. The cyclopropenones pass through asynchronous conical intersections and have dynamically concerted photodecarbonylation mechanisms. The COT-precursor has a higher quantum yield of 55% than cyclopropenone (28%) because these trajectories prefer to break a σCC bond to avoid the strained trans-cyclooctene geometries. Our solvated simulations show an increased quantum yield (58%) for the systems studied here.

Design rules for optimization of photophysical and kinetic properties of azoarene photoswitches.

Azoarenes are an important class of molecular photoswitches that often undergo E → Z isomerization with ultraviolet light and have short Z-isomer lifetimes. Azobenzene has been a widely studied photoswitch for decades but can be poorly suited for photopharmacological applications due to its UV-light absorption and short-lived Z-isomer half-life (t1/2). Recently, diazo photoswitches with one or more thiophene rings in place of a phenyl ring have emerged as promising candidates, as they exhibit a stable photostationary state (98% E → Z conversion) and E-isomer absorption (λmax) in the visible light range (405 nm). In this work, we performed density functional theory calculations [PBE0-D3BJ/6-31+G(d,p)] on 26 hemi-azothiophenes, substituted with one phenyl ring and one thiophene ring on the diazo bond. We calculated the E-isomer absorption (λmax) and Z-isomer t1/2 for a set of 26 hemi-azothiophenes. We compared their properties to thiophene-based photoswitches that have been studied previously. We separated the 26 proposed photoswitches into four quadrants based on their λmax and t1/2 relative to past generations of hemi-azothiophene photoswitches. We note 8 hemi-azothiophenes with redshifted λmax and longer t1/2 than previous systems. Our top candidate has λmax and a t1/2 approaching 360 nm and 279 years, respectively. The results here present a pathway towards leveraging and optimizing two properties of photoswitches previously thought to be inversely related.

Cross-conjugation controls the stabilities and photophysical properties of heteroazoarene photoswitches.

Azoarene photoswitches are versatile molecules that interconvert from their E-isomer to their Z-isomer with light. Azobenzene is a prototypical photoswitch but its derivatives can be poorly suited for in vivo applications such as photopharmacology due to undesired photochemical reactions promoted by ultraviolet light and the relatively short half-life (t1/2) of the Z-isomer (2 days). Experimental and computational studies suggest that these properties (λmax of the E isomer and t1/2 of the Z-isomer) are inversely related. We identified isomeric azobisthiophenes and azobisfurans from a high-throughput screening study of 1540 azoarenes as photoswitch candidates with improved λmax and t1/2 values relative to azobenzene. We used density functional theory to predict the activation free energies and vertical excitation energies of the E- and Z-isomers of 2,2- and 3,3-substituted azobisthiophenes and azobisfurans. The half-lives depend on whether the heterocycles are π-conjugated or cross-conjugated with the diazo π-bond. The 2,2-substituted azoarenes both have t1/2 values on the scale of 1 hour, while the 3,3-analogues have computed half-lives of 40 and 230 years (thiophene and furan, respectively). The 2,2-substituted heteroazoarenes have significantly higher λmax absorptions than their 3,3-substituted analogues: 76 nm for azofuran and 77 nm for azothiophene.

A Green Chemistry Approach toward the Stereospecific Synthesis of Densely Functionalized Cyclopropanes via the Solid-State Photodenitrogenation of Crystalline 1-Pyrazolines.

The cyclopropane ring features prominently in active pharmaceuticals, and this has spurred the development of synthetic methodologies that effectively incorporate this highly strained motif into such molecules. As such, elegant solutions to prepare densely functionalized cyclopropanes, particularly ones embedded within the core of complex structures, have become increasingly sought-after. Here we report the stereospecific synthesis of a set of cyclopropanes with vicinal quaternary stereocenters via the solvent-free solid-state photodenitrogenation of crystalline 1-pyrazolines. Density functional theory calculations at the M062X/6-31+G(d,p) level of theory were used to determine the origin of regioselectivity for the synthesis of the 1-pyrazolines; favorable in-phase frontier molecular orbital interactions are responsible for the observation of a single pyrazoline regioisomer. It was also shown that the loss of N2 may take place via a highly selective solid-state thermal reaction. Scalability of the solid-state photoreaction is enabled through aqueous nanocrystalline suspensions, making this method a "greener" alternative to effectively facilitate the construction of cyclopropane-containing molecular scaffolds.

Efficient Discovery of Visible Light-Activated Azoarene Photoswitches with Long Half-Lives Using Active Search.

Photoswitches are molecules that undergo a reversible, structural isomerization after exposure to certain wavelengths of light. The dynamic control offered by molecular photoswitches is favorable for materials chemistry, photopharmacology, and catalysis applications. Ideal photoswitches absorb visible light and have long-lived metastable isomers. We used high-throughput virtual screening to predict the absorption maxima (λmax) of the E-isomer and half-life (t1/2) of the Z-isomer. However, computing the photophysical and kinetic stabilities with density functional theory of each entry of a virtual molecular library containing thousands or millions of molecules is prohibitively time-consuming. We applied active search, a machine-learning technique, to intelligently search a chemical search space of 255 991 photoswitches based on 29 known azoarenes and their derivatives. We iteratively trained the active search algorithm on whether a candidate absorbed visible light (λmax > 450 nm). Active search was found to triple the discovery rate compared to random search. Further, we projected 1962 photoswitches to 2D using the Uniform Manifold Approximation and Projection algorithm and found that λmax depends on the core, which is tunable by substituents. We then incorporated a second stage of screening to predict the stabilities of the Z-isomers for the top candidates of each core. We identified four ideal photoswitches that concurrently satisfy the following criteria: λmax > 450 nm and t1/2 > 2 h.These candidates had λmax and t1/2 range from 465 to 531 nm and hours to days, respectively.

Machine-Learning Photodynamics Simulations Uncover the Role of Substituent Effects on the Photochemical Formation of Cubanes.

Photochemical [2 + 2]-cycloadditions store solar energy in chemical bonds and efficiently access strained organic molecular architectures. Functionalized [3]-ladderdienes undergo [2 + 2]-photocycloadditions to afford cubanes, a class of strained organic molecules. The substituents (e.g., methyl, trifluoromethyl, and cyclopropyl) affect the overall reactivities of these cubane precursors; the yields range from 1 to 48%. However, the origin of these substituent effects on the reactivities and chemoselectivities is not understood. We now integrate single and multireference calculations and machine-learning-accelerated nonadiabatic molecular dynamics (ML-NAMD) to understand how substituents affect the ultrafast dynamics and mechanism of [2 + 2]-photocycloadditions. Steric clashes between substituent groups destabilize the 4π-electrocyclic ring-opening pathway and minimum energy conical intersections by 0.72-1.15 eV and reaction energies by 0.68-2.34 eV. Noncovalent dispersive interactions stabilize the [2 + 2]-photocycloaddition pathway; the conical intersection energies are lower by 0.31-0.85 eV, and the reaction energies are lower by 0.03-0.82 eV. The 2 ps ML-NAMD trajectories reveal that closed-shell repulsions block a 6π-conrotatory electrocyclic ring-opening pathway with increasing steric bulk. Thirty-eight percent of the methyl-substituted [3]-ladderdiene trajectories proceed through the 6π-conrotatory electrocyclic ring-opening, whereas the trifluoromethyl- and cyclopropyl-substituted [3]-ladderdienes prefer the [2 + 2]-photocycloaddition pathways. The predicted cubane yields (H: 0.4% < CH3: 1% < CF3: 14% < cPr: 15%) match the experimental trend; these substituents predistort the reactants to resemble the conical intersection geometries leading to cubanes.

Benchmarking of Density Functionals for Z-Azoarene Half-Lives via Automated Transition State Search.

Molecular photoswitches use light to interconvert from a thermodynamically stable isomer into a metastable isomer. Photoswitches have been used in photopharmacology, catalysis, and molecular solar thermal (MOST) materials because of their spatiotemporal activation. Visible-light-absorbing photoswitches are especially attractive because low-energy light minimizes undesired photochemical reactions and enables biological applications. Ideal photoswitches require well-separated absorption spectra for both isomers and long-lived metastable states. However, predicting thermal half-lives with density functional theory is difficult because it requires locating transition structures and chosing an accurate model chemistry. We now report EZ-TS; by automatically calculating activation energies for the thermal Z → E isomerization. We used 28 density functionals [local spin density approximation, generalized gradient approximation, meta-GGA, hybrid GGA, and hybrid meta-GGA] and five basis sets [6-31G(d), 6-31+G(d,p), 6-311+G(d,p), cc-pVDZ, and aug-cc-pVDZ]. The hybrid GGA functionals performed the best among all tested functionals. We demonstrate that the mean absolute errors of 14 model chemistries approach chemical accuracy.