{"id":"6d396402-28b5-4e5c-9877-de6775a1667c","arxiv_id":"2411.15654","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Equilibrium 2D Gaussian-core systems exhibit string-like cooperative motion and persistent sub-diffusive dynamics in their solid and hexatic phases, most strongly under random pinning.","lead":"This simulation study shows that particles in two-dimensional Gaussian-core systems move cooperatively along string-like paths in low-temperature equilibrium phases, even without disorder. Random pinning amplifies this cooperative motion and produces sub-diffusive, non-Fickian dynamics that resemble glassy relaxation.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central claim depends on sub-diffusive γL(t) persisting as an equilibrium property; the evidence covers only about one decade and the authors concede they cannot establish longer-time persistence, so the contrast with earlier Fickian-crossover results is not yet supported.","rationale":"The reader identified the weakest assumption as the long-time persistence of sub-diffusive γL(t), and I agree: this is the load-bearing point because the paper's novelty claim is that cooperative motion causes a persistent departure from Fickian diffusion in equilibrium phases, in contrast with earlier colloidal studies that observed a crossover to Fickian behavior. All of the supporting evidence is limited to one system size (N=4096), one density, and a time window of about one decade for the exponent fit, with no error bars on β. The authors' explicit caveat that their data quality is inadequate to establish longer-time persistence confirms that the central interpretation is not yet secure. The qualitative observation of string-like cooperative motion and non-Gaussian tails is plausible and visually supported, so I would not reject the paper; however, the central claim should remain conditional until the finite-time-transient alternative is tested. A longer-time and larger-system simulation, plus a check of the lab-frame MSD, would directly determine whether the sub-diffusive regime is asymptotic or a pre-asymptotic crossover, and would also address whether the cage-relative coordinate choice creates the apparent slow relaxation. Therefore the reader's CONDITIONAL verdict is appropriate and no change is needed.","tokens_in":12717,"tokens_out":5257,"duration_ms":54965,"concrete_test":"Extend MD runs for the pure system at Γ−1 = 0.0112 and the RP system at Γ−1 = 0.0056 to at least t = 10^6 reduced units, with N = 4096 and one larger system N = 16384 at fixed density ρ = 0.628. Compute the local log-log slope of γL(t) in non-overlapping windows [10^3, 10^4], [10^4, 10^5], and [10^5, 10^6]. If the local slope increases monotonically toward 1, the sub-diffusive CMTR is a finite-time transient. Separately compute the lab-frame MSD; if it becomes Fickian while the cage-relative γL(t) remains sub-diffusive, the claimed 'slow relaxation' is coordinate-dependent rather than an intrinsic equilibrium transport property.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline claim is that in the CMTR, γL(t) ∼ t^β with β<1 persists up to t≈10^5 and that this 'departure from diffusive dynamics' is an equilibrium property of pure and disordered Gaussian-core systems, in contrast with earlier findings (Refs. 26, 27, 47). This claim rests on the assumption that the observed sub-diffusive window is not a finite-time transient that would eventually cross over to Fickian diffusion. The paper itself concedes the uncertainty: after reporting β<1 up to t≈10^5, it states 'our data quality is inadequate for conclusively establishing this' in the RP-system discussion and asks 'Will these non-Fickian dynamics survive for longer times?' without answering. For the pure system, the fitted regime shown in the Fig. 1(a) inset is roughly one decade (10^4 ≤ t ≤ 10^5), no error bars are provided for β, and only N=4096 is used with no system-size variation. Because cage-relative displacements suppress the dominant phonon background, the residual growth of γL(t) is controlled by rare, intermittent cage-breaking events whose statistics are poor at these times and with only 10 pinning realizations. A gradual upward curvature of γL(t) over this single decade could masquerade as β<1. In any finite equilibrium system with finite relaxation time, a cage-relative mean-square displacement must eventually become linear in t; the only question is the crossover time. The paper provides no evidence that the crossover time exceeds the simulation window, so the 'slow relaxation in contrast with earlier findings' is not established as an asymptotic equilibrium property.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"Using molecular dynamics simulations of the two-dimensional Gaussian-core model, the authors study dynamics across melting in pure systems and systems with random or commensurate pinning. They report a Cooperative Motion Temperature Regime (CMTR) overlapping the solid and hexatic phases in which a macroscopic fraction of particles move along string-like paths, the cage-relative Lindemann parameter grows as t^β with β<1 up to t≈10^5, the self-part of the van Hove function develops a non-Gaussian tail with an unshifted peak, and the non-Gaussian parameter remains nonzero at long times. The effect is stronger under random pinning, which suppresses true solidity, and weaker under commensurate pinning, which stabilizes the solid. The authors argue that this cooperative motion causes a departure from diffusive dynamics and slow relaxation of a glassy type.","tokens_in":13076,"tokens_out":4996,"duration_ms":43460,"significance":"If the result holds, the observation of string-like cooperative motion and non-Fickian cage-relative displacement in equilibrium crystalline and hexatic phases of a purely repulsive two-dimensional system is significant, as it connects glassy dynamics to equilibrium solid phases and challenges the usual association of such heterogeneous dynamics with supercooled liquids or strong disorder. The paper uses several complementary observables (cage-relative Lindemann parameter, van Hove functions in self and distinct parts, non-Gaussian parameter, and the stretching exponent k(t)), and the qualitative signatures are consistent across them. However, the central quantitative claim—persistence of β<1 to t≈10^5—rests on a single decade of data without statistical uncertainties or system-size checks, and the authors themselves acknowledge that the data quality is inadequate to establish longer-time persistence. The significance would be enhanced by a more careful treatment of the long-time limit and by a demonstration that the observed behavior is not a finite-time transient.","major_comments":[{"comment":"The central claim of the paper is that sub-diffusive growth of the cage-relative Lindemann parameter, γL(t)∼t^β with β<1, persists in the CMTR up to t≈10^5 and constitutes a departure from diffusive dynamics. As presented in Fig. 1(a) and its inset, the fitted regime spans roughly one decade (10^4≤t≤10^5), no error bars are provided for β, and only N=4096 is used with no system-size variation. The authors themselves state in the discussion of the RP system, 'our data quality is inadequate for conclusively establishing this', and pose the question 'Will these non-Fickian dynamics survive for longer times?' without answering it. In any finite equilibrium system, a cage-relative mean-square displacement must eventually become linear in t once relaxation is complete; the paper provides no evidence that the crossover time lies beyond the simulation window. The abstract's statement that cooperative motion 'causes a departure from the diffusive dynamics, causing slow relaxation' is therefore stronger than the data support, and the contrast drawn with Refs. [26,27,47] is not yet established. The authors should either extend the simulations (e.g., longer times, larger N, multiple realizations) to test for a crossover to β=1, or explicitly weaken the persistence claim.","section":"Model & method (Fig. 1 and following discussion)"},{"comment":"The exponent β(T) and the stretching exponent k(t) are extracted from fits in log-log or semi-log representations, but the manuscript does not describe the fitting procedure, the choice of fitting window, or the statistical uncertainties. For example, the claim β<1 within 10^4≤t≤10^5 is based on an inset without error bars, and the k(t) analysis in Fig. 2(d-f) is reported for individual time points; the exponential fit for k(t) is shown for only one case in Fig. S4. Without these details, the reader cannot assess whether a gradual upward curvature of γL(t) over the single decade could be misidentified as a persistent sub-diffusive exponent. Please provide error bars, define the fitting windows, and report the sensitivity of β and k to those choices.","section":"Measure of the non-Gaussian motion (Fig. 2 and Fig. S4)"},{"comment":"The comparison with earlier findings (Refs. [26,27,47]) is used to argue that the persistent sub-diffusive behavior and non-vanishing NGP at long times are 'in contrast with earlier findings.' However, the cited works study different systems (colloidal crystals, quasi-two-dimensional colloids) and mostly monitor standard mean-square displacements rather than cage-relative displacements. Since the use of cage-relative coordinates suppresses phonon contributions and changes the long-time behavior, the claimed contrast is not established unless the same observables are compared for those systems, or the coordinate dependence is explicitly addressed. The discussion in Section VI of the SM is a start, but it does not directly address the Fickian crossover in the earlier works.","section":"Non-Gaussian parameter (discussion of contrast with Refs. [26,27,47])"}],"minor_comments":[{"comment":"The sentence 'The resulting distribution of displacements, even in pure systems, becomes non-Gaussian and causes slow relaxation akin to glassy systems' presents a result before it has been derived; consider rephrasing as a preview or moving it to the conclusions.","section":"Introduction"},{"comment":"The dimensionless coupling Γ^{-1} is defined only in Section I of the SM; stating its definition (or at least a pointer) in the main text would help readers who do not consult the SM.","section":"Model & method"},{"comment":"There is a typo: 'where the the system is presumably a solid' should read 'where the system is presumably a solid.'","section":"Main text after Fig. 1(b)"},{"comment":"The figure captions are very long and contain much of the analysis narrative; the authors could streamline them and move interpretive statements to the text, which would improve readability.","section":"Figure captions"},{"comment":"The statement 'The effect of such atypical motion is amplified in the presence of uncorrelated disorders, where they are found down to the lowest temperatures, denying solidity' is a strong claim; it would benefit from a quantitative measure of 'denying solidity' (e.g., the vanishing of the shear modulus or the behavior of the Lindemann parameter in the low-T limit).","section":"Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The paper leans on the authors' prior static phase diagram (Ref. [16]) to locate the solid/hexatic/liquid boundaries, and the dynamic CMTR boundaries are presented as overlapping those static phases. If Ref. [16] or its interpretation is questioned, the mapping of CMTR to equilibrium phases would need to be revisited. The manuscript is letter-style; a revised version should include more details on the error analysis and simulation protocol in the SM to support the central quantitative claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe paper reports string-like cooperative motion in the equilibrium solid and hexatic phases of a pure two-dimensional Gaussian-core model, with sub-diffusive cage-relative mean-square displacement persisting up to t≈10^5. That observation, if it holds up, is genuinely new: prior work saw such cooperative dynamics only in supercooled liquids, near melting, or with disorder. The paper also shows random pinning extends this regime to the lowest temperatures, while commensurate pinning anchors the solid — a nice contrast.\n\nWhat the paper does well: the analysis uses several independent observables — cage-relative Lindemann parameter, van Hove function, non-Gaussian parameter, k(t) exponent, and the distinct van Hove peak at r→0. The visual evidence of string-like paths is compelling, and the use of cage-relative coordinates to suppress phonon background is appropriate. The authors are also honest: they explicitly say their data quality is inadequate to establish whether the sub-diffusion persists beyond their window.\n\nThe soft spot is exactly there. The claim of persistent sub-diffusion, and the contrast with earlier Fickian-crossover results, rests on about one decade of data, with no error bars on β, no system-size variation, and only 10 pinning realizations. A gradual upward curvature over that decade could masquerade as β<1. The authors themselves ask 'Will these non-Fickian dynamics survive for longer times?' and don't answer. So the asymptotic claim is not established. That said, this is a limitation on a quantitative extrapolation, not on the core qualitative observation: cooperative, string-like motion in an equilibrium pure solid is already interesting, and the non-Gaussian tails and NGP are clear.\n\nMinor issues: single density, single pinning fraction, and phase boundaries taken from the authors' own earlier paper. Acceptable, but a referee should ask for a second density or a check against another interaction.\n\nWho this is for: people working on 2D melting, colloidal dynamics, and glassy behavior in equilibrium soft matter. It deserves a serious referee — the observation is novel enough that the field needs to know about it, even if the long-time interpretation remains open.\n\nMy recommendation: send it to peer review, and let the referees push for longer trajectories or, at minimum, a softened claim about persistence. With that revision, it's a solid contribution.","headline":"A novel observation of string-like cooperative motion in pure equilibrium 2D crystals, with a real but contained caveat about the persistence of sub-diffusion.","tokens_in":13629,"tokens_out":2646,"would_cite":true,"duration_ms":25429,"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":"In equilibrium 2D solids, cooperative strings slow diffusion.","keywords":["two-dimensional melting","Gaussian-core model","cooperative motion","dynamical heterogeneity","sub-diffusion","cage-relative displacement","random pinning","hexatic phase"],"falsifier":"Run the same Gaussian-core simulations (pure and random-pinned) for times well beyond $t \\approx 10^5$, or in larger systems, and measure whether the exponent $\\beta(t)$ of $\\gamma_L(t)$ remains below 1; if $\\beta$ approaches 1 and the van Hove tails become Gaussian at longer times, the claimed CMTR is a transient. Alternatively, track individual particles in a colloidal monolayer experiment at the same effective temperatures and look for string-like cooperative displacements with a sub-diffusive Lindemann parameter.","tokens_in":12498,"feed_emoji":"🧊","tokens_out":5207,"duration_ms":43857,"temperature":0.7,"pith_summary":"This paper studies how particles with Gaussian-core repulsion move as a two-dimensional system melts from solid to liquid. It finds that in the solid and part of the hexatic phase, a macroscopic fraction of particles move together along long, tortuous string-like paths, producing sub-diffusive growth of the mean-square displacement and non-Gaussian displacement distributions. This cooperative motion is usually associated with glasses and supercooled liquids, but it appears here in ordinary equilibrium phases, even in a perfectly pure system. Randomly pinned impurities amplify the effect and extend it to the lowest temperatures studied, while impurities placed on the underlying lattice anchor the solid and suppress it. The authors argue this glass-like motion causes slow relaxation and marks a distinct dynamical regime, the Cooperative Motion Temperature Regime, whose boundaries differ from the thermodynamic phase boundaries.","feed_headline":"Cooperative strings slow diffusion in equilibrium 2D solids","feed_subtitle":"Pure 2D solid and hexatic phases show glass-like cooperative strings that slow diffusion down.","key_machinery":"The central object is the cage-relative displacement, $\\Delta r_i(t)$ minus the displacement of a particle's first neighbors, which removes long-wavelength collective motion and isolates local hopping. Normalized by the zero-temperature lattice spacing it defines the dynamic Lindemann parameter $\\gamma_L(t)$, whose long-time exponent $\\beta$ distinguishes solid ($\\beta = 0$), CMTR ($0 < \\beta < 1$), and liquid ($\\beta = 1$). The string-like cooperative paths are identified through trajectory inspection, the non-Gaussian parameter $\\alpha_2(t)$, the stretched-exponential decay exponent $k(t)$ of the van Hove tail, and the rising peak of the distinct van Hove function at $r \\to 0$.","core_discovery":"The central claim is that equilibrium, defect-free two-dimensional Gaussian-core systems pass through a Cooperative Motion Temperature Regime (CMTR) overlapping the solid and part of the hexatic phase, in which particles perform cage-breaking jumps along string-like paths, with a macroscopic fraction moving cooperatively. In this regime the cage-relative dynamic Lindemann parameter grows as $\\gamma_L(t) \\sim t^{\\beta}$ with $\\beta < 1$ persisting to $t \\approx 10^5$, the self-part of the van Hove function develops a long non-Gaussian tail with an unchanged peak position, and the distinct part $G_d(r,t)$ develops a rising peak at $r \\to 0$, the signature of particles occupying vacated sites. The sub-diffusive behavior and non-Gaussian tails persist well beyond the peak of the non-Gaussian parameter, in contrast with earlier claims of a crossover to Fickian diffusion in colloidal systems. Random pinning amplifies the regime, making it span the entire low-temperature phase and denying true solidity, while commensurate pinning anchors the solid to higher temperatures.","pith_inferences":["If the sub-diffusive regime is intrinsic, similar cooperative strings might appear in other soft-matter 2D systems with tunable interactions, not just Gaussian cores.","The apparent 'hexatic glass' in random-pinned systems could be tested by measuring hexatic order decay and comparing its relaxation time with the sub-diffusive exponent.","The authors' use of cage-relative coordinates suggests experiments should analyze relative, not absolute, displacements to detect these effects in colloidal monolayers.","The CMTR boundaries could shift if the interaction potential is steeper or softer; a systematic study across interaction softness would map how generic the phenomenon is."],"forward_implications":["In pure 2D Gaussian-core systems, the solid and part of the hexatic phase host glass-like cooperative motion, so slow relaxation is not exclusive to glasses.","Random pinning at 3.5% suppresses true solidity, extending the CMTR to the lowest temperatures and producing a regime the authors suggest may be a 'hexatic glass'.","Commensurate pinning anchors crystalline order, raising the melting temperature and narrowing the CMTR.","Dynamical boundaries (CMTR) do not coincide with thermodynamic phase boundaries, so dynamics and statics encode different melting pictures.","Cage-relative coordinates reveal these effects cleanly by removing long-wavelength motion; standard coordinates show the same qualitative trends with additional lattice peaks."],"supporting_citations":[{"why":"Supplies the static phase boundaries in pure and disordered Gaussian-core systems that the CMTR is compared against.","marker":"[16]"},{"why":"Defines the cage-relative displacement and dynamic melting criteria used throughout the study.","marker":"[22]"},{"why":"Establishes string-like cooperative motion as a hallmark of glassy dynamics, the phenomenon the paper finds in equilibrium phases.","marker":"[23]"},{"why":"Documents cooperative dynamics in confined colloidal systems, providing a baseline for comparing the present string-like motion.","marker":"[24]"},{"why":"Earlier simulation study of cooperative dynamics in two dimensions that the paper's persistent sub-diffusion contrasts with.","marker":"[26]"},{"why":"Earlier simulation finding of seemingly Fickian but heterogeneous dynamics in 2D colloids, the baseline the paper's long-time sub-diffusion contradicts.","marker":"[27]"},{"why":"Established the hexatic phase in the 2D Gaussian-core model, locating the phase region the CMTR overlaps.","marker":"[30]"},{"why":"Earlier study of dynamical heterogeneities in soft colloidal crystals reporting a Fickian crossover that the paper's long-time behavior challenges.","marker":"[47]"}],"fun_headline_variants":["Glass-like strings slow diffusion in equilibrium 2D solids","Equilibrium 2D solids show glassy cooperative strings","String motion in pure solids mimics glassy slowdown","Cooperative strings drive sub-diffusion in 2D melting","Pinned 2D systems extend glass-like cooperative motion"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the sub-diffusive growth seen up to $t \\approx 10^5$ is a genuine long-time property of the equilibrium phase, not a finite-time transient that would eventually cross over to ordinary diffusion; the authors note their data cannot conclusively rule out that crossover.","fun_headline_variants_meta":{"raw":{"variants":["Glass-like strings slow diffusion in equilibrium 2D solids","Equilibrium 2D solids show glassy cooperative strings","String motion in pure solids mimics glassy slowdown","Cooperative strings drive sub-diffusion in 2D melting","Pinned 2D systems extend glass-like cooperative motion"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000171,"raw_usage":{"total_tokens":1221,"prompt_tokens":845,"completion_tokens":376,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":461,"completion_tokens_details":{"reasoning_tokens":295}},"tokens_in":461,"tokens_out":376,"duration_ms":3638,"temperature":1.0,"reasoning_tokens":295,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:03:59.001524+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same Gaussian-core simulations (pure and random-pinned) for times well beyond $t \\approx 10^5$, or in larger systems, and measure whether the exponent $\\beta(t)$ of $\\gamma_L(t)$ remains below 1; if $\\beta$ approaches 1 and the van Hove tails become Gaussian at longer times, the claimed CMTR is a transient. Alternatively, track individual particles in a colloidal monolayer experiment at the same effective temperatures and look for string-like cooperative displacements with a sub-diffusive Lindemann parameter.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the static phase boundaries in pure and disordered Gaussian-core systems that the CMTR is compared against."},{"cited_title":"Zahn and G","cited_arxiv_id":null,"evidence_quote":"Defines the cage-relative displacement and dynamic melting criteria used throughout the study."},{"cited_title":"Donati, J","cited_arxiv_id":null,"evidence_quote":"Establishes string-like cooperative motion as a hallmark of glassy dynamics, the phenomenon the paper finds in equilibrium phases."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents cooperative dynamics in confined colloidal systems, providing a baseline for comparing the present string-like motion."},{"cited_title":"Zangi and S","cited_arxiv_id":null,"evidence_quote":"Earlier simulation study of cooperative dynamics in two dimensions that the paper's persistent sub-diffusion contrasts with."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Earlier simulation finding of seemingly Fickian but heterogeneous dynamics in 2D colloids, the baseline the paper's long-time sub-diffusion contradicts."},{"cited_title":"Prestipino, F","cited_arxiv_id":null,"evidence_quote":"Established the hexatic phase in the 2D Gaussian-core model, locating the phase region the CMTR overlaps."},{"cited_title":"conserva- tion of area under its trace, it should be multiplied by 2 πr, similar to the plots for Gs(r, t)","cited_arxiv_id":null,"evidence_quote":"Earlier study of dynamical heterogeneities in soft colloidal crystals reporting a Fickian crossover that the paper's long-time behavior challenges."}],"review_version":1}