{"id":"2bf6efeb-bbe3-40dd-ad62-6d1bf3097b62","arxiv_id":"2505.05254","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"For model PEG hydrogels, the equilibrium partition constant of linear polymer chains follows exp(-4(Rg/lcycle)^2), with lcycle the elastically effective network cycle length.","lead":"Researchers measured how polymer chains enter flexible hydrogels and found a simple law: the partition constant depends only on the squared ratio of chain size to the network's effective cycle length. The result gives a contactless method to measure partitioning and challenges the old rod-obstacle picture.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Most K0 values in Fig. 5 are single-point ratios measured above the stated dilute threshold, so the universal law is fit to concentration-dependent data.","rationale":"The reader's weakest_assumption concerned Eq. (2), the mixing osmotic pressure form, which is validated by direct drying on a subset. That is a real concern, but it is at least partially supported by the 8-point check in Fig. 2. The more load-bearing issue is that the K0 values used to establish Eq. (5) are, for most architectures, not true dilute-limit constants at all: they are single measurements at cext values the paper itself identifies as above the Nernst regime. This is not a matter of outside consensus; it is an internal inconsistency with the stated criterion in Fig. 3 and the text. For the one system with proper dilute data, the concentration dependence is visible (e.g., Msol=20 K rises from 0.03 to 0.08 between 20 and 30 g/L). If those biased points are removed, the universality claim rests on a very small fraction of the dataset, which would not justify a universal law. The reader did mention 'several K0 values taken above the stated dilute threshold' in the rationale, so there is partial agreement, but the formal 'weakest assumption' was placed elsewhere. The fix is straightforward: perform additional dilute-concentration measurements for the underrepresented architectures. Because the paper's central claim hinges on this, the CONDITIONAL verdict remains appropriate, but the required condition should explicitly include establishing genuine dilute-limit K0 values across the network parameter space.","tokens_in":14160,"tokens_out":9078,"duration_ms":87057,"concrete_test":"For each network architecture that currently has only one K0 point, measure additional dilute concentrations (cext/c*_ext < 0.8) to confirm a plateau in K vs cext, e.g., cext=20 and 40 g/L for Msol=5, 10 and 20 g/L for Msol=10, and 5 and 10 g/L for Msol=20. Then re-fit Eq. (5) using only well-defined dilute K0 values and check whether the collapse and the prefactor 4 remain. If fewer than five independent architectures retain a well-defined K0, or if the fit degrades or the prefactor shifts by more than a factor of two, the universal law is not supported by the current data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The universal law Eq. (5) is calibrated using K0 values from Tables S3 and S4. For all samples except the Mpre=20 kg/mol, cnet,0=60 g/L, f=4 series, each K0 is a single measured K at one external concentration: cext=80 g/L for Msol=5, 45 g/L for Msol=10, and 30 g/L for Msol=20. Table S1 gives c*_ext = 74.5, 43.9, and 25.8 g/L for these solutes, so these measurements fall at cext/c*_ext = 1.07, 1.03, and 1.16, all above the paper's own dilute limit. The text and Fig. 3 state that K~K0 only for cext/c*_ext < 1 and that above this threshold K > K0. The one system with multiple sub-c* points (Mpre=20, cnet,0=60) shows the bias is large for Msol=20: K rises from 0.03 at cext=20 g/L to 0.08 at 30 g/L and 0.16 at 60 g/L, so using the single 30 g/L point as K0 overestimates -lnK0 substantially. Since the Fig. 5 'various density' and 'various topology' datasets consist almost entirely of such single super-threshold points, the data collapse and fitted prefactor 4 do not demonstrate a dilute-limit universal law; they may instead represent a curve through concentration-dependent partition ratios.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports an experimental study of the equilibrium partitioning of linear poly(ethylene glycol) (PEG) chains into model PEG hydrogels with controlled network architecture. A label-free, osmotic-pressure-based method is introduced: the gel swelling ratio is used to infer the internal partitioned polymer concentration c_int through Eqs. (1) and (2), after validating Eq. (2) against an alternative additive decomposition by direct drying measurements. The paper then argues that the Nernst partition law holds in the dilute regime, that the conventional Ogston model is inadequate, and that the dilute partition constant follows the universal law K_0 = exp[-4(R_g/l_cycle)^2], where R_g is the solute gyration radius and l_cycle is the elastically effective cycle length computed from the Bethe approximation.","tokens_in":14541,"tokens_out":5552,"duration_ms":55640,"significance":"If the universal partition law were established, this would be a substantial contribution: it would replace the rod-based Ogston picture with a network-architecture-dependent length scale, and the osmotic method would be a useful label-free route to partition coefficients. The paper has notable strengths: the model networks are systematically varied in density, functionality, and connectivity; the central osmotic-pressure assumption Eq. (2) is directly cross-checked by drying for a subset of samples; and the collapse in Fig. 5 onto a single curve in R_g/l_cycle is a striking empirical observation. However, as detailed below, the evidence currently does not support the claimed dilute-limit universality: most K_0 values plotted in Fig. 5 are single-point measurements taken above the paper's own dilute threshold, the numerical prefactor and the choice of l_cycle are selected from the same data, and Eq. (2) is validated on only a subset of the architectures used in the universal-law fit.","major_comments":[{"comment":"Most K_0 values in Fig. 5 are not dilute-limit partition constants. For all samples except the M_pre = 20 kg/mol, c_net,0 = 60 g/L, f = 4 series, the plotted K_0 is a single measured K at c_ext = 80 g/L for M_sol = 5 kg/mol, c_ext = 45 g/L for M_sol = 10 kg/mol, or c_ext = 30 g/L for M_sol = 20 kg/mol. Table S1 gives c*_ext = 74.5, 43.9, and 25.8 g/L for these solutes, so these measurements are at c_ext/c*_ext = 1.07, 1.03, and 1.16, all above the dilute limit that the paper itself identifies. The text and Fig. 3 state that K > K_0 when c_ext/c*_ext > 1. The one system with multiple sub-c* points (M_pre = 20, c_net,0 = 60, f = 4, Table S3) shows a large bias for M_sol = 20: K rises from 0.03 at c_ext = 20 g/L (c/c* = 0.78) to 0.08 at c_ext = 30 g/L (c/c* = 1.16), so -ln K changes from 3.51 to 2.53. Using the single 30 g/L point as K_0 substantially overestimates -ln K_0 relative to the dilute value. Consequently, the 'various density' and 'various topology' datasets in Fig. 5 are largely concentration-dependent partition ratios, and the collapse and fitted prefactor 4 do not by themselves demonstrate a dilute-limit universal law.","section":"Section 'Nernst distribution law for polymer chain partitioning into polymer networks'; Tables S1, S3, S4"},{"comment":"The coefficient 4 in Eq. (5) is not derived from the theory; it is calibrated on the same data that are then shown as agreement, and l_cycle is selected among three candidate length scales (l_branch, l_cross, l_cycle) as the one that gives collapse in Fig. 5 while Fig. 6 shows that the other two fail. This is a model-selection procedure on the same dataset, not an independent test. The manuscript should state explicitly how the prefactor was obtained, report its uncertainty, and provide an out-of-sample test or a closed-form derivation; otherwise the 'universality' of Eq. (5) is a fitted empirical collapse with a selected axis rather than a validated prediction.","section":"Eq. (5) and Fig. 5; End Matter and Fig. 6"},{"comment":"The validation of Eq. (2) over the additive decomposition Eq. (3) is performed on only eight dried-gel samples (Table S5), all with f = 4, s = 0.5, c_net,0 = 60 g/L, and M_pre = 20 or 40 kg/mol. For all other densities and topologies in Fig. 5, c_int and hence K_0 are extracted from swelling measurements assuming Eq. (2) with the same prefactor A. If Eq. (3) or a different prefactor applied for f = 3, f = 8, s = 0.4, or c_net,0 = 30/90 g/L, every extracted K_0 would shift and the collapse could change. The authors should extend the direct drying check to at least one sample per topology class or provide a quantitative argument for transferability.","section":"Eq. (2), Fig. 2, and Table S5"},{"comment":"No replicate counts or error bars are reported for any K or K_0 value, and several entries are K_0 = 0 (e.g., M_pre = 10, M_sol = 20 rows in Table S3; M_pre = 40, f = 8, M_sol = 20 in Table S4). These zero values are excluded from the log-log inset without discussion, yet they are not small positive values as Eq. (5) would predict; they instead suggest a detection limit or a different regime. Without replicate measurements and a stated uncertainty propagation from the weights, osmotic pressures, and modulus data, the collapse in Fig. 5 cannot be distinguished from scatter of single-point measurements.","section":"Tables S3 and S4; Fig. 5"}],"minor_comments":[{"comment":"The term 'contactless' is used for a method that ultimately relies on weighing gels; the intended meaning appears to be 'label-free' or 'requiring no chemical modification.' Please clarify to avoid confusion.","section":"Abstract and introduction"},{"comment":"The solid lines for K_0 are drawn through data that include points near or above c*_ext; for M_sol = 20, the c_ext = 20 g/L point is at c_ext/c*_ext = 0.78, and the c_ext = 30 g/L point is above c*_ext. A statement of which points are included in each fit would make the determination of K_0 transparent.","section":"Fig. 3(b)"},{"comment":"Since the dilute threshold c*_ext is central to the interpretation, it would be helpful to state the experimental basis for Ψ* ≈ 0.24 and to mention how much c*_ext would vary with a reasonable uncertainty in Ψ*.","section":"Eq. (S2) and Table S1"},{"comment":"The table reports ν_0, μ_0, and ξ_0 without uncertainties; given that these feed directly into l_cycle and hence into the x-axis of Fig. 5, a sensitivity estimate would be useful.","section":"End Matter, Table I"}],"recommendation":"major_revision","confidential_remarks":"The manuscript comes from a group with a strong track record on model tetra-PEG networks, and the osmotic method is potentially valuable. My main concern is that the headline universal law is currently fitted to single-point, super-dilute-threshold measurements, and the prefactor and length-scale selection are made on the same dataset. These are correctable with additional dilute-regime measurements and a transparent fitting/validation protocol, so I do not see this as a reject; but the claims in the abstract and conclusions outpace the evidence as presented."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. The label-free osmotic method for measuring partition is genuinely new, and the validation of Eq. (2) by direct drying is careful and convincing. The scaling hypothesis, that K0 collapses on Rg/lcycle, is plausible and worth taking seriously. But the paper's central claim, the universal law with coefficient 4, is currently supported by data that mostly violate the paper's own dilute criterion. That is a load-bearing problem, not a cosmetic one.\n\nThe method: they infer cint from swelling against external PEG solutions, using Πmix = A(cnet + cint)^{3ν/(3ν−1)}. They check this against direct drying on a subset (Table S5, Fig. 2), and the agreement is good. They also try three network length scales and show that only lcycle collapses the data. That is a real observation if the data are trustworthy.\n\nThe soft spot: look at Tables S3 and S4. For the 'various density' and 'various topology' series, each K is a single measurement at cext = 80, 45, or 30 g/L for Msol = 5, 10, 20. Table S1 gives c* = 74.5, 43.9, 25.8 g/L. So nearly every point is at cext/c* > 1, where the paper itself says K > K0. For the one system with a full concentration sweep (Mpre = 20, cnet,0 = 60), K for Msol = 20 goes from 0.03 at 20 g/L to 0.08 at 30 g/L to 0.16 at 60 g/L. The single 30 g/L point used in Fig. 5 is roughly 2.7 times the dilute value. That overestimates −lnK0 and biases the fitted coefficient 4 upward. The collapse in Fig. 5 may therefore be a collapse of concentration-dependent K values, not of the dilute-limit constant.\n\nMinor related issues: no error bars or replicate counts anywhere, the coefficient is fit to the same data displayed, and lcycle is chosen among three candidates as the one that gives collapse. These would be survivable if the dilute-limit issue were fixed, but they compound the problem.\n\nThe paper is for people who work on gel mesh characterization, size-exclusion chromatography, and drug-delivery modeling. It deserves a serious referee: the method is novel and the scaling idea is interesting. But the universal law should not be accepted in this form. The authors need to remeasure K0 at genuinely dilute concentrations, or justify why the super-threshold points do not bias the fit, and they need replicates. I would ask for major revision, not desk reject.","headline":"A clever osmotic method and a plausible scaling idea, but most K0 values are measured above the paper's own dilute limit, so the universal law is not yet demonstrated.","tokens_in":15043,"tokens_out":2526,"would_cite":false,"duration_ms":25259,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["82.70.Gg","82.35.Lr"],"model":"deepseek-v4-flash","headline":"The equilibrium partition constant for linear polymer chains entering a flexible polymer network is governed by one dimensionless ratio, $K_0 = \\exp[-4(R_g/l_{\\mathrm{cycle}})^2]$.","keywords":["polymer partitioning","hydrogels","partition constant","gyration radius","cycle length","osmotic pressure","semidilute scaling","mesh-size probe"],"falsifier":"Direct dry-weight measurements, already performed for a subset of samples, should be extended to every network architecture used in the universal-law plot; if for any architecture the internal concentration inferred from swelling disagrees with direct drying under Eq. (2), the extracted $K_0$ shifts and the collapse in Fig. 5 would be violated.","tokens_in":13998,"feed_emoji":"💧","tokens_out":7084,"duration_ms":67258,"temperature":0.7,"pith_summary":"The paper claims that the equilibrium partition constant $K_0$ for linear polymer chains entering a flexible polymer network is set entirely by the network's elastically effective cycle length $l_{\\mathrm{cycle}}$ and the chain's gyration radius $R_g$, through $K_0 = \\exp[-4(R_g/l_{\\mathrm{cycle}})^2]$. To test this, the authors develop a label-free, contactless osmotic method: the extra osmotic pressure from partitioned chains makes the gel swell measurably, and from that swelling they read the internal concentration. They validate the underlying osmotic description by direct drying on a subset of gels, then show that all measured $K_0$ values, spanning three chain sizes and many network densities and topologies, collapse onto the single exponential law. If right, the law replaces the old rigid-rod picture of gel sieving and gives a simple one-length-scale description of partitioning that could be used to estimate mesh sizes from swelling alone.","feed_headline":"One ratio predicts how polymers enter a gel","feed_subtitle":"The partition constant collapses to a single exponential curve in chain size over network cycle length.","key_machinery":"The central object is $l_{\\mathrm{cycle}}$, the mean distance between the centroids of elastically effective cycles in the network, computed from a recursive branching calculation on the star-polymer precursors. It carries the argument because plotting $-\\ln K_0$ against $(R_g/l_{\\mathrm{cycle}})^2$ collapses all data onto a line of slope 4, while plotting against $R_g/l_{\\mathrm{branch}}$ or $R_g/l_{\\mathrm{cross}}$ does not. The measurement machinery is osmotic: the gel's swelling equilibrium, plus the semidilute scaling of mixing pressure with total polymer concentration, converts weight changes into internal concentrations.","core_discovery":"The central discovery is Eq. (5): $K_0 = \\exp[-4(R_g/l_{\\mathrm{cycle}})^2]$, where $R_g$ is the solute gyration radius and $l_{\\mathrm{cycle}}$ is the mean distance between centroids of elastically effective cycles in the network. The paper shows experimentally, on end-linked star-polymer PEG hydrogels, that this single ratio organizes every measured partition constant, for three solute molar masses and networks differing in precursor mass, polymer concentration, branching number $f$, and connectivity $p$. It also shows that the dilute-regime law holds ($K$ is constant at low $c_{\\mathrm{ext}}$) and that the conventional rigid-rod excluded-volume model badly underpredicts $K_0$ and misses its topology dependence. The implied physical picture is that flexible networks restrict chain entry through the conformational cost of threading elastically effective cycles, not through static excluded volume.","pith_inferences":["A natural next test is whether the numerical prefactor 4 is universal across polymer chemistries and solvents, or specific to PEG-water; measuring $K_0$ in another good-solvent flexible network would settle it.","Equation (5) suggests a reciprocal molecular-ruler application: with a calibrated set of chain sizes, a single swelling experiment maps $l_{\\mathrm{cycle}}$, so partitioning data could serve as a non-invasive structural assay for gels whose topology is unknown.","If the law extends to chemically distinct solutes, the prefactor would encode the solute-network interaction alone, providing a clean separation of size exclusion from interaction effects in membrane and chromatography design.","Because the collapse uses one universal curve, systematic deviations in a new dataset would flag either a breakdown of the osmotic scaling relation or a network defect population that the recursive branching cycle count misses."],"forward_implications":["Given $R_g$ of a solute and $l_{\\mathrm{cycle}}$ of a network, $K_0$ can be predicted with no adjustable parameters.","Swelling a gel in a dilute solution of chains with known $R_g$ becomes a contactless probe of $l_{\\mathrm{cycle}}$, that is, of the elastically effective mesh size.","Dilute-regime partition measurements can be summarized by a single constant $K_0$ rather than a concentration-dependent ratio.","The rod-obstacle picture, which predicts $K_0$ from excluded volume alone, is inadequate for flexible networks because it ignores the topology-dependent cycle structure.","Because $l_{\\mathrm{cycle}}$ depends on connectivity $p$ and functionality $f$, the law predicts that defects that break elastically effective cycles will raise $K_0$ (less exclusion) at fixed polymer concentration."],"supporting_citations":[{"why":"The rigid-rod exclusion model whose predictions the new law is tested against.","marker":"[1]"},{"why":"The gel-filtration extension of the rod model, another baseline for $K_0$.","marker":"[2]"},{"why":"Establishes the universal semidilute osmotic equation of state for gels used to infer $c_{\\mathrm{int}}$.","marker":"[26]"},{"why":"Extends the semidilute scaling principle to gels, supporting the total-concentration mixing pressure.","marker":"[27]"},{"why":"Provides the osmotic equation of state for external linear PEG solutions used to set $\\Pi_{\\mathrm{ext}}$.","marker":"[41]"},{"why":"Supplies shear modulus data used to evaluate the elastic contribution $\\Pi_{\\mathrm{el}}$.","marker":"[44]"},{"why":"Fixes the $R_g$ versus molar mass scaling used to compute solute gyration radii.","marker":"[45]"},{"why":"Branching-statistics recursion used to count elastically effective chains and cycles.","marker":"[46]"},{"why":"Post-gel properties calculation that yields the cycle count entering $l_{\\mathrm{cycle}}$.","marker":"[47]"}],"fun_headline_variants":["One ratio predicts polymer partitioning into gels","Universal law: chain-to-mesh ratio sets partition constant","Flexible network uptake: a single exponential law","Chain size over mesh size decides if a polymer enters a gel","New exponential law for polymer entry into flexible networks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the mixing osmotic pressure inside a gel that has absorbed partitioned chains is the same function of the total polymer concentration as in a chain-free solution, so that the internal concentration can be read off from swelling alone.","fun_headline_variants_meta":{"raw":{"variants":["One ratio predicts polymer partitioning into gels","Universal law: chain-to-mesh ratio sets partition constant","Flexible network uptake: a single exponential law","Chain size over mesh size decides if a polymer enters a gel","New exponential law for polymer entry into flexible networks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001036,"raw_usage":{"total_tokens":4318,"prompt_tokens":859,"completion_tokens":3459,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":475,"completion_tokens_details":{"reasoning_tokens":3385}},"tokens_in":475,"tokens_out":3459,"duration_ms":25867,"temperature":1.0,"reasoning_tokens":3385,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T23:09:03.827841+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Direct dry-weight measurements, already performed for a subset of samples, should be extended to every network architecture used in the universal-law plot; if for any architecture the internal concentration inferred from swelling disagrees with direct drying under Eq. (2), the extracted $K_0$ shifts and the collapse in Fig. 5 would be violated.","supporting_citations":[{"cited_title":"Universality of Osmotic Equation of State in Star Polymer Solutions","cited_arxiv_id":"2302.13669","evidence_quote":"Supplies shear modulus data used to evaluate the elastic contribution $\\Pi_{\\mathrm{el}}$."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The rigid-rod exclusion model whose predictions the new law is tested against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The gel-filtration extension of the rod model, another baseline for $K_0$."},{"cited_title":"Shibayama, Spatial inhomogeneity and dynamic fluc- tuations of polymer gels, Macromol","cited_arxiv_id":null,"evidence_quote":"Establishes the universal semidilute osmotic equation of state for gels used to infer $c_{\\mathrm{int}}$."},{"cited_title":"Yasuda, N","cited_arxiv_id":null,"evidence_quote":"Extends the semidilute scaling principle to gels, supporting the total-concentration mixing pressure."},{"cited_title":"Yoshikawa, N","cited_arxiv_id":null,"evidence_quote":"Provides the osmotic equation of state for external linear PEG solutions used to set $\\Pi_{\\mathrm{ext}}$."},{"cited_title":"Yoshikawa, N","cited_arxiv_id":null,"evidence_quote":"Fixes the $R_g$ versus molar mass scaling used to compute solute gyration radii."},{"cited_title":"Kawaguchi, G","cited_arxiv_id":null,"evidence_quote":"Branching-statistics recursion used to count elastically effective chains and cycles."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Post-gel properties calculation that yields the cycle count entering $l_{\\mathrm{cycle}}$."}],"review_version":1}