Divergence of detachment forces in the finite Voronoi model
Pith reviewed 2026-05-10 09:24 UTC · model grok-4.3
The pith
The finite Voronoi model develops diverging detachment forces that make tissue fracture times depend on the simulation time step.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
When the line tension on cell-medium interfaces exceeds the tension on cell-cell contacts, the finite Voronoi model exhibits a divergence of detachment forces. This produces strong time-step dependence in the fracture timescale of initially intact active clusters, with decreasing time step unphysically suppressing rupture events. The divergence arises from the model's use of straight Voronoi edges and fixed-radius circular arcs for boundaries. A simple regularization removes the artifact, after which the model is calibrated to a deformable polygon representation to study how physical parameters control fracture under two strategies.
What carries the argument
The geometric construction of cell boundaries from straight Voronoi edges and fixed-radius circular arcs, which produces diverging detachment forces as separation approaches the arc radius.
Load-bearing premise
The fixed-radius circular arcs and straight Voronoi edges remain an adequate geometric representation of cell boundaries all the way to the instant of detachment without additional physical regularization or cutoff.
What would settle it
Run simulations of active cell clusters with cell-medium tension larger than cell-cell tension and check whether fracture time increases without bound as the time step is made smaller; if rupture is increasingly suppressed, the divergence is confirmed.
Figures
read the original abstract
Detachment and fracture are central to many tissue-level processes, but they are challenging to simulate with Voronoi-type models that typically assume a confluent tissue. Here we analyze the finite Voronoi model, a nonconfluent extension of conventional Voronoi models, in which cell boundaries are composed of straight Voronoi edges and circular arcs of fixed radius $\ell$. When the line tension on cell-medium interfaces exceeds the tension on cell-cell contacts, we find that the model exhibits a strong time-step dependence in the fracture timescale of initially intact active clusters: decreasing $\Delta t$ can unphysically suppress cluster rupture events. We trace this behavior to a divergence of detachment forces in the finite Voronoi model and introduce a simple regularization. Finally, we calibrate the near-detachment mechanics against a deformable polygon model and examine how key physical parameters control the tissue fracture timescale under two different calibration strategies. Our results show that, for studies focused on fracture or intercellular adhesion in nonconfluent monolayers, a physically motivated calibration of near-detachment mechanics in the finite Voronoi model is essential.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes the finite Voronoi model for nonconfluent cell monolayers, in which boundaries consist of straight Voronoi edges and fixed-radius circular arcs of length ℓ. It reports that when cell-medium line tension exceeds cell-cell tension, detachment forces diverge, producing strong unphysical dependence of fracture timescale on integration time step Δt that can suppress rupture events. The authors trace the divergence to the model geometry, introduce a regularization, calibrate near-detachment mechanics against a deformable-polygon reference model, and examine how key parameters control fracture timescale under two calibration strategies.
Significance. If the central results hold, the work is significant for the cell-mechanics simulation community because it identifies and mitigates a numerical artifact that directly affects the reliability of finite Voronoi models when used for detachment and fracture studies. Explicit demonstration of the Δt-dependent suppression, together with the regularization and the calibration to an independent polygon model, supplies concrete guidance for model users. The authors are credited for keeping the argument self-contained, for grounding the fix in an external benchmark, and for exploring parameter dependence under multiple calibration choices.
minor comments (3)
- The abstract refers to 'two different calibration strategies' without naming them; the main text should introduce and label these strategies (e.g., in §3 or §4) so that readers can follow the subsequent parameter sweeps without backtracking.
- Notation for the fixed arc radius ℓ and the two line tensions should be defined once in the introduction and then used consistently; a short table of symbols with units would improve clarity.
- Figure captions for the time-step dependence plots should explicitly list the Δt values employed and state whether the regularization is active, to allow immediate visual assessment of the artifact.
Simulated Author's Rebuttal
We thank the referee for their positive assessment of our work on the finite Voronoi model and for recommending minor revision. The review accurately summarizes the key findings regarding detachment force divergence and the proposed regularization and calibration approaches.
Circularity Check
No significant circularity; derivation is self-contained from model geometry
full rationale
The paper derives the detachment force divergence mathematically from the finite Voronoi construction (straight edges + fixed-radius arcs) when cell-medium tension exceeds cell-cell tension; this is an explicit geometric consequence shown in the equations, not a fit or redefinition. Regularization is introduced as a fix, and near-detachment calibration uses an independent deformable-polygon reference model as benchmark. No self-citations are load-bearing for the central claim, no ansatz is smuggled, and no prediction reduces to its own input by construction. The argument about time-step dependence and calibration necessity stands on its own internal consistency.
Axiom & Free-Parameter Ledger
free parameters (2)
- arc radius ℓ
- line tensions (cell-cell vs cell-medium)
axioms (2)
- domain assumption Cell boundaries can be represented by straight Voronoi edges joined to circular arcs of constant radius without loss of essential mechanics near detachment.
- domain assumption Force balance on cell vertices remains valid even as the contact angle approaches the singular configuration.
Reference graph
Works this paper leans on
-
[1]
Details of finite Voronoi model force calculations In the standard/finite Voronoi model, as given in Refs. [5, 11], an inner vertex connecting three cells𝑖,𝑗, and𝑘is given by hin =𝛼 𝑖r𝑖 +𝛼 𝑗r 𝑗 +𝛼 𝑘r𝑘 .(A1) The three barycentric coordinates of the circumcenter are 𝛼𝑖 =|r 𝑗 −r 𝑘 |2 (r𝑖 −r 𝑗 ) · (r𝑖 −r 𝑘)/𝐷,(A2a) 𝛼 𝑗 =|r 𝑖 −r 𝑘 |2 (r 𝑗 −r 𝑖) · (r 𝑗 −r 𝑘)/𝐷,...
-
[2]
Details of deformable polygon model In our simulations of the deformable polygon model, each cell is represented by𝑀=100 vertices. Enumerating the vertices{h 𝑚}of cell𝑖in clockwise order, the force exerted on vertex𝑚is given by f𝑚 =− 𝜕𝐸 𝜕h𝑚 .(A13) We need the derivative of area𝐴 𝑖 with respect toh 𝑚 [11]: 𝜕 𝐴𝑖 𝜕h𝑚 = (h𝑚−1 −h 𝑚+1) ׈𝑧 2 ,(A14) and derivati...
-
[3]
K. K. Youssef and M. A. Nieto, Epithelial–mesenchymal transi- tion in tissue repair and degeneration, Nature Reviews Molecular Cell Biology25, 720 (2024)
work page 2024
-
[4]
S. Tripathi, H. Levine, and M. K. Jolly, The physics of cellu- lar decision making during epithelial–mesenchymal transition, Annual Review of Biophysics49, 1 (2020)
work page 2020
-
[5]
K. L. Harper, M. S. Sosa, D. Entenberg, H. Hosseini, J. F. Cheung, R. Nobre, A. Avivar-Valderas, C. Nagi, N. Girnius, R. J. Davis,et al., Mechanism of early dissemination and metastasis in her2+ mammary cancer, Nature540, 588 (2016)
work page 2016
-
[6]
K. J. Cheung and A. J. Ewald, A collective route to metastasis: Seeding by tumor cell clusters, Science352, 167 (2016)
work page 2016
-
[7]
D. Bi, X. Yang, M. C. Marchetti, and M. L. Manning, Motility- driven glass and jamming transitions in biological tissues, Phys. Rev. X6, 021011 (2016)
work page 2016
-
[8]
D. Bi, J. Lopez, J. M. Schwarz, and M. L. Manning, A density- independent rigidity transition in biological tissues, Nature Physics11, 1074 (2015)
work page 2015
-
[9]
Y. Chen, Q. Gao, J. Li, F. Mao, R. Tang, and H. Jiang, Activation of topological defects induces a brittle-to-ductile transition in epithelial monolayers, Phys. Rev. Lett.128, 018101 (2022)
work page 2022
-
[10]
F. Graner and Y. Sawada, Can surface adhesion drive cell rear- rangement? Part II: a geometrical model, Journal of theoretical biology164, 477 (1993)
work page 1993
-
[11]
M. Bock, A. K. Tyagi, J.-U. Kreft, and W. Alt, Generalized Voronoi tessellation as a model of two-dimensional cell tissue dynamics, Bulletin of mathematical biology72, 1696 (2010). 11
work page 2010
-
[12]
G. Schaller and M. Meyer-Hermann, Multicellular tumor spheroid in an off-lattice Voronoi-Delaunay cell model, Phys. Rev. E71, 051910 (2005)
work page 2005
- [13]
- [14]
-
[15]
Nonomura, Study on multicellular systems using a phase field model, PloS one7, e33501 (2012)
M. Nonomura, Study on multicellular systems using a phase field model, PloS one7, e33501 (2012)
work page 2012
-
[16]
B. Palmieri, Y. Bresler, D. Wirtz, and M. Grant, Multiple scale model for cell migration in monolayers: Elastic mismatch be- tween cells enhances motility, Scientific reports5, 11745 (2015)
work page 2015
- [17]
-
[18]
W. Wang, R. A. Law, E. Perez Ipi ˜na, K. Konstantopoulos, and B. A. Camley, Confinement, jamming, and adhesion in cancer cells dissociating from a collectively invading strand, PRX Life 3, 013012 (2025)
work page 2025
-
[19]
S. A. Sandersius and T. J. Newman, Modeling cell rheology with the subcellular element model, Physical biology5, 015002 (2008)
work page 2008
-
[20]
C. R. Sweet, S. Chatterjee, Z. Xu, K. Bisordi, E. D. Rosen, and M. Alber, Modelling platelet–blood flow interaction using the subcellular element langevin method, Journal of The Royal Society Interface8, 1760 (2011)
work page 2011
-
[21]
S. Sandersius, C. J. Weijer, and T. J. Newman, Emergent cell and tissue dynamics from subcellular modeling of active biome- chanical processes, Physical biology8, 045007 (2011)
work page 2011
-
[22]
A. Boromand, A. Signoriello, F. Ye, C. S. O’Hern, and M. D. Shattuck, Jamming of deformable polygons, Phys. Rev. Lett. 121, 248003 (2018)
work page 2018
- [23]
-
[24]
S. Weng, R. J. Huebner, and J. B. Wallingford, Convergent exten- sion requires adhesion-dependent biomechanical integration of cell crawling and junction contraction, Cell Reports39, 110666 (2022)
work page 2022
-
[25]
D. M. Sussman and M. Merkel, No unjamming transition in a Voronoi model of biological tissue, Soft matter14, 3397 (2018)
work page 2018
-
[26]
E. Lawson-Keister, T. Zhang, F. Nazari, F. Fagotto, and M. L. Manning, Differences in boundary behavior in the 3d vertex and Voronoi models, PLoS computational biology20, e1011724 (2024)
work page 2024
-
[27]
H. Yue, C. R. Packard, and D. M. Sussman, Scale-dependent sharpening of interfacial fluctuations in shape-based models of dense cellular sheets, Soft Matter20, 9444 (2024)
work page 2024
-
[28]
W. Wang and B. A. Camley, wwang721/pyafv, Zen- odo.18091659 (2026), [version to be pinned]
work page 2026
-
[29]
W. Wang and B. A. Camley, Controlling tissue size by active fracture, Phys. Rev. E113, 034405 (2026)
work page 2026
-
[30]
E. L. Kaplan and P. Meier, Nonparametric estimation from in- complete observations, Journal of the American Statistical As- sociation53, 457 (1958)
work page 1958
-
[31]
A. G. Fletcher, M. Osterfield, R. E. Baker, and S. Y. Shvartsman, Vertex models of epithelial morphogenesis, Biophysical journal 106, 2291 (2014)
work page 2014
-
[32]
P. Spahn and R. Reuter, A vertex model of Drosophila ventral furrow formation, PloS one8, e75051 (2013)
work page 2013
- [33]
-
[34]
J. Vangheel, H. Ramon, and B. Smeets, Rigidity transitions in a three-dimensional active foam model of cell monolayers with frictional contact interactions, Phys. Rev. Res.8, 013022 (2026)
work page 2026
-
[35]
Y.-S. Chu, S. Dufour, J. P. Thiery, E. Perez, and F. Pincet, Johnson-Kendall-Roberts theory applied to living cells, Phys. Rev. Lett.94, 028102 (2005)
work page 2005
-
[36]
S. W. Byers, C. L. Sommers, B. Hoxter, A. M. Mercurio, and A. Tozeren, Role of E-cadherin in the response of tumor cell aggregates to lymphatic, venous and arterial flow: measurement of cell-cell adhesion strength, Journal of cell science108, 2053 (1995)
work page 2053
- [37]
- [38]
-
[39]
P. E. Kloeden and E. Platen,Numerical Solution of Stochastic Differential Equations(Springer, 1992)
work page 1992
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.