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   
                    
                            
                            A Molecular-Scale Analysis of Pressure-Dependent Sliding Shear Stresses
                        
                    
    
            Stress-modified activated processes are analyzed using a model first proposed by Evans and Polanyi that uses transition-state theory to calculate the effect of some perturbation, described by an intensive variable, \(I\), on the reaction rate. They suggested that the rate constant depended primarily on the equilibrium between the transition state and the reactant, which, in turn, depends on the effect of the perturbation \(I\) on the Gibbs free energy, \(G=U-TS+IC\), where \(C\) is a variable conjugate to \(I\). For example, in the case of a hydrostatic pressure \(P\), the conjugate variable is the volume, \(-V\). This allows a pressure-dependent rate to be calculated from the equilibrium constant between the reactant and transition state. Advantages to this approach are that the analysis is independent of the pathway between the two states and can simultaneously include the effect of multiple perturbations. These ideas are applied to the Prandtl–Tomlinson model, which analyses the force-induced transition rate over a surface energy barrier. The Evans–Polanyi analysis is independent of the shape of the sliding potential and merely requires the locations of the initial and transition states. It also allows the effects of both normal and shear stresses to be analyzed to identify the molecular origins of the well-known pressure-dependent shear stress: \(\tau ={\tau }_{0}+{\mu }_{L}P\), where \({\tau }_{0}\) is a pressure-independent stress. The analysis reveals that \({\mu }_{L}\) depends on the molecular corrugation of the potential and that \({\tau }_{0}\) is velocity dependent, in accord with an empirical equation proposed by Briscoe and Evans. 
        more » 
        « less   
        
    
                            - Award ID(s):
- 2020525
- PAR ID:
- 10535335
- Publisher / Repository:
- Springer
- Date Published:
- Journal Name:
- Tribology Letters
- Volume:
- 71
- Issue:
- 4
- ISSN:
- 1023-8883
- Page Range / eLocation ID:
- 121
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
- 
            
- 
            Abstract The viscosity of fluids and their dependence on shear rate, known as shear thinning, plays a critical role in applications ranging from lubricants and coatings to biomedical and food-processing industries. Traditional models such as the Carreau and Eyring theories offer competing explanations for shear-thinning behavior. The Carreau model attributes viscosity reduction to molecular distortions, while the Eyring model describes shear thinning as a stress-induced transition over an activation energy barrier. This work proposes an extended-Eyring model that incorporates stress-dependent activation volumes, bridging key aspects of both theories. In modifying transition-state theory by using an Evans-Polanyi perturbation analysis, we derive a generalized viscosity equation that accounts for the molecular-scale rearrangements governing fluid flow. The model is validated against computational and experimental data, including shear-thinning behavior of pure squalane and polyethylene oxide (PEO) aqueous solutions. Comparative analysis with Carreau-Yasuda and conventional Eyring models demonstrates excellent accuracy in predicting viscosity trends over a wide range of shear rates. The introduction of stress-dependent activation volumes provides a description of molecular exchange kinetics accounting for structural reorganization under shear. These findings offer a unified framework for modeling shear thinning and have broad implications for designing advanced lubricants, polymer solutions, and complex fluids with tailored flow properties. Graphical Abstractmore » « less
- 
            Abstract 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
- 
            Abstract A rheological model for loose granular media is developed to capture both solid-like and fluid-like responses during shearing. The proposed model is built by following the mathematical structure of an extended Kelvin–Voigt model, where an elastic spring and plastic slider act in parallel to a viscous damper. This arrangement requires the partition of the total stress into rate-independent and rate-dependent stress components. To model the solid-like behavior, a simple frictional plasticity model is adopted without modifications, thus contributing to the rate-independent stress. Instead, the fluid-like or rate-dependent stress is further decomposed into deviatoric and volumetric parts, by proposing a new formulation based on a combination of the m(I) relation, originally developed under pressure-controlled shear, with a pressure-shear rate relation derived under volume-controlled shear. The proposed formulation allows the model to capture both the increase in the friction coefficient and the enhanced dilation at high shear rates. High-fidelity simulation data, obtained from discrete element method and multiscale modelling, are used to evaluate the performance of the proposed constitutive model. The model provides accurate results under both drained and undrained simple shear paths across a wide range of shear rates. Furthermore, it successfully reproduces at much lower computational cost the flowslide mobility computed through multiscale simulations, which is primarily regulated by the shear rate dependence of the material properties during the dynamic runout stage.more » « less
 An official website of the United States government
An official website of the United States government 
				
			 
					 
					
 
                                    