Todd J. Martínez

Force-induced activation of covalent bonds in mechanoresponsive polymeric materials

Mechanochemical transduction enables an extraordinary range of physiological processes such as the sense of touch, hearing, balance, muscle contraction, and the growth and remodelling of tissue and bone. Although biology is replete with materials systems that actively and functionally respond to mechanical stimuli, the default mechanochemical reaction of bulk polymers to large external stress is the unselective scission of covalent bonds, resulting in damage or failure. An alternative to this degradation process is the rational molecular design of synthetic materials such that mechanical stress favourably alters material properties. A few mechanosensitive polymers with this property have been developed; but their active response is mediated through non-covalent processes, which may limit the extent to which properties can be modified and the long-term stability in structural materials. Previously, we have shown with dissolved polymer strands incorporating mechanically sensitive chemical groups—so-called mechanophores—that the directional nature of mechanical forces can selectively break and re-form covalent bonds. We now demonstrate that such force-induced covalent-bond activation can also be realized with mechanophore-linked elastomeric and glassy polymers, by using a mechanophore that changes colour as it undergoes a reversible electrocyclic ring-opening reaction under tensile stress and thus allows us to directly and locally visualize the mechanochemical reaction. We find that pronounced changes in colour and fluorescence emerge with the accumulation of plastic deformation, indicating that in these polymeric materials the transduction of mechanical force into the ring-opening reaction is an activated process. We anticipate that force activation of covalent bonds can serve as a general strategy for the development of new mechanophore building blocks that impart polymeric materials with desirable functionalities ranging from damage sensing to fully regenerative self-healing.

Quantum Chemistry on Graphical Processing Units. 3. Analytical Energy Gradients, Geometry Optimization, and First Principles Molecular Dynamics.

We demonstrate that a video gaming machine containing two consumer graphical cards can outpace a state-of-the-art quad-core processor workstation by a factor of more than 180× in Hartree-Fock energy + gradient calculations. Such performance makes it possible to run large scale Hartree-Fock and Density Functional Theory calculations, which typically require hundreds of traditional processor cores, on a single workstation. Benchmark Born-Oppenheimer molecular dynamics simulations are performed on two molecular systems using the 3-21G basis set - a hydronium ion solvated by 30 waters (94 atoms, 405 basis functions) and an aspartic acid molecule solvated by 147 waters (457 atoms, 2014 basis functions). Our GPU implementation can perform 27 ps/day and 0.7 ps/day of ab initio molecular dynamics simulation on a single desktop computer for these systems.

Conical intersections and double excitations in time-dependent density functional theory

There is a clear need for computationally inexpensive electronic structure theory methods which can model excited state potential energy surfaces. Time-dependent density functional theory (TDDFT) has emerged as one of the most promising contenders in this context. Many previous tests have concentrated on vertical excitation energies, which can be compared to experimental absorption maxima. Here, we focus attention on more global aspects of the resulting potential energy surfaces, especially conical intersections which play a key role in photochemical mechanisms. We introduce a new method for minimal energy conical intersection (MECI) searches which does not require knowledge of the non-adiabatic coupling vector. Using this new method, we compute MECI geometries with multi-state complete active space perturbation theory (MS-CASPT2) and TDDFT. We show that TDDFT in the linear response and adiabatic approximations can predict MECI geometries and energetics quite accurately, but that there are a number of qualitative deficiencies which need to be addressed before TDDFT can be used routinely in photochemical problems. †Dedicated to Professor M. A. Robb on the occasion of his 60th birthday.

Quantum Chemistry on Graphical Processing Units. 1. Strategies for Two-Electron Integral Evaluation

Modern videogames place increasing demands on the computational and graphical hardware, leading to novel architectures that have great potential in the context of high performance computing and molecular simulation. We demonstrate that Graphical Processing Units (GPUs) can be used very efficiently to calculate two-electron repulsion integrals over Gaussian basis functions [Formula: see text] the first step in most quantum chemistry calculations. A benchmark test performed for the evaluation of approximately 10(6) (ss|ss) integrals over contracted s-orbitals showed that a naïve algorithm implemented on the GPU achieves up to 130-fold speedup over a traditional CPU implementation on an AMD Opteron. Subsequent calculations of the Coulomb operator for a 256-atom DNA strand show that the GPU advantage is maintained for basis sets including higher angular momentum functions.

Quantum Chemistry on Graphical Processing Units. 2. Direct Self-Consistent-Field Implementation

We demonstrate the use of graphical processing units (GPUs) to carry out complete self-consistent-field calculations for molecules with as many as 453 atoms (2131 basis functions). Speedups ranging from 28× to 650× are achieved as compared to a mature third-party quantum chemistry program (GAMESS) running on a traditional CPU. The computational organization used to construct the Coulomb and exchange operators is discussed. We also present results using three GPUs in parallel, combining coarse and fine-grained parallelism.

Trapping a Diradical Transition State by Mechanochemical Polymer Extension

Forced Open Traditionally, the study of reaction chemistry has relied on random encounters between molecules to initiate the proceedings. Heating and stirring increase the power and frequency of such encounters but provide little finer control. Very recently, chemists have learned how to initiate reactions more directly by embedding precursors in the backbone of a polymer large enough to manipulate with shear forces. Lenhardt et al. (p. 1057) applied this technique to a cyclopropyl ring-opening reaction. When the strained triangular carbon rings were embedded within a polymer, shear force applied by sonication ruptured their bonds as the polymer backbone stretched. The taut polymer then redirected the ring-opened intermediates toward a product differently arranged from that generated by simple heating. Strained carbon rings in the backbone of a polymer can be opened by external application of shear forces. Transition state structures are central to the rates and outcomes of chemical reactions, but their fleeting existence often leaves their properties to be inferred rather than observed. By treating polybutadiene with a difluorocarbene source, we embedded gem-difluorocyclopropanes (gDFCs) along the polymer backbone. We report that mechanochemical activation of the polymer under tension opens the gDFCs and traps a 1,3-diradical that is formally a transition state in their stress-free electrocyclic isomerization. The trapped diradical lives long enough that we can observe its noncanonical participation in bimolecular addition reactions. Furthermore, the application of a transient tensile force induces a net isomerization of the trans-gDFC into its less-stable cis isomer, leading to the counterintuitive result that the gDFC contracts in response to a transient force of extension.

Photodynamics of ethylene: ab initio studies of conical intersections

Abstract We have used extended basis sets and configuration interaction wave functions to systematically characterize the excited state potential energy surfaces of ethylene, including conical intersections between electronic states. The results are consistent with our previous ab initio multiple spawning simulations of ethylene photodynamics and electronic spectra. The C–C bond on the optically accessible V state is extended in planar geometries, suggesting a role for C–C stretching in the electronic absorption spectrum. A cascade of conical intersections connecting the V state in the Franck–Condon region to each of the low-lying Rydberg states has been identified, in addition to intersections connecting the excited state manifold back to the ground state. The D 2d twisted geometry of ethylene is found to be a saddle point, not a local minimum. Pyramidalization of one of the methylene units in twisted ethylene is found to be favorable, leading to a conical intersection. We have identified and characterized eight conical intersections involving the V state which are likely to be relevant in the photochemistry of ethylene.

First Principles Dynamics and Minimum Energy Pathways for Mechanochemical Ring Opening of Cyclobutene

We use ab initio steered molecular dynamics to investigate the mechanically induced ring opening of cyclobutene. We show that the dynamical results can be considered in terms of a force-modified potential energy surface (FMPES). We show how the minimal energy paths for the two possible competing conrotatory and disrotatory ring-opening reactions are affected by external force. We also locate minimal energy pathways in the presence of applied external force and show that the reactant, product, and transition state geometries are altered by the application of external force. The largest effects are on the transition state geometries and barrier heights. Our results provide a framework for future investigations of the role of external force on chemical reactivity.

Ab initio molecular dynamics study of cis–trans photoisomerization in ethylene

Abstract We have used ab initio multi-electronic state molecular dynamics to study the photoinduced cis–trans isomerization of ethylene. The initial motion on the excited state is a stretching of the CC bond and the photoisomerization begins within ∼50 fs of optical excitation. Quenching to the ground electronic state is found to be ultrafast and proceeds from an ionic state via a conical intersection. Accessing the conical intersection requires pyramidalization of one of the methylene groups and this can happen only after energy is funneled from the twisting mode into the pyramidalization mode.

Masked cyanoacrylates unveiled by mechanical force.

Mechanical damage of polymers is often a destructive and irreversible process. However, desirable outcomes may be achieved by controlling the location of chain cleavage events through careful design and incorporation of mechanically active chemical moieties known as mechanophores. It is possible that mechanophores can be used to generate reactive intermediates that can autopolymerize or cross-link, thus healing mechanically induced damage. Herein we report the generation of reactive cyanoacrylate units from a dicyanocyclobutane mechanophore located near the center of a polymer chain. Because cyanoacrylates (which are used as monomers in the preparation of superglue) autopolymerize, the generated cyanoacrylate-terminated polymers may be useful in self-healing polymers. Sonication studies of polymers with the mechanophore incorporated into the chain center have shown that selective cleavage of the mechanophore occurs. Trapping experiments with an amine-based chromophore support cyanoacrylate formation. Additionally, computational studies of small-molecule models predict that force-induced bond cleavage should occur with greater selectivity for the dicyanocyclobutane mechanophore than for a control molecule.