REVIEW 4 major objections 5 minor 25 references
Visualizing phonon edge states on molybdenum disulphide
T0 review · 4 major / 5 minor · reviewed 2026-08-15 · deepseek-v4-flash
Pith's one-line read Low-energy phonon modes pinned to the edges of MoS2 nanoparticles explain the enhanced atomic vibrations and the contrast fade seen in high-resolution electron microscopy near the boundary.
desk verdict A plausible and well-communicated edge-phonon explanation for enhanced vibrations at MoS2 edges, with a genuinely useful active-learning ENNP workflow, but the central mechanism rests on an unvalidated ML potential for low-frequency phonons. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the edge-localized phonon mode, identified by projecting the phonon density of states of a finite strip onto the formula units at the edges and coloring band-structure modes by the squared vibration amplitude weight on the edge atoms. The argument is carried by the frozen-phonon imaging model: atomic configurations sampled from molecular dynamics are used to generate and average exit waves, and the same machine-learned potential, an equivariant graph neural network trained on roughly 400 DFT configurations, is used for both the molecular dynamics and the phonon calculations. The mechanism is that low-energy edge modes, with energies below roughly $k_B T$, increase the mean-square vibrational amplitude near the edge, which smears the projected electrostatic potential and lowers the exit-wave peaks.
What would settle it
Compute the phonon band structure and edge-projected density of states of the same 20-unit-cell MoS2 strip using direct DFT force constants and compare with the machine-learned potential; if the low-frequency edge modes are not reproduced, the explanation fails. Also, measure the temperature dependence of the HR-TEM contrast decay and check that it follows the occupation of the predicted edge-mode energies.
Extended reading notes
Core claim
On the paper's own terms: the enhanced and correlated vibrational amplitudes near the edges of carbon-supported MoS2 nanoparticles, observed in molecular dynamics and in a recent HR-TEM experiment, are due to low-energy phonon modes localized at the nanoparticle edges. Projecting the phonon density of states onto edge formula units shows a strong shift toward low frequencies relative to the bulk sheet, and the band structure of a 20-unit-cell strip contains bands whose weight is concentrated on the Mo or S edge, some extending a few lattice constants into the interior. These modes suppress the imaginary part of the simulated electron exit wave at atomic columns up to a nanometer from the edge, matching the experimental contrast decay. The paper further claims that a workflow combining active learning with equivariant neural network potentials can deliver an accurate potential for this supported-nanoparticle system from only slightly more than 400 DFT configurations, with an ensemble of five potentials providing a posteriori error estimates.
Load-bearing premise
The machine-learned potential, trained on about 400 DFT configurations, reproduces the true low-energy vibrational modes of MoS2 edges even though its phonon accuracy is not separately validated against DFT.
Editorial extensions
If this is right
- HR-TEM image interpretation for MoS2 nanoparticles must include correlated atomic motion: the independent-atom Einstein model cannot reproduce the observed contrast decay.
- The vibrational amplitude enhancement extends roughly a nanometer from the edge, so quantitative image simulation must include atoms that are not under-coordinated.
- Edge phonon modes, not just the low coordination of the outermost atoms, are the cause of the enhanced amplitudes.
- The active-learning and equivariant-neural-network workflow can produce usable potentials from about 400 DFT configurations and can flag out-of-distribution configurations through ensemble disagreement.
Reading between the lines
- If edge-localized low-energy phonons are generic in finite 2D crystals, HR-TEM and STEM contrast simulations of other 2D materials will need the same correlated-molecular-dynamics treatment.
- The ensemble disagreement among the five potentials could be turned into a per-atom uncertainty map during molecular dynamics, giving a practical way to detect when the phonon picture itself is unreliable.
- A testable extension is temperature dependence: the edge modes' occupation grows with temperature, so the edge-to-bulk contrast decay should steepen at higher temperatures in a way a static disorder model would not predict.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents a combined molecular dynamics (MD) and phonon study of carbon-supported MoS2 nanoparticles, using an equivariant neural network potential (ENNP) trained on a limited set of DFT (PBE-D3) calculations via an active-learning workflow. The authors report enhanced vibrational amplitudes near the nanoparticle edges, both for Mo and S atoms, and show that this enhancement extends roughly a nanometer from the edge. They attribute this behavior to low-energy phonon modes localized at the edges, identified by projecting the phonon density of states (DOS) onto the edge formula units of a 20-unit-cell-wide strip. They further simulate HR-TEM exit waves using the frozen-phonon method with MD snapshots and observe a decay of the imaginary-part peak values near the edges, consistent with the experimental findings of Chen et al. The paper also emphasizes the value of correlated MD over independent-atom models for quantitative image interpretation and presents the active-learning workflow as a way to obtain accurate machine-learning potentials with built-in error estimation.
Significance. If the central claim holds, the paper makes a useful contribution by explaining the experimentally observed edge-induced contrast decay through phonon edge states, and by demonstrating the necessity of correlated atomic motion in HR-TEM simulations. The active-learning framework for training an ENNP with a few hundred DFT configurations is a practical step toward simulating realistic nanostructures beyond DFT. The manuscript is clearly written and the workflow is described in sufficient detail to be reproduced. However, the central explanation rests on the reliability of the ENNP for low-frequency phonons, a point that is not quantitatively validated. The paper therefore has clear potential but requires strengthening of the evidence before its main claims can be considered fully supported.
major comments (4)
- [Generation of the machine learning potentials and Phonon band structure simulations] The central claim that low-energy phonon edge states cause the enhanced vibrational amplitudes is not supported by any quantitative validation of the ENNP for phonon properties. The manuscript reports neither force/energy test-set errors nor a comparison of the ENNP phonon DOS, band structure, or force constants against DFT results. Because the same potential is used for the MD simulations, the exit-wave simulations, and the phonon analysis, any inaccuracy in the low-frequency force constants propagates directly into the edge-state conclusion. Please provide test-set error metrics and a direct DFT comparison of the low-energy phonon modes (for example, a DFT phonon calculation on a small strip or cluster) to demonstrate that the edge states are not an artifact of the machine-learning potential.
- [Phonon band structure simulations (Appendix)] The phonon calculations are performed on a free-standing 2D sheet, which possesses a quadratic transverse acoustic mode (the bending mode) that the authors acknowledge may distort the lowest parts of the DOS. This is a serious concern because the edge-induced DOS increase appears precisely at low energies. Without a correction, a constraining of the out-of-plane motion, or a comparison to a supported strip, it remains possible that the low-energy enhancement is dominated by the free-standing bending mode rather than by localized edge states. The authors should quantify the effect of this artifact on the edge-projected DOS, for instance by repeating the calculation with the out-of-plane directions constrained or by applying a weak restoring force.
- [Molecular Dynamics and Nanoparticle simulations] The vibrational amplitudes underlying the comparison with experiment are computed from a single 5 ps production trajectory per condition, following 5 ps of thermalization. This sampling may be insufficient to converge the variance of atomic positions, especially for the low-frequency edge modes whose periods are on the order of several picoseconds. The manuscript does not report block averaging, independent trajectories, or any convergence analysis with respect to simulation length. Given that the exit-wave comparison and the central interpretation rest on these amplitudes, the quantitative conclusions are not yet statistically supported. Additional analysis of statistical convergence is needed.
- [Abstract and Nanoparticle simulations] The claim of a 'built-in error estimation' is stronger than what is actually demonstrated. The ensemble of five potentials is used during active learning and for validation of shorter sequences, but the production MD and the phonon calculations use a single potential trained on the full generation-5 dataset. No uncertainty estimates or error bars are provided for the central results (vibrational amplitude profiles, exit-wave peak values, or phonon DOS). Please clarify the role of the ensemble in the final results and provide quantitative uncertainty measures for the main figures.
minor comments (5)
- [Abstract] The parenthetical citation 'by Chen et al (Nature Communications 12, 5007 (2021).' is missing a closing parenthesis.
- [Exit wave simulations] The text refers to 'the 3 nm wide MoS2 nanoparticle in Fig. 1(b)', but Fig. 1 shows a nanoparticle with radius 3 nm; the terminology is ambiguous and should be consistent.
- [Molecular Dynamics] The description of the trajectory collection could state explicitly that only a single trajectory per condition was used and that no block or ensemble averaging was performed, which would make the sampling limitations transparent.
- [Figure 3] The caption of Figure 3 mentions that bands are 'colored by their weight' but does not define the weight; the definition is given in the text, but adding it to the caption would improve readability.
- [Conclusion] The conclusion states that the potential 'shows that the vibrational amplitudes are significantly enhanced in a region near the edge' without acknowledging the potential limitations regarding trajectory length and phono validation; a brief caveat would be appropriate.
Circularity Check
No significant circularity: the edge-state explanation emerges from a DFT-trained ML potential and is not fitted to the target observation.
full rationale
The paper's derivation chain is: PBE-D3 DFT energies and forces feed an active-learning equivariant neural network potential (ENNP); the ENNP drives MD simulations from which vibrational amplitudes and frozen-phonon exit waves are computed; the same ENNP is then used for phonon band-structure and DOS calculations on a finite strip, revealing low-energy modes localized near the edges. None of these links is fitted to the experimental observation of Chen et al.; that experiment is used only as a qualitative comparison after the simulations are already produced. The ENNP is trained exclusively on DFT data, not on phonon spectra, edge-state weights, or the experimental exit-wave decay, so the enhanced edge vibrational amplitude and the edge-localized phonon modes are emergent properties of the fitted interatomic forces rather than quantities encoded by construction. Using the same potential for MD and phonon analysis is internally consistent, not circular, because the potential was not constructed to reproduce the target phenomenon. The only self-citations (refs. 15 and 16) support general methodological claims about active learning and the data efficiency of equivariant neural networks; they are not load-bearing for the edge-state mechanism and no uniqueness theorem is imported from the authors' prior work. The paper's main weakness is the absence of a quantitative phonon validation of the ENNP, such as force RMSEs or a DFT comparison of low-frequency modes; that is a correctness and robustness risk, not a circularity. No equation or fitted parameter is renamed as a prediction, and no reduction of the central claim to its inputs can be exhibited. Accordingly, no circular step is identified and the circularity score is 0.
Assumptions & free parameters
free parameters (4)
- ENNP model parameters (weights and biases of NequIP networks) =
trained on ~250-300 DFT configurations (generation 5); 5 ensemble members
- Active learning hyperparameters: force disagreement threshold and the 100 eV/A filter cutoff =
100 eV/A; configurations with larger forces discarded
- MD collection window =
5 ps after 5 ps thermalization
- Edge-cutoff difference for S-rich vs Mo-rich nanoparticles =
1 A smaller cutoff for Mo than S
assumptions (4)
- domain assumption PBE-D3 DFT is an accurate reference for the energetics and forces of MoS2/graphite interfaces and edges.
- domain assumption The NequIP equivariant neural network can extrapolate to the nanoparticle sizes and temperatures used in the MD after active learning on smaller clusters.
- domain assumption Phonon modes of a free-standing 20-unit-cell MoS2 strip are a valid proxy for the edge modes of the circular supported nanoparticles and of the experimental particles.
- domain assumption The frozen phonon approximation with 100 snapshots converges the simulated exit wave.
Cite this review
Pith. "Pith review of Visualizing phonon edge states on molybdenum disulphide." pith.science (2026). https://pith.science/paper/F33U3VCH
@misc{pith2026250908497,
author = {Pith},
title = {Pith review of: Visualizing phonon edge states on molybdenum disulphide},
year = {2026},
howpublished = {\url{https://pith.science/paper/F33U3VCH}},
note = {Machine review of arXiv:2509.08497}
}
read the original abstract
We employ Molecular Dynamics (MD) simulations to study atom vibrational amplitudes in carbon-supported Molybdenum Disulphide (MoS2) nanoparticles. Enhanced and correlated atom vibrational amplitudes are observed as the nanoparticle edges are approached from the bulk, consistent with recent experimental High-Resolution Transmission Electron Microscopy (HR-TEM) observations by Chen et al (Nature Communications 12, 5007 (2021). Analysis of phonon modes in finite systems explains the experimental observation by low-energy phonon modes confined at the nanoparticle edge, underscoring the need of full MD modeling for accurate HR-TEM image interpretation. Noticeably, we introduce a workflow for training Equivariant Neural Network-based machine learning potentials using limited Density Functional Theory (DFT) calculations. This approach effectively captures both covalent and van der Waals interactions, enabling accurate extrapolations of DFT calculations to larger systems with built-in error estimation.
Figures
Figures from the paper (5 more)
Reference graph
Works this paper leans on
-
[1]
C. Kisielowskiet al., Detection of single atoms and buried defects in three dimensions by aberration-corrected elec- tron microscope with 0.5- ˚A information limit, Microsc. Microanal.14, 469 (2008)
work page 2008
-
[2]
R. Erni, M. D. Rossell, C. Kisielowski, and U. Dah- men, Atomic-resolution imaging with a sub-50-pm elec- tron probe, Phys. Rev. Lett.102, 096101 (2009)
work page 2009
-
[3]
F.-R. Chen, D. V. Dyck, C. Kisielowski, L. P. Hansen, B. Barton, and S. Helveg, Probing atom dynamics of ex- cited Co-Mo-S nanocrystals in 3D, Nature Comm.12, 5007 (2021)
work page 2021
-
[4]
Z. Chen, Y. Jiang, Y.-T. Shao, M. E. Holtz, M. Odstrˇ cil, M. Guizar-Sicairos, I. Hanke, S. Ganschow, D. G. Schlom, and D. A. Muller, Electron ptychography achieves atomic-resolution limits set by lattice vibrations, Science372, 826 (2021)
work page 2021
-
[5]
R. F. Loane, P. Xu, and J. Silcox, Thermal vibrations in convergent-beam electron diffraction, Acta Crystallogr. A47, 267 (1991)
work page 1991
-
[6]
S. Hillyard and J. Silcox, Detector geometry, thermal dif- fuse scattering and strain effects in ADF STEM imaging, Ultramicroscopy58, 6 (1995)
work page 1995
-
[7]
D. Van Dyck, Persistent misconceptions about incoher- ence in electron microscopy, Ultramicroscopy111, 894 (2011)
work page 2011
-
[8]
J. S. Smith, B. Nebgen, N. Lubbers, O. Isayev, and A. E. Roitberg, Less is more: Sampling chemical space with active learning, J. Chem. Phys.148, 241733 (2018)
work page 2018
Show all 25 references
-
[9]
Thomas, T
N. Thomas, T. Smidt, S. Kearnes, L. Yang, L. Li, K. Kohlhoff, and P. Riley, Tensor field networks: Rotation- and translation-equivariant neural networks for 3D point clouds, arXiv.org , 1802.08219 (2018)
2018 arXiv
-
[10]
Batzner, A
S. Batzner, A. Musaelian, L. Sun, M. Geiger, J. P. Mailoa, M. Kornbluth, N. Molinari, T. E. Smidt, and B. Kozinsky, E(3)-equivariant graph neural networks for data-efficient and accurate interatomic potentials, Nature Comm.13, 2453 (2022)
2022
-
[11]
Madsen and T
J. Madsen and T. Susi, The abTEM code: transmission electron microscopy from first principles, Open Research Europe1, 24 (2021)
2021
-
[12]
See Supplemental Material at [URL will be inserted by publisher] for animated renderings of the phonon modes in Fig. 4. The files are named by the edge, the k-point and the energy in meV
-
[13]
K. T. Sch¨ utt, P.-J. Kindermans, H. E. Sauceda, S. Chmiela, A. Tkatchenko, and K.-R. M¨ uller, SchNet: A continuous-filter convolutional neural network for model- ing quantum interactions, Adv. Neural Inf. Process Syst. 30, 991 (2017)
2017
-
[14]
J. S. Smith, O. Isayev, and A. E. Roitberg, ANI-1: an ex- tensible neural network potential with DFT accuracy at force field computational cost, Chem. Sci.8, 3192 (2017)
2017
-
[15]
A. E. Mikkelsen, J. Schiøtz, T. Vegge, and K. W. Ja- cobsen, Is the water/Pt(111) interface ordered at room temperature?, J. Chem. Phys.155, 224701 (2021)
2021
-
[16]
Nu˜ nez Valencia, W
C. Nu˜ nez Valencia, W. B. Lomholdt, M. H. Leth Larsen, T. W. Hansen, and J. Schiøtz, Beam induced heating in electron microscopy modeled with machine learning interatomic potentials, Nanoscale16, 5750 (2024)
2024
-
[17]
Enkovaaraet al., Electronic structure calculations with GPAW: A real-space implementation of the projector augmented-wave method (2010)
J. Enkovaaraet al., Electronic structure calculations with GPAW: A real-space implementation of the projector augmented-wave method (2010)
2010
-
[18]
J. J. Mortensenet al., GPAW: An open Python package for electronic structure calculations, J. Chem. Phys.160, 092503 (2024)
2024
-
[19]
J. P. Perdew, K. Burke, and M. Ernzerhof, Generalized Gradient Approximation Made Simple, Phys. Rev. Lett. 77, 3865 (1996)
1996
-
[20]
Grimme, J
S. Grimme, J. Antony, S. Ehrlich, and H. Krieg, A con- sistent and accurate ab initio parametrization of density 6 functional dispersion correction (DFT-D) for the 94 ele- ments H-Pu, J. Chem. Phys.132, 154104 (2010)
2010
-
[21]
Hjorth Larsenet al., The atomic simulation environ- ment - A Python library for working with atoms, J
A. Hjorth Larsenet al., The atomic simulation environ- ment - A Python library for working with atoms, J. Phys.: Condens. Matter29, 273002 (2017)
2017
-
[22]
B. Deng, P. Zhong, K. J. Jun, J. Riebesell, K. Han, C. J. Bartel, and G. Ceder, CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling, Nat. Mach. Intell.5, 1031 (2023)
2023
-
[23]
Batatiaet al., A foundation model for atomistic mate- rials chemistry, arXiv.org , 2401.00096 (2024)
I. Batatiaet al., A foundation model for atomistic mate- rials chemistry, arXiv.org , 2401.00096 (2024)
2024 arXiv
-
[24]
Lobato and D
I. Lobato and D. Van Dyck, An accurate parameteriza- tion for scattering factors, electron densities and electro- static potentials for neutral atoms that obey all physical constraints, Acta Crystallogr. A70, 636 (2014). Appendix: computational details Generation of structures...
2014
-
[25]
The Kohn-Sham wave- functions were described in a plane wave basis with a cutoff of 600 eV
added to properly describe the van der Waals in- teractions between the layers. The Kohn-Sham wave- functions were described in a plane wave basis with a cutoff of 600 eV. As the supercell is large in the periodic directions, the Brillouin zone is sampled with a single k-point...
Reviewed August 15, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.