pith. machine review for the scientific record. sign in

arxiv: 2605.13458 · v1 · submitted 2026-05-13 · ❄️ cond-mat.mtrl-sci

Recognition: 2 theorem links

· Lean Theorem

Multiple Softening Q-vectors Driving a Cascade of CDW Phases in 1T-VSe₂

Zheng-Hong Li , Yung-Ting Lee , Yu-Chan Tai , Cheng-Tien Chiang , Chien-Cheng Kuo , Meng-Kai Lin , Chun-Liang Lin , Hung-Chung Hsueh
show 3 more authors
Ming-Chiang Chung Po-Tuan Chen Chi-Cheng Lee
Authors on Pith no claims yet

Pith reviewed 2026-05-14 18:29 UTC · model grok-4.3

classification ❄️ cond-mat.mtrl-sci
keywords charge density wave1T-VSe2phonon instabilitystructural relaxationmonolayer materialfirst-principles calculationphase transition
0
0 comments X

The pith

Multiple phonon instabilities at distinct Q-vectors in monolayer 1T-VSe2 drive separate CDW distortions that all converge on the same low-energy 2√3×4 superstructure.

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

The paper tracks the structural evolution of monolayer 1T-VSe2 by computing phonon spectra at successive stages of distortion. Several imaginary-frequency modes in the pristine lattice generate first-generation CDW phases. These intermediate structures remain unstable and undergo further phonon-driven relaxations that lower symmetry and enlarge the unit cell. All identified pathways, despite starting from different initial wave vectors, terminate at nearly degenerate 2√3×4 configurations. This pattern indicates that a single stable endpoint can be reached through multiple phonon-driven routes.

Core claim

Through iterative phonon-driven relaxations, multiple transformation pathways that originate from distinct imaginary-frequency modes at different Q_CDW vectors converge toward the same low-energy 2√3×4 CDW configuration. Although the pathways begin from separate intermediate CDW states, they reach nearly degenerate, energetically stable phases, demonstrating that different phonon-driven routes can lead to the same ground-state configuration.

What carries the argument

Iterative phonon analysis followed by structural relaxation, starting from imaginary-frequency modes at various Q_CDW vectors and mapping each step to larger superstructures.

If this is right

  • Distinct initial CDW phases can relax into one common ground-state structure.
  • The 2√3×4 configuration is the stable endpoint of the cascade.
  • Phonon instabilities at multiple wave vectors produce a hierarchical sequence of ordered phases.
  • Different transformation routes yield nearly degenerate energies for the final state.

Where Pith is reading between the lines

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

  • Analogous multi-step phonon cascades may govern CDW ordering in other layered transition-metal compounds where several competing wave vectors are present.
  • Experiments that resolve the sequence of lattice distortions as temperature is lowered could test whether all pathways are active.
  • Computational searches for CDW ground states in two-dimensional materials should include successive phonon relaxations to reduce the risk of reporting metastable intermediates.

Load-bearing premise

Standard first-principles phonon calculations capture every relevant lattice instability and the iterative relaxation procedure reaches the global energy minimum rather than a local one.

What would settle it

Identification of an energetically lower CDW superstructure in 1T-VSe2 whose atomic arrangement cannot be obtained by relaxing any of the reported intermediate phases would show that the convergence is incomplete.

Figures

Figures reproduced from arXiv: 2605.13458 by Cheng-Tien Chiang, Chi-Cheng Lee, Chien-Cheng Kuo, Chun-Liang Lin, Hung-Chung Hsueh, Meng-Kai Lin, Ming-Chiang Chung, Po-Tuan Chen, Yu-Chan Tai, Yung-Ting Lee, Zheng-Hong Li.

Figure 1
Figure 1. Figure 1: Schematic illustration of the phonon-driven cascade. The top panels show the phonon dispersions of the normal phase, intermediate CDW phase, and [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
read the original abstract

Charge density wave (CDW) formation in two-dimensional materials is governed by complex competing lattice instabilities that remain incompletely understood. Here, we investigate the structural evolution of monolayer $\mathrm{1T-VSe}_{2}$ using first-principles electronic and phonon calculations. The pristine phase exhibits several imaginary-frequency phonon modes associated with dominant instability wave vectors $\mathrm{Q}_{CDW}$, which generate the first-generation CDW phases. Subsequent phonon analyses reveal that several of these intermediate structures remain dynamically unstable and undergo further symmetry-lowering distortions into larger superstructures. Through iterative phonon-driven relaxations, we identify multiple transformation pathways that converge toward the same low-energy $2\sqrt{3}\times4$ CDW configuration. Although these pathways originate from distinct intermediate CDW states, they ultimately reach nearly degenerate energetically stable phases, demonstrating that different phonon-driven routes can lead to the same ground-state configuration. The results establish a unified phonon-driven cascade mechanism for hierarchical CDW formation in monolayer $\mathrm{1T-VSe}_{2}$ and provide a systematic framework for understanding competing ordered phases in low-dimensional quantum materials.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript uses first-principles electronic and phonon calculations to show that monolayer 1T-VSe₂ exhibits multiple imaginary-frequency modes at distinct Q_CDW vectors in the pristine phase. These drive first-generation CDW structures whose subsequent phonon instabilities trigger further symmetry-lowering distortions; iterative relaxations from several distinct intermediate states are reported to converge to the same low-energy 2√3×4 CDW configuration, establishing a unified phonon-driven cascade mechanism.

Significance. If the reported convergence is robust, the work supplies a concrete, phonon-based account of hierarchical CDW formation in a 2D material and demonstrates that multiple transformation pathways can terminate at the same ground-state superstructure. The direct use of successive phonon calculations to map competing instabilities is a methodological strength that could be extended to other low-dimensional systems.

major comments (1)
  1. [Results section on iterative relaxations and phonon analyses] The central claim that the iterative phonon-driven relaxations reach the global 2√3×4 minimum (rather than a local minimum) is load-bearing yet insufficiently demonstrated. Standard DFT phonon calculations are local by construction; without an independent global search, exhaustive enumeration of supercell distortions, or explicit checks for anharmonic or long-range instabilities, the possibility remains that lower-energy configurations exist outside the sampled pathways.
minor comments (2)
  1. [Abstract] The abstract states that the pathways 'converge toward the same low-energy 2√3×4 CDW configuration' but does not quantify the energy differences or degeneracy among the final structures; adding these values would strengthen the claim of near-degeneracy.
  2. [Introduction and Methods] Notation for the wave vectors (Q_CDW) and the final superstructure (2√3×4) is introduced without an explicit definition or reference to the Brillouin-zone path used; a short clarification or supplementary figure would improve readability.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the constructive report and positive assessment of the significance of our work. We address the single major comment below, clarifying the scope of our claims while strengthening the presentation of our results.

read point-by-point responses

  1. Referee: [Results section on iterative relaxations and phonon analyses] The central claim that the iterative phonon-driven relaxations reach the global 2√3×4 minimum (rather than a local minimum) is load-bearing yet insufficiently demonstrated. Standard DFT phonon calculations are local by construction; without an independent global search, exhaustive enumeration of supercell distortions, or explicit checks for anharmonic or long-range instabilities, the possibility remains that lower-energy configurations exist outside the sampled pathways.

    Authors: We agree that phonon-based relaxations are inherently local and that our study does not constitute an exhaustive global search. Our central result is the observation that multiple, independent transformation pathways—each initiated from a distinct imaginary phonon mode and first-generation CDW intermediate—converge to the identical 2√3×4 superstructure, which is dynamically stable (no imaginary frequencies). We have revised the manuscript to (i) explicitly state that we do not claim an absolute global minimum via exhaustive enumeration, (ii) add a dedicated paragraph discussing the local nature of the method and the computational impracticality of anharmonic or exhaustive supercell searches for this system, and (iii) include a direct energy comparison of the 2√3×4 phase against other CDW supercells reported in the VSe₂ literature. These additions make the evidential basis and its limitations transparent while preserving the demonstration that the phonon cascade robustly selects this structure from the sampled instabilities. revision: partial

Circularity Check

0 steps flagged

No circularity: direct first-principles phonon relaxations

full rationale

The paper performs standard DFT phonon calculations on the pristine structure to locate imaginary modes at specific Q-vectors, then relaxes the resulting distorted supercells and repeats the phonon analysis on the intermediates. This is an explicit computational search procedure whose outputs (the 2√3×4 ground-state configuration and the listed pathways) are generated by the simulation engine rather than being presupposed by any equation or fitted parameter. No self-definitional loops, no fitted inputs renamed as predictions, and no load-bearing self-citations or imported uniqueness theorems appear in the described workflow. The iterative cascade is therefore self-contained against external benchmarks and does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain assumption that density-functional-theory phonon calculations reliably identify the dominant instabilities and that iterative structural relaxations adequately sample the configuration space to locate the ground state.

axioms (1)
  • domain assumption Density functional theory approximations are sufficient to capture the phonon instabilities in 1T-VSe2
    Standard assumption invoked for all first-principles phonon work on this class of materials.

pith-pipeline@v0.9.0 · 5549 in / 1190 out tokens · 47151 ms · 2026-05-14T18:29:41.138482+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages

  1. [1]

    The Dynamics of Charge-Density Waves

    Grüner, G. The Dynamics of Charge-Density Waves. Rev. Mod. Phys. 1988, 60, 1129–1181. https://doi.org/10.1103/RevModPhys.60.1129

    work page doi:10.1103/revmodphys.60.1129 1988

  2. [2]

    Electronic Crystals: An Experimental Overview

    Monceau, P. Electronic Crystals: An Experimental Overview. Adv. Phys. 2012, 61 (4), 325– 581. https://doi.org/1.1080/00018732.2012.719674

    work page arXiv 2012

  3. [3]

    Doping Tunable CDW Phase Transition in Bulk 1T-ZrSe_2

    Ø rsted, A.; Scarfato, A.; Barreteau, C.; Giannini, E.; Renner, C. Doping Tunable CDW Phase Transition in Bulk 1T-ZrSe_2. Nano Lett. 2025, 25 (4), 1729–1735. https://doi.org/10.1021/acs.nanolett.4c06377

    work page doi:10.1021/acs.nanolett.4c06377 2025

  4. [4]

    Li, L.; O’Farrell, E. C. T.; Loh, K. P.; Eda, G.; Özyilmaz, B.; Castro Neto, A. H. Controlling Many-Body States by the Electric -Field Effect in a Two -Dimensional Material. Nature 2016, 529, 185–189. https://doi.org/10.1038/nature16175

    work page doi:10.1038/nature16175 2016

  5. [5]

    Functions of Bounded Variation and Free Discontinuity Problems

    Peierls, R. E., Quantum Theory of Solids, Oxford Classic Texts in the Physical Sciences. https://doi.org/10.1093/acprof:oso/9780198507819.001.0001

    work page doi:10.1093/acprof:oso/9780198507819.001.0001

  6. [6]

    W.; Guo, J

    Zhu, X.; Cao, Y.; Zhang, J.; Plummer, E. W.; Guo, J. Classification of Charge Density Waves Based on Their Nature. Proc. Natl. Acad. Sci. U.S.A. 2015, 112 (8), 2367–2371. https://doi.org/10.1073/pnas.1424791112

    work page doi:10.1073/pnas.1424791112 2015

  7. [7]

    Revealing the Peierls Instability Caused Charge Density Wave in Two-Dimensional NbSe2

    Lee, Y.-T.; Chen, P.-T.; Li, Z.-H.; Wu, J.-Y.; Kuo, C.-N.; Lue, C.-S.; Wu, C.-T.; Kuo, C.- C.; Chiang, C.-T.; Ozaki, T.; Lin, C.-L.; Lee, C.-C.; Hsueh, H.-C.; Chung, M.-C. Revealing the Peierls Instability Caused Charge Density Wave in Two-Dimensional NbSe2. ACS Mater. Lett. 2024, 6 (7), 2941–2947. https://doi.org/10.1021/acsmaterialslett.4c00142

    work page doi:10.1021/acsmaterialslett.4c00142 2024

  8. [8]

    M.; He, Y.; van Wezel, J.; Rahn, D

    Soumyanarayanan, A.; Yee, M. M.; He, Y.; van Wezel, J.; Rahn, D. J.; Rossnagel, K.; Hudson, E. W.; Norman, M. R.; Hoffman, J. E. Quantum Phase Transition from Triangular to Stripe Charge Order in NbSe2. Proc. Natl. Acad. Sci. U.S.A. 2013, 110, 1623–1627. https://doi.org/10.1073/pnas.1211387110

    work page doi:10.1073/pnas.1211387110 2013

  9. [9]

    Huang, Y.-T.; Li, Z.-Z. et al. Complex charge density waves in simple electronic systems of two-dimensional III2–VI3 materials. Nature Communications 2024, 15, 9983. https://doi.org/10.1038/s41467-024-54205-1

    work page doi:10.1038/s41467-024-54205-1 2024

  10. [10]

    Wang, Q. et al. Unraveling competing charge orders in monolayer VTe 2. Phys. Rev. B 2025, 112, 245410. https://doi.org/10.1103/ggfw-ypkp

    work page doi:10.1103/ggfw-ypkp 2025

  11. [11]

    Chen, P. et al. , Unique Gap Structure and Symmetry of the Charge Density Wave in Single-Layer VSe2. Phys. Rev. Lett. 2018, 121, 196402. https://doi.org/10.1103/PhysRevLett.121.196402

    work page doi:10.1103/physrevlett.121.196402 2018

  12. [12]

    et al., Dimensional crossover and symmetry transformation of charge density waves in VSe2

    Chen, P. et al., Dimensional crossover and symmetry transformation of charge density waves in VSe2. Phys. Rev. B 2022, 105, L161404. https://doi.org/10.1103/PhysRevB.105.L161404

    work page doi:10.1103/physrevb.105.l161404 2022

  13. [13]

    et al., Theoretical study of pseudo two-step metal–insulator transition in monolayer 1T-VSe2

    Du, H. et al., Theoretical study of pseudo two-step metal–insulator transition in monolayer 1T-VSe2. J. Appl. Phys. 2025, 137, 174305. https://doi.org/10.1063/5.0256435

    work page doi:10.1063/5.0256435 2025

  14. [14]

    Duvjir, G.; Choi, B. K. et al. Emergence of a Metal–Insulator Transition and High- Temperature Charge-Density Waves in VSe2 at the Monolayer Limit. Nano Lett. 2018, 18, 9, 5432–5438. https://doi.org/10.1021/acs.nanolett.8b01764

    work page doi:10.1021/acs.nanolett.8b01764 2018

  15. [16]

    Ozaki \ and\ author H

    Ozaki, T.; Kino, H. Numerical atomic basis orbitals from H to Kr. Phys. Rev. B 2004, 69, 195113. https://doi.org/10.1103/PhysRevB.69.195113

    work page doi:10.1103/physrevb.69.195113 2004

  16. [17]

    Generalized gradient approximation made simple

    Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77 (18), 3865–3868. https://doi.org/10.1103/PhysRevLett.77.3865

    work page doi:10.1103/physrevlett.77.3865 1996

  17. [18]

    Efficient Projector Expansion for the ab initio LCAO Method

    Ozaki, T.; Kino, H. Efficient Projector Expansion for the ab initio LCAO Method. Phys. Rev. B 2005, 72 (4), 045121. https://doi.org/10.1103/PhysRevB.72.045121

    work page doi:10.1103/physrevb.72.045121 2005

  18. [19]

    Ozaki ,\ https://doi.org/10.1103/PhysRevB.67.155108 journal journal Phys

    Ozaki, T. Variationally optimized atomic orbitals for large-scale electronic structures. Phys. Rev. B 2003, 67, 155108. https://journals.aps.org/prb/abstract/10.1103/PhysRevB.67.155108

    work page doi:10.1103/physrevb.67.155108 2003

  19. [20]

    Physical Review B , author =

    Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Phys. Rev. B 1976, 13, 5188–5192. https://doi.org/10.1103/PhysRevB.13.5188

    work page doi:10.1103/physrevb.13.5188 1976

  20. [21]

    Momma, K. et al. VESTA3 for Three-Dimensional Visualization of Crystal, Volumetric and Morphology Data. J. Appl. Crystallogr. 2011, 44, 1272–1276. https://doi.org/10.1107/S0021889811038970

    work page doi:10.1107/s0021889811038970 2011

  21. [22]

    Lee, Y.-T. et al. OpenMX Viewer: A web-based crystalline and molecular graphical user interface program. J. Mol. Graph. Model. 2019, 89, 192. https://doi.org/10.1016/j.jmgm.2019.03.013

    work page doi:10.1016/j.jmgm.2019.03.013 2019

  22. [23]

    Kresse, G. et al. Ab initio Force Constant Approach to Phonon Dispersion Relations of Diamond and Graphite. EPL. 1995, 32, 729. https://doi.org/10.1209/0295-5075/32/9/005

    work page doi:10.1209/0295-5075/32/9/005 1995

  23. [24]

    Parlinski, K. et al. First-principles determination of the soft mode in cubic ZrO2 V. Phys. Rev. Lett. 1997, 78, 4063. https://doi.org/10.1103/PhysRevLett.78.4063

    work page doi:10.1103/physrevlett.78.4063 1997

  24. [25]

    Pallikara, I. et al. The physical significance of imaginary phonon modes in crystals Electron. Struct. 2022, 4, 033002. https://doi.org/10.1088/2516-1075/ac78b3

    work page doi:10.1088/2516-1075/ac78b3 2022

  25. [26]

    Fumega, A. O. et al. Anharmonicity Reveals the Tunability of the Charge Density Wave Orders in Monolayer VSe2. Nano Lett. 2023, 23, 5, 1794–1800. https://pubs.acs.org/doi/10.1021/acs.nanolett.2c04584

    work page doi:10.1021/acs.nanolett.2c04584 2023

  26. [28]

    et al., Multiple charge density wave phases of monolayer VSe 2 manifested by graphene substrates

    Duvjir, G. et al., Multiple charge density wave phases of monolayer VSe 2 manifested by graphene substrates. Nanotechnology 2021, 32, 364002. https://doi.org/10.1088/1361-6528/ac06f3

    work page doi:10.1088/1361-6528/ac06f3 2021

  27. [29]

    et al., Charge Density Wave State Suppresses Ferromagnetic Ordering in VSe2 Monolayers

    Coelho, P.M. et al., Charge Density Wave State Suppresses Ferromagnetic Ordering in VSe2 Monolayers. J. Phys. Chem. C 2019, 123, 14089. https://doi.org/10.1021/acs.jpcc.9b04281

    work page doi:10.1021/acs.jpcc.9b04281 2019

  28. [30]

    Zhang, D.; Ha, J. et al. Strain engineering a 4𝑎×√3𝑎 charge-density-wave phase in transition-metal dichalcogenide 1T−VSe2. Phys. Rev. Materials 2017, 1, 024005. https://doi.org/10.1103/PhysRevMaterials.1.024005

    work page doi:10.1103/physrevmaterials.1.024005 2017

  29. [31]

    Chen, G.; Howard, S. T. et al. Correlating structural, electronic, and magnetic properties of epitaxial VSe2 thin films. Phys. Rev. B 2020, 102, 115149. https://doi.org/10.1103/PhysRevB.102.115149

    work page doi:10.1103/physrevb.102.115149 2020

  30. [32]

    Ma, W.; Yang, Y. et al. Nanoscale Determination of the Metal–Insulator Transition in Intercalated Bulk VSe2. Nano Lett. 2025, 25, 31, 11852–11859. https://doi.org/10.1021/acs.nanolett.5c02203

    work page doi:10.1021/acs.nanolett.5c02203 2025