pith. sign in

arxiv: 2606.18604 · v1 · pith:2EVHV5K5new · submitted 2026-06-17 · ⚛️ physics.bio-ph

A symmetric relaxation method for entire two-dimensional cellular networks and its implications

Kai Xu , Lifan Weng , Zihan Wang , Yuyang Lian , Bin Huang This is my paper

Pith reviewed 2026-06-26 18:54 UTC · model grok-4.3

classification ⚛️ physics.bio-ph
keywords symmetric relaxation2D cellular networksvon Neumann-Mullins lawT1 transitionsVoronoi networksAboav-Weaire lawLewis lawforce disequilibrium
0
0 comments X

The pith

Symmetric relaxation based on angle symmetries simulates entire 2D cellular networks, reproduces statistical laws, and stabilizes against T1 transitions.

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

This paper introduces a symmetric relaxation method for modeling the dynamics of both inner and marginal vertices in two-dimensional cellular networks. The method relies on central angle symmetry within cells and angle symmetry at vertices to determine how the network relaxes. When applied to irregular Voronoi networks generated via regular hexagon disordering, the simulations match the von Neumann-Mullins law for cell area change, with a geometric correction term enhancing the fit. The approach also recovers the Aboav-Weaire and Lewis laws, and indicates that symmetric relaxation prevents neighbor exchanges by lengthening short edges and increasing area differences between neighbors. This leads to the suggestion that T1 events require force imbalance to overcome the symmetry-driven stability.

Core claim

The symmetric relaxation method, which treats inner vertices via cell central angle symmetry and marginal vertices via vertex angle symmetry, enables full-network simulations on trimmed Voronoi initial conditions. These simulations confirm agreement with the von Neumann-Mullins law for both cell types, where a modified form with geometric correction improves accuracy, and reproduce the Aboav-Weaire law and Lewis law showing cells approach maximum inscribed polygons in ellipses. The method reveals that symmetric relaxation inhibits T1 transitions by reducing short edges and amplifying area disparity, implying T1 occurs when force disequilibrium surpasses this stabilization.

What carries the argument

The symmetric relaxation method for inner and marginal vertices, determined respectively by central angle symmetry of associated cells and angle symmetry at each vertex.

If this is right

  • Simulations of relaxed networks agree with the von Neumann-Mullins law for area evolution of both inner and marginal cells.
  • A modified von Neumann-Mullins equation including a geometric correction term significantly improves the quality of predictions.
  • The Aboav-Weaire law and Lewis law are reproduced in the relaxed networks.
  • Relaxed cells tend to approach the maximum inscribed polygons of ellipses according to the Lewis law.
  • Symmetric relaxation inhibits T1 neighbor exchange transitions by reducing the number of short edges and increasing area disparity among neighboring cells.

Where Pith is reading between the lines

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

  • The method could be tested on experimental images of real tissues to see if the same symmetry assumptions hold without additional forces.
  • Similar symmetry-based relaxation might apply to three-dimensional foams, though the geometry of vertices and cells would require adaptation.
  • The proposed trigger for T1 suggests that models without explicit force calculations may still capture transition statistics through effective rules derived from symmetry.
  • The conserved edge number distribution at irregularity value one indicates a possible universal feature of 2D networks that persists after relaxation.

Load-bearing premise

The assumption that central angle symmetry of associated cells and angle symmetry at each vertex are the dominant factors determining relaxation dynamics for both inner and marginal vertices, without other physical forces or asymmetries substantially altering the outcome.

What would settle it

If enforcing symmetric relaxation in simulations still allows frequent T1 events even at low force disequilibrium, or if the modified von Neumann-Mullins equation fails to improve predictions on independent network data, the central claims would be challenged.

read the original abstract

To simulate the relaxation of an entire 2D cellular network, this study proposes a symmetric relaxation method for both inner and marginal vertices. The relaxations of these two types of vertices are determined by the central angle symmetry of associated cells and the angle symmetry at each vertex, but with different major considerations. Trimmed Voronoi networks with varying irregularity are used as initial networks for the relaxation simulation. In particular, we propose a regular hexagon disordering method to generate Voronoi networks and find that the inner cells of networks with an irregularity value of one exhibit a conserved edge number distribution, as found in other 2D cellular networks. Simulation results agree with the von Neumann-Mullins law for both inner and marginal cells, and a modified equation including a geometric correction term significantly improves prediction quality. The Aboav-Weaire law and Lewis law are also reproduced, with the latter showing that relaxed cells tend to approach the ellipses' maximum inscribed polygons. Analysis of edge length, interior angle, and shape index reveals that symmetric relaxation inhibits T1 (neighbour exchange) topological transitions by reducing short edges while increasing area disparity among neighbouring cells. The findings suggest that T1 events may be triggered when force disequilibrium overcomes the stabilising effect of symmetric relaxation, providing a possible mechanistic explanation for T1 in 2D foams.

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 introduces a symmetric relaxation method for simulating the dynamics of entire 2D cellular networks, including marginal vertices. Inner vertices relax according to central-angle symmetry of adjacent cells, while marginal vertices use vertex-angle symmetry with different weighting. Initial networks are generated via a regular-hexagon disordering procedure applied to trimmed Voronoi tessellations at varying irregularity. Simulations are reported to recover the von Neumann-Mullins relation for both interior and boundary cells; a modified form that adds an empirical geometric correction term improves the fit. The Aboav-Weaire and Lewis laws are also recovered, and the authors interpret the suppression of T1 transitions under symmetric relaxation as evidence that topological changes occur when force imbalance overcomes the stabilizing symmetry.

Significance. If the symmetry rules can be shown to be equivalent to (or dominant over) mechanical force balance, the method would offer a computationally lightweight route to large-scale network evolution and a concrete mechanistic hypothesis for T1 onset in foams. The reproduction of three independent empirical laws on the same set of relaxed networks is a positive feature. However, the absence of a derivation linking the angle-symmetry update to dA/dt ∝ (n−6), the lack of any comparison against a standard tension-driven or curvature-driven solver, and the introduction of a free geometric correction term together limit the strength of the mechanistic claim.

major comments (3)
  1. [Results section on von Neumann-Mullins agreement] The central claim that the symmetric relaxation reproduces the von Neumann-Mullins law rests on the unshown step that the angle-symmetry update rule yields dA/dt = k(n−6). No derivation or perturbative test is supplied showing that the chosen update produces this proportionality rather than imposing it by construction.
  2. [Modified equation and associated figures/tables] The modified von Neumann-Mullins equation that includes a geometric correction term is reported to improve predictive quality, yet the term is listed among the free parameters and no first-principles derivation from the symmetry rules is given. This makes the reported improvement difficult to distinguish from a post-hoc fit.
  3. [Discussion of T1 mechanism] The mechanistic interpretation that T1 events are triggered when force disequilibrium overcomes symmetric-relaxation stabilization assumes that the angle-symmetry rule is equivalent to (or dominant over) the physical condition ΣT_i = 0. No perturbation test with unequal line tensions or explicit curvature terms, and no benchmark against a standard mechanical solver, is presented to support this dominance.
minor comments (2)
  1. [Methods] The regular-hexagon disordering procedure used to generate initial networks is described only in prose; a short pseudocode block or explicit formula for the irregularity parameter would improve reproducibility.
  2. [Method description] Notation for inner versus marginal vertices and the precise weighting between cell-angle and vertex-angle contributions is introduced without a compact summary table; readers must extract the rules from scattered paragraphs.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the detailed and constructive report. We address each major comment below, clarifying the scope of our claims while agreeing where revisions can strengthen the presentation.

read point-by-point responses

  1. Referee: [Results section on von Neumann-Mullins agreement] The central claim that the symmetric relaxation reproduces the von Neumann-Mullins law rests on the unshown step that the angle-symmetry update rule yields dA/dt = k(n−6). No derivation or perturbative test is supplied showing that the chosen update produces this proportionality rather than imposing it by construction.

    Authors: The symmetric relaxation method enforces local angle symmetries at vertices and is not constructed to impose the von Neumann-Mullins relation. The observed agreement with dA/dt ∝ (n−6) for both interior and marginal cells emerges from the collective dynamics in the simulations on disordered Voronoi networks. We do not claim a first-principles derivation from the update rule. We will add a clarifying paragraph in the revised manuscript stating that the reproduction is empirical and outlining why the symmetry rules are expected to produce area changes correlated with side number. revision: yes

  2. Referee: [Modified equation and associated figures/tables] The modified von Neumann-Mullins equation that includes a geometric correction term is reported to improve predictive quality, yet the term is listed among the free parameters and no first-principles derivation from the symmetry rules is given. This makes the reported improvement difficult to distinguish from a post-hoc fit.

    Authors: The geometric correction is introduced as an empirical adjustment to account for residual irregularities and boundary effects not fully captured by the basic symmetry rules. We agree that its status as a free parameter limits mechanistic interpretation. In revision we will relabel the term explicitly as empirical, move the associated fit statistics to a supplementary table, and discuss possible geometric origins tied to the angle-symmetry weighting without claiming a derivation. revision: partial

  3. Referee: [Discussion of T1 mechanism] The mechanistic interpretation that T1 events are triggered when force disequilibrium overcomes symmetric-relaxation stabilization assumes that the angle-symmetry rule is equivalent to (or dominant over) the physical condition ΣT_i = 0. No perturbation test with unequal line tensions or explicit curvature terms, and no benchmark against a standard mechanical solver, is presented to support this dominance.

    Authors: The manuscript presents the force-disequilibrium trigger as a possible mechanistic hypothesis suggested by the observation that symmetric relaxation strongly suppresses T1 transitions. It does not assert equivalence or dominance over the vector force-balance condition. We acknowledge that direct tests with unequal tensions or comparisons to curvature-driven or tension-driven solvers lie outside the present study. In revision we will expand the discussion section to state the assumptions and limitations of the interpretation more explicitly. revision: partial

Circularity Check

0 steps flagged

No significant circularity; derivation remains self-contained against external benchmarks

full rationale

The paper defines a symmetric relaxation update rule from central-angle and vertex-angle symmetry (distinct rules for inner vs. marginal vertices), applies it to trimmed Voronoi initial networks generated by a regular-hexagon disordering procedure, and then compares the resulting dynamics to the independently known von Neumann-Mullins, Aboav-Weaire and Lewis laws. No equation in the provided text shows the symmetry rule being defined in terms of dA/dt ∝ (n-6) or any other target law; the reported agreement is therefore an external check rather than a reduction by construction. The modified geometric-correction term is presented as an empirical improvement to an existing law, not as the input that was fitted. No self-citation chain is invoked to justify uniqueness or to forbid alternatives. The central claim therefore rests on independent simulation output compared with established empirical patterns.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

Only the abstract is available, so the ledger is inferred from stated elements; the geometric correction term is treated as a likely free parameter and the symmetry rules as a domain assumption.

free parameters (2)
  • geometric correction term
    Added to the von Neumann-Mullins equation to improve prediction quality; its value is not derived from first principles in the abstract.
  • irregularity value
    Parameter controlling disorder in the initial Voronoi networks; value of one is highlighted for conserved edge number distribution.
axioms (1)
  • domain assumption Relaxations of inner and marginal vertices are determined by central angle symmetry of associated cells and angle symmetry at each vertex.
    This is the foundational premise of the proposed symmetric relaxation method as described in the abstract.

pith-pipeline@v0.9.1-grok · 5773 in / 1465 out tokens · 36928 ms · 2026-06-26T18:54:28.490072+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

33 extracted references

  1. [1]

    The Perimeter-Area Relation

    Feder J. The Perimeter-Area Relation. In: Feder J, editor. Fractals. Boston, MA: Springer US; 1988. p. 200–11

    1988

  2. [2]

    Soap, cells and statistics —random patterns in two dimensions

    Weaire D, Riv ier N. Soap, cells and statistics —random patterns in two dimensions. Contemporary Physics. 1984;25(1):59–99

    1984

  3. [3]

    The kinetics of cellular patterns

    Glazier JA, Weaire D. The kinetics of cellular patterns. Journal of Physics: 11 Condensed Matter. 1992;4(8):1867–94

    1992

  4. [4]

    Equilibrium states and ground state of two-dimensional fluid foams

    Graner F, Jiang Y, Janiaud E, Flament C. Equilibrium states and ground state of two-dimensional fluid foams. Physical Review E. 2000;63:011402

    2000

  5. [5]

    Topological correlations in Bénard -Marangoni convective structures

    Cerisier P, Rahal S, Rivier N. Topological correlations in Bénard -Marangoni convective structures. Phys Rev E. 1996;54(5):5086–94

    1996

  6. [6]

    Boundary effects in biological planar networks: pentagons dominate Pyropia marginal cells

    Xu K, He F. Boundary effects in biological planar networks: pentagons dominate Pyropia marginal cells. Philosophical Magazine. 2026:1–13

    2026

  7. [7]

    Vertex models capturing subcellular scales in epithelial tissues

    Turlier H, Lange Z, Matthäus F, Qiu M. Vertex models capturing subcellular scales in epithelial tissues. PLOS Computational Biology. 2025;21(5):e1012993

    2025

  8. [8]

    Vertex stability and topological transitions in vertex models of foams and epithelia

    Spencer MA, Jabeen Z, Lubensky DK. Vertex stability and topological transitions in vertex models of foams and epithelia. The European Physical Journal E. 2017;40(1)

    2017

  9. [9]

    Vertex models of epithelial morphogenesis

    Fletcher AG, Osterfield M, Baker RE, Shvartsman SY. Vertex models of epithelial morphogenesis. Biophys J. 2014;106(11):2291–304

    2014

  10. [10]

    A curvature -driven vertex model for two-dimensional grain growth

    Marsh SP, Masumura RA, Pande CS. A curvature -driven vertex model for two-dimensional grain growth. Philosophical Magazine Letters. 1995;72(6):429–34

    1995

  11. [11]

    Scaling behavior of two -dimensional domain growth: Computer simulation of vertex models

    Nakashima K, Nag ai T, Kawasaki K. Scaling behavior of two -dimensional domain growth: Computer simulation of vertex models. Journal of Statistical Physics. 1989;57(3-4):759–87

    1989

  12. [12]

    A Nonlinear Vertex -based Model for Animation of Two-dimensional Dry Fo am

    Kelager M, Erleben K. A Nonlinear Vertex -based Model for Animation of Two-dimensional Dry Fo am. GRAPP 2010 - Proceedings of the International Conference on Computer Graphics Theory and Applications; Angers, France2010

    2010

  13. [13]

    Tissue fluidity promotes epithelial wound healing

    Tetley RJ, Staddon MF, Heller D, Hoppe A, Banerjee S, Mao Y. Tissue fluidity promotes epithelial wound healing. Nature Physics. 2019;15(11):1195–203

    2019

  14. [14]

    A geometry -based relaxation algorithm for equilibrating a trivalent polygonal network in two dimensions and its implications

    Xu K. A geometry -based relaxation algorithm for equilibrating a trivalent polygonal network in two dimensions and its implications. Philosophical Magazine. 2021;101(14):1632–53

    2021

  15. [15]

    Dynamic crushing of 2D cellular structures: A finite element study

    Zheng Z, Yu J, Li J. Dynamic crushing of 2D cellular structures: A finite element study. International Journal of Impact Engineering. 2005;32(1-4):650–64

    2005

  16. [16]

    Fitting quadratic curves to data points

    Chernov N, Huang Q, Ma H. Fitting quadratic curves to data points. British Journal of Mathematics & Computer Science. 2014;4(1):33–60

    2014

  17. [17]

    Building block analysis of 2D amorphous networks reveals medium range correlation

    Büchner C, Liu L, Stuckenholz S, Burson KM, Lichtenstein L, Heyde M, et al. Building block analysis of 2D amorphous networks reveals medium range correlation. Journal of Non-Crystalline Solids. 2016;435:40–7

    2016

  18. [18]

    Cell s tile a flat plane by controlling geometries during morphogenesis of Pyropia thalli

    Xu K, Xu Y, Ji D, Chen T, Chen C, Xie C. Cell s tile a flat plane by controlling geometries during morphogenesis of Pyropia thalli. PeerJ. 2017;5:e3314

    2017

  19. [19]

    The emergence of geometric order in proliferating metazoan epithelia

    Gibson MC, Patel AB, Nagpal R, Perrimon N. The emergence of geometric order in proliferating metazoan epithelia. Nature. 2006;442(7106):1038–41

    2006

  20. [20]

    Control of the Mitotic Cleavage Plane by Local Epithelial Topology

    Gibson WT, Veldhuis JH, Rubinstein B, Cartwright HN, Perrimon N, Brodland GW, et al. Control of the Mitotic Cleavage Plane by Local Epithelial Topology. Cell. 2011;144(3):427–38

    2011

  21. [21]

    EpiTools: 12 an open-source image analysis toolkit for quantifying epithelial growth dynamics

    Heller D, Hoppe A, Restrepo S, Gatti L, Tournier AL, Tapon N, et al. EpiTools: 12 an open-source image analysis toolkit for quantifying epithelial growth dynamics. Developmental cell. 2016;36(1):103–16

    2016

  22. [22]

    Hexagonal packing of Drosophila wing epithelial cells by the planar cell pol arity pathway

    Classen A -K, Anderson KI, Marois E, Eaton S. Hexagonal packing of Drosophila wing epithelial cells by the planar cell pol arity pathway. Developmental cell. 2005;9(6):805–17

    2005

  23. [23]

    Two‐dimensional motion of idealized grain boundaries

    Mullins WW. Two‐dimensional motion of idealized grain boundaries. Journal of Applied Physics. 1956;27(8):900–4

    1956

  24. [24]

    In: Herring C, editor

    Von Neumann J. In: Herring C, editor. Metal interfaces. Cleveland, Ohio: American Society for Metals; 1952. p. 108–10

    1952

  25. [25]

    Fluid Networks

    AREF H, HERDTLE T. Fluid Networks. In: MOFFATT HK, TSINOBER A, editors. Topological Fluid Mechanics: Cambridge University Press; 1990

    1990

  26. [26]

    Numerical studies of foam dynamics: University of Califo rnia, San Diego; 1991

    Herdtle T. Numerical studies of foam dynamics: University of Califo rnia, San Diego; 1991

    1991

  27. [27]

    Vertex models for two-dimensional grain growth

    Kawasaki K, Nagai T, Nakashima K. Vertex models for two-dimensional grain growth. Philosophical Magazine B. 1989;60(3):399–421

    1989

  28. [28]

    Ellipse packing in two -dimensional cell tessellation: a theoretical explanation for Lewis’s law and Aboav-Weaire’s law

    Xu K. Ellipse packing in two -dimensional cell tessellation: a theoretical explanation for Lewis’s law and Aboav-Weaire’s law. PeerJ. 2019;7:e6933

    2019

  29. [29]

    Geometric formulas of Lewis’s law and Aboav -Weaire’s law in two dimensions based on ellipse packing

    Xu K. Geometric formulas of Lewis’s law and Aboav -Weaire’s law in two dimensions based on ellipse packing. Philosophical Magazine Letters. 2019;99(9):317–25

    2019

  30. [30]

    The characteristics of maximum inscribed and minimum circumscribed polygons of ellipse (in Chinese)

    Su H. The characteristics of maximum inscribed and minimum circumscribed polygons of ellipse (in Chinese). Mathematics Teaching 1987;6:22–6

    1987

  31. [31]

    Structure-dependent mobility of a dry aqueous foam flowing along two parallel channels

    Jones SA, Dollet B, Méheust Y, Cox SJ, Cantat I. Structure-dependent mobility of a dry aqueous foam flowing along two parallel channels. Physics of Flu ids. 2013;25(6)

    2013

  32. [32]

    Viscous froth model applied to the motion and topological transformations of two-dimensional bubbles in a channel: three-bubble case

    Torres-Ulloa C, Grassia P. Viscous froth model applied to the motion and topological transformations of two-dimensional bubbles in a channel: three-bubble case. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences. 2022;478(2258)

    2022

  33. [33]

    Dynamics of two -dimensional soap froths

    Glazier JA, Gross SP, Stavans J. Dynamics of two -dimensional soap froths. Physical Review A. 1987;36(1):306–12. 13 Figure 1 Diagrams of symmetric relaxation for inner and marginal vertices. (A) Relaxation moves each inner vertex 𝑉 toward its target position 𝑃, which is the centroid of a triangle formed by three optimal rays corresponding to the vertex (1...

    1987