A Two-Step Baromechanical Cycle for Repeated Activation and Deactivation of Mechanophores

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.

Calculating high-pressure vibrational frequencies analytically with the extended hydrostatic compression force field approach

We report the implementation of the analytical Hessian for the mechanochemical extended hydrostatic compression force field method in the Q-Chem program package. To verify the implementation, the analytical Hessian was compared with finite difference calculations. In addition, we calculated the pressure dependency of the Raman active vibrational modes of methane, ethane, and hydrogen, as well as all IR and Raman active modes of Buckminsterfullerene, and compared the results with experimental and theoretical data. Our implementation paves the way for the analysis of geometric points on a pressure-deformed potential energy surface and provides a straightforward model to calculate the vibrational properties of molecules under high pressure.

Shear-activation of mechanochemical reactions through molecular deformation

Mechanical stress can directly activate chemical reactions by reducing the reaction energy barrier. A possible mechanism of such mechanochemical activation is structural deformation of the reactant species. However, the effect of deformation on the reaction energetics is unclear, especially, for shear stress-driven reactions. Here, we investigated shear stress-driven oligomerization reactions of cyclohexene on silica using a combination of reactive molecular dynamics simulations and ball-on-flat tribometer experiments. Both simulations and experiments captured an exponential increase in reaction yield with shear stress. Elemental analysis of ball-on-flat reaction products revealed the presence of oxygen in the polymers, a trend corroborated by the simulations, highlighting the critical role of surface oxygen atoms in oligomerization reactions. Structural analysis of the reacting molecules in simulations indicated the reactants were deformed just before a reaction occurred. Quantitative evidence of shear-induced deformation was established by comparing bond lengths in cyclohexene molecules in equilibrium and prior to reactions. Nudged elastic band calculations showed that the deformation had a small effect on the transition state energy but notably increased the reactant state energy, ultimately leading to a reduction in the energy barrier. Finally, a quantitative relationship was developed between molecular deformation and energy barrier reduction by mechanical stress.

Wandering through quantum-mechanochemistry: from concepts to reactivity and switches

Mechanochemistry has experienced a renaissance in recent years witnessing, at the molecular level, a remarkable interplay between theory and experiment. Molecular mechanochemistry has welcomed a broad spectrum of quantum-chemical methods to evaluate the influence of an external mechanical force on molecular properties. In this contribution, an overview is given on recent work on quantum mechanochemistry in the Brussels Quantum Chemistry group (ALGC). The effect of an external force was scrutinized both in fundamental topics, like reactivity descriptors in Conceptual DFT, and in applied topics, such as designing molecular force probes and tuning the stereoselectivity of certain types of reactions. In the conceptual part, a brief overview of the techniques introducing mechanical forces into a quantum-mechanical description of a molecule is followed by an introduction to conceptual DFT. The evolution of the electronic chemical potential (or electronegativity), chemical hardness and electrophilicity are investigated when a chemical bond in a series of diatomics is put under mechanical stress. Its counterpart, the influence of mechanical stress on bond angles, is analyzed by varying the strain present in alkyne triple bonds by applying a bending force, taking the strain promoted alkyne-azide coupling cycloaddition as an example. The increase of reactivity of the alkyne upon bending is probed by Fukui functions and the local softness. In the applied part, a new molecular force probe is presented based on an intramolecular 6π-electrocyclization in constrained polyenes operating under thermal conditions. A cyclic process is conceived where ring opening and closure are triggered by applying or removing an external pulling force. The efficiency of mechanical activation strongly depends on the magnitude of the applied force and the distance between the pulling points. The idea of pulling point distances as a tool to identify new mechanochemical processes is then tested in [28]hexaphyrins with an intricate equilibrium between Möbius aromatic and Hückel antiaromatic topologies. A mechanical force is shown to trigger the interconversion between the two topologies, using the distance matrix as a guide to select appropriate pulling points. In a final application, the Felkin-Anh model for the addition of nucleophiles to chiral carbonyls under the presence of an external mechanical force is scrutinized. By applying a force for restricting the conformational freedom of the chiral ketone, otherwise inaccessible reaction pathways are promoted on the force-modified potential energy surfaces resulting in a diastereoselectivity different from the force-free reaction.