Recognition: unknown
Linear and Non-Linear Rheology of Single and Double Cross-Linked Biopolymer Networks under Viscous Shear Flow
Pith reviewed 2026-05-07 13:12 UTC · model grok-4.3
The pith
Two-step yielding occurs in single biopolymer networks, and double-network stress curves match summed single networks only under large strains.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The two-step yielding (double peaks) reported for double-crosslinked networks at high strain amplitude is not caused by fiber alignment changes or seed numbers despite the presence of crosslinkers linking the subnetworks; the same peaks appear in single networks. In addition, the stress-strain response of a double network equals the sum of the two single-network responses only in the nonlinear regime, while the same superposition fails in the linear regime. Higher-order Fourier modes up to the fifth are prominent in double networks under nonlinear conditions, and nonlinearity strength in single networks depends on the initial structural seed number rather than flow conditions.
What carries the argument
Modified slender body theory that incorporates hydrodynamic interactions among fibers within single and intertwined double networks to compute stress and Fourier coefficients under oscillatory shear.
If this is right
- Two-step yielding is an intrinsic feature of these fiber networks and does not require inter-subnetwork crosslinkers.
- The nonlinear mechanical response of double networks can be approximated by adding the responses of the constituent single networks under large-amplitude shear.
- Controlling seed number during network formation offers a route to tune nonlinearity strength independently of applied flow conditions.
- Morphology-based design of initial structure can be used to engineer desired nonlinear rheology in fibrous biomaterials.
Where Pith is reading between the lines
- Single-network simulations may be sufficient to estimate double-network behavior at high strains, lowering computational cost for design studies.
- Experimental checks of fiber alignment under shear would clarify whether other mechanisms besides morphology produce the observed double peaks.
- The morphology dependence suggests that fabrication protocols that fix seed-like parameters could control nonlinear properties in tissue-engineering scaffolds.
Load-bearing premise
The chosen network generation rules and modified slender body theory correctly reproduce the hydrodynamic interactions and yielding processes that occur in real biopolymer networks.
What would settle it
An experiment on real double-crosslinked hydrogels at 100 percent strain amplitude whose measured stress-strain curve deviates from the sum of the two single-network curves, or whose single networks lack double peaks when fiber alignments are matched to the double-network case.
read the original abstract
In this research study, a numerical tool, which is based on a version of Slender Body theory, has been used and also modified to simulate the mechanical behaviour of single- and double-cross-linked biopolymer networks (hydrogel) under oscillatory shear flow. The hydrodynamic interactions among fibres of intertwined networks were considered. Then, the stress and Fourier coefficients (i.e. shear moduli) were evaluated for both linear and nonlinear regimes. It was found that the double peaks (two-step yielding) of two double network at 100% maximum strain amplitude (nonlinear regime) cannot happen due to changes in fibre alignments and seed numbers, although the crosslinkers between two subnetworks present, which was previously reported in the literature. In fact, we also observed two peaks for single network in nonlinear regime. Furthermore, it was shown that the stress-strain curve of double network is not predicted by just superimposing the results from the corresponding single networks at 5% maximum strain amplitude (linear regime), but this prediction can be provided at 100% maximum strain amplitude (nonlinear regime). The Fourier coefficients and corresponding amplitude (an indication of nonlinearity effects) for double network were quite considerable from zero to fifth modes in nonlinear regime, despite enough zero and first modes in linear regime. It was also shown that the nonlinearity effects can be related to the morphology of the initial structure, i.e. the seed number rather than the flow condition for the single network. These results can help scientists to better design enhance fibrous materials used in wound healing or tissue engineering.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents numerical simulations of single- and double-crosslinked biopolymer networks under oscillatory shear using a modified Slender Body theory that incorporates hydrodynamic interactions. It claims that two-step yielding (double peaks) does not appear in double networks at 100% strain amplitude due to fiber alignment changes and seed-number effects despite inter-subnetwork crosslinkers, that single networks also exhibit two peaks in the nonlinear regime, that stress-strain superposition of single-network results predicts double-network behavior only in the nonlinear regime (not at 5% linear strain), and that nonlinearity (via Fourier coefficients up to fifth mode) correlates more with initial morphology (seed number) than flow conditions.
Significance. If the model were validated, the findings on morphology-dependent yielding and regime-specific superposition would be relevant for understanding composite biopolymer hydrogel mechanics and could inform design of materials for tissue engineering. The inclusion of hydrodynamic interactions is a methodological strength. However, the absence of experimental calibration or benchmarks means the results currently offer limited insight beyond the specific simulation setup.
major comments (3)
- [Abstract and nonlinear regime results] Abstract and results on double networks: The central claim that double peaks 'cannot happen' at 100% strain due to fiber alignments and seed numbers (despite crosslinkers) directly contradicts prior literature on two-step yielding, yet the manuscript provides no quantitative alignment metrics (e.g., nematic order parameter), no systematic variation over multiple seeds, and no error bars or statistical robustness checks. This makes the explanation for the discrepancy unconvincing and load-bearing for the novelty assertion.
- [Methods] Methods and model description: All reported behaviors depend on the fidelity of the modified Slender Body theory and chosen network generation parameters (including seed numbers). No validation against experimental shear moduli, yielding strains, or known hydrodynamic effects in real biopolymer gels (e.g., actin or fibrin networks) is presented, nor is sensitivity to the theory modifications or parameter choices demonstrated. This is load-bearing for the reliability of both the no-two-step-yielding and superposition claims.
- [Results on superposition] Superposition analysis in linear vs. nonlinear regimes: The assertion that double-network stress-strain curves are not predicted by superposition at 5% strain but are at 100% strain lacks explicit description of the superposition procedure, quantitative error metrics (e.g., deviation between superimposed and simulated curves), and controls for how single-network results were scaled or combined. Without these, the regime-dependent prediction claim cannot be evaluated.
minor comments (3)
- [Abstract] Abstract contains grammatical issues ('two double network', 'enhance fibrous materials') and vague phrasing ('quite considerable from zero to fifth modes') that should be clarified for precision.
- [Introduction or Discussion] The manuscript should cite the specific prior literature on two-step yielding in double networks that is being contradicted, to allow readers to assess the discrepancy.
- [Results] No details are given on how fiber alignments were measured or visualized, nor on the number of independent realizations per seed number.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed comments on our manuscript. We address each major point below, clarifying our findings and outlining revisions where appropriate to improve clarity and robustness.
read point-by-point responses
-
Referee: [Abstract and nonlinear regime results] Abstract and results on double networks: The central claim that double peaks 'cannot happen' at 100% strain due to fiber alignments and seed numbers (despite crosslinkers) directly contradicts prior literature on two-step yielding, yet the manuscript provides no quantitative alignment metrics (e.g., nematic order parameter), no systematic variation over multiple seeds, and no error bars or statistical robustness checks. This makes the explanation for the discrepancy unconvincing and load-bearing for the novelty assertion.
Authors: Our simulations demonstrate that at 100% strain amplitude, fiber alignment under viscous shear suppresses the second yielding peak in double networks, even with inter-subnetwork crosslinkers, and this behavior is also seen in single networks. We tested multiple seed numbers for network generation and found consistent results, but we agree that quantitative support is needed. In the revision, we will compute and report the nematic order parameter to measure alignment changes, include error bars from multiple independent realizations, and add a discussion reconciling our hydrodynamic-inclusive results with prior literature on two-step yielding (which often lacks such interactions). This will strengthen the explanation without altering the core claim. revision: partial
-
Referee: [Methods] Methods and model description: All reported behaviors depend on the fidelity of the modified Slender Body theory and chosen network generation parameters (including seed numbers). No validation against experimental shear moduli, yielding strains, or known hydrodynamic effects in real biopolymer gels (e.g., actin or fibrin networks) is presented, nor is sensitivity to the theory modifications or parameter choices demonstrated. This is load-bearing for the reliability of both the no-two-step-yielding and superposition claims.
Authors: The modified Slender Body theory with hydrodynamic interactions builds on established approaches for semiflexible fiber networks. We will expand the methods section to include a sensitivity analysis varying seed numbers, crosslink densities, and hydrodynamic strength, showing that the reported behaviors (including absence of two-step yielding and regime-specific superposition) are robust. We will also compare linear-regime shear moduli to values reported in the literature for similar biopolymer systems. As this is a numerical study, full experimental calibration of the specific parameters is not feasible here but will be noted as a direction for future work. revision: partial
-
Referee: [Results on superposition] Superposition analysis in linear vs. nonlinear regimes: The assertion that double-network stress-strain curves are not predicted by superposition at 5% strain but are at 100% strain lacks explicit description of the superposition procedure, quantitative error metrics (e.g., deviation between superimposed and simulated curves), and controls for how single-network results were scaled or combined. Without these, the regime-dependent prediction claim cannot be evaluated.
Authors: We will add an explicit description of the superposition procedure, detailing the scaling (based on relative crosslink contributions and volume fractions) and combination of single-network stress-strain data. Quantitative metrics, such as root-mean-square deviation between the superimposed prediction and the directly simulated double-network response, will be included for both 5% and 100% strain amplitudes, along with controls for different scaling assumptions. This will allow readers to evaluate the regime dependence directly. revision: yes
- Direct experimental validation or calibration against specific biopolymer gel experiments (e.g., measured yielding strains in actin or fibrin networks), which would require new laboratory work beyond the scope of this simulation-based study.
Circularity Check
No circularity: results are direct simulation outputs with no self-referential definitions or load-bearing self-citations
full rationale
The paper's claims rest on numerical simulations of single- and double-crosslinked networks under oscillatory shear using a modified Slender Body theory that incorporates hydrodynamic interactions. The reported absence of two-step yielding in double networks at 100% strain (attributed to fiber alignment and seed-number effects) and the regime-dependent superposition of stress-strain curves are direct outputs of running the model on generated networks with varying seeds. No equations define a target quantity in terms of itself, no parameters are fitted to the reported peaks or superposition behavior before being called predictions, and the abstract and summary contain no self-citations that justify core premises. The chain is therefore self-contained: network generation plus hydrodynamic evolution produces the observed stress and Fourier coefficients without reducing the findings to tautologies.
Axiom & Free-Parameter Ledger
free parameters (1)
- seed number
axioms (1)
- domain assumption Slender Body theory provides a valid approximation for the hydrodynamics of slender fibers in viscous flow when modified for network interactions.
Reference graph
Works this paper leans on
-
[1]
J. S. Palmer, M. C. Boyce, Constitutive modeling of the stress –strain behavior of F -actin filament networks. Acta biomaterialia 4, 597-612 (2008)
2008
-
[2]
Ban et al., Strong triaxial coupling and anomalous Poisson effect in collagen networks
E. Ban et al., Strong triaxial coupling and anomalous Poisson effect in collagen networks. Proceedings of the National Academy of Sciences 116, 6790-6799 (2019)
2019
-
[3]
J. M. Jimenez et al. , Multiscale mechanical characterization and comput ational modeling of fibrin gels. Acta biomaterialia 162, 292-303 (2023)
2023
-
[4]
Kopnar, A
V. Kopnar, A. O’Connell, N. Shirshova, A. Aufderhorst -Roberts, Mechanistic Origins of Yielding in Hybrid Double-Network Hydrogels. Macromolecules 58, 8610-8621 (2025)
2025
-
[5]
M. L. Mug nai, R. Tchuenkam Batoum, E. Del Gado, Interspecies interactions in dual, fibrous gels enable control of gel structure and rheology. Proceedings of the National Academy of Sciences 122, e2423293122 (2025)
2025
-
[6]
Hajaliakbari, D
N. Hajaliakbari, D. Head, O. Harlen, The Effect s of Temperature on Hydrodynamic Interaction of Sedimenting Semi-Flexible Brownian Fibres
-
[7]
Hajaliakbari, D
N. Hajaliakbari, D. Head, O. G. Harlen, Hydrodynamic interactions between pairs of sedimenting semi-flexible Brownian fibers. Physics of Fluids 38 (2026)
2026
-
[8]
Maxian, A
O. Maxian, A. Donev, A simulation platform for slender, semiflexible, and inextensible fibers with Brownian hydrodynamics a nd steric repulsion. Physics of Fluids 36 (2024). 14
2024
-
[9]
Maxian, Hydrodynamics of transiently cross-linked actin networks: theory, numerics, and emergent behaviors (New York University, 2023)
O. Maxian, Hydrodynamics of transiently cross-linked actin networks: theory, numerics, and emergent behaviors (New York University, 2023)
2023
-
[10]
Maxian, R
O. Maxian, R. P. Peláez, A. Mogilner, A. Donev, Simulations of dynamically cross -linked actin networks: morphology, rheology, and hydrodynamic interactions. PLoS computational biology 17, e1009240 (2021)
2021
-
[11]
Maxian, A
O. Maxian, A. Mogilner, A. Donev, Integral -based spectral method for inextensible slender fibers in Stokes flow. Physical Review Fluids 6, 014102 (2021)
2021
-
[12]
J. M. Stockie (1997) Analysis and computation of immersed boundaries, with application to pulp fibres. (University of British Columbia)
1997
-
[13]
Livraghi, K
M. Livraghi, K. H öllring, C. R. Wick, D. M. Smith, A. -S. a. Smith, An exact algorithm to detect the percolation transition in molecular dynamics simulations of cross -linking polymer networks. Journal of chemical theory and computation 17, 6449-6457 (2021)
2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.