{"id":"34ec4bc9-4ec9-40df-b3a2-b8557be1e679","arxiv_id":"2608.02456","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"In 1T-TaS2, the 200 K transition randomizes interlayer stacking by coherently shifting Star-of-David clusters without destroying them, while the 350 K transition collapses the clusters themselves.","lead":"Using a machine-learned force field trained on first-principles data, the authors simulate nanosecond dynamics of the charge-density-wave material 1T-TaS2 and find that its two transitions are mechanistically distinct: the first re-arranges Star-of-David clusters while preserving them, the second melts them. This helps explain why ultra-thin films miss the first transition and gives a framework for light- and cavity-driven phase control.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"RDF first-peak shift to 3.28 Å at 250 K contradicts claimed preservation of SoD amplitude through the 200 K transition","rationale":"The reader identified the MLFF/PBE potential as the weakest assumption, which is a valid general concern. My reading found a more specific and immediate conflict in the presented data: the first RDF peak at 250 K coincides with the undistorted Ta–Ta distance, which contradicts the claim that the SoD amplitude is preserved through the first transition. This is not just a question of the functional's accuracy—it is a question of whether the simulations actually exhibit the claimed behavior. The proposed order-parameter check would settle the issue directly. The reader's conditionality is appropriate: the paper should be accepted only if the data are re-analyzed and the claim holds. Since my concern strengthens the reader's conditionality, the verdict should remain UNCHANGED.","tokens_in":9819,"tokens_out":11360,"duration_ms":105446,"concrete_test":"Compute an explicit SoD amplitude order parameter from the MLFF trajectories, e.g., the average projection of each Ta atom's displacement onto the ideal √13×√13 distortion pattern, as a function of temperature (100, 175, 200, 250, 300, 350, 400 K). If this order parameter remains at ≥80% of its 100 K value at 250 K, the RDF peak shift must have another explanation (e.g., thermal expansion or a change in which shell defines the 'first' peak), and the claim is supported. If it drops below ~50% already at 250 K, the central two-step mechanism is contradicted.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on the two transitions being mechanistically distinct: above 200 K, SoD clusters shift while preserving their local amplitude, and only near 350 K does the amplitude collapse. But Fig. 4e shows the first RDF peak for Ta–Ta distances shifts from 3.18 Å at 100 K to 3.28 Å at 250 K, which the authors themselves note corresponds to the Ta–Ta distance in the symmetric triangular lattice without SoD formation. This is above the 200 K transition but well below 350 K. If the SoD amplitude were preserved through the first transition, the intra-cluster short bonds (3.18 Å) should still dominate the first RDF peak; the shift to 3.28 Å suggests those short bonds are largely gone by 250 K. This is either a symptom of MLFF inaccuracy (if it destabilizes the SoD too readily) or a misreading of the RDF, but in either case the presented data undercut the distinction between stacking rearrangement and amplitude collapse. No explicit SoD amplitude order parameter is reported, so this inconsistency is not addressed.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports AIMD and machine-learning force-field molecular dynamics simulations of bulk 1T-TaS2, using up to 1404-atom cells and 5 ns trajectories, to characterize the microscopic mechanism of the two-step CDW transition. The central claim is that the two transitions are mechanistically distinct: the CCDW-NCCDW transition at ~200 K proceeds by coherent SoD shifts via a deformation-formation process that preserves the local SoD amplitude while randomizing interlayer stacking and nucleating domain walls; the NCCDW-ICCDW transition near 350 K involves collapse of the SoD amplitude itself. Evidence includes AIMD at 100 and 175 K, NEB barriers (108 meV for deformation-formation vs 207 meV for layer sliding; 79 meV for domain wall), MLFF-MD at 200-300 K showing SoD shifts with size-dependent transition times, and RDF/DOS analyses.","tokens_in":9934,"tokens_out":5047,"duration_ms":69323,"significance":"If the conclusions hold, this would resolve a long-standing debate on the 1T-TaS2 transition sequence, provide a microscopic explanation of the thickness dependence of the CCDW-NCCDW transition, and give a concrete framework for light-induced hidden phases and cavity-modified transitions. The paper's use of a first-principles-trained MLFF to reach previously inaccessible time/length scales is appropriate, and the distinction between deformation-formation and sliding mechanisms via NEB is a strong point. However, the absence of MLFF validation on extended cells and the inconsistency between the claimed preserved SoD amplitude and Fig. 4e currently prevent full confidence; the central claim is plausible but not yet fully demonstrated.","major_comments":[{"comment":"The central claim that the 200 K transition preserves the local SoD amplitude is difficult to reconcile with Fig. 4e, which shows the first Ta-Ta RDF peak shifting from 3.18 A at 100 K to 3.28 A at 250 K, the latter identified by the authors as the Ta-Ta distance of the undistorted triangular lattice. If the SoD short bonds were preserved, this shortest-peak position should remain near 3.18 A. The text attributes the shift to 'SoD shifts and thermal expansion,' but no direct SoD amplitude order parameter is reported to substantiate that the short-bond population survives. Please provide a quantitative mode-amplitude or bond-population analysis as a function of temperature, or reconcile the RDF shift with the proposed mechanism.","section":"Thermally driven SoD deformation and the NCCDW-ICCDW transition, Fig. 4e"},{"comment":"The MLFF is the workhorse for the central dynamics claims (5 ns, 1404 atoms), yet the manuscript reports no validation of the NequIP model: no training/validation force- or energy-RMSE, no held-out test set, and no direct comparison of MLFF-MD with AIMD on the extended cells or even on the smaller 2 sqrt(13) x sqrt(13) x 2 cell. The transition times and the RDF evolution are therefore only as reliable as the surrogate. Please report held-out errors, a comparison of MLFF and AIMD observables (e.g., RDF, transition events) in the training cell, and if possible a short MLFF-vs-AIMD check on a subset of extended-cell configurations.","section":"Methods: Machine-learning force field and large-scale molecular dynamics"},{"comment":"The claim that the SoD transition time 'increases exponentially with lattice size and temperature' is based on what appear to be single trajectories for a few superlattice sizes and temperatures. No error bars, replicate runs, or first-passage-time statistics are provided. Given that the size dependence is used to argue about experimental observability, this quantitative claim needs at least 3-5 independent runs per condition or an explicit statement that these are single-event estimates.","section":"Results: Stacking-order randomization, Fig. 3d"},{"comment":"The MLFF and NEB barriers are obtained with a PBE-type functional without Hubbard corrections, while the electronic-structure discussion uses on-site Coulomb interactions (Ref. [18]). In a strongly correlated system like 1T-TaS2, semi-local PBE can misestimate relative barriers; the ordering of 108 meV vs 0.35 eV/SoD is load-bearing for the two-step mechanism. Please provide PBE+U (or hybrid-functional) NEB barriers for the deformation-formation and SoD-collapse paths, or otherwise demonstrate that the barrier ordering is robust to the treatment of correlations.","section":"Methods: DFT calculations / NEB"}],"minor_comments":[{"comment":"There is a typographical artifact 'T¿200 K' in the Introduction; it should be 'T > 200 K'. Please check for other non-ASCII artifacts.","section":"Abstract/Introduction"},{"comment":"The text describes single SoD-shift events at specific times (5-9 ps). Please clarify in the caption or text that these are representative single trajectories, not ensemble averages, and state how many AIMD runs were performed.","section":"Results: AIMD, Fig. 2b"},{"comment":"The horizontal axis 'Transition time (ns)' is on a logarithmic scale; please make this explicit in the axis label or caption. Also specify whether 'exponential with lattice size' means with the number of SoDs or with the linear supercell dimension.","section":"Results: Fig. 3d"},{"comment":"Please state explicitly whether the VASP AIMD training data and the Quantum Espresso DFT/NEB calculations use the same functional, pseudopotential scheme, and dispersion corrections. If they differ, the MLFF training and the NEB barrier comparison may not be internally consistent.","section":"Methods: DFT calculations"},{"comment":"Provide a direct URL or accession link to the Figshare repository in addition to the DOI, for reader convenience.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The paper's mechanism is intriguing and likely to attract attention, but the reviewer concerns about MLFF validation and the Fig. 4e inconsistency are substantive; I would not accept in current form. If the authors can provide an explicit SoD amplitude order parameter and MLFF error metrics, the paper could become a strong candidate for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe short version: this is a serious computational kinetics study with a plausible mechanism, but the RDF data as presented seem to undercut the paper's central claim that the 200 K transition preserves SoD amplitude. I'd send it out for review, but with major revision.\n\nWhat is new: they actually resolve the two-step dynamics with an MLFF trained on PBE AIMD, reaching 1404-atom cells over nanoseconds. The NEB barriers (108 meV for deformation-formation vs 207 meV for sliding, 79 meV for domain wall formation) are a concrete, useful result. The thickness-dependent interpretation—that stacking randomization, not domain walls alone, drives metallization—is worth considering. I give them credit for extracting a coherent mechanistic picture from simulation.\n\nThe soft spots: first, the RDF contradiction. Fig 4e shows the first Ta-Ta RDF peak shifts to 3.28 Å by 250 K—exactly the undistorted triangular lattice value. The authors claim this is due to \"SoD shifts and thermal expansion,\" but 0.1 Å is far too large for thermal expansion, and if the SoD amplitude were preserved, the short intra-cluster bonds should keep the peak near 3.18 Å. This looks like the MLFF is dissolving the SoD amplitude too readily, or the authors are misreading a time-averaged RDF. No direct SoD amplitude order parameter is reported, so we can't distinguish. This is not a minor issue; it goes to the core of the mechanism.\n\nSecond, the MLFF validation is thin. Trained on small AIMD cells at 100-500 K; no held-out test error, no validation on the extended cells. The transition-time scaling (Fig 3d) comes from single trajectories, and the claim of exponential scaling with cell size is an inference from a few points with no error bars. The sentence about temperature dependence of transition time appears garbled.\n\nThat said, the high-level picture—first transition rearranges stacking/domain walls, second melts the distortion—is consistent with earlier stacking-driven models and with the thickness dependence, so the overall direction is probably right. But the paper's own RDF data don't yet support it cleanly.\n\nWho's it for: TaS2 specialists and anyone using MLFF for correlated materials. It deserves a serious referee because the question is important and the NEB barriers are a real contribution, but it needs major revision: direct amplitude order parameter, proper MLFF validation, replicate trajectories, and a reconciliation of the RDF.\n\nI'd accept it for peer review but not as it stands.","headline":"A serious computational kinetics study with a clean mechanistic story, but the RDF data as presented seem to undermine the central claim that the 200 K transition preserves SoD amplitude; needs direct order parameter and MLFF validation before I'd trust the distinction.","tokens_in":10563,"tokens_out":4415,"would_cite":false,"duration_ms":38993,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["71.45.Lr","71.30.+h","63.22.-m"],"model":"deepseek-v4-flash","headline":"The 200 K step of the CDW transition in 1T-TaS2 shifts Star-of-David clusters while preserving their distortion; the 350 K step melts that distortion, and the paper proves this with machine-learning molecular dynamics.","keywords":["charge density wave","1T-TaS2","Star-of-David distortion","machine learning force field","molecular dynamics","stacking order","domain walls","metal-insulator transition"],"falsifier":"Measure the intensity of the sqrt13 x sqrt13 superlattice reflection and the interlayer stacking correlation in the same time-resolved diffraction experiment while heating through 200 K: if the superlattice peak weakens or the peak width broadens at 200 K, the claim that the SoD amplitude is preserved would be falsified; only a loss of stacking correlation without superlattice attenuation would support the paper's mechanism.","tokens_in":1316,"feed_emoji":"🔬","tokens_out":3871,"duration_ms":60534,"temperature":0.7,"pith_summary":"This paper uses machine-learning-force-field molecular dynamics to watch the 200 K and 350 K transitions in bulk 1T-TaS2 in real time. It claims the two steps are mechanistically distinct: above 200 K, Star-of-David clusters translate by transiently dissolving and re-forming around shifted centers, preserving their local amplitude while randomizing interlayer stacking and nucleating domain walls. Only near 350 K does the Star-of-David distortion itself collapse. This resolves a long-standing question about whether the transitions are driven by amplitude collapse or lattice rearrangement, and it explains why the 200 K transition disappears in thin films.","feed_headline":"Star-of-David clusters shift first, melt at 350 K","feed_subtitle":"Simulations reveal the 200 K step preserves the local distortion while randomizing stacking; the 350 K step destroys it.","key_machinery":"The central mechanism is the deformation-formation process: a low-barrier path for SoD cluster translation (108 meV, versus 207 meV for layer sliding) where the cluster's internal distortion weakens just enough for the cluster to re-center, then re-forms. The authors access this dynamics with a NequIP machine-learning force field trained on PBE-based AIMD trajectories, enabling 5-ns simulations of 1404-atom supercells that capture collective coherent breathing motions and domain-wall nucleation.","core_discovery":"The paper establishes that the commensurate-to-near-commensurate transition (200 K) is a lattice-rearrangement event: Star-of-David clusters shift coherently through a deformation-formation process in which the local displacement pattern temporarily dissolves and re-forms around a new center, leaving the SoD amplitude intact. This shifting randomizes the interlayer stacking order and creates metastable in-plane domain walls, both of which contribute to metallization. The near-commensurate-to-incommensurate transition (350 K) is instead an amplitude-collapse event, where the SoD distortion fully melts into a symmetric triangular lattice. The two transitions therefore have different order para","pith_inferences":["A direct experimental test could use time-resolved electron or X-ray diffraction just above 200 K: if the superlattice peak intensity stays constant while the interlayer stacking correlation decays, the deformation-formation picture is confirmed.","The reliance on a PBE-trained force field implies that strong-correlation renormalization of the 0.35 eV/SoD collapse barrier might shift the absolute temperature of the second transition, even though the qualitative two-step ordering is robust.","The exponential size dependence of the shifting time suggests that in real crystals, the CCDW-NCCDW transition might be governed by rare, locally coherent events rather than global cooperative motion, which could be probed by carefully measuring pre-transitional dynamics in nanoscale samples."],"forward_implications":["If correct, the metallization at the 200 K transition is driven primarily by randomization of interlayer stacking, not by in-plane domain walls alone, explaining why thin films (which lack out-of-plane dispersion) do not show the transition.","The 350 K transition is a genuine melting of the SoD distortion, independent of stacking configuration, consistent with its appearance in all sample thicknesses.","The transition time for SoD shifting grows exponentially with lattice size because coherence of the breathing motion becomes harder to achieve, suggesting that experimentally observed domain-wall dynamics are nucleation-limited.","Light-induced hidden phases should retain the intralayer SoD amplitude while losing interlayer correlations, pointing to a coherent deformation-formation process as the underlying route.","Cavity-modified transition temperatures can be interpreted as renormalization of the 108 meV shifting barrier rather than a change in the electronic correlation strength."],"fun_headline_variants":["Two-step CDW melt: first shift, then dissolve","Simulation: SoD clusters shift, then amplitude collapses","200 K step moves clusters; 350 K step melts them","Real-time simulation reveals CDW two-step: shift then melt","Star-of-David clusters slide first, then melt away"],"cache_read_input_tokens":11776,"weakest_assumption_plain":"The entire dynamical picture rests on the machine-learning force field faithfully reproducing the density-functional-theory potential-energy surface, especially the relative heights of the 108 meV shifting barrier and the 0.35 eV/SoD collapse barrier; if strong electron correlations alter these barriers, the mechanistic ordering of the two transitions could change.","fun_headline_variants_meta":{"raw":{"variants":["Two-step CDW melt: first shift, then dissolve","Simulation: SoD clusters shift, then amplitude collapses","200 K step moves clusters; 350 K step melts them","Real-time simulation reveals CDW two-step: shift then melt","Star-of-David clusters slide first, then melt away"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000877,"raw_usage":{"total_tokens":3637,"prompt_tokens":757,"completion_tokens":2880,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":501,"completion_tokens_details":{"reasoning_tokens":2797}},"tokens_in":501,"tokens_out":2880,"duration_ms":18173,"temperature":1.0,"reasoning_tokens":2797,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T06:57:41.499926+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the intensity of the sqrt13 x sqrt13 superlattice reflection and the interlayer stacking correlation in the same time-resolved diffraction experiment while heating through 200 K: if the superlattice peak weakens or the peak width broadens at 200 K, the claim that the SoD amplitude is preserved would be falsified; only a loss of stacking correlation without superlattice attenuation would support the paper's mechanism.","supporting_citations":[],"review_version":1}