3.1 MSD and relative concentration
The relationship between the MSD of n-hexadecane molecules and time is shown in Fig. 3, which is used to characterize the diffusion of lubricant molecules during the sliding process of two friction pairs from the initial position to complete separation. Compared with the MSD value of the n-hexadecane molecule in the model of the substrate without ZnO layer, the MSD value in the model of the substrate with ZnO layer is smaller under the same simulation conditions. It means that the n-hexadecane molecules in the model of the substrate without ZnO layer show a more intense degree of diffusion than that of the model of the substrate with ZnO layer under the same conditions. Although the n-hexadecane molecules still exhibit diffusivity characteristic in the model of the substrate with ZnO layer, it is obvious that the ZnO layer limits the diffusion of n-hexadecane molecules. In addition, during the sliding process, the n-hexadecane molecule maintains a high degree of order and a low degree of dispersion under the internal action of ZnO layer.
The dispersion characteristics of lubricants molecules, which are in the two models of the substrates with and without ZnO layer, can be obtained through the relative concentration curve. Figure 4 and Fig. 5 exhibit the relative concentration changes of n-hexadecane lubricants molecules in the two models of substrates with and without ZnO layer during the sliding process. The relative concentration curves of n-hexadecane molecules in two Figs both have obvious peaks close to the surfaces of friction pairs and fluctuate slightly during sliding processes. These phenomenons indicate that the n-hexadecane molecules aggregate and form lubricating films on the upper and lower surfaces of friction pairs. 7 different values of relative concentration, which are the peaks in the curve, are obtained in the model of the substrate without ZnO layer in Fig. 4, while only 6 different values of relative concentration are obtained in the model of the substrate with ZnO layer in Fig. 5. These phenomenons mean that, during the sliding process, the lubricating fluids are divided into 7 and 6 lubricating films in two models, respectively.
As shown in Fig. 5, the thickness of lubricating films, which is about 25 Å, in the model of the substrate with ZnO layer is thinner than that of, which is about 31 Å, in the model of the substrate without ZnO layer. Furthermore, in the model of the substrate with ZnO layer, the relative concentration of lubricating films close to the surfaces of two friction pairs is about 4.3 and 4.7, respectively. The relative concentration of the other four lubricating films, which are in the middle of the curve, remains stable at around 3.5 and oscillates between 3.2 and 4.2 with distance. Comparing the curve of lubricating films in the model of the substrate without ZnO layer, in the model of the substrate with ZnO layer, the values of peak maintain higher, and meanwhile have smaller fluctuations during the sliding process. It illustrates that the distribution of n-hexadecane molecules in the model of the substrate with ZnO layer is relatively uniform and stable with considerably smaller scatter in the value of migration in the z-direction with time. In Fig. 4, the values of the peak are lower in the middle of the curve with fluctuating around three insignificant peaks. These results indicate that the formed lubricating films are not obvious and unstable in the model of the substrate without ZnO layer.
It is observed that the relative concentration close to the upper and lower surfaces in Fig. 5 fluctuates smaller with time. The phenomenons mean that the lubricating films, which are close to the upper and lower surfaces, spread and slid at the same speed with the friction pair sliding. These two lubricating films are called adsorption films, which adsorb on the surface of friction pair. Based on the curves in Fig. 4 and Fig. 5, the lubricating films maintains stability during the sliding. Actually, during the sliding process, the interface molecular film has three different physical statuses, which are solid-like, amorphous, and liquid-like[40]. Based on classical fluid mechanics theory the strong adsorption force of alkane molecules on the upper and lower surfaces can make lubricating film showing solidity, forming a thick solid adsorption film to participate in shear lubrication, and meanwhile, the adsorption force limits the spread of alkanes[41]. Therefore, to analyze the interaction effect between the lubricants molecules and ZnO layer, the interaction distance between them needs to be investigated.
3.2 RDF and coordination number
Figure 6 draws the RDF curves of the interaction distance between the lubricating molecules and the surface molecules of the friction pair in the models of the substrate with and without ZnO layer, respectively. The abscissa r represents the interaction distance between the surface molecules and the lubricants molecules. In the two curves, the peaks and valleys are observed significantly. Initially, the RDF values are all zero with r < 0.75Å. Afterward the curves steeply sharp to the maximum values at the first peak, r = 2.25 Å. It implies that a large number of n-hexadecane molecules gather near the surfaces of friction pairs at r = 2.25 Å, meanwhile, ZnO molecules on the surface of the friction pair polarize the n-hexadecane molecules close to the surface. The density of lubricants molecules near the peaks is much higher than the average density of the curves. The deep valleys, which represent the weak interaction areas, between the two peaks are observed in the curves. The second peaks of both curves represent the medium degree of connection strength. For the model of the substrate without ZnO layer, the values of the RDF curve are generally smaller than the values of the RDF curve obtained in the model of the substrate with ZnO layer.
Molecular aggregations are caused by the distance of molecules reducing, and moreover, a smaller distance from the reference atom means a stronger aggregating degree of molecules. Compared with copper and zinc atoms, ZnO molecules have stronger interaction with lubricants molecules, enabling the lubricants molecules to exist in a more compact form on the surfaces during the sliding process. The curves abscissa corresponds to the interaction distance between the surface molecules and the lubricants molecules. In Fig. 6, the first peaks of the curves indicate that the distance between the lubricating films near the surface and the surface is all 2.25 Å. For the model of the substrate without ZnO layer, the second peak in the curve is smooth, furthermore, to the end of the curve, the curve shows the weak interaction between the Cu-Zn alloy molecules and the lubricants molecules.
The coordination number is used to measure the number of n-hexadecane molecules surrounding the surface of the friction pair. The coordination number not only reflects the binding ability and coordination relationship between lubricants molecules and the surface molecules, but also describes the tightness of the lubricants molecules during the sliding process. The larger the coordination number means the tighter molecules arrangement. According to formula (3), the position of the first peak of the RDF in Fig. 6 is Rmin=2.25 Å, furthermore, the number of carbon atoms in the adsorption films, which are on the two kinds of friction pair surfaces, was calculated to characterize the tightness of n-hexadecane molecules. For the substrates with and without ZnO layer, the coordination numbers are 119.4 and 91.4 respectively. It can be obtained that the carbon atoms of n-hexadecane gathered in the lubricant films, which are the closest to the surfaces of substrate with ZnO layer, is about 1/4 of the total number of the simulated carbon atoms in the adsorption film. Furthermore, this coordination number is about 1.3 times that of the model of the substrate without ZnO layer. According to the RDF image and coordination number results, the ZnO layer leads to more n-hexadecane molecules gathering near the friction surface to form adsorption film, which proves the result of Fig. 5. Due to the hydrogen bond and van der Waals interaction with the ZnO molecules, adsorption film adsorbs on the surface of ZnO layer, moreover, the adsorption capacity decreases with the increasing distance from the surfaces. Figure 7 shows the snapshot of the sliding process in the model of the substrate with ZnO layer and the sketch of films movement during the sliding process, respectively. Figure 7(a) exhibits that, in the model of the substrate without ZnO layer, the lubricants protect the friction surface and are subjected to shear stress during the sliding process. Figure 7(b) is the schematic diagram of lubricating films movement during the sliding process. The lubricating films are arranged in parallel in the film thickness direction and maintain the same movement direction. Therefore, the transmitting motion behavior of lubricating films, which include adsorption film, needs to be investigated in detail.
3.3 Shear stress
Based on the above results, the interaction between the n-hexadecane molecules and the surface of the friction pair belongs to physical adsorption behavior. During the sliding process, the lubricating films separate the two surfaces of the friction pair, meanwhile, the friction force is caused by conquering the shear stress between the two surfaces. Adsorption film adsorbs on the surface of friction pair to participate in the sliding. Figure 8 shows the shear stress on the upper and lower surfaces of the friction pair of the model without ZnO layer during the sliding process. It is observed that the stresses on the upper and lower surfaces maintain and fluctuate near their own fixed values of the vertical axis, respectively. Similar phenomenons are shown in Fig. 9 (a), which is the shear stress change on the upper and lower surfaces of the friction pair of the model with ZnO layer during the sliding process. However, the overall shapes of stress curves in Fig. 9 (a) show more regular than these in Fig. 8. By observing the enlarged curves in Fig. 9 (b), during the sliding process, the shear stresses fluctuate regularly within a certain range and exhibit periodic characteristics. The stick-slip effect, which is caused by the interaction of molecules, between the friction pairs and the lubricating films leads to a regular curve shown in Fig. 9 (a), and moreover, the hilly bulge area on the shear stress in Fig. 9 (a) is the significant feature of the stick-slip effect. During the process of the lubricating films from adhesive contact to the beginning of sliding, the shear stress gradually increases from the minimum dynamic shear stress to the maximum static shear stress, then the lubricating film molecules begin to slide and enter a new adhesive contact with the friction pair surface, and meanwhile, the shear stress gradually decreases to the minimum dynamic shear stress. During the sliding process, the alternate adhesion and sliding of the lubricating films reduce the wear of the friction pair. Comparing the shear stress curves in Fig. 8 and Fig. 9(a), it is inferred that the ZnO layer on the friction surface causes periodic shear stress and maintains a stable stick-slip effect, which results in a firm adsorption molecular film structure and a strong interaction between the surface and lubricants.
Based on the result of Fig. 6, the strong polarity, which is between the ZnO lubricants on the surface of the friction pair and the n-hexadecane molecules, leads to the strong interaction between the surface and lubricants. With the combined effect of Van der Waals force and electrostatic force, lubricants molecules form lubricating films during the sliding, which are shown in Fig. 10(a). Between these lubricating films, low concentration cavities are generated due to the lubricants molecules gathering up or down to the lubricating films. The polar molecules in the lubricating film are arranged in regular order on the surface of the film. During the sliding process, the polar molecules on two adjacent lubricating films contact and slide with each other to produce the stick-slip effect. With the action of the stick-slip effect, the polar molecules in one lubricating film transfer the motion to the next lubricating film, resulting in periodic shear stress as shown in Fig. 9 (b).
Figure 10 (b) shows the arrangement state of n-hexadecane molecules in the model of the substrate with ZnO layer. The results of the coordination number indicate a large number of lubricants molecules in adsorption film are aggregated on the surface of the substrate with ZnO layer. Therefore, compared with lubricants molecules adsorbed the surface of the substrate without ZnO layer (Fig. 10 (a)), more lubricants molecules are aggregated and bear the load on the surface of the substrate with ZnO layer. Figure 11 describes the change of pressure on the surface of the substrate. The pressure on the surface of the substrate with ZnO layer drops at the sliding start, and then stabilizes at about 2.4Gpa finally, and meanwhile, the pressure values on the surface of the substrate with ZnO layer are higher than that on the surface of the substrate without ZnO layer. Under the same sliding conditions, the ZnO layer causes lubricants molecules to be arranged regularly on the surface into a stable structure, which sustains greater loads.