{"id":"c7a2c291-1481-46f1-8e00-eeac105e4ac7","arxiv_id":"2608.06797","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The amplitude of the t^(3/2) short-time relaxation of density correlations is exactly set by equilibrium contact geometry in reversible Brownian hard-particle systems.","lead":"This paper derives a formula for the strength of the first fractional (t to the 3/2) decay of density correlations in dense Brownian liquids, tying it to hard particle contacts. If correct, it gives experimentalists and simulators a boundary condition to separate short-time collision physics from later glassy relaxation.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Spurious factor 1/2 in Eq. (C4) makes the printed derivation of the hard-sphere formula (12) off by a factor 4; the central prediction is not actually derived as written.","rationale":"The reader's conditional verdict is reasonable but pointed at the boundary-layer asymptotic argument. My independent check found a more immediate and fully concrete problem in the hard-sphere specialization. The spurious 1/2 in Eq. (C4) changes the predicted amplitude by a factor 4; since Eq. (12) is the central experimentally testable result, the derivation as printed is internally inconsistent. This is not an objection to the physical claim itself—the corrected derivative reproduces the known prefactor—but it means the manuscript must be revised before the 'exact, fit-free' claim is supported. The SQS soft-interface check and the projection hierarchy are independent and unaffected. The boundary-layer asymptotics in Appendix C should still be made rigorous, but the factor error is the sharper concern. No ad hominem; this is a checkable algebraic point. I therefore keep the reader's CONDITIONAL verdict, with an additional specific required correction.","tokens_in":9880,"tokens_out":32722,"duration_ms":287782,"concrete_test":"Take the hard-sphere reduction in Appendix C and recompute Eq. (12) from Eq. (11) using the directly differentiated normal derivative \\partial_{h_{ij}}\\rho_k = i k\\cdot\\hat n_{ij}(e^{ik\\cdot r_j}-e^{ik\\cdot r_i}) (no 1/2), keeping D_{ij}=2D_0 and the pair identity (C5). If the resulting prefactor is (4\\sqrt{\\pi}/9)\\rho\\sigma^2 g(\\sigma^+)(2D_0)^{3/2} k^2[1-j_0(k\\sigma)+2j_2(k\\sigma)]/S(k), then Eq. (C4) contains a typographical factor-1/2 error and Eq. (12) is correct after correction; if the prefactor differs, Eq. (12)'s amplitude is wrong by a factor and the central prediction fails as printed. A second check: evaluate the same substitution with the printed C4 and confirm the prefactor is one quarter of Eq. (12), isolating the inconsistency.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The manuscript's central checkable prediction is Eq. (12), the hard-sphere t^{3/2} amplitude. Appendix C is supposed to derive it from Eq. (11). Eq. (C4) states \\partial_{h_{ij}} \\rho_k = (i k\\cdot\\hat n_{ij}/2)(e^{ik\\cdot r_j}-e^{ik\\cdot r_i}). Direct differentiation of h_{ij}=|r_j-r_i|-\\sigma gives \\partial h_{ij}/\\partial r_j=\\hat n_{ij} and \\partial h_{ij}/\\partial r_i=-\\hat n_{ij}, while \\partial\\rho_k/\\partial r_j = i k e^{ik\\cdot r_j} and \\partial\\rho_k/\\partial r_i = i k e^{ik\\cdot r_i}; hence \\partial_{h_{ij}}\\rho_k = i k\\cdot\\hat n_{ij}(e^{ik\\cdot r_j}-e^{ik\\cdot r_i}), with no factor 1/2. Substituting the printed C4 into Eq. (11) and using the pair identity (C5) and integral (C6) yields B_k^{HS} = (2/9)\\sqrt{2\\pi}\\,\\rho\\sigma^2 g(\\sigma^+) D_0^{3/2} k^2[1-j_0(k\\sigma)+2j_2(k\\sigma)]/S(k), exactly one quarter of Eq. (12). Since Eq. (12) is the paper's advertised fit-free prediction, the derivation as written does not produce the stated amplitude. This is a concrete algebraic inconsistency, independent of the boundary-layer asymptotics question; it must be corrected (remove the 1/2 in C4) and the final formula re-verified before the exactness claim can be accepted. The independent literature agreement of the angular factor does not resolve the prefactor ambiguity.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper addresses the origin of the t^{3/2} term in the short-time decay of density-density correlations in reversible overdamped Brownian systems with hard contacts. It derives a surface formula (Eq. 11) expressing the amplitude B_k as an equilibrium average over contact faces of the squared conormal density flux divided by the square root of the normal diffusivity. For monodisperse hard spheres this reduces to a fit-free prediction (Eq. 12) in terms of S(k), g(σ+), and D0. The authors also derive the corresponding normalization for a soft internal interface, verify the normalization against an exactly solvable one-gap SQS model, and prove via a Gram-Schur projection argument that adding regular collective modes does not change the leading contact coefficient. The derivation is self-contained: the inputs are independent equilibrium and short-time transport quantities, and the SQS check is evaluated exactly from the stated expressions.","tokens_in":10221,"tokens_out":14819,"duration_ms":123848,"significance":"If the central formulas are correct, the paper provides a useful microscopic boundary condition for memory-kernel reconstructions and connects contact physics to mode-coupling and collective relaxation. The general surface formula (11) and the projection invariance are conceptually attractive, and the exactly solvable SQS test with 80-digit evaluation is a concrete strength. The hard-sphere angular factor agrees with earlier classical calculations, and the derivation uses only independent measurable inputs. The main reservation is that the printed derivation of Eq. (12) contains an algebraic factor error in Appendix C that must be fixed; because this error affects the headline prediction, the paper requires a major revision rather than minor polishing.","major_comments":[{"comment":"The derivative of the density phase with respect to the contact coordinate is misstated. For h_ij=|r_j-r_i|-σ, one has ∂h_ij/∂r_j = n̂_ij and ∂h_ij/∂r_i = -n̂_ij, and with ρ_k = Σ_l e^{ik·r_l} the chain rule gives ∂_{h_ij}ρ_k = i k·n̂_ij (e^{ik·r_j} - e^{ik·r_i}), with no factor 1/2. Substituting the printed Eq. (C4) into Eq. (11) and using Eqs. (C5)-(C6) produces exactly one quarter of Eq. (12). Removing the spurious 1/2 restores Eq. (12). Because Eq. (12) is the advertised fit-free prediction, this is a load-bearing algebraic error that must be corrected and the resulting formula re-verified before the central claim is accepted.","section":"Appendix C, Eq. (C4)"},{"comment":"The paper states that freezing the mobility tensor and contact geometry over a boundary layer of thickness O(s^{-1/2}) and discarding tangential gradients, curvature, and D-variation introduces errors only at O(s^{-1}), and on this basis calls Eq. (11) exact. This separation of scales is plausible for smooth coefficients, but it is asserted rather than proved. Since the exactness of the leading s^{-1/2} coefficient is one of the paper's central claims, the authors should either supply a short argument that these terms contribute only at O(s^{-1}) or O(s^{-3/2}) under the stated piecewise-smoothness hypotheses, or explicitly reformulate the claim as an asymptotic exactness statement under that separation assumption.","section":"Appendix C, Eqs. (C1)-(C3)"}],"minor_comments":[{"comment":"The definition 'A = -Ω^†' uses a symbol Ω that is never defined; presumably it is a typo for the generator of the diffusion or for L^†, and it should be corrected.","section":"Section II, Eq. (3)"},{"comment":"The acronym SQS is introduced without expansion; either define it explicitly or point the reader to the infinite-dimensional gap model of Refs. [37-39] when the abbreviation first appears.","section":"Section IV"},{"comment":"The label R0,reg(s) appears in the caption but is only defined in the text around Eq. (25); the caption should state that this is the residual after subtraction of the contact block.","section":"Figure 4 caption"},{"comment":"The parabolic-cylinder function Dν is invoked without a definition or reference; since the appendix claims an exact resolvent evaluation, the notation should be specified or a standard reference should be given.","section":"Appendix D, Eq. (D3)"}],"recommendation":"major_revision","confidential_remarks":"The factor 1/2 in Eq. (C4) appears to be a typo, but it sits at the center of the paper's main quantitative claim, which is why I recommend major revision rather than minor. The paper should also make explicit how Eq. (12) relates to the prior hard-sphere short-time results [3,14], since the angular factor itself is not new; the novelty lies in the general surface formula and in the soft-interface and projection-hierarchy results. I do not see grounds for concern about the independence of the inputs or about improper citation behavior."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's real content is the general contact-surface formula (11) for reversible overdamped Brownian systems, plus the soft-interface normalization (16) and the projection-invariance argument. Those are worth a look. The hard-sphere result (12) is not new—the paper honestly credits the angular factor to the classical calculations—so the value added is the general framework and the SQS exact check, which are clean and plausible.\n\nBut there is a concrete problem the reader missed. The stress-test note is right: Eq. (C4) has a factor 1/2 in the density derivative, and that factor is not there. Differentiating h_ij directly gives ∂h_ij ρ_k = i k·n̂(e^{ik·r_j} − e^{ik·r_i}), no 1/2. With the printed C4, substitution into (11) yields exactly one quarter of Eq. (12). So as written, the appendix does not derive the paper's advertised fit-free amplitude. This is a fixable typo, not a conceptual flaw—remove the 1/2 and the derivation goes through—but it is load-bearing because Eq. (12) is the headline checkable prediction. The paper needs this corrected before the exactness claim is taken seriously.\n\nThe other soft spot is exactly what the reader flagged: the boundary-layer asymptotic in Appendix C is a sketch. The assertion that tangential gradients, curvature, and D-variation enter only at O(s^{-1}) is plausible but not rigorously established. That is a lesser concern than the prefactor, and the SQS check lends support, but a serious referee would want a cleaner statement or a proof.\n\nThe citation pattern and inputs are honest. S(k), g(σ+), and D0 are independent measurable quantities; there is no circularity. Self-citations are unrelated to the derivation. The SQS resolvent values are reproduced from stated expressions, which is good practice.\n\nWho gets value from this? People working on memory-kernel reconstructions and MCT-like short-time boundary conditions. The general surface formula, once the prefactor bug is fixed, is a useful organizing result and a legitimate reference for the t^{3/2} amplitude. As is, I would not cite Eq. (12) without running the algebra myself, and I would not pass it to a student without a warning about C4.\n\nRecommendation: send it to peer review, but make the author fix C4, re-verify Eq. (12), and tighten the boundary-layer argument. The core idea deserves referee time; the current draft does not justify its central formula as printed.","headline":"The general surface formula is interesting and likely right, but the printed derivation of the headline hard-sphere prediction is off by a factor 4 because of a spurious 1/2 in Eq. (C4).","tokens_in":10765,"tokens_out":3934,"would_cite":false,"duration_ms":32785,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Hard-sphere contacts exactly determine the leading $t^{3/2}$ density-relaxation amplitude in Brownian liquids.","keywords":["density correlation","hard-sphere contacts","fractional relaxation","memory kernel","Brownian dynamics","intermediate scattering function","contact boundary layer","short-time diffusion"],"falsifier":"Run overdamped Brownian dynamics of monodisperse hard spheres over a short-time window, independently measure $S(k)$, $g(\\sigma^+)$, and $D_0$, and plot $[\\phi_k(t)-1+\\Gamma_k t]/t^{3/2}$; if this ratio does not approach the predicted $B_k$ from Eq. (12) across a range of wave numbers, the claimed exactness fails.","tokens_in":9640,"feed_emoji":"🔬","tokens_out":7908,"duration_ms":62649,"temperature":0.7,"pith_summary":"This paper claims that the leading fractional term in the short-time decay of density correlations in dense Brownian liquids---the $t^{3/2}$ correction that appears before collective relaxation---is exactly fixed by hard particle contacts. Treating each contact as a reflecting boundary in configuration space, it derives a surface formula for the amplitude $B_k$ involving the equilibrium contact probability, the mobility normal to the contact surface, and the response of the density wave to contact motion. For monodisperse hard spheres, the formula becomes an explicit, fit-free prediction in terms of the static structure factor $S(k)$, the radial distribution at contact $g(\\sigma^+)$, and the short-time diffusion coefficient $D_0$. If correct, the result supplies an exact microscopic boundary condition connecting collision kinetics to later caging and structural relaxation within the same memory equation.","feed_headline":"Hard-sphere collisions fix the t^{3/2} decay of density correlations","feed_subtitle":"A fit-free formula from structure, contact value, and diffusion links collision kinetics to collective relaxation.","key_machinery":"The central object is the exact variational Schur complement $\\Sigma_k(s)=\\sup_{v\\in Q\\mathcal{D}(E)}\\{2\\operatorname{Re}\\ell_k(v)-s\\|v\\|^2-\\mathcal{E}(v,v)\\}$, which supplies the Laplace-space memory kernel in the density correlation equation. At a hard contact, the density mode carries a finite boundary flux $q_{k,c}$, so the generator action becomes a surface functional; evaluating the local variational problem on a half-line gives the $s^{-1/2}$ term whose coefficient is the contact-surface formula. The soft-interface calibration uses the corresponding full-line response function, yielding half the reflecting normalization. The projection hierarchy then shows, through a first-column identity for the representing vector of the contact functional, that regular modes contribute only at $O(s^{-1})$ and leave the $s^{-1/2}$ contact term unchanged.","core_discovery":"The paper proves that for reversible overdamped Brownian systems with piecewise-smooth regular contact faces and positive normal mobility, the coefficient $B_k$ of the $t^{3/2}$ term in $\\phi_k(t)=1-\\Gamma_k t+B_k t^{3/2}+O(t^2)$ is exactly $B_k=\\frac{4}{3\\sqrt{\\pi}}\\sum_c\\langle\\delta(h_c)|q_{k,c}|^2/\\sqrt{D_c}\\rangle_\\mu$, an equilibrium average over every contact surface. Mechanistically, diffusion samples a layer of thickness $O(\\sqrt{t})$ at each reflecting contact, and the normal displacement contributes a factor $t$, producing the fractional power; the amplitude collects the many-body weight. For identical hard spheres with scalar diffusivity $D_0$, this reduces to a closed expression in $S(k)$, $g(\\sigma^+)$, $D_0$, and $\\sigma$ whose angular factor $1-j_0(k\\sigma)+2j_2(k\\sigma)$ matches classical short-time hard-sphere calculations and respects density conservation with $B_k=O(k^4)$ as $k\\to 0$. The same construction calibrates soft internal interfaces, with half the reflecting-boundary normalization, and proves via a projection hierarchy built on the shifted Gram matrix that regular collective variables do not alter the leading contact coefficient.","pith_inferences":["If the surface-integral formula holds for a broad class of mobility tensors, the same construction could be extended to anisotropic or position-dependent diffusivity in mixtures and colloidal systems with hydrodynamic interactions, since $D_c$ already enters as a local quantity.","The factor-of-two half-line versus full-line normalization suggests a general dictionary between reflecting confining boundaries and internal interfaces for any observable with a derivative jump, which could be tested in soft-potential or trap experiments beyond the exactly solvable one-gap model.","A natural stress test is the predicted conservation law $B_k\\sim k^4$ as $k\\to 0$; deviations would signal either a non-conserved density response or contact-geometry contributions beyond regular faces.","The invariance result implies that any collective closure respecting the boundary-flux pairing is automatically consistent with the exact short-time coefficient, so a discrepancy in $B_k$ would point to the short-time dynamics itself rather than the later closure."],"forward_implications":["Experimenters can compute $B_k$ entirely from independently measured equilibrium and transport inputs---$S(k)$, $g(\\sigma^+)$, and $D_0$---and predict the short-time plateau of $[\\phi_k(t)-1+\\Gamma_k t]/t^{3/2}$ without any fitting.","Memory-kernel reconstructions and collective closure schemes gain an independently fixed short-time boundary condition, so later-time theory only has to model caging and structural escape.","The contact coefficient carries a universal wave-number dependence $G(x)=1-j_0(x)+2j_2(x)$, so different wave vectors provide multiple independent tests of the same physical inputs.","For soft interfaces, the leading fractional coefficient is exactly half the reflecting-boundary value with the corresponding one-sided slope, giving a separate calibration for penetrable particles or effective potentials.","Adding regular collective variables, however many, leaves the leading $s^{-1/2}$ contact term invariant; modes carrying their own surface distributions must be assigned to the contact sector."],"supporting_citations":[{"why":"Identifies the nonanalytic contact contribution to Brownian density correlations that the present $t^{3/2}$ coefficient generalizes.","marker":"[1]"},{"why":"Formulates the dynamic structure factor through a memory function, establishing the short-time branch the paper extends.","marker":"[2]"},{"why":"Provides explicit short-time hard-sphere results whose angular contact factor Eq. (12) reproduces.","marker":"[3]"},{"why":"Connects hard-sphere short-time corrections with colloidal diffusion and structure, supporting the contact-factor comparison.","marker":"[14]"},{"why":"Introduces the ensemble projection underlying the exact memory equation used throughout the paper.","marker":"[4]"},{"why":"Formulates generalized Langevin equations and memory kernels that the projection sector builds on.","marker":"[5]"},{"why":"Supplies the projection identity that the paper's weak form Schur complement extends to surface functionals.","marker":"[35]"},{"why":"Provides the exactly solvable one-gap model used to calibrate the soft-interface normalization.","marker":"[37]"}],"fun_headline_variants":["Contacts set the exact t^{3/2} decay in dense liquids","Fit-free law: collisions dictate density relaxation amplitude","Reflecting contacts force the t^{3/2} density signal","Hard-sphere contacts pin down fractional relaxation amplitude"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation assumes that within the thin diffusion layer where contacts act, the mobility, contact geometry, and density response are effectively frozen, with curvature and tangential variation entering only one order later; if that separation of scales is violated, the leading coefficient is approximate rather than exact.","fun_headline_variants_meta":{"raw":{"variants":["Contacts set the exact t^{3/2} decay in dense liquids","Fit-free law: collisions dictate density relaxation amplitude","Reflecting contacts force the t^{3/2} density signal","Hard-sphere contacts pin down fractional relaxation amplitude"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000285,"raw_usage":{"total_tokens":1747,"prompt_tokens":1084,"completion_tokens":663,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":700,"completion_tokens_details":{"reasoning_tokens":595}},"tokens_in":700,"tokens_out":663,"duration_ms":5487,"temperature":1.0,"reasoning_tokens":595,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T14:30:26.923522+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run overdamped Brownian dynamics of monodisperse hard spheres over a short-time window, independently measure $S(k)$, $g(\\sigma^+)$, and $D_0$, and plot $[\\phi_k(t)-1+\\Gamma_k t]/t^{3/2}$; if this ratio does not approach the predicted $B_k$ from Eq. (12) across a range of wave numbers, the claimed exactness fails.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the nonanalytic contact contribution to Brownian density correlations that the present $t^{3/2}$ coefficient generalizes."},{"cited_title":"Cichocki and W","cited_arxiv_id":null,"evidence_quote":"Formulates the dynamic structure factor through a memory function, establishing the short-time branch the paper extends."},{"cited_title":"Cichocki and B","cited_arxiv_id":null,"evidence_quote":"Provides explicit short-time hard-sphere results whose angular contact factor Eq. (12) reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Connects hard-sphere short-time corrections with colloidal diffusion and structure, supporting the contact-factor comparison."},{"cited_title":"Mori, Prog","cited_arxiv_id":null,"evidence_quote":"Formulates generalized Langevin equations and memory kernels that the projection sector builds on."}],"review_version":2}