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.

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.

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.

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.

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.

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.

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.

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.

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.

Asymmetric Ketone Reduction by Imine Reductases.

The rapidly growing area of asymmetric imine reduction by imine reductases (IREDs) has provided alternative routes to chiral amines. Here we report the expansion of the reaction scope of IREDs by showing the stereoselective reduction of 2,2,2-trifluoroacetophenone. Assisted by an in silico analysis of energy barriers, we evaluated asymmetric hydrogenations of carbonyls and imines while considering the influence of substrate reactivity on the chemoselectivity of this novel class of reductases. We report the asymmetric reduction of C=N as well as C=O bonds catalysed by members of the IRED enzyme family.

High Oxidation State Molybdenum N-Heterocyclic Carbene Alkylidyne Complexes: Synthesis, Mechanistic Studies, and Reactivity.

The first synthetic protocol to high oxidation state molybdenum(VI) N-heterocyclic carbene (NHC) alkylidyne complexes (NHC=1,3-diisopropylimidazol-2-ylidene, 1,3-dimethyl-4,5-R2 -imidazol-2-ylidene, R2 =H, Cl, CN) is reported. Steric limitations of the NHCs and the benzylidyne are described. All novel complexes were characterized by single crystal X-ray diffraction and solution NMR techniques. It was shown that all complexes presented here show activity in the self-metathesis of 1-phenyl-1-propyne at room temperature. To identify mechanistic differences, an experimental sequence to detect dissociation of ligands was developed. Results reveal dissociation of less electron-donating NHCs in course of the reaction. Mechanistic and reactivity differences were attributed to electronic and steric effects through Tolman's electronic parameter and the percentage of buried volume. Furthermore, Mo-1 containing the 1,3-dimethylimidazol-2-ylidene ligand showed good activity in self-metathesis reactions of p-substituted 1-phenyl-1-propynes with electron-donating moieties at room temperature.

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.

Comparison of classical reaction paths and tunneling paths studied with the semiclassical instanton theory.

Atom tunneling in the hydrogen atom transfer reaction of the 2,4,6-tri-tert-butylphenyl radical to 3,5-di-tert-butylneophyl, which has a short but strongly curved reaction path, was investigated using instanton theory. We found the tunneling path to deviate qualitatively from the classical intrinsic reaction coordinate, the steepest-descent path in mass-weighted Cartesian coordinates. To perform that comparison, we implemented a new variant of the predictor-corrector algorithm for the calculation of the intrinsic reaction coordinate. We used the reaction force analysis method as a means to decompose the reaction barrier into structural and electronic components. Due to the narrow energy barrier, atom tunneling is important in the abovementioned reaction, even above room temperature. Our calculated rate constants between 350 K and 100 K agree well with experimental values. We found a H/D kinetic isotope effect of almost 106 at 100 K. Tunneling dominates the protium transfer below 400 K and the deuterium transfer below 300 K. We compared the lengths of the tunneling path and the classical path for the hydrogen atom transfer in the reaction HCl + Cl and quantified the corner cutting in this reaction. At low temperature, the tunneling path is about 40% shorter than the classical path.

An Aluminum Fluoride Complex with an Appended Ammonium Salt as an Exceptionally Active Cooperative Catalyst for the Asymmetric Carboxycyanation of Aldehydes.

Al-F bonds are among the most stable σ bonds known, exhibiting an even higher bond energy than Si-F bonds. Despite a stability advantage and a potentially high Lewis acidity of Al-F complexes, they have not been described as structurally defined catalysts for enantioselective reactions. We show that Al-F salen complexes with appended ammonium moieties give exceptional catalytic activity in asymmetric carboxycyanations. In addition to aromatic aldehydes, enal and aliphatic substrates are well accepted. Turnover numbers up to around 104 were achieved, whereas with previous catalysts 101 -102 turnovers were typically attained. In contrast to Al-Me and Al-Cl salen complexes, the analogous Al-F species are remarkably stable towards air, water, and heat, and can be recovered unchanged after catalysis. They possess a considerably increased Lewis acidity as shown by DFT calculations.

Atom tunnelling in the reaction NH3+ + H2 → NH4+ + H and its astrochemical relevance.

The title reaction is involved in the formation of ammonia in the interstellar medium. We have calculated thermal rates including atom tunnelling using different rate theories. Canonical variational theory with microcanonically optimised multidimensional tunnelling was used for bimolecular rates, modelling the gas-phase reaction and also a surface-catalysed reaction of the Eley-Rideal type. Instanton theory provided unimolecular rates, which model the Langmuir-Hinshelwood type surface reaction. The potential energy was calculated on the CCSD(T)-F12 level of theory on the fly. We report thermal rates and H/D kinetic isotope effects. The latter have implications for observed H/D fractionation in molecular clouds. Tunnelling causes rate constants to be sufficient for the reaction to play a role in interstellar chemistry even at cryogenic temperature. We also discuss intricacies and limitations of the different tunnelling approximations to treat this reaction, including its pre-reactive minimum.

Reaction rates and kinetic isotope effects of H2 + OH → H2O + H.

We calculated reaction rate constants including atom tunneling of the reaction of dihydrogen with the hydroxy radical down to a temperature of 50 K. Instanton theory and canonical variational theory with microcanonical optimized multidimensional tunneling were applied using a fitted potential energy surface [J. Chen et al., J. Chem. Phys. 138, 154301 (2013)]. All possible protium/deuterium isotopologues were considered. Atom tunneling increases at about 250 K (200 K for deuterium transfer). Even at 50 K the rate constants of all isotopologues remain in the interval 4 ⋅ 10(-20) to 4 ⋅ 10(-17) cm(3) s(-1), demonstrating that even deuterated versions of the title reaction are possibly relevant to astrochemical processes in molecular clouds. The transferred hydrogen atom dominates the kinetic isotope effect at all temperatures.

Atom Tunneling in Chemistry.

Quantum mechanical tunneling of atoms is increasingly found to play an important role in many chemical transformations. Experimentally, atom tunneling can be indirectly detected by temperature-independent rate constants at low temperature or by enhanced kinetic isotope effects. In contrast, the influence of tunneling on the reaction rates can be monitored directly through computational investigations. The tunnel effect, for example, changes reaction paths and branching ratios, enables chemical reactions in an astrochemical environment that would be impossible by thermal transition, and influences biochemical processes.