Mechanochemistry of phosphate esters confined between sliding iron surfaces.

The molecular structure of lubricant additives controls not only their adsorption and dissociation behaviour at the nanoscale, but also their ability to reduce friction and wear at the macroscale. Here, we show using nonequilibrium molecular dynamics simulations with a reactive force field that tri(s-butyl)phosphate dissociates much faster than tri(n-butyl)phosphate when heated and compressed between sliding iron surfaces. For both molecules, dissociative chemisorption proceeds through cleavage of carbon-oxygen bonds. The dissociation rate increases exponentially with temperature and stress. When the rate-temperature-stress data are fitted with the Bell model, both molecules have similar activation energies and activation volumes and the higher reactivity of tri(s-butyl)phosphate is due to a larger pre-exponential factor. These observations are consistent with experiments using the antiwear additive zinc dialkyldithiophosphate. This study represents a crucial step towards the virtual screening of lubricant additives with different substituents to optimise tribological performance.

Adsorption of Organic Friction Modifier Additives.

Organic friction modifier additives (OFMs) are surfactant molecules added to engine oils to reduce friction in the boundary lubrication regime. They are thought to work by forming an absorbed layer, which provides low friction. This paper studied the relationship between the adsorption of OFMs and their friction and wear-reducing properties in a rubbing contact formed by a stationary glass ball and a rotating silicon disk under the boundary lubrication regime. The effect of molecular structure was investigated by using OFMs of various tail saturation and head group chemistry. OFMs tested were oleic acid, octadecylamine, oleylamine, and glycerol monooleate. The thickness of an OFM-adsorbed layer in hexadecane, examined in situ by spectroscopic ellipsometry and quartz crystal microbalance (QCM), depends on the molecular structure and the concentration of the OFM. As expected, the saturated, linear chain gives the thickest film. A critical OFM layer thickness of about 0.6 nm is necessary to achieve low initial and maximum friction. The thicker OFM layers are accompanied by narrower wear tracks, which are rougher than the wider, smoother wear tracks formed with thinner OFM layers. The interplay between the thickness of the OFM layer and wear track surface roughness results in all OFM layers having similar steady state friction. This shows that the apparent effect of OFM depends on the stage of the rubbing test: initially on friction, and then subsequently on surface damage. Despite OFMs and the base oil having similar refractive indices, ellipsometry was found to be a suitable technique for examining the adsorption of OFM additives from an oil-based solution and showed reasonable correlation with QCM results.

Mechanochemistry of Zinc Dialkyldithiophosphate on Steel Surfaces under Elastohydrodynamic Lubrication Conditions.

Zinc dialkyldithiophosphate (ZDDP) is added to engine lubricants to reduce wear and ensure reliable operation. ZDDP reacts under rubbing conditions to form protective zinc/iron phosphate tribofilms on steel surfaces. Recently, it has been demonstrated that this process can be promoted by applied stresses in lubricated contacts, as well as temperature, and is thus mechanochemical in origin. In this study, a tribology test rig, capable of applying very high loads, has been developed to generate ZDDP tribofilms under full-film elastohydrodynamic lubrication (EHL) conditions in steel/steel ball-on-disk contacts. This provides a well-defined temperature and stress environment with negligible direct asperity contact in which to study mechanochemical processes. ZDDPs with branched primary and secondary alkyl substituents have been studied in three base oils, two with high EHL friction and one with low EHL friction. In the high EHL friction base oils, the tribofilm growth rate increases exponentially with shear stress and temperature for both ZDDPs, as predicted by a stress-augmented thermal activation model. Conversely, under otherwise identical conditions, negligible ZDDP tribofilm formation takes place in the low EHL friction base oil. This confirms that the ZDDP reaction is driven by macroscopic shear stress rather than hydrostatic pressure. The secondary ZDDP forms tribofilms considerably faster than the primary ZDDP under equivalent conditions, suggesting that the initial decomposition reaction is the rate-determining step for tribofilm formation. The rate of tribofilm growth is independent of ZDDP concentration over the range studied, indicating that this process follows zero-order kinetics. Under full-film EHL conditions, ZDDP tribofilm formation is promoted by macroscopic shear stress applied through the base oil molecules, which induces asymmetric stress on adsorbed ZDDP molecules to promote their decomposition and initiate rapid phosphate polymerization.

Nonequilibrium Molecular Dynamics Simulations of Organic Friction Modifiers Adsorbed on Iron Oxide Surfaces.

For the successful development and application of lubricants, a full understanding of the nanoscale behavior of complex tribological systems is required, but this is difficult to obtain experimentally. In this study, we use nonequilibrium molecular dynamics (NEMD) simulations to examine the atomistic structure and friction properties of commercially relevant organic friction modifier (OFM) monolayers adsorbed on iron oxide surfaces and lubricated by a thin, separating layer of hexadecane. Specifically, acid, amide, and glyceride OFMs, with saturated and Z-unsaturated hydrocarbon tail groups, are simulated at various surface coverages and sliding velocities. At low and medium coverage, the OFMs form liquidlike and amorphous monolayers, respectively, which are significantly interdigitated with the hexadecane lubricant, resulting in relatively high friction coefficients. At high coverage, solidlike monolayers are formed for all of the OFMs, which, during sliding, results in slip planes between well-defined OFM and hexadecane layers, yielding a marked reduction in the friction coefficient. When present at equal surface coverage, OFMs with saturated and Z-unsaturated tail groups are found to yield similar structure and friction behavior. OFMs with glyceride head groups yield significantly lower friction coefficients than amide and particularly carboxylic acid head groups. For all of the OFMs and coverages simulated, the friction coefficient is found to increase linearly with the logarithm of sliding velocity; however, the gradient of this increase depends on the coverage. The structure and friction details obtained from these simulations agree well with experimental results and also shed light on the relative tribological performance of these OFMs through nanoscale structural variations. This has important implications in terms of the applicability of NEMD to aid the development of new formulations to control friction.

A Comparison of Classical Force-Fields for Molecular Dynamics Simulations of Lubricants.

For the successful development and application of lubricants, a full understanding of their complex nanoscale behavior under a wide range of external conditions is required, but this is difficult to obtain experimentally. Nonequilibrium molecular dynamics (NEMD) simulations can be used to yield unique insights into the atomic-scale structure and friction of lubricants and additives; however, the accuracy of the results depend on the chosen force-field. In this study, we demonstrate that the use of an accurate, all-atom force-field is critical in order to; (i) accurately predict important properties of long-chain, linear molecules; and (ii) reproduce experimental friction behavior of multi-component tribological systems. In particular, we focus on n-hexadecane, an important model lubricant with a wide range of industrial applications. Moreover, simulating conditions common in tribological systems, i.e., high temperatures and pressures (HTHP), allows the limits of the selected force-fields to be tested. In the first section, a large number of united-atom and all-atom force-fields are benchmarked in terms of their density and viscosity prediction accuracy of n-hexadecane using equilibrium molecular dynamics (EMD) simulations at ambient and HTHP conditions. Whilst united-atom force-fields accurately reproduce experimental density, the viscosity is significantly under-predicted compared to all-atom force-fields and experiments. Moreover, some all-atom force-fields yield elevated melting points, leading to significant overestimation of both the density and viscosity. In the second section, the most accurate united-atom and all-atom force-field are compared in confined NEMD simulations which probe the structure and friction of stearic acid adsorbed on iron oxide and separated by a thin layer of n-hexadecane. The united-atom force-field provides an accurate representation of the structure of the confined stearic acid film; however, friction coefficients are consistently under-predicted and the friction-coverage and friction-velocity behavior deviates from that observed using all-atom force-fields and experimentally. This has important implications regarding force-field selection for NEMD simulations of systems containing long-chain, linear molecules; specifically, it is recommended that accurate all-atom potentials, such as L-OPLS-AA, are employed.