1. Introduction
Lubrication is crucial for the functional performance of a wide range of mechanical and biological interfaces. Most interfaces operate under a mixed regime of lubrication, meaning that a certain degree of direct interaction of stochastic solid contacting surfaces occurs. Therefore, the conditions endured by many interfaces require the formation of an intervening fluid film to extend beyond its mere bulk microscale rheology. Protective boundary-active additives are used to form surface-adhered/bonded tribofilms. These can contribute to a reduction in contact friction and prevent wear in the nanoscale interactions of surface asperities in mechanical as well as biological interfaces. Hitherto, despite significant improvements in lubricant technology, it is estimated that 23% of the world’s total energy consumption is dissipated in tribological contacts [
1].
Typically, additives form around 5–24 wt% of fully formulated lubricants for engineering applications [
2,
3]. Depending on application, lubricants comprise a blend of hydrocarbon base stock, as well as detergents, dispersants, antioxidants, friction and viscosity modifiers, anti-foaming agents, pour point depressants, corrosion inhibitors and anti-wear additives [
3]. In the formulation, the constituent physical chemistry enhances certain aspects of the lubricant system performance and can be adjusted to suit specific applications.
Despite their crucial intended functions, the mechanisms through which most additives adsorb or bond to surfaces remain somewhat empirical, particularly due to their potential synergistic as well as antagonistic competition in occupying attachment sites on boundary surfaces [
4]. The lack of a fundamental understanding is exacerbated by the significant difficulty associated with taking in situ measurements of additives’ activation and tribofilm characterisation.
Molybdenum-based additives were first patented in 1939 in the form of colloidal molybdenum disulphide (MoS
2) as an anti-wear and friction modifying additive [
5]. The use of oil soluble molybdenum dialkyldithiophosphate (MoDDP) was reported in 1955 as an alternative to zinc dialkyldithiophosphate (ZDDP) which was already in use [
6]. Other forms of oil soluble molybdenum-based friction modifiers have also been used commercially, including molybdenum dithiocarbamate (MoDTC). The reaction process that transforms an oil soluble molybdenum-based molecule into a functional friction-reducing layer of MoS
2 is chemical and requires energy and suitable reactants for initiation.
The Arrhenius equation has long been applied in chemistry to model the rate-limiting step in a reaction process [
7]. Through modifying the Arrhenius equation, found in the seminal works of Prandtl [
8] and Eyring [
9], the mechanical contribution to reaction rates can also be considered as an energy input. Tysoe [
10] and Spikes [
11] have recently provided detailed summaries of the evolution of these concepts and their application in tribology through tribochemistry or mechanochemisty.
The shear-induced thermal activation model has recently gained prominence following the work of Gosvami et al. [
12] on in situ activation of ZDDP via scanning Atomic Force Microscopy (AFM). Using a range of contact pressures and temperatures, Gosvami et al. [
12] recorded the change in surface height at intervals after scanning with an AFM tip over a surface submerged in a lubricant containing ZDDP. The tip was scanned parallel to the cantilever support in order to promote additive activation, and then operated in tapping mode over the scanned area and its surroundings to assess the formed tribofilm height. By evaluating the increased film height after intervals of time, they were able to determine an exponential film growth rate which they linked to an Arrhenius model using the generated contact pressure as the source of supplied mechanical energy.
Zhang and Spikes [
13] investigated the formation of ZDDP tribofilm using a Mini Traction Machine (MTM) to provide the necessary input mechanical energy and interferometry to measure the tribofilm growth. They adjusted the slide-roll ratio and applied contact load, thus the resulting shear and contact pressure, respectively. They concluded that shear, not normal pressure should be used as the main source of supplied mechanical energy, also when using the Arrhenius equation as a predictive tool. Zhang and Spikes [
13] used the Eyring’s reaction path concept [
14] to demonstrate that induced shear not only lowers the activation energy barrier for the chemical reaction of boundary active additives, but also increases the potential barrier for the reversal of the reaction. Adams et al. [
15] showed that for the case of dimethyl disulfide on a copper substrate in Ultra-High Vacuum (UHV), a Surface Force Apparatus (SFA) can be used to activate the molecules through shear and bond them to the surface. They made full use of the UHV environment by employing angle-resolved X-ray Photoelectron Spectroscopy (XPS) to track the chemical composition of the surface following the SFA activated chemical reaction.
The Prandtl [
8] and Eyring [
9,
14] models apply a force-distance term (work done:
) as the mechanical input energy, whilst Gosvami et al. [
12] used the pressure volume term (
) and Zhang and Spikes [
13] used a shear stress-volume term (
). The pressure and shear stress terms effectively represent the characteristic forces applied in the contact, whilst the volume terms (referred to as the activation volume) are necessary to provide a dimensionless exponential term in the Arrhenius equation. Zhang and Spikes [
13] and Tysoe [
10] provided experimental evidence that shear stress is the most pertinent characteristic force expression for the mechanical contribution to activation. The activation volume can be considered for vertically orientated molecules to be the area of the surface occupied by the molecule multiplied by the activation length [
11]. Previous studies have suggested that a molecular area between 20 and 400
and an activation volume between 7.8 and 180
are applicable for various molecules, including ZDDP, fatty acids and base polymers [
13,
16,
17,
18,
19,
20].
A number of studies have measured the rate of growth of tribofilms with both fully formulated and partially formulated lubricants [
21,
22,
23,
24,
25,
26,
27]. These studies have shown that there is a slow initial growth period, followed by an increased growth rate and finally a steady state region with a minimal increase or decrease in tribofilm thickness with time. This is particularly noteworthy in the study of Gosvami et al. [
12] who fitted a linear relationship to the initial growth stage and modelled the subsequent faster growth rate with the Arrhenius model. They noted a “limiting value”, beyond which the growth appears to cease and some small patches of the tribofilm are removed from the edges of the scan area under extreme AFM contact loads (600–1700 nN). An important point to note is that neither the slow initial growth rate nor the limiting value can be explained by the Arrhenius model alone.
Akchurin and Bosman [
28] used the data from the experimental work of Gosvami et al. [
12] and applied a simple wear term to explain the observed plateauing of the growth rate. This wear rate term was included with the Arrhenius growth model and showed good agreement with the experimental results, although the underlying physical nature of the wear mechanism was not discussed in any detail. Many papers [
29,
30,
31] have shown that wear at the nanoscale does not follow the behaviour demonstrated at larger scales described by classical wear models [
32]. Therefore, despite the agreement between the models presented for AFM (nanoscale) [
28] and the tribometric activated film growth (microscale) [
24], it is unlikely that the same wear mechanism would be dominant at both physical scales. This is particularly evident as the film height is not seen to reduce in AFM studies, where the rate of growth merely decreases, and the film height stabilises. This has been hypothesised to be as the result of combined stable growth and wear rates acting in concert [
24].
AFM is well-suited to the investigation of activation of lubricant additives. However, its use in the simultaneous measurement of surface height and friction in Lateral Force Microscopy (LFM) has not hitherto been reported in the literature. Additionally, there are some gaps in the current understanding of film growth characteristics, particularly in relation to the initial and limiting film growth mechanisms. This paper presents an advance in the understanding of tribofilm growth through novel simultaneous in situ nanoscopic measurement of friction and film growth height. In addition, previously presented models are adopted and further developed in order to provide better predictions of the early stage tribofilm growth through the inclusion of an adsorption model. It is shown that the inclusion of an adsorption model provides for better prediction of early stage tribofilm growth. It explains the observed limiting maximum tribofilm thickness and allows for the inclusion of additive presence at the adsorption sites.
2. Experimental Methodology
The current investigation uses a Group III high viscosity base oil with a 0.5 wt% concentration (approximately equating to 500 ppm) of MoDTC. The concentration MoDTC used here guards against interruption in tribofilm growth on account of a limited supply of the additive. MoDTC is a typical engine oil additive for generation of low friction MoS
2 films and provides anti-oxidation properties [
33]. The MoS
2 films are generally produced in tandem with ZDDP films for combined wear protection and low frictional performance where direct boundary interactions occur [
34]. The base oil used here has been investigated in detail for its shear characteristics as well as for boundary interactions of this additive at microscale tribometry, using a s range of engineering surface materials [
35]. The investigation considered the interaction of a reciprocating sliding tribometer, intended to replicate a piston ring to cylinder liner conjunction, with different lubricant-material combinations. These included Nickle-Silicon-carbide, Borided Steel and EN14 steel (for the sliding strip). These materials were tested for frictional losses with the use of base oil containing only a single friction modifier species. The base oil, without friction modifier, was also investigated in isolation using LFM to determine the Eyring shear stress of the lubricant. The substrates used in the current investigation are super-finished EN31 steel specimens with isotropic topography and roughness of approximately 0.05 µm (Sa). These plates were washed in an ultrasonic bath with petroleum ether (40–60) for a period of 10 min and air dried prior to experimentation.
To determine a base line coefficient of friction for the test specimens, friction was measured in the absence of any lubricant for a range of AFM probe tip contact loads. Bruker DNP-10 twin-arm silicon nitride cantilever with a nominal tip radius of 20 nm was used for these measurements. In order to characterise the cantilever a “blind” calibration technique was employed in line with that proposed by Buenviaje et al. [
36] using a TGF 11 monocrystalline silicon grating. Following the cantilever calibration, the tip geometry was determined by taking a topographical measurement of a TGT 1 silicon wafer surface with calibrated peaks. Off-line deconvolution of the topographical data enables the determination of tip geometry. The measured tip radii were within 1 nm of those specified by the manufacturer, when averaged over 20 peaks from the calibration sample’s scanned area.
Measurements were taken from 3 areas of all the EN31 steel samples, with a new tip calibrated and used for each test. The steel sample was a 50 mm × 20 mm × 5 mm plate with two counter-sunk mounting holes. A range of loads were applied to the sample and a 1 × 2 µm
2 area scanned (2 µm in the fast scan direction) at a sliding speed of 2 µm/s, whilst simultaneously recording the applied load, friction and the height of the cantilever. The results are shown in
Figure 1, indicating a coefficient of friction of 0.27 for the nominally dry tested surface.
Figure 1 shows very good agreement between the results of the three tests, indicating that the procedure and sample surface topography are both repeatable and precisely controlled. The coefficient of friction is found as the gradient of the line and is in good agreement with similar tests on EN14 steel [
35,
37]. The contact is taken as nominally dry since the tests were conducted in a controlled atmosphere with 50 ± 2% relative humidity resulting in a slight offset to the intercept, in agreement with water-surface condensation and meniscus growth models [
38,
39]. This is assumed to produce a constant offset to the friction trace, providing a small positive intercept to the graph in
Figure 1, but insufficient to cause any hydrodynamic lift at the slow sliding speeds investigated.
To investigate the in-situ growth of the tribofilm, a tip load of 160 nN was applied which produces a mean contact pressure of approximately 5 GPa. The scan rate was set at 29.6 Hz in order to allow 5000 full scans of a 1 × 2 μm
2 (2 µm in the fast scan direction) region over a period of 6 h. These values are similar to those employed by Gosvami et al. [
12] and have been shown to be sufficient for the formation of ZDDP-based tribofilms. These conditions result in a very high tip sliding speed of 118.4 μm·s
−1 which is found to result in non-Newtonian shear behaviour of the lubricant used in previous investigations [
35].
The experimental procedure for additive activation is commenced by following the same calibration procedure described for the dry case, followed by measuring the surface in a nominally dry state. The initial measurements were made to ensure that measured friction was in line with the aforementioned initial dry testing condition and provide a baseline surface height. The lubricant was then injected into the contact while the tip was loaded against the surface, thereby restricting any drift. Sufficient lubricant was added to ensure that the generated meniscus was far removed from the scanning location. Continuous friction and height measurements were conducted for 5000 surface scans over a 6 h test duration. This procedure was repeated three times under the same scan conditions with a different tip and test area in each case. The bulk material and ambient temperatures during the testing was controlled to 20 ± 0.5 °C.
3. Experimental Results
The results of three separate tests are provided in
Figure 2 and
Figure 3, showing the cumulative tribofilm height relative to the initial scan and the evolution of coefficient of friction with time.
An overlay of the mean measured coefficient of friction and the mean surface height for the three scans can be seen in
Figure 4.
Figure 2,
Figure 3 and
Figure 4 show several distinct trends. Firstly, there is a clear initial drop in friction relative to the nominally dry contact measurement which is attributed to the reduction in adhesion through wetting of the contact under partial (mixed) regime of lubrication in line with that proposed by Bowden and Tabor [
40], Tambe and Bhushan [
41] and Chong and Rahnejat [
42]. Secondly, following this initial drop relative to dry friction, there is a significant, but unstable further reduction in friction over the first 45 min of scanning. This coincides with an unstable film growth rate. Referring to
Figure 2 and
Figure 3, for this region of variations an unstable film growth can be observed for tests 1 and 3, but not for test 2 where a steady reduction in friction with a steady film growth is observed. In the case of tests 1 and 3, an increase in friction is observed prior to a stable coefficient of friction being reached after 120 min of scanning. Test 2 displays a slower frictional response rate, also reaching a stable coefficient of friction after approximately 120 min.
Contrary to the evolving frictional characteristics, there is almost negligible change in the tribofilm height for the initial 30 min scanning duration, followed by an unsteady film growth for the next 30 min interval for tests 1 and 3, whilst in test 2 a slower, but smoother change in the growth of film height is observed. Finally, for all tests, the film growth stabilises after an hour and appears to plateau near the completion of test duration. The film growth in
Figure 2, between approximately 120 and 300 min, appears to follow a near linear growth rate which is in agreement with the shear-promoted thermal activation model as the contact conditions and coefficient of friction remain almost unaltered. However, it is essential to modify the shear-induced thermal activation model in order to account for the reduced growth rate nearing the end of these tests and take into account the reduced film formation rate observed at the beginning of the tests, despite relatively high friction.
Several studies have considered the wear of AFM tips through either scanning electron microscopy [
43,
44], or measurement of calibration samples [
45] or even by measuring a nominally worn tip using a second fresh tip [
29]. These techniques are not possible in the current study due to the application of the lubricant to the sample surface and the fragility of the components. The tip wear reported elsewhere has been rather slight and it is unclear what effects the lubricant and the formed tribofilm would have on the tip wear rate. Therefore, the tip wear is neglected in the current study.
4. Discussion
The stress-dependent film growth rate is predicted by the Arrhenius shear-promoted thermal activation model as:
where
is the growth rate with the unit:
,
is referred to as the pre-exponential term,
is the internal activation energy, given as the energy barrier in the absence of any applied stress,
is the Boltzmann’s constant (
),
is the absolute temperature and
is the expended mechanical energy as:
where
is the mean value of stress component affecting the activation barrier, in this case the contact shear stress and
is the activation volume.
where
and
are described by Gosvami et al. [
12] as the attempt frequency (i.e., the rate of interactions between an additive molecule and the target surface) and the molar volume (the volume occupied by one mole of the growth species) respectively.
In the current application, and should be defined rather differently as the mechanical input energy is produced through shear rather than application of pressure. Therefore, the attempt frequency, , becomes the number of molecules sheared by the AFM tip per second. This can be estimated by the area occupied by a single molecule on the surface, the footprint tip contact area and its sliding speed. The volume, , is therefore no longer the molar volume, but the molecular volume (i.e., the volume of the product produced by the reaction of one molecule). is, therefore, regarded as the product of the area occupied by a molecule and the molecular-chain length.
Applying these definitions and the conditions of the experiment to the Arrhenius shear-promoted thermal activation model and adjusting for the best fit to the experimental data, yields the results shown in
Figure 5. The best fit was found by applying the following conditions:
,
,
and
.
However, as
Figure 5 shows, the Arrhenius shear-promoted thermal activation model clearly does not predict a limiting value or any initial state of transience. Therefore, it is rather inadequate for the evaluation of full tribofilm forming characteristics.
MoDTC and ZDDP are polar molecules and can adsorb to metallic surfaces through van der Waals interactions [
3]. This provides for a high viscosity weakly adsorbed near surface layer [
34], which can be considered as the mechanism for the initial stages of tribofilm formation, whereas there would not yet be a significant bond or performance benefit to the additive film formation in the case of MoDTC and ZDDP. This is thought to be the form an additive film takes through application of thermal energy alone. Gosvami et al. [
12] showed that for ZDDP such a layer is formed to a depth of approximately 10 nm but was easily removed using the scanning conditions that were later shown to induce film growth.
From the state of an adsorbed layer, sufficient energy input is required to form a more permanent structure and a stronger bond to form with the target surface (a process commonly referred to as ‘activation’). Adams et al. showed that, in UHV tribometry with a SFA, molecules bonded in a similar structure which can be sheared into activation [
15]. This proposition is supported by the work of Zhang and Spikes [
13], among others, at the microscale of tribometry. It is hypothesised that with adequate shear the polar attachment of the additives forms an anchor-point while the free ends are forced over and laid in the direction of shear [
46]. Essentially, it is thought that the shearing action in combination with the available thermal energy would force the free ends of molecules into proximity with each other or with the solid surface, forming covalent bonds and thus promote the formation of a tribofilm.
The requirement for a molecule to be adsorbed to the surface prior to activation by the sliding tip provides a physical restriction to the chemical reaction rate predicted by the Arrhenius model. The finite coverage of molecules adsorbed to the surface should also be considered unless a large period of time is elapsed between successive scanning passes for a film covering the surface to stabilise. Therefore, the rate of adsorption must be taken into account in the modelling of the activation process.
Langmuir is credited with the first mathematical model for surface adsorption [
47,
48,
49]. The Langmuir model considers the formation of an adsorbed monolayer onto a surface from a gas as a function of the adsorption, desorption rates and surface coverage through a simple expression as:
where
is the rate of adsorption,
is the rate of desorption,
is the fraction of the area covered by a deposited layer and
is time.
A Langmuir adsorption model has been applied previously to the case of lubricant additive film formation by Dacre and Bovington [
50,
51]. They investigated the effect of concentration, temperature and stirring on the adsorption of ZDDP on steel and determined values of
and
for a range of cases.
Although the Langmuir model allows for a limiting growth rate at the onset of the scanning process, it does not provide an explanation for the limiting growth observed at the culmination of the experiments. A possible reason for the reduction in the growth rate and for the initial transience is provided through inclusion of a multi-molecular layer adsorption model. The first multi-molecular layer adsorption model was provided by Brunauer, Emmett and Teller [
52] and have since been named after them as the so-called BET model. Under the proposed conditions, additive molecules would have the ability to adsorb in layers upon each other in the form of a stack. Desorption can also occur at all points within the stack such that a steady state condition is always reached with a finite number of layers. Under these conditions the shear experienced by each molecule at each adsorption site may be considered as inversely proportional to the number of molecular layers resident there [
53,
54].
A BET multi-molecular layer adsorption model can be applied and solved analytically for a finite maximum number of layers. For the current case, an initial condition at which there would be no adsorbent coverage on the surface was applied. In the time between successive passes of the AFM tip, the growth on any layer
can be described as:
where
denotes the descrete time step,
accounts for the probability that any molecule in the stack can desorb to reduce the coverage at layer
and the growth rate would be limited by the fraction coverage of the layer
, multiplying through by
[
55].
Applying the multi-molecular layer adsorption model to the Arrhenius equation yields:
where
is the maximum number of molecules in a stack accounted for in the simulation study.
Equation (6) shows that the effective shear stress term in the mechanical energy is reduced as the number of adsorbed layers increases. The surface monolayer adsorption and the shear-promoted thermal activation can therefore be modelled as a symbiotic loop in which molecules adsorb over a given period of time, with a scanning pass activating some proportion of these. Therefore, a new surface coverage is determined, and adsorption is allowed over the next time period before another scanning pass occurs.
Another factor that should be considered concerns the removal of the adhered molecules which are unsuccessful in the bonding reaction. If the expended shear energy is sufficient, these molecules may be removed from the surface and return to their dissolved state in the bulk oil. As the desorption rate is energy dependent, it would be logical to relate the desorption rate to the generated thermal and mechanical energies, while only using the former for desorption occurring outside the contact domain. As a result, a process may be conceived as illustrated in
Figure 6. A similar process was described by Heuberger et al. [
56] for the tribochemical film formation of ZDDP using MTM.
Applying the multi-molecular layer model yields the results shown in
Figure 7, where the average depth of the adsorbed film is added to the activated tribofilm height predicted by the Arrhenius equation. An estimated layer height of 2.4 nm has been applied here as this has been reported as the mean molecular length of free ZDDP molecules determined by Georges et al. [
57] using dynamic light scattering.
Figure 8 shows the predicted fraction of the surface covered with bonded molecules 1, 2, 3, 4 and 5 layers deep within the test period. It can be seen that the mono-layer coverage reduces after 60 min as it becomes increasingly covered with further molecular layers.
Figure 9 shows that the total predicted surface coverage stabilises at 85% of the surface after 180 min, although there remains slight fluctuation due to changes in shear stress and activation rate.
Referring back to
Figure 7, it is shown that the multi-molecular layer adsorption model not only provides a reduced initial growth rate in comparison with the use of the Arrhenius model alone (
Figure 5), but also provides an explanation for the reducing film growth at the end of the test period. With additional testing for further understanding of the adsorption and desorption rates for each of the multi-molecular layers, it is expected that even better conformance of predictions to the measured data may be achieved.
In micro-scale contacts, some change would have to be made to the cyclic process shown in
Figure 6 as the adsorbed multi-molecular layers, predicted through BET, may be sheared-off under the contact conditions, particularly in asperity-pair ploughing. The unusual case of the AFM scanning probe removes this issue in the current investigation as the fixed end of the cantilever holding the tip is a reactive loading system, which adjusts the tip height in order to maintain a constant load over varying topography. This is of course, not the case in micro-scale contacts.
Figure 7 and
Figure 8 show that the rate of growth of the tribofilm predicted through the new combined BET and Arrhenius model begins to fall as the adsorbed film tends to an average depth of 3 molecules. This can be seen to occur even though there is ample adsorbed film coverage as shown in
Figure 9.
The inclusion of an adsorption model in the prediction of tribofilm growth offers an opportunity for the development of competitive/antagonistic additive technology, where boundary active additives compete for adsorption sites. It has been observed that certain additive combinations provide harmonic/synergistic or disruptive/antagonistic relationships. In the case of MoDTC and ZDDP, it has been noted that there is a harmonic/synergistic relationship in automotive lubricant applications where low shear strength MoS
2 layers form within the ZDDP-based tribofilm, offering a slip plane which reduces friction, whilst ZDDP offers exceptional wear protection [
34,
58]. However, if the ratio of adsorption rates between the ZDDP and MoDTC are altered, then the ratio of the tribofilm composition can be tailored accordingly. MoDTC is one of several lubricant-soluble polar molecules capable of producing a MoS
2 tribofilm. Therefore, it may be possible to make a selection based on the predicted performance by the expounded model. Alternatively, detergents such as calcium sulfonate have been shown to disrupt the growth of ZDDP, leading to essential compromises in lubricant formulation [
59,
60]. These relationships are well-suited for predictions with a competitive adsorption model as they allow a limited number of surface sites to be accounted for prior to the additive activation event. As a result, it is possible to determine a fractional content of subsequent layers in combined MoDTC and ZDDP competition and the growth of ZDDP in competition with a detergent molecule as well as other surfactant combinations.
Ostwald and de Izaguirre [
61] developed an early competitive adsorption model, based on the work of Freundlich [
62] and a similar model also exists for Langmuir-based adsorption [
63].
Similarly, the current approach allows for discrepancies between additives’ performances observed for different surface materials or coatings. In a typical fully formulated lubricant, there are many species competing for bonding/adsorbing to a finite number of surface sites, and variation in the surface forces is likely to cause significant differences in tribofilm growth and subsequent lubricant system performance in mixed and boundary regimes of lubrication. The analytical model proposed shows reasonable correlation to the test results in
Figure 7. However, the average results were generated from a limited number of tests supported by observed trends [
12].