Mechanochemical activation has created new opportunities for applications such as solvent-free chemical synthesis, polymer processing, and lubrication. However, mechanistic understanding of these processes is still limited because the mechanochemical response of a system is a complex function of many variables, including the direction of applied stress and the chemical features of the reactants in non-equilibrium conditions. Here, we studied shear-activated reactions of simple cyclic organic molecules to isolate the effect of chemical structure on reaction yield and pathway. Reactive molecular dynamics simulations were used to model methylcyclopentane, cyclohexane, and cyclohexene subject to pressure and shear stress between silica surfaces. Cyclohexene was found to be more susceptible to mechanochemical activation of oxidative chemisorption and subsequent oligomerization reactions than either methylcyclopentane or cyclohexane. The oligomerization trend was consistent with shear-driven polymerization yield measured in ball-on-flat sliding experiments. Analysis of the simulations showed the distribution of carbon atom sites at which oxidative chemisorption occurred and identified the double bond in cyclohexene as being the origin of its shear susceptibility. Lastly, the most common reaction pathways for association were identified, providing insight into how the chemical structures of the precursor molecules determined their response to mechanochemical activation.
more »
« less
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.
more » « less- Award ID(s):
- 2038499
- PAR ID:
- 10489569
- Publisher / Repository:
- Nature Publishing Group
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 14
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Electrocatalysis plays a key role in the development of renewable energy technologies. Along with the design of electrocatalysts, the microenvironment around catalytic sites has received increasing attention because it affects the distribution and mass transport of reaction species and impacts the reaction kinetics. In this Concept article, we highlight some mechanistic insights into the effect of microenvironment on gas‐involving electrocatalytic reactions, including CO2reduction, 2‐electron oxygen reduction, and hydrazine oxidation, demonstrating their sensitivity to the wetting properties of microenvironment. For reactions with a gaseous reactant, a moderately hydrophobic microenvironment can greatly enhance the mass transport of gaseous species to accelerate the reaction kinetics while improving the stability of gas‐diffusion electrodes. In contrast, for reactions with a liquid reactant and gaseous product, a hydrophilic microenvironment improves the exposure of catalytic sites to the reactant, while a hydrophobic microenvironment benefits the reaction on the other end by accelerating the diffusion and detachment of generated gas bubbles, which would otherwise block the catalytic sites from the reactant. These understandings and insights can provide important guidelines on the control and optimization of catalyst microenvironment for the development of efficient electrolyzers.
-
Tribochemistry, which is another name for mechanochemistry driven by shear, deals with complex and dynamic interfacial processes that can lead to surface wear or formation of beneficial tribofilms. For better mechanistic understanding of these processes, we investigated the reactivity of tribopolymerization of organic molecules with different internal ring strain (methylcyclopentane, cyclohexane, and cyclohexene) on a stainless steel (SS) surface in inert (N2), oxidizing (O2), and reducing (H2) environments at room temperature. On the clean stainless steel surface, precursor molecules were found to physisorb with a broad range of molecular orientations. In inert and reducing environments, the strain-free cyclohexane showed the lowest tribochemical activity among the three molecules tested. Compared to the N2 environment, the tribochemical activity in H2 was suppressed. In the O2 environment, only cyclohexene produced tribofilms and methylcyclopentane while cyclohexane did not. When tribofilms were analyzed with Raman spectroscopy, the spectral features of diamond-like carbon (DLC) or amorphous carbon (a-C) were observed due to photochemical degradation of triboproducts. Based on infrared spectroscopy, tribofilms were found to be organic polymers containing oxygenated groups. Whenever polymeric tribrofilms were produced, wear volume was suppressed by orders of magnitudes but not completely to zero. These results support previously suggested mechanisms which involve surface oxygen as a reactant species in the tribopolymerization process.more » « less
-
more » « less
Computational fluid dynamics models often employ the free shear boundary condition at free surfaces, a result from the continuity of the stress and the large viscosity contrast at liquid–gas interfaces. This study leverages nonequilibrium molecular dynamics simulations to investigate the validity of the free shear boundary condition on the exposed surface of a liquid meniscus at the nanoscale. The primary objective is elucidating the fundamental mechanisms and behavior of fluid interactions within a capillary meniscus formed between two carbon nanotubes (CNTs) in shear-driven flow. Shear-driven flow simulations were conducted by varying the velocity of a solid slab to induce different shear rates in the adjacent water molecules. The results demonstrate, for the first time, negligible shear at the free surface, supporting the free shear assumption from the nanoscale point of view. A force balance analysis reveals that capillary and surface tension forces dominate within the meniscus, dictating its shape and stability. Meniscus deformation was observed and primarily attributed to interatomic interactions between water molecules and CNTs, driven by a combination of short-range repulsive forces and van der Waals attractions. The minimal contribution from shear forces suggests that interatomic forces, rather than applied shear stress, are the primary drivers of the meniscus deformation. These findings offer valuable insights into fluid behavior and a sound fundamental analysis of the free shear boundary condition at the nanoscale.
-
This paper shows how the effect of combined normal and shear stresses on the rates of tribochemical reactions can be calculated using Evans-Polanyi (E-P) perturbation theory. The E-P approach is based on transition-state theory, where the rate of reaction is taken to be proportional to the concentration of activated complex. The equilibrium constant depends on the molar Gibbs free energy change between the initial- and transition-states, which, in turn, depends on the stresses. E-P theory has been used previously to successfully calculate the effects of normal stresses on reaction rates. In this case, ln(Rate) varies linearly with stress with a slope given by an activation volume, which broadly corresponds to the volume difference between the reactant and activated complex. An advantage of E-P theory is that it can calculate the influence of several perturbations, for example, the normal stress dependence of the shear stress during sliding. In this paper, E-P theory is used to calculate shear-induced, tribochemical reaction rates. The results depend on four elementary activation volumes for different contributions to the Gibbs free energy: two of them due to normal and shear stresses for sliding over the surface and two more for the surface reaction. The results of the calculations show that there is a linear dependence of ln(Rate) on the normal stress but that the coefficient of proportionality between the ln(Rate) and the normal stress now has contributions from all elementary-step activation volumes. Counterintuitively, the analysis predicts that the ln(Rate)-normal stress evolution tends, at zero normal stress, to an asymptotic rate constant that depends on sliding velocity and differs from the thermal reaction rate. The theoretical prediction is verified for the shear-induced decomposition of ethyl thiolate species adsorbed on a Cu(100) single crystal substrate that decomposes by C‒S bond cleavage. The theoretical analyses show that tribochemical reactions can be influenced by either just normal stresses or by a combination of normal and shear stresses, but that the latter effect is much greater. Finally, it is predicted that there should be a linear relationship between the activation energy and the logarithm of the pre-exponential factor of the asymptotic rate constant.more » « less
