A Look Inside the Black Box of Machine Learning Photodynamics Simulations.

Photochemical reactions are of great importance in chemistry, biology, and materials science because they take advantage of a renewable energy source, mild reaction conditions, and high atom economy. Light absorption can excite molecules to a higher energy electronic state of the same spin multiplicity. The following nonadiabatic processes induce molecular transformations that afford exotic molecular architectures and high-energy-isomers that are inaccessible by thermal means. Computational simulations now complement time-resolved instrumentation to reveal ultrafast excited-state mechanistic information for photochemical reactions that is essential in disentangling elusive spectroscopic features, excited-state lifetimes, and excited-state mechanistic critical points. Nonadiabatic molecular dynamics (NAMD), powered by surface hopping techniques, is among the most widely applied techniques to model the photochemical reactions of medium-sized molecules. However, the computational efficiency is limited because of the requisite thousands of multiconfigurational quantum-chemical calculations multiplied by hundreds of trajectories. Machine learning (ML) has emerged as a revolutionary force in computational chemistry to predict the outcome of the resource-intensive multiconfigurational calculations on the fly. An ML potential trained with a substantial set of quantum-chemical calculations can predict the energies and forces with errors under chemical accuracy at a negligible cost. The integration of ML potentials in NAMD dramatically extends the maximum simulation time scale by ∼10 000-fold to the nanosecond regime.In this Account, we present a comprehensive demonstration of ML photodynamics simulations and summarize our most recent applications in resolving complex photochemical reactions. First, we address three fundamental components of ML techniques for photodynamics simulations: the quantum-chemical data set, the ML potential, and NAMD. Second, we describe best practices in building training data and our procedure toward training the ML photodynamics model with our recent literature contributions. We introduce a convenient training data generation scheme combining Wigner sampling and geometrical interpolation. It trains reliable and effective ML potentials suitable for subsequent active learning to detect undersampled data. We demonstrate how active learning automatically discovers new mechanistic pathways and reproduces experimental results. We point out that atomic permutation is an essential data augmentation approach to improve the learnability of distance-based molecular descriptors for highly symmetric molecules. Third, we demonstrate the utility of ML-photodynamics by showing the results of ML photodynamics simulations of (1) photo-torquoselective 4π disrotatory electrocyclic ring closing of norbornyl cyclohexadiene, which reveals a thermal conversion from experimentally unobserved intermediates to the reactant in 1 ns; (2) [2 + 2] photocycloaddition of substituted [3]-syn-ladderdienes in competition with 4π and 6π electrocyclic ring-opening reactions, uncovering substituent effects to explain the reported increased quantum yield of substituted cubane precursors; and (3) photochemical 4π disrotatory electrocyclic reactions of fluorobenzenes in nanoseconds with XMS-CASPT2-level training data. We expect this Account to broaden understanding of ML photodynamics and inspire future developments and applications to increasingly large molecules within complex environments on long time scales.

Multiconfigurational Calculations and Photodynamics Describe Norbornadiene Photochemistry.

Storing solar energy is a vital component of using renewable energy sources to meet the growing demands of the global energy economy. Molecular solar thermal (MOST) energy storage is a promising means to store solar energy with on-demand energy release. The light-induced isomerization reaction of norbornadiene (NBD) to quadricyclane (QC) is of great interest because of the generally high energy storage density (0.97 MJ kg-1) and long thermal reversion lifetime (t1/2,300K = 8346 years). However, the mechanistic details of the ultrafast excited-state [2 + 2]-cycloaddition are largely unknown due to the limitations of experimental techniques in resolving accurate excited-state molecular structures. We now present a full computational study on the excited-state deactivation mechanism of NBD and its dimethyl dicyano derivative (DMDCNBD) in the gas phase. Our multiconfigurational calculations and nonadiabatic molecular dynamics simulations have enumerated the possible pathways with 557 S2 trajectories of NBD for 500 fs and 492 S1 trajectories of DMDCNBD for 800 fs. The simulations predicted the S2 and S1 lifetimes of NBD (62 and 221 fs, respectively) and the S1 lifetime of DMDCNBD (190 fs). The predicted quantum yields of QC and DCQC are 10 and 43%, respectively. Our simulations also show the mechanisms of forming other possible reaction products and their quantum yields.

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.

Photochemically Mediated Polymerization of Molecular Furan and Pyridine: Synthesis of Nanothreads at Reduced Pressures.

Nanothreads are emerging one-dimensional sp3-hybridized materials with high predicted tensile strength and a tunable band gap. They can be synthesized by compressing aromatic or nonaromatic small molecules to pressures ranging from 15-30 GPa. Recently, new avenues are being sought that reduce the pressure required to afford nanothreads; the focus has been placed on the polymerization of molecules with reduced aromaticity, favorable stacking, and/or the use of higher reaction temperatures. Herein, we report the photochemically mediated polymerization of pyridine and furan aromatic precursors, which achieves nanothread formation at reduced pressures. In the case of pyridine, it was found that a combination of slow compression/decompression with broadband UV light exposure yielded a crystalline product featuring a six-fold diffraction pattern with similar interplanar spacings to previously synthesized pyridine-derived nanothreads at a reduced pressure. When furan is compressed to 8 GPa and exposed to broadband UV light, a crystalline solid is recovered that similarly demonstrates X-ray diffraction with an interplanar spacing akin to that of the high-pressure synthesized furan-derived nanothreads. Our method realizes a 1.9-fold reduction in the maximum pressure required to afford furan-derived nanothreads and a 1.4-fold reduction in pressure required for pyridine-derived nanothreads. Density functional theory and multiconfigurational wavefunction-based computations were used to understand the photochemical activation of furan and subsequent cascade thermal cycloadditions. The reduction of the onset pressure is caused by an initial [4+4] cycloaddition followed by increasingly facile thermal [4+2]-cycloadditions during polymerization.

Excited-State Distortions Promote the Photochemical 4π-Electrocyclizations of Fluorobenzenes via Machine Learning Accelerated Photodynamics Simulations.

Benzene fluorination increases chemoselectivities for Dewar-benzenes via 4π-disrotatory electrocyclization. However, the origin of the chemo- and regioselectivities of fluorobenzenes remains unexplained because of the experimental limitations in resolving the excited-state structures on ultrafast timescales. The computational cost of multiconfigurational nonadiabatic molecular dynamics simulations is also currently cost-prohibitive. We now provide high-fidelity structural information and reaction outcome predictions with machine-learning-accelerated photodynamics simulations of a series of fluorobenzenes, C6 F6-n Hn , n=0-3, to study their S1 →S0 decay in 4 ns. We trained neural networks with XMS-CASPT2(6,7)/aug-cc-pVDZ calculations, which reproduced the S1 absorption features with mean absolute errors of 0.04 eV (<2 nm). The predicted nonradiative decay constants for C6 F4 H2 , C6 F6 , C6 F3 H3 , and C6 F5 H are 116, 60, 28, and 12 ps, respectively, in broad qualitative agreement with the experiments. Our calculations show that a pseudo Jahn-Teller distortion of fluorinated benzenes leads to an S1 local-minimum region that extends the excited-state lifetimes of fluorobenzenes. The pseudo Jahn-Teller distortions reduce when fluorination decreases. Our analysis of the S1 dynamics shows that the pseudo-Jahn-Teller distortions promote an excited-state cis-trans isomerization of a πC-C bond. We characterized the surface hopping points from our NAMD simulations and identified instantaneous nuclear momentum as a factor that promotes the electrocyclizations.

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.

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.

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.

Role of the Perfluoro Effect in the Selective Photochemical Isomerization of Hexafluorobenzene.

Hexafluorobenzene and many of its derivatives exhibit a chemoselective photochemical isomerization, resulting in highly strained, Dewar-type bicyclohexenes. While the changes in absorption and emission associated with benzene hexafluorination have been attributed to the so-called "perfluoro effect", the resulting electronic structure and photochemical reactivity of hexafluorobenzene is still unclear. We now use a combination of ultrafast time-resolved spectroscopy, multiconfigurational computations, and non-adiabatic dynamics simulations to develop a holistic description of the absorption, emission, and photochemical dynamics of the 4π-electrocyclic ring-closing of hexafluorobenzene and the fluorination effect along the reaction coordinate. Our calculations suggest that the electron-withdrawing fluorine substituents induce a vibronic coupling between the lowest-energy 1B2u (ππ*) and 1E1g (πσ*) excited states by selectively stabilizing the σ-type states. The vibronic coupling occurs along vibrational modes of e2u symmetry which distorts the excited-state minimum geometry resulting in the experimentally broad, featureless absorption bands, and a ∼100 nm Stokes shift in fluorescence-in stark contrast to benzene. Finally, the vibronic coupling is shown to simultaneously destabilize the reaction pathway toward hexafluoro-benzvalene and promote molecular vibrations along the 4π ring-closing pathway, resulting in the chemoselectivity for hexafluoro-Dewar-benzene.

A Theoretical Stereoselectivity Model of Photochemical Denitrogenations of Diazoalkanes Toward Strained 1,3-Dihalogenated Bicyclobutanes.

Photochemical reactions exemplify "green" chemistry and are an essential tool for synthesizing highly strained molecules under mild conditions with light. The light-promoted denitrogenation of bicyclic azoalkanes affords functionalized, stereoenriched bicyclo[1.1.0]butanes. These reactions were revisited with multireference calculations and non-adiabatic molecular dynamics (NAMD) simulations to provide a detailed analysis of the photophysics, reactivities, and unexplained stereoselectivity of a series of diazabicyclo[2.1.1]hexenes. We used complete active space self-consistent field (CASSCF) calculations with an (8,8) active space and ANO-S-VDZP basis set; the CASSCF energies were corrected with CASPT2 (8,8)/ANO-S-VDZP. The nature of the electronic excitation is n → π* and ranges from 3.77 to 3.91 eV for the diazabicyclo[2.1.1]hexenes reported here. Minimum energy path calculations showed stepwise C-N bond breaking and led directly to a minimum energy crossing point, corresponding to a stereochemical "double inversion" product. Wigner sampling of diazabicyclo[2.1.1]hexene provided initial conditions for 692 NAMD trajectories. We identified competing complete stereoselective and stereochemical scrambling pathways. The stereoselective pathways feature concerted bicyclobutane inversion and N2 extrusion. The stereochemical scrambling pathways involve N2 extrusion followed by bicyclobutane planarization, leading to stereochemical scrambling. The predicted diastereomeric excess (d.e.) almost exactly matches the experiment (calc.d.e. = 46% vs exp.d.e. = 47%). Our NAMD simulations with 672, 568, and 596 trajectories for 1-F, 1-Cl, and 1-Br predicted a d.e. of 94-97% for the double inversion products. Halogenation significantly perturbs the potential energy surface (PES) toward the retention products due to hyperconjugative interactions. The nC → σ*C-X, X = F, Cl, Br hyperconjugative effect leads to a broader shoulder region on the PES for double inversion.

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.

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.

Multiconfigurational Calculations and Nonadiabatic Molecular Dynamics Explain Tricyclooctadiene Photochemical Chemoselectivity.

Sunlight is a renewable energy source that can be stored in chemical bonds using photochemical reactions. The synthesis of exotic and strained molecules is especially attractive with photochemical techniques because of the associated efficient and mild reaction conditions. We have understood the photophysics and subsequent photochemistry of a possible cubane precursor, tricyclo[4,2,0,02,5]octa-3,7-diene with complete active space self-consistent field (CASSCF) calculations with an (8,7) active space and the ANO-S-VDZP basis set to. The CASSCF energies were corrected with a second-order perturbative correction CASPT2(8,7)/ANO-S-VDZP. The S0 → S1 vertical excitation energy of 1 is 6.25 eV, which is a π → π* excitation. The minimum energy path from the S1 Franck-Condon point leads to a 4π-disrotatory electrocyclic ring-opening reaction to afford bicyclo[4,2,0]octa-2,4,7-triene. The 2D potential energy surface scan located a rhomboidal S1/S0 minimum energy crossing point connecting 1 and cubane, suggesting that a cycloaddition is theoretically possible. We used the fewest switches surface hopping to study the photodynamics of this cycloaddition: 85% of 1722 trajectories relaxed to eight products; the major products are bicyclo[4,2,0]octa-2,4,7-triene (30%) and cycloocta-1,3,5,7-tetraene (32%). Only 0.4% of trajectories undergo a [2 + 2] cycloaddition to form cubane.

Invasive mechanical ventilation using a bilevel PAP ST device in a healthy swine model.

PURPOSE: The Coronavirus Disease 2019 (COVID-19) pandemic may cause an acute shortage of ventilators. Standard noninvasive bilevel positive airway pressure devices with spontaneous and timed respirations (bilevel PAP ST) could provide invasive ventilation but evidence on their effectiveness in this capacity is limited. We sought to evaluate the ability of bilevel PAP ST to effect gas exchange via invasive ventilation in a healthy swine model. METHODS: Two single limb respiratory circuits with passive filtered exhalation were constructed and evaluated. Next, two bilevel PAP ST devices, designed for sleep laboratory and home use, were tested on an intubated healthy swine model using these circuits. These devices were compared to an anesthesia ventilator. RESULTS: We evaluated respiratory mechanics, minute ventilation, oxygenation, and presence of rebreathing for all of these devices. Both bilevel PAP ST devices were able to control the measured parameters. There were noted differences in performance between the two devices. Despite these differences, both devices provided effective invasive ventilation by controlling minute ventilation and providing adequate oxygenation in the animal model. CONCLUSIONS: Commercially available bilevel PAP ST can provide invasive ventilation with a single limb respiratory circuit and in-line filters to control oxygenation and ventilation without significant rebreathing in a swine model. Further study is needed to evaluate safety and efficacy in clinical disease models. In the setting of a ventilator shortage during the COVID-19 pandemic, and in other resource-constrained situations, these devices may be considered as an effective alternative means for invasive ventilation.

Cyclic Thiosulfinates and Cyclic Disulfides Selectively Cross-Link Thiols While Avoiding Modification of Lone Thiols.

This work addresses the need for chemical tools that can selectively form cross-links. Contemporary thiol-selective cross-linkers, for example, modify all accessible thiols, but only form cross-links between a subset. The resulting terminal "dead-end" modifications of lone thiols are toxic, confound cross-linking-based studies of macromolecular structure, and are an undesired, and currently unavoidable, byproduct in polymer synthesis. Using the thiol pair of Cu/Zn-superoxide dismutase (SOD1), we demonstrated that cyclic disulfides, including the drug/nutritional supplement lipoic acid, efficiently cross-linked thiol pairs but avoided dead-end modifications. Thiolate-directed nucleophilic attack upon the cyclic disulfide resulted in thiol-disulfide exchange and ring cleavage. The resulting disulfide-tethered terminal thiolate moiety either directed the reverse reaction, releasing the cyclic disulfide, or participated in oxidative disulfide (cross-link) formation. We hypothesized, and confirmed with density functional theory (DFT) calculations, that mono- S-oxo derivatives of cyclic disulfides formed a terminal sulfenic acid upon ring cleavage that obviated the previously rate-limiting step, thiol oxidation, and accelerated the new rate-determining step, ring cleavage. Our calculations suggest that the origin of accelerated ring cleavage is improved frontier molecular orbital overlap in the thiolate-disulfide interchange transition. Five- to seven-membered cyclic thiosulfinates were synthesized and efficiently cross-linked up to 104-fold faster than their cyclic disulfide precursors; functioned in the presence of biological concentrations of glutathione; and acted as cell-permeable, potent, tolerable, intracellular cross-linkers. This new class of thiol cross-linkers exhibited click-like attributes including, high yields driven by the enthalpies of disulfide and water formation, orthogonality with common functional groups, water-compatibility, and ring strain-dependence.

Nanosecond laser flash photolysis of a 6-nitroindolinospiropyran in solution and in nanocrystalline suspension under single excitation conditions.

Nanosecond transient absorption spectroscopy was used to study the photochemical ring-opening reaction for a 6-nitroindolinospiropyran (SP1) in solution and in nanocrystalline (NC) suspension at 298 K. We measured the kinetics in argon purged and air saturated acetonitrile and found that the presence of oxygen affected two out of the three components of the kinetic decay at 440 nm. These are assigned to the triplet excited states of the Z- and E-merocyanines (3Z-MC* and 3E-MC*). In contrast, a long-lived growth component at 550 nm and the decay of a band centered at 590 nm showed no dependence on oxygen and are assigned, respectively, to the ground state Z- and E-merocyanines (Z-MC0 and E-MC0). Laser flash photolysis studies performed in NC suspensions initially showed a very broad, featureless absorption spectrum that decayed uniformly for ca. 70 ns before revealing a more defined spectrum that persisted for greater than 4 ms and is consistent with a mixture of the more stable Z- and E-MC0 structures. We performed quantum mechanical calculations on the interconversion of E- and Z-MCs on the S0 and S1 potential energy surfaces. The computed UV-vis spectra for a scan along the Z → E interconversion reaction coordinate show substantial absorptivity from 300-600 nm, which suggests that the broad, featureless transient absorption spectrum results from the contribution of the transition structure and other high-energy species during the Z to E isomerization.

Understanding dispersive charge-transport in crystalline organic-semiconductors.

The effect of short-range order and dispersivity on charge-transport for organic crystalline semiconductors are important and unresolved questions. This exploration is the first to discern the role of short-range order on charge-transport for crystalline organic semiconductors. A multimode computational approach (including Molecular Dynamics and kinetic Monte Carlo simulations) is employed to understand the hole mobility dispersivity of crystalline organic semiconductors. Crystalline organic solids feature a mesoscale region where dispersive charge-transport dominates; our calculations show a clear transition of charge-mobility from non-dispersive to dispersive. An empirical relationship between the dispersive and non-dispersive transport transition region and ideal simulation box thickness is put forth. The dispersive to non-dispersive transition region occurs when energetic disorder approaches 72 meV. Non-dispersive transport is observed for simulation box sizes greater than 3.7 nm, which corresponds to approximately 12 π-stacked layers in typical π-stacked organic solids. A qualitative relationship is deduced between the variability of measured dispersive hole and variability of computed dispersive hole mobilities and system size. This relationship will guide future charge-transport investigations of condensed-phase organic systems.

1,3-Dipolar cycloaddition reactivities of perfluorinated aryl azides with enamines and strained dipolarophiles.

The reactivities of enamines and predistorted (strained) dipolarophiles toward perfluoroaryl azides (PFAAs) were explored experimentally and computationally. Kinetic analyses indicate that PFAAs undergo (3 + 2) cycloadditions with enamines up to 4 orders of magnitude faster than phenyl azide reacts with these dipolarophiles. DFT calculations were used to identify the origin of this rate acceleration. Orbital interactions between the cycloaddends are larger due to the relatively low-lying LUMO of PFAAs. The triazolines resulting from PFAA-enamine cycloadditions rearrange to amidines at room temperature, while (3 + 2) cycloadditions of enamines and phenyl azide yield stable, isolable triazolines. The 1,3-dipolar cycloadditions of norbornene and DIBAC also show increased reactivity toward PFAAs over phenyl azide but are slower than enamine-azide cycloadditions.

Mono-, Di-, and Trifluoroalkyl Substituent Effects on the Torquoselectivities of Cyclobutene and Oxetene Electrocyclic Ring Openings.

The reactivities and torquoselectivities of electrocyclic ring opening reactions of fluoromethyl-substituted cyclobutenes and oxetenes were studied with M06-2X density functional theory. The torquoselectivities of a series of mono-, di-, and trifluoromethylcyclobutenes and oxetenes result from the interplay of favorable orbital interactions and closed-shell repulsions. When the substituent rotates inward, there can be a favorable interaction between the breaking σ(CO) bond and the σ(CF)* orbital (σ(CO) → σ(CF)*) of the fluoromethyl group in fluoromethyloxetenes. The preference for rotation of a fluoromethyl group is decreased in trifluoromethyloxetenes because closed-shell repulsions between the breaking σ(CO) bond and trifluoromethyl substituent orbitals compete with the σ(CO) → σ(CF)* interaction.

Substituent effects on rates and torquoselectivities of electrocyclic ring-openings of N-substituted 2-azetines.

Transition structures for the conrotatory electrocyclic ring-opening reactions of N-substituted 2-azetines were computed with the density functional M06-2X/6-31+G(d,p). A wide range of substituents from π acceptors (e.g., CHO, CN) to π donors (NMe2, OMe) was explored. Acceptor substituents delocalize the nitrogen lone pair and stabilize the reactant state of 2-azetines, while donors destabilize the 2-azetine reactant state. The conrotatory ring-opening is torquoselective, and the transition state for the outward rotation of the N-substituent and inward rotation of the nitrogen lone pair is preferred. This transition structure is stabilized by an interaction between the nitrogen lone pair and the vacant π* orbital. The activation free energies are linearly related to the reaction free energies and the Taft σR(0) parameter.