pith. sign in

arxiv: 2604.23642 · v1 · submitted 2026-04-26 · ❄️ cond-mat.soft · cond-mat.mtrl-sci

Constitutive relations for colloidal gel

Saikat Roy , Yezaz Ahmed Gadi Man This is my paper

Pith reviewed 2026-05-08 05:05 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.mtrl-sci
keywords colloidal gelsconstitutive relationsdepletion gelsfrictional gelsfree energy expansionstress-strain responseamorphous materialscontinuum mechanics
0
0 comments X

The pith

Colloidal gels require free-energy expansions around stressed reference states to derive accurate constitutive relations, unlike traditional stress-free models.

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

The paper establishes that standard continuum theories for depletion gels expand free energy around a stress-free reference state and assume affine deformations, but these premises fail for amorphous gels that never reach global energy minima and instead contain local compressive, tensile, and zero-stress regions. Expansion around a stressed reference state therefore generates qualitatively different scalar terms in the harmonic expansion, leading the authors to propose simpler constitutive relations that relate global stress and strain more faithfully. Large-scale simulations of both depletion and frictional gels over wide parameter ranges confirm that the new relations capture mechanical response where prior models do not. A sympathetic reader cares because reliable constitutive laws are needed to predict how soft, disordered materials deform, flow, and yield in practical settings.

Core claim

Because gels contain local regions of significant compressive and tensile stress interspersed with zero-stress regions and never attain a true global minimum in the potential energy landscape, the free-energy expansion around a stressed reference state produces scalar terms in the harmonic expansion whose effects differ from those obtained around an unstressed state; the authors therefore replace traditional continuum constitutive relations with simple new ones and demonstrate their robustness through extensive numerical simulations.

What carries the argument

Harmonic expansion of free energy around a stressed reference state, which supplies the distinct scalar contributions that enter the proposed constitutive relations for gel stress-strain response.

If this is right

  • The new constitutive relations correctly reproduce the mechanical response of both depletion and frictional gels over wide ranges of particle interactions and densities.
  • Local stress heterogeneity and non-affine particle motion are directly incorporated into the global stress-strain law without additional fitting parameters.
  • Qualitatively different harmonic terms appear in the free-energy expansion, altering predicted moduli, yield stresses, and strain responses compared with prior models.
  • The approach remains robust when applied to gels prepared under different interaction strengths and friction coefficients.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same stressed-state expansion could be tested on other amorphous soft solids such as glasses or emulsions where local stresses are known to persist.
  • Direct imaging of local stress distributions inside gels could be used to calibrate the scalar coefficients in the new relations.
  • Preparation protocols that reduce or increase internal stress heterogeneity might be predicted to shift the measured constitutive curves in specific ways.
  • The framework suggests that rate-dependent or history-dependent effects in gels may also require expansion around non-equilibrium stressed states.

Load-bearing premise

Gels never reach a true global minimum in the potential energy landscape and therefore always contain local regions of significant compressive and tensile stress interspersed with zero-stress regions.

What would settle it

A large-scale simulation or rheology experiment in which the stress-strain curves of depletion or frictional gels match the predictions of traditional stress-free constitutive relations more closely than the stressed-state relations proposed here, across the same parameter space.

Figures

Figures reproduced from arXiv: 2604.23642 by Saikat Roy, Yezaz Ahmed Gadi Man.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Contact orientation distribution is plotted on a view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. (a) Evolution of contact anisotropy, view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Comparison between the prediction of the proposed view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. (a) view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. (a) Illustration of typical depletion gel structure at view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Contact orientation probability density function, view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Directional variation of the different components view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Angular variation of the pairwise central interaction view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Angular variation of view at source ↗
read the original abstract

The theoretical treatment of depletion gels with central interactions often involves expanding the free energy around a stress-free reference state to derive a constitutive relation between global stress and strain. The premise upon which the previous continuum theories are based, i.e., the stress-free reference state and the affine deformation, both of which do not hold in the context of amorphous gel materials. Gels never reach a true global minimum in the potential energy landscape and contain local regions of significant compressive and tensile stress, interspersed with zero-stress regions. Hence, expansion of free energy around a stressed reference state will produce scalar terms in harmonic expansion, the effects of which are qualitatively different from the terms appearing in the expansion around an unstressed reference state. In this study, we demonstrate the limitations of traditional continuum theories and propose simple constitutive relations that better capture the mechanical response of gel materials. The robustness of the proposed relations is established through large-scale numerical simulations of depletion and frictional gels across a vast parameter space.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 2 minor

Summary. The paper argues that traditional continuum theories for colloidal gels rely on free-energy expansions around a stress-free reference state and affine deformations, both of which fail for amorphous gels that never reach global energy minima and contain interspersed regions of compressive, tensile, and zero stress. It proposes new simple constitutive relations obtained by expanding around a stressed reference state, which are claimed to generate qualitatively different scalar terms, and asserts that these relations better capture gel mechanics as demonstrated by large-scale simulations of depletion and frictional gels over a broad parameter space.

Significance. If the new relations can be derived independently, shown to introduce necessary terms not capturable by adjustable parameters in existing frameworks, and validated with quantitative metrics, the work would address a genuine limitation in modeling non-affine response in soft amorphous materials. The focus on prestressed local states offers a potentially useful conceptual shift for colloidal gels and related systems, though its impact depends on explicit demonstration that the claimed qualitative differences improve predictive power.

major comments (3)
  1. Abstract and theoretical section: the central claim that expansion around a stressed reference state produces scalar terms qualitatively different from those in the stress-free case is asserted without any explicit derivation, expansion steps, or identification of the additional invariants arising from the prestress tensor; no equations are provided showing how the local stress heterogeneity alters the harmonic terms.
  2. Validation section: the robustness of the proposed relations is established via simulations, yet no quantitative comparison metrics (e.g., R², error norms, or cross-validation scores), error analysis, or direct head-to-head fits against traditional models are reported; this leaves open whether the improvement is independent or partly by construction.
  3. Simulation analysis: the premise that gels contain local stressed and zero-stress regions requiring a stressed reference (rather than averaging to an effective stress-free state at continuum scale) is not isolated or tested; no evidence is given that zero-stress regions do not dominate or that the new terms are necessary rather than reparametrizations of existing constitutive forms.
minor comments (2)
  1. Clarify notation for any new scalar invariants or terms introduced in the constitutive relations to allow direct comparison with classical affine models.
  2. Include tables or figures with explicit fit statistics and parameter counts for both traditional and proposed relations across the simulated parameter space.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the constructive and detailed report. We address each major comment below and outline the revisions we will make to the manuscript.

read point-by-point responses

  1. Referee: Abstract and theoretical section: the central claim that expansion around a stressed reference state produces scalar terms qualitatively different from those in the stress-free case is asserted without any explicit derivation, expansion steps, or identification of the additional invariants arising from the prestress tensor; no equations are provided showing how the local stress heterogeneity alters the harmonic terms.

    Authors: We agree that the theoretical section requires an explicit derivation to substantiate the central claim. In the revised manuscript we will insert a step-by-step expansion of the free-energy functional around a prestressed reference state, explicitly identifying the additional invariants generated by the prestress tensor and demonstrating how local stress heterogeneity modifies the harmonic terms relative to the stress-free expansion. revision: yes

  2. Referee: Validation section: the robustness of the proposed relations is established via simulations, yet no quantitative comparison metrics (e.g., R², error norms, or cross-validation scores), error analysis, or direct head-to-head fits against traditional models are reported; this leaves open whether the improvement is independent or partly by construction.

    Authors: We acknowledge that quantitative metrics are needed to demonstrate the improvement. We will add R² values, error norms, cross-validation scores, and direct head-to-head comparisons against traditional constitutive models in the revised validation section, thereby clarifying that the reported gains are independent of parameter adjustment. revision: yes

  3. Referee: Simulation analysis: the premise that gels contain local stressed and zero-stress regions requiring a stressed reference (rather than averaging to an effective stress-free state at continuum scale) is not isolated or tested; no evidence is given that zero-stress regions do not dominate or that the new terms are necessary rather than reparametrizations of existing constitutive forms.

    Authors: Our simulations of depletion and frictional gels across a wide parameter space already show that the new relations capture mechanical responses that traditional models cannot reproduce, consistent with the presence of mixed local stress states. To isolate this effect we will add local stress-distribution analysis and conditional averaging in the revised manuscript, demonstrating that zero-stress regions do not dominate the continuum response and that the new scalar terms are not simple reparametrizations. revision: yes

Circularity Check

0 steps flagged

No significant circularity detected; derivation remains self-contained.

full rationale

The paper asserts that gels contain local stressed regions (unlike stress-free affine models), motivating expansion around a stressed reference state that yields qualitatively different scalar terms, then proposes simple constitutive relations validated on independent large-scale simulations of depletion and frictional gels across wide parameter space. No equations or text in the provided material reduce the proposed relations to simulation outputs by construction, nor do any self-citations or fitted inputs appear as load-bearing steps. The simulations function as external falsification rather than parameter tuning that forces the claimed improvement.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the domain assumption that amorphous gels possess persistent local stresses and lack a global stress-free minimum; no free parameters or new entities are introduced in the abstract.

axioms (2)
  • domain assumption Gels never reach a true global minimum in the potential energy landscape
    Stated explicitly in the abstract as the reason traditional expansions fail.
  • domain assumption Gels contain local regions of significant compressive and tensile stress, interspersed with zero-stress regions
    Invoked to argue that expansion around a stressed reference state is necessary.

pith-pipeline@v0.9.0 · 5461 in / 1440 out tokens · 74655 ms · 2026-05-08T05:05:53.003121+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages · 1 internal anchor

  1. [1]

    hint towards a possible discrepancy between the measured data and the theoretical prediction [3]: upon increasing the dimensionless attraction strength, a sevenfold rise is observed for the shear modulus whereas the product of the cluster number density,N c, and coordination number,z c displays only a threefold rise. Gels also exhibit intriguing power law...

    work page internal anchor Pith review Pith/arXiv arXiv 2026

  2. [2]

    fails miserably to predict the power-law growth of stress and the evolution of NSD. We also propose simple constitutive relations that relate the components of the stress tensor to the appropriate internal state parameters, correctly encoding the non-linear effects of the quenched stresses. The generality of the proposed re- lations is rigorously tested t...

    work page 2023

  3. [3]

    A linear correc- tion term is also introduced to the potential to ensure that both the potential energy and force become zero smoothly at the cutoff

    r≥r c (12) where,β= 1/k BT,βU 0 denotes the dimensionless strength of the potential minimum andκ −1 dictates 6 the range of attraction.r 0 =a i +a j represents the equilibrium separation distance between a particle pair, wherea i anda j are the radii of thei th and thej th particles.r c sets the potential cutoff. A linear correc- tion term is also introdu...

    work page

  4. [4]

    J. N. Nampoothiri, M. D’Eon, K. Ramola, FIG. 10. Directional variation of the different components (fn(θ, ϕ), ftθ(θ, ϕ), ftϕ(θ, ϕ)) of the pairwise interaction force normalized by the mean normal force over all interactions is shown for the frictional gel system across different volume fractions. Here,βU 0 = 50. FIG. 11. Angular variation of fn(θ, ϕ), ftθ...

    work page 2022

  5. [5]

    Vinutha, F

    H. Vinutha, F. D. Diaz Ruiz, X. Mao, B. Chakraborty, and E. Del Gado, Stress–stress correlations reveal force chains in gels, The Journal of chemical physics158 (2023)

    work page 2023

  6. [6]

    Zaccone, H

    A. Zaccone, H. Wu, and E. Del Gado, Elasticity of ar- rested short-ranged attractive colloids:¡? format?¿ homo- geneous and heterogeneous glasses, Physical review let- 10 ters103, 208301 (2009)

    work page 2009

  7. [7]

    Ramakrishnan, Y.-L

    S. Ramakrishnan, Y.-L. Chen, K. Schweizer, and C. Zukoski, Elasticity and clustering in concentrated de- pletion gels, Physical Review E—Statistical, Nonlinear, and Soft Matter Physics70, 040401 (2004)

    work page 2004

  8. [8]

    K. A. Whitaker, Z. Varga, L. C. Hsiao, M. J. Solomon, J. W. Swan, and E. M. Furst, Colloidal gel elasticity arises from the packing of locally glassy clusters, Nature communications10, 2237 (2019)

    work page 2019

  9. [9]

    Roy and M

    S. Roy and M. S. Tirumkudulu, Yielding in a strongly aggregated colloidal gel. part i: 2d simulations, Journal of Rheology60, 559 (2016)

    work page 2016

  10. [10]

    Roy and M

    S. Roy and M. S. Tirumkudulu, Yielding in a strongly aggregated colloidal gel. part ii: Theory, Journal of Rhe- ology60, 575 (2016)

    work page 2016

  11. [11]

    Roy and M

    S. Roy and M. S. Tirumkudulu, Universality in consoli- dation of colloidal gels, Soft matter12, 9402 (2016)

    work page 2016

  12. [12]

    Roy and M

    S. Roy and M. S. Tirumkudulu, Micro-mechanical theory of shear yield stress for strongly flocculated colloidal gel, Soft Matter16, 1801 (2020)

    work page 2020

  13. [13]

    Buscall and L

    R. Buscall and L. R. White, The consolidation of concen- trated suspensions. part 1.—the theory of sedimentation, Journal of the Chemical Society, Faraday Transactions 1: Physical Chemistry in Condensed Phases83, 873 (1987)

    work page 1987

  14. [14]

    Tateno, Y

    M. Tateno, Y. Wang, and H. Tanaka, Mechanical self- organization of particle networks during uniaxial com- pression yielding, Physical Review X14, 011035 (2024)

    work page 2024

  15. [15]

    Tateno, Y

    M. Tateno, Y. Wang, and H. Tanaka, Void connectiv- ity and criticality in the compression-induced gel-to-glass transition of short-range attractive colloids, Physical Re- view Letters134, 048201 (2025)

    work page 2025

  16. [16]

    Y. A. Gadi Man, D. S. Dagur, and S. Roy, Origin of normal stress difference in consolidating strong depletion gels, The Journal of Chemical Physics163(2025)

    work page 2025

  17. [17]

    Y. A. Gadi Man, D. S. Dagur, and S. Roy, Micro- mechanical insights into the stress transmission in strongly aggregating colloidal gel, The Journal of Chem- ical Physics158(2023)

    work page 2023

  18. [18]

    M. M. Islam and D. R. Lester, Normal stress differences in the consolidation of strong colloidal gels, Rheologica Acta60, 59 (2021)

    work page 2021

  19. [19]

    R. Seto, R. Botet, M. Meireles, G. K. Auernhammer, and B. Cabane, Compressive consolidation of strongly aggregated particle gels, Journal of rheology57, 1347 (2013)

    work page 2013

  20. [20]

    C. J. Dibble, M. Kogan, and M. J. Solomon, Structure and dynamics of colloidal depletion gels: Coincidence of transitions and heterogeneity, Physical Review E74, 041403 (2006)

    work page 2006

  21. [21]

    P. A. Cundall and O. D. Strack, A discrete numeri- cal model for granular assemblies, geotechnique29, 47 (1979)

    work page 1979

  22. [22]

    See supplemental material for additional information

    work page

  23. [23]

    L. D. Landau, E. M. Lifˇ sic, E. M. Lifshitz, A. M. Kose- vich, and L. P. Pitaevskii,Theory of elasticity: volume 7, Vol. 7 (Elsevier, 1986)

    work page 1986

  24. [24]

    Ken-Ichi, Distribution of directional data and fabric tensors, International Journal of Engineering Science22, 149 (1984)

    K. Ken-Ichi, Distribution of directional data and fabric tensors, International Journal of Engineering Science22, 149 (1984)

    work page 1984

  25. [25]

    Plimpton, Fast parallel algorithms for short-range molecular dynamics, Journal of computational physics 117, 1 (1995)

    S. Plimpton, Fast parallel algorithms for short-range molecular dynamics, Journal of computational physics 117, 1 (1995)

    work page 1995

  26. [26]

    Varga and J

    Z. Varga and J. Swan, Hydrodynamic interactions en- hance gelation in dispersions of colloids with short-ranged attraction and long-ranged repulsion, Soft matter12, 7670 (2016)

    work page 2016

  27. [27]

    De Graaf, W

    J. De Graaf, W. C. Poon, M. J. Haughey, and M. Hermes, Hydrodynamics strongly affect the dynamics of colloidal gelation but not gel structure, Soft Matter15, 10 (2019)

    work page 2019

  28. [28]

    C. P. Royall, M. A. Faers, S. L. Fussell, and J. E. Hal- lett, Real space analysis of colloidal gels: Triumphs, chal- lenges and future directions, Journal of Physics: Con- densed Matter33, 453002 (2021)

    work page 2021

  29. [29]

    ´ aP Meeker, R

    S. ´ aP Meeker, R. ´ a. L. Evans, M. ´ aM Robins,et al., De- layed sedimentation of transient gels in colloid–polymer mixtures: dark-field observation, rheology and dynamic light scattering studies, Faraday Discussions112, 143 (1999)

    work page 1999

  30. [30]

    Shah, Y.-L

    S. Shah, Y.-L. Chen, K. Schweizer, and C. Zukoski, Vis- coelasticity and rheology of depletion flocculated gels and fluids, The Journal of chemical physics119, 8747 (2003)

    work page 2003

  31. [31]

    Buscall, P

    R. Buscall, P. D. Mills, J. W. Goodwin, and D. Lawson, Scaling behaviour of the rheology of aggregate networks formed from colloidal particles, Journal of the Chemical Society, Faraday Transactions 1: Physical Chemistry in Condensed Phases84, 4249 (1988)

    work page 1988

  32. [32]

    Buscall, P

    R. Buscall, P. Mills, R. Stewart, D. Sutton, L. White, and G. Yates, The rheology of strongly-flocculated sus- pensions, Journal of Non-Newtonian Fluid Mechanics24, 183 (1987)

    work page 1987

  33. [33]

    Bonacci, X

    F. Bonacci, X. Chateau, E. M. Furst, J. Fusier, J. Goyon, and A. Lemaˆ ıtre, Contact and macroscopic ageing in col- loidal suspensions, Nature Materials19, 775 (2020)

    work page 2020

  34. [34]

    Bonacci, X

    F. Bonacci, X. Chateau, E. M. Furst, J. Goyon, and A. Lemaˆ ıtre, Yield stress aging in attractive colloidal sus- pensions, Physical Review Letters128, 018003 (2022)

    work page 2022

  35. [35]

    R. Seto, R. Mari, J. F. Morris, and M. M. Denn, Discon- tinuous shear thickening of frictional hard-sphere suspen- sions, Physical review letters111, 218301 (2013)

    work page 2013

  36. [36]

    C.-P. Hsu, S. N. Ramakrishna, M. Zanini, N. D. Spencer, and L. Isa, Roughness-dependent tribology effects on dis- continuous shear thickening, Proceedings of the National Academy of Sciences115, 5117 (2018)

    work page 2018

  37. [37]

    R. Mari, R. Seto, J. F. Morris, and M. M. Denn, Discon- tinuous shear thickening in brownian suspensions by dy- namic simulation, Proceedings of the National Academy of Sciences112, 15326 (2015)

    work page 2015

  38. [38]

    N. Y. Lin, B. M. Guy, M. Hermes, C. Ness, J. Sun, W. C. Poon, and I. Cohen, Hydrodynamic and contact contri- butions to continuous shear thickening in colloidal sus- pensions, Physical review letters115, 228304 (2015)

    work page 2015

  39. [39]

    Singh, C

    A. Singh, C. Ness, R. Seto, J. J. de Pablo, and H. M. Jaeger, Shear thickening and jamming of dense suspen- sions: the “roll” of friction, Physical Review Letters124, 248005 (2020)

    work page 2020

  40. [40]

    Silbert, D

    L. Silbert, D. Erta¸ s, G. Grest, T. Halsey, D. Levine, and S. Plimpton, Granular flow down an inclined plane: Bagnold scaling and rheology, Phys. Rev. E64, 051302 (2001)

    work page 2001

  41. [41]

    Luding, Cohesive, frictional powders: contact models for tension, Granular matter10, 235 (2008)

    S. Luding, Cohesive, frictional powders: contact models for tension, Granular matter10, 235 (2008)

    work page 2008