- 1AC2T Research GmbH, Wiener Neustadt, Austria
- 2Institute of Engineering Design and Product Development, Technische Universität Wien, Vienna, Austria
- 3Material Engineering Division, Toyota Motor Europe NV/SA, Zaventem, Belgium
- 4Institute of Applied Physics, Technische Universität Wien, Vienna, Austria
Among the so-called analytical models of friction, the most popular and widely used one, the Prandtl-Tomlinson model in one and two dimensions is considered here to numerically describe the sliding of the tip within an atomic force microscope over a periodic and atomically flat surface. Because in these PT-models, the Newtonian equations of motion for the AFM-tip are Langevin-type coupled stochastic differential equations the resulting friction and reaction forces must be statistically correctly determined and interpreted. For this, it is firstly shown that the friction and reaction forces as averages of the time-resolved ones over the sliding part, are normally (Gaussian) distributed. Then based on this, an efficient numerical scheme is developed and implemented to accurately estimate the means and standard deviations of friction and reaction forces without performing too many repetitions for the same sliding experiments. The used corrugation potential is the simplest one obtained from the Fourier series expansion of the two-dimensional (2D) periodic potential, e.g., for an fcc(111) surface, which permits sliding on both commensurate and incommensurate paths. In this manner, it is proven that the PT-models predict both frictional regimes, namely the structural superlubricity and stick-slip along (in)commensurate sliding paths, if the ratio of mean corrugation and elastic energies is properly set.
1 Introduction
Frictional experiments in an atomic force microscope can be straightforwardly simulated by using either molecular dynamics (MD) or ab-initio (first-principles) methods, see for example, Refs. (Kobayashi et al., 2016; Wolloch et al., 2015). From a computational point of view, however, these numerical calculations are extremely demanding and time consuming. Therefore, the so-called one-dimensional (1D) and two-dimensional (2D) analytical models of friction, such as the Prandtl-Tomlinson and Frenkel-Kontorova models, (Dong et al., 2011) are still remaining to be highly efficient alternatives to the atomistic models and hence very popular because of their simplicity with which they are capturing the main undergoing mechanisms within nanotribological AFM-experiments. (Pawlak et al., 2016)
The Prandtl-Tomlinson model (PT-model), in its original formulation, (Prandtl, 1928) was introduced by Prandtl while developing a kinetic theory of solids assuming a periodic interaction between the bodies. Curiously, in the other paper ascribed to the PT-model, (Tomlinson, 1929) Tomlinson was not dealing at all with these concepts, since he was outlining a molecular theory of friction based on the findings by Lennard-Jones. Despite this historical inaccuracy, (Popov and Gray, 2012) within a 1D or 2D PT-model, one can clearly distinguish between various frictional regimes by considering a single parameter, namely the ratio between the average interaction and elastic energies. (Socoliuc et al., 2004)
Traditionally, if the coefficient-of-friction (CoF) observed for a tribological system is less than or equal to 0.01 that system is said to evolve within a superlubric frictional regime. (Martin and Erdemir, 2018) Superlubricity was predicted theoretically in 1990 by Hirano and Shinjo, (Hirano and Shinjo, 1990) and its occurrence was shown experimentally 1 year later. (Hirano et al., 1991) Nowadays, beyond this structural superlubricity discovered in the nineties as caused by the incommensurability of interacting surfaces, one speaks also about thermolubricity (Krylov et al., 2005) and chemolubricity (Erdemir and Martin, 2007) which are two other forms of superlubricity caused by the temperature acting as lubricant and due to chemistry, respectively. For further information on superlubricity, the reader is recommended to consult the excellent reviews in the newest book edited by J.-M. Martin and A. Erdemir. (Erdemir et al., 2021)
In this contribution, the structural superlubricity will be addressed numerically by using the most popular and hence widely applied analytical model of friction, namely the PT-model in one-dimension and two-dimensions, in its well-known form from the kinetic theory of solids and solid mechanics, see Section 2. For this, the periodic corrugation of interest is explicitly derived in Section 4 and directly introduced into the corresponding 1D/2D PT-model. Since the resulting Newtonian equations of motion are forming a set of stochastic differential equations (SDEs), they are solved numerically by applying an adequate fourth-order Runge-Kutta method as presented in Section 3. Due to the randomness of the thermally induced force mimicking the impact of temperature on the sliding, both frictional regimes, either stick-slip or superlubric, must be statistically evaluated and interpreted, for example, in terms of distributions as shown in Secion 5. Beyond the development of a statistically proper numerical scheme for comparing the calculated frictional performance along various paths on periodic atomic surfaces, the ultimate goal of this contribution is to elucidate in which circumstances, i.e., range of various parameters such as materials, elastic coupling of AFM-tip to the drive, sliding velocity or absolute temperature, one could detect the superlubric frictional regime with high accuracy, see Section 6.
2 Prandtl-Tomlinson Models
In an atomic force microscope, the movement of the tip on a corrugated surface of a substrate in dry contact conditions is described within the kinetic theory of solids by the Langevin equation of motion, (Filippov et al., 2010)
with r(t) being the time-dependent position of the AFM-tip of mass m and γ is characterizing the damping proportional to the velocity of the AFM-tip. Here the first force on the left-hand side is that which accelerates the AFM-tip, the second one acts against the former and it is proportional to the velocity of the AFM-tip, whereas the last, third force on the left-hand side is due to the interaction of the AFM-tip with the surface and with the drive.
If the total potential energy consists of the time independent corrugation potential U(r) and the elastic potential of the tip-drive subsystem, respectively, then
where kc stands for the effective stiffness of the linear coupling (e.g., via spring) between the AFM-tip and drive. Note that the drive is supposed to move with a constant sliding velocity vs, such that
with Φ defining the orientation of the sliding path with respect to the Cartesian x-axis, whereas i and j being the in-plane unity vectors along the Cartesian axes. The relative direction of vs also sets the commensurability/incommensurability of the sliding path versus the periodicity of corrugation U(r), for more details see Secion 4 below. For example, for a 2D corrugation one can consider the same functional form as for the 1D, such that the latter 1D corrugation provides a periodic atomic chain according to Φ, whereas in the 2D case it represents the simplest 2D periodic potential corresponding to the crystalline orientation of substrate’s surface.
In Eq. 1, a normal (Gaussian) distribution is assumed for the amplitude of the thermal-noise-induced random force ζ(t), thus for its ensemble averaged auto-correlation function holds that
where δ denotes the Dirac delta-function, kB is the Boltzmann constant and T stands for the absolute temperature, and its random orientation is given by
with θ being a uniformly distributed random angle.
3 Runge-Kutta Methods
In the following, some details on the implementation are given, which can be directly used for a 1D PT-model as well as for a 2D PT-model where ξ(t) is measured along the sliding path. In this case, the second-order stochastic differential equation in Eq. 1 is equivalent to a system of two coupled first-order differential equations,
which, once numerically solved, yields simultaneously the position ξ(t) and velocity v(t) of the AFM-tip along the sliding path for ∀t ≥ t0, if the initial position ξ0 and velocity v0 of the AFM-tip are both known at the starting time t0, i.e.,
The Runge-Kutta (RK) methods for solving differential equations numerically were independently developed by Runge (Runge, 1895) and Kutta (Kutta, 1901) more than a century ago. During their long history, (Butcher and Wanner, 1996) they were continuously improved, (Kalogiratou et al., 2014) such that also RK-methods exist for solving SDEs, one of which will be here applied as follows. The fourth-order Runge-Kutta method applied for the coupled SDEs in Eq. 6, i.e., when T > 0, means that one knowing ξn = ξ(tn) and vn = v(tn) for tn = n Δt (∀n = 0, 1, 2, …), computes for i = 1, 2, 3, 4
to finally get
where all the involved real-valued coefficients as derived in Ref. (Kasdin, 1995), Gi (i = 1, …, 4) are Gaussian random numbers and W = 2mγ kBT in accordance with Eq. 5. Due to this latter randomness at finite absolute temperature, Eq. 6 admits not only a single exact solution, but a set of differently valued solutions which, however, are all statistically equivalent down to T = 1 mK–as we found. On the other hand, in accordance with Eq. 5, the thermal-noise-induced random force ζ(t) is vanishing at absolute zero temperature, since W = 0 when T = 0 K, and hence the first-order SDE in Eq. 6 turns at T = 0 K into an ordinary differential equation (ODE) and therefore the ‘classical’ fourth-order Runge-Kutta method for ODEs can be immediately applied. Formally, this means that one considers the former expressions with wi = 0 (i = 1, …, 4) and the following tableaux of Runge-Kutta coefficients, (Press et al., 2007)
4 Corrugation Potential and Frictional Regimes
In the following, the corrugation potential of fcc(111)–the conceptual model considered in this contribution–is derived from the Fourier series expansion of the 2D translational invariant potential. For this, consider the 2D primitive lattice vectors for an fcc(111) crystalline surface as being
with
and hence the corresponding translation vector within the 2D reciprocal space reads Km = m1b1 + m2b2 (
Now, from the Fourier series expansion of the 2D periodic potential U(r) = U(r + Tn),
written for
can be immediately considered as a proper model of the interaction potential between the AFM-tip and fcc(111) surface, e.g., by assuming for sake of simplicity that U±10 = U0±1 = U00. In this manner, one could have any desired sinusoidal corrugation as defined by φ, for example,
providing the temperature independent force acting on the AFM-tip in the form of
Indeed, when φ = 0, the most applied and simplest sinusoidal corrugation of period equal to the lattice constant a results,
Considering in these latter expressions ξ = Γ cos Φ, via Φ one can easily select any commensurate/incommensurate path for sliding, by letting Φ either to show or not to point along a high-symmetry direction on fcc(111), for example, when Φ = 0 or Φ = π/12 = 15◦.
The values for all the parameters entering Eqs 1–4, 13 are from Ref. (Dong et al., 2011), e.g., a = 2.88 Å, and U0 = 0.6 eV, and hence correspond to fcc Au(111). However, for the large-parametric study presented in Seciton 6, both U0 and kc were continuously varied to determine the impact of the substrate’s material and that of the coupling strength on the friction. When kc is kept unaltered, then its constant value is set to 1 N/m. As already mentioned, it turned out that independently of the sliding path, whether it is commensurate or incommensurate, the realization of a normal stick-slip as illustrated in Figure 1, or a superlubric frictional regime provided by Figure 2 indeed is only depending on a single quantity,
namely on the ratio between the mean interaction energy of the AFM-tip with the surface and the elastic energy of spring coupling the AFM-tip to the drive, (Socoliuc et al., 2004) see also Section 6.
FIGURE 1. Stick-slip movement of the AFM-tip at RT of 300 K along an incommensurate sliding path Φ = 15 deg for an η set accordingly, when U0 is kept constant and kc is adjusted, recall Eqs 3, 14. The trajectories of drive (black) and AFM-tip (green) are shown in the top panel, whereas the friction and reaction forces (green) at each position of the drive are given in the middle and bottom panels together with their means (blue) over the entire sliding path.
FIGURE 2. As in Figure 1, for an η-value leading to a superlubric sliding of the AFM-tip on an fcc(111) substrate.
5 Statistics of Forces
From Eq. 13 it becomes immediately clear that the average of the lateral force
is a proper measure for the friction force and completely in accordance with the experimental practice. (Schwarz and Hölscher, 2016; Peng et al., 2020) Note, however, that accordingly the calculated friction force could be either positive or negative valued, depending on whether the drive is more times in advance with respect to the AFM-tip or the other way around during the entire sliding period. In this Eq. 15,
either frictional or reaction one (the latter only within the 2D PT-model), and its standard deviation
are both calculated using numerical results from n = 104 times repeated “identical” sliding experiments, since it is assumed that after such an amount of repetitions
As can be directly seen from Figure 3A, the distribution of time-averaged friction force within a PT-model is a normal (Gaussian) one,
of an almost perfect shape in view of its third (skewness) and fourth (kurtosis) central moments,
FIGURE 3. In the top panel, the histogram of the time-averaged friction force (blue) for the same sliding experiment repeated 104 times together with the corresponding normal (Gaussian) distribution (solid red line) and the mean (dashed red line). Here the skewness (sk) and kurtosis (ku), i.e., the third and fourth central moments, are measuring the high quality of the normal distribution for the time-averaged friction force. This is additionally also quantified by Pearson’s linear correlation coefficient r2 in the bottom panel, where the probability plot of the time-averaged friction force is given.
Furthermore, by inspecting the normal probability plot given in Figure 3B which per definition identifies graphically the departure of random data from normality, i.e., from the Gaussian behavior, it can be observed that the ordered time-averaged friction force shows against the theoretical quantiles qi (i = 1, …, n − 1) indeed a perfect linear correlation in terms of the Pearson’s coefficient,
where
FIGURE 4. As in Figure 3, for the reaction force.
Since our main goal is to numerically study the occurrence of structural superlubricity as function of various parameters, such as temperature and velocity, for example, a 104 times repetition of the sliding for a single parametric configuration is too demanding computationally and therefore one should develop a numerical scheme to achieve accurate approximations of
In a next attempt, it was iteratively attempted to achieve the convergence of the mean force. 1) In a first step, one increases step-by-step the number of repetitions n until the mean
Finally, it was decided to realize the convergence of the mean friction and/or reaction force by setting the accuracy ϵG > 0 with which the normality of forces should be achieved. Namely, increasing either step-by-step or blockwise the number of repetitions, it is continuously verified whether
If these criteria are satisfied, the mean of friction and reaction forces are considered to be converged and the normality of their distributions is additionally quantified by Pearson’s linear correlation coefficient, recall Eq. 21. For the case that one of the criteria in Eq. 22 and/or
FIGURE 5. Convergence of the friction force skewness (blue) and kurtosis (green) within the 1D PT-model, together with the number n of repetitions needed (red dashed line) to achieve a given accuracy ɛG, and to this accuracy corresponding Pearson’s linear correlation coefficient
FIGURE 6. As in Figure 5, for the friction (top row) and reaction force (bottom row) within a stick-slip (left panels) and a superlubric frictional regime (right panels) obtained by applying the 2D PT-model.
6 Large-Parametric Studies
Having fixed the convergence criteria with ϵG = 0.02, large parametric studies of sliding along (in)commensurate paths were performed first using the 1D PT-model. By varying the system and environment properties, we aim to gain understanding on the dependence of the friction forces and thus the frictional regime on any of these parameters. As illustrated in Figure 7, the temperature values of T = 0 K and T = 300 K were considered when solving the corresponding ODEs and SDEs respectively. In these calculations, the η-dependence of structural superlubricity was studied by keeping all parameters entering Eq. 1 constant, except the corrugation potential amplitude U0 and the spring constant of the coupling of the AFM-tip to the drive kc, respectively, which were varied independently. A different behaviour is observed in the case of each variation.
FIGURE 7. η-dependence of structural superlubricity as predicted by the 1D PT-model at 0 K (solid line) and at RT of 300 K (full circles) along various incommensurate sliding paths Φ = 5, 10, 15, 20, 25 and 30 deg (depicted in violet as specified in the legend), recall Figures 1, 2, when either U0(top row) or kc(bottom row) is kept constant in Eq. 14.
When U0 is kept constant, one can observe a short region of very low friction forces, followed by a quick increase with a diminishing slope as the ratio η increases. The introduction of thermal noise at non-zero temperatures leads to the reduction of the friction forces in the stick-slip regime, however, the region of the superlubric regime is not significantly affected. At very low η values, which correspond to very high spring constants kc, negative friction force values of low magnitude are observed. This can be due to the non-convergence of the average force, as it was already observed that a larger number of statistically independent realisations is typically required in the superlubric regime. On the other hand, it could pinpoint to an incompatibility of the very high elastic coupling energy for the given interaction energy of the AFM-tip with the surface, where the assumptions of the PT-model may not be valid.
In the opposite case where the spring constant kc is kept fixed and the interaction potential is varied, only positive friction forces are observed. The superlubric regime holds for larger values of η, while the increase of forces is less abrupt and shows a linear trend. Here, the introduction of thermal noise also leads to a reduction of the friction forces, however, it also leads to a significant expansion of the range of η values which lead to a superlubric regime. It becomes clear, therefore, from the analysis of Figure 7 that the value of η uniquely determines the type of frictional regime, either stick-slip or superlubric between the AFM-tip and any semi-infinite crystalline substrate.
Note that all these numerical studies evidence significantly larger intervals for η where the superlubric frictional regime could occur than that of η < 1 assumed in the literature for both U0 and kc kept constant, recall for example Ref. (Socoliuc et al., 2004). This means that by varying the material and system properties, the occurrence of superlubricity could be further tuned and even extended on (in)commensurate sliding paths. Furthermore, it seems also that the superlubricity for η < 1, recall Figure 7, it is more pronounced than over its extension for 1 ≤ η < 20, probably because the former is energetically and the latter symmetrically more favorable. Even small negative values for friction force are obtained when η < 1 and U0 is varied, which shows that probably a sliding distance of only 3.5 nm which was considered here, is not enough to get in average the drive in front of the AFM-tip.
All the same holds for the friction forces calculated using the 2D PT-model, cf. Figure 7 with the left panels of Figure 8. What concerns the η-dependence of the reaction force, this differs not too much from that of the friction force, as one observes on the right column of Figure 8, and shows that keeping the sliding on an incommensurate path needs in magnitude higher reaction forces, when η increases. Interestingly, the incommensurate paths are not equivalent, since the η-dependence of the friction and reaction forces is Φ-dependent too. Furthermore, both forces, i.e., the frictional in both 1D and 2D PT-models, and the reaction force within the 2D PT-model, feature a strong temperature dependence for a given sliding velocity.
FIGURE 8. As in Figure 7, for three two-dimensional incommensurate sliding paths Φ = 5, 10 and 15 deg, when using the 2D PT-model.
Apart from identifying the frictional regime according to the material and system properties, the 1D and 2D PT-models also enable us to study the temperature and velocity dependence of friction and reaction forces in both stick-slip and superlubric frictional regimes. In Figures 9, 10 only the calculated friction force for a large range of temperature and sliding velocity values is given, because in magnitude much smaller reaction force behaves pretty similarly. It can be readily observed that the increase of temperature leads to a monotonic decrease of friction forces within both frictional regimes, as the increasing thermal motion aids the AFM-tip to overcome the interaction energy barriers and slide with a lower mechanical energy. On the other hand, the increase of sliding velocity results in increased friction forces strongly depending on the frictional regime, probably due to different impact of the damping. A linear dependence of the friction and reaction force (although the latter is not shown explicitly here) on the logarithm of the sliding velocity can be observed within the stick-slip regime,
of which slope a1kBT + a0 is linearly depending on the absolute temperature T. In this Eq. 23kB stands for the Boltzmann constant and
provides the temperature and sliding velocity dependence of the friction and reaction forces. For the values of the polynomial coefficients ap (p = 0, 1, 2) and further computational details, see the Supplementary Material. All these are suggesting at least a power-law dependence of friction reaction force on the sliding velocity with a temperature-dependent exponent.
FIGURE 9. Temperature and sliding velocity dependence of the friction force as predicted by the 2D PT-model within the stick-slip regime along the incommensurate path for Φ = 15 deg, when either U0(top row) or kc(bottom row) is kept constant in Eq. 14.
FIGURE 10. As in Figure 9, but within the superlubric frictional regime.
7 Summary and Conclusions
In summary, the PT-model was considered in one and two dimensions for the numerical description of an AFM-tip sliding over a periodic and atomically flat surface. A system of SDEs was considered and implemented for the inclusion of the effects of temperature on the motion of the tip. The conditions for the derivation of statistically accurate friction and reaction forces in a numerically efficient way were defined. The proposed statistically proper scheme can be applied for sliding on both commensurate and incommensurate paths, while it allows the study of the impact of both system and environmental conditions on the resulting friction forces.
By the numerical study of the AFM-tip sliding using the PT-models, it was concluded that the ratio between the mean interaction energy of the AFM-tip with the surface and the elastic energy of spring coupling the AFM-tip to the drive, can uniquely determine the type of frictional regime, either stick-slip or superlubric between the AFM-tip and any semi-infinite crystalline substrate. In addition, the friction force in both regimes is predicted to be a function of temperature and sliding velocity.
All these findings can support the design of appropriate experimental studies for the realization of structural superlubricity, while they can also shed light into the design of novel systems showing ultra-low friction which can be beneficial in several engineering applications.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
MS: made the calculations and representation of data KG: wrote the manuscript and interpretation of data CG: interpretation of data and discussions AV: developed the theory, made the programming, wrote the manuscript.
Funding
This work was funded by the Austrian COMET-Project INTRIBOLOGY (FFG-No. 872176) and carried out at the “Excellence Centre of Tribology” (AC2T research GmbH). MS, KG, and AV have received funding also from the European Union’s Horizon 2020 research and innovation prgramme under grant agreement No. 814494, project i-TRIBOMAT.
Conflict of Interest
Authors MS and AV were employed by company AC2T Research GmbH, whereas author KG was employed by Toyota Motor Europe NV/SA.
The remaining author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors acknowledge TU Wien Bibliothek too for financial support through its Open Access Funding Program.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmech.2021.742684/full#supplementary-material
References
Butcher, J. C., and Wanner, G. (1996). Runge-kutta Methods: Some Historical Notes. Appl. Numer. Maths. 22, 113–151. doi:10.1016/s0168-9274(96)00048-7
Dong, Y., Vadakkepatt, A., and Martini, A. (2011). Analytical Models for Atomic Friction. Tribol. Lett. 44, 367–386. doi:10.1007/s11249-011-9850-2
Filippov, A. E., Vanossi, A., and Urbakh, M. (2010). Origin of Friction Anisotropy on a Quasicrystal Surface. Phys. Rev. Lett. 104, 074302. doi:10.1103/PhysRevLett.104.074302
Hirano, M., and Shinjo, K. (1990). Atomistic Locking and Friction. Phys. Rev. B 41, 11837–11851. doi:10.1103/physrevb.41.11837
Hirano, M., Shinjo, K., Kaneko, R., and Murata, Y. (1991). Anisotropy of Frictional Forces in Muscovite Mica. Phys. Rev. Lett. 67, 2642–2645. doi:10.1103/physrevlett.67.2642
Kalogiratou, Z., Monovasilis, T., Psihoyios, G., and Simos, T. E. (2014). Runge-kutta Type Methods with Special Properties for the Numerical Integration of Ordinary Differential Equations. Phys. Rep. 536, 75–146. doi:10.1016/j.physrep.2013.11.003
Kasdin, N. J. (1995). Runge-kutta Algorithm for the Numerical Integration of Stochastic Differential Equations. J. Guidance, Control Dyn. 18 (1), 114–120. doi:10.2514/3.56665
Kobayashi, K., Liang, Y., Amano, K.-i., Murata, S., Matsuoka, T., Takahashi, S., et al. (2016). Molecular Dynamics Simulation of Atomic Force Microscopy at the Water-Muscovite Interface: Hydration Layer Structure and Force Analysis. Langmuir 32 (15), 3608–3616. doi:10.1021/acs.langmuir.5b04277
Krylov, S. Y., Jinesh, K. B., Valk, H., Dienwiebel, M., and Frenken, J. W. (2005). Thermally Induced Suppression of Friction at the Atomic Scale. Phys. Rev. E Stat. Nonlin Soft Matter Phys. 71, 065101. doi:10.1103/PhysRevE.71.065101
Kutta, W. (1901). Beitrag zur näherungsweisen Integration totaler Differentialgleichungen. Z. Math. Phys. 46, 435.
Martin, J. M., and Erdemir, A. (2018). Superlubricity: Friction's Vanishing Act. Phys. Today 71 (4), 40–46. doi:10.1063/pt.3.3897
Pawlak, R., Ouyang, W., Filippov, A. E., Kalikhman-Razvozov, L., Kawai, S., Glatzel, T., et al. (2016). Single-Molecule Tribology: Force Microscopy Manipulation of a Porphyrin Derivative on a Copper Surface. ACS Nano 10 (1), 713–722. doi:10.1021/acsnano.5b05761
Peng, Y., Zeng, X., Yu, K., and Lang, H. (2020). Deformation Induced Atomic-Scale Frictional Characteristics of Atomically Thin Two-Dimensional Materials. Carbon 163, 186–196. doi:10.1016/j.carbon.2020.03.024
Popov, V. L., and Gray, J. A. T. (2012). Prandtl-Tomlinson Model: History and Applications in Friction, Plasticity, and Nanotechnologies. Z. Angew. Math. Mech. 92 (9), 683–708. doi:10.1002/zamm.201200097
Prandtl, L. (1928). Ein Gedankenmodell zur kinetischen Theorie der festen Körper. Z. Angew. Math. Mech. 8, 85–106. doi:10.1002/zamm.19280080202
Press, W. H., Teukolsky, S. A., Vetterling, W. T., and Flannery, B. P. (2007). Numerical Recipes – the Art of Scientific Computing. Cambridge: University Press.
Runge, C. (1895). Über die numerische Auflösung von Differentialgleichungen. Math. Ann. 46, 167–178. doi:10.1007/bf01446807
Schwarz, U. D., and Hölscher, H. (2016). Exploring and Explaining Friction with the Prandtl-Tomlinson Model. ACS Nano 10 (1), 38–41. doi:10.1021/acsnano.5b08251
Socoliuc, A., Bennewitz, R., Gnecco, E., and Meyer, E. (2004). Transition from Stick-Slip to Continuous Sliding in Atomic Friction: Entering a New Regime of Ultralow Friction. Phys. Rev. Lett. 92, 134301. doi:10.1103/physrevlett.92.134301
Tomlinson, G. A. (1929). CVI.A Molecular Theory of Friction. Lond. Edinb. Dublin Philos. Mag. J. Sci. 7, 905–939. doi:10.1080/14786440608564819
Keywords: Prandtl-Tomlinson model, Langevin dynamical systems, statistics of forces, stick-slip, superlubricity
Citation: Srbulovic M, Gkagkas K, Gachot C and Vernes A (2021) Statistics of Sliding on Periodic and Atomically Flat Surfaces. Front. Mech. Eng 7:742684. doi: 10.3389/fmech.2021.742684
Received: 16 July 2021; Accepted: 15 September 2021;
Published: 06 October 2021.
Edited by:
Dirk Spaltmann, Federal Institute for Materials Research and Testing (BAM), GermanyReviewed by:
Yitian Peng, Donghua University, ChinaFeodor M. Borodich, Cardiff University, United Kingdom
Copyright © 2021 Srbulovic, Gkagkas, Gachot and Vernes. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: András Vernes, vernes@ac2t.at