REVIEW 3 major objections 4 minor 44 references
Free Growth under Tension
T0 review · 3 major / 4 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read A cell colony can develop tensile stress with no motility force at all, driven by growth plus substrate friction.
desk verdict A genuinely new mechanism for tensile stress in growing colonies without motility, with a clean analytic model and honest calibration; the expansion-speed prediction is partly self-consistency but the existence claim holds. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the growth law $k = \kappa(P_H - P) + \Delta k\Delta x\,\delta(x - x_0)$, which couples proliferation to mechanical stress through a bulk response coefficient $\kappa$ and a surface term $\Delta k\Delta x$: cells near the free boundary proliferate faster because growth there is mechanically cheaper. From that law, the continuity equation $\partial_x v = k$, and force balance with linear substrate friction, the argument reduces the stress field to the screened Poisson equation $\partial_x^2 P_x = (P_x - P_H)/\lambda^2$, whose solution is the hyperbolic-cosine pressure profile of Eq. (7). The same growth law, with a motility traction term added to the force balance, yields the motile pressure profile and the generalized expansion speed. The crucial entity is therefore the combination of negative homeostatic pressure and rim growth: the rim produces an inward flux, friction converts that flux into tension, and the tension suppresses growth in the bulk.
What would settle it
Measure the velocity field inside a non-motile expanding colony that is under tension: the model predicts a retrograde inward flow of cells throughout the interior, so the absence of such a flow would falsify the mechanism.
Extended reading notes
Core claim
The central discovery is that tension arises from the interplay of a negative homeostatic pressure and boundary-enhanced growth. The growth rate is written as $k = \kappa(P_H - P) + \Delta k \Delta x\,\delta(x - x_0)$, with the surface term representing faster growth in a thin layer at the free edge. Combining this with force balance, where substrate friction is the only external force, $\partial_x\sigma_{xx} = 2\rho\gamma v$, gives a screened Poisson equation for the pressure, $P_x = P_H\bigl(1 - \cosh(x/\lambda)/\cosh(x_0/\lambda)\bigr)$ with $\lambda^2 = (\rho\kappa\gamma)^{-1}$. For homeostatic pressures between a critical negative value $P_H^c = -\Delta k\Delta x\sqrt{4\rho\gamma/\kappa}$ and zero, colonies expand indefinitely at constant speed $v(x_0) = P_H\sqrt{\kappa/(4\rho\gamma)} + \Delta k\Delta x$ while their interior is tensile; the tension is balanced by friction against an inward retrograde flow of cells. Motility adds a boundary-localized traction that generates additional tension and extends the tensile-growth phase, and the two-dimensional circular case reduces to the same physics with growth proportional to the perimeter.
Load-bearing premise
The whole mechanism depends on the assumption that cells in a thin rim at the colony boundary proliferate faster than cells in the bulk; if that surface growth advantage disappeared, no inward flux would form and tension would not arise.
Editorial extensions
If this is right
- Non-motile tensile colonies should show a retrograde flow: cells move inward from the proliferating rim toward the center, except in a thin layer at the boundary; this is directly measurable by particle image velocimetry.
- The expansion speed of an indefinitely growing colony is constant and independent of colony size; it is linear in the growth force for non-motile colonies and quadratic in the motility speed for motile ones.
- Motility generates its own boundary tension $T$ that is quadratic in $v_0$ and linear in $G$, and it widens the range of negative homeostatic pressures for which a colony can grow indefinitely under tension.
- In two dimensions, growth is perimeter-controlled: cell number grows as $t^2$, radius grows linearly, and fingers at the frontier do not lead to fractal boundaries.
- The analytical profiles allow experimental stress profiles to be decomposed into friction and motility contributions, providing a procedure for inferring the underlying mechanics from traction measurements.
Reading between the lines
- The same mechanism should apply to other proliferating tissues or colonies on frictional substrates, wherever growth is enhanced at free edges and friction resists inward flow; this is an extension the paper does not develop.
- Because the surface growth term $\Delta k\Delta x$ is partly calibrated from the expansion speed it later predicts, a full test of the theory would require measuring $\Delta k$ and $\Delta x$ independently; absent that, the quantitative match is less decisive.
- A simultaneous measurement of stress and velocity fields in an expanding monolayer would separate this growth-friction mechanism from motility-alignment mechanisms: the former predicts inward flow correlated with tension, while alignment models predict outward-oriented traction.
- The predicted exponential buildup length $\lambda$ could be extracted experimentally from stress and velocity profiles and compared with the bulk response coefficient $\kappa$ measured by independent deformation experiments.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript combines 2PG-type agent-based simulations and a one-dimensional continuum theory to argue that cell colonies expanding on a substrate can be under tensile stress even in the absence of motility forces. The proposed mechanism is that cells near the free boundary grow faster than the bulk, producing an inward cell flux; substrate friction opposes this flux and generates tension when the homeostatic pressure is negative. The authors derive exponential pressure profiles, a constant expansion speed v(x0) = PH sqrt(kappa/(4 rho gamma)) + Delta k Delta x, and a phase boundary between finite and indefinitely growing colonies. They extend the model to motile colonies, where boundary polarization adds a motility tension term, and compare the resulting profiles and speeds with simulations.
Significance. If the central claim holds, the paper offers a mechanistically distinct route to tensile stress in freely expanding monolayers, complementing the established motility-alignment and confinement-sorting explanations. The analytical model is transparent, the simulation method is standard, and the predicted retrograde flow in non-motile tensile colonies is a concrete, falsifiable experimental target. The comparison between analytic pressure profiles and simulations in Fig. 3 is largely convincing for the shape of the profiles. However, the quantitative predictive claim for the expansion speed is weakened by the in-sample calibration of the surface-growth parameter Delta k Delta x from the very speed the theory is said to predict, as discussed below. The existence of tension in non-motile colonies is a direct simulation observation and is not undermined by this calibration issue.
major comments (3)
- [Section I, after Eq. (11) and SI Fig. S7] The central quantitative claim v(x0) = PH sqrt(kappa/(4 rho gamma)) + Delta k Delta x is presented as a prediction, but the text states 'We use this constant expansion speed to obtain Delta k Delta x.' SI Fig. S7(b) confirms that the surface-growth contribution is obtained as the difference between the total measured front speed and the independently computed bulk term. Therefore the agreement between Eq. (11) and the simulated front speed used for calibration is a consistency check, not an independent prediction. The linearity in G of the simulated speed is also not an independent test, because Delta k Delta x is derived from that speed and could absorb the same linear dependence. Please either measure Delta k locally in the simulations, provide a microscopic derivation of its value, or explicitly re-frame Eq. (11) as a parameterized relation rather than a parameter-free prediction.
- [Main text Fig. 3 caption vs. Section I, after Eq. (11)] The Fig. 3 caption claims that the theory curves are 'calculated from the piecewise solution with all parameters measured independently (i.e. not fitted)' and that only Delta x is estimated. This conflicts with the main-text statement that Delta k Delta x is obtained from the constant expansion speed. Since Delta k enters the piecewise pressure profile in Eq. (S6), the pressure and velocity comparisons in Fig. 3 are partially constrained by the fitted quantity, not fully independent. The manuscript should clearly list, for every parameter entering Eqs. (S6)-(S7), whether it is measured in a separate simulation, estimated from a structural assumption, or calibrated from the expansion speed.
- [Section I, Eq. (1) and phase-II mechanism] The existence of indefinite tensile growth for PH < 0 rests entirely on the surface-growth term Delta k Delta x in Eq. (1): without it, Eq. (11) gives a negative expansion speed and no phase II. The physical motivation for boundary-enhanced growth is plausible, but the manuscript does not provide a direct measurement of the local growth rate k(x) in the simulation rim, nor a derivation of Delta k from the 2PG microscopic rules. Given that this term is load-bearing, please report the measured k(x) profile in the boundary layer and compare it with the assumed step/delta form, or state clearly that the surface growth amplitude is an inferred effective parameter.
minor comments (4)
- [Throughout] The manuscript contains many typographical errors, including 'optianed', 'boudnary', 'therory', 'shortcommings', 'enabeled', 'reaveals', 'arrise', 'mechanims', 'accrding', 'consitituting', and 'partcile'. A thorough proofreading pass is needed.
- [SI Section S2] The text refers to 'Main text Fig.7(a)' when comparing quasi-1d and 2d expansion speeds, but the main text has Fig. 6(a) for this comparison; please correct the cross-reference.
- [SI Fig. S10 caption] The caption uses 'B^* = 11.4' where the context indicates this should be 'G^* = 11.4'.
- [Main text Fig. 3(b) and Section I, paragraph after Eq. (11)] The description of the pressure peak for expanding fronts with PH < 0 is qualitatively clear, but the figure inset is too small to see the predicted boundary pressure peak; a larger inset or an explicit zoom would help the reader verify this nontrivial prediction.
Circularity Check
Expansion-speed prediction is calibrated in-sample: Eq. (11)'s ΔkΔx is obtained from the measured front speed it then claims to reproduce; the motile Eq. (14) uses the same subtraction.
-
fitted input called prediction
[Section I, after Eq. (11)]
"Eq.(11) predicts a constant expansion speed for all colonies in phase II or III and is linear in PH− PC H . Indeed, the simulations display a constant expansion speed linear in G (see Fig.S7). We use this constant expansion speed to obtain ∆k∆x. Obtaining the homeostatic pressure and other bulk tissue properties from bulk simulations (i.e. without a fit, see SI), and estimating ∆x∗ = 0.7 reproduces the simulation data remarkably well (Fig.3)."
The claimed prediction v(x0)=PH sqrt(κ/(4ργ))+ΔkΔx is not out-of-sample: ΔkΔx is obtained from the very measured expansion speed by subtracting the independently computed bulk term PH sqrt(κ/(4ργ)) (SI Fig. S7 shows the surface contribution as the difference between total speed and bulk contribution). The subsequent agreement with the expansion speeds used for the subtraction is therefore a consistency check. Because ΔkΔx also enters the piecewise pressure and velocity profiles (Eq. S6), the profile comparisons in Fig. 3 are partly constrained by this calibration, although PH, κ, ρ, γ, and Δx are obtained separately.
-
fitted input called prediction
[SI Fig. S10 caption]
"The surface growth contribution (green) is obtained from the difference between the total expansion speed (blue) and the bulk contribution (orange)."
For motile colonies, the same surface-growth parameter ΔkΔx is calibrated from the measured total expansion speed and then inserted into Eq. (14), v(x0) = (PH + T/(1+λm/λ)) sqrt(κ/(4ργ)) + ΔkΔx. The agreement between Eq. (14) and the expansion-speed data is therefore not an independent test of the theory; it is guaranteed by construction for the surface-growth part. The motility contributions T and λm are independently measured from polarization and traction profiles, so the motile tension mechanism retains independent content, but the quantitative speed prediction is in-sample.
full rationale
The paper's central existence claim—that a non-motile expanding colony can be under tension—is a direct simulation observation and does not reduce to a fit; the analytic mechanism (friction balancing the inward flux produced by boundary-enhanced growth) is a coherent derivation from stated mechanical assumptions. Homeostatic pressure PH, friction γ, compressibility κ, and density ρ are measured independently, and Δx is estimated rather than fitted. However, the load-bearing quantitative prediction of the expansion speed is circular in an explicit, quotable way: Eq. (11) contains the parameter ΔkΔx, and the paper states 'We use this constant expansion speed to obtain ΔkΔx.' SI Fig. S7 confirms that ΔkΔx is the difference between the measured front speed and the computed bulk term, so the agreement in Fig. 3 and Fig. S7 is a consistency check rather than a prediction. The same calibration is reused in the motile expansion speed Eq. (14) via the caption of Fig. S10. The phase boundary between finite and infinite tensile colonies also inherits this calibration through PC_H = −ΔkΔx sqrt(4ργ/κ), so the predicted transition is not fully out-of-sample. An independent measurement of the surface-growth contribution (e.g., from local division-rate profiles in the rim) would settle the quantitative claim; without it, the central speed prediction is forced by construction. Self-citations to Refs. [2] and [25] motivate the surface-growth ansatz but are not the main circular step; the circularity is the in-sample calibration of ΔkΔx. Overall score 6: partial circularity of the central quantitative prediction, while the qualitative tensile mechanism remains independently supported by simulation.
Assumptions & free parameters
free parameters (3)
- ΔkΔx (surface growth amplitude) =
calibrated from expansion speed
- Δx* (surface growth layer width) =
0.7 in units of σ
- λm (motility force decay length) =
estimated from polarization/force density profile (Fig. S9)
assumptions (6)
- domain assumption Growth rate is linear in pressure: k = κ(PH - P).
- domain assumption The colony interior has constant density ρ.
- domain assumption In the quasi-1D colony the y-direction is homeostatically balanced: σyy = -PH.
- domain assumption Substrate friction is linear: f_ext = -2ργv.
- ad hoc to paper Surface growth is confined to a layer of width Δx, or to a delta function.
- ad hoc to paper The average cell polarization decays exponentially from the boundary.
Cite this review
Pith. "Pith review of Free Growth under Tension." pith.science (2026). https://pith.science/paper/WVORXACY
@misc{pith2026250203927,
author = {Pith},
title = {Pith review of: Free Growth under Tension},
year = {2026},
howpublished = {\url{https://pith.science/paper/WVORXACY}},
note = {Machine review of arXiv:2502.03927}
}
read the original abstract
Ever since the ground breaking work of Trepat et al. in 2009, we know that cell colonies growing on a substrate can be under tensile mechanical stress. The origin of tension has so far been attributed to cellular motility forces being oriented outward of the colony. Works in the field mainly revolve around how this orientation of the forces can be explained, ranging from velocity alignment, self-sorting due to self-propulsion, to kenotaxis. In this work, we demonstrate that tension in growing colonies can also be explained without cellular motility forces! Using a combination of well established tissue growth simulation technique and analytical modelling, we show how tension can arise as a consequence of simple mechanics of growing tissues. Combining these models with a minimalistic motility model shows how colonies can expand while under even larger tension. Furthermore, our results and analytical models provide novel analysis procedures to identify the underlying mechanics.
Figures
Reference graph
Works this paper leans on
-
[1]
Montel, F. et al. Stress Clamp Experiments on Multi- cellular Tumor Spheroids. Phys. Rev. Lett. 107, 188102 (2011)
work page 2011
-
[2]
Delarue, M. et al. Mechanical Control of Cell flow in Multicellular Spheroids. Phys. Rev. Lett. 110, 138103 (2013)
work page 2013
-
[3]
Hallatschek, O. et al. Proliferating active matter. Nat. Rev. Phys. 1–13 (2023)
work page 2023
-
[4]
Shraiman, B. I. Mechanical feedback as a possible regu- lator of tissue growth. Proc. Natl. Acad. Sci. 102, 3318– 3323 (2005)
work page 2005
-
[5]
Belteton, S. A. et al. Real-time conversion of tissue-scale mechanical forces into an interdigitated growth pattern. Nat. Plants 7, 826–841 (2021)
work page 2021
-
[6]
Basan, M., Risler, T., Joanny, J.-F., Sastre-Garau, X. & Prost, J. Homeostatic competition drives tumor growth and metastasis nucleation. HFSP J. 3, 265–272 (2009)
work page 2009
-
[7]
Trepat, X. et al. Physical forces during collective cell migration. Nat. Phys. 5, 426–430 (2009)
work page 2009
-
[8]
Basan, M., Prost, J., Joanny, J.-F. & Elgeti, J. Dissipa- tive particle dynamics simulations for biological tissues: Rheology and competition. Phys. Biol. 8, 026014 (2011)
work page 2011
Show all 44 references
-
[9]
& Levine, H
Basan, M., Elgeti, J., Hannezo, E., Rappel, W.-J. & Levine, H. Alignment of cellular motility forces with tissue flow as a mechanism for efficient wound healing. Proc. Natl. Acad. Sci. 110, 2452–2459 (2013)
2013
-
[10]
E., Hannezo, E., Trepat, X., Fredberg, J
Angelini, T. E., Hannezo, E., Trepat, X., Fredberg, J. J. & Weitz, D. A. Cell Migration Driven by Cooperative Substrate Deformation Patterns. Phys. Rev. Lett. 104, 168104 (2010)
2010
-
[11]
Angelini, T. E. et al. Glass-like dynamics of collective cell migration. Proc. Natl. Acad. Sci. 108, 4714–4719 (2011)
2011
-
[12]
Poujade, M. et al. Collective migration of an epithelial monolayer in response to a model wound. Proc. Natl. Acad. Sci. 104, 15988–15993 (2007)
2007
-
[13]
Ravasio, A. et al. Gap geometry dictates epithelial clo- sure efficiency. Nat. Commun. 6, 7683 (2015)
2015
-
[14]
Reffay, M. et al. Interplay of RhoA and mechanical forces in collective cell migration driven by leader cells. Nat. Cell Biol. 16, 217–223 (2014)
2014
-
[15]
M., Gelfand, I
Omelchenko, T., Vasiliev, J. M., Gelfand, I. M., Feder, H. H. & Bonder, E. M. Rho-dependent formation of ep- ithelial “leader” cells during wound healing. Proc. Natl. Acad. Sci. 100, 10788–10793 (2003)
2003
-
[16]
& Elgeti, J
Sarkar, D., Gompper, G. & Elgeti, J. A minimal model for structure, dynamics, and tension of monolayered cell colonies. Commun. Phys. 4, 1–8 (2021)
2021
-
[17]
Elgeti, J., Winkler, R. G. & Gompper, G. Physics of microswimmers—single particle motion and collective behavior: A review. Rep. Prog. Phys. 78, 056601 (2015)
2015
-
[18]
& Marchetti, M
Fily, Y. & Marchetti, M. C. Athermal Phase Separation of Self-Propelled Particles with No Alignment.Phys. Rev. Lett. 108, 235702 (2012)
2012
-
[19]
S., Hagan, M
Redner, G. S., Hagan, M. F. & Baskaran, A. Struc- ture and Dynamics of a Phase-Separating Active Col- loidal Fluid. Phys. Rev. Lett. 110, 055701 (2013)
2013
-
[20]
& Gompper, G
Elgeti, J. & Gompper, G. Wall accumulation of self- propelled spheres. EPL 101, 48003 (2013)
2013
-
[21]
Fily, Y., Baskaran, A. & F. Hagan, M. Dynamics of self- propelled particles under strong confinement. Soft Matter 10, 5609–5617 (2014)
2014
-
[22]
L., Garc´ ıa-Asenjo, J
Br´ u, A., Albertos, S., Subiza, J. L., Garc´ ıa-Asenjo, J. L. & Br´ u, I. The Universal Dynamics of Tumor Growth. Biophysical Journal 85, 2948–2961 (2003)
2003
-
[23]
Heinrich, M. A. et al. Size-dependent patterns of cell proliferation and migration in freely-expanding epithelia. eLife 9, e58945 (2020)
2020
-
[24]
Vedula, S. R. K. et al. Emerging modes of collective cell migration induced by geometrical constraints.Proc. Natl. Acad. Sci. 109, 12974–12979 (2012)
2012
-
[25]
& Elgeti, J
Podewitz, N., Delarue, M. & Elgeti, J. Tissue homeosta- sis: A tensile state. EPL 109, 58005 (2015)
2015
-
[26]
& Elgeti, J
Podewitz, N., J¨ ulicher, F., Gompper, G. & Elgeti, J. Interface dynamics of competing tissues. New J. Phys. 18, 083020 (2016)
2016
-
[27]
Tambe, D. T. et al. Collective cell guidance by coopera- tive intercellular forces. Nat. Mater. 10, 469–475 (2011)
2011
-
[28]
Sunyer, R. et al. Collective cell durotaxis emerges from long-range intercellular force transmission. Science 353, 1157–1161 (2016)
2016
-
[29]
Kim, J. H. et al. Propulsion and navigation within the advancing monolayer sheet. Nat. Mater. 12, 856–863 (2013)
2013
-
[30]
Marel, A.-K., Podewitz, N., Zorn, M., R¨ adler, J. O. & Elgeti, J. Alignment of cell division axes in directed ep- ithelial cell migration. New J. Phys. 16, 115005 (2014)
2014
-
[31]
Zisis, T. et al. Disentangling cadherin-mediated cell-cell interactions in collective cancer cell migration. Biophys- ical Journal 121, 44–60 (2022)
2022
-
[32]
& Trepat, X
Alert, R. & Trepat, X. Physical Models of Collective Cell Migration. Annu. Rev. Condens. Matter Phys. 11, 77–101 (2020)
2020
-
[33]
Vercruysse, E. et al. Geometry-driven migration effi- ciency of autonomous epithelial cell clusters. Nat. Phys. 20, 1492–1500 (2024)
2024
-
[34]
Williamson, John. J. & Salbreux, G. Stability and Roughness of Interfaces in Mechanically Regulated Tis- sues. Phys. Rev. Lett. 121, 238102 (2018). 8
2018
-
[35]
L., Gompper, G
B¨ uscher, T., Diez, A. L., Gompper, G. & Elgeti, J. Insta- bility and fingering of interfaces in growing tissue. New J. Phys. 22, 083005 (2020)
2020
-
[36]
Vedula, S. R. K. et al. Mechanics of epithelial closure over non-adherent environments. Nat. Commun. 6, 6111 (2015)
2015
-
[37]
& Vattulainen, I
Nikunen, P., Karttunen, M. & Vattulainen, I. How would you integrate the equations of motion in dissipative par- ticle dynamics simulations? Comput. Phys. Commun. 153, 407–423 (2003). Free Growth under Tension Chenyun Yao 1 and Jens Elgeti 1, ∗ 1Theoretical Physics of Living M...
2003
-
[38]
x∈ [0,x 0− ∆x] P =PH +D1ex/λ +D2e−x/λ
-
[39]
The boundary conditions read:
x∈ (x0− ∆x,x 0] P =PH + 2∆k κ +C1ex/λ +C2e−x/λ As in the main text, we assume mirror symmetry around x = 0 and thus only treat the positive half space here. The boundary conditions read:
-
[40]
Zero pressure at the boundary P (x0) = 0
-
[41]
Pressure and its derivative are continuous, in particular at x0− ∆x
-
[42]
Mirror symmetry at x = 0 implies dP dx = 0 at x = 0 Using the boundary conditions we fix the constants of integration: P = { PH(1−e x−x0 λ ) + 2∆k κ (cosh ∆x λ − 1)e x−x0 λ for x∈ [0,x 0− ∆x] (PH + 2∆k κ )(1−e x−x0 λ ) + 2∆k κ e− ∆x λ sinh (x−x0 λ ) for x∈ (x0− ∆x,x 0] (S6) ∗ ...
2025 arXiv
-
[43]
Local stress and pressure in an inhomogeneous system of spherical active Brownian particles
Shibananda Das, Gerhard Gompper, and Roland G Winkler. Local stress and pressure in an inhomogeneous system of spherical active Brownian particles. Scientific reports, 9(1):1–11, 2019. Publisher: Nature Publishing Group
2019
-
[44]
Multidimensional binary search trees used for associative searching
Jon Louis Bentley. Multidimensional binary search trees used for associative searching. Commun. ACM , 18(9):509–517, September 1975. 6 FIG. S6. Comparison of pressure and velocity profiles between different colonies. (a) The pressure profiles of non-motile colonies with differ...
1975
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.