{"id":"b6ca9a7b-386e-4a0b-9a47-69f7a19e4de6","arxiv_id":"1908.01484","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"The mean pore radius from Torquato's and Gubbins's pore size distributions can serve as the geometrical mesh size of polymer solutions across dilute to concentrated densities.","lead":"This paper shows that the average size of the pores between polymer chains, measured with two standard pore-size definitions, follows the same scaling laws as the mesh size of a polymer solution. This gives researchers a direct way to compute the mesh size and its full distribution from particle coordinates in simulations.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The identification <r> = xi is underdetermined: scaling fixes only the exponent, <r>_G and <r>_T differ by a factor near 2, and absolute values depend on the arbitrary monomer radius r_m. At most one estimator can equal xi.","rationale":"The paper does several things well: the PSD computation is clearly described, the two definitions are correctly contrasted, the model gel in Fig. 12 provides a direct test that favors PG over PT, and the OC expression in Eq. (46) is a genuine parameter-free prediction for <r>_T that works over a wide range despite the poor shape agreement in Fig. 16. My concern is not that the method is useless; it is that the headline identification 'xi can be identified with either <r>_T or <r>_G' is stronger than the evidence. Scaling collapse fixes only the exponent and is insensitive to O(1) prefactors, while the two means differ by a factor of about 2. At most one of them can be the actual mesh size in absolute terms, and the text ultimately selects Gubbins on the basis of a single idealized gel. The r_m dependence noted by the reader is part of the same underdetermination: because the pore space is defined through a hard-sphere cutoff for a soft WCA potential, absolute pore radii are convention-dependent. This does not invalidate the qualitative and methodological contributions, but it should be stated as a calibration problem rather than as an unambiguous identification, which is what the conditional verdict captures.","tokens_in":32908,"tokens_out":10167,"duration_ms":117017,"concrete_test":"Run a controlled calibration on model networks with known mesh size h: simulate diamond or cubic-lattice networks with strand lengths h = 5, 10, 20, 50 and monomer hard-core radius sigma/2 at several densities spanning the semidilute window; compute PT(r) and PG(r) with r_m = 0.4, 0.5, 0.6. Check whether <r>_G or <r>_T, with or without a fixed r_m-shift, reproduces h within 10% across all strand lengths and densities. If neither does, or if the correct estimator or prefactor changes with h, the identification xi = <r> in solutions cannot be taken as quantitative; if Gubbins recovers h robustly, the conditional concern reduces to a calibration statement.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The claim that xi can be identified with either <r>_T or <r>_G is underdetermined at the quantitative level. The comparison in Fig. 10 tests only the semidilute exponent: Eq. (19)-(20) fix the prefactor via the approximate overlap density rho*, but the scaling line is an estimate, not a ground truth, and any length with the same rho-dependence would collapse. The two estimators are not interchangeable: Fig. 11b gives <r>_G/<r>_T ~ 2, so at most one of them can equal the same xi; the paper itself later concludes that PG, and hence <r>_G, is the better pore-size indicator, using the model gel in Fig. 12 as the only case with an independently known mesh. The absolute values additionally depend on the arbitrary hard-sphere monomer radius r_m = sigma/2 (Eq. (34) and footnote [66]); when <r> is comparable to r_m, as it is near the operational rho** = 0.3, a change of r_m by 0.1 shifts <r> by a non-negligible fraction. The parameter-free OC model (Eq. (46)) supplies an absolute prefactor, but its large-eta limit is built from the scaling of R_g and eta_c (Eqs. (47)-(51)), so it does not independently certify that either mean is the true mesh size. Thus the evidence supports 'same scaling behavior' but not the stronger quantitative identification, especially with both estimators simultaneously.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper uses molecular dynamics simulations of coarse-grained Kremer-Grest polymer chains (N = 50, 200, 1000; monomer densities from 0.001 to 1.0) to characterize the geometrical mesh size ξ of polymer solutions. It measures the monomer density fluctuation correlation length ξ_c from both the static structure factor and the radial distribution function, and computes pore size distributions (PSDs) using the definitions of Torquato and of Gubbins. The authors report that the mean pore sizes ⟨r⟩_T and ⟨r⟩_G both follow the scaling form predicted for ξ in the semidilute regime, and argue that ξ can be identified with either quantity. They also introduce a parameter-free 'overlapping chains' (OC) model, in which the polymer solution is mapped onto overlapping spheres with radius R_g and density ρ/N, and show that the resulting analytical mean ⟨r⟩_T^{OC} reproduces ⟨r⟩_T accurately for ρ < ρ**. A model cubic polymer gel is used as an independent test to argue that Gubbins's PSD is the more direct indicator of the true pore size.","tokens_in":33251,"tokens_out":4732,"duration_ms":52614,"significance":"If the quantitative identification ⟨r⟩ = ξ holds, the paper would provide a direct, distribution-level route to the mesh size, which is more informative than scaling estimates or correlation-length measurements and is relevant for models of nanoparticle diffusion in polymer solutions. The simulation work is careful and standard, and the paper is strong on several counts: it compares two PSD definitions on the same systems, validates Gubbins's PSD against a model gel with independently known mesh size, and derives an elegant parameter-free approximation for ⟨r⟩_T from R_g and ρ. The scaling-exponent agreement in Figs. 10 and 11 is convincing. However, the quantitative identification of both ⟨r⟩_T and ⟨r⟩_G with the same ξ is not established: the two means differ by a factor near 2, the absolute values depend on the arbitrary monomer radius r_m, and the scaling comparison fixes only an exponent with an approximate prefactor.","major_comments":[{"comment":"The central claim that ξ can be identified with either ⟨r⟩_T or ⟨r⟩_G is not supported by the evidence presented. The comparison in Fig. 10 tests essentially the semidilute exponent: Eq. (19) carries an unknown O(1) prefactor and uses the approximate overlap density ρ* from Eq. (20), so any length with the same χ-dependence would collapse onto the same line. Fig. 11b shows ⟨r⟩_G/⟨r⟩_T ≈ 2, so at most one of the two means can equal the same numerical ξ. The model-gel test in Fig. 12 independently favors PG rather than PT. The conclusion should therefore be restricted to proportionality (⟨r⟩_T ≈ c_T ξ, ⟨r⟩_G ≈ c_G ξ with c_G/c_T ≈ 2), unless a calibration against an independently known mesh size is supplied.","section":"IV C, Figs. 10 and 11"},{"comment":"The absolute values of ⟨r⟩_T and ⟨r⟩_G depend on the arbitrary hard-sphere monomer radius r_m = σ/2 used to define the pore-solid interface. The authors acknowledge this in footnote [66], but they do not quantify the sensitivity. Near the operational concentrated-regime onset ρ** = 0.3 the mean pore sizes become comparable to the monomer size (for N = 50, R_g0 ≈ 4.8, so scaling values of ⟨r⟩ are of order unity), and a change of r_m by 0.1σ would shift ⟨r⟩ by a non-negligible fraction. Since the literal identification ⟨r⟩ = ξ is central to the paper, a sensitivity study over r_m is required, at least for the densities where ⟨r⟩ ≈ σ.","section":"IV B, Eq. (34) and footnote [66]"},{"comment":"The correlation-length data in Fig. 5 are obtained from fits that the authors themselves describe as not satisfactory in the intermediate density range. In Sec. IV A they state that for densities 0.11 < ρ < 0.64 neither the low-density form Eq. (9) nor the high-density form Eq. (10) gives a really satisfactory description, and that they nevertheless fit with an exponential 'to have at least an estimate'. These ξ_c values are then used in Fig. 10 in comparisons that support the paper's negative conclusions about ξ = ξ_c outside the semidilute regime. The paper should either restrict the ξ_c analysis to the ranges where the functional forms are justified or report fit quality measures that allow the reader to judge the uncertainty.","section":"IV A, Fig. 4"},{"comment":"The OC model is genuinely parameter-free, but the agreement shown in Fig. 17 does not independently certify that ⟨r⟩_T equals ξ. Equation (51) demonstrates that, for long chains, the OC expression reproduces exactly the scaling form of ξ in Eq. (18), because the inputs R_g and η_c are themselves chosen from the same scaling theory. The agreement is therefore partly a consistency check of the scaling ansatz rather than an external validation of the absolute pore size. The text should be revised to present the OC result as a practical estimator of ⟨r⟩_T, not as evidence for the literal identification ⟨r⟩_T = ξ.","section":"IV D, Eqs. (46)-(51), Fig. 17"}],"minor_comments":[{"comment":"The caption contains a typo: 'Torquato's PDS' should read 'Torquato's PSD'.","section":"Fig. 13 caption"},{"comment":"The caption reads 'contineous lines'; this should be 'continuous lines'.","section":"Fig. 7 caption"},{"comment":"No statistical uncertainties are reported for ⟨r⟩_T, ⟨r⟩_G, or ξ_c. Given that the paper makes quantitative claims about 'high accuracy', error bars or at least a statement of the statistical error from block averaging should be added.","section":"Figs. 5, 10, 11, 17"},{"comment":"The logarithmic fit line is quoted as 2.18 + 0.0458 ln(χ), but the fit range, the uncertainty of the coefficients, and the goodness of fit are not given. This should be specified.","section":"Fig. 11b"},{"comment":"The derivation of the Ornstein-Zernike form from Eq. (11) is only sketched in footnote [56]; a brief sentence in the main text would make the step clearer.","section":"II B, Eq. (12)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this paper gives you a practical way to measure mesh size distributions in polymer simulations, not just an average. Second: the quantitative identification of the mean pore size with xi is weaker than the abstract claims, but that doesn't sink the method.\n\nThe genuinely new piece is the application of two standard porous-media pore size distributions (Torquato, Gubbins) to polymer solutions, and the overlapping-chains (OC) estimator, Eq. (46), which reproduces <r>_T from Rg and density with no fit parameters. I checked the logic: it uses the analytical Torquato PSD for overlapping spheres and plugs in Rg and chain density. That it works across the semidilute range is a real result. The gel test in Fig. 12 is also a good idea—a system with a known mesh—and it supports the earlier conclusion that Gubbins's PG is the more interpretable distribution. The simulations are standard Kremer-Grest, with proper equilibration checks and several chain lengths; the scaling collapse in Fig. 10 is convincing.\n\nNow the soft spots. The paper's abstract says xi can be identified with either <r>_T or <r>_G, but these two differ by a factor near 2 (Fig. 11b), so at most one can equal xi. The scaling comparison in Fig. 10 fixes only the exponent; any length with the same rho-dependence collapses onto that line. The absolute values also depend on the hard-sphere monomer radius r_m = sigma/2, and near rho** the pores are not much larger than r_m, so the shift is not negligible. The paper does flag this in footnote 66, but the abstract overstates. There are also no error bars on the reported <r> values, and the xi_c fits in the intermediate density range are made with a functional form the authors say is not justified—they are honest about it, but it is a weak point. None of this kills the method; it does mean the quantitative claim should be softened to \"both estimators have the same scaling as xi, and PG gives the better absolute estimate,\" with the gel simulation as the real benchmark.\n\nWho is it for: anyone simulating or interpreting diffusion of nanoparticles or proteins in polymer solutions, and anyone who needs a mesh size distribution rather than a single correlation length. The citation pattern looks appropriate. It deserves a serious referee; I would send it out with a request for error bars and a careful rewrite of the identification claim.","headline":"A genuinely useful methods paper: PSD-based mesh size distributions plus a working parameter-free estimator, though the abstract overstates the quantitative identification of <r> with xi.","tokens_in":33768,"tokens_out":2448,"would_cite":true,"duration_ms":25888,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["61.25.Hq"],"model":"deepseek-v4-flash","headline":"The paper argues that the mean values of the two standard pore size distributions both follow the scaling predicted for the polymer mesh size, so the mesh size can be identified with the average pore size, and a parameter-free formula…","keywords":["polymer solutions","mesh size","pore size distribution","scaling theory","correlation length","coarse-grained simulation","radius of gyration","semidilute regime"],"falsifier":"Compute the two pore-size distributions at a density just below the concentrated-regime onset, $\\rho\\approx 0.3$, using two monomer radii, e.g. $r_m=0.4\\sigma$ and $r_m=0.6\\sigma$. If the resulting means $\\langle r\\rangle_T$ or $\\langle r\\rangle_G$ shift by more than the scatter that scaling theory permits, or if $\\langle r\\rangle_G/\\langle r\\rangle_T$ departs significantly from $\\approx 2$, then the claimed identification of the mesh size with a pore-size mean depends too strongly on the hard-sphere monomer model to be the geometric mesh size.","tokens_in":32714,"feed_emoji":"🫧","tokens_out":7802,"duration_ms":73314,"temperature":0.7,"pith_summary":"This paper tries to establish a direct geometric way to measure the \"mesh size\" of a polymer solution—the typical gap through which dissolved particles must move. Simulating coarse-grained chains from dilute to concentrated densities, the authors compute two standard pore-size distributions and show that their mean values both follow the scaling behavior predicted for the mesh size by polymer scaling theory. This allows the mesh size to be identified with a directly measurable quantity, and gives access to the full distribution of mesh sizes rather than only an average. In the semidilute regime, a parameter-free \"overlapping chains\" formula, using only the radius of gyration and the monomer density, reproduces the Torquato average pore size to high accuracy.","feed_headline":"Mean pore size follows polymer mesh scaling across densities","feed_subtitle":"The mesh that gates nanoparticle diffusion can be read from pore sizes using only Rg and density.","key_machinery":"The load-bearing object is the pore size distribution (PSD), computed by two definitions. Torquato's PSD, $P_T(r)$, is the probability density for a random solvent point to lie at distance $r$ from the nearest point on the monomer surface; it is obtained from the fraction of pore volume accessible to the center of a hard probe sphere of radius $r$. Gubbins's PSD, $P_G(r)$, instead assigns to $r$ the volume coverable by spheres of radius $r$ or smaller (the Connolly/reentrant surface), so its mean gives the pore size directly—for an isolated spherical pore $\\langle r\\rangle_G$ equals the pore radius while $\\langle r\\rangle_T=R/4$. The argument is carried by the analytic PSD of overlapping spheres, whose mean yields the parameter-free formula Eq. (46); the Puiseux asymptotics of the exponential integral convert this formula into the scaling powers $\\chi^{-\\nu+1/3}$ and $\\chi^{-\\nu}$.","core_discovery":"The central claim is that the geometrical mesh size $\\xi$ of a polymer solution equals the mean pore size defined by either of the two standard pore-size distributions. Torquato's mean $\\langle r\\rangle_T$ (distance from a random solvent point to the nearest monomer surface) and Gubbins's mean $\\langle r\\rangle_G$ (radius of the largest sphere coverable at that point) both collapse onto the scaling prediction $\\xi/R_{g0}\\propto \\chi^{-\\nu}$ in the semidilute regime and onto the dilute-limit behavior $\\chi^{-\\nu+1/3}$, while the density-fluctuation correlation length $\\xi_c$ follows this only in the semidilute regime and rises again at high density. The identification is supported by a model polymer gel with known mesh size, where Gubbins's distribution peaks at the true mesh size. The paper further claims that the parameter-free overlapping-chains expression $\\langle r\\rangle_T^{\\mathrm{OC}}=(R_g/3)e^{\\eta_c}E_{2/3}(\\eta_c)$, with $\\eta_c=4\\pi R_g^3\\rho/(3N)$, reproduces the simulated $\\langle r\\rangle_T$ in the semidilute regime and asymptotically yields the scaling law for long chains.","pith_inferences":["Because Gubbins's mean gives the true pore radius in idealized pores, $\\langle r\\rangle_G$ is arguably the better proxy for the particle radius separating free from obstructed diffusion; the systematic factor $\\langle r\\rangle_G/\\langle r\\rangle_T\\approx 2$ should be checked against tracer-diffusion measurements in the same systems.","The overlapping-chains formula implies that the mean pore size is controlled mainly by the chain-level packing variable $R_g^3\\rho/N$; in real systems this could let one map measured radii of gyration and concentrations directly onto mesh sizes without any scattering-derived correlation length.","A direct numerical test of the proposed identification would be to compare $\\langle r\\rangle_G$ with the tracer radius at which the diffusivity crossover occurs in the same simulated solutions; the paper does not perform this comparison."],"forward_implications":["The mesh size becomes a directly measurable geometric quantity at any density, with no unknown multiplicative prefactor and no dependence on an assumed overlap concentration.","The full distribution of mesh sizes, not just its mean, is accessible; this enables direct study of how spatial heterogeneity of the mesh affects tracer-particle diffusion.","In semidilute solutions, $\\langle r\\rangle_T$ (and hence, up to the factor $\\approx 2$, $\\langle r\\rangle_G$) can be estimated without simulation or scattering by inserting the radius of gyration and density into Eq. (46).","Using $\\xi_c$ as a proxy for the mesh size is ruled out except in the semidilute regime, since $\\xi_c$ increases with density in concentrated solutions while pore sizes decrease.","The PSD method extends to other polymer topologies and to networks and gels, where the mesh-size distribution is currently hard to access."],"supporting_citations":[{"why":"defines Torquato's pore size distribution and its formalism for porous media","marker":"[1]"},{"why":"defines Gubbins's pore size distribution via volumes coverable by spheres of radius r or smaller","marker":"[2]"},{"why":"supplies the blob model and scaling concepts underlying the predicted mesh-size behavior","marker":"[3]"},{"why":"provides the scaling expressions for correlation length, radius of gyration, and semidilute/concentrated regimes","marker":"[4]"},{"why":"derives the analytic pore-size distributions for hard and overlapping spheres used for the mapping","marker":"[38]"},{"why":"gives the exact hard-sphere PSD expression with density-dependent coefficients used in fits","marker":"[39]"},{"why":"provides the computational algorithm and further analytic results used for Torquato's PSD","marker":"[40]"},{"why":"supplies the optimization-based algorithm used to compute Gubbins's PSD","marker":"[43]"},{"why":"provides the bead-spring polymer model used in the molecular dynamics simulations","marker":"[57]"}],"fun_headline_variants":["Pore size reveals polymer mesh across densities","Mesh size from pore distributions: scaling matched","Polymer mesh equals mean pore size, no fit needed","Pore size distributions pin down mesh scaling","Rg and density predict mesh size in semidilute"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The entire pore geometry is built by treating each monomer as a hard sphere of radius $\\sigma/2$; whenever the average pore size approaches that monomer radius, the computed means—and hence the identification $\\langle r\\rangle = \\xi$—depend on this modeling choice.","fun_headline_variants_meta":{"raw":{"variants":["Pore size reveals polymer mesh across densities","Mesh size from pore distributions: scaling matched","Polymer mesh equals mean pore size, no fit needed","Pore size distributions pin down mesh scaling","Rg and density predict mesh size in semidilute"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000172,"raw_usage":{"total_tokens":1353,"prompt_tokens":1098,"completion_tokens":255,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":714,"completion_tokens_details":{"reasoning_tokens":196}},"tokens_in":714,"tokens_out":255,"duration_ms":3060,"temperature":1.0,"reasoning_tokens":196,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T15:10:25.882429+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the two pore-size distributions at a density just below the concentrated-regime onset, $\\rho\\approx 0.3$, using two monomer radii, e.g. $r_m=0.4\\sigma$ and $r_m=0.6\\sigma$. If the resulting means $\\langle r\\rangle_T$ or $\\langle r\\rangle_G$ shift by more than the scatter that scaling theory permits, or if $\\langle r\\rangle_G/\\langle r\\rangle_T$ departs significantly from $\\approx 2$, then the claimed identification of the mesh size with a pore-size mean depends too strongly on the hard-sphere monomer model to be the geometric mesh size.","supporting_citations":[{"cited_title":"Koshy, T","cited_arxiv_id":null,"evidence_quote":"derives the analytic pore-size distributions for hard and overlapping spheres used for the mapping"},{"cited_title":"Torquato, Phys","cited_arxiv_id":null,"evidence_quote":"gives the exact hard-sphere PSD expression with density-dependent coefficients used in fits"},{"cited_title":"Prager, Chem","cited_arxiv_id":null,"evidence_quote":"provides the computational algorithm and further analytic results used for Torquato's PSD"},{"cited_title":"Torquato, Random heterogeneous materials: mi- crostructure and macroscopic properties (Springer Sci- ence & Business Media, 2013)","cited_arxiv_id":null,"evidence_quote":"supplies the optimization-based algorithm used to compute Gubbins's PSD"}],"review_version":1}