Pith. sign in

REVIEW 3 major objections 6 minor 30 references

Rectification of stress by fiber networks: Manifestation of non-linear screening through self-organized buckling

T0 review · 3 major / 6 minor · reviewed 2026-08-03 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The full nonlinear response of a buckling fiber network is captured by a renormalized Poisson ratio, which rises from ~1/3 to ~1 and explains stress rectification.

desk verdict A credible simulation study showing a renormalized Poisson ratio in buckled networks, but the central claim is fitted in-sample under an untested isotropy assumption — promising, not proven. read the letter →

arxiv 2511.20754 v1 pith:7EJ6XQ6Z submitted 2025-11-25 physics.bio-ph cond-mat.soft

Kanaya Malakar , Albert Countryman , Bulbul Chakraborty This is my paper · ORCID
classification physics.bio-phcond-mat.soft
keywords nonlinearelasticitybucklingPoissonratiomechanicalscreeningprestressforcedipolefibernetworkKagomelattice
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper establishes that a mildly buckling fiber network, driven by a microscopic extensile force dipole, can be described at large scales by ordinary isotropic elasticity with a single renormalized parameter: the Poisson ratio. As the dipole strength grows, the fitted Poisson ratio rises from the linear-network value of about 1/3 to essentially 1, the stability limit in two dimensions. That single number reproduces the simulated far-field stress in all three mechanical regimes the authors identify—linear, nonlinear, and rectified—including the striking rectified phase in which a contractile far-field response appears despite an extensile applied force. The mechanism is 'mechanical screening': the force dipole polarizes the network through self-organized buckling patterns (twisted units on the midpoints' Kagome lattice), and the prestressed network responds as an effective medium that can mask or even overshoot the sign of the source. If correct, this gives a single-parameter route to predicting force transmission in biological fiber networks without invoking strain or a reference state, and it identifies the rectified phase as a precursor to shear instability and flow.

What carries the argument

The carrying object is the force-dependent Poisson ratio ν(|F|), obtained by fitting the angular profile of the far-field force-moment tensor D (measured on an annulus at radius 20) to the VCT prediction for an isotropic elastic modulus tensor K. VCT is a continuum theory that recasts prestressed elasticity in the language of electrostatics: the external force dipole acts as a vector 'charge,' the network's buckling deformations generate a polarization field P, and the screened response is governed by the effective modulus K = I + χ, with susceptibility χ. In 2D, an isotropic K leaves only ν as the response parameter. The microscopic generator of the nonlinearity is the bucklable spring, whi

What would settle it

In the rectified phase (|F|=30), measure all three independent components of the far-field stress tensor as a function of polar angle, and simultaneously apply a small pure-shear perturbation to independently extract the shear modulus. The isotropic single-ν reduction predicts (i) all components fit one ν≈1 and (ii) the shear modulus vanishes; a deviation in either—directional anisotropy or a finite shear modulus—would falsify the central claim.

Watch

Extended reading notes

Core claim

The paper's central claim is that the full nonlinear response of the elastic network in all three mechanical regimes is captured by just a renormalization of the Poisson ratio, ν(|F|), which the authors extract by fitting simulation data to a vector-charge-theory (VCT) prediction for an isotropic elastic medium. The fit yields ν≈0.33 at dipole force 5, rising through 0.62 at 15, 0.90 at 25, and ≈1 at 30 (Table 1), with r²≥0.99 at each force. The authors show that the sign changes in the far-field tensor D, including the rectified phase in which D_yy and D_xy are opposite to the applied stress S, are consequences of this single parameter crossing the 2D stability limit, implying a vanishing s

Load-bearing premise

The entire reduction to one scalar Poisson ratio presumes that the far-field response is governed by an isotropic elastic modulus tensor, so that a single ν(|F|) captures the material; the paper's own buckling patterns are organized into anisotropic domains, and if the effective tensor is anisotropic, the fitted ν is a directional average and the implication that ν→1 means a vanishing shear modulus does not follow.

Editorial extensions

If this is right

  • In the rectified phase (|F|≈30), ν reaches ≈1, which in 2D implies the shear modulus μ→0; the paper interprets the contractile far-field response to an extensile dipole as the precursor to an activity-driven shear instability, plasticity, or flow.
  • Because ν is dimensionless and extracted from stress measurements alone, the VCT screening description applies to adaptable or renewable networks (such as cytoskeletal or extracellular-matrix assemblies) where strain and reference states are ill-defined.
  • The prediction of the full angular stress profile in every phase means that any future experiment measuring the far-field stress around a local force in a buckling network can be directly compared with a one-parameter theory; deviations would signal anisotropy or higher-order effects.
  • Disorder up to 10% bond dilution preserves the rectification curve, indicating the screening mechanism is robust to structural heterogeneity, not a lattice artifact.
  • Phase N corresponds to nearly complete screening (the far field shows almost no net response), while phase R corresponds to overscreening (the response has the opposite sign); this maps the three phases onto the dielectric analogy of a screening transition.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the isotropic single-parameter description is exact, then a second, weak force dipole placed far away should interact with the first through a screened potential that changes sign between phase N and phase R; this is a testable prediction for paired-dipole experiments.
  • The discontinuous jumps in ν at the phase boundaries, together with the ordered buckling domains, suggest a first-order transition in the effective elastic constants; cycling the force up and down should reveal hysteresis, which the paper's quasi-static annealing does not address.
  • A direct independent shear measurement on the network in phase R would settle the claim that ν→1 corresponds to μ→0; if the measured shear modulus remains finite, the isotropic mapping breaks down at large force.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper studies a triangular lattice of non-linear springs that can buckle at their midpoints, subject to a local extensile force dipole. Using energy minimization with periodic boundary conditions, it identifies three mechanical regimes in the far-field force-moment response: a linear regime (L), a nonlinear regime (N) with a sharp drop in Tr(D) and sign reversal of D_yy, and a rectified regime (R) where the far-field response is contractile to an extensile dipole. The paper connects these regimes to organized patterns of Kagome-lattice soft modes ('twisted units') that appear at the phase boundaries. It then invokes vector charge theory (VCT) of prestressed elasticity, generalizes it to a nonlinear dielectric with a field-dependent susceptibility, and claims that the full nonlinear response in all three regimes is captured by a single renormalized Poisson ratio ν(|F|) that rises from ~1/3 at small force to ~1 at |F|=30. The screening analogy is proposed as the mechanism for rectification. Supporting results include fitted ν values with high r², a linear-network baseline, disorder checks with 1% and 10% bond dilution, and visual comparisons of predicted and simulated stress fields.

Significance. If the main claim holds, the paper provides an unusually simple reduction: a strongly nonlinear, internally patterned fiber network responds at large scales as an isotropic elastic medium with a single force-dependent Poisson ratio, and the rectified phase is a precursor to shear instability. The empirical package is valuable: the simulations are clearly specified, the three-phase characterization is based on simultaneous fits to all three stress components with r² near 0.99, the linear-network control gives ν=1/3, and disorder robustness is checked. These are real strengths. However, the central quantitative claim is only as strong as the isotropy assumption on which the single-ν reduction rests, and that assumption is asserted rather than tested. The paper also explicitly acknowledges that the shear modulus cannot be deduced from the stress response, which qualifies the ν→1 instability interpretation. The significance is therefore conditional: the paper demonstrates a striking empirical correlation and a plausible theoretical framework, but it does not yet establish that the renormalized Poisson ratio is a well-defined material property of the effective medium.

major comments (3)
  1. [Major] The single-ν reduction is obtained by fitting the simulated stress field to the VCT Green function derived from an isotropic elastic modulus tensor, described as 'the simplest possible elastic modulus'. This is an assumption, not a test. The paper's own Fig. 4 shows 'distinctive anisotropic spatial patterns' in σ_yy and σ_θθ in all three phases. If the effective modulus tensor is anisotropic, the fitted ν is a directional average and the relation ν=(λ−μ)/(λ+μ) used to conclude ν→1 implies μ→0 does not follow. The reported r² values (Table 2) are in-sample and do not discriminate between the isotropic form and an anisotropic null model. I request that the authors fit an anisotropic modulus tensor (e.g., orthotropic or the full fourth-rank K) to the same data and compare goodness of fit, or otherwise provide a quantitative test of isotropy (e.g., angle-resolved residuals). Without this, th
  2. [Major] The 'predictions' in Fig. 5 are generated using the same values of ν that were fitted to the same simulated stress data from which they are compared. This is an in-sample check, not an out-of-sample validation. The high r² in Table 2 shows only that the isotropic form can represent the data within the fitting range, not that the fitted ν has predictive power. To validate the reduction, the authors should fit ν on one subset of the data (e.g., one stress component, one angular sector, or one q-range) and predict the remaining components, or use cross-validation. As written, the comparison in Fig. 5 does not provide independent evidence for the central claim.
  3. [Major] The phase boundaries at |F|=15 and |F|=26 are described as discontinuous, and the buckling patterns at these forces show large ordered domains. However, the simulation protocol is a quasi-static ramp with monotonically increasing |F|, starting from an unperturbed lattice. No hysteresis or protocol dependence is reported. Because the transitions are associated with self-organized domain patterns, metastability is a real possibility: the system might remain in a metastable branch when |F| is increased, and the sharp changes in Tr(D) and ν(F) could be artifacts of the annealing direction. I request checks with decreasing |F|, different Δ|F|, and/or multiple random initial configurations. This is particularly important because the claim of 'discontinuous' changes underlies the phase classification.
minor comments (6)
  1. [Minor] Typo: 'reponse' should be 'response' in the sentence 'the rectification is a manifestation of nonlinear elastic screening'.
  2. [Minor] The notation '|F|_0 = Δ|F|, and Δ|F| = 2κ/l_0' is unclear: it suggests the increment equals 2κ/l_0, but the text earlier says |F|_0 = Δ|F|. Please define the initial force and the increment explicitly.
  3. [Minor] The value ν≈1 at |F|=30 is reported without uncertainty; the fitting constrains ν to [−1,1] in 2D, so a fit at the boundary is likely censored. Reporting a confidence interval or the fact that the optimizer hit the bound would be informative.
  4. [Minor] The caption says 'semi quantitatively' for the agreement; the main text says 'reproduce the stress patterns ... semi quantitatively'. Please specify what 'semi-quantitative' means (e.g., relative error, angular range, q-range) and state the limitations explicitly.
  5. [Minor] The caption uses 'Fig 1B and fig 3C' with inconsistent abbreviations; please standardize. Also, the color scale for the stress fields is not defined.
  6. [Minor] The statement 'we cannot deduce the shear modulus from the stress response' is in tension with the earlier claim that ν→1 implies μ→0. Please clarify that the μ→0 inference relies on the isotropic relation and is therefore conditional on the isotropy assumption.

Circularity Check

1 steps flagged · score 6.0 of 10

The central ν(|F|) renormalization is a fitted parameter, and the Fig. 5 'predictions' re-evaluate the same fit on the same data; the isotropy assumption that makes ν the only parameter is asserted, not tested.

  1. fitted input called prediction [Table 1; Fig. 5; 'Fitting Procedure' and 'Comparing simulations to VCT predictions' (Supplemental); Eq. (4)]
    "We extract the Poisson ratio (ν(|F|)) by fitting the simulation results for the stress response to the theoretical predictions. Table 1 shows the values of the Poisson ratio as a function of|F|. ... In fig 5, we test the validity of this approach by comparing the stress responses predicted using the fitted values of ν(|F|)to the simulated ones. ... The full nonlinear response of the elastic network in all three mechanical regimes, L (linear), N (non-linear), and R (rectified), is captured by just a renormalization of the Poisson ratio."

    ν(|F|) is not derived from the buckling patterns or from a parameter-free theory; it is the best-fit parameter of the VCT isotropic Green function (Eq. 4) to the very simulated stress components whose 'prediction' it is then used to demonstrate. The fitting procedure states: 'The Python package SymFit is used to simultaneously fit all components of σ_ij(q) to the same value of ν.' Fig. 5 therefore re-evaluates the fitted model on the same dataset; no held-out data or independent prediction is involved. 'Captured by just a renormalization of the Poisson ratio' is thus a restatement of the fit, not a test that the framework predicted the response.

full rationale

The paper is transparent that ν is a fitted quantity, and the simulations are the primary data product. The circularity is in the validation language: Fig. 5 is called a test/prediction, but it uses the same fitted ν on the same stress data used to obtain ν. Since the central claim ('full nonlinear response ... captured by just a renormalization of the Poisson ratio') is literally the fitted parameter, that claim reduces to the fit by construction. There is residual non-circular content: the reported r² values (0.9887–0.9948) show that an assumed isotropic, single-ν form can represent all three stress components simultaneously, and the buckling-pattern correlation is a separate observation. However, the isotropy of K is assumed ('Assuming an isotropic form for K'), not demonstrated against an anisotropic modulus tensor; the paper's own Fig. 4 shows anisotropic stress patterns, and the paper admits 'we cannot deduce the shear modulus from the stress response,' making the ν→1 ⇒ μ→0 instability interpretation conditional on the assumed isotropic relation. I do not count the VCT self-citations as independently circular: VCT is an external published formalism being tested against simulation here, and author overlap is not itself the load-bearing step. Overall: one central in-sample 'prediction' reduces by construction, so score 6.

Assumptions & free parameters 3 free parameters · 5 assumptions · 2 invented entities

The quantitative payload reduces to one fitted parameter ν(|F|) plus the hand-chosen κ/µ. The framework (VCT), the D↔P bridge, and the 'response depends only on ν' reduction are imported from prior work, partly the authors' own; the new diagnostics (Kagome twisted units) are presented without falsifiable statistics. This makes the ledger lean: the paper is an interpretive application of an established theory to one lattice model.

free parameters (3)
  • Renormalized Poisson ratio ν(|F|) = 0.33, 0.34, 0.62, 0.72, 0.90, 1.00 at |F|=5,10,15,20,25,30; fit constrained to [−1,1]; 1.00 at F=30 is at the constraint
    The single free parameter of the VCT Green function fit (Table 1; Fitting Procedure). The paper's central quantitative claim is literally the fitted values.
  • Buckling/stretch energy ratio κ/µ = 1/1000
    Hand-chosen so buckling is the preferred mode; all reported phases and ν values are conditional on this ratio (Eq. 1).
  • Fitting q-range and angular binning = unspecified q-range; 30 angular bins
    Fits are restricted to q where σ_ij(q) ∼ sin(d q_x/2)/q already holds; the range boundaries are not stated, so r²≈0.99 is partially a selection outcome.
assumptions (5)
  • domain assumption In 2D, the VCT (prestressed-elasticity) response of an isotropic medium depends only on the Poisson ratio ν, not on the shear modulus
    Foundation of the single-parameter fit; adopted from the authors' prior VCT papers (Refs 18, 26, 27), §'The macroscopic measures...' ('In 2D, the response depends only on ν 18').
  • domain assumption The measured far-field force-moment D equals a boundary integral of the polarization: D_αβ = −∮ n_λ P_αλ r_β
    Imported from Ronceray & Lenz (Ref 23); this is the bridge that gives the D-measurements their 'screening' meaning. Stated, not derived, in §'The macroscopic measures...'.
  • ad hoc to paper The buckled network's far-field response is governed by an isotropic elastic modulus tensor K, so a single scalar ν(|F|) can describe all three stress components simultaneously
    Chosen for simplicity ('the simplest possible elastic modulus') in the Fitting Procedure; contradicted in spirit by the anisotropic buckling domains (Fig 3C, 4).
  • domain assumption The stepwise energy-minimization ('annealing') protocol reaches the mechanically relevant equilibrium; results are independent of step size and initial state (no hysteresis, no metastable trapping)
    Supplemental 'Simulation details': F increased by ∆F from the previous minimum; discontinuous transitions at |F|=15, 26 are claimed without protocol checks.
  • ad hoc to paper Kagome-lattice soft modes (twisted units) are the correct basis for reading out organized buckling
    Supplemental 'Kagome construction': the Kagome network 'is a theoretical construct... there are no springs connecting the Kagome points'; no derivation that these modes carry the response.
invented entities (2)
  • Twisted units (chirality-paired Kagome triangles)
    purpose: Diagnostic to visualize organized buckling and to correlate domain structure with the L/N/R phase transitions
    Introduced in 'Kagome construction'; correlation with phases is asserted from visually selected force values (Fig 1B, 3C) without a quantitative domain-size or order-parameter statistic.
  • Mechanical screening (force-dipole screening via buckling, with the polarization field P)
    purpose: Explains rectification as dielectric-analog screening of the extensile dipole, the paper's headline interpretation
    The VCT polarization formalism predates this paper (Refs 18, 26, 27); what is new is its application to this network. The only quantitative support is the in-sample ν fits; the polarization interpretation of the measured D is inherited from Ref 23.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Rectification of stress by fiber networks: Manifestation of non-linear screening through self-organized buckling." pith.science (2026). https://pith.science/paper/7EJ6XQ6Z

@misc{pith2026251120754,
  author       = {Pith},
  title        = {Pith review of: Rectification of stress by fiber networks: Manifestation of non-linear screening through self-organized buckling},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/7EJ6XQ6Z}},
  note         = {Machine review of arXiv:2511.20754}
}
read the original abstract

Force transmission at large length scales is crucial for such biological functions as cell motility and morphogenesis. The networks that transmit these forces are malleable, patterned by active forces generated at the microscale by biological motors. In this paper we explore a simple model of a non-linear fiber network which has only two modes of deformation, but exhibits diverse mechanical phases with distinct large-scale response, tuned by the strength of a microscopic force dipole. We demonstrate, via numerical simulations, that the network is remodeled by organized patterns of buckling, which lead to a renormalization of the Poisson ratio. Finally, we show that the emergent behavior at large length scales can be ascribed to "mechanical screening" of the force dipole, analogous to dielectric screening of charges in electrostatics.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

30 extracted references

  1. [1]

    Abhishek, A

    M. Abhishek, A. Dhanuka, D. S. Banerjee and M. Rao, Phys. Rev. E, 2025, 112, 045433

    2025

  2. [2]

    Roychowdhury, S

    A. Roychowdhury, S. Dasgupta and M. Rao, Phys. Rev. E, 2025, 112, 044414

    2025

  3. [3]

    D. A. Gagnon, C. Dessi, J. P. Berezney, R. Boros, D. T.-N. Chen, Z. Dogic and D. L. Blair, Physical Review Letters, 2020, 125, 178003

    2020

  4. [4]

    C. P. Broedersz and F. C. MacKintosh, Soft Matter, 2011, 7, 3186--3191

    2011

  5. [5]

    Broedersz and F

    C. Broedersz and F. MacKintosh, Reviews of Modern Physics, 2014, 86, 995--1036

    2014

  6. [6]

    Fürthauer, D

    S. Fürthauer, D. J. Needleman and M. J. Shelley, New Journal of Physics, 2021, 23, 013012

    2021

  7. [7]

    Wyse Jackson, J

    T. Wyse Jackson, J. Michel, P. Lwin, L. A. Fortier, M. Das, L. J. Bonassar and I. Cohen, Science Advances, 2022, 8, eabk2805

    2022

  8. [8]

    Xu and S

    X. Xu and S. A. Safran, Phys. Rev. E, 2015, 92, 032728

    2015

Show all 30 references
  1. [9]

    Murrell and M

    M. Murrell and M. Gardel, Proc. Natl. Acad. Sci, 2012, 109, 20820--20825

    2012

  2. [10]

    H. Wang, A. Abhilash, C. S. Chen, R. G. Wells and V. B. Shenoy, Biophysical Journal, 2014, 107, 2592--2603

    2014

  3. [11]

    Najma, W

    B. Najma, W. Wei, A. Baskaran, P. Foster and G. Duclos, Proc. Natl. Acad. Sci., 2024, 121, year

    2024

  4. [12]

    Zakharov, M

    A. Zakharov, M. Awan, A. Gopinath, S. Lee, A. Ramasubramanian and K. Dasbiswas, Science Advances, 2024, 10, year

    2024

  5. [13]

    Farhadi, S

    L. Farhadi, S. N. Ricketts, M. J. Rust, M. Das, R. M. Robertson-Anderson and J. L. Ross, Soft Matter, 2020, 16, 7191--7201

    2020

  6. [14]

    Burla, J

    F. Burla, J. Tauber, S. Dussi, J. van der Gucht and G. H. Koenderink, Nature Physics, 2019, 15, 549--553

    2019

  7. [15]

    Fürthauer, B

    S. Fürthauer, B. Lemma, P. J. Foster, S. C. Ems-McClung, C.-H. Yu, C. E. Walczak, Z. Dogic, D. J. Needleman and M. J. Shelley, Nature Physics, 2019, 15, 1295--1300

    2019

  8. [16]

    Tanaka and S

    F. Tanaka and S. F. Edwards, Macromolecules, 1992, 25, 1516--1523

    1992

  9. [17]

    Ronceray, C

    P. Ronceray, C. P. Broedersz and M. Lenz, Proceedings of the National Academy of Sciences, 2016, 113, 2827--2832

    2016

  10. [18]

    J. N. Nampoothiri, M. D'Eon, K. Ramola, B. Chakraborty and S. Bhattacharjee, Physical Review E, 2022, 106, 065004

    2022

  11. [19]

    T. C. Lubensky, C. L. Kane, X. Mao, A. Souslov and K. Sun, Reports on Progress in Physics, 2015, 78, 073901

    2015

  12. [20]

    Mao and T

    X. Mao and T. C. Lubensky, Annual Review of Condensed Matter Physics, 2018, 9, 413--433

    2018

  13. [21]

    Malakar, A

    K. Malakar, A. Countryman and B. Chakraborty, Supplemental Information, 2025

    2025

  14. [22]

    Benoist, G

    F. Benoist, G. Saggiorato and M. Lenz, Soft Matter, 2023, 19, 2970--2976

    2023

  15. [23]

    Ronceray and M

    P. Ronceray and M. Lenz, Soft Matter, 2015, 11, 1597--1605

    2015

  16. [24]

    X. Mao, N. Xu and T. C. Lubensky, Physical Review Letters, 2010, 104, 085504

    2010

  17. [25]

    Zangwill, Modern Electrodynamics, Cambridge University Press, 2013

    A. Zangwill, Modern Electrodynamics, Cambridge University Press, 2013

    2013

  18. [26]

    Chakraborty, J

    S. Chakraborty, J. N. Nampoothiri, S. Bhattacharjee, B. Chakraborty and K. Ramola, Stress Response of Jammed Solids: Prestress and Screening, 2025, https://arxiv.org/abs/2509.14336

    2025

  19. [27]

    J. N. Nampoothiri, Y. Wang, K. Ramola, J. Zhang, S. Bhattacharjee and B. Chakraborty, Physical Review Letters, 2020, 125, 118002

    2020

  20. [28]

    T. C. T. Ting and T. Chen, The Quarterly Journal of Mechanics and Applied Mathematics, 2005, 58, 73--82

    2005

  21. [29]

    J. L. Shivers, S. Arzash and F. MacKintosh, Physical Review Letters, 2020, 124, 038002

    2020

  22. [30]

    k P 5h572|aDph88 „ j2oJj*muozo L\|Â

    J. Berezney, B. L. Goode, S. Fraden and Z. Dogic, Proceedings of the National Academy of Sciences, 2022, 119, e2115895119 mcitethebibliography communicationsArXiv.tex0000664000000000000000000016635015111372745014323 0ustar rootroot [twoside,twocolumn,9pt] article extsizes [sup...

    2022

Pith tools

Reviewed August 3, 2026 · model on record in the stance chip above.