{"id":"c859ac7b-7f86-4e3a-802c-9161e8d7e85e","arxiv_id":"1908.03655","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Stripe and quasi-stripe CDW phases in monolayer MX2 compounds emerge as local minima of an isotropic McMillan-Nakanishi-Shiba free energy once the 120-degree wave-vector constraint is relaxed.","lead":"A numerical study of a classic Ginzburg-Landau model finds many unexpected local minima in the free energy of charge-density-wave materials, including states with striped and quasi-striped patterns. The result suggests that anisotropic stripe patterns in monolayer tantalum dichalcogenides can appear without invoking anisotropic interactions.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The T/stripe claims rest on minima in a 2D slice with experimental angles fixed; full 4D stability is unverified.","rationale":"The reader's verdict of CONDITIONAL is appropriate, and the reader's rationale already notes that fixing the T-phase angles lowers circularity. However, the reader's stated weakest assumption is transferability of bulk parameters to monolayers, whereas I regard the constrained-slice issue as more directly load-bearing: even with perfect monolayer parameters, the reported T/stripe 'local minima' are not shown to be stationary in the full Q-space. The two-parameter-set comparison and reproduction of Nakanishi-Shiba results in the 120°-constrained case are genuine supporting evidence, but they do not cover the unconstrained anisotropic search, which is performed with only one parameter set and N = 1. The concern is therefore not that the theory is wrong, but that the central claim has not yet been demonstrated as stated. A full-space minimization is a straightforward computational check and would settle the issue.","tokens_in":11742,"tokens_out":7053,"duration_ms":83597,"concrete_test":"Recompute the T and stripe searches of Section 1.3 by minimizing F over all four independent components of (Q(1), Q(2)) with Q(3) = −Q(1) − Q(2), without fixing δφ2, δφ3, and without imposing Q(1) = Q_C(1). Initialize from the Figure 4 minima and from random seeds, and compute the 4x4 Hessian at the reported points. If the relaxed stationary points move away from the experimental angles or have negative Hessian eigenvalues, the claimed anisotropic minima are artifacts of the imposed constraints. Repeating the same unconstrained search with N = 2 and N = 3 would additionally test whether the low-order truncation used in Figure 4 preserves the minima.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that the isotropic McMillan-Nakanishi-Shiba free energy has local minima corresponding to the T and stripe phases without imposing anisotropic interactions. The search in Section 1.3 does not establish this. For the T phase, the authors state: 'From the triple-Q condition there are four independent degrees of freedom. Here, we fix two degrees of freedom, namely, the angles between domain walls, to visualize the free energy. The angles between domain walls are known from experiments [14]', fixing δφ2 = 180° − φC and δφ3 = 150°. For the stripe phase, they impose Q(1) = Q_C(1) based on reference [11]. The minima shown in Figure 4 are therefore minima of the free energy restricted to the remaining 2D slice of Q-space. A point that minimizes F on a codimension-2 surface need not be stationary in the full 4D space; if the gradient along the two fixed directions is nonzero, these are not free-energy local minima of the unconstrained problem. Since the abstract claims the phases appear 'naturally' as local minima of the isotropic free energy, the experimentally known angles cannot be inputs at this stage. The paper reports no full 4D minimization, no Hessian check, and no unconstrained search over all four independent components of the wavevectors. This is the most load-bearing gap because the key conclusion depends on these points being genuine minima, not artifacts of the imposed experimental constraints.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper claims to explain the origin of stripe and quasi-stripe (T) CDW domain-wall phases in monolayer 1T-TaS2 and 2H-TaSe2 within the isotropic McMillan-Nakanishi-Shiba Ginzburg-Landau free energy, without invoking anisotropic Coulomb interactions or interlayer coupling. The authors first map the free energy over wave-vector space under a 120° constraint, finding a multivalley landscape with many local minima. Removing the 120° constraint, they report free-energy minima whose wave-vector triples form regular triangles and identify these with the experimentally observed T and stripe phases. They conclude that anisotropic domain walls can appear naturally from isotropic interactions, with Coulomb and interlayer effects demoted to secondary factors.","tokens_in":12083,"tokens_out":5037,"duration_ms":50335,"significance":"If fully established, the claim would be significant: it would provide an alternative explanation for stripe and T phases in monolayers, predict new CDW phases (e.g., √7×√7 TaSe2 with T domain walls), and reframe prior Coulomb/interlayer explanations as secondary. The paper's strength is the extension of Nakanishi-Shiba to a full two-dimensional wave-vector landscape, with two parameter sets giving qualitatively similar multivalley structures. The central difficulty is that the key unconstrained claim is not actually demonstrated by the reported constrained minimizations, so the significance is currently conditional on additional numerical evidence.","major_comments":[{"comment":"The claimed T-phase minima are found on a codimension-2 slice of Q-space. The text states that 'from the triple-Q condition there are four independent degrees of freedom. Here, we fix two degrees of freedom, namely, the angles between domain walls, to visualize the free energy. The angles between domain walls are known from experiments [14].' Fixing δφ2 and δφ3 to the experimental values means the displayed minima are not necessarily stationary points of the full free energy: the gradient along the two fixed directions is not checked. Since the abstract claims the stripe and T phases appear 'naturally' as local minima of the isotropic free energy, the experimentally known angles cannot be used as inputs at this stage. The authors should perform an unconstrained minimization over all four independent wave-vector components (with a Hessian verification) and report whether the previously found points are genuine stationary points. The same concern applies to the stripe phase, where Q(1)=Q_C(1) is imposed from ref [11] without checking full 4D stationarity.","section":"§1.3, Figure 4(a)"},{"comment":"The agreement with experiment for the T phase is partly circular. Because the domain-wall angles δφ2 and δφ3 are fixed to the values from ref [14] before minimization, the observed agreement in those angles is built in; the only nontrivial comparison is the domain size, and the paper states that 'the domain size in our calculation is smaller than bulk crystal.' The central claim that the T phase appears naturally from isotropic interactions would require the minimization to reproduce the experimental angles without imposing them.","section":"§1.3, Figure 4(a)"},{"comment":"The free-energy coefficients a0, a1, b, c, d and the gradient operator e(-i∇) are taken from fits to bulk 1T-TaS2 and 2H-TaSe2 (refs [30,31,32]) and used without modification for free-standing monolayers. The paper provides no monolayer-specific validation, sensitivity analysis, or error estimates for these parameters. Because the multivalley landscape and the positions of the T/stripe minima are quantitative outputs, the possibility that these minima are artifacts of bulk-derived parameters is a load-bearing uncertainty. The authors should either provide monolayer-specific parameter estimates, show that the qualitative conclusions are robust to parameter variations, or clearly state this as a limitation.","section":"§1.1 and Methods"},{"comment":"Numerical results are reported without convergence analysis, error bars, or a description of the numerical tolerances used to identify local minima. The claim that type-1 and type-2 free energies give 'almost identical local minima' is not quantified, and no code or data repository is provided. To make the multivalley landscape reproducible, the authors should specify the grid, the convergence criteria for solving ∂F/∂Δ=0, and the criterion for when a point is called a local minimum.","section":"Numerical results, Figures 2-4"}],"minor_comments":[{"comment":"Reference [33] is listed as 'To be published'; this is not an accessible citation and should be replaced or removed.","section":"References"},{"comment":"The statement that 'the only input from experimental data is the incommensurate wave vectors Q_IC' is contradicted by §1.3, where experimental domain-wall angles from ref [14] are used as inputs; the text should be revised to acknowledge this.","section":"Introduction and §1.3"},{"comment":"The distinction between type-1 and type-2 free energies is described only as 'different forms of e(-i∇)'; the explicit forms or equations should be given.","section":"Methods"},{"comment":"Figure 4 shows free-energy landscapes without contour levels or a color scale, making it difficult to judge the depth and sharpness of the claimed minima; adding contours and marking the constrained slice would improve clarity.","section":"Figure 4"},{"comment":"The caption states 'in this article we focus on |q(1)|=|q(2)|=|q(3)| and each |q(i)| are separated by a 120°,' which is inconsistent with §1.3 where the 120° constraint is removed; the scope should be clarified.","section":"Figure 6 caption"},{"comment":"The entropy/Kosterlitz-Thouless analogy in the Discussion is qualitative and not developed; no entropy contribution is computed. The paragraph should be framed as speculation or supported by a quantitative estimate.","section":"Discussion"}],"recommendation":"major_revision","confidential_remarks":"The paper addresses a significant question and the multivalley idea is attractive. The main technical gap—constrained versus full minimization—is fixable in principle but requires substantial new numerical work, so major revision is appropriate. I would also urge the editor to require a code/data availability statement, given that the numerical results are central to the claims."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Good to have something concrete to argue with. The genuinely new piece is the exhaustive 2D sweep of the McMillan-Nakanishi-Shiba free energy for 1T-TaS2 and 2H-TaSe2, both under and beyond the 120-degree Q-vector constraint. The multivalley landscape with new local minima near Q_C is a visible extension of Nakanishi-Shiba's 1D paths, and the fact that type-1 and type-2 parameter sets give nearly identical minima is real supporting evidence. Their reading of the stripe and T phases as natural consequences of commensurability and higher harmonics, without invoking Coulomb or inter-layer interactions, is a clean idea worth taking seriously.\n\nThe soft spots are mostly about how much the numerics actually prove. The T-phase section fixes two of the four wave-vector degrees of freedom to the experimental domain-wall angles from ref [14], then reports minima in the remaining 2D plane. That is not a demonstration that those points are stationary in the full 4D Q-space; a point can minimize F on a codimension-2 slice and still have nonzero gradient along the fixed directions. Same for the stripe phase, where Q(1) is pinned to Q_C(1) on experimental grounds. Repeating the search without pre-fixing the angles is not a cosmetic improvement; it is the difference between finding local minima of F and finding minima of a restricted function. The paper claims the phases 'appear naturally,' so the absence of the unconstrained calculation is the load-bearing gap. Minor but worth saying: no code, no convergence analysis, no Hessian check, and the coefficients are taken from bulk fits with no monolayer-specific validation. The authors are upfront that they used ref [30] parameters, so it's not hidden, but it still limits quantitative force. The experimental comparison in Fig 4 is only about domain size, and they acknowledge the size mismatch, so the agreement is qualitative.\n\nI also want to note what is not wrong. The citation pattern is reasonable: they are explicit about what Nakanishi-Shiba did and what they extend. The claim about entropy and domain-wall formation is speculative but labeled as such. I don't think the paper is sloppy; it just stops short of the calculation that would close the argument.\n\nFor a referee: I would send it out, but with a request for either a full 4D unconstrained minimization on a reasonable grid or a Hessian check at the reported minima, plus code and parameter disclosure. If those come back clean, the paper would be a solid contribution; right now it's a strong proposal with a missing eigenvalue.","headline":"Worth a serious referee, but the central 'natural minima' claim needs an unconstrained 4D check before I believe it.","tokens_in":12562,"tokens_out":2049,"would_cite":false,"duration_ms":22417,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Stripe and quasi-stripe CDW phases in monolayer MX2 compounds can arise from an isotropic free-energy landscape alone, with no anisotropic interactions assumed.","keywords":["charge density wave","stripe phase","quasi-stripe phase","McMillan-Nakanishi-Shiba model","multivalley free energy landscape","1T-TaS2","2H-TaSe2","topological defects"],"falsifier":"Compute the same free-energy landscape with parameters obtained directly from monolayer 1T-TaS2 or 2H-TaSe2, for example from first-principles phonon calculations or measured monolayer CDW wave vectors, and check whether the regular-triangle local minima for the stripe and T phases survive when the 120° constraint is removed. A free-standing monolayer 1T-TaS2 that shows no quasi-stripe phase on heating would also contradict the central prediction.","tokens_in":11543,"feed_emoji":"⚡","tokens_out":9173,"duration_ms":87329,"temperature":0.7,"pith_summary":"This paper aims to explain why stripe-shaped and quasi-stripe (triclinic) charge-density-wave (CDW) patterns appear in monolayers of the transition-metal dichalcogenides 1T-TaS2 and 2H-TaSe2. It argues that these anisotropic domain-wall structures do not require anisotropic interactions, such as Coulomb repulsion between domain walls or inter-layer stacking. Instead, they emerge as local minima of an isotropic Ginzburg-Landau free-energy landscape, the McMillan-Nakanishi-Shiba model, once the search over CDW wave vectors is broadened beyond the usual symmetry-constrained paths. If correct, the long-standing puzzle of anisotropic CDW phases in ultrathin crystals becomes a problem of topological defect formation and harmonic interference, and the model predicts new CDW phases that experiments could look for.","feed_headline":"Stripe CDW phases arise without anisotropic forces","feed_subtitle":"A multivalley energy landscape in monolayer TaS2 and TaSe2 creates stripe and quasi-stripe order by itself.","key_machinery":"The machinery is the McMillan-Nakanishi-Shiba free energy, a Ginzburg-Landau functional with three complex CDW order parameters $\\psi_i(\\mathbf{r})$ whose wave vectors $\\mathbf{Q}^{(i)}$ obey the triple-Q condition $\\mathbf{Q}^{(1)}+\\mathbf{Q}^{(2)}+\\mathbf{Q}^{(3)}=0$. The paper searches this free energy over the full two-dimensional $\\mathbf{Q}^{(1)}$ plane, including higher harmonics of the order parameter (indexed by $N$), which are essential for creating the multivalley structure: with $N=0$ the landscape is flat, while $N=1,2,3$ develop many local minima. After removing the 120° separation constraint, the local minima of the three $\\mathbf{Q}^{(i)}$ form regular triangles, and that geometric condition selects anisotropic domain-wall phases (stripe vs. quasi-stripe) depending on the commensurability indices $(\\mu,\\nu)$ of the material.","core_discovery":"The central claim is that anisotropic stripe and quasi-stripe CDW domain walls appear naturally, not accidentally: they are local minima of the isotropic McMillan-Nakanishi-Shiba free energy for monolayer 1T-TaS2 and 2H-TaSe2. For CDW wave vectors separated by 120°, the free-energy landscape already contains many metastable minima beyond the known incommensurate, nearly commensurate, and commensurate states. When the 120° constraint is removed, triplets of these minima form regular triangles in wave-vector space; the corresponding real-space domain walls have the angles and quasi-stripe character of the experimentally observed T phase in 1T-TaS2 and stripe phase in 2H-TaSe2. The paper therefore concludes that Coulomb domain-wall repulsion and inter-layer coupling are secondary factors that adjust domain sizes, not the origin of the anisotropic phase, and that the same mechanism predicts additional CDW states.","pith_inferences":["Beyond the paper: the same mechanism should operate in any van der Waals material with three coexisting CDW wave vectors satisfying the triple-Q condition, so the model's predictions could be tested in other MX2 compounds with different $(\\mu,\\nu)$.","Beyond the paper: if hidden CDW states are metastable minima of this landscape, then the path taken in wave-vector space during a quench should determine which minimum is reached; momentum-resolved pump-probe experiments could test this by mapping the transient wave vectors.","Beyond the paper: the analogy with Abrikosov vortex formation suggests a quantitative prediction, namely that the density of domain walls should grow continuously with temperature in a single monolayer, and the heating/cooling asymmetry of domain size should be measurable in transport or scanning tunneling microscopy."],"forward_implications":["Monolayer 1T-TaS2 and 2H-TaSe2 should show stripe or quasi-stripe CDW phases even in the complete absence of inter-layer stacking, consistent with observed monolayer CDW textures.","Coulomb interaction between domain walls and inter-layer interaction are demoted to secondary factors that change domain size but are not needed to create anisotropic domain walls.","The domain-wall angle is set by the commensurability indices $(\\mu,\\nu)$ through $\\phi_C = \\cos^{-1}((\\mu+\\nu/2)/\\sqrt{\\mu^2+\\mu\\nu+\\nu^2})$: both $\\mu$ and $\\nu$ nonzero gives quasi-stripe T walls, while $\\nu=0$ gives stripe walls.","The multivalley landscape predicts new CDW phases and offers a route to explaining 'hidden' CDW states as metastable local minima.","On heating, anisotropic domain walls can form one by one, breaking the threefold rotational symmetry and increasing entropy, which explains why stripe and T phases appear only on heating from the commensurate state."],"supporting_citations":[{"why":"supplies the original McMillan free energy with triple-Q CDW order parameters and the discommensuration concept the paper extends.","marker":"[2]"},{"why":"supplies the 1T-TaS2 free-energy parameters and the earlier domain-like NC-phase analysis that the paper generalizes to a two-dimensional wave-vector search.","marker":"[30]"},{"why":"supplies the 2H-TaSe2 free-energy parameters and the IC-C transition analysis that the paper generalizes.","marker":"[31]"},{"why":"provides the ring-like diffuse-scattering data used to fix the type-1 free-energy coefficients for 1T-TaS2.","marker":"[32]"},{"why":"is the inter-layer stacking explanation of stripe and T phases that the paper argues is not necessary in monolayers.","marker":"[4]"},{"why":"provides the experimental T-phase domain-wall angles and domain sizes in 1T-TaS2 used to compare with the calculated quasi-stripe structure.","marker":"[14]"},{"why":"provides the transmission-electron-microscopy observation of multiple T domain walls in a free-standing monolayer 1T-TaS2 that motivates the monolayer calculation.","marker":"[22]"},{"why":"provides the experimental broken-hexagonal stripe structure of 2H-TaSe2 used as the stripe-phase comparison.","marker":"[11]"}],"fun_headline_variants":["Stripe CDWs arise from free energy minima, not extra forces","Anisotropic CDW phases need no anisotropic interactions","Free energy landscape spawns stripe and quasi-stripe CDW","Stripe phases in monolayer TMDs are ground-state minima, not accidents","Multivalley free energy explains stripe CDW without anisotropic forces"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing assumption is that the free-energy coefficients fitted to bulk 1T-TaS2 and 2H-TaSe2 transfer unchanged to a single free-standing monolayer; the paper uses the bulk parameters without monolayer-specific validation and relies on them in Section 1.1 and the Methods.","fun_headline_variants_meta":{"raw":{"variants":["Stripe CDWs arise from free energy minima, not extra forces","Anisotropic CDW phases need no anisotropic interactions","Free energy landscape spawns stripe and quasi-stripe CDW","Stripe phases in monolayer TMDs are ground-state minima, not accidents","Multivalley free energy explains stripe CDW without anisotropic forces"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000807,"raw_usage":{"total_tokens":3621,"prompt_tokens":1104,"completion_tokens":2517,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":720,"completion_tokens_details":{"reasoning_tokens":2429}},"tokens_in":720,"tokens_out":2517,"duration_ms":17820,"temperature":1.0,"reasoning_tokens":2429,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:06:39.554669+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the same free-energy landscape with parameters obtained directly from monolayer 1T-TaS2 or 2H-TaSe2, for example from first-principles phonon calculations or measured monolayer CDW wave vectors, and check whether the regular-triangle local minima for the stripe and T phases survive when the 120° constraint is removed. A free-standing monolayer 1T-TaS2 that shows no quasi-stripe phase on heating would also contradict the central prediction.","supporting_citations":[{"cited_title":"Theory of discommensurations and the commensurate-incommensurate charge- density-wave phase transition,","cited_arxiv_id":null,"evidence_quote":"supplies the original McMillan free energy with triple-Q CDW order parameters and the discommensuration concept the paper extends."},{"cited_title":"Domain-like Incommensurate Charge-Density-Wave States and the First- Order Incommensurate-Commensurate Transitions in Layered Tantalum Dichalcogenides. I. 1T- Polytype,","cited_arxiv_id":null,"evidence_quote":"supplies the 1T-TaS2 free-energy parameters and the earlier domain-like NC-phase analysis that the paper generalizes to a two-dimensional wave-vector search."},{"cited_title":"Domain-like Incommensurate Charge-Density-Wave States and the First- Order Incommensurate-Commensurate Transitions in Layered Tantalum Dichalcogenides. II. 2H- Polytype,","cited_arxiv_id":null,"evidence_quote":"supplies the 2H-TaSe2 free-energy parameters and the IC-C transition analysis that the paper generalizes."},{"cited_title":"ORIGIN OF THE STABILIZATION OF THE NEARLY COMMENSURATE PHASE IN 1T-TaS2,","cited_arxiv_id":null,"evidence_quote":"provides the ring-like diffuse-scattering data used to fix the type-1 free-energy coefficients for 1T-TaS2."},{"cited_title":"Theory of Three-Dimensional Orderings of Charge-Density Waves in 1T- TaX2 (X: S, Se),","cited_arxiv_id":null,"evidence_quote":"is the inter-layer stacking explanation of stripe and T phases that the paper argues is not necessary in monolayers."},{"cited_title":"X-ray study of the new charge-density-wave phase in 1T-TaS2,","cited_arxiv_id":null,"evidence_quote":"provides the experimental T-phase domain-wall angles and domain sizes in 1T-TaS2 used to compare with the calculated quasi-stripe structure."},{"cited_title":"Direct observation of monolayer, bi- layer, and tri-layer charge density waves in 1T-TaS2 by transmission electron microscopy without a substrate,","cited_arxiv_id":null,"evidence_quote":"provides the transmission-electron-microscopy observation of multiple T domain walls in a free-standing monolayer 1T-TaS2 that motivates the monolayer calculation."},{"cited_title":"Broken Hexagonal Symmetry in the Incommensurate Charge-Density Wave Structure of 2H-TaSe2,","cited_arxiv_id":null,"evidence_quote":"provides the experimental broken-hexagonal stripe structure of 2H-TaSe2 used as the stripe-phase comparison."}],"review_version":1}