Tribology of Sliding Elastic Media

Anomalous phase diagram of the elastic interface with nonlocal hydrodynamic interactions in the presence of quenched disorder

We investigate the influence of quenched disorder on the steady states of driven systems of the elastic interface with nonlocal hydrodynamic interactions. The generalized elastic model (GEM), which has been used to characterize numerous physical systems such as polymers, membranes, single-file systems, rough interfaces, and fluctuating surfaces, is a standard approach to studying the dynamics of elastic interfaces with nonlocal hydrodynamic interactions. The criticality and phase transition of the quenched generalized elastic model are investigated numerically and the results are presented in a phase diagram spanned by two tuning parameters. We demonstrate that in the one-dimensional disordered driven GEM, three qualitatively different behavior regimes are possible with a proper specification of the order parameter (mean velocity) for this system. In the vanishing order parameter regime, the steady-state order parameter approaches zero in the thermodynamic limit. A system with a nonzero mean velocity can be in either the continuous regime, which is characterized by a second-order phase transition, or the discontinuous regime, which is characterized by a first-order phase transition. The focus of this research is to investigate the critical scaling features near the pinning-depinning threshold. The behavior of the quenched generalized elastic model at the critical depinning force is explored. Near the depinning threshold, the critical exponent is obtained numerically.

Detecting depinning and nonequilibrium transitions with unsupervised machine learning

Using numerical simulations of a model disk system, we demonstrate that a machine learning generated order-parameter-like measure can detect depinning transitions and different dynamic flow phases in systems driven far from equilibrium. We specifically consider monodisperse passive disks with short range interactions undergoing a depinning phase transition when driven over quenched disorder. The machine learning derived order-parameter-like measure identifies the depinning transition as well as different dynamical regimes, such as the transition from a flowing liquid to a phase separated liquid-solid state that is not readily distinguished with traditional measures such as velocity-force curves or Voronoi tessellation. The order-parameter-like measure also shows markedly distinct behavior in the limit of high density where jamming effects occur. Our results should be general to the broad class of particle-based systems that exhibit depinning transitions and nonequilibrium phase transitions.

Pinned-to-sliding transition and structural crossovers for helically confined charges

We explore the nonequilibrium dissipative dynamics of a system of identical charged particles trapped on a closed helix. The particles are subject to an external force accelerating them along the underlying structure. The effective interactions between the charges induce a coupling of the center of mass to the relative motion which in turn gives rise to a pinned-to-sliding transition with increasing magnitude of the external force. In the sliding regime we observe an Ohmic behavior signified by a constant mobility. Within the same regime a structural transition of the helical particle chain takes place with increasing the helix radius leading to a global change of the crystalline arrangement. The resulting crystal is characterized by the existence of multiple defects whose number increases with the helix radius.

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.

Seismic cycles, size of the largest events, and the avalanche size distribution in a model of seismicity

We address several questions on the behavior of a numerical model recently introduced to study seismic phenomena, which includes relaxation in the plates as a key ingredient. First, we make an analysis of the scaling of the largest events with system size and show that, when parameters are appropriately interpreted, the typical size of the largest events scale as the system size, without the necessity to tune any parameter. Second, we show that the temporal activity in the model is inherently nonstationary and obtain from here justification and support for the concept of a "seismic cycle" in the temporal evolution of seismic activity. Finally, we ask for the reasons that make the model display a realistic value of the decaying exponent b in the Gutenberg-Richter law for the avalanche size distribution. We explain why relaxation induces a systematic increase in b from its value b≃0.4 observed in the absence of relaxation. However, we have not been able to justify the actual robustness of the model in displaying a consistent b value around the experimentally observed value b≃1.

Friction phenomena and phase transition in the underdamped two-dimensional Frenkel-Kontorova model

Locked-to-sliding phase transition has been studied in the driven two-dimensional Frenkel-Kontorova model with the square symmetric substrate potential. It is found that as the driving force increases, the system transfers from the locked state to the sliding state where the motion of particles is in the direction different from that of driving force. With the further increase in driving force, at some critical value, the particles start to move in the direction of driving force. These two critical forces, the static friction or depinning force, and the kinetic friction force for which particles move in the direction of driving force have been analyzed for different system parameters. Different scenarios of phase transitions have been examined and dynamical phases are classified. In the case of zero misfit angle, the analytical expressions for static and kinetic friction force have been obtained.

Coarse-graining the dynamics of a driven interface in the presence of mobile impurities: effective description via diffusion maps

Developing effective descriptions of the microscopic dynamics of many physical phenomena can both dramatically enhance their computational exploration and lead to a more fundamental understanding of the underlying physics. Previously, an effective description of a driven interface in the presence of mobile impurities, based on an Ising variant model and a single empirical coarse variable, was partially successful [M. Haataja, Phys. Rev. Lett. 92, 160603 (2004)]; yet it underlined the necessity of selecting additional coarse variables in certain parameter regimes. In this paper we use a data mining approach to help identify the coarse variables required. We discuss the implementation of this diffusion map approach, the selection of a similarity measure between system snapshots required in the approach, and the correspondence between empirically selected and automatically detected coarse variables. We conclude by illustrating the use of the diffusion map variables in assisting the atomistic simulations and we discuss the translation of information between fine and coarse descriptions using lifting and restriction operators.

Dry friction avalanches: experiment and theory

Experimental evidence and theoretical models are presented supporting the conjecture that dry friction stick-slip is described by self-organized criticality. We use the data, obtained with a pin-on-disk tribometer set to measure lateral force, to examine the variation of the friction force as a function of time. We study nominally flat surfaces of matching aluminum and steel. The probability distribution of force drops follows a negative power law with exponents mu in the range 3.2-3.5. The frequency power spectrum follows a 1/f alpha pattern with alpha in the range 1-1.8. We first compare these experimental results with the well-known Robin Hood model of self-organized criticality. We find good agreement between theory and experiment for the force-drop distribution but not for the power spectrum. We explain this on a physical basis and propose a model which takes explicitly into account the stiffness and inertia of the tribometer. Specifically, we numerically solve the equation of motion of a block on a friction surface pulled by a spring and show that for certain spring constants the motion is characterized by the same power law spectrum as in experiments. We propose a physical picture relating the fluctuations of the force drops to the microscopic geometry of the surface.

Dynamic phases 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. When driving force increases, the system transfers from a disorder locked state to an ordered sliding state corresponding to a moving crystal. By varying the values of the effective elastic constant, damping, and temperature, we found different scenarios and intermediate phases during the ordering transition. For a soft atomic layer, the system passes through a plastic-channel regime that appears as a steady-state regime at higher values of the damping coefficient. For high values of the effective elastic constant, when the atomic layer is stiff, the intermediate plastic phase corresponds to a traffic-jam regime with immobile islands in the sea of running atoms. At a high driving of the stiff layer, a solitonlike elastic flow of atoms has been observed.

Apparent hysteresis in a driven system with self-organized drag

The motion of extended defects in materials is often resisted by interaction with diffusing impurities. The defect undergoes a transition from slow to fast migration, or vice versa. This transition is commonly hysteretic, with the defect jumping back and forth between these kinetic states. We explore such hysteresis within a kinetic Monte Carlo simulation. After identifying the slow variables, we construct an effective potential that quantitatively describes the stable and metastable states of the system (a bifurcation diagram) and the kinetics of the transitions between these. This description provides a means to determine whether hysteresis will be observed in a particular simulation and the detailed nature of this hysteresis (switching times).

Motional ordering of a charge-density wave in the sliding state

We have used high-resolution x-ray scattering, in the presence of an applied direct current, for studying the correlation lengths in the sliding charge-density wave (CDW) state. Transport properties were simultaneously measured in situ during the experiment. We find that, while the transverse correlation is reduced when the CDW moves, the CDW becomes more ordered in the direction of motion. This is the first report of a motional ordering process in a periodic system other than a vortex lattice.