Shear-activated chemisorption and association of cyclic organic molecules

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.

Reactive Simulations of Silica Functionalization with Aromatic Hydrocarbons

Reactive molecular dynamics simulations are used to model the covalent functionalization of amorphous silica with aromatic hydrocarbons. Simulations show that the surface density of silanol-terminated phenyl, naphthyl, and anthracenyl molecules is lower than the maximum value calculated based on molecule geometry, and the simulation densities decrease faster with the number of aromatic rings than the geometric densities. The trends are analyzed in terms of the surface-silanol bonding configurations, tilt angles, local conformational ordering, and aggregation of surface-bound molecules under steady-state conditions. Results show that the surface density is affected by both the size and symmetry of the aromatic hydrocarbons. The correlations among bonding, orientation, and surface density identified here may guide the selection or design of molecules for functionalized surfaces.