Pith. sign in

REVIEW 3 major objections 5 minor 27 references

This paper establishes that in misfit layer compounds the La/Pb ratio acts as a chemical gate: it continuously dopes a NbSe2 monolayer and switches its charge density wave from a 3x3 to a 2x2 pattern, with long-range order vanishing at high

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

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

T0 review · deepseek-v4-flash

2026-08-01 13:26 UTC pith:ZHNW5WLG

load-bearing objection Genuinely new mapping of the misfit doping series, with a credible 3x3->2x2 transition; the QPI-derived doping axis is the weak point, but not a fatal one. the 3 major comments →

arxiv 2607.19095 v1 pith:ZHNW5WLG submitted 2026-07-21 cond-mat.mtrl-sci

Doping tunable charge density waves in misfit layer compounds

Hugo Le Du , Ludovica Zullo , Justine Cordiez , Robin Salvatore , Daniel Schmieg , Arindam Mukherjee , Giovanni Marini , Dominik Volavka
show 9 more authors
Francois Debontridder Marie Herve Tomas Samuely Shunsuke Sasaki Florent Pawula Etienne Janod Matteo Calandra Laurent Cario Tristan Cren
This is my paper
classification cond-mat.mtrl-sci
keywords charge density wavesmisfit layer compoundsNbSe2doping controlscanning tunneling microscopyquasiparticle interferencedensity functional theory2x2 CDW
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

This paper tries to show that a misfit layer compound—a naturally stacked sandwich of a rocksalt layer and a NbSe2 layer—can act as a chemical gate for the NbSe2 sheet, with the La/Pb ratio in the rocksalt layer setting the electron doping. Using density-functional phonon calculations and low-temperature scanning tunneling microscopy, the authors find that increasing doping drives the charge density wave (CDW) of NbSe2 from the usual 3x3 pattern to a 2x2 pattern, through a regime where both coexist, and eventually suppresses long-range order altogether. The result matters because uniform, sample-scale doping of a single TMD layer has been hard to achieve by gating or intercalation; here stoichiometry alone selects the CDW ordering vector. If correct, misfits offer a synthesis-level knob for engineering quantum phases in two-dimensional materials.

Core claim

The central claim is that in (LaxPb1−xSe)1.14(NbSe2)2 the electron count on the NbSe2 monolayer is set by the La/Pb ratio—La3+ donates charge, Pb2+ does not—and that this continuous doping axis controls the CDW ground state. STM images at 4.2 K show a pure 3x3 CDW at x=0, coexisting 3x3 and 2x2 domains at x=0.15–0.3, a pure 2x2 phase at x=0.4 (≈0.28 e−/Nb), and a short-range 3x3 resurgence at x=0.6 (≈0.35 e−/Nb), while DFPT and SSCHA calculations place the 3x3→2x2 crossover near 0.2–0.3 e−/Nb and the complete collapse of the CDW near 0.4–0.5 e−/Nb. The authors also show that the 2x2 order is not driven by Fermi-surface nesting but likely by electron–phonon coupling, and that misfit strain br

What carries the argument

The load-bearing object is the misfit heterostructure itself: the rocksalt (LaxPb1−xSe) layer acts as a built-in gate that transfers electrons into the adjacent NbSe2 monolayer, and alloying La for Pb sets the transferred charge per Nb atom. Two probe elements carry the experimental argument: low-temperature STM topographs reveal the CDW wavevector directly, and quasiparticle-interference (QPI) contours—hexagonal ring sizes fit to T-matrix simulations on a DFT band structure—convert each image into a local doping value. On the theory side, DFPT phonon calculations on a field-effect model of the misfit predict the CDW instabilities, with SSCHA anharmonic corrections fixing the critical doping

Load-bearing premise

The entire doping axis—and with it the placement of the 3x3→2x2→collapse sequence—rests on the assumption that the size of the hexagonal quasiparticle-interference contour, matched to T-matrix simulations built on the authors' DFT band structure, faithfully reports the electron count in the exact region where each CDW is imaged.

What would settle it

Take the same crystals used for STM and measure the Fermi-surface area directly by angle-resolved photoemission (or by the dI/dV onset of the conduction band); if the spectroscopically measured doping disagrees with the QPI-inferred values, then the CDW phase boundaries attributed to specific electron counts would need to be re-anchored.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Tuning the La/Pb ratio provides uniform, sample-scale electron doping of the NbSe2 monolayer, reaching levels (≈0.6 e−/Nb at x=1) inaccessible to electric-field gating.
  • The CDW ordering vector responds continuously to that doping: pure 3x3 at x=0, coexisting 3x3 and 2x2 around x=0.15–0.3, pure 2x2 at x≈0.4, and a short-range 3x3 remnant at x≈0.6.
  • DFPT and SSCHA calculations locate the 3x3→2x2 crossover near 0.2–0.3 e−/Nb and the collapse of the CDW near 0.4–0.5 e−/Nb, meaning long-range order vanishes before the fully La end-member is reached.
  • The 2x2 order is not driven by Fermi-surface nesting; the electronic susceptibility shows no nesting at M, pointing to electron–phonon coupling as the mechanism.
  • Because misfit strain breaks the threefold symmetry and becomes decisive near the compositional transition, the CDW phase is simultaneously a sensitive probe of strain in this platform.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The doping calibration assumes the QPI contour size maps rigidly to band filling; since the x=0 compound already shows PbSe-derived states that distort the QPI pattern, a photoemission or transport cross-check on identical crystals would show whether the quoted electron counts at each CDW boundary are accurate.
  • If strain is what stabilizes the 2x2 phase only when the system is near the doping-induced transition, then applying controlled uniaxial strain to undoped NbSe2 should induce a 2x2 CDW at lower electron counts—an experiment the paper does not report but its logic implies.
  • The short-range 3x3 order at x=0.6 could be disorder-stabilized patches rather than a true bulk phase; temperature-dependent measurements would tell whether they are fluctuating remnants of the CDW or pinned impurity states.
  • The same misfit architecture has been used to tune superconductivity in related compounds, so the CDW phase diagram mapped here invites a search for a doping-controlled CDW–superconductivity interplay in NbSe2 itself.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 5 minor

Summary. The paper reports a combined density-functional and scanning tunneling microscopy study of the misfit layer compounds (LaxPb1−xSe)1.14(NbSe2)2 for x = 0, 0.15, 0.3, 0.4, and 0.6. The authors argue that increasing the La/Pb ratio electron-dopes the NbSe2 layers and drives the CDW from the canonical 3×3 order of undoped monolayer NbSe2, through a coexistence region, to a pure 2×2 order, and finally to a suppression of long-range CDW order at high doping. The theoretical part uses a field-effect model for the misfit, harmonic DFPT, and SSCHA to map phonon instabilities versus doping; the experimental part uses STM topographies and FFTs to identify the CDW order, and QPI contour sizes to estimate local carrier densities. The paper also attributes the anisotropic 2×2 component at low/intermediate doping to the uniaxial misfit strain.

Significance. If validated, the central claim is significant: it would provide a continuous, sample-scale chemical doping axis for CDWs in TMD monolayers, reaching electron densities far beyond electrostatic gating, and would establish misfit layer compounds as a designable platform for CDW engineering. The paper combines multiple methods—DFPT/SSCHA, STM imaging, QPI simulations, and a susceptibility check—and the CDW ordering itself is directly observed rather than inferred. The identification of a 2×2 CDW in NbSe2 is a notable result. However, the quantitative doping axis that links experiment to theory is under-validated, and the high-doping endpoint is not cleanly confirmed by the data as presented.

major comments (3)
  1. [Sec. III, Table I and Fig. S6] The QPI-based doping values (0.12, 0.21, 0.28, 0.35 e−/Nb) are extracted by fitting the size of the hexagonal QPI contour to rigid-band T-matrix simulations built on the authors' own DFT band structure. The authors state in Sec. III that for x=0 the rigid-band QPI description is invalid because PbSe-derived states modify the spectral weight, yet x=0 serves as the 0.0 e−/Nb reference and no error bars or cross-check against the ARPES-determined global doping of Ref. [14] are provided. The extracted values also differ from the nominal 0.57x relation used in the Introduction (e.g., 0.12 vs 0.086 for x=0.15; 0.28 vs 0.228 for x=0.4), and this discrepancy is not discussed. Since the paper's thesis is that doping, not La content per se, controls the CDW, this calibration uncertainty directly affects the claimed agreement with the DFT phase diagram in Table I.
  2. [Sec. II SSCHA paragraph; Introduction] The paper states that 'the compound with 50% of the La content' has 'a corresponding charge transfer of 0.4 electrons per Nb atoms,' but the Introduction defines the charge transfer as 0.57x electrons per Nb (Ref. [14]), which gives 0.285 e/Nb for x=0.5. Similarly, '60–80% La content' is said to give '≈0.5 e/Nb', while 0.57×0.6 = 0.342 and 0.57×0.8 = 0.456. This internal inconsistency in the x-to-doping mapping affects the claimed critical doping for CDW collapse and makes it difficult to compare the DFT phase diagram with the STM doping values in Table I, which do not follow the 0.57x relation either. The authors should clarify which calibration is used where.
  3. [Sec. III, high-doping regime; Table I] At x=0.6, the QPI-derived doping in Table I is 0.35 e−/Nb, below the calculated collapse threshold of 0.4 e−/Nb. The observed short-range 3×3 order is therefore not the predicted 2×2 (expected near 0.3 e−/Nb) nor the collapsed state (expected for ≥0.4 e−/Nb). The text quotes a nominal '≈0.5 e−/Nb' for this sample to support the collapse interpretation, contradicting the QPI value in the same table. If the local doping is 0.35 e−/Nb, the presence of 3×3 order is a discrepancy with the DFT phase diagram, not a confirmation of collapse. A quantitative explanation—e.g., disorder-broadened transitions or a refined charge-transfer model—is required.
minor comments (5)
  1. [Sec. III] In the sentence 'the observed CDW evolve systematically with doping,' the verb should agree with the singular subject 'CDW' ('evolves').
  2. [Fig. 4 caption] The caption states that Fermi levels show systematic shifts but does not specify that these shifts are extracted from QPI contour sizes, not directly from topography. This should be stated explicitly to avoid confusion.
  3. [Table I] The DFT and STM rows both use the heading 'Doping (e−/Nb)' but contain values on different scales (DFT: 0.0–0.5; STM: 0.0–0.35). Add a note or subheading indicating the different provenance of these numbers, and list the corresponding La content x for the DFT entries.
  4. [Sec. II and Sec. III] The phrase 'the compound with 50% of the La content' and '60% to 80% La content' should consistently use the x notation defined in the Introduction, since 'La content' could be confused with the doping level (e-/Nb).
  5. [Sec. II] The term 'quasi-commensurate' is used without a definition. Define it in terms of the observed incommensurate lattice peaks and CDW peak positions in the FFT.

Circularity Check

1 steps flagged

Doping axis relies on authors' own prior charge-transfer calibration and DFT-based QPI fitting, but CDW observations are direct and not forced.

specific steps
  1. self citation load bearing [Section I (Introduction) and Section II (QPI doping extraction)]
    "Based on our DFT+ARPES study [12,14], it is possible to increase the charge transfer to NbSe2 by increasing the percentage of La in the (LaxPb1−xSe)1.14(NbSe2)2 samples... a typical charge transfer of 0.57x electrons per niobium atom in the TMD layers [14]."

    The quantitative doping axis of the phase diagram is imported from the authors' own prior preprint (ref [14]) rather than derived in this work. The local doping values are then estimated by matching QPI contours to T-matrix simulations built on the same group's DFT band structure (refs [13,20]) and compared with DFT phonon predictions at those same doping values. This makes the experiment-theory agreement in Table I partly a consistency check within the authors' own theoretical framework rather than an independent confirmation of the doping values. The observed CDW patterns themselves are direct STM data and are not fitted, so the central claim is not forced by construction.

full rationale

The central observation—different CDW orders (3x3, 2x2, coexistence, short-range order) at different La/Pb ratios—is direct STM topography and is not derived from the doping calibration. The QPI-based doping estimates are an independent observable (Fermi-surface contour size) and are not identical to the nominal 0.57x values, so the mapping is not a tautology. However, the paper leans heavily on self-citations: the charge-transfer calibration (0.57x e/Nb) comes from the same authors' unpublished work (ref [14]), and the QPI simulations use the same group's DFT band structure. The paper itself concedes that the rigid-band QPI analysis fails at the x=0 reference point due to PbSe-derived states, which undermines confidence in the absolute doping axis. This is a correctness/validation weakness rather than a logical circularity, but it lowers the independence of the claimed DFT-experiment agreement.

Axiom & Free-Parameter Ledger

3 free parameters · 3 axioms · 0 invented entities

The paper introduces no new physical entities. It relies on the authors' previously developed charged-gate model for misfits, a fitted charge-transfer coefficient, and structural parameters from their prior work; the QPI doping estimates are fit outputs, not measured inputs.

free parameters (3)
  • Charge transfer coefficient (0.57 e-/Nb per x) = 0.57 electrons per Nb per unit La content x
    Used to map composition to doping; taken from the authors' earlier ARPES/DFT study (ref 14); not re-derived here.
  • Gate geometry in the 1L-FET model = optimized gate positions listed in Tab. S1
    Free parameters tuned to reproduce the charge transfer of full misfit calculations (refs 13,15).
  • Strain lattice parameters = a1=3.457 Å, a2=3.437 Å, angle=59.6°
    Taken from experimental structural study (ref 12); chosen to reflect the misfit compression and central to the symmetry-breaking explanation for the 2x2 modulation direction.
axioms (3)
  • domain assumption The misfit effect on the TMD layer can be modeled as a uniform positively charged gate plus a repulsive potential barrier (single-gate field-effect setup).
    Invoked in Sec. II ('modeling the misfit as a collection of field-effect transistors') and validated only against the authors' own prior calculations (refs 13,15).
  • domain assumption Rigid-band approximation: doping only shifts the Fermi level of monolayer NbSe2 without modifying band structure or electronic wavefunctions relevant to QPI.
    Used in Sec. II QPI simulations; explicitly acknowledged to fail for x=0 (PbSe) due to PbSe-derived bands crossing the Nb-d bands.
  • domain assumption DFPT harmonic phonons plus SSCHA anharmonicity at T=0 reliably capture the CDW instability and its critical doping.
    The SSCHA correction is only computed for 50% La; other compositions rely on harmonic calculations (Fig. 2, Sec. II), so the anharmonic boundary at 0.4 e-/Nb is an extrapolation.

pith-pipeline@v1.3.0-alltime-deepseek · 10632 in / 14775 out tokens · 131216 ms · 2026-08-01T13:26:55.432354+00:00 · methodology

0 comments
read the original abstract

The ability to tune charge density waves (CDWs) through external control knobs, such as doping, pressure or strain is crucial for exploring the phase diagram of two dimensional (2D) or quasi-2D materials. Yet, controlling CDWs critical temperature and ordering vector remains a challenge for current experimental techniques. In this work, we establish misfit layer compound heterostructures as a reliable platform to manipulate CDWs in transition metal dichalcogenides. By combining ab initio calculations with low-temperature scanning tunneling microscopy, we show how to achieve doping tunable control over NbSe2 CDW by chemically alloying in the rocksalt subunit. Crucially, we prove that tuning the La Pb ratio in the misfit family (LaxPb1xSe)1.14(NbSe2)2 enables stabilization of different CDW orders, such as 2x2 or 3x3 patterns, and even coexisting phases. This work paves the way for engineering transition metal dichalcogenides with tailored charge density waves within misfit heterostructures.

Figures

Figures reproduced from arXiv: 2607.19095 by Arindam Mukherjee, Daniel Schmieg, Dominik Volavka, Etienne Janod, Florent Pawula, Francois Debontridder, Giovanni Marini, Hugo Le Du, Justine Cordiez, Laurent Cario, Ludovica Zullo, Marie Herve, Matteo Calandra, Robin Salvatore, Shunsuke Sasaki, Tomas Samuely, Tristan Cren.

Figure 1
Figure 1. Figure 1: a) Side-view of (LaxPb1−xSe)1.14(NbSe2)2. As an example, the compound with x = 0.5 is shown. b) The field￾effect modeling scheme for the surface of misfit layer compounds. In this scheme, MLC surfaces can be modeled in a single layer field effect transistor setup (1L FET) [16, 17] by replacing the rocksalt with a positively charged gate, a negatively charged monolayer TMD then corresponds to the misfit ter… view at source ↗
Figure 2
Figure 2. Figure 2: Theoretical predictions for the doping evolution of the CDW of NbSe [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: a) Fourier filtered STM topography of (PbSe) [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: a-c) Fourier filtered STM topographies of different (La [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

27 extracted references · 1 linked inside Pith

  1. [1]

    Revolinsky, G

    E. Revolinsky, G. Spiering, and D. Beerntsen, Journal of Physics and Chemistry of Solids26, 1029 (1965)

    1965

  2. [2]

    Wilson, F

    J. Wilson, F. Di Salvo, and S. Mahajan, Advances in Physics24, 117 (1975)

    1975

  3. [3]

    J. M. Lu, O. Zheliuk, I. Leermakers, N. F. Q. Yuan, U. Zeitler, K. T. Law, and J. T. Ye, Science350, 1353 (2015)

    2015

  4. [4]

    Saito, Y

    Y. Saito, Y. Nakamura, M. S. Bahramy, Y. Kohama, J. Ye, Y. Kasahara, Y. Nakagawa, M. Onga, M. Toku- naga, T. Nojima, Y. Yanase, and Y. Iwasa, Nature Physics12, 144 (2016)

    2016

  5. [5]

    S. Gao, F. Flicker, R. Sankar, H. Zhao, Z. Ren, B.Rachmilowitz, S.Balachandar, F.Chou, K.S.Burch, Z. Wang, J. van Wezel, and I. Zeljkovic, Proceedings of the National Academy of Sciences115, 6986 (2018)

    2018

  6. [6]

    A. K. Kundu, A. Rajapitamahuni, E. Vescovo, I. I. Klimovskikh, H. Berger, and T. Valla, Communica- tions Materials5, 208 (2024)

    2024

  7. [7]

    Ørsted, A

    A. Ørsted, A. Scarfato, C. Barreteau, E. Giannini, and C. Renner, Nano Letters25, 1729 (2025)

    2025

  8. [8]

    S. Das, A. Sebastian, E. Pop,et al., Nature Electronics 4, 786 (2021)

    2021

  9. [9]

    J. Liu, B. Li, and Q. Li, Magnetochemistry8, 172 (2022)

    2022

  10. [10]

    Li and F

    Z. Li and F. Chen, Applied Physics Reviews4, 011103 (2017)

    2017

  11. [11]

    Rouxel and A

    J. Rouxel and A. Meerschaut, Molecular Crystals and Liquid Crystals244(1994)

    1994

  12. [12]

    R. T. Leriche, A. Palacio-Morales, M. Campetella, C. Tresca, S. Sasaki, C. Brun, F. Debontridder, 9 P. David, I. Arfaoui, O. Šofranko, T. Samuely, G. Kre- mer, C. Monney, T. Jaouen, L. Cario, M. Calan- dra, and T. Cren, Advanced Functional Materials31, 2007706 (2021)

    2021

  13. [13]

    Zullo, G

    L. Zullo, G. Marini, T. Cren, and M. Calandra, Nano Letters23, 6658 (2023)

    2023

  14. [14]

    Tuning the charge transfer of transition metal dichalcogenides via misfit layer compounds,

    H. Le Du, L. Zullo, J. Cordiez, R. Salvatore, G. Marini, M. Hervé, D. Pierucci, S. Sasaki, F. Pawula, E. Janod, C. Bigi, M. Zonno, F. Bertran, T. Jaouen, P. L. Fèvre, M. Calandra, L. Cario, T. Cren, and M. D. An- gelo, “Tuning the charge transfer of transition metal dichalcogenides via misfit layer compounds,” (2026), arXiv:2605.18317 [cond-mat.mtrl-sci]

    Pith/arXiv arXiv 2026

  15. [15]

    Zullo, G

    L. Zullo, G. Setnikar, A. Pawbake, T. Cren, C. Brun, J. Cordiez, S. Sasaki, L. Cario, G. Marini, M. Calandra, and M.-A. Méasson, Physical Review B110, 075430 (2024)

    2024

  16. [16]

    Brumme, M

    T. Brumme, M. Calandra, and F. Mauri, Physical Re- view B91, 155436 (2015)

    2015

  17. [17]

    Sohier, M

    T. Sohier, M. Calandra, and F. Mauri, Physical Review B96, 075448 (2017)

    2017

  18. [18]

    Bianco, L

    R. Bianco, L. Monacelli, M. Calandra, F. Mauri, and I. Errea, Physical Review Letters125, 106101 (2020)

    2020

  19. [19]

    Monacelli, R

    L. Monacelli, R. Bianco, M. Cherubini, M. Calandra, I. Errea, and F. Mauri, Journal of Physics: Condensed Matter33, 363001 (2021)

    2021

  20. [20]

    Zullo,Emergent properties of misfit layer compounds from first principles, Ph.D

    L. Zullo,Emergent properties of misfit layer compounds from first principles, Ph.D. thesis, Sorbonne Université and Università degli Studi di Trento (2024)

    2024

  21. [21]

    W.-C. Chiu, S. Mardanya, R. Markiewicz, J. Niemi- nen, B. Singh, T. Hakioglu, A. Agarwal, T.-R. Chang, H. Lin, and A. Bansil, ACS Nano19, 18108 (2025)

    2025

  22. [22]

    Luque Merino, F

    R. Luque Merino, F. Carrascoso, E. Henríquez-Guerra, M. R. Calvo, R. Frisenda, and A. Castellanos-Gomez, Nano Letters25, 16379 (2025)

    2025

  23. [23]

    Zhang, J

    D. Zhang, J. Ha, H. Baek, Y.-H. Chan, F. D. Natterer, A. F. Myers, J. D. Schumacher, W. G. Cullen, A. V. Davydov, Y. Kuk, M. Y. Chou, N. B. Zhitenev, and J.A.Stroscio,PhysicalReviewMater.1,024005(2017)

    2017

  24. [24]

    Moiré-induced symmetry breaking of charge order in van der waals heterostructures,

    S. Sajan, L. Pätzold, T. Agarwal, C. Pfister, H. Guo, S. Duan, P. V. Sruthibhai, M. Rossi, M. N. Gas- tiasoro, S. Barja, R. P. Singh, T. Wehling, and M. M. Ugeda, “Moiré-induced symmetry breaking of charge order in van der waals heterostructures,” (2026), arXiv:2603.05759 [cond-mat.mes-hall]

    arXiv 2026

  25. [25]

    C. J. Arguello, S. P. Chockalingam, E. P. Rosenthal, L. Zhao, C. Gutiérrez, J. H. Kang, W. C. Chung, R. M. Fernandes, S. Jia, A. J. Millis, R. J. Cava, and A. N. Pasupathy, Physical Review B89, 235115 (2014)

    2014

  26. [26]

    Weber, R

    F. Weber, R. Hott, R. Heid, K.-P. Bohnen, S. Rosenkranz, J.-P. Castellan, R. Osborn, A. H. Said, B. M. Leu, and D. Reznik, Physical Review B87, 245111 (2013)

    2013

  27. [27]

    M. D. Johannes and I. I. Mazin, Physical Review B77, 165135 (2008)

    2008