Next Article in Journal
Synergistic Effect of Spark Plasma Sintering Driven Solid-Solution Phases on Scratch Resistance in Two-Dimensional Materials
Previous Article in Journal
Wear-Resistant Elastomeric Composites Based on Unvulcanized Rubber Compound and Recycled Polytetrafluoroethylene
Previous Article in Special Issue
Predicting Wear under Boundary Lubrication: A Decisive Statistical Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Coarse-Grained Molecular Dynamics Simulations of Organic Friction Modifier Adsorption on Rough Surfaces under Shear

1
Department of Complex System Science, Graduate School of Informatics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan
2
Faculty of Mechanical Engineering, University Teknologi Malaysia (UTM), Johor Bahru 81310, Malaysia
3
Future Value Creation Research Center (FV-CRC), Graduate School of Informatics, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8601, Japan
*
Author to whom correspondence should be addressed.
Lubricants 2024, 12(2), 30; https://doi.org/10.3390/lubricants12020030
Submission received: 22 December 2023 / Revised: 19 January 2024 / Accepted: 22 January 2024 / Published: 24 January 2024
(This article belongs to the Special Issue Advances in Boundary Lubrication)

Abstract

:
Reducing friction energy losses is crucial in mechanical systems, often achieved through lubrication strategies employing friction modifiers. These additives adsorb onto surfaces, forming boundary film to prevent solid–solid contacts. However, atomistic simulation techniques used to study these additives often ignore surface roughness due to high computational cost. This study addresses this gap by employing Coarse-Grained Molecular Dynamics (CG MD) to investigate the impact of surface roughness on the adsorption of Organic Friction Modifiers (OFMs) under shear. Traditional self-diffusion methods prove inadequate for determining the damping coefficients in CG models because of strong OFM adsorption effects. Therefore, shear-induced motion is introduced for the coefficient determination. The simulation reveals that a symmetrical model (identical opposing surfaces) shows OFM slip, desorption, and re-adsorption trends on rough surfaces, while an asymmetrical model (smooth cylinder on a flat surface) demonstrates increased adsorption on rough flat surfaces (up to 60.9%) compared to smooth flat surfaces under similar shearing conditions. However, rough flat surfaces with a smaller wavelength (6 nm) exhibit faster OFM desorption along the asperity region, up to four times more than a 24 nm wavelength surface. This research emphasizes the importance of considering surface roughness in simulating OFM behavior for lubrication applications.

1. Introduction

Enhancing lubricity for mating surfaces is imperative in mechanical systems to minimize their energy losses due to friction. This is often achieved by incorporating small quantities of friction modifier additives into low-viscosity base oils [1,2,3,4,5]. In contrast to commonly used friction modifiers that contain sulfur and phosphorus, such as Molybdenum Dithiocarbamate (MoDTC) and Molybdenum Dithiophosphate (MoDTP), Organic Friction Modifiers (OFMs) are composed solely of carbon, hydrogen, oxygen, and nitrogen, thereby garnering significant attention and usage in pursuing greener lubrication [2,6,7]. However, OFMs are generally not as effective as molybdenum compounds in reducing friction [8,9]. Therefore, it is crucial to further improve the performance of OFMs.
Typical OFM molecules comprise at least one polar head group and one alkyl chain tail. Previous studies have suggested that OFMs form adsorbed films on sliding surfaces, preventing direct solid–solid contacts, thereby reducing friction in the boundary and mixed lubrication regime [6,10,11,12]. The chemistry of OFMs and solid surfaces determines their interaction strength. Additionally, the morphology of solid surfaces, including surface roughness and textures, determines the pressure distribution in the contact region. Consequently, these factors affect the structural and mechanical properties of OFM adsorption films under shear. A thorough understanding of these correlations is essential for improving the friction reduction performance of OFMs.
Several advanced experimental technologies have been employed to investigate the surface adsorption of OFMs. Atomic Force Microscopy (AFM) has been utilized to measure the morphology and frictional properties of adsorbed films [13,14,15,16,17]. Quartz Crystal Microbalance (QCM) has been adopted to assess not only adsorption mass but also adsorption kinetics [18,19,20,21]. Additionally, Neutron Reflectometry (NR) has been used to analyze the structure of adsorption films in more detail, including the number of layers and the thickness and density of each layer [22,23,24]. These experimental approaches significantly contribute to understanding the effect of OFM chemistry on surface adsorption, primarily under static conditions. Regarding the polar head groups, which govern the surface adsorption strength, studies have indicated that carboxyl and amine groups are more effective than alcohol, ester, and nitrile groups [2,6,20,25]. Additionally, increasing the number of polar groups can enhance adsorption strength through multi-site adsorption or chelate effect [20,26,27,28,29]. Studies have also shown that straight tails, such as saturated or trans-unsaturated alky chains, form more densely packed adsorption films than kink tails, such as cis-unsaturated alky chains, resulting in better friction-reducing performance [15,30]. Nevertheless, measuring adsorbed OFM films on rough surfaces becomes challenging due to their typical thickness of approximately 1 to 2 nm [14,15]. Therefore, there remains a scarcity of research on the correlation between surface roughness and the adsorption efficacy of OFMs. Furthermore, the limitations of experimental technologies become apparent in the lack of thorough analyses concerning the distribution and mechanisms of adsorption and desorption on rough surfaces, especially in the presence of shearing [31].
To address this gap, Molecular Dynamics (MD) simulations are often employed to fundamentally elucidate the adsorption of molecules like friction modifiers at the atomic level. Shi et al. [32] investigated the impact of molecular polarity and temperature on the physisorption behavior of OFMs. Ewen et al. [33] focused on the slip and friction behavior of alkanes on monolayer OFM films adsorbed on solid surfaces during shear motion. Chen et al. [34] recently analyzed the adsorption behavior of N-(2,2,6,6-tetramethyl-1-oxyl-4-piperidinyl)-dodecaneamide (TEMPO-based) OFM on iron oxide surfaces, revealing the mechanism behind its experimentally proven superior lubrication performance. Concerning surface roughness, Eder et al. [35] studied different molecular surface coverages and rough substrates, focusing on the Derjaguin offset. Ewen et al. [36] demonstrated that OFMs can strongly adsorb onto surfaces with random nanoscale roughness (0.2–0.8 nm RMS) to prevent direct contact with solids under high-pressure conditions. Their study highlighted the correlation between nanoscale roughness and stearic acid coverage, indicating lower friction forces in systems with higher coverage.
While the above studies contributed valuable insights, it is crucial to note that the most widely used atomistic simulation techniques have size limitations, where the surface roughness dimensions were either minimized or the system was completely immersed in a lubricating solution [37,38,39]. One significant concern with this approach is the potential oversight of the space for liquid molecule dispersion out of confinement during shearing. This assumption may impact the number of molecules available for adsorption within the confinement.
To address the limitation and simulate larger systems that better mimic realistic contact conditions, we propose using Coarse-Grained (CG) models to simulate the adsorption behavior of OFMs on rough surfaces. CG models allow for the simulation of larger systems with significant surface roughness. However, there remains a challenge in developing CG models that can accurately reproduce the static and dynamic characteristics of target systems [40,41]. Expanding on our previous research [42], we introduce an enhanced approach to replicate the dynamic characteristics of CG models against All-atom (AA) models. This improved methodology is applied to larger-scale rough surface structures, enabling the study of the motion and distribution of OFMs during shearing processes in a more representative manner.

2. Simulation Details

2.1. Simulation Method

In this study, MD simulations were performed using CG models and Dissipative Particle Dynamics (DPD) [43,44,45]. The equation governing the motion of CG beads is expressed as follows:
M i r ¨ i = F i C + F i D + F i R ,
Here, M i and r i represent the mass and position of CG bead i , and F i C , F i D , and F i R denote the conservative, damping, and random forces exerted on CG bead i , respectively. F i C is derived from the negative derivative of interactive potentials between CG beads (CG potentials), and F i D and F i R are calculated as follows:
F i D = j f i j D = j ζ w D r i j e i j · v i j e i j ζ w D r i j v i j e i j · v i j e i j ,
F i R = j f i j R = j σ w R r i j e i j · θ i j e i j + σ w R r i j θ i j e i j · θ i j e i j ,
Here, e i j = r i r j / r i j denotes the unit vector from bead j to i , v i j = v i v j represent the relative velocity of bead i to j , and θ i j is a time-varying three-dimensional vector with normally distributed random numbers as its elements. The superscripts and indicate directions parallel and perpendicular to e i j , respectively. The parameters ζ (referred to as damping coefficient) and σ and weight functions w D and w R serve to adjust the strength of damping and random forces. To achieve the canonical ensemble, the following relationships must hold for both the parallel and perpendicular components:
σ 2 = 2 k B T ζ ,
w R r i j 2 = w D r i j = m a x 0 , 1 r i j / r c D P D 2 ,
where k B is the Boltzmann constant, T is the temperature, and r c D P D is the cutoff distance for the damping and random forces. For simplicity, we used ζ = ζ and w D r i j = w D r i j in this study.
The CG potentials and damping coefficient must be determined a priori. The former affects the static and dynamic properties of simulation systems, whereas the latter solely affects the dynamic properties of simulation systems. Therefore, this study initially derived the CG potentials by matching static properties obtained from AA simulations and then determined the damping coefficient by aligning with dynamic properties from AA simulations.
All AA and CG MD simulations were conducted using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [46], specifically with an in-house extended DPD package for the CG simulations [45]. The simulation results were visualized using the Open Visualization Tool (OVITO) [47].

2.2. Materials and Coarse-Grained Models

Dodecane (C12H26), stearic acid (C18H36O2), and hematite plates with (100) surface orientation were used as base oil, OFM additive, and solid surfaces, respectively. As illustrated in Figure 1, dodecane and stearic acid molecules were coarse-grained by mapping three consecutive carbon atoms and their bonding atoms into a single bead, with the interaction center defined as the center of mass of the constituent atoms. Thus, each dodecane and stearic acid were modeled as four and six beads, respectively, connected by springs. The bead containing the carboxyl group was polar, while the others were nonpolar. For simplicity, the constituent beads of dodecane and stearic acid are collectively referred to as “liquid beads” hereafter. The hematite plates were coarse-grained by arranging solid beads in a simple tetragonal lattice with lattice constants a = b = 5.714 Å and c = 5.814 Å (see Figure 1).
As depicted in Figure 2, CG potentials included bond–length potentials ( U b l n p and U b l p ), bond–angle potentials ( U b a n p and U b a p ), nonbonded liquid–liquid interaction potentials ( U l l n p and U l l p ), and nonbonded liquid–solid interaction potentials ( U l s n p and U l s p ). The subscripts “np” and “p” indicate involvement without and with polar beads, respectively, with U l l n p being the only exception, which was also applied to polar and nonpolar liquid bead pairs. Notably, U l l p and U l l n p were applied to liquid–bead pairs not involved in the bond–length and bond–angle interactions. These potentials were derived through the iterative Boltzmann inversion method, aligning with reference distribution functions characterizing static structural properties obtained from AA simulations, the same as in our previous studies [42,48]. Specifically, our previous work [42] extensively covered the derivation and validation of U l l n p , U b l n p , and U b a n p for dodecane, demonstrating their high accuracy and transferability across a wide pressure range of up to 1 GPa. The remaining CG potentials, U b l p , U b a p , U l l p , U l s p , and U l s n p , alongside their validation, are detailed in Appendix A.

2.3. Determination of Damping Coefficients

The damping coefficients included ζ l l p , ζ l l n p , ζ l s p , and ζ l s n p for bead pairs of polar liquid, nonpolar liquid, polar liquid–solid, and nonpolar liquid–solid. Note that ζ l l n p was also used for bead pairs of polar and nonpolar liquids. Conventionally, damping coefficients are adjusted by matching the self-diffusion coefficient of thermal motion obtained from AA simulations. However, due to the strong interactions of the polar end groups, stearic acid molecules form dimers in the bulk and adsorb on solid surfaces [49,50]. Therefore, thermally induced motion rarely occurs between the polar end groups and between them and solid surfaces, making it challenging to adjust ζ l l p , and ζ l s p by the conventional method, as evident from Equation (2). In this study, we addressed this challenge by employing shear-induced motion rather than thermal motion [51].
In brief, ζ l l n p was initially determined by matching the shear viscosities from non-equilibrium AA and CG MD simulations of bulk dodecane systems under shear. Subsequently, ζ l l p was determined similarly using bulk stearic acid systems and the previously determined ζ l l n p . For adjustment of ζ l s n p and ζ l s p , we matched the Coefficients of Friction (CoF) from non-equilibrium AA and CG MD simulations, using a single layer of dodecane and stearic acid, respectively, lying flat between two hematite (100) surfaces, and applying a normal pressure and moving the upper solid surface horizontally. Appendix B provides details of the simulation setups and calculation of viscosities and CoF.
Figure 3 and Figure 4 show the variation of bulk shear viscosities and CoF with the liquid–liquid and liquid–solid damping coefficients in the CG simulations, respectively. The values obtained from the parent AA simulations are also plotted in the figures with horizontal lines for comparison. These shear-related properties increase monotonically with the damping coefficients, showcasing their effectiveness in fine-tuning the damping coefficients. Table 1 summarizes the final values of ζ l l n p , ζ l l p , ζ l s n p , and ζ l s p , along with the corresponding shear-related properties from the CG and AA simulations. The good agreement validates the high accuracy of our CG models, even for reproducing dynamic properties.

2.4. Models and Conditions of CG Shear Simulations

This study used two types of models, symmetric and asymmetric, for the CG shear simulations, as shown in Figure 5 and Figure 6. The symmetric model had identical smooth or rough solid surfaces on the top and bottom, while the asymmetric model had a smooth cylindrical solid surface at the top and a smooth or rough flat surface at the bottom. In the symmetric model, the entire region between the solid surfaces was pressurized, whereas, in the asymmetric model, the pressure was only present in the region below the flat apex (24 nm in length) of the upper cylindrical surface, allowing liquid films to disperse from the pressurized region to surrounding spaces.
Both models had dimensions of 48 and 4 nm in the x and y directions, respectively, with periodic boundary conditions applied. A one-dimensional sinusoidal function z = A sin 2 π x / λ was used to describe the roughness of the upper and lower solid surfaces in the symmetric model and the lower solid surface in the asymmetric model. To capture the essence of surface roughness effect at reasonable computational costs, the amplitude A was fixed at 1.7 nm, and the wavelength λ was varied to 24, 48 nm, and infinity in the symmetric model, and 6, 12, 24, 48 nm, and infinity in the asymmetric model. It is to note that λ = means smooth surfaces. The radius of curvature of the cylindrical upper surface was 40 nm.
A monolayer physically adsorbed film was formed by placing the polar liquid beads of 168 stearic acid molecules on each solid surface, with the molecular chains perpendicular to the surface. This corresponded to an areal density of 0.875 molecules per nm2. Then, a mixed solution of 300 dodecane and 300 stearic acid molecules was inserted into the central region between the two solid surfaces. The liquid films were then compressed by applying a constant pressure of 0.5 GPa on the upper solid surface for 0.5 ns. The shear simulation was performed by moving the upper solid surface at 10.0 m/s in the x direction while keeping the lower solid surface stationary. The system temperature was maintained at 353.15 K, and the simulation timestep was set to 5.0 fs. Movie S1 in the Supplementary Materials demonstrates the shear motion in the asymmetric model with a wavelength of 6 nm.

3. Results and Discussion

3.1. Symmetric Models

3.1.1. Adsorption Behavior

In the symmetrical model, Figure 7a depicts the time-dependent changes in OFM adsorption during the shearing of identical opposing surfaces. The adsorption criterion is that the distance between the OFM polar beads and their nearest solid beads is less than 0.7 nm, corresponding to the equilibrium distance of interaction between the polar liquid bead and the solid bead ( U l s p ). The adsorption fraction is calculated using n(t)/N, where n(t) represents the number of adsorbed OFMs at time t, and N = 636 is the total number of OFMs. It is noted that the number of adsorbed OFM molecules presented is only along the bottom surface. Figure 7a reveals that the adsorption of OFM reaches a steady state after 20 ns; however, contrasting the stable saturation on the smooth surfaces, fluctuations in OFM adsorption are observed on the rough surfaces. The mean adsorption fraction is calculated from the final 30 ns of the simulations, as shown in Figure 7b. The adsorption fraction decreases slightly with increasing roughness wavelength, by 3.2% (24 nm) and 6.2% (48 nm) compared with the smooth surfaces. The error bars indicate that the variations in adsorption fraction on the rough surfaces are approximately twice as large as those on the smooth surfaces, emphasizing higher fluctuations on rough surfaces. These findings suggest that, in comparison to the stable adsorption of OFMs on smooth surfaces, as also confirmed in previous simulation studies [7,32], the adsorption of OFMs on rough surfaces is less stable.

3.1.2. Interpretation of Fluctuating OFM Adsorption Behavior

Figure 8a examines the fluctuating adsorption behavior of the OFM molecules on the rough surfaces in the symmetric model. Although the fluctuation period for the surfaces with a 48 nm wavelength is half that for the surfaces with a 24 nm wavelength, the reasons for the fluctuations are the same, as explained below. Therefore, the subsequent discussion will focus on the 24 nm wavelength. As observed at 18.8, 19.4, and 19.9 ns, there are three significant states of the fluctuations: local maximum adsorbed OFM molecules, transitioning state, and local minimum adsorbed OFM molecules. Snapshots of these states are provided in Figure 8b–d.
To comprehend the movement of OFMs during shearing, attention is focused on a cluster of polar beads (depicted in red) already adhered to the convex areas of the lower corrugated rough surface in Figure 8b. Under sufficiently high pressure during shearing, this cluster of OFMs (in red) slips from the convex region as opposing peaks align, as in Figure 8c. As the shearing progresses, these OFM molecules slide away from the convex region, detaching and migrating to the concave region. Concurrently, another cluster of OFMs (depicted in blue) is sheared onto the vacant sites along the convex region. In Figure 8d, these OFM molecules (in blue) re-adsorb onto the convex region of the opposing surface. Meanwhile, the disengaged molecules (in red) are either carried away until they coincide with another convex region or re-adsorb along the transitional slope region of the surface feature. This slip, desorption, and re-adsorption cycle persists throughout the simulation, contributing to the fluctuating adsorption trend illustrated in Figure 7a. The cyclic behavior is attributed to the identical opposing rough surfaces (i.e., roughness amplitude and wavelength), potentially resulting in a less stable layer of adsorbed OFM.

3.2. Asymmetric Models

3.2.1. Adsorption Behavior

In the asymmetrical model, Figure 9a illustrates the variation in OFM adsorption over a 300 ns simulation as a smooth cylinder slides across a smooth, flat surface or a flat surface, with sinusoidal roughness wavelengths ranging from 6 to 48 nm. On the smooth surface, the plot indicates an initial increase followed by a gradual decrease in adsorbed OFM, finally leveling off into a steady state. This initial rise is attributed to the wedge effect induced by the cylindrical indenter, drawing in additional molecules, including those not initially adsorbed. Some of these molecules may adhere to available sites, increasing adsorption coverage within the confinement. However, sufficiently high pressure and shear forces could also lead to the slip and subsequent removal of adsorbed OFMs within the confined space. This cyclic process continues until a steady state is reached, typically around 200 ns in the simulation.
Similarly, simulations involving flat bottom surfaces with varying roughness reach a steady state at around 200 ns, displaying a consistent increase before saturation. The total number of adsorbed OFM molecules increases with smaller roughness wavelength values. Comparing this behavior with sliding on the smooth flat surface, the contact experiences an increase in mean OFM adsorption sites ranging between 9.9% (48 nm) and 60.9% (6 nm), as depicted in Figure 9b. The impact of surface roughness in the asymmetric model is clearly more pronounced than in the symmetric model shown in Figure 7b. Notably, the curvature of the indenter closely resembles the surface asperity curvature with a 48 nm wavelength, leading to a similar fluctuation as observed in Figure 7a. This behavior reflects the cyclic slip, desorption, and re-adsorption pattern of OFM molecules observed previously. The fluctuations diminish when the curvature of the indenter deviates from the surface asperity curvature.
The study further investigates specific locations on the rough surface where OFM molecules adhere, including the concave, convex, and transitional slope regions within the confinement, as shown in Figure 10. For all the rough surfaces, the adsorption fraction is the highest on the slope region. This is because that the total area of slope regions is the largest, providing more adsorption sites than the convex and concave regions. Moreover, due to the corner effect, more solid beads are in close contact with the polar beads adsorbed on the slope regions compared to the other regions, thereby resulting in stronger adsorption in the slope regions.
With decreasing wavelength, i.e., increasing surface roughness, the adsorption fraction increases on the slope and convex regions but decreases on the concave regions. Particularly, the change induced by surface roughness is more pronounced in the slope regions than in the convex and concave regions, which dictates the trend of increasing adsorption fraction with decreasing wavelength as observed in Figure 9b. The reasons are discussed below, based on the motions of OFMs as shown in Movie S2 in the Supplementary Materials. Due to the shear motion of the upper surface in the horizontal direction, OFMs possess horizontal rather than vertical velocity. As seen in Movie S2, OFM molecules not adsorbed onto the convex regions tend to traverse horizontally in the shallower cavity regions instead of diffusing vertically into the deeper cavity regions. If there are unoccupied slopes in the shallower regions, OFMs are likely to adsorb onto these slopes; otherwise, they move to the next convex region. As the wavelength decreases, although the total area of the slopes remains nearly constant, the area of slopes in the shallower region increases. This provides more favorable adsorption sites, resulting in a higher number of adsorbed OFMs on the slope regions. Additionally, with decreasing wavelength, the distance and travelling time of OFMs between neighboring convexes decrease. The shortened time make it more difficult for OFMs to diffuse to the deeper cavity regions, thereby leading to a decreased adsorption in the concave regions. Moreover, compared to OFMs in larger cavities of surfaces with longer wavelengths, OFMs in smaller cavities of surfaces with shorter wavelengths are more likely to entrain into the contacting interface, consequently resulting in an increased adsorption fraction in the convex regions.

3.2.2. Interpretation of OFM Desorption in Shearing

Focusing on the convex region, where boundary friction is expected to occur, the desorption time (td) for OFMs on a convex region is calculated as td = tend − t0, where tend is the time at which the last initially adsorbed OFM leaves the original convex region, and t0 is the time at which shearing commenced. Figure 11 presents the desorption time, indicating the duration needed for OFM, initially pre-adsorbed in convex regions, to completely desorb from their original locations. The plot illustrates a noticeable increase in the desorption time with larger roughness wavelength values. The variations in desorption time can be attributed to the increased number of surface peaks or asperities, with a smaller wavelength value correlating to an expedited OFM desorption process.
To further interpret OFM desorption time, the current study employs the Greenwood and Williamson (GW) rough surface contact model [52]. The model assumes only one rough surface, similar to the asymmetric model simulated in the present study. Following this model, predictions for asperity contact area and contact load for the asymmetrical model, between the flat apex of the upper cylindrical surface and the rough flat surface, are made based on the following equations [52]:
T o t a l   a s p e r i t y   c o n t a c t   a r e a m 2 : A a = π ξ β σ F 1 h / σ A ,
T o t a l   a s p e r i t y   c o n t a c t   l o a d   ( N ) : P a = 4 3 ξ β 1 / 2 σ 3 / 2 F 3 / 2 h / σ E A ,
where E , A , and h refers to the reduced modulus of elasticity, apparent contact area, and gap between the upper and lower surfaces, in the units of Pa, m2, and m, respectively. Equations (6) and (7) highlight that the primary roughness geometrical factors influencing asperity area and load are ξ β σ , and ξ β 1 / 2 σ 3 / 2 , respectively.
Utilizing the model, the rough surfaces generated with different wavelength values are translated into GW parameters, as detailed in Table 2. Since the amplitude remains constant, the variations in the adopted configurations of rough surfaces result in the same RMS roughness of 1.204 nm. The distinction lies in the mean asperity peak radius ( β ) and density ( ξ ). Following the GW parameters, a surface with a smaller wavelength value yields higher asperity density and smaller-sized asperities (based on curvature radius).
Table 2 reveals a consistent total asperity contact area across the various generated rough surfaces, while the load supported by surface asperities exhibits a significant difference. Surfaces with a 6 nm wavelength value are expected to carry an asperity load up to three times larger than those with a 24 nm wavelength value. The faster desorption rate observed for smaller wavelength values in Figure 11 is likely attributed to the increased load on surface asperities or the convex region, resulting in quicker desorption of OFM molecules. The faster desorption could also encourage more OFM molecules to migrate to the transitional slope region during shear, further increasing the adsorbed molecules in this region, as shown in Figure 10.
Additionally, Table 2 highlights the association between smaller wavelength values and higher mean asperity slopes. Sliding on such slopes often aligns with Euler’s interlocking mechanism, where shear is proportionate to the tangent of the asperity slope. In scenarios involving potential solid wear particles (of similar length scale as the surface roughness) in the concave region, the surplus of adsorbed OFM molecules along the transitional slope between the concave and convex regions, particularly with smaller wavelength values, offers a more convenient pathway for these particles to slide out of the concave region with greater ease. Consequently, an enhanced number of adsorbed OFMs in this region is a mitigating factor, reducing the potential for erosive wear.

4. Conclusions

In this study, CG MD was employed to model OFM adsorption under shear in two distinct configurations: (1) symmetric model (identical opposing surfaces) and (2) asymmetric model (smooth cylinder on flat surface). The construction of the CG model faced challenges due to the strong effect of OFM adsorption and dimer cluster, making the traditional self-diffusion methods unsuitable for determining the damping coefficients of DPD. To address this challenge, shear-induced motion was employed, allowing the CG model to accurately match the viscosity of the AA model in liquid–liquid situations and the CoF of the AA model in liquid–solid situations. This approach yielded a more appropriate determination of the damping coefficients, improving the relevance of the CG model in simulating OFM behavior under shear conditions.
The symmetric model, considering identical opposing surfaces in shear, reveals that, despite a higher number of adsorbed OFMs by up to 6.2%, a distinct cyclic slip, desorption, and re-adsorption process occurs when opposing surfaces share similar roughness wavelength and amplitude during shear. This cyclic behavior, absent in the atomically smooth flat-on-flat configuration, emphasizes the complex nature of OFM adsorption behavior under shear conditions in the presence of surface roughness.
Simulations with a more realistic contact scenario with the asymmetric model involving a smooth cylinder sliding on a rough flat surface also show an overall increase in adsorbed OFMs, reaching a substantial 60.9%. This increase is observed primarily in a rough flat surface with a 6 nm wavelength compared to sliding on a smooth flat surface. However, it is important to note that this configuration’s surplus of adsorbed molecules mainly accumulates along the transitional slope between convex and concave regions. This rough flat surface also exhibits faster OFM desorption along the convex or asperity region, attributed to the larger load carried by the asperity along the convex regions, reaching up to four times compared to a rough surface with a 24 nm wavelength.
In summary, this study observed the varying OFM adsorption and desorption behavior under diverse shear conditions, particularly concerning surface roughness. Unlike observations on atomically smooth surfaces, the cyclic patterns of OFM desorption and re-adsorption between rough surfaces, alongside enhanced adsorption in realistic contact scenarios, underscore the pivotal role of surface roughness in dictating the dynamic behavior of OFM adsorption under shear.
The study highlights the potential ramifications of overlooking surface roughness when simulating or designing OFM molecular structures. Such oversights may compromise the reliability of friction-reducing performance. Striking a balance between simulation details and complexity while accurately modeling rough surfaces is essential. This ensures a comprehensive understanding of OFM interactions under varying shear conditions, facilitating reliable predictions and informed design choices for OFMs in lubrication applications.

Supplementary Materials

The following supplementary information can be downloaded at: https://www.mdpi.com/article/10.3390/lubricants12020030/s1, Movie S1: Shear motion in the asymmetric model with a wavelength of 6 nm. Movie S2: Representative motions of OFMs observed in the asymmetric model with a wavelength of 12 nm. The other OFM and dodecane molecules are omitted for clarity.

Author Contributions

Conceptualization, J.T., W.W.F.C. and H.Z.; methodology, J.T. and H.Z.; software, J.T.; validation, J.T.; formal analysis, J.T.; investigation, J.T., W.W.F.C. and H.Z.; resources, J.T. and H.Z.; data curation, J.T.; writing—original draft preparation, J.T. and W.W.F.C.; writing—review and editing, J.T., W.W.F.C. and H.Z.; visualization, J.T.; supervision, W.W.F.C. and H.Z.; funding acquisition, H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by JSPS KAKENHI Grants (Nos. 21H01238 and 23H05448).

Data Availability Statement

The data presented in this study are available in the article.

Acknowledgments

J.T. acknowledges the financial support provided by the THERS Interdisciplinary Frontier Next Generation Researcher Program (No. RD191018).

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Coarse-Grained Potentials and Validation

Appendix A.1. Liquid–Liquid Potentials

The figures below present the bonded potentials for bond–length ( U b l p ) and bond–angle ( U b a p ), as well as the nonbonded potential ( U l l p ) for polar beads, alongside their corresponding distribution functions. In all the distribution functions, the coarse-grained (CG) results reasonably closely match the all-atom (AA) results, validating the accuracy of the corresponding CG potentials.
Figure A1. (a) Coarse-grained bond-length potential involving polar bead, U b l p . (b) Bond length distribution functions derived from simulations of an isolated stearic acid using all-atom and our coarse-grained models.
Figure A1. (a) Coarse-grained bond-length potential involving polar bead, U b l p . (b) Bond length distribution functions derived from simulations of an isolated stearic acid using all-atom and our coarse-grained models.
Lubricants 12 00030 g0a1
Figure A2. (a) Coarse-grained bond-angle potential involving polar bead, U b a p . (b) Bond angle distribution functions derived from simulations of an isolated stearic acid using all-atom and our coarse-grained models.
Figure A2. (a) Coarse-grained bond-angle potential involving polar bead, U b a p . (b) Bond angle distribution functions derived from simulations of an isolated stearic acid using all-atom and our coarse-grained models.
Lubricants 12 00030 g0a2
Figure A3. (a) Coarse-grained nonbonded polar liquid–liquid potential, U l l p . (b) Radial distribution functions of polar beads derived from simulation of bulk stearic acid systems using all-atom and our coarse-grained models.
Figure A3. (a) Coarse-grained nonbonded polar liquid–liquid potential, U l l p . (b) Radial distribution functions of polar beads derived from simulation of bulk stearic acid systems using all-atom and our coarse-grained models.
Lubricants 12 00030 g0a3

Appendix A.2. Liquid–Solid Potentials

The figures below present the nonpolar and polar liquid–solid potentials, U l s n p and U l s p , alongside the number density distributions (NDDs) of nonpolar and polar liquid beads versus their distance from the solid surface. In Figure A4b, the CG result strongly aligns with the AA one. Similarly, Figure A5b demonstrates a close match between the first peak in the CG and AA results. These validate the effectiveness of the CG liquid–solid potentials in reproducing the structures of dodecane and adsorbed stearic acid films on hematite surfaces.
Figure A4. (a) Coarse-grained nonbonded nonpolar liquid–solid potential, U l s n p . (b) Number density distribution functions of nonpolar liquid beads derived from simulation of a dodecane film on a hematite plate using all-atom and our coarse-grained models.
Figure A4. (a) Coarse-grained nonbonded nonpolar liquid–solid potential, U l s n p . (b) Number density distribution functions of nonpolar liquid beads derived from simulation of a dodecane film on a hematite plate using all-atom and our coarse-grained models.
Lubricants 12 00030 g0a4
Figure A5. (a) Coarse-grained nonbonded polar liquid–solid potential, U l s p . (b) Number density distribution functions of polar liquid beads derived from simulation of a stearic acid film on a hematite plate using all-atom and our coarse-grained models.
Figure A5. (a) Coarse-grained nonbonded polar liquid–solid potential, U l s p . (b) Number density distribution functions of polar liquid beads derived from simulation of a stearic acid film on a hematite plate using all-atom and our coarse-grained models.
Lubricants 12 00030 g0a5

Appendix B. Simulation Setup for Fine-Tuning Damping Coefficients

Appendix B.1. Liquid–Liquid Damping Coefficients

As described in the main text, the liquid–liquid damping coefficients ζ l l n p and ζ l l p were determined by matching the shear viscosities from non-equilibrium AA and CG simulations of bulk dodecane and stearic acid systems under shear, respectively. All simulations maintained a temperature of 353.15 K. The bulk dodecane system comprised 405 molecules in a box of 6.25 nm × 6.25 nm × 4.17 nm in the x, y, and z directions, yielding a density of 0.69 g/cm3 at 0.1 MPa, which was the same as in our previous study [42]. The bulk stearic acid system contained 283 molecules in a box of 6.26 nm × 6.26 nm × 4.17 nm in the x, y, and z directions, producing a density of 0.83 g/cm3, which approximates the experimental value under identical temperature and pressure conditions [55]. Periodic boundary conditions were applied to all the three directions.
To derive bulk viscosities, a constant engineering strain rate was applied to deform the simulation boxes and the motion of the systems was calculated using the SLLOD algorithm [51]. As depicted in Figure A6, shear was applied to the x–z plane in the x direction at a velocity of 1 m/s, generating a velocity gradient in the y direction. Viscosities were calculated using the equation η = P x y / s , where P x y is the pressure tensor and s   ( = v x / l y ) is the shear rate with v x being the shear velocity and l y being the box length in the y direction. Each shear simulation spanned 3.1 ns. The individual viscosity values shown in Figure 3 and Table 1 in the main text are the average over the last 1 ns of three independent simulations.
Figure A6. Snapshot from coarse-grained shear simulation of bulk stearic acid. The green and light blue represent the polar and nonpolar beads, respectively.
Figure A6. Snapshot from coarse-grained shear simulation of bulk stearic acid. The green and light blue represent the polar and nonpolar beads, respectively.
Lubricants 12 00030 g0a6

Appendix B.2. Liquid–Solid Damping Coefficients

As described in the main text, the liquid-solid damping coefficients ζ l s n p and ζ l s p were determined by matching the coefficients of friction (CoF) from non-equilibrium AA and CG simulations, which involved the shearing of a monolayer layer of dodecane and stearic acid, respectively, between two hematite substrates. As depicted in Figure A7, the solid substrates measured 4.03 nm and 4.12 nm in the x and y directions, respectively, with periodic boundary conditions applied. The confined liquid film comprised 50 dodecane or 30 stearic acid molecules.
A normal pressure of 0.5 GPa and a shear velocity of 10 m/s in the x direction were applied to the top layer of the upper substrate, while the bottom layer of the lower substrate remained fixed. CoFs were computed by dividing the shear stress in the x direction of the x-y plane by the normal pressure. Each shear simulation spanned 5 ns. The individual CoF values shown in Figure 4 and Table 1 in the main text are the average over the last 2 ns of three independent simulations.
Figure A7. Snapshot from coarse-grained shear simulation of confined monolayer stearic acid film between solid substrates. The green, light blue, and gray represent the polar and nonpolar liquid beads and solid beads, respectively.
Figure A7. Snapshot from coarse-grained shear simulation of confined monolayer stearic acid film between solid substrates. The green, light blue, and gray represent the polar and nonpolar liquid beads and solid beads, respectively.
Lubricants 12 00030 g0a7

References

  1. Holmberg, K.; Erdemir, A. Influence of Tribology on Global Energy Consumption, Costs and Emissions. Friction 2017, 5, 263–284. [Google Scholar] [CrossRef]
  2. Spikes, H. Friction Modifier Additives. Tribol. Lett. 2015, 60, 5. [Google Scholar] [CrossRef]
  3. McQueen, J.S.; Gao, H.; Black, E.D.; Gangopadhyay, A.K.; Jensen, R.K. Friction and Wear of Tribofilms Formed by Zinc Dialkyl Dithiophosphate Antiwear Additive in Low Viscosity Engine Oils. Tribol. Int. 2005, 38, 289–297. [Google Scholar] [CrossRef]
  4. Zhang, J.; Meng, Y. Boundary Lubrication by Adsorption Film. Friction 2015, 3, 115–147. [Google Scholar] [CrossRef]
  5. Vaitkunaite, G.; Espejo, C.; Wang, C.; Thiébaut, B.; Charrin, C.; Neville, A.; Morina, A. MoS2 Tribofilm Distribution from Low Viscosity Lubricants and Its Effect on Friction. Tribol. Int. 2020, 151, 106531. [Google Scholar] [CrossRef]
  6. Cyriac, F.; Tee, X.Y.; Poornachary, S.K.; Chow, P.S. Influence of Structural Factors on the Tribological Performance of Organic Friction Modifiers. Friction 2021, 9, 380–400. [Google Scholar] [CrossRef]
  7. Ewen, J.P.; Gattinoni, C.; Morgan, N.; Spikes, H.A.; Dini, D. Nonequilibrium Molecular Dynamics Simulations of Organic Friction Modifiers Adsorbed on Iron Oxide Surfaces. Langmuir 2016, 32, 4450–4463. [Google Scholar] [CrossRef]
  8. Wang, W.; Li, C.; Yang, J.; Shen, Y.; Xu, J. Friction Performance of MoDTP and Ester-containing Lubricants between CKS Piston Ring and Cast Iron Cylinder Liner. Lubr. Sci. 2018, 30, 33–43. [Google Scholar] [CrossRef]
  9. Cyriac, F.; Yi, T.X.; Poornachary, S.K.; Chow, P.S. Boundary Lubrication Performance of Polymeric and Organic Friction Modifiers in the Presence of an Anti-Wear Additive. Tribol. Int. 2022, 165, 107256. [Google Scholar] [CrossRef]
  10. Tang, Z.; Li, S. A Review of Recent Developments of Friction Modifiers for Liquid Lubricants (2007–Present). Curr. Opin. Solid State Mater. Sci. 2014, 18, 119–139. [Google Scholar] [CrossRef]
  11. Cañellas, G.; Emeric, A.; Combarros, M.; Navarro, A.; Beltran, L.; Vilaseca, M.; Vives, J. Tribological Performance of Esters, Friction Modifier and Antiwear Additives for Electric Vehicle Applications. Lubricants 2023, 11, 109. [Google Scholar] [CrossRef]
  12. Shi, J.; Li, H.; Lu, Y.; Sun, L.; Xu, S.; Fan, X. Synergistic Lubrication of Organic Friction Modifiers in Boundary Lubrication Regime by Molecular Dynamics Simulations. Appl. Surf. Sci. 2023, 623, 157087. [Google Scholar] [CrossRef]
  13. Hou, J.; Tsukamoto, M.; Zhang, H.; Fukuzawa, K.; Itoh, S.; Azuma, N. Characterization of Organic Friction Modifiers Using Lateral Force Microscopy and Eyring Activation Energy Model. Tribol. Int. 2023, 178, 108052. [Google Scholar] [CrossRef]
  14. Campen, S.; Green, J.H.; Lamb, G.D.; Spikes, H.A. In Situ Study of Model Organic Friction Modifiers Using Liquid Cell AFM: Self-Assembly of Octadecylamine. Tribol. Lett. 2015, 58, 39. [Google Scholar] [CrossRef]
  15. Campen, S.; Green, J.H.; Lamb, G.D.; Spikes, H.A. In Situ Study of Model Organic Friction Modifiers Using Liquid Cell AFM: Saturated and Mono-Unsaturated Carboxylic Acids. Tribol. Lett. 2015, 57, 18. [Google Scholar] [CrossRef]
  16. Sahoo, R.R.; Biswas, S.K. Frictional Response of Fatty Acids on Steel. J. Colloid Interf. Sci. 2009, 333, 707–718. [Google Scholar] [CrossRef]
  17. Hamdan, S.H.; Lee, C.T.; Lee, M.B.; Chong, W.W.F.; Chong, C.T.; Sanip, S.M. Synergistic Nano-Tribological Interaction between Zinc Dialkyldithiophosphate (ZDDP) and Methyl Oleate for Biodiesel-Fueled Engines. Friction 2021, 9, 612–626. [Google Scholar] [CrossRef]
  18. Zachariah, Z.; Nalam, P.C.; Ravindra, A.; Raju, A.; Mohanlal, A.; Wang, K.; Castillo, R.V.; Espinosa-Marzal, R.M. Correlation Between the Adsorption and the Nanotribological Performance of Fatty Acid-Based Organic Friction Modifiers on Stainless Steel. Tribol. Lett. 2020, 68, 11. [Google Scholar] [CrossRef]
  19. Gmür, T.A.; Mandal, J.; Cayer-Barrioz, J.; Spencer, N.D. Towards a Polymer-Brush-Based Friction Modifier for Oil. Tribol. Lett. 2021, 69, 124. [Google Scholar] [CrossRef]
  20. Nalam, P.C.; Pham, A.; Castillo, R.V.; Espinosa-Marzal, R.M. Adsorption Behavior and Nanotribology of Amine-Based Friction Modifiers on Steel Surfaces. J. Phys. Chem. C 2019, 123, 13672–13680. [Google Scholar] [CrossRef]
  21. Fry, B.M.; Moody, G.; Spikes, H.A.; Wong, J.S.S. Adsorption of Organic Friction Modifier Additives. Langmuir 2020, 36, 1147–1155. [Google Scholar] [CrossRef] [PubMed]
  22. Shen, W.; Hirayama, T.; Yamashita, N.; Adachi, M.; Oshio, T.; Tsuneoka, H.; Tagawa, K.; Yagishita, K.; Yamada, N.L. Relationship between Interfacial Adsorption of Additive Molecules and Reduction of Friction Coefficient in the Organic Friction Modifiers-ZDDP Combinations. Tribol. Int. 2022, 167, 107365. [Google Scholar] [CrossRef]
  23. Cyriac, F.; Yamashita, N.; Hirayama, T.; Yi, T.X.; Poornachary, S.K.; Chow, P.S. Mechanistic Insights into the Effect of Structural Factors on Film Formation and Tribological Performance of Organic Friction Modifiers. Tribol. Int. 2021, 164, 107243. [Google Scholar] [CrossRef]
  24. Yamashita, N.; Hirayama, T.; Hino, M.; Yamada, N.L. Neutron Reflectometry under High Shear in Narrow Gap for Tribology Study. Sci. Rep. 2023, 13, 18268. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, X.; Tsukamoto, M.; Zhang, H.; Mitsuya, Y.; Itoh, S.; Fukuzawa, K. Experimental Study of Application of Molecules with a Cyclic Head Group Containing a Free Radical as Organic Friction Modifiers. J. Adv. Mech. Des. Syst. Manuf. 2020, 14, JAMDSM0044. [Google Scholar] [CrossRef]
  26. Hou, J.; Tsukamoto, M.; Hor, S.; Chen, X.; Yang, J.; Zhang, H.; Koga, N.; Yasuda, K.; Fukuzawa, K.; Itoh, S.; et al. Molecules with a TEMPO-Based Head Group as High-Performance Organic Friction Modifiers. Friction 2023, 11, 316–332. [Google Scholar] [CrossRef]
  27. Hu, W.; Xu, Y.; Zeng, X.; Li, J. Alkyl-Ethylene Amines as Effective Organic Friction Modifiers for the Boundary Lubrication Regime. Langmuir 2020, 36, 6716–6727. [Google Scholar] [CrossRef]
  28. Desanker, M.; He, X.; Lu, J.; Liu, P.; Pickens, D.B.; Delferro, M.; Marks, T.J.; Chung, Y.-W.; Wang, Q.J. Alkyl-Cyclens as Effective Sulfur- and Phosphorus-Free Friction Modifiers for Boundary Lubrication. ACS Appl. Mater. Inter. 2017, 9, 9118–9125. [Google Scholar] [CrossRef]
  29. Desanker, M.; He, X.; Lu, J.; Johnson, B.A.; Liu, Z.; Delferro, M.; Ren, N.; Lockwood, F.E.; Greco, A.; Erdemir, A.; et al. High-Performance Heterocyclic Friction Modifiers for Boundary Lubrication. Tribol. Lett. 2018, 66, 50. [Google Scholar] [CrossRef]
  30. Ouyang, C.; Bai, P.; Wen, X.; Zhang, X.; Meng, Y.; Ma, L.; Tian, Y. Effects of Conformational Entropy on Antiwear Performances of Organic Friction Modifiers. Tribol. Int. 2021, 156, 106848. [Google Scholar] [CrossRef]
  31. Jacobs, T.D.B.; Ryan, K.E.; Keating, P.L.; Grierson, D.S.; Lefever, J.A.; Turner, K.T.; Harrison, J.A.; Carpick, R.W. The Effect of Atomic-Scale Roughness on the Adhesion of Nanoscale Asperities: A Combined Simulation and Experimental Investigation. Tribol. Lett. 2013, 50, 81–93. [Google Scholar] [CrossRef]
  32. Shi, J.; Zhou, Q.; Sun, K.; Liu, G.; Zhou, F. Understanding Adsorption Behaviors of Organic Friction Modifiers on Hydroxylated SiO2 (001) Surfaces: Effects of Molecular Polarity and Temperature. Langmuir 2020, 36, 8543–8553. [Google Scholar] [CrossRef] [PubMed]
  33. Ewen, J.P.; Kannam, S.K.; Todd, B.D.; Dini, D. Slip of Alkanes Confined between Surfactant Monolayers Adsorbed on Solid Surfaces. Langmuir 2018, 34, 3864–3873. [Google Scholar] [CrossRef] [PubMed]
  34. Chen, X.; Yang, J.; Yasuda, K.; Koga, N.; Zhang, H. Adsorption Behavior of TEMPO-Based Organic Friction Modifiers during Sliding between Iron Oxide Surfaces: A Molecular Dynamics Study. Langmuir 2022, 38, 3170–3179. [Google Scholar] [CrossRef] [PubMed]
  35. Eder, S.J.; Vernes, A.; Betz, G. On the Derjaguin Offset in Boundary-Lubricated Nanotribological Systems. Langmuir 2013, 29, 13760–13772. [Google Scholar] [CrossRef] [PubMed]
  36. Ewen, J.P.; Echeverri Restrepo, S.; Morgan, N.; Dini, D. Nonequilibrium Molecular Dynamics Simulations of Stearic Acid Adsorbed on Iron Surfaces with Nanoscale Roughness. Tribol. Int. 2017, 107, 264–273. [Google Scholar] [CrossRef]
  37. Gao, J.; Luedtke, W.D.; Landman, U. Structures, Solvation Forces and Shear of Molecular Films in a Rough Nano-Confinement. Tribol. Lett. 2000, 9, 3–13. [Google Scholar] [CrossRef]
  38. Bhushan, B.; Israelachvili, J.N.; Landman, U. Nanotribology: Friction, Wear and Lubrication at the Atomic Scale. Nature 1995, 374, 607–616. [Google Scholar] [CrossRef]
  39. Math, S.; Gao, J.; Landman, U. Interfacial Segregation, Structure, and Diffusion of n -Alkane Mixture Films Adsorbed on Smooth and Rough Gold Surfaces. J. Phys. Chem. C 2022, 126, 4209–4219. [Google Scholar] [CrossRef]
  40. Zhang, H.; Fukuda, M.; Washizu, H.; Kinjo, T.; Yoshida, H.; Fukuzawa, K.; Itoh, S. Shear Thinning Behavior of Nanometer-Thick Perfluoropolyether Films Confined between Corrugated Solid Surfaces: A Coarse-Grained Molecular Dynamics Study. Tribol. Int. 2016, 93, 163–171. [Google Scholar] [CrossRef]
  41. Wang, H.; Junghans, C.; Kremer, K. Comparative Atomistic and Coarse-Grained Study of Water: What Do We Lose by Coarse-Graining? Eur. Phys. J. E 2009, 28, 221–229. [Google Scholar] [CrossRef] [PubMed]
  42. Tang, J.; Kobayashi, T.; Zhang, H.; Fukuzawa, K.; Itoh, S. Enhancing Pressure Consistency and Transferability of Structure-Based Coarse-Graining. Phys. Chem. Chem. Phys. 2023, 25, 2256–2264. [Google Scholar] [CrossRef] [PubMed]
  43. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular Dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef]
  44. Junghans, C.; Praprotnik, M.; Kremer, K. Transport Properties Controlled by a Thermostat: An Extended Dissipative Particle Dynamics Thermostat. Soft Matter 2008, 4, 156–161. [Google Scholar] [CrossRef] [PubMed]
  45. Kobayashi, T.; Zhang, H.; Fukuzawa, K.; Itoh, S. Effect of Transverse Dissipative Particle Dynamics on Dynamic Properties of Nanometer-Thick Liquid Films on Solid Surfaces. Mol. Simulat. 2020, 46, 1281–1290. [Google Scholar] [CrossRef]
  46. Thompson, A.P.; Aktulga, H.M.; Berger, R.; Bolintineanu, D.S.; Brown, W.M.; Crozier, P.S.; In ’T Veld, P.J.; Kohlmeyer, A.; Moore, S.G.; Nguyen, T.D.; et al. LAMMPS—A Flexible Simulation Tool for Particle-Based Materials Modeling at the Atomic, Meso, and Continuum Scales. Comput. Phys. Commun. 2022, 271, 108171. [Google Scholar] [CrossRef]
  47. Stukowski, A. Visualization and Analysis of Atomistic Simulation Data with OVITO–the Open Visualization Tool. Model. Simul. Mater. Sci. Eng. 2010, 18, 015012. [Google Scholar] [CrossRef]
  48. Fukuda, M.; Zhang, H.; Ishiguro, T.; Fukuzawa, K.; Itoh, S. Structure-Based Coarse-Graining for Inhomogeneous Liquid Polymer Systems. J. Chem. Phys. 2013, 139, 054901. [Google Scholar] [CrossRef]
  49. Yamamoto, S.; Matsuda, H.; Kasahara, Y.; Iwahashi, M.; Takagi, T.; Baba, T.; Kanamori, T. Dynamic Molecular Behavior of Semi-Fluorinated Oleic, Elaidic and Stearic Acids in the Liquid State. J. Oleo Sci. 2012, 61, 649–657. [Google Scholar] [CrossRef]
  50. Wheeler, D.H.; Potente, D.; Wittcoff, H. Adsorption of Dimer, Trimer, Stearic, Oleic, Linoleic, Nonanoic and Azelaic Acids on Ferric Oxide. J. Am. Oil Chem. Soc. 1971, 48, 125–128. [Google Scholar] [CrossRef]
  51. Hess, B. Determining the Shear Viscosity of Model Liquids from Molecular Dynamics Simulations. J. Chem. Phys. 2002, 116, 209–217. [Google Scholar] [CrossRef]
  52. Greenwood, J.A.; Williamson, J.B.P. Contact of Nominally Flat Surfaces. Proc. R. Soc. A Math. Phys. Eng. Sci. 1966, 295, 300–319. [Google Scholar] [CrossRef]
  53. Tomanik, E.; Chacon, H.; Teixeira, G. A Simple Numerical Procedure to Calculate the Input Data of Greenwood-Williamson Model of Asperity Contact for Actual Engineering Surfaces. Tribol. Ser. 2003, 41, 205–215. [Google Scholar] [CrossRef]
  54. Bhushan, B. Surface Roughness Analysis and Measurement Techniques. In Modern Tribology Handbook, Two Volume Set; CRC press: Boca Raton, FL, USA, 2000; pp. 79–150. [Google Scholar]
  55. Fisher, C.H. N-fatty Acids: Comparison of Published Densities and Molar Volumes. J. Am. Oil Chem. Soc. 1995, 72, 681–685. [Google Scholar] [CrossRef]
Figure 1. Schematic of coarse-grained models (right) mapped from all-atom models (left). a, b, and c are the lattice constants in the x, y, and z directions, respectively.
Figure 1. Schematic of coarse-grained models (right) mapped from all-atom models (left). a, b, and c are the lattice constants in the x, y, and z directions, respectively.
Lubricants 12 00030 g001
Figure 2. Schematic of all potentials between CG beads.
Figure 2. Schematic of all potentials between CG beads.
Lubricants 12 00030 g002
Figure 3. Variation of viscosity with the damping coefficients between liquid beads. Markers and horizontal lines indicate the results from CG and AA simulations, respectively.
Figure 3. Variation of viscosity with the damping coefficients between liquid beads. Markers and horizontal lines indicate the results from CG and AA simulations, respectively.
Lubricants 12 00030 g003
Figure 4. Variation of CoF with the damping coefficients between liquid and solid beads. Markers and horizontal lines indicate the results from CG and AA simulations, respectively.
Figure 4. Variation of CoF with the damping coefficients between liquid and solid beads. Markers and horizontal lines indicate the results from CG and AA simulations, respectively.
Lubricants 12 00030 g004
Figure 5. Front view of symmetric models. The gray, light blue, and green represent the solid, nonpolar liquid, and polar liquid beads, respectively. The top arrow indicates the shear direction of the upper solid.
Figure 5. Front view of symmetric models. The gray, light blue, and green represent the solid, nonpolar liquid, and polar liquid beads, respectively. The top arrow indicates the shear direction of the upper solid.
Lubricants 12 00030 g005
Figure 6. Front view of asymmetric models. The gray, light blue, and green represent the solid, nonpolar liquid, and polar liquid beads, respectively. The top arrow indicates the shear direction of the upper solid.
Figure 6. Front view of asymmetric models. The gray, light blue, and green represent the solid, nonpolar liquid, and polar liquid beads, respectively. The top arrow indicates the shear direction of the upper solid.
Lubricants 12 00030 g006
Figure 7. (a) Adsorption trend throughout the simulation duration in the symmetric model. (b) Mean and standard deviation (indicated by error bars) of adsorption fraction during the last 30 ns of the simulations for surfaces with different wavelengths.
Figure 7. (a) Adsorption trend throughout the simulation duration in the symmetric model. (b) Mean and standard deviation (indicated by error bars) of adsorption fraction during the last 30 ns of the simulations for surfaces with different wavelengths.
Lubricants 12 00030 g007
Figure 8. (a) Fluctuations in OFM adsorption on rough surfaces in the symmetric model. (bd) Snapshots showing three characteristic states for adsorption of OFM molecules on the surfaces with 24 nm wavelength. Dodecane molecules are omitted for clarity. The gray, light blue, and green represent the solid, nonpolar liquid, and polar liquid beads, respectively. The red and blue indicate polar liquid beads that initially adsorbed and re-adsorbed onto the convex region of the lower surface, respectively. The arrows at the top of the snapshots indicate the shear direction of the upper solids.
Figure 8. (a) Fluctuations in OFM adsorption on rough surfaces in the symmetric model. (bd) Snapshots showing three characteristic states for adsorption of OFM molecules on the surfaces with 24 nm wavelength. Dodecane molecules are omitted for clarity. The gray, light blue, and green represent the solid, nonpolar liquid, and polar liquid beads, respectively. The red and blue indicate polar liquid beads that initially adsorbed and re-adsorbed onto the convex region of the lower surface, respectively. The arrows at the top of the snapshots indicate the shear direction of the upper solids.
Lubricants 12 00030 g008
Figure 9. (a) Adsorption trend throughout the simulation duration in the asymmetric model. (b) Mean and standard deviation (indicated by error bars) of adsorption fraction during the last 100 ns of the simulations for surfaces with different wavelengths.
Figure 9. (a) Adsorption trend throughout the simulation duration in the asymmetric model. (b) Mean and standard deviation (indicated by error bars) of adsorption fraction during the last 100 ns of the simulations for surfaces with different wavelengths.
Lubricants 12 00030 g009
Figure 10. Average adsorption fraction over the last 100 ns in different regions of rough surfaces with varying wavelengths.
Figure 10. Average adsorption fraction over the last 100 ns in different regions of rough surfaces with varying wavelengths.
Lubricants 12 00030 g010
Figure 11. Time for OFM desorption from convex regions on rough surfaces with varying wavelengths.
Figure 11. Time for OFM desorption from convex regions on rough surfaces with varying wavelengths.
Lubricants 12 00030 g011
Table 1. Determined values of damping coefficients and comparison of viscosity and CoF between coarse-grained and all-atom simulations.
Table 1. Determined values of damping coefficients and comparison of viscosity and CoF between coarse-grained and all-atom simulations.
ζ [(kcal·fs/(mol·Å2)]Coarse-GrainedAll-Atom
Viscosity, η [mPa∙s]
ζ l l n p = 4.1 1.341.18
ζ l l p = 36 5.986.24
CoF
ζ l s n p = 18 0.09060.0912
ζ l s p = 132 0.2110.196
Table 2. Surface properties based on Greenwood and Williamson’s rough surface contact model.
Table 2. Surface properties based on Greenwood and Williamson’s rough surface contact model.
Wavelength [nm]RMS Roughness, σ [nm]Mean Asperity Peak Radius, β [nm] [53]Mean Asperity Peak Density, ξ [nm−2] [53]Mean Asperity Slope (rad) [54] ξ β σ
(×10−2)
ξ β 1 / 2 σ 3 / 2
(×10−3)
61.2040.543127.77781.13041.817027.0585
122.14776.94440.56541.796413.4524
248.58341.73610.28271.79486.7232
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tang, J.; Chong, W.W.F.; Zhang, H. Coarse-Grained Molecular Dynamics Simulations of Organic Friction Modifier Adsorption on Rough Surfaces under Shear. Lubricants 2024, 12, 30. https://doi.org/10.3390/lubricants12020030

AMA Style

Tang J, Chong WWF, Zhang H. Coarse-Grained Molecular Dynamics Simulations of Organic Friction Modifier Adsorption on Rough Surfaces under Shear. Lubricants. 2024; 12(2):30. https://doi.org/10.3390/lubricants12020030

Chicago/Turabian Style

Tang, Jiahao, William Woei Fong Chong, and Hedong Zhang. 2024. "Coarse-Grained Molecular Dynamics Simulations of Organic Friction Modifier Adsorption on Rough Surfaces under Shear" Lubricants 12, no. 2: 30. https://doi.org/10.3390/lubricants12020030

APA Style

Tang, J., Chong, W. W. F., & Zhang, H. (2024). Coarse-Grained Molecular Dynamics Simulations of Organic Friction Modifier Adsorption on Rough Surfaces under Shear. Lubricants, 12(2), 30. https://doi.org/10.3390/lubricants12020030

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop