Bicyclo[2.2.0]hexene: A Multicyclic Mechanophore with Reactivity Diversified by External Forces.

Polymer mechanochemistry has been established as an enabling tool in accessing chemical reactivity and reaction pathways that are distinctive from their thermal counterparts. However, eliciting diversified reaction pathways by activating different constituent chemical bonds from the same mechanophore structure remains challenging. Here, we report the design of a bicyclo[2.2.0]hexene (BCH) mechanophore to leverage its structural simplicity and relatively low molecular symmetry to demonstrate this idea of multimodal activation. Upon changing the attachment points of pendant polymer chains, three different C-C bonds in bicyclo[2.2.0]hexene are specifically activated via externally applied force by sonication. Experimental characterization confirms that in different scenarios of polymer attachment, the regioisomers of BCH undergo different activation reactions, entailing retro-[2+2] cycloreversion, 1,3-allylic migration, and retro-4π ring-opening reactions, respectively. Control experiments with small-molecule analogues reveal that the observed diversified reactivity of BCH regioisomers is possible only with mechanical force. Theoretical studies further elucidate that the differences in the positions of substitution between regioisomers have a minimal impact on the potential energy surface of the parent BCH scaffold. The mechanochemical selectivity between different C-C bonds in each constitutional isomer is a result of selective and effective coupling of force to the aligned C-C bond in each case.

Mechanochemistry of Pterodactylane.

Pterodactylane is a [4]-ladderane with substituents on the central rung. Comparing the mechanochemistry of the [4]-ladderane structure when pulled from the central rung versus the end rung revealed a striking difference in the threshold force of mechanoactivation: the threshold force is dramatically lowered from 1.9 nN when pulled on the end rung to 0.7 nN when pulled on the central rung. We investigated the bicyclic products formed from the mechanochemical activation of pterodactylane experimentally and computationally, which are distinct from the mechanochemical products of ladderanes being activated from the end rung. We compared the products of pterodactylane's mechanochemical and thermal activation to reveal differences and similarities in the mechanochemical and thermal pathways of pterodactylane transformation. Interestingly, we also discovered the presence of elementary steps that are accelerated or suppressed by force within the same mechanochemical reaction of pterodactylane, suggesting rich mechanochemical manifolds of multicyclic structures. We rationalized the greatly enhanced mechanochemical reactivity of the central rung of pterodactylane and discovered force-free ground state bond length to be a good low-cost predictor of the threshold force for cyclobutane-based mechanophores. These findings advance our understanding of mechanochemical reactivities and pathways, and they will guide future designs of mechanophores with low threshold forces to facilitate their applications in force-responsive materials.

Efficient Acceleration of Reaction Discovery in the Ab Initio Nanoreactor: Phenyl Radical Oxidation Chemistry.

Over the years, many computational strategies have been employed to elucidate reaction networks. One of these methods is accelerated molecular dynamics, which can circumvent the expense required in dynamics to find all reactants and products (local minima) and transition states (first-order saddle points) on a potential energy surface (PES) by using fictitious forces that promote reaction events. The ab initio nanoreactor uses these accelerating forces to study large chemical reaction networks from first-principles quantum mechanics. In the initial nanoreactor studies, this acceleration was done through a piston periodic compression potential, which pushes molecules together to induce entropically unfavorable bimolecular reactions. However, the piston is not effective for discovering intramolecular and dissociative reactions, such as those integral to the decomposition channels of phenyl radical oxidation. In fact, the choice of accelerating forces dictates not only the rate of reaction discovery but also the types of reactions discovered; thus, it is critical to understand the biases and efficacies of these forces. In this study, we examine forces using metadynamics, attractive potentials, and local thermostats for accelerating reaction discovery. For each force, we construct a separate phenyl radical combustion reaction network using solely that force in discovery trajectories. We elucidate the enthalpic and entropic trends of each accelerating force and highlight their efficiency in reaction discovery. Comparing the nanoreactor-constructed reaction networks with literature renditions of the phenyl radical combustion PES shows that a combination of accelerating forces is best suited for reaction discovery.

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis.

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

3-Ligated Phenothiazinyl-terephthalonitrile Dyads and Triads - Synthesis, Electronic Properties, Delayed Fluorescence and Electronic Structure.

The bromine-lithium exchange-borylation-Suzuki sequence efficiently furnishes phenothiazine-terephthalonitrile donor-acceptor dyads and triads in high yields. In contrast to most phenothiazine-acceptor conjugates the title compounds are ligated in p-position to the phenothiazine nitrogen atom. Moreover, the acceptors are either directly linked or ligated by an arylene bridge and p-anisyl N-substituents on the phenothiazine are chosen to lock the tricycle into an intra-configuration. Cyclic voltammetry reveals effects of bridging and ligation of the N-substituent. Optical spectroscopy likewise displays similar band gaps, large Stokes shifts and substantial to high quantum yields in solution, in the solid state and in PMMA matrix. Time-resolved fluorescence spectroscopy indicates quite long fluorescence decay times in solution and emission components in the microsecond time range. TADF properties are further assessed by fluorescence increase in deoxygenated solution, gated emission spectroscopy and temperature-dependent determination of phosphorescence. The nature of the electronically excited states is investigated by DFT/MRCI. While for the directly ligated dyad a singlet-triplet energy gap Δ E ( S 1 - T 1 ) ${{E}_{{({\rm S}}_{1}-{{\rm T}}_{1})}{\rm \ }}$ of 0.24 eV can be estimated and is consistently confirmed by quantum chemical calculations on the lowest energy conformer, even lower Δ E ( S 1 - T 1 ) ${{\rm \Delta }{E}_{{(S}_{1}-{T}_{1})}{\rm \ }}$ of 0.029 and 0.008 eV are estimated for the investigated dyads and the triad in the solid state and in PMMA matrix.

Mechanoresponsive Metal-Organic Cage-Crosslinked Polymer Hydrogels.

We report the formation of metal-organic cage-crosslinked polymer hydrogels. To enable crosslinking of the cages and subsequent network formation, we used homodifunctionalized poly(ethylene glycol) (PEG) chains terminally substituted with bipyridines as ligands for the Pd6 L4 corners. The encapsulation of guest molecules into supramolecular self-assembled metal-organic cage-crosslinked hydrogels, as well as ultrasound-induced disassembly of the cages with release of their cargo, is presented in addition to their characterization by nuclear magnetic resonance (NMR) techniques, rheology, and comprehensive small-angle X-ray scattering (SAXS) experiments. The constrained geometries simulating external force (CoGEF) method and barriers using a force-modified potential energy surface (FMPES) suggest that the cage-opening mechanism starts with the dissociation of one pyridine ligand at around 0.5 nN. We show the efficient sonochemical activation of the hydrogels HG3 -6 , increasing the non-covalent guest-loading of completely unmodified drugs available for release by a factor of ten in comparison to non-crosslinked, star-shaped assemblies in solution.

A multi-stage single photochrome system for controlled photoswitching responses.

The ability of molecular photoswitches to convert on/off responses into large macroscale property change is fundamental to light-responsive materials. However, moving beyond simple binary responses necessitates the introduction of new elements that control the chemistry of the photoswitching process at the molecular scale. To achieve this goal, we designed, synthesized and developed a single photochrome, based on a modified donor-acceptor Stenhouse adduct (DASA), capable of independently addressing multiple molecular states. The multi-stage photoswitch enables complex switching phenomena. To demonstrate this, we show spatial control of the transformation of a three-stage photoswitch by tuning the population of intermediates along the multi-step reaction pathway of the DASAs without interfering with either the first or final stage. This allows for a photonic three-stage logic gate where the secondary wavelength solely negates the input of the primary wavelength. These results provide a new strategy to move beyond traditional on/off binary photochromic systems and enable the design of future molecular logic systems.

Proton Transfer from a Photoacid to a Water Wire: First Principles Simulations and Fast Fluorescence Spectroscopy.

Proton transfer reactions are ubiquitous in chemistry, especially in aqueous solutions. We investigate photoinduced proton transfer between the photoacid 8-hydroxypyrene-1,3,6-trisulfonate (HPTS) and water using fast fluorescence spectroscopy and ab initio molecular dynamics simulations. Photoexcitation causes rapid proton release from the HPTS hydroxyl. Previous experiments on HPTS/water described the progress from photoexcitation to proton diffusion using kinetic equations with two time constants. The shortest time constant has been interpreted as protonated and photoexcited HPTS evolving into an "associated" state, where the proton is "shared" between the HPTS hydroxyl and an originally hydrogen bonded water. The longer time constant has been interpreted as indicating evolution to a "solvent separated" state where the shared proton undergoes long distance diffusion. In this work, we refine the previous experimental results using very pure HPTS. We then use excited state ab initio molecular dynamics to elucidate the detailed molecular mechanism of aqueous excited state proton transfer in HPTS. We find that the initial excitation results in rapid rearrangement of water, forming a strong hydrogen bonded network (a "water wire") around HPTS. HPTS then deprotonates in ≤3 ps, resulting in a proton that migrates back and forth along the wire before localizing on a single water molecule. We find a near linear relationship between the emission wavelength and proton-HPTS distance over the simulated time scale, suggesting that the emission wavelength can be used as a ruler for the proton distance. Our simulations reveal that the "associated" state corresponds to a water wire with a mobile proton and that the diffusion of the proton away from this water wire (to a generalized "solvent-separated" state) corresponds to the longest experimental time constant.

Flyby reaction trajectories: Chemical dynamics under extrinsic force.

Dynamic effects are an important determinant of chemical reactivity and selectivity, but the deliberate manipulation of atomic motions during a chemical transformation is not straightforward. Here, we demonstrate that extrinsic force exerted upon cyclobutanes by stretching pendant polymer chains influences product selectivity through force-imparted nonstatistical dynamic effects on the stepwise ring-opening reaction. The high product stereoselectivity is quantified by carbon-13 labeling and shown to depend on external force, reactant stereochemistry, and intermediate stability. Computational modeling and simulations show that, besides altering energy barriers, the mechanical force activates reactive intramolecular motions nonstatistically, setting up "flyby trajectories" that advance directly to product without isomerization excursions. A mechanistic model incorporating nonstatistical dynamic effects accounts for isomer-dependent mechanochemical stereoselectivity.

Understanding the Mechanochemistry of Ladder-Type Cyclobutane Mechanophores by Single Molecule Force Spectroscopy.

We have recently reported a series of ladder-type cyclobutane mechanophores, polymers of which can transform from nonconjugated structures to conjugated structures and change many properties at once. These multicyclic mechanophores, namely, exo-ladderane/ene, endo-benzoladderene, and exo-bicyclohexene-peri-naphthalene, have different ring structures fused to the first cyclobutane, significantly different free energy changes for ring-opening, and different stereochemistry. To better understand their mechanochemistry, we used single molecule force spectroscopy (SMFS) to characterize their force-extension behavior and measure the threshold forces. The threshold forces correlate with the activation energy of the first bond, but not with the strain of the fused rings distal to the polymer main chain, suggesting that the activation of these ladder-type mechanophores occurs with similar early transition states, which is supported by force-modified potential energy surface calculations. We further determined the stereochemistry of the mechanically generated dienes and observed significant and variable contour length elongation for these mechanophores both experimentally and computationally. The fundamental understanding of ladder-type mechanophores will facilitate future design of multicyclic mechanophores with amplified force-response and their applications as mechanically responsive materials.

Electrostatic Control of Photoisomerization in Channelrhodopsin 2.

Channelrhodopsin 2 (ChR2) is the most commonly used tool in optogenetics. Because of its faster photocycle compared to wild-type (WT) ChR2, the E123T mutant of ChR2 is a useful optogenetic tool when fast neuronal stimulation is needed. Interestingly, in spite of its faster photocycle, the initial step of the photocycle in E123T (photoisomerization of retinal protonated Schiff base or RPSB) was found experimentally to be much slower than that of WT ChR2. The E123T mutant replaces the negatively charged E123 residue with a neutral T123 residue, perturbing the electric field around the RPSB. Understanding the RPSB photoisomerization mechanism in ChR2 mutants will provide molecular-level insights into how ChR2 photochemical reactivity can be controlled, which will lay the foundation for improving the design of optogenetic tools. In this work, we combine ab initio nonadiabatic dynamics simulation, excited state free energy calculation, and reaction path search to comprehensively characterize the RPSB photoisomerization mechanism in the E123T mutant of ChR2. Our simulation agrees with previous experiments in predicting a red-shifted absorption spectrum and significant slowdown of photoisomerization in the E123T mutant. Interestingly, our simulations predict similar photoisomerization quantum yields for the mutant and WT despite the differences in excited-state lifetime and absorption maximum. Upon mutation, the neutralization of the negative charge on the E123 residue increases the isomerization barrier, alters the reaction pathway, and changes the relative stability of two fluorescent states. Our findings provide new insight into the intricate role of the electrostatic environment on the RPSB photoisomerization mechanism in microbial rhodopsins.

Substituent Effects in Mechanochemical Allowed and Forbidden Cyclobutene Ring-Opening Reactions.

Woodward and Hoffman once jested that a very powerful Maxwell demon could seize a molecule of cyclobutene at its methylene groups and tear it open in a disrotatory fashion to obtain butadiene (Woodward, R. B.; Hoffmann, R. The Conservation of Orbital Symmetry. Angew. Chem., Int. Ed. 1969, 8, 781-853). Nearly 40 years later, that demon was discovered, and the field of covalent polymer mechanochemistry was born. In the decade since our demon was befriended, many fundamental investigations have been undertaken to build up our understanding of force-modified pathways for electrocyclic ring-opening reactions. Here, we seek to extend that fundamental understanding by exploring substituent effects in allowed and forbidden ring-opening reactions of cyclobutene (CBE) and benzocyclobutene (BCB) using a combination of single-molecule force spectroscopy (SMFS) and computation. We show that, while the forbidden ring-opening of cis-BCB occurs at a lower force than the allowed ring-opening of trans-BCB on the time scale of the SMFS experiment, the opposite is true for cis- and trans-CBE. Such a reactivity flip is explained through computational analysis and discussion of the so-called allowed/forbidden gap.

Proton Transfer Dynamics in the Aprotic Proton Accepting Solvent 1-Methylimidazole.

The dynamics of proton transfer to the aprotic solvent 1-methylimidazole (MeIm, proton acceptor) from the photoacid 8-hydroxypyrene-1,3,6-trisulfonic acid trisodium salt (HPTS) was investigated using fast fluorescence measurements. The closely related molecule, 8-methoxypyrene-1,3,6-trisulfonic acid trisodium salt (MPTS), which is not a photoacid, was also studied for comparison. Following optical excitation, the wavelength-dependent population dynamics of HPTS in MeIm resulting from the deprotonation process were collected over the entire fluorescence emission window. Analysis of the time-dependent fluorescence spectra revealed four distinct fluorescence bands that appear and decay on different time scales. We label these four states as protonated (P), associated I (AI), associated II (AII), and deprotonated (D). We find that the simple kinetic scheme of P → AI → AII → D is not consistent with the data. Instead, the kinetic scheme that describes the data has P decaying into AI, which mainly goes on to deprotonation (D), but AI can also feed into AII. AII can return to AI or decay to the ground state, but does not deprotonate within experimental error. Quantum chemistry and excited state QM/MM Born-Oppenheimer molecular dynamics simulations indicate that AI and AII are two H-bonding conformations of MeIm to the HPTS hydroxyl, axial, and equatorial, respectively.

Vibrational analysis of methyl cation-Rare gas atom complexes: CH3 +-Rg (Rg = He, Ne, Ar, Kr).

The vibrational spectra of simple CH3 +-Rg (Rg = He, Ne, Ar, Kr) complexes have been studied by vibrational configuration interaction theory relying on multidimensional potential energy surfaces (PESs) obtained from explicitly correlated coupled cluster calculations, CCSD(T)-F12a. In agreement with experimental results, the series of rare gas atoms leads to rather unsystematic results and indicates huge zero point vibrational energy effects for the helium complex. In order to study these sensitive complexes more consistently, we also introduce configuration averaged vibrational self-consistent field theory, which is a generalization of standard vibrational self-consistent field theory to several configurations. The vibrational spectra of the complexes are compared to that of the methyl cation, for which corrections due to scalar-relativistic effects, high-order coupled-cluster terms, e.g., quadruple excitations, and core-valence correlation have explicitly been accounted for. The occurrence of tunneling splittings for the vibrational ground-state of CH3 +-He has been investigated on the basis of semiclassical instanton theory. These calculations and a direct comparison of the energy profiles along the intrinsic reaction coordinates with that of the hydronium cation, H3O+, suggest that tunneling effects for vibrationally excited states should be very small.

Asymmetric Carboxycyanation of Aldehydes by Cooperative AlF/Onium Salt Catalysts: from Cyanoformate to KCN as Cyanide Source.

Asymmetric 1,2-additions of cyanide yield enantioenriched cyanohydrins as versatile chiral building blocks. Next to HCN, volatile organic cyanide sources are usually used. Among them, cyanoformates are more attractive on technical scale than TMSCN for cost reasons, but catalytic productivity is usually lower. Here, the development of a new strategy for cyanations is described, in which this activity disadvantage is overcome. A Lewis acidic Al center cooperates with an aprotic onium moiety within a remarkably robust bifunctional Al-F-salen complex. This allowed for unprecedented turnover numbers of up to 104 . DFT studies suggest an unexpected unique trimolecular pathway in which the ammonium bound cyanide attacks the aldehyde, which itself is activated by the carbonyl group of the cyanoformate binding to the Al center. In addition, a novel practical carboxycyanation method was developed that makes use of KCN as the sole cyanide source. The use of a pyrocarbonate as carboxylating reagent provided the best results.

Dual-Level Approach to Instanton Theory.

Instanton theory is an established method to calculate rate constants of chemical reactions including atom tunneling. Technical and methodological improvements increased its applicability. Still, a large number of energy and gradient calculations is necessary to optimize the instanton tunneling path, and second derivatives of the potential energy along the tunneling path have to be evaluated, restricting the range of suitable electronic structure methods. To enhance the applicability of instanton theory, we present a dual-level approach in which instanton optimizations and Hessian calculations are performed using an efficient but approximate electronic structure method, and the potential energy along the tunneling path is recalculated using a more accurate method. This procedure extends the applicability of instanton theory to high-level electronic structure methods for which analytic gradients may not be available, like local linear-scaling approaches. We demonstrate for the analytical Eckart barrier and three molecular systems how the dual-level instanton approach corrects for the largest part of the error caused by the inaccuracy of the efficient electronic structure method. This reduces the error of the calculated rate constants significantly.

The Lewis Pair Polymerization of Lactones Using Metal Halides and N-Heterocyclic Olefins: Theoretical Insights.

Lewis pair polymerization employing N-Heterocyclic olefins (NHOs) and simple metal halides as co-catalysts has emerged as a useful tool to polymerize diverse lactones. To elucidate some of the mechanistic aspects that remain unclear to date and to better understand the impact of the metal species, computational methods have been applied. Several key aspects have been considered: (1) the formation of NHO-metal halide adducts has been evaluated for eight different NHOs and three different Lewis acids, (2) the coordination of four lactones to MgCl2 was studied and (3) the deprotonation of an initiator (butanol) was investigated in the presence and absence of metal halide for one specific Lewis pair. It was found that the propensity for adduct formation can be influenced, perhaps even designed, by varying both organic and metallic components. Apart from the NHO backbone, the substituents on the exocyclic, olefinic carbon have emerged as interesting tuning site. The tendency to form adducts is ZnCl2 > MgCl2 > LiCl. If lactones coordinate to MgCl2, the most likely binding mode is via the carbonyl oxygen. A chelating coordination cannot be ruled out and seems to gain importance upon increasing ring-size of the lactone. For a representative NHO, it is demonstrated that in a metal-free setting an initiating alcohol cannot be deprotonated, while in the presence of MgCl2 the same process is exothermic with a low barrier.

Double Regioselective Asymmetric C-Allylation of Isoxazolinones: Iridium-Catalyzed N-Allylation Followed by an Aza-Cope Rearrangement.

Isoxazolinones are biologically and synthetically interesting densely functionalized heterocycles, which for a long time were not accessible in enantioenriched form by asymmetric catalysis. Next to the deficit of enantioselective methods, the functionalization of isoxazolinones is often plagued by regioselectivity issues due to the competition of various nucleophilic centers within the heterocycles. We report the first regio- and enantioselective C-allylations of isoxazolinones. These occur with high regioselectivity in favor of the linear allylation products, although Ir phosphoramidite catalysts were used, which commonly results in branched isomers. Our studies suggest that this outcome is the result of a reaction cascade via an initial regio- and enantioselective N-allylation to provide a branched allyl intermediate, followed by a spontaneous [3,3]-rearrangement resulting in chirality transfer.

Support for a Dioxyallyl Cation in the Mechanism Leading to (-)-Levoglucosenone.

Levoglucosenone (LGO) is the major product formed when cellulose is pyrolyzed in the presence of acid at temperatures between 170 and 350 °C. The current intense interest in biomass conversion has led to a number of reports on its preparation; however, there is still uncertainty on the mechanism leading to LGO. We propose a new mechanism which involves a C2-C1 hydride shift followed by intramolecular trapping of a dioxyallyl cation. The reaction has been modeled using DFT calculations from the known LGO precursors levoglucosan and 1,4:3,6-dianhydro-α-D-glucopyranose to a common intermediate with calculated barriers of 10.6 and 13.5 kcal·mol-1, respectively. A discussion of the literature on the formation of LGO from late pathway intermediates is also provided.