Pith. sign in

REVIEW 2 major objections 5 minor 8 references

Effect of crystallinity on the frictional and wear performance of molybdenum disulfide: A molecular dynamics study

T0 review · 2 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Crystalline MoS2 slides with near-zero friction; friction climbs steeply as structural order falls.

desk verdict Graded crystallinity MD series gives a clean friction trend, but the wear claim rests on a single run and contradicts the authors' own cited experiment. read the letter →

arxiv 2506.00690 v1 pith:5ECEBTUF submitted 2025-05-31 cond-mat.mtrl-sci

Abhiram B R , Ilia Ponomarev , Tomas Polcar This is my paper
classification cond-mat.mtrl-sci
keywords molybdenumdisulfidecrystallinityfrictionwearmoleculardynamicssuperlubricityreactiveforcefieldsolidlubricant
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper uses reactive molecular dynamics to ask whether the internal structural order of MoS2, not just its surface chemistry, controls how it behaves as a solid lubricant. It claims that as crystallinity is reduced, the coefficient of friction rises monotonically from 0.00014 for the crystal to 0.28 for the amorphous form, and that wear depth under a 24.85 nN load grows from 5.32 Å to 14.8 Å. The proposed cause is that a crystalline surface stays atomically flat and lets the slider take a zig-zag minimum-energy path, while disordered surfaces introduce roughness and dissipate energy. If true, crystallinity becomes a direct design lever for ultralow-friction, wear-resistant MoS2 coatings, and superlubricity can occur even when a rigid slider, not MoS2 itself, moves across the material.

What carries the argument

The machinery is a suite of reactive molecular dynamics simulations in which a spherical rigid carbon body slides over MoS2 samples with controlled crystallinity (crystalline, polycrystalline PC75/PC25/PC10, and amorphous) under normal loads from 1.49 to 9.94 nN for friction, and 24.85 nN for wear. Crystallinity is quantified by the percentage of Mo atoms whose second-neighbor shell forms a hexagonal 2H pattern (Appendix A), and friction is read from the slope of frictional force versus normal load across three independent runs. The explanatory mechanisms are surface flatness, assessed by a Gaussian density mesh that maps out-of-plane roughness, and frictional anisotropy, inferred from the observed lateral zig-zag motion of the slider.

What would settle it

Take a highly crystalline and a fully amorphous MoS2 film and measure their coefficients of friction and wear depths with an AFM tip under identical load, speed, and dry environment; if the amorphous film does not show a substantially higher COF and deeper wear than the crystalline film, the claimed ordering fails. A complementary computational test is to rerun the same sliding protocol with an independently fitted reactive or machine-learned potential and check whether the COF ordering across crystallinity is preserved.

Watch

Extended reading notes

Core claim

On its own terms, the paper establishes a monotonic relationship: lower crystallinity means higher friction and lower wear resistance in MoS2. Crystalline MoS2 reaches a coefficient of friction of 0.00014±0.00005 under a rigid carbon slider, which qualifies as superlubricity, whereas polycrystalline samples with 75, 25, and 10 percent crystallinity give 0.05, 0.11, and 0.17, and amorphous MoS2 gives 0.28. The explanation the authors offer is twofold: the crystalline surface remains flat under load, minimizing mechanical interlocking, and frictional anisotropy lets the slider follow a lateral zig-zag path of least resistance rather than a straight line. The same structural order is claimed to protect against wear, since crystalline MoS2 shows a wear depth of only 5.32 Å compared with 14.8 Å for amorphous MoS2.

Load-bearing premise

The entire friction and wear comparison rests on the Mo-S-C reactive force field being accurate for sliding contact, but the paper validates it only against the energetics of crystalline and amorphous phases, not against any friction or wear measurement.

Editorial extensions

If this is right

  • Crystalline MoS2 can exhibit superlubricity against a rigid carbon slider, not only in self-mated MoS2 contacts, broadening the conditions under which superlubricity is expected.
  • The coefficient of friction is continuously tunable by crystallinity, from roughly 0.00014 for the crystal to 0.28 for the amorphous form, offering a structural route to control friction.
  • Wear resistance follows the same order: higher crystallinity gives shallower wear tracks under the same high load, so coating design should preserve crystalline order for durability.
  • The combination of flat surface and frictional anisotropy provides a concrete atomic-scale picture of why ordered layered materials slide easily, which can guide the search for other ultralow-friction solids.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The three-order spread in COF across crystallinity suggests that even a small amorphous fraction creates local energy barriers; a testable extension is to map COF as a function of crystallinity percentage with finer resolution than the 75/25/10 points used here.
  • The rigid-body slider assumption neglects deformation of the counterbody; allowing the carbon tip to be flexible might change the superlubricity value, since deformable contacts can increase contact area and energy dissipation.
  • The wear direction here runs opposite to the experimental comparison cited in the paper's own introduction (amorphous MoS2 reported to wear four times better than crystalline), so a direct AFM wear experiment across controlled crystallinities would settle which trend is physical.
  • The reported COF of 0.00014 is close to the numerical floor of the simulation; verifying it with a different reactive potential or a DFT-informed machine-learned potential would check that it is not an artifact of a single force field.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. This paper uses reactive molecular dynamics (ReaxFF) simulations to study how the degree of crystallinity of MoS2 affects its frictional and wear behavior. The authors construct crystalline, polycrystalline (PC75, PC25, PC10), and amorphous MoS2 substrates, then slide a rigid spherical carbon body over each surface. The coefficient of friction (COF) is extracted from the slope of friction force versus normal load over five loads, with three independent runs per structure. The reported COF increases monotonically as crystallinity decreases: 0.00014 for crystalline, 0.05 for PC75, 0.11 for PC25, 0.17 for PC10, and 0.28 for amorphous MoS2. The authors attribute the ultralow friction of crystalline MoS2 to a flat sliding surface and frictional anisotropy that lets the rigid body follow a zig-zag, minimum-energy path. For wear, a single high-load (24.85 nN) sliding run per structure is performed, and wear depth is measured from the initial and final vertical position of the rigid body. Crystalline MoS2 shows a depth of 5.32 Å versus 14.8 Å for amorphous MoS2, leading the authors to conclude that higher crystallinity enhances wear resistance. The paper frames the work as extending MoS2 superlubricity studies to rigid-body sliding conditions and as evidence that structural order, not just surface chemistry, controls tribological performance.

Significance. The friction half of this study is a useful contribution: the COF trend with crystallinity is large, monotonic, and based on repeated runs with reported error bars, and the crystallinity metric is clearly defined in Appendix A. If the wear claim were equally well supported, the paper would offer a clean computational demonstration that structural order controls both friction and wear in MoS2. However, the wear conclusion rests on a single simulation per structure and a penetration-depth proxy rather than a material-loss measure, and it directly contradicts the authors' own cited experiment [29]. The paper also uses a ReaxFF potential that is validated only on energetics, with no direct validation for sliding contact or wear. These issues do not undermine the COF trend, but they prevent the paper's combined friction-and-wear message from being accepted as established.

major comments (2)
  1. [Sec. 3.3] The wear-resistance claim is not supported by the evidence as presented. The 'depth of wear' is measured from the initial and final position of the rigid body, a metric that conflates elastic and plastic substrate deformation, rigid-body sinking, and actual material removal. Only one sliding run at 24.85 nN is performed per structure, so no statistical uncertainty is available. The reported amorphous wear depth of 14.8 Å exceeds the rigid-body radius (10 Å), indicating that the quantity measured is largely indentation rather than wear. Moreover, the paper's own introduction (Section 1) cites Hesam et al. [29], who found amorphous MoS2 to be approximately four times more wear resistant than its crystalline counterpart; the opposite ordering in Figure 6 is never reconciled with that experiment. The authors should either replace the penetration-depth metric with a direct measure of material loss (for example, counting atoms detached from the substrate), repeat the wear simulations to provide error bars, and explicitly discuss the discrepancy with [29], or reframe the conclusion as a statement about deformation resistance under the specific rigid-body, high-load conditions simulated here.
  2. [Sec. 2] The ReaxFF Mo-S-C parameters from Ponomarev and Polcar [62] are described as validated to reproduce the energetics of crystalline and amorphous MoS2, but the manuscript provides no validation for sliding contact, adhesion, or wear. Because the quantitative claims (COF = 0.00014, wear depths of 5.32 Å versus 14.8 Å) and the qualitative ordering depend on the accuracy of interlayer shear and carbon-MoS2 interactions, the transferability of the force field to tribological loading is load-bearing. I request that the authors either supply additional validation (for example, comparison with experimental friction coefficients or with independent interlayer interaction data) or explicitly state this limitation and discuss how it might affect the magnitude, and possibly the ordering, of the simulated friction and wear results. This is a request for transparency and evidence, not an assertion that the force field is wrong.
minor comments (5)
  1. [Abstract] The sentence 'Our results reveal a pronounced reduction in COF with decreasing crystallinity' should read 'increase in COF with decreasing crystallinity' (or 'reduction with increasing crystallinity') to be consistent with Table 1.
  2. [Sec. 3.3 and Table 1] The amorphous COF is given as 0.2736 in the text and 0.28 in the abstract and Table 1; please use consistent significant figures.
  3. [Figure 5 caption] The caption states 'red and blue colored regions indicating atoms at the bottom at top surface'; this is unclear and should read 'highest and lowest positions among the selected atoms.'
  4. [Sec. 3.3] The phrase 'complete structural disintegration for both for PC10 and amorphous' contains a typo ('both for').
  5. [Appendix A] The crystallinity criteria are described as minimizing false positives, but no sensitivity analysis is given; a brief test of how the reported percentages change with the chosen angle thresholds would strengthen the metric.

Circularity Check

0 steps flagged · score 1.0 of 10

No circular derivation: the COF and wear trends are direct simulation outputs, and the only self-citations provide force-field provenance or supportive context rather than forcing the reported trends.

full rationale

The paper's central derivation chain is: construct MoS2 structures with crystallinity percentages defined by an independent structural criterion (Appendix A), perform ReaxFF MD sliding simulations (Section 2), record normal and lateral forces, define COF as the slope of the force-load regression (Section 3.1), and define wear depth from the rigid body's initial and final vertical position (Section 3.3). None of these steps fits the crystallinity-friction or crystallinity-wear trend into the data. The crystallinity labels are structural (Mo-Mo-Mo angle criteria), the COF values are averaged simulation outputs with error bars, and the wear depths are direct measurements from high-load runs. The ReaxFF parameters are cited to the authors' prior force-field papers [62,63], but that citation is provenance for the interaction model, not a derivation of the friction or wear ordering; the cited work validates energetics, not the tribological outputs claimed here, so the central claims do not reduce to the citation. The comparison to experimental work, including the contradictory wear result of Hesam et al. [29], is a validation and correctness question rather than a circularity: no equation in the paper defines crystallinity in terms of friction or wear, and no fitted parameter is renamed as a prediction. The self-citations (e.g., [46,52,65]) are supportive or contextual, and the simulation results stand independently of them. I therefore find no significant circularity, while noting that the force field and single-run wear metric are legitimate reproducibility and validity concerns.

Assumptions & free parameters 2 free parameters · 4 assumptions · 0 invented entities

The central results depend on an empirical force field from the same group and on hand-chosen structural criteria. No new entities are introduced. The ledger is small, but the ReaxFF parameters carry a lot of weight, and the crystallinity thresholds define the independent variable without a sensitivity test.

free parameters (2)
  • ReaxFF Mo-S-C force field parameters (Ponomarev and Polcar, 2024) = See Ref. [62]
    All COF and wear values are outputs of this empirical force field, which was fit to energetics in prior work, not to friction data. The central results inherit its accuracy.
  • Crystallinity classification thresholds = Angle windows 50-70 deg, 110-130 deg, and over 170 deg; minimum counts 5, 5, 2
    Appendix A defines crystalline Mo by these hand-chosen thresholds. The degree of crystallinity, the paper's main independent variable, depends on this choice. The authors claim it minimizes false positives, but no sensitivity test is shown.
assumptions (4)
  • domain assumption ReaxFF accurately describes Mo-S-C interactions in both crystalline and amorphous MoS2 under sliding contact.
    Invoked in Section 2 when selecting ReaxFF over REBO; the force field is validated only against energetics and crystallization, not against tribological quantities.
  • domain assumption The rigid, non-deformable carbon sphere is a valid counterbody for studying friction and wear of MoS2.
    Section 2: the spherical rigid body is fixed via 'fix rigid'; this suppresses deformation and wear of the counterbody, which may bias wear mechanisms.
  • ad hoc to paper Mo-Mo-Mo angular criteria in Appendix A faithfully quantify crystallinity independent of the friction outcome.
    The thresholds are introduced ad hoc in Appendix A to define the independent variable; no external calibration is provided.
  • domain assumption Finite simulation size (63 x 40 x 35.5 A^3, periodic X/Y) and 1 ns sliding time are sufficient to reach steady-state friction.
    Section 2 describes the cell and 1 ns runs; no convergence checks versus cell size or sliding length are reported.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Effect of crystallinity on the frictional and wear performance of molybdenum disulfide: A molecular dynamics study." pith.science (2026). https://pith.science/paper/5ECEBTUF

@misc{pith2026250600690,
  author       = {Pith},
  title        = {Pith review of: Effect of crystallinity on the frictional and wear performance of molybdenum disulfide: A molecular dynamics study},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5ECEBTUF}},
  note         = {Machine review of arXiv:2506.00690}
}
read the original abstract

The frictional and wear performance of molybdenum disulfide (MoS2) is significantly influenced by its intrinsic arrangement of crystals or crystallinity. In this study, we investigate the effect of crystallinty on coefficient of friction (COF) and wear in MoS2 using a suite of reactive molecular dynamics (MD) simulations. A range of configurations, from amorphous to crystalline, is modeled to capture the effect of structural order on the tribological behavior. To study friction and wear, we simulate the sliding of a spherical rigid carbon body over the MoS2 surface under varying crystallinity conditions. Our results reveal a pronounced reduction in COF with decreasing crystallinity, with crystalline MoS2 exhibiting superlubricity. This behavior is attributed to the preservation of a flat sliding surface and frictional anisotropy, which enables lateral movement along low-resistance paths. In contrast, amorphous and polycrystalline MoS2 with lower degrees of crystallinity displays a substantially higher COF, driven by increased surface roughness and atomic-scale energy dissipation. Furthermore, we examine the wear mechanisms under high normal loads, demonstrating that crystallinity enhances wear resistance by mitigating material deformation. These findings provide atomic-scale insights into the tribological performance of MoS2, emphasizing the critical role of structural order in achieving ultralow friction. Our work corroborates with previous studies on superlubricity in MoS2 and extends this understanding to rigid-body sliding conditions, offering valuable implications for designing low-friction and wear resistant solid lubricants.

Figures

Figures reproduced from arXiv: 2506.00690 by the authors.

Figure 1
Figure 1. Sequence of annealing and quenching temperatures used to create polycrystalline models with different percentages of crystals achieve the desired polycrystalline structures, a sequence of annealing and quenching processes were employed, allowing control over the percentage of crystallinity. Detailed information regarding the annealing temperatures and the resulting crystallinity percentages are provided in [PITH_FU… view at source ↗
Figure 2
Figure 2. Atomic conJiguration of PC25, rigid body and [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. Plots of normal force with respect to frictional force for MoS [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Zig-zag movement of rigid body during sliding simulation in crystalline MoS2 and lowest positions among the group of selected atoms. A constant color in the plot implies that all the atoms are positioned at the same level. This is evident in Figures 5a & 5b (crystallin…
Figure 5
Figure 5. Figure 5: Geometric surface meshing for MoS2 with different levels of crystallinity. Color gradients are used to represent the surface roughness with red and blue colored regions indicating atoms at the bottom at top surface. Crystalline, PC75, and PC25 con=igurations (Figures 7…
Figure 6
Figure 6. Figure 6: Atomic representation of wear in MoS2 with different levels of crystallinity [PITH_FULL_IMAGE:figures/full_fig_p010_6.png]
Figure 7
Figure 7. Figure 7: Color coded displacement of atoms at the top surface of MoS [PITH_FULL_IMAGE:figures/full_fig_p011_7.png]

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

8 extracted references · 8 canonical work pages

  1. [2]

    4 We modeled MoS2 in various con=igurations, including crystalline, polycrystalline (PC), and amorphous structures

    Modeling and simulation In this section, we describe the modeling strategy used to create MoS2 with varying degrees of crystallinity, along with the details of the MD simulations. 4 We modeled MoS2 in various con=igurations, including crystalline, polycrystalline (PC), and amorphous structures. The PC models were created with 75%, 25%, and 10% crystallini...

    work page

  2. [3]

    First, we compute the COF for different degrees of crystallinity, ranging from amorphous to crystalline structures

    Results and discussion In this section, we present and discuss the key =indings from our study, focusing on the in=luence of crystallinity on the frictional and wear behavior of MoS2. First, we compute the COF for different degrees of crystallinity, ranging from amorphous to crystalline structures. Second, we investigate the underlying friction mechanisms...

    work page

  3. [4]

    crystalline Mo

    Conclusion This study systematically investigates the effect of crystallinity on the frictional and wear performance of MoS2 using rigid body sliding simulations. Our =indings reveal that crystallinity plays a crucial role in determining the COF and wear resistance. Crystalline MoS2 exhibits superlubricity, with an exceptionally low COF of 0.00014, due to...

    work page 2017

  4. [8]

    Renevier, J

    N. Renevier, J. Hamphire, V. Fox, J. Witts, T . Allen, D. Teer, Advantages of using self-lubricating, hard, wear-resistant MoS2-based coatings, Surface and Coatings Technology 142 (2001) 67–77. [9] M. Buscema, M. Barkelid, V . Zwiller , H. S. van der Zant, G. A. Steele, A. Castellanos-Gomez, Large and tunable photothermoelectric effect in single-layer MoS...

    work page 2001

  5. [22]

    Wa n g , L

    P. Wa n g , L . Q i a o , J . X u , W. L i , W. L i u , E r o s i o n m e c h a n i s m o f M o S2-based =ilms exposed to atomic oxygen environments, ACS Applied Materials & Interfaces 7 (23) (2015) 12943–12950. [23] T. Ku b a r t , T. Po l c a r, L . Ko p e c k y ` , R . N o va k , D. N o v ´ a k o v ´ a , Te m p e r a t u r e d e p e n d e n c e o f t r...

    work page 2015

  6. [36]

    H. P . Lal, B. R. Abhiram, D. Ghosh, Prediction of nonlocal elasticity parameters using high-throughput molecular dynamics simulations and machine learning, European Journal of Mechanics-A/Solids 103 (2024) 105175. [37] B. R. Abhiram, D. Ghosh, In=luence of nano=iller agglomeration on mechanical properties of nanocomposites: a multiscale study, in: AIAA S...

    work page 2024

  7. [49]

    S h i , Z

    Y. S h i , Z . C a i , J . P u , L . Wa n g , Q. Xu e , I n t e r fa c i a l m o l e c u l a r d e fo r m a t i o n m e c h a n i s m fo r l o w f r i c t i o n o f MoS2 determined using ReaxFF-MD simulation, Ceramics International 45 (2) (2019) 2258–2265. [50] A. Bondarev, I. Ponomarev, R. Muydinov, T . Polcar , Friend or foe? revising the role of oxygen...

    work page 2019

  8. [61]

    H. M. Aktulga, J. C. Fogarty, S. A. Pandit, A. Y . Grama, Parallel reactive molecular dynamics: Numerical methods and algorithmic techniques, parallel computing 38 (4-5) (2012) 245–259. [62] I. Ponomarev, T . Polcar , A new reactive Mo-SC force =ield to explore structure and properties of MoS2-C 2D materials and =ilms, Computational Materials Science 241 ...

    work page 2012

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.