Stiff-stilbene photoswitch ruptures bonds not by pulling but by local heating

JEDI: A versatile code for strain analysis of molecular and periodic systems under deformation

Stretching or compression can induce significant energetic, geometric, and spectroscopic changes in materials. To fully exploit these effects in the design of mechano- or piezo-chromic materials, self-healing polymers, and other mechanoresponsive devices, a detailed knowledge about the distribution of mechanical strain in the material is essential. Within the past decade, Judgement of Energy DIstribution (JEDI) analysis has emerged as a useful tool for this purpose. Based on the harmonic approximation, the strain energy in each bond length, bond angle, and dihedral angle of the deformed system is calculated using quantum chemical methods. This allows the identification of the force-bearing scaffold of the system, leading to an understanding of mechanochemical processes at the most fundamental level. Here, we present a publicly available code that generalizes the JEDI analysis, which has previously only been available for isolated molecules. Now, the code has been extended to two- and three-dimensional periodic systems, supramolecular clusters, and substructures of chemical systems under various types of deformation. Due to the implementation of JEDI into the Atomic Simulation Environment, the JEDI analysis can be interfaced with a plethora of program packages that allow the calculation of electronic energies for molecular systems and systems with periodic boundary conditions. The automated generation of a color-coded three-dimensional structure via the Visual Molecular Dynamics program allows insightful visual analyses of the force-bearing scaffold of the strained system.

ReaxFF-based nonadiabatic dynamics method for azobenzene derivatives

ReaxFF reactive force fields have been parameterized for the ground and first excited states of azobenzene and its derivatives. In addition, an extended set of ab initio reference data ensures wide applicability, including to azosystems in complex environments. Based on the optimized force fields, nonadiabatic surface hopping simulations produce photoisomerization quantum yields and decay times of azobenzene, both in the gas phase and in n-hexane solution, in reasonable agreement with higher level theory and experiment. The transferability to other azo-compounds is illustrated for different arylazopyrazoles as well as ethylene-bridged azobenzene. Moreover, it has been shown that the model can be easily extended to adsorbates on metal surfaces. The simulation of the ring-opening of cyclobutene triggered by the photoisomerization of azobenzene in a macrocycle highlights the advantages of a reactive force field model.

Mechanochemically Triggered Topology Changes in Expanded Porphyrins

A hitherto unexplored class of molecules for molecular force probe applications are expanded porphyrins. This work proves that mechanical force is an effective stimulus to trigger the interconversion between Hückel and Möbius topologies in [28]hexaphyrin, making these expanded porphyrins suitable to act as conformational mechanophores operating at mild (sub-1 nN) force conditions. A straightforward approach based on distance matrices is proposed for the selection of pulling scenarios that promote either the planar Hückel topology or the three lowest lying Möbius topologies. This approach is supported by quantum mechanochemical calculations. Force distribution analyses reveal that [28]hexaphyrin selectively allocates the external mechanical energy to molecular regions that trigger Hückel-Möbius interconversions, explaining why certain pulling scenarios favor the Hückel two-sided topology and others favor Möbius single-sided topologies. The meso-substitution pattern on [28]hexaphyrin determines whether the energy difference between the different topologies can be overcome by mechanical activation.

Effect of ring rize on photoisomerization properties of stiff stilbene macrocycles

A series of stiff stilbene macrocycles have been studied to investigate the possible impact of the macrocycle ring size on their photodynamic properties. The results show that reducing the ring size counteracts the photoisomerization ability of the macrocycles. However, even the smallest macrocycle studied (stiff stilbene subunits linked by a six carbon chain) showed some degree of isomerization when irradiated. DFT calculations of the energy differences between the E- and Z-isomers show the same trend as the experimental results. Interestingly the DFT study highlights that the energy difference between the E- and Z-isomers of even the largest macrocycle (linked by a twelve carbon chain) is significantly higher than that of the stiff stilbene unit itself. In general, it is indicated that addition of even a flexible chain to the stiff stilbene unit may significantly affect its photochemical properties and increase the photostability of the resulting macrocycle.

Franck-Condon Theory of Quantum Mechanochemistry

The spectroscopic Franck-Condon (FC) principle is extended to mechanochemistry. If the external force is applied rapidly (the sudden-force regime), then the transition between the potential energy surface and the force-modified potential energy surface is analogous to the optical electronic transition. Such a transition produces a nonequilibrium ensemble of vibrationally excited molecules. This excess of vibrational energy is another activation source in addition to the well-known reaction barrier modulation by the external force. In the same time, the nonequilibrium vibrational distribution implies nonstatistical kinetics of a mechanochemical transformation. Mechanochemical FC principle thus provides a conceptual picture for the sudden-force mechanochemistry and opens possibilities for quantitative calculations of the mechanochemical rates and mechanisms. Here we use it to compute the dissociation rates of a model diatomic molecule and to explain the selectivity in mechanochemical bond breaking in n-butane. The approach is predicted to be relevant for large-magnitude external forces, applied instantaneously. Otherwise, the excess vibrational energy will dissipate due to intramolecular vibrational redistribution and interaction with environment.

Force-induced retro-click reaction of triazoles competes with adjacent single-bond rupture

The highly controversial force-induced cycloreversion of 1,2,3-triazole, its well-known retro-click reaction, is shown to be possible only for 1,5-substituted triazoles, but competes with rupture of an adjacent single-bond. We draw this conclusion from both static and dynamic calculations under external mechanical forces applied to unsubstituted and 1,4- and 1,5-substituted triazoles. The JEDI (Judgement of Energy DIstribution) analysis, a quantum chemical tool quantifying the distribution of strain energy in mechanically deformed molecules, is employed to identify the key factors facilitating the force-induced retro-click reaction in these systems. For 1,4-substituted triazoles it is shown to be impossible, but the parallel alignment of the scissile bond in 1,5-substituted triazoles with the acting force makes it generally feasible. However, the weakness of the carbon-nitrogen bond connecting the triazole ring to the linker prevents selective cycloreversion.

Quantum Chemical Strain Analysis For Mechanochemical Processes

The use of mechanical force to initiate a chemical reaction is an efficient alternative to the conventional sources of activation energy, i.e., heat, light, and electricity. Applications of mechanochemistry in academic and industrial laboratories are diverse, ranging from chemical syntheses in ball mills and ultrasound baths to direct activation of covalent bonds using an atomic force microscope. The vectorial nature of force is advantageous because specific covalent bonds can be preconditioned for rupture by selective stretching. However, the influence of mechanical force on single molecules is still not understood at a fundamental level, which limits the applicability of mechanochemistry. As a result, many chemists still resort to rules of thumb when it comes to conducting mechanochemical syntheses. In this Account, we show that comprehension of mechanochemistry at the molecular level can be tremendously advanced by quantum chemistry, in particular by using quantum chemical force analysis tools. One such tool is the JEDI (Judgement of Energy DIstribution) analysis, which provides a convenient approach to analyze the distribution of strain energy in a mechanically deformed molecule. Based on the harmonic approximation, the strain energy contribution is calculated for each bond length, bond angle and dihedral angle, thus providing a comprehensive picture of how force affects molecules. This Account examines the theoretical foundations of quantum chemical force analysis and provides a critical overview of the performance of the JEDI analysis in various mechanochemical applications. We explain in detail how this analysis tool is to be used to identify the "force-bearing scaffold" of a distorted molecule, which allows both the rationalization and the optimization of diverse mechanochemical processes. More precisely, we show that the inclusion of every bond, bending and torsion of a molecule allows a particularly insightful discussion of the distribution of mechanical strain in deformed molecules. We illustrate the usefulness of the JEDI analysis by rationalizing the finding that a knot tremendously weakens a polymer strand via a "choking" motion of the torsions in the curved part of the knot, thus leading to facilitated bond rupture in the immediate vicinity of the knot. Moreover, we demonstrate that the JEDI analysis can be exploited to devise methods for the stabilization of inherently strained molecules. In addition to applications in the electronic ground state, the JEDI analysis can also be used in the electronically excited state to determine the mechanical energy that a molecular photoswitch can release into its environment during photoisomerization. This approach allows the quantification of the mechanical efficiency of a photoswitch, i.e., the part of the energy that becomes available for the motion into a specific direction, which enables us to judge whether a photoswitch is capable of performing a desired switching function.

Advances in Quantum Mechanochemistry: Electronic Structure Methods and Force Analysis

In quantum mechanochemistry, quantum chemical methods are used to describe molecules under the influence of an external force. The calculation of geometries, energies, transition states, reaction rates, and spectroscopic properties of molecules on the force-modified potential energy surfaces is the key to gain an in-depth understanding of mechanochemical processes at the molecular level. In this review, we present recent advances in the field of quantum mechanochemistry and introduce the quantum chemical methods used to calculate the properties of molecules under an external force. We place special emphasis on quantum chemical force analysis tools, which can be used to identify the mechanochemically relevant degrees of freedom in a deformed molecule, and spotlight selected applications of quantum mechanochemical methods to point out their synergistic relationship with experiments.

Can Strained Hydrocarbons Be "Forced" To Be Stable?

Many strained hydrocarbons are prone to isomerization, dimerization, and trimerization under normal laboratory conditions. Here we investigate a method to stabilize angle-strained cycloalkynes by applying a mechanical pulling force to the carbon atoms adjacent to the triple bond, which partially linearizes the C≡C-C bond angles. We discuss various methods of applying such pulling forces, including photoswitches and incorporation into additional strained macrocycles. We use the computational JEDI (Judgement of Energy DIstribution) analysis to quantify the distribution of energy in strained cycloheptyne and judge the change in stability upon application of an external force via isodesmic and homodesmotic reactions. We find that cycloheptyne can indeed be stabilized by external forces. However, the force generated by photoswitches during isomerization is too low to lead to a significant stabilization of the molecule. Hence, stronger forces are needed, which can be achieved by incorporating cycloheptyne into a second strained macrocycle.

Response to Comment on T. Stauch, A. Dreuw, “Stiff-stilbene photoswitch ruptures bonds not by pulling but by local heating”, Phys. Chem. Chem. Phys., 2016, 18, 15848

Computational analysis of mechanical stress energy revealed that photoisomerization of stiff-stilbene cannot initiate the ring opening of cyclobutene via a purely mechanical pathway.

Comment on T. Stauch, A. Dreuw, “Stiff-stilbene photoswitch ruptures bonds not by pulling but by local heating”, Phys. Chem. Chem. Phys., 2016,18, 15848

We explain why the key premise of the above-cited paper (“local heating”) is based on misunderstanding of how fast solutes thermalize, of how molecular strain affects reaction rates and of the role of force in the modern models of mechanochemical kinetics.