Why liquids can appear to solidify during squeeze-out – Even when they don’t

Combining Molecular Dynamics and Machine Learning to Analyze Shear Thinning for Alkane and Globular Lubricants in the Low Shear Regime

Lubricants with desirable frictional properties are important in achieving an energy-saving society. Lubricants at the interfaces of mechanical components are confined under high shear rates and pressures and behave quite differently from the bulk material. Computational approaches such as nonequilibrium molecular dynamics (NEMD) simulations have been performed to probe the molecular behavior of lubricants. However, the low-shear-velocity regions of the materials have rarely been simulated owing to the expensive calculations necessary to do so, and the molecular dynamics under shear velocities comparable with that in the experiments are not clearly understood. In this study, we performed NEMD simulations of extremely confined lubricants, i.e., two molecular layers for four types of lubricants confined in mica walls, under shear velocities from 0.001 to 1 m/s. While we confirmed shear thinning, the velocity profiles could not show the flow behavior when the shear velocity was much slower than thermal fluctuations. Therefore, we used an unsupervised machine learning approach to detect molecular movements that contribute to shear thinning. First, we extracted the simple features of molecular movements from large amounts of MD data, which were found to correlate with the effective viscosity. Subsequently, the extracted features were interpreted by examining the trajectories contributing to these features. The magnitude of diffusion corresponded to the viscosity, and the location of slips that varied depending on the spherical and chain lubricants was irrelevant. Finally, we attempted to apply a modified Stokes-Einstein relation at equilibrium to the nonequilibrium and confined systems. While systems with low shear rates obeyed the relation sufficiently, large deviations were observed under large shear rates.

History-Dependent Stress Relaxation of Liquids under High-Confinement: A Molecular Dynamics Study

When liquids are confined into a nanometer-scale slit, the induced layering-like film structure allows the liquid to sustain non-isotropic stresses and thus be load-bearing. Such anisotropic characteristics of liquid under confinement arise naturally from the liquids’ wavenumber dependent compressibility, which does not need solidification to take place as a prerequisite. In other words, liquids under confinement can still retain fluidity with molecules being (sub-)diffusive. However, the extensively prolonged structural relaxation times can cause hysteresis of stress relaxation of confined molecules in response to the motions of confining walls and thereby rendering the quasi-static stress tensors history-dependent. In this work, by means of molecular dynamics, stress tensors of a highly confined key base-oil component, i.e., 1-decene trimer, are calculated after its relaxation from being compressed and decompressed. A maximum of 77.1 MPa normal stress discrepancy has been detected within a triple-layer boundary film. Analyses with respect to molecular morphology indicate that among the effects (e.g., confinement, molecular structure, and film density) that can potentially affect confined stresses, the ordering status of the confined molecules plays a predominant role.

Scale-Dependent Friction-Coverage Relations and Nonlocal Dissipation in Surfactant Monolayers

Surfactant molecules, known as organic friction modifiers (OFMs), are routinely added to lubricants to reduce friction and wear between sliding surfaces. In macroscale experiments, friction generally decreases as the coverage of OFM molecules on the sliding surfaces increases; however, recent nanoscale experiments with sharp atomic force microscopy (AFM) tips have shown increasing friction. To elucidate the origin of these opposite trends, we use nonequilibrium molecular dynamics (NEMD) simulations and study kinetic friction between OFM monolayers and an indenting nanoscale asperity. For this purpose, we investigate various coverages of stearamide OFMs on iron oxide surfaces and silica AFM tips with different radii of curvature. We show that the differences between the friction-coverage relations from macroscale and nanoscale experiments are due to molecular plowing in the latter. For our small tip radii, the friction coefficient and indentation depth both have a nonmonotonic dependence on OFM surface coverage, with maxima occurring at intermediate coverage. We rationalize the nonmonotonic relations through a competition of two effects (confinement and packing density) that varying the surface coverage has on the effective stiffness of the OFM monolayers. We also show that kinetic friction is not very sensitive to the sliding velocity in the range studied, indicating that it originates from instabilities. Indeed, we find that friction predominately originates from plowing of the monolayers by the leading edge of the tip, where gauche defects are created, while thermal dissipation is mostly localized in molecules toward the trailing edge of the tip, where the chains return to a more extended conformation.

Computational Simulations of Nanoconfined Argon Film through Adsorption–Desorption in a Uniform Slit Pore

We performed hybrid grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations to investigate the adsorption-desorption isotherms of argon molecules confined between commensurate and incommensurate contacts in nanoscale thickness. The recently proposed mid-density scheme was applied to the obtained hysteresis loops to produce a realistic equilibrium phase of nanoconfined fluids. The appropriate chemical potentials can be determined if the equilibrium structures predicted by GCMC/MD simulations are consistent with those observed in previously developed liquid-vapor molecular dynamics (LVMD) simulations. With the chemical potential as input, the equilibrium structures obtained by GCMC/MD simulations can be used as reasonable initial configurations for future metadynamics free energy calculations.

Experimental and simulation study of the high-pressure behavior of squalane and poly-α-olefins

The equation of state, dynamical properties, and molecular-scale structure of squalane and mixtures of poly-α-olefins at room temperature are studied with a combination of state-of-the-art, high-pressure experiments and molecular-dynamics simulations. Diamond-anvil cell experiments indicate that both materials are non-hydrostatic media at pressures above ∼1 GPa. The equation of state does not exhibit any sign of a first-order phase transition. High-pressure x-ray diffraction experiments on squalane show that there are no Bragg peaks, and hence, the apparent solidification occurs without crystallization. These observations are complemented by a survey of the equation of state and dynamical properties using simulations. The results show that molecular diffusion is essentially arrested above about 1 GPa, which supports the hypothesis that the samples are kinetically trapped in metastable amorphous-solid states. The shear viscosity becomes extremely large at very high pressures, and the coefficient governing its increase from ambient pressure is in good agreement with the available literature data. Finally, simulated radial distribution functions are used to explore the evolution of the molecular-scale structure with increasing pressure. Subtle changes in the short-range real-space correlations are related to a collapse of the molecular conformations with increasing pressure, while the evolution of the static structure factor shows excellent correlation with the available x-ray diffraction data. These results are of indirect relevance to oil-based lubricants, as the pressures involved are comparable to those found in engines, and hence, the ability of lubricating thin films to act as load-bearing media can be linked to the solidification phenomena studied in this work.