Pith. sign in

REVIEW 4 major objections 7 minor 42 references

Ultrafast atomic dimerization of Peierls distortion in semimetal molybdenum ditelluride

T0 review · 4 major / 7 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 Under an 800 nm laser pulse, few-layer 1T'-MoTe2 shears and shortens its Mo-Mo bonds within picoseconds, then holds a distorted structure for nanoseconds.

desk verdict A genuinely new ultrafast observation in few-layer 1T'-MoTe2, but the dimerization claim is supported by one unmodeled spacing and a thermal calibration that is partly circular. read the letter →

arxiv 2505.21096 v1 pith:N7JYKSWH submitted 2025-05-27 cond-mat.mtrl-sci

Zhong Wang , Chunlong Hu , Changchang Gong , Fuyong Hua , Qian You , Wenxi Liang This is my paper
classification cond-mat.mtrl-sci PACS 87.15.ht61.43.-j78.47.J
keywords 1T'-MoTe2ultrafastelectrondiffractionPeierlsdistortiondimerizationsheardisplacementDebye-Wallereffectpairdistributionfunctiondensityfunctionaltheory
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

The paper tries to establish that illuminating few-layer 1T'-MoTe2 with an 800 nm laser pulse does not just heat the crystal: within a few picoseconds the layers shear and the Mo-Mo bonds shorten, and under strong excitation a distorted structure persists into nanoseconds. It reaches this picture by combining ultrafast electron diffraction with density functional theory (DFT), reading Bragg-peak intensities, interplanar spacings, and pair distribution functions. The authors argue that both motions are driven by photodoped electrons occupying antibonding states near the Fermi level, and that the shear is a partial step toward the Td phase. If correct, this shows light can drive the Peierls distortion in a semimetal rather than only melt it, and it gives a time-resolved view of the first steps of a 1T'-to-Td transition.

What carries the argument

The argument runs through the kinematic structure factor $F(hkl)=\sum_j f_j \exp[-i2\pi(hx_j+ky_j+lz_j)]$, which converts changes in atomic fractional coordinates into changes in Bragg intensities. The authors compare measured intensity drops with structure-factor predictions from five distortion models and use the same models to compute interplanar-spacing changes; only shear plus Mo-Mo bond shortening matches the picosecond data. The Debye-Waller contribution from random atomic vibrations is calibrated with the (020) peak, and the pair distribution function of reduced diffraction elements is used to follow atomic-pair separations. On the electronic side, the crystal orbital Hamilton population (COHP), a per-bond decomposition of band-structure energy, carries the bonding/antibonding argument: the near-Fermi Mo d states are antibonding, so photoexcited electrons change the bond order and drive dimerization.

What would settle it

Measure the (020) intensity drop at 3 ps while independently determining the lattice temperature from high-angle diffuse scattering or a phonon-mode thermometer; if the (020) drop exceeds the Debye-Waller prediction for that temperature, the thermal calibration is contaminated. A second check is fluence scaling at low pump power, where the purely thermal (020) drop should be linear in fluence; any nonlinear early-time deviation signals a nonthermal (020) response.

Watch

Extended reading notes

Core claim

The central claim is that the measured lattice response in few-layer 1T'-MoTe2 after 800 nm excitation is a superposition of three motions rather than a single thermal expansion. At 3 ps, the interplanar spacing data match a shear model in which the monoclinic angle $\beta$ is reduced by roughly 0.5°, while the (020) spacing contracts rapidly; the authors attribute the contraction to Mo-Mo bond shortening, that is, enhanced dimerization. Structure-factor calculations show that no single distortion model reproduces all Bragg intensities at 3 ps, but the combination of shear and Mo-Mo bond shortening is the one consistent with the best-measured peak. At 3 ns, all planes expand nearly uniformly while some intensities keep changing, which the authors interpret as a metastable structure with additional b-c plane distortion. DFT and crystal orbital Hamilton population (COHP) calculations supply the mechanism: photodoped electrons occupy interlayer antibonding Te p states, driving the shear, and antibonding Mo d states near the Fermi level, so the relaxed carrier population lowers the COHP integral and shortens the Mo-Mo bond.

Load-bearing premise

The load-bearing assumption is that the (020) Bragg peak's intensity drop is caused entirely by random atomic vibrations (the Debye-Waller effect), so it can be used to subtract thermal motion from every other peak; yet at 3 ns the paper itself invokes an additional nonthermal (020) drop from b-c plane distortions, so any early-time nonthermal (020) response would corrupt the subtraction.

Editorial extensions

If this is right

  • Photoexcitation can drive the Mo-Mo dimerization of the Peierls distortion instead of melting it; the dimer bond shortens within a few picoseconds.
  • The roughly 0.5° reduction of $\beta$ at 3 ps is a partial step along the 1T'-to-Td transition; completing the transition would require more electron doping or fewer compensating holes.
  • The nanosecond-lived metastable structure, seen only above about 1 mJ/cm2 fluence, implies a long-lived photoinduced lattice state distinct from ordinary thermal expansion.
  • Using one Bragg peak as an internal Debye-Waller thermometer can isolate nonthermal lattice motions in layered materials, provided that peak is truly unaffected by structural distortions.

Reading between the lines

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

  • The authors themselves note that the observations give no direct evidence of the distortions; an experiment with higher PDF resolution that directly resolves the nearest Mo-Mo distance contracting would confirm the dimerization claim without relying on structure-factor modeling.
  • If the mechanism is electronic, electrostatic gating or chemical intercalation that dopes electrons into the same antibonding states should reproduce the shear-then-dimerize sequence at lower power, a testable prediction the paper does not make.
  • The nanosecond metastable state, if it is a distinct structural phase rather than slow thermal relaxation, could act as a bistable lattice switch; checking whether it reverts to 1T' or anneals into Td would map the free-energy landscape.
  • The residual mismatch among Bragg peaks at 3 ps hints at defect or polaron contributions; repeating the measurement on a thicker or higher-quality crystal would show whether that mismatch shrinks.
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

4 major / 7 minor

Summary. The paper reports ultrafast electron diffraction (UED) measurements on few-layer 1T'-MoTe2 following 800 nm photoexcitation, together with DFT, DOS, and COHP calculations. The authors propose that the observed lattice response within a few picoseconds is best described by a combination of interlayer shear displacement (a decrease of the monoclinic angle β by about 0.5°) and intralayer Mo-Mo bond shortening (dimerization), followed by a long-lived metastable structure on the nanosecond timescale. Evidence is drawn from Bragg peak intensity changes, interplanar spacing changes, a qualitative pair distribution function (PDF) analysis, and electronic-structure calculations. The paper explicitly acknowledges that none of the five distortion models reproduces the measured nonthermal intensity changes for all Bragg peaks at 3 ps, and that the PDF analysis is only qualitative.

Significance. If validated, the claimed observation of ultrafast photoinduced Mo-Mo dimerization — rather than the bond elongation reported in some previous work — would be a significant contribution to the understanding of Peierls distortion dynamics in 1T'-MoTe2 and related TMDCs. The paper has clear strengths: sub-picosecond UED with multi-fluence control, a transparent structure-factor framework, explicit three-stage temporal evolution, and DFT/COHP calculations that connect the proposed motions to the electronic structure. The shear-displacement component is reasonably supported by the interplanar spacing data. However, the dimerization component, which is the most novel part of the claim, is not supported at the same level: the intensity fit fails globally, the bond-shortening interpretation rests on a single (020) spacing measurement with no quantitative displacement model, and the PDF resolution is admitted to be insufficient to resolve individual Mo-Mo pairs. The paper's own stated limitations therefore place the central claim only at the level of a plausible interpretation, not a quantitative demonstration.

major comments (4)
  1. [Results and discussions, Fig. 2(c) and the text following Eq. (1)] The paper states that none of the five distortion models reproduces the measured nonthermal intensity changes for all Bragg peaks at 3 ps, and that the combination of shear displacement and Mo-Mo bond shortening is selected on the basis of the (310) peak alone. Because the intensity analysis is presented as a primary support for the two-motion model, the failure of the global fit means that alternative combinations of thermal and nonthermal contributions cannot be excluded. This weakens the central structural assignment unless the interplanar-spacing evidence is independently quantitative.
  2. [Results and discussions, Fig. 2(e) and Fig. 3 (PDF analysis)] The Mo-Mo bond-shortening claim rests on the rapid shrink of the (020) interplanar spacing plus a qualitative PDF. No calculated Δd/d for the proposed bond-shortening displacement, no explicit bond-length change, and no quantitative comparison with alternatives such as anisotropic thermal strain or defect-induced local distortions are provided. The text itself states that the PDF has about 1 Å resolution so that no single atom pair separation is resolved, and that only qualitative discussion is possible. Consequently, the dimerization component of the central claim is not quantitatively established.
  3. [Results and discussions, thermal-calibration paragraph and Fig. 4] The (020) intensity drop is used to calibrate the thermal (Debye-Waller) contribution for all other Bragg peaks, based on the assumption that the (020) structure factor is unaffected by shear displacement and Mo-Mo bond shortening. However, later in the same section, the paper invokes a b-c plane distortion at 3 ns that also depresses the (020) intensity. If the (020) peak possesses any nonthermal response at early times, the subtracted nonthermal intensities in Fig. 2(c) and all model comparisons derived from them are miscalibrated. This is a partial circularity that the manuscript does not resolve.
  4. [Fig. 5(c) and the COHP discussion] The COHP analysis as stated predicts that photoexcited electrons occupying antibonding Mo d orbitals should elongate the Mo-Mo metallic bond, which is the opposite of the observed shortening. The rebuttal invokes subtle COHP features near the Fermi level, a small photoexcited density of about 0.015 e/cell, and the instrument's temporal resolution, but no quantitative calculation of the change in integrated COHP or of the photoexcited orbital occupancy is presented. The proposed electronic mechanism for dimerization is plausible but is not demonstrated to the same standard as the shear-displacement mechanism.
minor comments (7)
  1. [Eq. (1) in the main text] Equation (1) appears garbled in the manuscript ('*f!exp' and '4hx!+ky!+lz!9'); it should be typeset with a proper summation and correctly placed brackets for the fractional coordinates.
  2. [Supplementary Note 1, Eq. (2)] Equation (2) in the supplementary information has unbalanced parentheses and an unclear average structure; please correct the expression for S(Q).
  3. [Abstract and Introduction] The abstract contains grammatical issues such as '1T'-MoTe2 possess' and 'The desired details of pathway and time span are still insufficient'; these should be revised for clarity.
  4. [Fig. 2(c) caption and main text] The atomic displacement used in the distortion models is quoted as '~1‱' in the caption; please specify the value in absolute units (Å or fractional coordinates) and state whether it is a rigid displacement or a per-atom parameter.
  5. [Fig. 2(e)] The decrease of β by about 0.5° is reported without an uncertainty estimate; please provide the fit uncertainty or confidence interval for this angle.
  6. [Fig. 3(b) caption] The arrows in Fig. 3(b) are described as guides to the eye; please state explicitly what feature each arrow marks (e.g., peak shift or valley shift) so that the three-stage evolution can be followed by the reader.
  7. [References] Reference [31] is an arXiv preprint; if a published version exists, it should be cited. The supplementary reference numbering also restarts at [37] and is disconnected from the main-text list; please unify the numbering.

Circularity Check

1 steps flagged · score 4.0 of 10

Partial circularity in the nonthermal intensity analysis: the (020) thermal calibration assumes the shear + Mo-Mo bond-shortening model that the paper concludes, and the paper later invokes (020) nonthermal response at 3 ns; d-spacing and PDF evidence remain independent.

  1. self definitional [Results and discussions, paragraph following Eq. (1) and Fig. 2(c)]
    "Given that the structure factor of (020) crystal plane is not affected by these two motions, the intensity drop of (020) peak is considered constituting only of the contribution of increased MSD, which is better demonstrated by the measurement with varied fluences, as discussed later."

    The nonthermal intensity changes used to identify the shear + Mo-Mo bond-shortening model are defined by subtracting a thermal contribution calibrated on the assumption that these two motions do not affect the (020) plane. That assumption is exactly the defining property of the concluded model: the selected distortions leave the (020) structure factor unchanged. Therefore the intensity dataset on which the model is selected is constructed from the model's key premise. The paper later undermines this premise at 3 ns by interpreting extra (020) intensity drop from a b-c plane distortion, showing (020) is not intrinsically insensitive to nonthermal motion.

full rationale

The paper's central interpretation rests on three lines: nonthermal Bragg intensity analysis, interplanar-spacing kinetics, and qualitative PDF. The intensity line is partially circular: the nonthermal changes are defined by subtracting a thermal contribution scaled to the (020) peak, under the explicit assumption that shear and Mo-Mo bond motions do not affect (020). That assumption is the defining property of the combined model selected from these very nonthermal intensities, so the input is constructed from the conclusion. The paper itself later states that a b-c plane distortion brings additional intensity drop to the (020) peak, abandoning the assumption at 3 ns without establishing that it holds at 3 ps. The shear displacement is also supported by the d-spacing pattern at 3 ps, but the magnitude (beta reduced by ~0.5 degrees) is a fitted parameter, not a derived prediction. The PDF discussion is explicitly qualitative. The DFT/COHP section is an ex post rationalization rather than a parameter-free prediction, but it is not circular. Self-citations (refs [17] and [31]) are minor and not load-bearing. Overall, the central claim is not forced by definition or by a self-citation chain, and the d-spacing and PDF observations provide independent content, so the circularity is partial rather than total.

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

The central interpretation rests on several fitted or chosen parameters (shear angle, atomic displacement, U, carrier density) and on the assumption that the (020) peak is purely thermal. The DFT provides an independent but qualitative narrative. None of these are derived from first principles in this paper.

free parameters (4)
  • shear angle reduction beta = ~0.5 degrees at 3 ps
    The monoclinic angle between a- and c-axis is reduced by about 0.5 degrees to make the shear displacement model match the observed interplanar spacing changes at 3 ps (Fig. 2(e)). This value is chosen to fit the data, not derived from an independent calculation.
  • atomic displacement magnitude for distortion models = ~1e-4 (relative units)
    The five distortion models in Fig. 2(c) are calculated with an atomic displacement of approximately 1e-4 of the lattice parameter. This amplitude is chosen by hand to compare relative intensity changes, not determined from the data.
  • Hubbard U for Mo 4d = 5.0 eV
    The DFT+U calculation sets U=5.0 eV for Mo 4d orbitals, following ref [7]. This choice affects the band structure, DOS, and COHP used to explain the observed distortions, but it is taken from prior literature, not derived in this work.
  • photoexcited carrier density = ~0.015 e/cell
    The paper estimates the photodoped electron concentration (about 0.015 e/cell) in the discussion of Mo-Mo bond shortening. The estimate is used to argue that the excitation is weak compared with the ~0.4 e/cell needed for bond elongation, but the conversion from fluence to carrier density is not detailed.
assumptions (5)
  • standard math The structure factor depends only on fractional atomic coordinates, and the Debye-Waller model describes thermal intensity loss as a linear function of squared scattering vector.
    Used in Eq. (1) and in the analysis of Fig. 2(b) to separate thermal from nonthermal intensity changes. These are standard results in diffraction physics.
  • ad hoc to paper The (020) Bragg peak intensity drop is caused solely by increased atomic mean-square displacement (thermal effect) and is unaffected by shear displacement or Mo-Mo bond shortening.
    Stated in the paragraph after Eq. (1): 'the intensity drop of (020) peak is considered constituting only of the contribution of increased MSD.' This assumption is load-bearing for the thermal subtraction, and the paper itself later suggests (020) may have nonthermal contributions on the nanosecond timescale.
  • domain assumption The five distortion models considered (combinations of shear and Mo-Mo bond changes) span the relevant structural degrees of freedom.
    The comparison in Fig. 2(c) assumes that the true distortion is among the tested models. The failure to match all peaks is attributed to defects/polarons, so the model set is taken as sufficient for the central interpretation.
  • domain assumption GGA-PBE+U with U=5.0 eV and DFT-D3 gives a reliable description of the band structure, DOS, and COHP of 1T'-MoTe2.
    The DFT results (band structure with partially occupied Te p band near the Fermi level, positive COHP-Mo above the Fermi level) are used to explain the observed shear and bond shortening. This is a standard but approximate method, and the conclusions depend on its accuracy.
  • domain assumption The UED probe along the c-axis and the quasi-single-crystal morphology allow PDF analysis in the a-b plane only, and the reduced-diffraction-element construction is valid.
    The PDF model reduces atoms to columns projected in the a-b plane (Eq. (2), Fig. S1), which is necessary to interpret the measured PDF. Any error in this reduction affects the qualitative PDF interpretation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Ultrafast atomic dimerization of Peierls distortion in semimetal molybdenum ditelluride." pith.science (2026). https://pith.science/paper/N7JYKSWH

@misc{pith2026250521096,
  author       = {Pith},
  title        = {Pith review of: Ultrafast atomic dimerization of Peierls distortion in semimetal molybdenum ditelluride},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/N7JYKSWH}},
  note         = {Machine review of arXiv:2505.21096}
}
abstract

Semimetal molybdenum ditelluride (1T'-MoTe$_2$) possess diverse phase transitions enriching its application prospects. The structural response during these transitions is crucial to understanding the underlying mechanisms, but the desired details of pathway and time span are still insufficient. Here, we investigate the lattice evolution in few-layer 1T'-MoTe$_2$ after photoexcitation, using ultrafast electron diffraction and density functional theory (DFT) calculations. The observed complex lattice responses with unintuitively evolving Bragg peak intensity and interplanar spacing, are best interpreted as the combination of shear displacement and Mo-Mo bond shortening in a few picoseconds, and a metastable structure in nanoseconds, basing on the analyses of structure factor and pair distribution function. The DFT calculations reveal that, the photodoped electrons induced population change of the antibonding states close to Fermi level, lead to the shear displacement and the dimerization of Mo pairs. Our findings present new insights for elucidating the picture of Peierls distortion in 1T'-MoTe$_2$.

Figures

Figures reproduced from arXiv: 2505.21096 by the authors.

Figure 1
Figure 1. Sample characterization, lattice structure, and UED measurement. (a) Monoclinic lattice of 1T’-MoTe2 with layered structure bonding by van der Waals force. Dashed frames, unit cell. Lower right panel, high-resolution transmission electron microscopy image with a crystalline grain outlined in white, and the schematic of grain orientations. (b) Raman spectrum. Inset, selected area electron diffraction pattern. (c) Sch… view at source ↗
Figure 5
Figure 5. Calculated band structure, DOS, and COHP. (a) Orbital-projected electronic band. The size of marker denotes the contribution strength of each orbital. The s orbitals of Mo are mainly located below -1.5 eV, thus hardly visible. (b) DOS of bulk 1T′-MoTe2. (c) COHP of Mo-Mo bond. The d orbitals of Mo with energy higher than -1.15 eV are antibonding states. Dashed line, the Fermi Level. We also calculated the COHP, whic… view at source ↗

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

42 extracted references · 39 canonical work pages

  1. [1]

    Manzeli S et al 2017 Nat. Rev. Mater. 2 17033

    work page 2017

  2. [2]

    Keum D H et al 2015 Nat. Phys. 11 482

    work page 2015

  3. [3]

    Zhu M et al 2019 Small 15 1903159

    work page 2019

  4. [4]

    Tang F, Po H C, Vishwanath A and Wan X 2019 Nat. Phys. 15 470

    2019

  5. [5]

    Wang Z et al 2019 Phys. Rev. Lett. 123 186401

    work page 2019

  6. [6]

    Su B et al 2022 Adv. Sci. 9 2101532

    work page 2022

  7. [8]

    Sakai H et al 2016 Sci. Adv. 2 e1601378

    work page 2016

  8. [9]

    Paul S et al 2020 Phys. Rev. B 102 054103

    work page 2020

Show all 42 references
  1. [10]

    Zhuang M Y et al 2019 Phys. Rev. X 9 021036

    2019

  2. [11]

    Fukuda T et al 2020 Appl. Phys. Lett. 116 093103

    2020

  3. [12]

    Qi Y et al 2022 ACS Nano 16 11124

    2022

  4. [13]

    Sie E J et al 2019 Nature 565 61

    2019

  5. [14]

    Xu C et al 2023 Nat. Commun. 14 1256

    2023

  6. [15]

    Deng K et al 2016 Nat. Phys. 12 1105

    2016

  7. [16]

    Chen S et al 2016 Nano Lett. 16 5852

    2016

  8. [18]

    Wisser J J et al 2025 Phys. Rev. Applied 23 034075

    2025

  9. [19]

    Schäfer S, Liang W and Zewail A H 2011 Chem. Phys. Lett. 515 278

    2011

  10. [20]

    Hu J et al 2023 Adv. Funct. Mater. 33 2206395

    2023

  11. [21]

    Dai Y et al 2021 Nano Lett. 21 9903

    2021

  12. [22]

    Chen J et al 2022 J. Phys. Chem. C 126 9407

    2022

  13. [23]

    Phys.: Condens

    Choi J and Jhi S 2018 J. Phys.: Condens. Matter 30 105403

    2018

  14. [24]

    Sahu R et al 2017 J. Appl. Phys. 121 105101

    2017

  15. [25]

    Yang J et al 2021 Nature 596 531

    2021

  16. [26]

    Morrison V R et al 2014 Science 346 445

    2014

  17. [27]

    Wu X et al 2017 Sci. Adv. 3 e1602388

    2017

  18. [29]

    Kresse G and Furthmüller J 1996 Comput.l Mater. Sci. 6 15

    1996

  19. [30]

    Ferry D K 2021 Appl. Phys. Rev. 8 021324

    2021

  20. [31]

    Wang Z et al 2025 arXiv:2502.16145 [cond-mat.mtrl-sci] 9

    2025 arXiv

  21. [32]

    Han X et al 2021 ACS Appl. Mater. Interfaces 13, 44703

    2021

  22. [34]

    Canadell E et al 1989 J. Am. Chem. Soc. 111 3778

    1989

  23. [35]

    Shi H et al 2013 ACS Nano 7 1072

    2013

  24. [36]

    #(!$%𝑓&-(𝑄.)⟩𝑑𝑄. ∫𝐼/.(𝑄.)𝑑𝑄.(!

    Britz A et al 2021 Struct. Dynam. 8 014501 10 Supplementary Information Ultrafast atomic dimerization of Peierls distortion in semimetal molybdenum ditelluride Zhong Wang(王众)1,2, Chunlong Hu(胡春龙)1,2, Changchang Gong(巩畅畅)1,2, Fuyong Hua(华傅甬)1,2, Qian You(游骞)1,2,†, Wenxi Liang(梁...

    2021

  25. [37]

    B 30 090701

    Hu C et al 2021 Chinese Phys. B 30 090701

    2021

  26. [38]

    Kresse G and Hafner J 1993 Phys. Rev. B 47 558

    1993

  27. [39]

    Blöchl P E 1994 Phys. Rev. B 50 17953

    1994

  28. [40]

    Grimme S et al 2010 J. Chem. Phys. 132 154104

    2010

  29. [41]

    Kim H, Kang S, Hamada I and Son Y 2017 Phys. Rev. B 95 180101

    2017

  30. [42]

    Nelson R et al 2020 J. Comput. Chem. 41 1931

    2020

  31. [43]

    Shi H, Luo M and Wang W 2019 Comput. Phys. Commun. 238 295

    2019

  32. [44]

    Tran D T, Svensson G and Tai C 2017 J. Appl. Crystallogr. 50 304

    2017

  33. [45]

    A 52 257

    Peng L, Ren G, Dudarev S L and Whelan M J 1996 Acta Cryst. A 52 257

    1996

  34. [46]

    Optical Mater

    Shi Y et al 2025 Adv. Optical Mater. https://doi.org/10.1002/adom.202403141

    2025 doi

Pith tools

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