{"id":"7de8dfd4-b1b4-4fc6-9be4-29d4fe1a2dac","arxiv_id":"2608.07195","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"A seven-parameter free energy formula, fitted to Monte Carlo data, reproduces the radius-of-gyration distribution of 3D interacting self-avoiding walks across the coil-globule crossover, except near the theta point.","lead":"This paper proposes a formula for the full distribution of a polymer's size, the radius of gyration, for interacting polymer chains of any length in good or poor solvent. If the formula holds, it gives researchers a fast way to interpret single-molecule or chromatin imaging data, not just the average size.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The central free-energy formula inherits its t^2 dependence from the unverified extensivity/analyticity ansatz for all contact cumulants (Eq. 37); the paper's own Θ-point consistency note (n≈1.53) and undefined exponent q leave Eq. 56 as an in-sample fit rather than an explicit derived free energy.","rationale":"The reader's weakest assumption is also the load-bearing one for the central claim. A seven-parameter function of ε can fit a wide family of histograms, so in-sample agreement with the very data used for Bayesian inference is weak evidence that the specific exponents in Eq. 56 are correct. The paper is commendably honest about the near-Θ deviations and about the n≈1.53 inconsistency, but those admissions support the concern rather than resolving it. The stretched-chain term (N t)^(-q) in Eq. 51 is non-analytic in t and is not accompanied by a defined value of q; re-derivation from Sec. III gives q=(3ν-1)/(3(1-ν))≈0.618, but as printed the formula is incomplete. This does not change the appropriate verdict: the paper deserves a CONDITIONAL status pending an independent test of Eq. 37 or an out-of-sample validation, not rejection. The proposed conditional-cumulant measurement would settle whether the assumed analytic t-expansion is real or merely an artifact of the fitting procedure.","tokens_in":19354,"tokens_out":15267,"duration_ms":150574,"concrete_test":"Use the paper's ε=0 SAW simulation code to generate contact-number statistics conditioned on R for at least N=512, 1024, 2048, and 4096. For each N and R bin, compute κ_1, κ_2, κ_3 from the conditional cumulant-generating function K_ε(R)=ln⟨e^{εm}⟩_{0,N,R}; then plot κ_n(t)/N against t=(N/R^3)^g with g=1/(3ν-1)≈1.309. Check three things: (i) collapse of κ_n/N for different N (extensivity); (ii) after subtracting the constant, a leading term linear in t; (iii) the next term scaling as t^2 rather than t^p with p≈1.53. If any of these fail, Eq. 37 fails and Eq. 56 loses its derivation; if they hold, the formula's basis is substantially strengthened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. 56 is the central claim, but its functional form is not derived; it is the direct product of Eq. 37, the assertion that all fixed-R contact cumulants are extensive and analytic, κ_n(t)=N(a_n+B_n t+B'_n t^2+o(t^2)), for every n. That assertion is never tested directly. The only evidence offered is the final seven-parameter Bayesian fit to the same MC histograms used for comparison, a procedure that can compensate for a wrong t-dependence by adjusting a_i and c_i, and hence cannot validate the assumed exponents. The paper itself flags the inconsistency at the Θ-point: consistency with ν_θ=1/2 would replace the t^2 term by t^n with n≈1.53, so the t^2 ansatz cannot be universal. In addition, the SAW entropy contains a non-analytic stretched-chain term (N t)^(-q) (Eq. 51); q is never defined in the text, and no analogous interaction-dependent singular term appears in Eq. 56. Until Eq. 37 is checked against conditional contact-number statistics, the N/ε factorization and the identification of a_1(ε) with the second virial coefficient remain phenomenological fits, not an explicit validated free energy.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper claims to derive an explicit finite-size free energy for a three-dimensional interacting self-avoiding walk (ISAW) as a function of the radius of gyration, chain length N, and interaction energy ε. The central result is Eq. (56), βF_N(t|ε) = a1(ε) N t + a2(ε) N t^2 + a3(ε) (N t)^{-q} + a4(ε) (N t^2)^{2/3} + c_N(ε) ln N t, with t a renormalized density and c_N(ε) given by Eq. (53). The derivation proceeds through scaling arguments for the SAW free energy, a cumulant expansion of the contact-number distribution (Section IV), and ad hoc additions of logarithmic, stretched-chain, and surface terms (Section V). The seven ε-dependent coefficients are inferred by Bayesian MCMC from Monte Carlo simulations of lattice ISAWs, and the resulting distributions are compared with the same simulated histograms. The paper also uses the fitted formula to construct a finite-size coil-globule phase diagram.","tokens_in":19756,"tokens_out":3912,"duration_ms":39190,"significance":"If the central formula were established as a genuine free-energy expression, it would provide a compact, transferable description of the entire radius-of-gyration distribution across coil and globule regimes, with practical value for interpreting single-molecule and super-resolution experiments. The manuscript's strengths are its extensive lattice simulations (chain lengths up to ~23,000, many ε values), the transparency of the Bayesian inference procedure, and the provided reproducible code. However, the claimed derivation rests on unverified assumptions about the cumulants, the functional form is partly the outcome of trial-and-error model selection, and the validation is in-sample. These issues currently limit the result to a well-tested phenomenological fit rather than the explicit derived free energy advertised in the abstract and conclusions.","major_comments":[{"comment":"The central step of the derivation is the assertion that every contact-number cumulant κ_n(t) is extensive and has the analytic asymptotic form N(a_n + B_n t + B'_n t^2 + o(t^2)). This is stated for all orders n without proof or direct numerical check. Since Eqs. (41) and (56) inherit their t and t^2 dependence from this assumption, the claimed generality of the free-energy formula is not established. The paper itself notes in Section VI.D that consistency with the known Θ-point exponent ν_θ = 1/2 would replace the t^2 term by t^n with n ≈ 1.53, explicitly acknowledging that the exponent 2 is not universal. The authors should either derive the cumulant expansion from a controlled approximation, or directly measure κ_n(t) from the simulations at fixed R and test Eq. (37).","section":"Section IV.B, Eq. (37)"},{"comment":"The logarithmic term c_N(ε) ln N t, the stretched-chain term (N t)^{-q}, and the surface term (N t^2)^{2/3} are introduced through \"careful trial and error\" and \"several trials\" after poor fits, rather than derived. The exponent q is never defined, estimated, or reported, so Eq. (56) is not an explicit closed-form free energy in the sense claimed. At minimum, q must be specified and the derivation of the surface term's exponent should be made independent of the fitting data, or the manuscript should be reframed as presenting a phenomenological free-energy form.","section":"Section V.A–V.B, Eqs. (51)–(56)"},{"comment":"The validation is circular. The parameter vector θ(ε) = (a1, a2, a3, a4, c0, c1, c2) is obtained by maximizing the likelihood in Eq. (57) over the same {R^2_i,N} histograms that are later compared with the theory in Figures 2–5. With seven parameters per ε and additional piecewise polynomial regressions in ε (Section VI.C), the agreement measures the quality of the fit, not the predictive validity of Eq. (56). To support the abstract's claim that this \"validates our resulting free energy expression,\" the authors should perform a held-out validation: fit θ(ε) on one subset of the (N, ε) grid and compare predictions on the complementary subset, or use k-fold cross-validation.","section":"Section VI.B–VI.D"},{"comment":"The authors report that a4 unreasonably increases as ε → 0 and attribute this to \"overfitting in these conditions.\" This admission indicates that the fitted coefficients can compensate for missing or misspecified terms in Eq. (56). Consequently, the physical interpretation of a1(ε) as the second virial coefficient and of other coefficients as specific free-energy contributions is not uniquely supported: equally good fits might be obtained with different term balances. The claim that the parameters have \"a clear interpretation\" (Section VII) therefore needs additional support, for example by checking the stability of each a_i(ε) when terms are removed one at a time.","section":"Section VI.C"}],"minor_comments":[{"comment":"The exponent q in the term (N t)^{-q} is left undefined in the text and never estimated; please state its value, or at minimum how it is determined, since it is part of the central formula.","section":"Section V.A, Eq. (51)"},{"comment":"The phrase \"infection point\" should be \"inflection point\".","section":"Section VI.E"},{"comment":"Reference [21] gives the publisher of de Gennes' book as \"Cornel University Press\"; this should be \"Cornell University Press\".","section":"References"},{"comment":"The explanation that the discrepancy at ε = 0.28 is \"numerical rather than theoretical\" is speculative and not tested (e.g., by longer simulations or different equilibration protocols); since the parameters are fitted to the same data, the deviation could equally indicate missing terms in Eq. (56).","section":"Section VI.D"},{"comment":"The matching conditions are presented as if they determine the free-energy exponents uniquely, but the derivation assumes power-law forms for S, G, and C without justifying why sums of only two power laws are sufficient; a sentence acknowledging this modeling assumption would improve clarity.","section":"Section III, Eqs. (8)–(11)"}],"recommendation":"major_revision","confidential_remarks":"The central issue is that the paper's headline claim—an explicit, validated free-energy expression—is undercut by the in-sample fitting procedure and by the paper's own admission that the t^2 exponent is inconsistent with ν_θ = 1/2. This is fixable if the manuscript is reframed as a phenomenological model with held-out validation and with the cumulant ansatz tested directly. The authors are not claiming anything outside consensus; the problem is internal consistency between the derivation narrative and what is actually demonstrated. The paper may be better received if the abstract and conclusions are moderated and the code is emphasized as a practical tool."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my take. This paper gives an explicit finite-size free energy for the radius-of-gyration distribution of a 3D ISAW, Eq. (56), with a new surface term (N t^2)^(2/3) and a modified logarithmic correction c_N(ε) ln N t. It builds on the Lhuillier–Victor–Imbert framework and extends it to the interacting case. The fits to the full R^2 histograms across the coil–globule crossover, for N up to ~23000 and a wide ε range, are genuinely good, and the code for the interpolation functions is released. That is real work and a useful tool for fitting single-molecule and chromatin data.\n\nWhere it is soft: the 'validation' is in-sample. The seven parameters θ(ε) are fitted to the same Monte Carlo histograms that are later shown as agreement. That does not invalidate the formula as a fit, but it means the agreement is not independent confirmation. The functional form itself was assembled partly by trial and error, which the paper admits. More specifically, the cumulant expansion rests on Eq. (37), the assertion that all fixed-R contact-number cumulants are extensive and analytic as κ_n(t)=N(a_n+B_n t+B'_n t^2+o(t^2)) for every n. That assertion is never tested directly and is not obvious. The paper itself flags the related inconsistency at the Θ-point: consistency with ν_θ=1/2 would replace t^2 by t^n with n≈1.53. And the exponent q in the stretched-chain term (N t)^(-q) is never defined in the text. These are real gaps, but they are stated honestly, and the near-theta deviation is acknowledged rather than hidden.\n\nWeighing it up: the central claim holds as a phenomenological model, not as a first-principles derivation. The reader's conditional verdict is about right. The paper deserves a serious referee, because the question matters, the formula is explicit and testable, and code and data are available. The referee should ask for out-of-sample cross-validation (e.g., fit on a subset of N, predict the rest), explicit parameter values with uncertainties for each ε, and a clear statement of what q is and what evidence supports Eq. (37). If those are addressed, it becomes a genuinely useful reference for polymer and chromatin experiments. If not, it remains an over-claimed fit.\n\nRecommendation: send to peer review, but with the expectation of significant revision.","headline":"A useful, honest phenomenological fit of the ISAW radius-of-gyration distribution, but the central formula is not derived and its validation is in-sample; worth refereeing with demands for out-of-sample tests.","tokens_in":20275,"tokens_out":2476,"would_cite":true,"duration_ms":22351,"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":"The paper derives an explicit free energy as a function of the radius of gyration for an interacting self-avoiding walk, showing that the resulting distribution matches Monte Carlo histograms from coil to globule except at the immediate…","keywords":["interacting self-avoiding walk","radius of gyration","coil-globule transition","finite-size scaling","cumulant expansion","renormalized density","Monte Carlo simulation","Bayesian inference"],"falsifier":"Measure, from the same fixed-$R$ contact-count histograms used in the paper, the third cumulant $\\kappa_3(t)$ for a range of $N$ and $\\varepsilon$ and test whether it follows $N(a_3+B_3 t+B_3' t^2)$ across the crossover. If it shows non-analytic behaviour in $t$ or an extra $N$-dependence, the general form Eq. (56) is falsified. A second, direct test is to fit the exponent of the quadratic term in $t$ from histograms taken in a narrow band around $\\varepsilon_\\Theta$; consistency with $\\nu_\\theta=1/2$ would demand an effective exponent of about 1.53 instead of 2, settling whether the $t^2$ term is truly universal.","tokens_in":19084,"feed_emoji":"📏","tokens_out":11887,"duration_ms":96468,"temperature":0.7,"pith_summary":"The paper's aim is an explicit, closed-form free energy for a single interacting self-avoiding walk as a function of its instantaneous radius of gyration $R$, chain length $N$, and monomer-monomer interaction $\\varepsilon$. Exponentiating that free energy gives the full probability distribution $P(R|\\varepsilon,N)$ of the radius of gyration. The proposed form uses one renormalized density variable $t = (N/R^d)^{1/(\\nu d-1)}$, combines entropic, surface-energy and logarithmic finite-size terms, and depends on $N$ and $\\varepsilon$ through seven coefficients that are functions of $\\varepsilon$ only. The coefficients are fixed by Bayesian inference on extensive Monte Carlo histograms, and the resulting distributions agree with simulation across the coil and globule regimes, with deviations only very close to the $\\Theta$-point. If correct, the formula makes the measured radius-of-gyration distribution a direct probe of solvent quality and of the coil-globule crossover.","feed_headline":"Explicit free energy formula predicts polymer radius-of-gyration distributions","feed_subtitle":"One seven-parameter free energy reproduces full R_g histograms in Monte Carlo, failing only at the theta point.","key_machinery":"The load-bearing object is the renormalized density $t = \\rho^{g}$ with $g=1/(\\nu d-1)$, built from the monomer density $\\rho=N/R^d$ and the Flory exponent $\\nu$. Matching the free energies of stretched, globular, and scale-invariant coil conformations fixes the exponents of the terms, and $t$ is the variable in which the low-density free energy is assumed analytic. The second piece is the cumulant expansion that expresses the interacting free energy as the non-interacting one minus the generating function of the contact-number cumulants $\\kappa_n(t)$; the assumption that each $\\kappa_n(t)=N(a_n+B_n t+B'_n t^2+o(t^2))$ is extensive and analytic in $t$ turns the infinite sum into the quadratic structure of Eq. (56). The surface term $(N t^2)^{(d-1)/d}$ is fixed by requiring both the crossover scaling $\\hat t=N^{1/2}t$ and the $N^{(d-1)/d}$ dependence of a globule surface.","core_discovery":"The central claim is that for a three-dimensional lattice ISAW the finite-size free energy as a function of the renormalized density $t=(N/R^3)^{1/(3\\nu-1)}$ is $$\\$\\beta$ F_N(t|\\varepsilon) = a_1(\\varepsilon) N t + a_2(\\varepsilon) N $t^{2}$ + a_3(\\varepsilon) (N t)^{-q} + a_4(\\varepsilon) (N $t^{2}$)^{2/3} + c_N(\\varepsilon) \\ln(N t),$$ with $c_N(\\varepsilon)=c_0(\\varepsilon)+c_1(\\varepsilon)/N+c_2(\\varepsilon)\\ln N$. The radius-of-gyration distribution then follows from $P(R|\\varepsilon,N) \\propto \\exp[-\\beta F_N(t|\\varepsilon)]$. The first two terms are the low-density virial-like expansion in $t$; the $(N t)^{-q}$ term and the logarithm encode entropic and surface counting corrections; the $(N t^2)^{2/3}$ term is a new surface-energy contribution chosen by scaling arguments. The paper validates the formula against Monte Carlo data for chain lengths from about 8 to 23,000 monomers and interaction strengths from 0 to 0.39, reporting agreement for the full histograms, means, and medians except in a narrow window at the immediate vicinity of the $\\Theta$-point, and uses it to build a finite-$N$ coil-globule phase diagram.","pith_inferences":["If Eq. (56) is correct, the distribution itself supplies a likelihood function for estimating $\\varepsilon$ from a single-molecule histogram of $R$, which is the only order parameter accessible in super-resolution imaging of chromatin or intrinsically disordered proteins.","The fitted $a_4(\\varepsilon)$ increasing toward $\\varepsilon=0$ hints at parameter compensation; extrapolating the formula to $\\varepsilon$ values outside the fitted range, or using the individual coefficients as physical virial coefficients, may therefore be unsafe.","The unresolved contradiction between the quadratic $t^2$ term and the required $\\nu_\\theta=1/2$ at the $\\Theta$-point suggests the true crossover variable may involve an effective exponent $n\\simeq 1.53$; simulations at larger $N$ just above $\\varepsilon_\\Theta$ could settle which form is fundamental.","The surface term $(N t^2)^{2/3}$ predicts a specific dependence on spatial dimension $d$ through $(N t^2)^{(d-1)/d}$; testing the formula in $d=2$ would discriminate it from earlier surface-energy forms."],"forward_implications":["For a lattice ISAW, the full finite-size distribution $P(R^2|\\varepsilon,N)$ is available in closed form once the seven $\\varepsilon$-dependent coefficients are known, without further simulation.","A two-dimensional finite-$N$ phase diagram can be constructed from the inflection of the mean radius as a function of $N$; the phase boundary is consistent with $\\varepsilon_N = a N^{-b} + \\varepsilon_\\Theta$ or with the classical $N^{-1/2}+N^{-1}$ form, with $\\varepsilon_\\Theta \\approx 0.27$.","The formula reproduces the mean, median, and entire histograms of $R^2$ for $N$ from roughly 8 to 23,000 and $\\varepsilon$ from 0 to 0.39, which the paper takes as validating the factorized $N$- and $\\varepsilon$-dependence.","Because the $\\varepsilon$-dependence is contained entirely in the coefficients, the same free energy can be used to compare different chain lengths at a given solvent quality without refitting.","The authors propose the expression as a basis for mapping off-lattice or coarse-grained polymer models onto the lattice ISAW description."],"supporting_citations":[{"why":"Starts the framework where the SAW free energy is written as $N$ times a power of density, whose matching fixes the shape of the final expression.","marker":"[13]"},{"why":"Establishes the three conformation classes (stretched, coil, globule) and the crossover scale $\\hat t = N^{1/2}t$ that the new surface term is built to satisfy.","marker":"[12]"},{"why":"Supplies the scaling exponents ($\\nu$, $\\gamma$) and the dilute-solution virial picture that justify analyticity in $t$.","marker":"[21]"},{"why":"Reference simulations of three-dimensional $\\Theta$ polymers, setting the expected $\\Theta$-point and the finite-size transition behaviour the model must reproduce.","marker":"[8]"},{"why":"Provides the fitted phase-boundary form $\\varepsilon_N = A N^{-1/2} + B N^{-1} + \\varepsilon_\\Theta$ used for comparison in the phase diagram.","marker":"[18]"},{"why":"Documents the lower density at the globule surface, the empirical motivation for including a distinct surface-energy term.","marker":"[9]"},{"why":"Classical source of the surface-energy form and of the $N^{-1/2}$ finite-size shift of the transition point that the paper extends.","marker":"[25]"},{"why":"The reptation-move Monte Carlo algorithm used to sample ISAW conformations in dense globular states.","marker":"[30]"},{"why":"The affine-invariant ensemble sampler used for the Bayesian inference of the seven coefficients.","marker":"[32]"}],"fun_headline_variants":["Seven-parameter free energy reproduces full polymer size distributions","Unified formula captures coil-globule transition and finite-size effects","Explicit finite-size free energy for interacting self-avoiding walks","Polymer radius of gyration distributions predicted by a single formula","From coils to globules: one free energy formula predicts polymer size"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation's load-bearing premise is that every cumulant of the contact-number distribution at fixed radius is extensive and analytic in the renormalized density $t$, with only linear and quadratic corrections to the leading $N$ term; anything more than that cannot be absorbed into the seven fitted coefficients without making Eq. (56) a fit rather than the claimed general expression.","fun_headline_variants_meta":{"raw":{"variants":["Seven-parameter free energy reproduces full polymer size distributions","Unified formula captures coil-globule transition and finite-size effects","Explicit finite-size free energy for interacting self-avoiding walks","Polymer radius of gyration distributions predicted by a single formula","From coils to globules: one free energy formula predicts polymer size"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000533,"raw_usage":{"total_tokens":2624,"prompt_tokens":1063,"completion_tokens":1561,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":679,"completion_tokens_details":{"reasoning_tokens":1473}},"tokens_in":679,"tokens_out":1561,"duration_ms":14618,"temperature":1.0,"reasoning_tokens":1473,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T13:05:38.975236+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure, from the same fixed-$R$ contact-count histograms used in the paper, the third cumulant $\\kappa_3(t)$ for a range of $N$ and $\\varepsilon$ and test whether it follows $N(a_3+B_3 t+B_3' t^2)$ across the crossover. If it shows non-analytic behaviour in $t$ or an extra $N$-dependence, the general form Eq. (56) is falsified. A second, direct test is to fit the exponent of the quadratic term in $t$ from histograms taken in a narrow band around $\\varepsilon_\\Theta$; consistency with $\\nu_\\theta=1/2$ would demand an effective exponent of about 1.53 instead of 2, settling whether the $t^2$ term is truly universal.","supporting_citations":[{"cited_title":"Lhuillier, A simple model for polymeric fractals in a good solvent and an improved version of the flory ap- proximation, Journal de Physique49, 705 (1988)","cited_arxiv_id":null,"evidence_quote":"Starts the framework where the SAW free energy is written as $N$ times a power of density, whose matching fixes the shape of the final expression."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the three conformation classes (stretched, coil, globule) and the crossover scale $\\hat t = N^{1/2}t$ that the new surface term is built to satisfy."},{"cited_title":"de Gennes,Scaling concepts in polymer physics (Cornel University Press, 1979)","cited_arxiv_id":null,"evidence_quote":"Supplies the scaling exponents ($\\nu$, $\\gamma$) and the dilute-solution virial picture that justify analyticity in $t$."},{"cited_title":"Grassberger and R","cited_arxiv_id":null,"evidence_quote":"Reference simulations of three-dimensional $\\Theta$ polymers, setting the expected $\\Theta$-point and the finite-size transition behaviour the model must reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Documents the lower density at the globule surface, the empirical motivation for including a distinct surface-energy term."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Classical source of the surface-energy form and of the $N^{-1/2}$ finite-size shift of the transition point that the paper extends."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The reptation-move Monte Carlo algorithm used to sample ISAW conformations in dense globular states."},{"cited_title":"Foreman-Mackey, D","cited_arxiv_id":null,"evidence_quote":"The affine-invariant ensemble sampler used for the Bayesian inference of the seven coefficients."}],"review_version":1}