{"id":"132b0b49-0a98-4335-a67e-beccd984bf5b","arxiv_id":"1908.07768","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Fluid-immersed disordered spring networks exhibit a new intermediate viscoelastic plateau, absent in vacuum, in which both network and fluid respond non-affinely.","lead":"This paper simulates disordered spring networks immersed in a viscous fluid and finds a previously unseen response regime where both the network and the fluid deform non-affinely, producing an extra plateau in the elastic modulus. The result suggests immersed biopolymer networks and hydrogels have a richer frequency-dependent mechanical response than vacuum models predict.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (4) defines the stress with f_αβ = f_β − f_α (difference of total node forces), not the bond force; this is not a valid virial expression and the reported G*(ω) spectra are suspect.","rationale":"The reader identified the node-only point-force coupling as the weakest assumption, and that is indeed a legitimate limitation the authors themselves flag. However, the more pressing issue is the internal inconsistency in Eq. (4), which defines the very observable from which the central plateau is extracted. The expression as written uses the difference of total node forces rather than the bond force; this is not the standard virial stress and can fail qualitatively even for a uniformly strained linear chain. If the code follows Eq. (4), all quantitative results are invalid. If it is merely a typographical error, the paper is still unreproducible because the correct formula is not supplied. This concern is more direct than the point-coupling assumption because it questions the measured G*(ω) itself, not the physical range of applicability. The concrete test of comparing the two stress expressions will settle whether the manuscript is fundamentally flawed or just miswritten. Until that is resolved, the paper should not be accepted in any form, not even conditionally, since the central data cannot be trusted as reported. Hence the verdict should move from CONDITIONAL to UNVERDICTED.","tokens_in":7750,"tokens_out":17043,"duration_ms":180884,"concrete_test":"Implement both stress expressions for the same p = 0.8 network at ζ/η = 10^5 and a single frequency (e.g., ωη/k = 10^−3): (a) the formula in Eq. (4) exactly as written, and (b) the standard virial Σ_bonds F_αβ^x r_αβ^y with F_αβ = k[(u_β − u_α)·t̂]t̂. Also test a three-node chain under uniform strain. If the two results differ (they must in general), Eq. (4) is incorrect and the manuscript requires major revision; if they coincide, the notation in Eq. (4) has a different meaning that must be clarified.","verdict_should_be":"UNVERDICTED","load_bearing_attack":"Equation (4) is the central observable: every plateau and peak is read from G*(ω). As written, f_αβ is defined as the difference of total node forces f_β − f_α, where f_α is the sum of all spring forces at node α (Eq. 2). The virial shear stress of a central-force network requires the bond force F_αβ = k[(u_β − u_α)·t̂]t̂ contracted with the bond vector r_β − r_α. The difference f_β − f_α includes forces from all other bonds incident on α and β; it is not the force in bond αβ. Counterexample: a uniformly strained one-dimensional chain has f_α = 0 at every interior node, so Eq. (4) returns zero, while the true stress is nonzero. Unless the numerical code actually used the correct bond-force virial and Eq. (4) is a typo, the intermediate plateau and both loss peaks could be artifacts of an incorrect stress evaluation. This is an internal inconsistency in the paper's central calculation, not a boundary-condition or model-choice issue.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript describes a matrix-based numerical solver for the linear oscillatory rheology of two-dimensional disordered spring networks immersed in a Stokes fluid, with network–fluid coupling through point drag at the nodes. The authors scan the drag coefficient ζ, fluid viscosity η, and frequency ω and identify four response regimes from affinity and coupling metrics. Their main claim is a previously unreported intermediate regime at ζ≫η in which both the fluid and the network respond non-affinely and remain coupled, producing an additional storage-modulus plateau G_NA between the low-frequency G_0 and the high-frequency affine G_aff, with two peaks in the network contribution to the loss modulus at the plateau boundaries. They also report linear dependence of the plateau moduli on bond density p and note that finite-frequency driving rigidifies sub-isostatic networks. The central results are presented as averaged spectra over ten 100×100 networks, supplemented by checks of mesh insensitivity and incompressibility.","tokens_in":7959,"tokens_out":12875,"duration_ms":198821,"significance":"If correct, the claimed fluid-coupled non-affine plateau would be a genuinely new element in the rheology of immersed fiber networks, distinct from the vacuum response and from one-way coupled models; it is also potentially relevant to biopolymer networks and hydrogels. The systematic parameter scan and the matrix-based solver are useful technical contributions, and the authors are careful to check mesh insensitivity and incompressibility and to average over multiple realizations. The significance, however, cannot be fully assessed from the manuscript as written: the stress formula in Eq. (4) is questionable, and the plateau characterization relies on a Maxwell-mode fit without error bars or residuals. These issues must be resolved before the new regime can be considered established.","major_comments":[{"comment":"The stress formula in Eq. (4) uses f_αβ ≡ f_β − f_α, the difference of the total nodal forces, rather than the force carried by the bond αβ. The standard pair virial for a central-force network requires F_αβ = k[(u_β − u_α)·t̂]t̂ contracted with (x_β − x_α). As written, Eq. (4) fails even in a uniformly strained one-dimensional chain, where interior nodes have f_α = 0 and the expression returns zero while every bond carries a nonzero force. Since all plateau and loss-peak claims are read from G*(ω) computed by this formula, this is a load-bearing issue. If the code uses the correct bond-force virial, Eq. (4) and the definition of f_αβ must be corrected; if the code implements Eq. (4) as written, the reported spectra do not measure the shear modulus.","section":"Viscoelastic spectra, Eq. (4)"},{"comment":"The identification of the intermediate plateau and the two loss peaks is obtained by fitting G′ and G″ to a two-Maxwell-unit ansatz whose two relaxation rates are by construction the fitted parameters. No error bars are shown for the modulus curves and no fit residuals are reported, and the finding ζ_eff≈ζ and η_eff≈η is therefore a consistency check built into the fit form rather than independent evidence. Please report per-frequency standard errors (the text says each curve is an average over 10 networks), the fit residuals, and a model-independent plateau criterion, such as the frequency range over which d log G′/d log ω stays within a small tolerance of zero and over which two resolved maxima in G″−ωη are visible.","section":"Viscoelastic spectra, Eqs. (5)–(6) and Fig. 4"},{"comment":"The new regime is identified from smooth crossovers in the metrics, and the scaling argument in the text explicitly stops short of predicting it: the text states that 'no terms in (3) dominate, so a limiting solution for v cannot be inferred' and that 'Fig. 3 suggests a smooth crossover'. To support the claim that this is a distinct regime rather than a crossover, please provide an independent diagnostic of the regime boundaries, for example a collapse of Γ_net, Γ_fl, and Γ_dc on ωζ/k and ωη/k, or a direct demonstration that the plateau width scales as ζ/η. Such a test would tie the plateau in G*(ω) to the regime identified from the metrics without relying on the fitted Maxwell modes.","section":"Response regimes and Fig. 3"}],"minor_comments":[{"comment":"The axis labels in Fig. 4 are garbled: '10 3' should be '10^{-3}', and '(G′′ i ) / k' should be '(G′′ − ωη)/k' or similar.","section":"Fig. 4"},{"comment":"The name 'Tildesley' is misspelled as 'Tildedsley'.","section":"Reference [28]"},{"comment":"Fig. 5 shows no error bars for G_NA, G_0, or Γ_net even though each data point is an average over 10 networks; the uncertainties on the fitted p* values are given, but the underlying scatter should be visible.","section":"Fig. 5"},{"comment":"The Supplemental Material (Figs. S1–S4) is cited but was not included with the manuscript; it should be submitted for review because the mesh-insensitivity and incompressibility checks and the p<p_c spectra are used to support claims in the main text.","section":"Supplemental Material"},{"comment":"In the Methodology paragraph, 'The fluid velocity vα(t) at mesh nodes' should read 'at network nodes', since the staggered mesh refers to the fluid grid points.","section":"Methodology"}],"recommendation":"major_revision","confidential_remarks":"The decisive question is whether Eq. (4) is a typographical error or a faithful description of the computation. If it is the latter, the central claim of the paper is unsupported and rejection would be appropriate; if it is a typo, the paper may be salvageable with a careful revision. I recommend asking the authors to supply a verification test for the stress calculation (e.g., an affine deformation with known modulus) as part of the revision. Also, the Supplemental Material should be required for review, since the main text relies on Figs. S1–S4."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThe thing to know: this paper reports a genuinely new intermediate plateau in the storage modulus of fluid-immersed spring networks, in a regime where both network and fluid deform non-affinely but stay coupled. The regime sits between the usual low-frequency non-affine plateau and the high-frequency affine plateau, so it expands the known rheology of immersed fiber networks. The claim looks plausible: the numerics are standard two-fluid Stokes plus spring network, the regime diagram in Fig. 1 is well thought out, and the comparison with the no-hydrodynamic-interactions case (insets in Fig. 4) shows the plateau is tied to full coupling. The p-dependence of the intermediate modulus G_NA differs from G_0, which supports the claim that this is a distinct elastic state, not just a rescaling of the vacuum response.\n\nWhat the paper does well: clean numerical study, honest about the point-force drag assumption (flagged in the Discussion), and the scaling arguments for the regimes are reasonable. The fits give ζ_eff ≈ ζ and η_eff ≈ η, which is reassuring.\n\nNow the soft spots, in descending severity.\n\nFirst, Eq. (4) is not a valid virial stress as written. The text says f_αβ ≡ f_β − f_α, where f_α is the total spring force at node α (Eq. 2). The true bond stress requires the force in bond αβ, i.e., k[(u_β−u_α)·t̂]t̂, contracted with r_β−r_α. The difference of total node forces includes contributions from all other bonds and is not the bond force. A uniform strain on a chain gives zero net forces at interior nodes, so Eq. (4) would return zero stress while the true stress is nonzero. Unless the code used the correct bond-force virial and Eq. (4) is just a typo, the plateau and the loss peaks could be artifacts of a wrong stress evaluation. The paper includes no code or data to check. This is the biggest issue and needs to be fixed before I'd trust the spectra.\n\nSecond, there are no error bars on the modulus curves. They average over 10 networks, but the spread is not shown. The plateau is identified by fitting G′ and G″ to a two-Maxwell-unit model; without error bars and without independent verification, the plateau location and boundaries are less solid than the plots suggest.\n\nThird, the node-only coupling approximation is acknowledged but not tested against distributed or rod-based coupling. The authors argue it should not change the picture, but that is an assertion, not a demonstration.\n\nThe physics is interesting and the authors are clearly competent. If Eq. (4) is corrected and the spectrum data are made available, this would be a solid paper. As it stands, it deserves a careful referee, but the referee should ask for the corrected stress formula and the code. I would not cite it in its current form.\n\nRecommended for peer review: yes, with the expectation of revision.","headline":"A plausible new intermediate plateau in immersed-network rheology, but Eq. (4) as printed is not a valid stress expression and needs fixing before the spectra can be trusted.","tokens_in":8513,"tokens_out":6863,"would_cite":false,"duration_ms":66832,"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":"A fluid-coupled non-affine plateau appears in immersed fiber network rheology.","keywords":["immersed fiber networks","non-affine deformation","viscoelastic plateau","hydrodynamic interactions","Stokes flow","spring networks","two-fluid model","complex shear modulus"],"falsifier":"Replace the nodal point drag in the same coupled equations with a drag distributed continuously along each spring segment (for instance by discretizing each fiber into beads that all experience solvent drag), and compute the complex modulus across ωη/k ≪ 1 ≪ ωζ/k; if the intermediate plateau and its twin loss peaks vanish, the central claim is false. A bulk rheology experiment on a well-characterized immersed fiber network with independently known ζ/η would also test it: the predicted plateau and two loss peaks must appear between the low- and high-frequency limits.","tokens_in":7504,"feed_emoji":"🕸️","tokens_out":6069,"duration_ms":225472,"temperature":0.7,"pith_summary":"The paper asks how a disordered network of Hookean springs responds to oscillatory shear when it sits in a viscous fluid, with drag forces coupling the two at every network node. It claims that sweeping the drag coefficient and frequency reveals four response regimes, one of which has no counterpart in a dry network: at high drag and intermediate frequencies both the network displacements and the fluid flow are non-affine, yet strongly coupled. In that regime the storage modulus develops an extra plateau at a value between the low-frequency modulus and the high-frequency affine modulus, with two peaks in the loss modulus at the plateau's boundaries. If this is right, immersed biopolymer networks and hydrogels have a richer linear viscoelastic response than vacuum-based models suggest, and fluid coupling becomes a design handle.","feed_headline":"Immersed fiber networks gain a fluid-driven stiffness plateau","feed_subtitle":"When drag dominates viscosity, both network and fluid deform non-affinely, adding a storage plateau between two known limits.","key_machinery":"The machinery is a matrix-based solver for the linear oscillatory response of a bond-diluted triangular spring network coupled to a viscous solvent through point drag forces at network nodes: each node obeys f_α = ζ(∂_t u_α − v_α), and the fluid satisfies Stokes' equations with the nodal forces as delta-function sources. The coupled equations are discretized on staggered meshes and solved as one sparse matrix equation for the complex amplitudes of displacements, fluid velocities, and pressures under oscillatory shear. Three dimensionless metrics quantify non-affinity of the network, non-affinity of the fluid, and decoupling between the two, and the complex modulus is fit to a parallel spring plus two Maxwell units to extract the plateau's characteristic rates.","core_discovery":"The central discovery is that a fluid-immersed, disordered spring network can store elastic energy in an intermediate-frequency regime in which neither the network displacement field nor the fluid velocity field is affine, while the two fields remain strongly coupled. For drag coefficient ζ much larger than the solvent viscosity η and for frequencies satisfying ωη/k ≪ 1 ≪ ωζ/k, the storage modulus G′(ω) flattens into a plateau G_NA lying between the low-frequency modulus G0 and the high-frequency affine modulus G_aff = pk√3/4, and the network contribution to G″(ω) shows two peaks at the plateau edges. The plateau is absent when hydrodynamic interactions are removed by replacing the true fluid velocity with its affine prediction. Fitting the spectra to a spring plus two viscoelastic Maxwell elements gives effective rates ζeff ≈ ζ and ηeff ≈ η, and the plateau value scales linearly with p − p∗ in a way that cannot be mapped onto G0 by changing network density.","pith_inferences":["If node-only drag is replaced by distributed drag along fiber segments, the plateau may persist with only rescaled relaxation rates; a clean way to check is a bead-spring or immersed-boundary simulation where the fluid force is applied along the whole spring, not just at endpoints.","The plateau's high and low edges are set by η and ζ separately, so varying solvent viscosity at fixed network density should shift only one edge; measuring that would show whether the two-Maxwell description is physically meaningful.","The paper's finite-frequency rigidification of sub-isostatic networks suggests that many 'soft' fiber materials may appear rigid in practice, with the fluid-controlled plateau providing a second stiffness scale that cells or other sensors could in principle detect."],"forward_implications":["For p > p_c, strong drag produces an intermediate plateau G_NA in the storage modulus, bracketed by two peaks in the network contribution to the loss modulus.","The plateau disappears if hydrodynamic interactions are switched off, so any experiment that sees it is direct evidence of full fluid–network coupling.","The fitted plateau frequencies scale with the bare parameters, ζeff ≈ ζ and ηeff ≈ η, so the plateau location is set by the same drag and viscosity that define the regimes.","Below the rigidity threshold, any finite frequency gives terminal scaling G′ ∼ ω² and G″ ∼ ω, meaning the solvent rigidifies floppy networks at all accessible frequencies.","Network configurations in the intermediate plateau deform differently from vacuum non-affinity at the same modulus, so tuning p cannot reproduce the plateau by rescaling a vacuum network."],"supporting_citations":[{"why":"Supplies the two-fluid model of network and solvent that the paper discretizes.","marker":"[25]"},{"why":"Provides the bond-diluted triangular spring-network construction and previous spring-network rheology.","marker":"[16]"},{"why":"Previous treatment of hydrodynamic interactions in such networks, extended here to strong drag.","marker":"[18]"},{"why":"Rod-based network-fluid coupling with the same high-frequency affine limit, used to argue node-only coupling suffices.","marker":"[15]"},{"why":"Origin of the non-affinity metric generalized to complex displacement fields.","marker":"[10]"},{"why":"Supplemental checks of mesh insensitivity and the fitted plateau rates.","marker":"[27]"},{"why":"Sparse direct solver used to invert the coupled system matrix.","marker":"[29]"},{"why":"Bead-spring polymer dynamics formalism invoked to extrapolate the node-coupling results to other coupling schemes.","marker":"[4]"}],"fun_headline_variants":["Fluid drag adds a non-affine stiffness plateau to fiber networks","Immersed networks gain a fluid-coupled elastic plateau","Hydrodynamics yields an extra plateau in immersed fiber networks","Non-affine fluid coupling creates a storage plateau in networks","Fluid-coupled non-affine plateau in immersed fiber networks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the fluid–network interaction can be represented entirely as point drag forces acting at the discrete network nodes; if realistic coupling along the full length of each fiber changes the physics, the intermediate plateau may disappear.","fun_headline_variants_meta":{"raw":{"variants":["Fluid drag adds a non-affine stiffness plateau to fiber networks","Immersed networks gain a fluid-coupled elastic plateau","Hydrodynamics yields an extra plateau in immersed fiber networks","Non-affine fluid coupling creates a storage plateau in networks","Fluid-coupled non-affine plateau in immersed fiber networks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000242,"raw_usage":{"total_tokens":1467,"prompt_tokens":831,"completion_tokens":636,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":447,"completion_tokens_details":{"reasoning_tokens":555}},"tokens_in":447,"tokens_out":636,"duration_ms":6113,"temperature":1.0,"reasoning_tokens":555,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:56:38.453011+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Replace the nodal point drag in the same coupled equations with a drag distributed continuously along each spring segment (for instance by discretizing each fiber into beads that all experience solvent drag), and compute the complex modulus across ωη/k ≪ 1 ≪ ωζ/k; if the intermediate plateau and its twin loss peaks vanish, the central claim is false. A bulk rheology experiment on a well-characterized immersed fiber network with independently known ζ/η would also test it: the predicted plateau and two loss peaks must appear between the low- and high-frequency limits.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the two-fluid model of network and solvent that the paper discretizes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the bond-diluted triangular spring-network construction and previous spring-network rheology."},{"cited_title":"Dennison and H","cited_arxiv_id":null,"evidence_quote":"Previous treatment of hydrodynamic interactions in such networks, extended here to strong drag."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Rod-based network-fluid coupling with the same high-frequency affine limit, used to argue node-only coupling suffices."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Origin of the non-affinity metric generalized to complex displacement fields."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplemental checks of mesh insensitivity and the fitted plateau rates."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Sparse direct solver used to invert the coupled system matrix."}],"review_version":1}