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.

For the successful development and application of novel lubricant additives, a full understanding of their tribological behaviour at the nanoscale is required, but this can be difficult to obtain experimentally. In this study, nonequilibrium molecular dynamics simulations are used to examine the friction and wear reduction mechanisms of promising carbon nanoparticle friction modifier additives. Specifically, the friction and wear behaviour of carbon nanodiamonds (CNDs) and carbon nano-onions (CNOs) confined between a-iron slabs is probed at a range of coverages, pressures, and sliding velocities. At high coverage and low pressure, the nanoparticles do not indent into the a-iron slabs during sliding, leading to zero wear and a low friction coefficient. At low coverage and high pressure, the nanoparticles indent into, and plough through the slabs during sliding, leading to atomic-scale wear and a much higher friction coefficient. This contribution to the friction coefficient is well predicted by an expression developed for macroscopic indentation by Bowden and Tabor. Even at the highest pressures and lowest coverages simulated, both nanoparticles were able to maintain separation of the opposing slabs and reduce friction by approximately 75 % compared to when no nanoparticle was present, which agrees well with experimental observations. CNO nanoparticles yielded a lower indentation (wear) depth and lower friction coefficients at equal coverage and pressure with respect to CND, making them more attractive friction modifier additives. Potential changes in behaviour on harder and softer surfaces are also discussed, together with the implications that these results have in terms of the application of the studied nanoparticles as lubricants additives.

Molecular dynamics simulations of bilayers in a surfactant/co-surfactant/water system with explicit solvent molecules show formation of topologically distinct gel phases depending upon the bilayer composition. At low temperatures, the bilayers transform from the tilted gel phase, Lβ', to the one dimensional (1D) rippled, Pβ' phase as the surfactant concentration is increased. More interestingly, we observe a two dimensional (2D) square phase at higher surfactant concentration which, upon heating, transforms to the gel Lβ' phase. The thickness modulations in the 1D rippled and square phases are asymmetric in two surfactant leaflets and the bilayer thickness varies by a factor of ∼2 between maximum and minimum. The 1D ripple consists of a thinner interdigitated region of smaller extent alternating with a thicker non-interdigitated region. The 2D ripple phase is made up of two superimposed square lattices of maximum and minimum thicknesses with molecules of high tilt forming a square lattice translated from the lattice formed with the thickness minima. Using Voronoi diagrams we analyze the intricate interplay between the area-per-head-group, height modulations and chain tilt for the different ripple symmetries. Our simulations indicate that composition plays an important role in controlling the formation of low temperature gel phase symmetries and rippling accommodates the increased area-per-head-group of the surfactant molecules.

Grafted polymers on the surface of lipid membranes have potential applications in liposome-based drug delivery and supported membrane systems. The effect of polymer grafting on the phase behavior of bilayers made up of single-tail lipids is investigated using dissipative particle dynamics. The bilayer is maintained in a tensionless state using a barostat. Simulations are carried out by varying the grafting fraction, G(f), defined as the ratio of the number of polymer molecules to the number of lipid molecules, and the length of the lipid tails. At low G(f), the bilayer shows a sharp transition from the gel (L(beta)) to the liquid-crystalline (L(alpha)) phase. This main melting transition temperature is lowered as G(f) is increased, and above a critical value of G(f), the interdigitated L(betaI) phase is observed prior to the main transition. The temperature range over which the intermediate phases are observed is a function of the lipid tail length and G(f). At higher grafting fractions, the presence of the L(betaI) phase is attributed to the increase in the area per head group due to the lateral pressure exerted by the polymer brush. The areal expansion and decrease in the melting temperatures as a function of G(f) were found to follow the scalings predicted by the self-consistent mean field theories for grafted polymer membranes. Our study shows that the grafted polymer density can be used to effectively control the temperature range and occurrence of a given bilayer phase.

No abstract