Theory and simulation of sliding friction

Friction and nonlinear dynamics

The nonlinear dynamics associated with sliding friction forms a broad interdisciplinary research field that involves complex dynamical processes and patterns covering a broad range of time and length scales. Progress in experimental techniques and computational resources has stimulated the development of more refined and accurate mathematical and numerical models, capable of capturing many of the essentially nonlinear phenomena involved in friction.

Onset of sliding in amorphous films triggered by high-frequency oscillatory shear

We investigate the change of the static friction threshold of weakly adhesive amorphous interfaces in the presence of the shear ultrasonic oscillation. Prior to sliding, a softening of the shear interfacial stiffness is observed under either static or high-amplitude oscillatory shear. We find that the nonlinear shear ultrasound, regardless of its polarization, triggers the macroscopic sliding at these interfaces far below the static threshold. Such unjamming transition is due to the vibration-induced decrease of the apparent coefficient of static friction, which provides a mechanism for understanding the reduction of the yielding threshold of granular media by the acoustic fluidization.

Measurement of the friction between single polystyrene nanospheres and silicon surface using atomic force microscopy

In the present work, the individual nanoparticles have been manipulated on a silicon surface, using atomic force microscope (AFM) techniques. As a model system, near-spherical polystyrene nanoparticles with radii from 28.85 nm to 228.2 nm were deposited on a nanosmooth silicon wafer. Experiments demonstrated that when the normal force is above a threshold load, nanoparticles could steadily be pushed by the tip of the AFM along the defined pathway. The tests allow us to quantitatively study the interfacial friction between the nanoparticle and the surface. It was found that the friction could be affected by various factors such as the load, the particle size, and the surface treatment. The results showed that the friction between particles and substrate is proportional to the two-third power of the radius, which is in agreement with the Hertzian theory. It can also be seen that the ratio between the kinetic and the static friction was slightly changed from 0.3 to 0.6, depending on the size of the particles. However, the value of the ratio was little affected by other factors such as the particles' location, the tip normal force and the surface modification. The results provided new insights into the intriguing friction phenomenon on the nanoscale.

Onset of frictional slip by domain nucleation in adsorbed monolayers

It has been known for centuries that a body in contact with a substrate will start to slide when the lateral force exceeds the static friction force. Yet the microscopic mechanisms ruling the crossover from static to dynamic friction are still the object of active research. Here, we analyze the onset of slip of a xenon (Xe) monolayer sliding on a copper (Cu) substrate. We consider thermal-activated creep under a small external lateral force, and observe that slip proceeds by the nucleation and growth of domains in the commensurate interface between the film and the substrate. We measure the activation energy for the nucleation process considering its dependence on the external force, the substrate corrugation, and particle interactions in the film. To understand the results, we use the classical theory of nucleation and compute analytically the activation energy which turns out to be in excellent agreement with numerical results. We discuss the relevance of our results to understand experiments on the sliding of adsorbed monolayers.

Dynamical transitions and sliding friction of the phase-field-crystal model with pinning

We study the nonlinear driven response and sliding friction behavior of the phase-field-crystal (PFC) model with pinning including both thermal fluctuations and inertial effects. The model provides a continuous description of adsorbed layers on a substrate under the action of an external driving force at finite temperatures, allowing for both elastic and plastic deformations. We derive general stochastic dynamical equations for the particle and momentum densities including both thermal fluctuations and inertial effects. The resulting coupled equations for the PFC model are studied numerically. At sufficiently low temperatures, we find that the velocity response of an initially pinned commensurate layer shows hysteresis with dynamical melting and freezing transitions for increasing and decreasing applied forces at different critical values. The main features of the nonlinear response in the PFC model are similar to the results obtained previously with molecular dynamics simulations of particle models for adsorbed layers.

Nonlinear driven response of a phase-field crystal in a periodic pinning potential

We study numerically the phase diagram and the response under a driving force of the phase field crystal model for pinned lattice systems introduced recently for both one- and two-dimensional systems. The model describes the lattice system as a continuous density field in the presence of a periodic pinning potential, allowing for both elastic and plastic deformations of the lattice. We first present results for phase diagrams of the model in the absence of a driving force. The nonlinear response to a driving force on an initially pinned commensurate phase is then studied via overdamped dynamic equations of motion for different values of mismatch and pinning strengths. For large pinning strength the driven depinning transitions are continuous, and the sliding velocity varies with the force from the threshold with power-law exponents in agreement with analytical predictions. Transverse depinning transitions in the moving state are also found in two dimensions. Surprisingly, for sufficiently weak pinning potential we find a discontinuous depinning transition with hysteresis even in one dimension under overdamped dynamics. We also characterize structural changes of the system in some detail close to the depinning transition.

Transition from smooth sliding to stick-slip motion in a single frictional contact

We show that the transition from smooth sliding to stick-slip motion in a single planar frictional junction always takes place at an atomic-scale relative velocity of the substrates.

Boundary lubrication with a glassy interface

Recently introduced constitutive equations for the rheology of dense, disordered materials are investigated in the context of stick-slip experiments in boundary lubrication. The model is based on a generalization of the shear transformation zone (STZ) theory, in which plastic deformation is represented by a population of mesoscopic regions which may undergo nonaffine deformations in response to stress. The generalization we study phenomenologically incorporates the effects of aging and glassy relaxation. Under experimental conditions associated with typical transitions from stick-slip to steady sliding and stop-start tests, these effects can be dominant, although the full STZ description is necessary to account for more complex, chaotic transitions.

Collective stochastic resonance in shear-induced melting of sliding bilayers

The far-from-equilibrium dynamics of two crystalline two-dimensional monolayers driven past each other is studied using Brownian dynamics simulations. While at very high and low driving rates the layers slide past one another retaining their crystalline order, for intermediate range of drives the system alternates irregularly between the crystalline and fluidlike phases. A dynamical phase diagram in the space of interlayer coupling and drive is obtained. A qualitative understanding of this stochastic alternation between the liquidlike and crystalline phases is proposed in terms of a reduced model within which it can be understood as a stochastic resonance for the dynamics of collective order parameter variables. This remarkable example of stochastic resonance in a spatially extended system should be seen in experiments which we propose in the paper.

Transition from stick-slip to smooth sliding: an earthquakelike model

We present a detailed study of an earthquakelike model that exhibits a "transition" from stick-slip motion to smooth sliding at a velocity of the order of those observed in experiments. This contrasts with the many previous microscopic models in which the transition velocity is many orders of magnitude too large. The results show that experimentally observed smooth sliding at the macroscopic scale must correspond to microscopic-scale stick-slip motion.

Driven dynamics: a photodriven Frenkel-Kontorova model

In this study, we examine the dynamics of a one-dimensional Frenkel-Kontorova chain consisting of nanosize clusters (the "particles") and photochromic molecules (the "bonds"), also being subjected to a periodic substrate potential. Whether the whole chain should be running or be locked depends on both the frequency and the wavelength of the light (keeping the other parameters fixed), as observed through numerical simulation. In the locked state, the particles are bound at the bottom of the external potential and vibrate backwards and forwards at a constant amplitude. In the running state, the initially fed energy is transformed into directed motion as a whole. It is of interest to note that the driving energy is introduced to the system by the irradiation of light, and the driven mechanism is based on the dynamical competition between the inherent lengths of the moving object (the chain) and the supporting carrier (the isotropic surface). However, the most important feature is that the light-induced conformational changes of the chromophore lead to the time-and-space dependence of the rest lengths of the bonds.

Friction in a solid lubricant film

Molecular-dynamics study of a thin (one to five layers) lubricant film between two substrates in moving contact are performed using Langevin equations with an external damping coefficient depending on distance and velocity of atoms relative the substrates, motivated by microscopic configurations. They show that the minimal friction coefficient is obtained for the solid-sliding regime. A detailed analysis of the results, the comparison with other microscopic modeling approaches of friction, and the evaluation of quantities that can be compared to experiments, such as the velocity of the transition from stick slip to smooth sliding, are used to discuss the relevance of the microscopic simulations of friction.

Locked-to-running transition in the two-dimensional underdamped driven Frenkel-Kontorova model

We study the nonlinear dc response of a two-dimensional underdamped system of interacting atoms subject to an isotropic periodic external potential with triangular symmetry. We consider various values of the effective elastic constant of the system, two different atomic interaction potentials, and different concentrations of atoms. In the case of a closely packed layer, when its structure is commensurate with the substrate, there is a locked-to-running transition as a function of the driving force, whose mechanism depends on the effective elastic constant. For a low elastic constant, where the layer is weakly coupled, the transition is achieved via the creation of an avalanche of moving particles that leaves a depleted region in its wake. On increasing the effective elastic constant the depleted region becomes less marked and there is a crossover to a scenario in which an island of moving particles nucleates the transition. In the case of a partially filled atomic layer, several dynamical phase transitions between states with different atomic mobility are observed. The mobility of atoms as a function of the external force can vary nonmonotonically with increasing force. For the case of a small external damping, the system can be trapped at a large force in an immobile metastable state, thus demonstrating a "fuse-safety device" on an atomic scale.

Transverse thermal depinning and nonlinear sliding friction of an adsorbed monolayer

We study the response of an adsorbed monolayer under a driving force as a model of sliding friction phenomena between two crystalline surfaces with a boundary lubrication layer. Using Langevin-dynamics simulation, we determine the nonlinear response in the direction transverse to a high symmetry direction along which the layer is already sliding. We find that below a finite transition temperature there exist a critical depinning force and hysteresis effects in the transverse response in the dynamical state when the adlayer is sliding smoothly along the longitudinal direction.

Collective topological dynamics in the frenkel-kontorova chains

Topological dynamics of an array of harmonically coupled damped dc-driven nonlinear oscillators are studied by introducing a dynamical contraction factor and a deviation factor. Different dynamical transitions are identified, and topological changes for these transitions are studied. A bifurcation from the kink state to the kink-antikink-pair state is found, which relates the topological change to the spatiotemporal dynamics of the system. The presence of antikinks leads to the extension of the localized kink, and collisions of kinks and antikinks induce strong oscillations of the topology of the array.

Driven kink in the frenkel-kontorova model

The dynamics of dc driven chain of harmonically interacting atoms in the external sinusoidal potential (the Frenkel-Kontorova model) is studied. It is shown that in the underdamped case the motion of the topological soliton (kink) becomes unstable at a high velocity due to excitation of the localized intrinsic kink mode (the discrete shape mode, or discrete breather) in the kink tail. When the amplitude of the breather's oscillation becomes large enough, it decays into a kink-antikink pair. The subsequent collision of newly created kink and antikink leads to a sharp transition to the running state, where all atoms of the chain slide over the external potential almost freely.

Cellular automaton models of driven diffusive Frenkel-Kontorova-type systems

Three cellular automaton models of increasing complexity are introduced to model driven diffusive systems related to the generalized Frenkel-Kontorova (FK) models recently proposed by Braun et al. [Phys. Rev. E 58, 1311 (1998)]. The models are defined in terms of parallel updating rules. Simulation results are presented for these models. The features are qualitatively similar to those models defined previously in terms of sequentially updating rules. Essential features of the FK model such as phase transitions, jamming due to atoms in the immobile state, and hysteresis in the relationship between the fraction of atoms in the running state and the bias field are captured. Formulating in terms of parallel updating rules has the advantage that the models can be treated analytically by following the time evolution of the occupation on every site of the lattice. Results of this analytical approach are given for the two simpler models. The steady state properties are found by studying the stable fixed points of a closed set of dynamical equations obtained within the approximation of retaining spatial correlations only up to two nearest-neighboring sites. Results are found to be in good agreement with numerical data.

Structural rearrangement of solid surfaces due to competing adsorbate-substrate interactions

Employing a generalized lattice gas theory and the Brownian dynamics simulation, we show that the competing displacive interaction in an adsorbate may cause a continuous distortive transition in the underlying substrate. The threshold for the transition is determined by the competition of the substrate rigidity and the quasielastic energy induced by the adsorbate. In the presence of a strong pinning and repulsive lateral interaction, the resulting structure appears as a compromise between the square lattice of the substrate and the hexagonal arrangement of the adsorbate. For hexagonal substrate lattices the simulation demonstrates that various adsorbate structures (from honeycomb lattices to quasicrystalline pentagonal configurations) may be observed, depending on the effective radii of interaction. Due to the long-ranged coupling the substrate may acquire a substructure induced by the adsorbate. This paper represents a generalization of the work published in Phys. Rev. Lett. 81, 3904 (1998).

From Static to Kinetic Friction in Confined Liquid Films

The transition from rest to sliding contact of atomically smooth solids separated by molecularly thin liquid films was studied. The films could be deformed nearly reversibly to a large fraction of the film thickness. The modulus of elasticity and yield stress were low, considerably less than for a molecular crystal or glass in the bulk. The transition to dissipative sliding was typically (but not always) discontinuous. The dissipative stress was then nearly velocity-independent. The similar response of monolayers strongly attached to the solid surfaces, presenting a well-defined interface for sliding, suggests that the physical mechanism of sliding may involve wall slip.