## Abstract

Mechanical equilibrium states of cellular matter are overwhelmingly metastable and separated from each other by topology changes. Using theory and simulations, it is shown that for a wide class of energy functionals in 2D, including those describing tissue cell layers, local energy differences between neighboring metastable states as well as global energy differences between initial states and ground states are governed by simple, universal relations. Knowledge of instantaneous length of an edge undergoing a T1 transition is sufficient to predict local energy changes, while the initial edge length distribution yields a successful prediction for the global energy difference. An analytical understanding of the model parameters is provided.

Original language | English (US) |
---|---|

Article number | 248001 |

Journal | Physical review letters |

Volume | 120 |

Issue number | 24 |

DOIs | |

State | Published - Jun 11 2018 |

## ASJC Scopus subject areas

- Physics and Astronomy(all)

## Fingerprint

Dive into the research topics of 'Universal Features of Metastable State Energies in Cellular Matter'. Together they form a unique fingerprint.## Cite this

- APA
- Standard
- Harvard
- Vancouver
- Author
- BIBTEX
- RIS

*Physical review letters*,

*120*(24), Article 248001. https://doi.org/10.1103/PhysRevLett.120.248001

**Universal Features of Metastable State Energies in Cellular Matter.**/ Kim, Sangwoo; Wang, Yiliang; Hilgenfeldt, Sascha.

In: Physical review letters, Vol. 120, No. 24, 248001, 11.06.2018.

Research output: Contribution to journal › Article › peer-review

*Physical review letters*, vol. 120, no. 24, 248001. https://doi.org/10.1103/PhysRevLett.120.248001

**Universal Features of Metastable State Energies in Cellular Matter**. In: Physical review letters. 2018 ; Vol. 120, No. 24.

}

TY - JOUR

T1 - Universal Features of Metastable State Energies in Cellular Matter

AU - Kim, Sangwoo

AU - Wang, Yiliang

AU - Hilgenfeldt, Sascha

N1 - Funding Information: We thank Ken Brakke for valuable advice on surface evolver code implementation, and the NSF under Grant No. 1504301 for support. Funding Information: Kim Sangwoo * Wang Yiliang * Hilgenfeldt Sascha Mechanical Sciences and Engineering, University of Illinois , Urbana-Champaign, Illinois 61801, USA * These authors contributed equally to this work. 11 June 2018 15 June 2018 120 24 248001 26 March 2018 31 January 2018 © 2018 American Physical Society 2018 American Physical Society Mechanical equilibrium states of cellular matter are overwhelmingly metastable and separated from each other by topology changes. Using theory and simulations, it is shown that for a wide class of energy functionals in 2D, including those describing tissue cell layers, local energy differences between neighboring metastable states as well as global energy differences between initial states and ground states are governed by simple, universal relations. Knowledge of instantaneous length of an edge undergoing a T 1 transition is sufficient to predict local energy changes, while the initial edge length distribution yields a successful prediction for the global energy difference. An analytical understanding of the model parameters is provided. National Science Foundation 10.13039/100000001 1504301 In interacting many-particle systems, energy landscapes are complex and hard to analyze, in particular when disorder precludes symmetries. Considerable effort has focused on particle aggregates with short-range interactions (hard core or soft) in the context of granular media [1–4] , optimal packings [5–8] , or the description of jammed states [9–12] . In cellular matter , on the other hand, the main energy contributions result from the shape and properties of the interfaces between deformable domains that fill available space (with a negligible continuous phase), making the interfaces surfaces of polygons (in two dimensions) or polyhedra (in three dimensions) [13,14] . The exclusion of bulk energy contributions generally means that the areas (2D) or volumes (3D) of individual domains remain constant, while their shape and relative positioning is variable; the simplest physical example is a dry soap froth [15] . Cellular matter also includes large classes of systems considered in the context of modeling biological tissues, with energy contributions from elasticity and cell-cell adhesion [16–18] , bulk elasticity [19,20] , or viscous effects [21,22] . Recent work has focused on low-energy states of cellular matter, which we will call “ground states,” although the global lowest-energy state is in general unknown and may not be unique [23] . Ground states in this sense are found through a variety of protocols and annealing strategies, and their energies are typically insensitive to the method [23,24] . Tissuelike 2D systems show a qualitative transition of the ground state: For low values of interdomain adhesion energy (relative to elastic deformation penalties), the material ground state retains rigidity (finite resistance to external forces) [18] , while for higher adhesion it becomes degenerate [25,26] with individual domains minimizing their energy separately (the material becomes “floppy”). This “loss of rigidity transition” [19] occurs for static as well as for fluctuating systems, where it resembles a solid-fluid transition [19,27] . The present work, by contrast, focuses on metastable states (local energy minima) significantly above the ground state energy. These are common in cellular systems in nature: If the energy barriers exceed thermal energies, as expected for domains above colloidal size, a system needs induced stimuli to evolve towards the ground state. We show how limited information on the geometry of a generic rigid 2D cellular system quantitatively predicts the energy of the metastable states, governs individual topological transitions, and describes an efficient pathway of lowering energy towards a ground state. Cellular matter domains (identified by index i ) interact with nearest neighbors only, each contributing to a total energy E = ∑ i E i . Restricting ourselves to 2D systems, the requirement of dominant interfacial energy means that we can generally write E i = ∫ P i u P , i ( s ) d s = c 0 i P i + 1 2 c 1 i P i 2 + ⋯ , (1) expanding the general energy per length u P , i ( s ) in the perimeter lengths P i . The shown truncation after the second term is representative of the generic class of 2D tissue model studied in the recent literature [16,18,19,25,28] . For domains or cells of the same type, the coefficients are uniform ( c 0 , c 1 ). Setting c 1 = 0 describes a 2D foam, identifying the interfacial tension with c 0 ( = 1 without loss of generality). Including the second order term in P i recovers the general case of “tissue” energy. By rearranging terms, we obtain the following two functionals, E f = 1 2 L 0 ∑ i P i , E t = 1 2 L 0 ∑ i ( ( P i - P i , 0 ) 2 P i , 0 - γ P i ) , (2) where all lengths are normalized by L 0 , the edge length of a regular hexagon of area A tot / N for a system of N domains covering an area A tot . In E t , the first and second terms can be interpreted as perimeter elasticity and adhesion energy, respectively [16,18] . P i , 0 is the mechanical equilibrium perimeter of cell i in isolation, here chosen as the perimeter of a circle with the same area as domain i (other choices of P i , 0 merely rescale relevant energy differences [18] ). The dimensionless adhesion strength γ is normalized by the perimeter elastic modulus. It was shown [17–19,25] that loss of rigidity occurs when γ > γ c ≈ 0.12 . Below this value, γ is a nongeometric determinant of system energy. We will show that energies can, nevertheless, be inferred from geometry alone. Metastable states are separated from each other by T 1 topological transitions [29] , where a single edge of length L reorients (it “flips”) to change topology of four neighboring cells [Fig. 1(a) ]. For rigid or solid states (foams or tissues with γ < γ c ), there are metastable states on both sides of the transition, while the intermediate configuration of four-way-connected edges is a local maximum (the energy barrier) [30] . 1 10.1103/PhysRevLett.120.248001.f1 FIG. 1. (a) A T 1 transition flipping an edge of initial length L , going through an unstable four-way intermediate configuration and ending up with altered topologies of adjacent cells. Colors represent polygonal edge number: 5 (yellow), 6 (gray), 7 (blue). (b) Binned plot of edge length change | δ L | as a function of distance r from the center of the T 1 edge; r s (vertical dashed line) is the distance beyond which the analytical model assumes stochastic length changes. (c) The standard deviations of absolute (open symbols) and relative (filled symbols) edge length changes for r > r s decay as exp ( - κ r ) (black dashed line). We evaluate metastable states in surface evolver (SE) [31] with the quadratic or circular arc vertex models (i.e., edges between domains contain additional vertices and are generally not straight), on rectangles with periodic boundary conditions containing typically N = 400 or 900 domains. Initial patterns are Voronoi constructs from various point distributions. SE fixes domain areas to match a desired area distribution (its coefficient of variation c A quantifies polydispersity) and finds a local energy minimum with the given topology and energy functional. We analyze the metastable states after and before the T 1 . Geometrically, T 1 s are local events—the edge length changes | δ L | and their standard deviation σ δ L decay exponentially with distance from the flipping edge. Figures 1(b) and 1(c) identify a characteristic decay scale κ - 1 ≈ 2.8 , in quantitative agreement with earlier findings [32] . This study focuses not on the energy barrier height between the states [28] , but on the distribution of the energy differences Δ E between the system energies after and before the T 1 . We find that the expectation value of Δ E has a strong linear correlation with the initial length L of the flipping edge: Δ E = α ( L - L c ) . (3) As seen in Fig. 2(a) , the scatter around this linear relation is particularly small for energetically favorable T 1 s ( Δ E < 0 ). For this range, data from ≳ 30 000 T 1 transitions for systems of various polydispersities and various energy functionals were analyzed. It is a surprising fact that α and L c are found to be system independent and robust against protocol changes: (i) different methods of domain preparation—see Supplemental Material [33] —have no perceptible effect on Eq. (3) or on the scatter of the data; (ii) the order of T 1 s is irrelevant; (iii) widely different polydispersities result in the same relation [Fig. 2(b) ]; (iv) even simulations using E t are in quantitative agreement with those using E f : If the perimeter lengths of the tissue in mechanical equilibrium are P i * , computing E f * = ( 1 / 2 L 0 ) ∑ i P i * yields an equivalent foam energy whose correlation with L is quantitatively the same [Fig. 2(b) ], even though the energetics of the T 1 processes that yield the configurations are quite different, and the configurations are not metastable states under E f . The linear correlation Eq. (3) remains unchanged for all tense tissues ( γ < γ c ). Beyond γ c , system states lose rigidity and all Δ E are trivially zero. Very recent work on three-dimensional epithelial sheets [27] likewise finds a linear relation between flipping-face area and energy differences in agreement with Eq. (3) . 2 10.1103/PhysRevLett.120.248001.f2 FIG. 2. (a) Energy change by T 1 transition Δ E versus initial edge length L in monodisperse foam samples (blue) with binned data for Δ E < 0 (orange circles) and the overall linear relation Eq. (3) (dashed line); (b) the same in polydisperse foams and polydisperse tissue systems with different adhesion; (c) simplified local configurations for analytical calculations: (left) monodisperse hexagonal pattern, (center) one quadruple defect with equal cell areas, (right) a quadruple defect with area polydispersity (at fixed total area). Vertex positions indicated in red are fixed, the others represent optimization variables. (d) Combined Δ E ( L ) data from all simulations, binned (circles) with best linear fit using α , L c (black) and analytical linear relation using α a , L c , a (red); (e) like (d), for Δ E ( ℓ ) data and correlations using β , ℓ c and their analytical analogs β a , ℓ c , a . An average over all data is described well by Eq. (3) with a universal critical edge length L c ≈ 0.611 and a universal slope α ≈ 0.827 . Beyond empirical data, we can obtain analytical approximations to α and L c from the simple elementary T 1 transition between a honeycomb pattern and a quadruple defect [two neighboring dislocations, Fig. 2(c) ]. Changing the areas of the pentagons and heptagons generates elementary polydisperse configurations. All vertices not belonging to cells participating in the T 1 are fixed at their honeycomb positions [red in Fig. 2(c) ]. Minimizing E f with respect to the remaining degrees of freedom yields analytical metastable state geometries as solutions to a system of algebraic equations (see Supplemental Material [33] ). A linear fit to the resulting Δ E ( L ) values obtains α a ≈ 0.791 and L c , a ≈ 0.627 , in very good agreement with data [Fig. 2(d) ]. Now we use this information about energetic effects of (spatially local) T 1 s to infer the global energy of a given metastable state, not only for purposes of easy general diagnostics, but in order to assess whether metastability interferes with the ability to detect the loss-of-rigidity transition mentioned above. Both simulations and analytical computations are used. The simulations should reflect processes of mechanical excitation overcoming energy barriers, e.g., by shearing foams [34,35] , agitating emulsions [36,37] , or by cell mobility in tissues [38] , so that the system energy approaches a ground state through successive T 1 s. Equation (3) suggests flipping short edges ( L < L c ) will selectively lower the system energy. However, simulations may miss energetically favorable edge flips if these edges are surrounded by large-area cells (they are relatively short, but absolutely longer than L c ). Therefore, we shall focus on relative edge length ℓ , ℓ = L / min ( L 0 i , L 0 j ) , (4) where L 0 i = 2 1 / 2 3 - 3 / 4 A i 1 / 2 is the edge length of a regular hexagon of area A i , and the domains i , j share the edge. As Fig. 2(e) shows, Δ E ( ℓ ) is still a linear function, Δ E = β ( ℓ - ℓ c ) , (5) and can still be described with system-independent parameters ℓ c ≈ 0.654 and β ≈ 0.791 (see Supplemental Material for data [33] ). These parameters can be understood by analyzing the extreme cases: Near the ground state, the domain shapes do not deviate much from regular polygons, for which the ratio of perimeter to A 1 / 2 is essentially constant [39] . An average polygon undergoing a T 1 with L ≈ L c then has a shorter perimeter by the factor ( 5 + L c ) / 6 ; with Eq. (4) this leads to the estimate ℓ c , a = 6 L c / ( 5 + L c ) ≈ 0.654 . Conversely, any T 1 with ℓ → 0 must have the same result as L → 0 , so that α L c = β ℓ c , resulting in β a ≈ 0.772 . The analytical estimate again proves very accurate [Fig. 2(e) ]. Our SE simulations establish metastable states after every T 1 of a selected edge. Figure 3(a) compares the energy reduction Δ E tot ( n ) = E ( n ) - E ( 0 ) after n T 1 s using different selection strategies: systematic cycling through a complete list of edges, random selection, and the “greedy” algorithm suggested by Eq. (5) , which always flips the edge with the shortest current ℓ . All algorithms reverse T 1 s with Δ E > 0 , and try a different edge next; they all asymptote to very similar energies (supporting the notion of a well-defined ground state energy), but the greedy algorithm needs much less computational effort (and its final energy is slightly lower). These findings are independent of polydispersity or energy functional. 3 10.1103/PhysRevLett.120.248001.f3 FIG. 3. (a) Decrease of energy with the number of T 1 transitions for three different algorithms. The “greedy” algorithm always flips the current shortest edge and yields a very good approximation to the ground state energy in a small number of steps. (b) Evolution of the relative edge length probability distribution p n ( ℓ ) for a polydisperse ( c A = 0.4 ) foam sample of N = 900 , showing distributions at n = 0 [ p ( ℓ ) ], n = 88 , and n = n fin = 324 . Analytically, a total energy drop as in Fig. 3(a) can be predicted under the assumption that the effects of the (spatially localized) T 1 s are independent. Then, E tot ( n ) can be inferred from the initial probability distribution p ( ℓ ) only—precisely those edges with ℓ < ℓ c should flip. In Fig. 3(b) , a typical development of p ( ℓ ) with n in simulations is shown—indeed, the probability weight below ℓ c becomes negligible towards the end. Then, an edge of length ℓ flips after n ( ℓ ) T 1 s, such that n ( ℓ ) = 3 N ∫ 0 ℓ p ( ℓ ′ ) d ℓ ′ , (6) where 3 N is the total number of edges, and the predicted final number of T 1 s is n fin = n ( ℓ c ) . It follows that Δ E tot ( n ) = ∫ 0 ℓ ( n ) 3 N Δ E ( ℓ ) p ( ℓ ) d ℓ , (7) where ℓ ( n ) is given by inverting Eq. (6) . Taking into account Eq. (5) and using integration by parts, it is easy to show that Δ E tot ( n ) 3 N β ℓ c = ℓ ( n ) - ℓ c ℓ c P ( ℓ ( n ) ) - 1 ℓ c ∫ 0 ℓ ( n ) P ( ℓ ) d ℓ , (8) with P ( ℓ ) ≡ ∫ 0 ℓ p ( ℓ ′ ) d ℓ ′ . Equation (8) gives the predicted energy decrease as a fraction of a hypothetical maximum; note that - β ℓ c = Δ E ( 0 ) according to Eq. (5) . The prediction Eq. (8) only needs the initial distribution p ( ℓ ) for ℓ < ℓ c ; any integrable fit to p ( ℓ ) yields an explicit analytical expression for Δ E tot . Figure 4 compares greedy simulation results of different foams and tissues (only energy-lowering steps are accepted) with Eq. (8) ; for tissue systems, equivalent foam energies E f * are again used. The agreement is good, but | Δ E tot | is systematically underestimated by typically 5%–15%. 4 10.1103/PhysRevLett.120.248001.f4 FIG. 4. Energy decrease with T 1 number n comparing the simulation results (red; greedy algorithm) with theoretical predictions from Eq. (8) (orange, dashed line) and the refined theory Eq. (9) (blue, solid line). (a) Foams of different polydispersity, using E f ; (b) examples of tissues ( E f * , c A = 0.4 ) with γ = 0 and γ = 0.1 . This bias can be eliminated by modeling the shape changes in p n ( ℓ ) shown in Fig. 3(b) ; these come about because T 1 s induce exponentially decaying fluctuations in the absolute or relative lengths of edges beyond a characteristic distance r s [cf. Figs. 1(b) and 1(c) ]. This stochastic fluctuation of width σ ℓ acts as a convolution on p ( ℓ ) , increasing its width and lowering the value of ℓ ( n ) to ℓ ( n ) - Δ ℓ , so that the currently shortest edges become slightly shorter and their T 1 lowers the energy slightly more. The system independence of the features seen in Figs. 1(c) and 2(e) allows for an analytical computation of this convolution in the limit of Gaussian distributions, from which Δ ℓ ≈ 0.037 follows. The details are found in the Supplemental Material [33] . Accordingly, we modify Eq. (8) to Δ E tot ( n ) 3 N β ℓ c = ℓ eff ( n ) - ℓ c ℓ c P ( ℓ ( n ) ) - 1 ℓ c ∫ 0 ℓ ( n ) P ( ℓ ) d ℓ , (9) with ℓ eff = max ( ℓ - Δ ℓ , 0 ) , to avoid negative edge lengths. The systematic error in the comparisons to simulation results is largely eliminated (see Fig. 4 ), though a statistical error of a few percent remains (see Supplemental Material [33] ). The predicted | Δ E tot | is still obtained from the initial distribution only, and thus the asymptotic ground state energy is accurately predicted from just a snapshot of an initial metastable state. We stress that the simulations employ a variety of strategies for annealing to the ground state [33] , which can lead to a larger empirical n fin , but nevertheless, this “single-shot” prediction of Δ E tot is in good agreement. Also note that tissue samples with γ = 0.1 are much closer to the critical γ c than those with γ = 0 , but the quality of the prediction is unchanged. We have demonstrated that the geometry of 2D metastable states quantitatively determines their energy both locally and globally, beyond the trivial summing of edge lengths to obtain a foam energy: locally, the Δ E of a T 1 is predicted by its edge length. Globally, T 1 energies integrate to approximate the metastable state energy above the ground state Δ E tot . Energy-lowering T 1 transitions are almost exclusively confined to edges with relative length ℓ < ℓ c , and the critical value is universal across polydispersities and energy functional forms. Only these edges “store” the structural energy above the ground state, and they are relatively few [we did not find metastable states with P ( ℓ c ) > 0.18 ]. Apart from the foam and tissue models discussed here, we have conducted less extensive simulations with energy functionals including area elasticity, with altered boundary conditions, and even with springlike interactions, without changes to the reported relations. The remarkable simplicity and generality of these findings is reminiscent of the classification of rigid and floppy ground states by the purely geometric shape index p ¯ = ∑ i P i / A i 1 / 2 / N of the domains, independent of energy functionals [18,19,40] . Likewise, we rationalize the universal nature of our results by the strong geometric constraints imposed by a space-filling 2D structure with rigid domain boundaries, making all possible T 1 energy changes perturbative. Beyond the loss of rigidity transition, domains acquire different geometric degrees of freedom, and the relations cease to be valid. Importantly, however, our results show that metastable states of rigid or solid systems can have values of p ¯ (i.e., equivalent foam energies) significantly larger than the critical p ¯ c for loss of rigidity in the ground state . While there are other indicators of floppy or fluid systems, this illustrates that the diagnostic meaning of p ¯ depends on whether the system is close to the ground state. According to our results, a simple snapshot of any metastable 2D sample (a tissue, an emulsion, a polycrystal) in a rigid or solid state suffices to classify it in terms of its distance from the type of ground state analyzed in previous work [19,25] . Short edges are weak spots favoring T 1 transitions, and a concentration of short edges indicates mechanically weak regions. The diagnostics of material properties and their spatial distribution (in industrial applications) or the occurrence and distribution of pathological changes (in biological tissues) is aided by these findings. The geometric information used here can be further combined with topological statistics [34,41,42] , which is the subject of ongoing work [43] . We thank Ken Brakke for valuable advice on surface evolver code implementation, and the NSF under Grant No. 1504301 for support. [1] 1 C. h. Liu , S. R. Nagel , D. A. Schecter , S. N. Coppersmith , S. Majumdar , O. Narayan , and T. A. Witten , Science 269 , 513 ( 1995 ). SCIEAS 0036-8075 10.1126/science.269.5223.513 [2] 2 H. M. Jaeger , S. R. Nagel , and R. P. Behringer , Rev. Mod. Phys. 68 , 1259 ( 1996 ). RMPHAT 0034-6861 10.1103/RevModPhys.68.1259 [3] 3 D. Lohse , R. Rauhé , R. Bergmann , and D. van der Meer , Nature (London) 432 , 689 ( 2004 ). NATUAS 0028-0836 10.1038/432689a [4] 4 T. S. Majmudar and R. P. Behringer , Nature (London) 435 , 1079 ( 2005 ). NATUAS 0028-0836 10.1038/nature03805 [5] 5 A. Donev , I. Cisse , D. Sachs , E. A. Variano , F. H. Stillinger , R. Connelly , S. Torquato , and P. M. Chaikin , Science 303 , 990 ( 2004 ). SCIEAS 0036-8075 10.1126/science.1093010 [6] 6 W. Man , A. Donev , F. H. Stillinger , M. T. Sullivan , W. B. Russel , D. Heeger , S. Inati , S. Torquato , and P. M. Chaikin , Phys. Rev. Lett. 94 , 198001 ( 2005 ). PRLTAO 0031-9007 10.1103/PhysRevLett.94.198001 [7] 7 T. Aste and D. Weaire , The Pursuit of Perfect Packing ( Institute of Physics , London, 2000 ). [8] 8 S. Torquato , T. M. Truskett , and P. G. Debenedetti , Phys. Rev. Lett. 84 , 2064 ( 2000 ). PRLTAO 0031-9007 10.1103/PhysRevLett.84.2064 [9] 9 S. Torquato and F. H. Stillinger , Rev. Mod. Phys. 82 , 2633 ( 2010 ). RMPHAT 0034-6861 10.1103/RevModPhys.82.2633 [10] 10 S. Atkinson , F. H. Stillinger , and S. Torquato , Proc. Natl. Acad. Sci. U.S.A. 111 , 18436 ( 2014 ). PNASA6 0027-8424 10.1073/pnas.1408371112 [11] 11 A. J. Liu and S. R. Nagel , Annu. Rev. Condens. Matter Phys. 1 , 347 ( 2010 ). ARCMCX 1947-5454 10.1146/annurev-conmatphys-070909-104045 [12] 12 C. E. Zachary , Y. Jiao , and S. Torquato , Phys. Rev. Lett. 106 , 178001 ( 2011 ). PRLTAO 0031-9007 10.1103/PhysRevLett.106.178001 [13] 13 S. Hilgenfeldt , A. M. Kraynik , S. A. Koehler , and H. A. Stone , Phys. Rev. Lett. 86 , 2685 ( 2001 ). PRLTAO 0031-9007 10.1103/PhysRevLett.86.2685 [14] 14 S. Hilgenfeldt , A. M. Kraynik , D. A. Reinelt , and J. M. Sullivan , Europhys. Lett. 67 , 484 ( 2004 ). EULEEJ 0295-5075 10.1209/epl/i2003-10295-7 [15] 15 D. Weaire and S. Hutzler , The Physics of Foams ( Oxford University Press , Oxford, 2000 ). [16] 16 S. Hilgenfeldt , S. Erisken , and R. W. Carthew , Proc. Natl. Acad. Sci. U.S.A. 105 , 907 ( 2008 ). PNASA6 0027-8424 10.1073/pnas.0711077105 [17] 17 M. L. Manning , R. A. Foty , M. S. Steinberg , and E.-M. Schoetz , Proc. Natl. Acad. Sci. U.S.A. 107 , 12517 ( 2010 ). PNASA6 0027-8424 10.1073/pnas.1003743107 [18] 18 S. Kim and S. Hilgenfeldt , Soft Matter 11 , 7270 ( 2015 ). SMOABF 1744-683X 10.1039/C5SM01563D [19] 19 D. Bi , J. H. Lopez , J. M. Schwarz , and M. L. Manning , Nat. Phys. 11 , 1074 ( 2015 ). NPAHAX 1745-2473 10.1038/nphys3471 [20] 20 A. Fletcher , M. Osterfield , R. Baker , and S. Shvartsman , Biophys. J. 106 , 2291 ( 2014 ). BIOJAU 0006-3495 10.1016/j.bpj.2013.11.4498 [21] 21 G. W. Brodland and J. H. Veldhuis , J. Biomech. 35 , 673 ( 2002 ). JBMCB5 0021-9290 10.1016/S0021-9290(02)00006-4 [22] 22 G. W. Brodland , D. I.-L. Chen , and J. H. Veldhuis , Int. J. Plast. 22 , 965 ( 2006 ). IJPLER 0749-6419 10.1016/j.ijplas.2005.05.002 [23] 23 G. Zhang , F. H. Stillinger , and S. Torquato , Phys. Rev. E 96 , 042146 ( 2017 ). PRESCM 2470-0045 10.1103/PhysRevE.96.042146 [24] 24 M. Alava , P. Duxbury , C. Moukarzel , and H. Rieger , Phase Transitions Crit. Phenom. 18 , 143 ( 2001 ). PTCPEJ 1062-7901 10.1016/S1062-7901(01)80009-4 [25] 25 D. Staple , R. Farhadifar , J. C. Röper , B. Aigouy , S. Eaton , and F. Jülicher , Eur. Phys. J. E 33 , 117 ( 2010 ). EPJSFH 1292-8941 10.1140/epje/i2010-10677-0 [26] 26 R. Farhadifar , J.-C. Röper , B. Aigouy , S. Eaton , and F. Jülicher , Curr. Biol. 17 , 2095 ( 2007 ). CUBLE2 0960-9822 10.1016/j.cub.2007.11.049 [27] 27 M. Krajnc , S. Dasgupta , P. Ziherl , and J. Prost , arXiv:1805.06500 . [28] 28 D. Bi , J. H. Lopez , J. M. Schwarz , and M. L. Manning , Soft Matter 10 , 1885 ( 2014 ). SMOABF 1744-683X 10.1039/c3sm52893f [29] 29 D. Weaire and N. Rivier , Contemp. Phys. 25 , 59 ( 1984 ). CTPHAF 0010-7514 10.1080/00107518408210979 [30] 30 J. E. Taylor , Ann. Math. 103 , 489 ( 1976 ). ANMAAH 0003-486X 10.2307/1970949 [31] 31 K. Brakke , Exp. Math. 1 , 141 ( 1992 ). 1058-6458 10.1080/10586458.1992.10504253 [32] 32 S. Cox , F. Graner , and M. F. Vaz , Soft Matter 4 , 1871 ( 2008 ). SMOABF 1744-683X 10.1039/b802792g [33] 33a See Supplemental Material at http://link.aps.org/supplemental/10.1103/PhysRevLett.120.248001 for details; 33b see also H. X. Zhu , S. M. Thorpe , and A. H. Windle , Philos. Mag. A 81 , 2765 ( 2001 ); PMAADG 0141-8610 10.1080/01418610010032364 33c S. Lloyd , IEEE Trans. Inf. Theory 28 , 129 ( 1982 ). IETTAW 0018-9448 10.1109/TIT.1982.1056489 [34] 34 C. Quilliet , S. A. Talebi , D. Rabaud , S. C. J. Käfer , and F. Graner , Philos. Mag. Lett. 88 , 651 ( 2008 ). PMLEEG 0950-0839 10.1080/09500830802334249 [35] 35 A. Wyn , I. Davies , and S. Cox , Eur. Phys. J. E 26 , 81 ( 2008 ). EPJSFH 1292-8941 10.1140/epje/i2007-10286-0 [36] 36 D. Chen , K. W. Desmond , and E. R. Weeks , Soft Matter 8 , 10486 ( 2012 ). SMOABF 1744-683X 10.1039/c2sm26023a [37] 37 Y. Gai , C. M. Leong , W. Cai , and S. K. Tang , Proc. Natl. Acad. Sci. U.S.A. 113 , 12082 ( 2016 ). PNASA6 0027-8424 10.1073/pnas.1606601113 [38] 38 C. Collinet , M. Rauzi , P.-F. Lenne , and T. Lecuit , Nat. Cell Biol. 17 , 1247 ( 2015 ). NCBIFN 1465-7392 10.1038/ncb3226 [39] 39 F. Graner , Y. Jiang , E. Janiaud , and C. Flament , Phys. Rev. E 63 , 011402 ( 2000 ). 10.1103/PhysRevE.63.011402 [40] 40 J.-A. Park , J. H. Kim , D. Bi , J. A. Mitchel , N. T. Qazvini , K. Tantisira , C. Y. Park , M. McGill , S.-H. Kim , B. Gweon , Nat. Mater. 14 , 1040 ( 2015 ). NMAACR 1476-1122 10.1038/nmat4357 [41] 41 S. Kim , J. J. Cassidy , B. Yang , R. W. Carthew , and S. Hilgenfeldt , Biophys. J. 111 , 2735 ( 2016 ). BIOJAU 0006-3495 10.1016/j.bpj.2016.11.004 [42] 42 M. P. Miklius and S. Hilgenfeldt , Phys. Rev. Lett. 108 , 015502 ( 2012 ). PRLTAO 0031-9007 10.1103/PhysRevLett.108.015502 [43] 43 S. Kim and S. Hilgenfeldt , Publisher Copyright: © 2018 American Physical Society.

PY - 2018/6/11

Y1 - 2018/6/11

N2 - Mechanical equilibrium states of cellular matter are overwhelmingly metastable and separated from each other by topology changes. Using theory and simulations, it is shown that for a wide class of energy functionals in 2D, including those describing tissue cell layers, local energy differences between neighboring metastable states as well as global energy differences between initial states and ground states are governed by simple, universal relations. Knowledge of instantaneous length of an edge undergoing a T1 transition is sufficient to predict local energy changes, while the initial edge length distribution yields a successful prediction for the global energy difference. An analytical understanding of the model parameters is provided.

AB - Mechanical equilibrium states of cellular matter are overwhelmingly metastable and separated from each other by topology changes. Using theory and simulations, it is shown that for a wide class of energy functionals in 2D, including those describing tissue cell layers, local energy differences between neighboring metastable states as well as global energy differences between initial states and ground states are governed by simple, universal relations. Knowledge of instantaneous length of an edge undergoing a T1 transition is sufficient to predict local energy changes, while the initial edge length distribution yields a successful prediction for the global energy difference. An analytical understanding of the model parameters is provided.

UR - http://www.scopus.com/inward/record.url?scp=85048566390&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85048566390&partnerID=8YFLogxK

U2 - 10.1103/PhysRevLett.120.248001

DO - 10.1103/PhysRevLett.120.248001

M3 - Article

C2 - 29957000

AN - SCOPUS:85048566390

SN - 0031-9007

VL - 120

JO - Physical review letters

JF - Physical review letters

IS - 24

M1 - 248001

ER -