We find a peculiar velocity quantization phenomenon in the classical motion of an idealized solid lubricant, consisting of a harmonic chain interposed between two incommensurate sliders, as sketched at the right. The ratio v_{cm}/v_{ext} of the mean chain center-of-mass velocity to the externally imposed relative velocity of the sliders remains pinned to exact ``plateau'' values for wide ranges of parameters, such as sliders corrugation amplitudes, external velocity, chain stiffness and dissipation, and is strictly determined by one length ratio only.
The relevant lengths are:
a_{0} = lubricant rest periodicity
a_{+} = static substrate periodicity (taken as unit length)
a_{-} = periodicity of substrate sliding at v_{ext} (=0.1 in all reported simulations)
These lengths define two lengthscale ratios r_{+}=a_{+}/a_{0} and r_{-}=a_{-}/a_{0}. We assume, without loss of generality, r_{+} < r_{-} and 1/r_{+} < r_{-}.
The plateau dynamics of the model is illustrated in the animation at the right, where the (+) potential is drawn in red, and the (-) potential is drawn in blue. The solid lubricant chain trapped between two sliding periodic substrates is represented by circles. The two sliders act as shifting cosine potentials on every particle in the chain. The resulting model is a generalization of the classical Frenkel-Kontorova model. The animation has r_{+}=φ and r_{-}=r_{+}/(r_{+}-1)=φ², with φ = the golden mean (GM) = (1+5^{½})/2≈1.61.
v_{cm}/v_{ext} is constant for variations in particular of the spring constant K of the springs joining the lubricant particles (of unit mass). K=2 in the figure.
Phenomenological dissipation produces a force on particle i equaling - γ (2 v_{i} - v_{ext}), with γ=0.1.
Any of several frames of reference can be chosen according to convenience
to examine this dynamics.
The figure example uses the frame of reference where the
potential (+) is static
(v_{+}=0
v_{-}=v_{ext}).
Here
the plateau velocity is
v_{cm}/v_{ext} = 1 - 1/r_{+}.
The velocity plateaus are understood in terms of (soliton) density waves produced by density mismatch to the (+) potential and driven by the (-) potential. This point is clearest when r_{+} is near unity, as in the animation at the right, with r_{+}=σ, with σ = 1+(π/310)^{½}≈1.10. Once every 1/(r_{+}-1) minima, an extra lubricant particle must "squeeze in" to preserve the extra mean density given by the slightly smaller a_{0} < a_{+}. This extra density lump is precisely a soliton or kink: in this example a soliton occurs every about 10 lubricant particles.
One could visualize solitons more clearly by plotting again the particle position x_{i} on the horizontal scale, and moving the particles vertically proportionally to the distance x_{i} - x_{i-1} between neighboring particles (which is also the inverse of the local particle density). In this scheme, and in the (-) frame of reference, the previous animation becomes the animation at the right. One sees that neighboring particles move consistently closer in the correspondence of the minima of the upper (-) potential. These closeups are precisely the solitons (or kinks). Solitons pin to the (-) substrate, at its sliding velocity v_{ext}: in this pinned state the (-) frame of reference is the best one to visualize solitons.
In this animation, the dashed line represents the distance a_{+} at which the lubricant particles match the minima of the (+) substrate. The solid line represents the equilibrium interparticle separation a_{0}, i.e. the inverse mean density of the lubricant.
As the lubricant becomes stiffer (larger spring K), solitons gradually expand in size (involving more particles), interact more and more strongly, and decrease in amplitude. Please see the animation for K=100 to clarify this point. Eventually, when the chain becomes too stiff, the solitons density amplitude is so small that it ceases to pin to the (-) substrate: the plateau is then lost, and the particle speed quits the perfect plateau value v_{cm}/v_{ext} = 1 - 1/r_{+} (in the + reference). The animation of the non-plateau state K=300 shows an alternating formation and flattening out of the kinks lattice as it forces its way through the (-) corrugation.
On the opposite, soft-springs side, kinks melt away due to retardation: particles are just too slow to react and arrange their relative distance as required by the movement of the potentials, thus kinks get swept out, and the plateau ends. Of course, if v_{ext} is lowered from the 0.1 used in all calculations to, e.g. 0.01, then the plateau extends to softer chains. Additionally, the kinks formation mechanism (particles tend to stick to the minima of the (+) potential) is also altered if the amplitude of the (-) potential is larger than that of the (+) potential, and this can induce irregular motion and particles clustering, as illustrated in the animation of the non-plateau state K=0.5.
Like for kinks, the anti-kink quantized motion is best seen in the (-) frame of reference, again putting the inter-particle separation along the vertical direction. The animation at the right illustrates the anti-kink movement: "peaks" of increased interparticle separation (less dense lubricant) remain pinned to the maxima of the (-) potential.
In all previous examples we selected r_{-}=r_{+}/(r_{+}-1) as a convenient choice, since the distance between solitons matches the distance between minima of the (-) potential. In the present example instead the choice of r_{-} is generic, which means that the antikink densities and (-)-potential minima density are incommensurate, so that the anti-kinks arrange irregularly and do not form a perfect lattice. With this generic r_{-}, the density (1-r_{+})/a_{+}=0.11/a_{+} of antikinks exceeds the density 1/a_{-} = r_{+}/(r_{-}a_{+})=0.0995../a_{+} of (-) minima. Accordingly, occasionally a pair of antikinks pins to a single (-) maximum, as e.g. near the center of the figure at the right.
Simulations show that this soliton physics is so robust that it survives even to regimes r_{+} > 2 and r_{+} < 0.5 where one could not tell individual solitons apart. Also, finite size effect do not destroy this solitonic dynamics, as the finite (and very small!) size of the chain at the right demonstrates.
All nice things eventually come to an end, of course! The plateau state is lost as K becomes either too small or too large, as discussed above. The plateau is also destroyed by a very large sliding velocity v_{ext}, which does not leave the kinks time to form and pin to the (-) substrate. Also the application of an external force (acting on all particles) can lead to an unpinned state. If this force F_{ext} is large enough it makes the particles slide quickly, thus melting the kinks lattice, or at least making it slip against the (-) substrate. Because a force only produces an acceleration, which does not imply an immediate increase in velocity, an hysteretical behavior is observed as the force is slowly ramped up to the unpinned state and back down to the dynamical plateau state. The model shows an intermediate-force region of coexistence of pinned and unpinned behavior for the same F_{ext}, depending on the initial state. For example, see the animation of the pinned state obtained increasing F_{ext} to 0.012 and the animation of unpinned state obtained with the same F_{ext}=0.012, but obtained by decreasing F_{ext} from the large-F_{ext} unpinned state. Upon further decrease of the force to F_{ext}=0.0068828, the solitons remain unpinned, but intermittencies occur less frequently, and the speed approaches slowly the plateau one. These 3 calculations have K=100 r_{+}=σ, r_{-}=r_{+}/(r_{+}-1) as above. For more details about this hysteresis, see [cond-mat/0609117]
The oscillatory motions of the particles about their average drift can even be heard as audio files. The smoother, more harmonic oscillation for K=10 results in a smoother sound than the sharper bouncier K=1 motion, which is richer in harmonic contents. As the auditive sensation, and a more detailed Fourier analysis demonstrates, the oscillations of the particles around their average drift is periodic!
Nonperiodic (or rather quasi periodic) oscillations of the individual lubricant particles are found for other incommensuracy ratios. Hear for example the nonperiodic sounds of π/2.
We have found evidence of the quantized plateau dynamics also in a 2D model where the lubricant atoms can move also transversally to the sliding direction. The plateau velocity is supported by a lubricant as many as 6-layers thick. Preliminary calculations show evidence of the quantized lubricant motion also in a fully 3D model, based on Lennard-Jones interactions. Incommensuracy could even arise through rotation, as in this movie.
The impressive robustness of the velocity plateaus, due to their topological nature should be realized in nature, but need not be especially easy to observe, especially since it is a property of the lubricant dynamics, which is usually hard to probe. However, the plateau dynamics has tribological implications: frictional forces, those opposing the advancement of the (-) substrate decrease sharply as the quantized-sliding regime starts.
On the other hand, irregular asymmetric sliding of the lubricant can have potentially disruptive effects on lubricated friction in a macroscopical multi-grain context.
To date, we are aware of no clear experimental evidence of the velocity-plateau phenomenon.
A number of local works and publications report different sides of this kink-dragging mechanism, including:
http://materia.fisica.unimi.it/manini/theses/distante.pdf
http://materia.fisica.unimi.it/manini/theses/castelli.pdf
http://materia.fisica.unimi.it/manini/theses/caruso.pdf
http://materia.fisica.unimi.it/manini/theses/cesaratto.pdf
http://materia.fisica.unimi.it/manini/theses/divitini.pdf
created: 07 Mar 2006 | last modified: 27 Apr 2012 by Nicola Manini |